Mafic Plinian volcanism and ignimbrite emplacement at Tofua volcano, Tonga

Tofua Island is the largest emergent mafic volcano within the Tofua arc, Tonga, southwest Pacific. The volcano is dominated by a distinctive caldera averaging 4 km in diameter, containing a freshwater lake in the south and east. The latest paroxysmal (VEI 5-6) explosive volcanism includes two phases of activity, each emplacing a high-grade ignimbrite. The products are basaltic andesites with between 52 wt.% and 57 wt.% SiO(2). The first and largest eruption caused the inward collapse of a stratovolcano and produced the 'Tofua' ignimbrite and a sub-circular caldera located slightly northwest of the island's centre. This ignimbrite was deposited in a radial fashion over the entire island, with associated Plinian fall deposits up to 0.5 m thick on islands > 40 km away. Common sub-rounded and frequently cauliform scoria bombs throughout the ignimbrite attest to a small degree of marginal magma-water interaction. The common intense welding of the coarse-grained eruptive products, however, suggests that the majority of the erupted magma was hot, water-undersaturated and supplied at high rates with moderately low fragmentation efficiency and low levels of interaction with external water. We propose that the development of a water-saturated dacite body at shallow (<6 km) depth resulted in failure of the chamber roof to cause sudden evacuation of material, producing a Plinian eruption column. Following a brief period of quiescence, largescale faulting in the southeast of the island produced a second explosive phase believed to result from recharge of a chemically distinct magma depleted in incompatible elements. This similar, but smaller eruption, emplaced the 'Hokula' Ignimbrite sheet in the northeast of the island. A maximum total volume of 8 km(3) of juvenile material was erupted by these events. The main eruption column is estimated to have reached a height of similar to 12 km, and to have produced a major atmospheric injection of gas, and tephra recorded in the widespread series of fall deposits found on coral islands 40-80 km to the east (in the direction of regional upper-tropospheric winds). Radiocarbon dating of charcoal below the Tofua ignimbrite and organic material below the related fall units imply this eruption sequence occurred post 1,000 years BP. We estimate an eruption magnitude of 2.24x10(13) kg, sulphur release of 12 Tg and tentatively assign this eruption to the AD 1030 volcanic sulphate spike recorded in Antarctic ice sheet records.

Regulation of arsenic mobility on basaltic glass surfaces by speciation and pH

The importance of geothermal energy as a source for electricity generation and district heating has increased over recent decades. Arsenic can be a significant constituent of the geothermal fluids pumped to the surface during power generation. Dissolved As exists in different oxidation states, mainly as As(III) and As(V), and the charge of individual species varies with pH. Basaltic glass is one of the most important rock types in many high-temperature geothermal fields. Static batch and dynamic column experiments were combined to generate and validate sorption coefficients for As(III) and As(V) in contact with basaltic glass at pH 3-10. Validation was carried out by two empirical kinetic models and a surface complexation model (SCM). The SCM provided a better fit to the experimental column data than kinetic models at high pH values. However, in certain circumstances, an adequate estimation of As transport in the column could not be attained without incorporation of kinetic reactions. The varying mobility with pH was due to the combined effects of the variable charge of the basaltic glass with the pH point of zero charge at 6.8 and the individual As species as pH shifted, respectively. The mobility of As(III) decreased with increasing pH. The opposite was true for As(V), being nearly immobile at pH 3 to being highly mobile at pH 10. Incorporation of appropriate sorption constants, based on the measured pH and Eh of geothermal fluids, into regional groundwater-flow models should allow prediction of the As(III) and As(V) transport from geothermal systems to adjacent drinking water sources and ecosystems.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Petrogenesis of Basaltic Andesites at Tofua Volcano

Tofua volcano is situated midway along the Tonga oceanic arc and has undergone two phases of ignimbrite-forming activity. The eruptive products are almost entirely basaltic andesites (52 center dot 5-57 wt % SiO2) with the exception of a volumetrically minor pre-caldera dacite. The suite displays a strong tholeiitic trend with K2O <1 wt %. Phenocryst assemblages typically comprise plagioclase + clinopyroxene +/- orthopyroxene with microlites of Ti-magnetite. Olivine (Fo(83-88)) is rare and believed to be dominantly antecrystic. An increase in the extent and frequency of reverse zoning in phenocrysts, sieve-textured plagioclase and the occurrence of antecrystic phases in post-caldera lavas record a shift to dynamic conditions, allowing the interaction of magma batches that were previously distinct. Pyroxene thermobarometry suggests crystallization at 950-1200 degrees C and 0 center dot 8-1 center dot 8 kbar. Volatile measurements of glassy melt inclusions indicate a maximum H2O content of 4 center dot 16 wt % H2O, and CO2-H2O saturation curves indicate that crystallization occurred at two levels, at depths of 4-5 center dot 5 km and 1 center dot 5-2 center dot 5 km. Major and trace element models suggest that the compositions of the majority of the samples represent a differentiation trend whereby the dacite was produced by 65% fractional crystallization of the most primitive basaltic andesite. Trace element models suggest that the sub-arc mantle source is the residuum of depleted Indian mid-ocean ridge basalt mantle (IDMM-1% melt), whereas radiogenic isotope data imply addition of 0 center dot 2% average Tongan sediment melt and a fluid component derived from the subducted altered Pacific oceanic crust. A horizontal array on the U-Th equiline diagram and Ra excesses of up to 500% suggest fluid addition to the mantle wedge within the last few thousand years. Time-integrated (Ra-226/Th-230) vs Sr/Th and Ba/Th fractionation models imply differentiation timescales of up to 4500 years for the dacitic magma compositions at Tofua.

