Volcanic emissions of metals and halogens from White Island (New Zealand) and Erebus volcano (Antarctica) determined with chemical traps

Volcanic emission rates of As, Sb, Pb, Hg, Se, Cl, and F were determined at Erebus volcano, Antarctica and White Island, New Zealand, using chemical traps. The trace metal fluxes were determined by combining the species to S ratios in the solutions with SO₂ emission rates measured by correlation spectrometry at the two volcanoes. At Erebus volcano, fluxes for the metals Pb and Hg were 2.0 × 10^− 4 and 8.1 × 10^− 6 kg s^− 11, respectively. Fluxes for Cl, F, As, Sb and Se (0.35, 0.15, 2.5 × 10^− 4, 1.2 × 10^− 5, and 4.5 × 10^− 6 kg s^− 1, respectively) agreed within error limits for values determined previously by the LiOH impregnated filter method [Zreda-Gostynska, G., Kyle, P., Finnegan, D., Prestbo, K., 1997. Volcanic gas emissions from Mount Erebus and their impact on the Antarctic environment. Journal of Geophysical Research, 102(B7): 15039–15055.], demonstrating the utility of the chemical trap method. A fall in the As/S ratio from 7 × 10^− 4 in 1997/1999 to 3 × 10^− 4 in 2000 at Erebus coincided with a change in the frequency and style of eruptive activity that may have been due to injection of magma into the system. At White Island, chemical trap data indicated fluxes of Cl = 0.90, F = 0.0079, Pb = 2.7 × 10^− 4, Hg = 1.1 × 10^− 5, As = 1.3 × 10^− 4, Sb = 1.9 × 10^− 5 and Se = 1.5 × 10^− 5 kg s^− 1. Samples collected 600 m downwind of the active crater were comparable to samples collected adjacent to the main gas vent, showing that this method can still be used at some distance from a degassing vent.     

Refinement of the late Quaternary geologic history of Erebus volcano,Antarctica using Ar/Ar and Cl age determinations

Six new ⁴⁰Ar/³⁹Ar and three cosmogenic ³⁶Cl age determinations provide new insight into the late Quaternary eruptive history of Erebus volcano. Anorthoclase from 3 lava flows on the caldera rim have ⁴⁰Ar/³⁹Ar ages of 23 ± 12, 81 ± 3 and 172 ± 10 ka (all uncertainties 2σ). The ages confirm the presence of a second, younger, superimposed caldera near the southwestern margin of the summit plateau and show that eruptive activity has occurred in the summit region for 77 ± 13 ka longer than previously thought. Trachyte from “Ice Station” on the eastern flank is 159 ± 2 ka, similar in age to those at Bomb Peak and Aurora Cliffs. The widespread occurrences of trachyte on the eastern flank of Erebus suggest a major previously unrecognized episode of trachytic volcanism. The trachyte lavas are chemically and isotopically distinct from alkaline lavas erupted contemporaneously in the summit region < 5 km away.     

Eruption episodes and magma recharge events in andesitic systems: Mt Taranaki,New Zealand

Eruption episodes, where a series of eruption events are generically related, can include the eruption of a wide spectrum of volcanic activity over decadal periods. This paper concentrates on the opening phases of an eruption episode which occurred approximately 1800 yrs BP from Mt Taranaki, New Zealand. These events spanned the eruption of differing bulk compositions and styles from two distinct vent locations; an andesitic sub-plinian eruption from the summit vent and a scoria cone-building eruption of basaltic magma from a satellite vent. Compositional profiles and zoning textures of plagioclase, amphibole and clinopyroxene phenocrysts from the opening andesitic event show evidence of magma mixing and subsequent crystallisation just prior to the initiation of the eruption episode. Titanomagnetite grain morphology and Ti variation suggest that the magma mixing event occurred within a few days to weeks before the eruption acting as a trigger for it. We present a magmatic model which is constrained by the petrological observations and eruptions of the episode. In this model magma differentiation at depth causes its rise and recharging of a mid-crustal magma storage area at 5–7 km. Although the recharging magma differed slightly in oxygen fugacity and temperature, it was compositionally and physically similar enough to the residing andesitic magma to allow efficient mixing. The petrological characteristics described here can be readily observed and enable identification of mixing events in other recent eruption episodes.     

