Volatile fluxes from volcanoes

Volatile fluxes from Mid Ocean Ridge (MOR) and subaerial volcanism have been estimated or re-evaluated using several natural tracers-³He, ²¹⁰Po, SO₂-and chemical ratios of volatile species in lavas and volcanic gases.
These estimates confirm the net predominance of anthropogenic fluxes over volcanic fluxes for CO₂, SO₂ and trace metals.
They also suggest that, while most of the volatiles transferred during MOR volcanism come from the mantle, volatiles stored at the surface of the Earth supply an appreciable fraction of subaerial fluxes and can be the dominant source for some of them.
The surface inventory of volatile species cannot result from steady-state degassing with constant rate and needs much greater fluxes in the past or other volatile supply processes. This inventory is the result of several of the following processes: capture of the solar nebula and its subsequent partial escape, impact degassing of accreting bodies, and, from Archean to present mantle, degasssing through volcanism and associated phenomena, with recycling into the mantle through subduction.     

Downstream composition changes of acidic volcanic waters discharged into the Banyupahit stream, Ijen caldera, Indonesia

The crater lake of Kawah Ijen volcano contains extremely low pH (<0.4) waters with high SO₄ (70000 mg/kg), Cl (21000 mg/kg), F (1500 mg/kg), Al (5000 mg/kg), Fe (2000 mg/kg) and trace metal (Cu 0.5, Zn 4, Pb 3 mg/kg) contents. These brines seep outward through the western crater rim and reappear on the other side as streamlets, which form the headwaters of the Banyupahit stream. The Banyupahit first mixes with fresh rivers and thermal springs in the Ijen caldera and then irrigates a coastal agricultural plain which is 30 km from the summit crater.We discuss the downstream composition changes affecting the Banyupahit waters by using stable isotope, chemical and mineralogical data collected from sites along the stream length. The saturation of the stream waters with respect to minerals was evaluated with SOLVEQ and WATEQ4F and compared with the geochemical observations. An aluminous mineralogy (alunogen, pickeringite, tamarugite and kalinite) develops in the upper part of the Banyupahit due to concentration of the headwaters by evaporation. Downstream attenuation of dissolved element concentrations results principally from dilution and from mineral precipitation. The stream pH changes from 0 at the source to >4 close to the mouth. The δD and δ¹⁸O values and the relative SO₄–Cl–F contents of the Banyupahit waters indicate that the tributaries are mostly meteoric. Dissolved SO₄ in the acidic stream come only from the crater lake seepages and are not involved later in microbially mediated reactions, as shown by their δ³⁴S and δ¹⁸O values. Re-equilibration of the stream SO₄ oxygen-isotope composition with H₂O from tributaries does not occur.Calcium, SiO₂, Al, Fe, K and SO₄ behave non-conservatively in the stream waters. Gypsum, silica (amorphous or poorly ordered), a basic aluminum hydroxysulfate (basaluminite?), K-jarosite and amorphous ferric hydroxide may exert a solubility control on these elements along the entire stream length, or in certain stream sections, consistent with the thermochemical model results. Downstream concentration trends and mineral saturation levels suggest that precipitation of Sr-, Pb-rich barite and celestite consume Ba, Sr and Pb, whereas dissolved Cu, Pb and Zn may adsorb onto solid particles, especially after the junctions of the acidic stream with non-acidic rivers. We calculated that significant fluxes of SO₄, F, Cl, Al, SiO₂, Ti, Mn and Cu may reach the irrigation system, possibly causing serious environmental impacts such as soil acidification and induration.     

Practical method for drawing a VGP path

A statistical method is proposed for the smoothing of polar wander paths and for giving their confidence limits by the mean of successive ellipses. The method is fully parametrical and is based on the relations between the inertia matrix and the parameters of the Fisher distribution, from which a bivariate form is deduced. An elementary tensorial calculation gives the parameters of the confidence ellipse around a vectorial weighted mean, for any unimodal set of vectors. This model can also be used for other statistical tests, wherever the rotational symmetry hypothesis is not consistent (i.e. the fold test).     

The synergistic powers of AEM and GIS geodatabase models in water resources studies