A field and reactive transport model study of arsenic in a basaltic rock aquifer

The use of geothermal energy as a source for electricity and district heating has increased over recent decades. Dissolved As can be an important constituent of the geothermal fluids brought to the Earth's surface. Here the field application of laboratory measured adsorption coefficients of aqueous As species on basaltic glass surfaces is discussed. The mobility of As species in the basaltic aquifer in the Nesjavellir geothermal system, Iceland was modelled by the one-dimensional (1D) reactive transport model PHREEQC ver. 2, constrained by a long time series of field measurements with the chemical composition of geothermal effluent fluids, pH, Eh and, occasionally, Fe- and As-dissolved species measurements. Di-, tri- and tetrathioarsenic species (As(OH)S22-, AsS₃H^2-, AsS33- and As(SH)4-) were the dominant form of dissolved As in geothermal waters exiting the power plant (2.556μM total As) but converted to some extent to arsenite (H₃AsO₃) and arsenate HAsO42- oxyanions coinciding with rapid oxidation of S2- to S2O32- and finally to SO42- during surface runoff before feeding into a basaltic lava field with a total As concentration of 0.882μM following dilution with other surface waters. A continuous 25-a data set monitoring groundwater chemistry along a cross section of warm springs on the Lake Thingvallavatn shoreline allowed calibration of the 1D model. Furthermore, a series of ground water wells located in the basaltic lava field, provided access along the line of flow of the geothermal effluent waters towards the lake. The conservative ion Cl^- moved through the basaltic lava field (4100m) in less than10a but As was retarded considerably due to surface reactions and has entered a groundwater well 850m down the flow path as arsenate in accordance to the prediction of the 1D model. The 1D model predicted a complete breakthrough of arsenate in the year 2100. In a reduced system arsenite should be retained for about 1ka. © 2011 Elsevier Ltd.

The spatial distribution of Holocene cryptotephras in north-west Europe since 7 ka: implications for understanding ash fall events from Icelandic eruptions

We present distribution maps for all cryptotephras (distal volcanic ash layers) younger than 7 ka that have been reported from three or more lakes or peatlands in north-west Europe. All but one of the tephras originates from Iceland; the exception has been attributed to Jan Mayen. We find strong spatial patterning in tephra occurrence at the landscape scale; most, but not all of the tephra occurrences are significantly spatially clustered, which likely reflects atmospheric and weather patterns at the time of the eruptions. Contrary to expectations based on atmospheric modelling studies, tephras appear to be at least as abundant in Ireland and northern Scotland as in Scandinavia. Rhyolitic and other felsic tephras occur in lakes and peatlands throughout the study region, but andesitic and basaltic tephras are largely restricted to lake sites in the Faroe Islands and Ireland. Explanations of some of these patterns will require further research on the effects of different methodologies for locating and characterizing cryptotephras. These new maps will help to guide future investigations in tephrochronology and volcanic hazard analysis.

Middle Stone Age (MSA) site distributions in eastern Africa and their relationship to Quaternary environmental change, refugia and the evolution of Homo sapiens

This paper considers the evolution of Homo sapiens in eastern Africa in relation to refugia and bottlenecks around ~200 ka BP, at a macro scale. Middle Stone Age (MSA) lithics, site distributions and locations are analysed in relation to palaeovegetation maps of the last glacial/interglacial cycle, which are used as a proxy for earlier climate cycles. A ‘‘push and pull’’ model is then postulated for the spread of Homo sapiens out of refugia in eastern Africa, involving both volcanism (push) and habitat availability (pull). A date within OIS 5 is suggested for this expansion to other parts of the continent, and potentially further a?eld, contrary to a frequently proposed expansion within OIS 3. ©2008 Elsevier Ltd. All rights reserved.

Early Holocene M similar to 6 explosive eruption from Plosky volcanic massif (Kamchatka) and its tephra as a link between terrestrial and marine paleoenvironmental records