Hydrogen emissions from Erebus volcano, Antarctica

The continuous measurement of molecular hydrogen (H₂) emissions from passively degassing volcanoes has recently been made possible using a new generation of low-cost electrochemical sensors. We have used such sensors to measure H₂, along with SO₂, H₂O and CO₂, in the gas and aerosol plume emitted from the phonolite lava lake at Erebus volcano, Antarctica. The measurements were made at the crater rim between December 2010 and January 2011. Combined with measurements of the long-term SO₂ emission rate for Erebus, they indicate a characteristic H₂ flux of 0.03?kg s^–1 (2.8?Mg? day^–1). The observed H₂ content in the plume is consistent with previous estimates of redox conditions in the lava lake inferred from mineral compositions and the observed CO₂/CO ratio in the gas plume (～0.9 log units below the quartz–fayalite–magnetite buffer). These measurements suggest that H₂ does not combust at the surface of the lake, and that H₂ is kinetically inert in the gas/aerosol plume, retaining the signature of the high-temperature chemical equilibrium reached in the lava lake. We also observe a cyclical variation in the H₂/SO₂ ratio with a period of ～10?min. These cycles correspond to oscillatory patterns of surface motion of the lava lake that have been interpreted as signs of a pulsatory magma supply at the top of the magmatic conduit.     

Investigating potential future changes in surface water flooding hazard and impact

Surface water flooding (SWF) is a recurrent hazard that affects lives and livelihoods. Climate change is projected to change the frequency of extreme rainfall events that can lead to SWF. Increasingly, data from Regional Climate Models (RCMs) are being used to investigate the potential water-related impacts of climate change; such assessments often focus on broad-scale fluvial flooding and the use of coarse resolution (>12 km) RCMs. However, high-resolution (<4 km) convection-permitting RCMs are now becoming available that allow impact assessments of more localised SWF to be made. At the same time, there has been an increasing demand for more robust and timely real-time forecast and alert information on SWF. In the UK, a real-time SWF Hazard Impact Model framework has been developed. The system uses 1-km gridded surface runoff estimates from a hydrological model to simulate the SWF hazard. These are linked to detailed inundation model outputs through an Impact Library to assess impacts on property, people, transport, and infrastructure for four severity levels. Here, a set of high-resolution (1.5 km and 12 km) RCM data has been used as input to a grid-based hydrological model over southern Britain to simulate Current (1996–2009) and Future (~2100s; RCP8.5) surface runoff. Counts of threshold-exceedance for surface runoff and precipitation (at 1-, 3- and 6-hr durations) are analysed. Results show that the percentage increases in surface runoff extremes, are less than those of precipitation extremes. The higher-resolution RCM simulates the largest percentage increases, which occur in winter, and the winter exceedance counts are greater than summer exceedance counts. For property impacts, the largest percentage increases are also in winter; however, it is the 12-km RCM output that leads to the largest percentage increase in impacts. The added-value of high-resolution climate model data for hydrological modelling is from capturing the more intense convective storms in surface runoff estimates.     

The variation and controls of mean transit times in Australian headwater catchments

