Evolution,transfer and release of magmas and volcanic gases: An introduction

In this introduction, recent progress in the study of the evolution, transfer and release of magmas and volcanic gases is briefly reviewed, based on discussions at the MAG2007 international conference held in Taipei, Taiwan, in April 2007. The meeting pooled the diverse expertise of igneous and volcanic gas geochemists. This special issue, which presents six case studies covering a range of topics from variations in mantle source compositions to the nature of volcanic degassing, reflects that diversity.     

Forced imbibition into a limestone: measuring P‐wave velocity and water saturation dependence on injection rate

Quantitative interpretation of time‐lapse seismic data requires knowledge of the relationship between elastic wave velocities and fluid saturation. This relationship is not unique but depends on the spatial distribution of the fluid in the pore‐space of the rock. In turn, the fluid distribution depends on the injection rate. To study this dependency, forced imbibition experiments with variable injection rates have been performed on an air‐dry limestone sample. Water was injected into a cylindrical sample and was monitored by X‐Ray Computed Tomography and ultrasonic time‐of‐flight measurements across the sample. The measurements show that the P‐wave velocity decreases well before the saturation front approaches the ultrasonic raypath. This decrease is followed by an increase as the saturation front crosses the raypath. The observed patterns of the acoustic response and water saturation as functions of the injection rate are consistent with previous observations on sandstone. The results confirm that the injection rate has significant influence on fluid distribution and the corresponding acoustic response. The complexity of the acoustic response —‐ that is not monotonic with changes in saturation, and which at the same saturation varies between hydrostatic conditions and states of dynamic fluid flow – may have implications for the interpretation of time‐lapse seismic responses.     

Seismic attenuation and velocity dispersion in fractured rocks: The role played by fracture contact areas

The presence of fractures in fluid‐saturated porous rocks is usually associated with strong seismic P‐wave attenuation and velocity dispersion. This energy dissipation can be caused by oscillatory wave‐induced fluid pressure diffusion between the fractures and the host rock, an intrinsic attenuation mechanism generally referred to as wave‐induced fluid flow. Geological observations suggest that fracture surfaces are highly irregular at the millimetre and sub‐millimetre scale, which finds its expression in geometrical and mechanical complexities of the contact area between the fracture faces. It is well known that contact areas strongly affect the overall mechanical fracture properties. However, existing models for seismic attenuation and velocity dispersion in fractured rocks neglect this complexity. In this work, we explore the effects of fracture contact areas on seismic P‐wave attenuation and velocity dispersion using oscillatory relaxation simulations based on quasi‐static poroelastic equations. We verify that the geometrical and mechanical details of fracture contact areas have a strong impact on seismic signatures. In addition, our numerical approach allows us to quantify the vertical solid displacement jump across fractures, the key quantity in the linear slip theory. We find that the displacement jump is strongly affected by the geometrical details of the fracture contact area and, due to the oscillatory fluid pressure diffusion process, is complex‐valued and frequency‐dependent. By using laboratory measurements of stress‐induced changes in the fracture contact area, we relate seismic attenuation and dispersion to the effective stress. The corresponding results do indeed indicate that seismic attenuation and phase velocity may constitute useful attributes to constrain the effective stress. Alternatively, knowledge of the effective stress may help to identify the regions in which wave induced fluid flow is expected to be the dominant attenuation mechanism.     

Gully catchments as a sediment sink,not just a source: Results from a long‐term (~12 500 year) sediment budget

Sediment delivery from hillslopes to trunk streams represents a significant pathway of mass transfer in the landscape, with a large fraction facilitated by gully systems. The internal gully geomorphic dynamics represent a considerable gap in many landscape and empirical erosion models, therefore a better understanding of these processes over longer timescales (10–10⁴ years) is needed. This study analyses the sediment mass balance and storage dynamics within a headwater gully catchment in central Europe over the last ~12 500 years. Human induced erosion resulted in hillslope erosion rates ~2.3 times higher than under naturally de‐vegetated conditions (during the Younger Dryas), however the total sediment inputs to the gully system (and therefore gully aggradation), were similar. Net gully storage has consistently increased to become the second largest term in the sediment budget after hillslope erosion (storage is ~45% and ~73% of inputs during two separate erosion and aggradation cycles). In terms of the depletion of gully sediment storage, the sediment mass balance shows that export beyond the gully fan was not significant until the last ~500 years, due to reduced gully fan accommodation space. The significance of storage effects on the gully sediment mass balance, particularly the export terms, means that it would be difficult to determine the influences of human impact and/or climatic changes from floodplain or lake sedimentary archives alone and that the sediment budgets of the headwater catchments from which they drain are more likely to provide these mechanistic links. Copyright © 2015 John Wiley & Sons, Ltd.     

Retreating ice: research in pro‐glacial areas matters