We report tephrochronological and geochemical data on early Holocene activity from Plosky volcanic massif in the Kliuchevskoi volcanic group, Kamchatka Peninsula. Explosive activity of this volcano lasted for similar to 1.5 kyr, produced a series of widely dispersed tephra layers, and was followed by profuse low-viscosity lava flows. This eruptive episode started a major reorganization of the volcanic structures in the western part of the Kliuchevskoi volcanic group. An explosive eruption from Plosky (M similar to 6), previously unstudied, produced tephra (coded PL2) of a volume of 10-12 km(3) (11-13 Gt), being one of the largest Holocene explosive eruptions in Kamchatka. Characteristic diagnostic features of the PL2 tephra are predominantly vitric sponge-shaped fragments with rare phenocrysts and microlites of plagioclase, olivine and pyroxenes, medium- to high-K basaltic andesitic bulk composition, high-K, high-Al and high-P trachyandesitic glass composition with SiO2 = 57.5-59.5 wt%, K2O = 2.3-2.7 wt%, Al2O3 = 15.8-16.5 wt%, and P2O5 = 0.5-0.7 wt%. Other diagnostic features include a typical subduction-related pattern of incompatible elements, high concentrations of all REE (> 10x mantle values), moderate enrichment in LREE (La/Yb similar to 5.3), and non-fractionated mantle-like pattern of LILE. Geochemical fingerprinting of the PL2 tephra with the help of EMP and LA-ICP-MS analyses allowed us to map its occurrence in terrestrial sections across Kamchatka and to identify this layer in Bering Sea sediment cores at a distance of > 600 km from the source. New high-precision C-14 dates suggest that the PL2 eruption occurred similar to 10,200 cal BP, which makes it a valuable isochrone for early Holocene climate fluctuations and permits direct links between terrestrial and marine paleoenvironmental records. The terrestrial and marine C-14 dates related to the PL2 tephra have allowed us to estimate an early Holocene reservoir age for the western Bering Sea at 1,410 +/- A 64 C-14 years. Another important tephra from the early Holocene eruptive episode of Plosky volcano, coded PL1, was dated at 11,650 cal BP. This marker is the oldest geochemically characterized and dated tephra marker layer in Kamchatka to date and is an important local marker for the Younger Dryas-early Holocene transition. One more tephra from Plosky, coded PL3, can be used as a marker northeast of the source at a distance of similar to 110 km.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Rapid Petrogenesis of Dacites at Fonualei Volcano

Fonualei is unusual amongst subaerial volcanoes in the Tonga arc because it has erupted dacitic vesicular lavas, tuffs and phreomagmatic deposits for the last 165 years. The total volume of dacite may approach 5 km(3) and overlies basal basaltic andesite and andesite lavas that are constrained to be less than a few millennia in age. All of the products are crystal-poor and formed from relatively low-viscosity magmas inferred to have had temperatures of 1100-1000 degrees C, 2-4 wt % H2O and oxygen fugacities 1-2 log units above the quartz-fayalite-magnetite buffer. Major and trace element data, along with Sr-Nd-Pb and U-Th-Ra isotope data, are used to assess competing models for the origin of the dacites. Positive correlations between Sc and Zr and Sr rule out evolution of the within-dacite compositional array by closed-system crystal fractionation of a single magma batch. An origin by partial melting of lower crustal amphibolites cannot reproduce these data trends or, arguably, any of the dacites either. Instead, we develop a model in which the dacites reflect mixing between two dacitic magmas, each the product of fractional crystallization of basaltic andesite magmas formed by different degrees of partial melting. Mixing was efficient because the two magmas had similar temperatures and viscosities. This is inferred to have occurred at shallow (2-6 km) depths beneath the volcano. U-Th-Ra disequilibria in the basaltic andesite and andesite indicate that the parental magmas had fluids added to their mantle source regions less than 8 kyr ago and that fractionation to the dacitic compositions took less than a few millennia. The 165 year eruption period for the dacites implies that mixing occurred on a similar timescale, possibly during ascent in conduits. The composition of the dacites renders them unsuitable candidates as contributors to average continental crust.

Recent contribution of sediments and fluids to the mantle's volatile budget

Subduction modifies the cycling of Earth's volatile elements. Fluid-rich sediments and hydrated oceanic lithosphere enter the convecting mantle at subduction zones. Some of the sediments and volatile components are released from the subducting slab, promote mantle melting and are returned to the surface by volcanism. The remainder continue into the deeper mantle. Quantification of the fate of these volatiles requires an understanding of both the nature and timing of fluid release and mantle melting(1). Here we analyse the trace element and isotopic geochemistry of fragments of upper mantle rocks that were transported to the surface by volcanic eruptions above the Batan Island subduction zone, Philippines. We find that the mantle fragments exhibit extreme disequilibrium between their U-Th-Ra isotopic ratios, which we interpret to result from the interaction of wet sediment melts and slab-derived fluids with rocks in the overlying mantle wedge. We infer that wet sediments were delivered from the slab to the mantle wedge between 8,000 and 10,000 years ago, whereas aqueous fluids were delivered separately much later. We estimate that about 625 ppm of water is retained in the wedge. A significant volume of water could therefore be delivered to the mantle transition zone at the base of the upper mantle, or even to the deeper mantle.