Comprehensive studies of water resources systems require integration of modeling tools and data associated with individual processes. An object-oriented approach is presented here that associates ground water models based upon the analytic element method (AEM) with geographic information system (GIS) geodatabase features using an AEM Model Interface. Each aquifer object contains a prescribed geometry, a mathematical representation in the AEM, and GIS hydrogeologic data. The synergistic understanding inherent in such an approach is illustrated by a study linking local AEM model predictions of water elevation with ground water geodatabase objects. This AEM Model Interface provides a key component in establishing a common object-oriented geodatabase modeling approach linking ground water to a variety of natural and social processes.     

Distinguishing volcanic debris avalanche deposits from their reworked products: the Perrier sequence (French Massif Central)

Debris avalanches associated with volcanic sector collapse are usually high-volume high-mobility phenomena. Debris avalanche deposit remobilisation by cohesive debris flows and landslides is common, so they can share textural characteristics such as hummocks and jigsaw cracks. Distinguishing original deposits from reworked products is critical for geological understanding and hazard assessment because of their different origin, frequency and environmental impact. We present a methodology based on field evidence to differentiate such epiclastic breccias. Basal contact mapping constrained by accurate altitude and location data allows the reconstruction of deposit stratigraphy and geometry. Lithological analysis helps to distinguish the different units. Incorporation structures, kinematic indicators and component mingling textures are used to characterise erosion and transport mechanisms. We apply this method to the enigmatic sequence at Perrier (French Massif Central), where four units (U1–U4) have been interpreted either as debris flow or debris avalanche deposits. The sequence results from activity on the Monts Dore Volcano about 2 Ma ago. The epiclastic units are matrix supported with an almost flat top. U2 and U3 have clear debris flow deposit affinities such as rounded clasts and intact blocks (no jigsaw cracks). U1 and U4 have jigsaw cracked blocks with matrix injection and stretched sediment blocks. U1 lacks large blocks (>10 m wide) and has a homogenous matrix with an upward increase of trapped air vesicle content and size. This unit is interpreted as a cohesive debris flow deposit spawned from a debris avalanche upstream. In contrast, U4 has large mega-blocks (up to 40 m wide), sharp contacts between mixed facies zones with different colours and numerous jigsaw fit blocks (open jigsaw cracks filled by monogenic intra-clast matrix). Mega-blocks are concentrated near the deposit base and are spatially associated with major substratum erosion. This deposit has a debris avalanche distal facies with local debris flow affinities due to partial water saturation. We also identify two landslide deposits (L1 and L2) resulting from recent reworking that has produced a similar facies to U1 and U4. These are distinguishable from the original deposits, as they contain blocks of mixed U1/U4 facies, a distinctly less consolidated and more porous matrix and a fresh hummocky topography. This work shows how to differentiate epiclastic deposits with similar characteristics, but different origins. In doing so, we improve understanding of present and past instability of the Monts Dore and identify present landslide hazards at Perrier.     

Melting of enriched mantle beneath Pitcairn seamounts: Unusual U–Th–Ra systematics provide insights into melt extraction processes

U-series systematics as well as Sr isotopes were measured on young seamount lavas from the Pitcairn hotspot collected during the Polynaut cruise. The combined U-series and Sr isotope data reveal typical mixing relationships between two endmembers. One typical ‘plume’ endmember with radiogenic ⁸⁷Sr/⁸⁶Sr and relatively low ²³⁰Th/²³⁸U and a ‘lithosphere’ endmember with less radiogenic ⁸⁷Sr/⁸⁶Sr and relatively larger ²³⁰Th/²³⁸U. Remarkably, all the lavas, except for a few arguably older samples, are characterized by ²²⁶Ra deficits relative to ²³⁰Th. On the basis of water content and trace element systematics, we argue that this is due to melting in the presence of phlogopite, which is only stable at lithosphere temperatures. A melting model including the diffusive exchange of elements among phlogopite, garnet and melt is used to constrain melting conditions of the lithosphere. These unusual ²²⁶Ra–²³⁰Th signatures can be explained by relatively slow melting rates at low matrix porosity. Our model also demonstrates that the effective partitioning behavior is dependent on the melting rate. A simple thermal model for lithosphere heating and melting is in good agreement with predicted melting rates.