Glacier forefields are landscapes in transition from glacial to non‐glacial conditions; this implies intense geomorphic, hydrological and ecological dynamics with important on‐ and off‐site effects. This special issue collects 13 papers covering recent research in both (sub‐)polar and alpine pro‐glacial environments that focus on (i) pro‐glacial sediment sources, (ii) pro‐glacial rivers, (iii) pro‐glacial lakes, (iv) ground water and ice, and (v) the development of soil and vegetation in its interplay with morphodynamics. Advances in mapping, surveying and geophysical techniques form the basis for research perspectives related to the historical evolution of pro‐glacial areas, the understanding of complex interactions of multiple processes, and the effects of continued glacier recession. Copyright © 2015 John Wiley & Sons, Ltd.     

In situ pore-pressure evolution during dynamic CPT measurements in soft sediments of the western Baltic Sea

We present in situ strength and pore-pressure measurements from 57 dynamic cone penetration tests in sediments of Mecklenburg (n?=?51), Eckernförde (n?=?2) and Gelting (n?=?4) bays, western Baltic Sea, characterised by thick mud layers and partially free microbial gas resulting from the degradation of organic material. In Mecklenburg and Eckernförde bays, sediment sampling by nine gravity cores served sedimentological characterisation, analyses of geotechnical properties, and laboratory shear tests. At selected localities, high-resolution echo-sounder profiles were acquired. Our aim was to deploy a dynamic cone penetrometer (CPT) to infer sediment shear strength and cohesion of the sea bottom as a function of fluid saturation. The results show very variable changes in pore pressure and sediment strength during the CPT deployments. The majority of the CPT measurements (n?=?54) show initially negative pore-pressure values during penetration, and a delayed response towards positive pressures thereafter. This so-called type B pore-pressure signal was recorded in all three bays, and is typically found in soft muds with high water contents and undrained shear strengths of 1.6–6.4 kPa. The type B signal is further affected by displacement of sediment and fluid upon penetration of the lance, skin effects during dynamic profiling, enhanced consolidation and strength of individual horizons, the presence of free gas, and a dilatory response of the sediment. In Mecklenburg Bay, the remaining small number of CPT measurements (n?=?3) show a well-defined peak in both pore pressure and cone resistance during penetration, i.e. an initial marked increase which is followed by exponential pore-pressure decay during dissipation. This so-called type A pore-pressure signal is associated with normally consolidated mud, with indurated clay layers showing significantly higher undrained shear strength (up to 19 kPa). In Eckernförde and Gelting bays pore-pressure response type B is exclusively found, while in Mecklenburg Bay types A and B were detected. Despite the striking similarities in incremental density increase and shear strength behaviour with depth, gas occurrence and subtle variations in the coarse-grained fraction cause distinct pore-pressure curves. Gaseous muds interbedded with silty and sandy layers are most common in the three bays, and the potential effect of free gas (i.e. undersaturated pore space) on in situ strength has to be explored further.     

A web‐based visualization tool for exploring stakeholder conflicts in land‐use planning

This article describes an open source web‐based visualization tool for exploring stakeholder conflicts in land‐use planning. It implements a multi‐criteria, decision‐analytic framework which solves a conflict‐constrained knapsack problem in order to find Pareto efficient combinations of actions which maximize public value, given the conflict constraint. The user interface is centered around an interactive map in which geographical subareas are color‐shaded in proportion to the mean strength of preference for a particular chosen action within the population in the selected area. The aim of this work is to investigate the use of web technologies and formal problem‐solving methods for improving public participation in land‐use planning and support local governments in the pursuit of improved decision‐making. The tool was developed using the free statistical programming language R and the web application framework Shiny. Its usefulness is demonstrated through a case study in Upplands Väsby municipality in Stockholm County, Sweden. The article ends with a summary together with some ideas for future research and development.     

Walker circulation controls ENSO atmospheric feedbacks in uncoupled and coupled climate model simulations

Many climate models strongly underestimate the two most important atmospheric feedbacks operating in El Niño/Southern Oscillation (ENSO), the positive (amplifying) zonal surface wind feedback and negative (damping) surface-heat flux feedback (hereafter ENSO atmospheric feedbacks, EAF). This hampers a realistic representation of ENSO dynamics in these models. Here we show that the atmospheric components of climate models participating in the 5th phase of the Coupled Model Intercomparison Project (CMIP5) when forced by observed sea surface temperatures (SST), already underestimate EAF on average by 23%, but less than their coupled counterparts (on average by 54%). There is a pronounced tendency of atmosphere models to simulate stronger EAF, when they exhibit a stronger mean deep convection and enhanced cloud cover over the western equatorial Pacific (WEP), indicative of a stronger rising branch of the Pacific Walker Circulation (PWC). Further, differences in the mean deep convection over the WEP between the coupled and uncoupled models explain a large part of the differences in EAF, with the deep convection in the coupled models strongly depending on the equatorial Pacific SST bias. Experiments with a single atmosphere model support the relation between the equatorial Pacific atmospheric mean state, the SST bias and the EAF. An implemented cold SST bias in the observed SST forcing weakens deep convection and reduces cloud cover in the rising branch of the PWC, causing weaker EAF. A warm SST bias has the opposite effect. Our results elucidate how biases in the mean state of the PWC and equatorial SST hamper a realistic simulation of the EAF.