Mantle flow, volatiles, slab-surface temperatures and melting dynamics in the north Tonga arc-Lau back-arc basin

The Fonualei Spreading Center affords an excellent opportunity to evaluate geochemical changes with increasing depth to the slab in the Lau back-arc basin. We present H2O and CO2 concentrations and Sr, Nd, Pb, Hf and U-Th-Ra isotope data for selected glasses as well as new Hf isotope data from boninites and seamounts to the north of the Tonga arc. The Pb and Hf isotope data are used to show that mantle flow is oriented to the southwest and that the tear in the northern end of the slab may not extend east as far as the boninite locality. Along the Fonualei Spreading Center, key geochemical parameters change smoothly with increasing distance from the arc front and increasing slab surface temperatures. The latter may range from 720 to 866 degrees C, based on decreasing H2O/Ce ratios. Consistent with experimental data, the geochemical trends are interpreted to reflect changes in the amount and composition of wet pelite melts or super-critical fluids and aqueous fluids derived from the slab. With one exception, all of the lavas preserve both U-238 excesses and Ra-226 excesses. We suggest that lavas from the Fonualei Spreading Center and Valu Fa Ridge are dominated by fluid-fluxed melting whereas those from the East and Central Lau Spreading Centers, where slab surface temperatures exceed similar to 850-900 degrees C, are largely derived through decompression. A similar observation is found for the Manus and East Scotia back-arc basins and may reflect the expiry of a key phase such as lawsonite in the subducted basaltic crust.

Distal volcanic ash layers in the Lateglacial Interstadial (GI-1): problems of stratigraphic discrimination

A new Icelandic ash layer has been detected in mid-Interstadial sediments in a number of Scottish Lateglacial sequences and has been named the Penifiler Tephra. It is rhyolitic in composition and possesses a chemistry, which is similar to the Borrobol Tephra of early Lateglacial Interstadial age, which also occurs in a number of these same sequences. Where the Borrobol Tephra has been identified in these sequences it consistently exhibits a diffuse distribution accompanied in some cases by stratigraphic bimodality. A number of sedimentological and taphonomic factors are considered in order to account for this distribution. One possibility is that these distributions are produced by taphonomic factors. Another possibility is that the Borrobol Tephra may not be the product of a single Icelandic eruption, but of two events closely spaced in time. In at least two of the sequences investigated in this study, basaltic shards were found in association with the Penifiler and Borrobol tephras, suggesting either a basaltic phase associated with these eruptions, or coincident eruptions from a separate basaltic volcanic centre. The discovery of the new Penifiler Tephra makes a contribution to the regional tephrostratigraphic framework, and provides an additional isochron for assessing the synchroneity of palaeoenvironmental changes during the Interstadial. The true stratigraphic nature and age of the Borrobol Tephra, however, remains unresolved and, therefore, its use as an isochron is more problematic. The possible occurrence of basaltic populations may strengthen correlations with basaltic tephras recently detected in the NGRIP ice-core.

Three new distal tephras in sediments spanning the Last Glacial-Interglacial Transition in Scotland

Three new microtephras are reported from a number of lake sites from the Inner Hebrides and Scottish mainland. One occurs stratigrapically in the middle of Greenland Interstadial 1 (GI-1) and has been named the Penifiler Tephra. It is rhyolitic and possesses a geochemical signature that is very similar to that of the Borrobol Tephra, which also occurs in three of the sequences reported here, but which lies close to the lower boundary of GI-1. The second occurs stratigraphically in the early Holocene below the Saksunarvatn Ash and is named the Ashik Tephra. This tephra is geochemically bimodal, with a rhyolitic component comparable to the An Druim Tephra that occurs later in the Holocene, and a basaltic component which is similar to the Saksunarvatn Ash. A third tephra occurs stratigraphically above the Saksunarvatn Ash and is provisionally named the Breakish Tephra. The consistent inter-site correlation demonstrated for these new tephras at several sites enhances the regional tephrostratigraphic framework, and increases the potential for correlating palaeoenvironmental events during GI-1 and the early Holocene. However, the occurrence of multiple tephras with similar geochemistry in close stratigraphic and temporal proximity has implications for the rigour with which tephrostratigraphic investigations must be performed.