Determining mean transit times in headwater catchments is critical for understanding catchment functioning and understanding their responses to changes in landuse or climate. Determining whether mean transit times (MTTs) correlate with drainage density, slope angle, area, or land cover permits a better understanding of the controls on water flow through catchments and allows first-order predictions of MTTs in other catchments to be made. This study assesses whether there are identifiable controls on MTTs determined using ³H in headwater catchments of southeast Australia. Despite MTTs at baseflow varying from a few years to >100 years, it was difficult to predict MTTs using single or groups of readily-measured catchment attributes. The lack of readily-identifiable correlations hampers the prediction of MTTs in adjacent catchments even where these have similar geology, land use, and topography. The long MTTs of the Australian headwater catchments are probably in part due to the catchments having high storage volumes in deeply-weathered regolith, combined with low recharge rates due to high evapotranspiration. However, the difficulty in estimating storage volumes at the catchment scale hampers the use of this parameter to estimate MTTs. The runoff coefficient (the fraction of rainfall exported via the stream) is probably also controlled by evapotranspiration and recharge rates. Correlations between the runoff coefficient and MTTs in individual catchments allow predictions of MTTs in nearby catchments to be made. MTTs are shorter in high rainfall periods as the catchments wet up and shallow water stores are mobilized. Despite the contribution of younger water, the major ion geochemistry in individual catchments commonly does not correlate with MTTs, probably reflecting heterogeneous reactions and varying degrees of evapotranspiration. Documenting MTTs in catchments with high storage volumes and/or low recharge rates elsewhere is important for understanding MTTs in diverse environments.     

Century‐long records reveal slight,ecoregion‐localized changes in Athabasca River flows

Reports of abruptly declining flows of Canada's Athabasca River have prompted concern because this large, free‐flowing river could be representative for northern North America, provides water for the massive Athabasca oil‐sands projects and flows to the extensive and biodiverse Peace–Athabasca, Slave and Mackenzie River deltas. To investigate historic hydrology along the river and its major tributaries, we expanded the time series with interpolations for short data gaps; calculations of annual discharges from early, summer‐only records; and by splicing records across sequential hydrometric gauges. These produced composite, century‐long records (1913–2011) and trend detection with linear Pearson correlation provided similar outcomes to nonparametric Kendall τ‐b tests. These revealed that the mountain and foothills reaches displayed slight increases in winter discharges versus larger declines in summer discharges and consequently declining annual flows (~0.16% per year at Hinton; p < 0.01). Conversely, with contrasting boreal contributions, the Athabasca River at Athabasca displayed no overall trend in monthly or annual flows, but there was correspondence with the Pacific Decadal Oscillation that contributed to a temporary flow decline from 1970 to 2000. These findings from century‐long records contrast with interpretations from numerous shorter‐term studies and emphasize the need for sufficient time series for hydrologic trend analyses. For Northern Hemisphere rivers, the study interval should be at least 80 years to span two Pacific Decadal Oscillation cycles and dampen the influence from phase transitions. Most prior trend analyses considered only a few decades, and this weakens interpretations of the hydrologic consequences of climate change. Copyright © 2014 John Wiley & Sons, Ltd.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

An experimental investigation of volatile exsolution in evolving magma chambers

Previous laboratory experiments investigating the fluid dynamics of replenished magma chambers have been extended to model effects resulting from the release of gas. Turbulent transfer of heat between a layer of dense, hot and volatile-rich mafic magma overlying cooler more evolved magma can lead to crystallization and exsolution of volatiles in the lower layer. Small gas bubbles can cause the bulk density to decrease to that of the upper layer and thus produce sudden overturning and initiate mixing, followed by further exsolution of gas and explosive eruption. These processes have been modelled in the laboratory using a chemical reaction between sodium or potassium carbonate and nitric acid to release small bubbles of CO₂. We have investigated both the initial overturning produced by gas release in the lower layer, and the subsequent evolution of gas due to intimate mixing of the two layers. The latter experiments, in which the reactants remained isolated in the two layers until overturning occurred, demonstrated unambiguously that the fluxes of chemical components across the interfaces between convecting layers are very slow compared to the flux of heat. This shows that the evolution of layers of magma of different origins and composition can take place nearly independently of each other. The magmas can coexist in the same stratified chamber, until their bulk densities become equal and they mix together. The processes illustrated in these experiments could occur in H₂O-bearing magmas such as in the calcalkaline association and in CO₂-bearing mafic magmas such as in silica undersaturated suites.