Open access peer-reviewed chapter

Studying Volcanic Plumbing Systems – Multidisciplinary Approaches to a Multifaceted Problem

Written By

Steffi Burchardt and Olivier Galland

Reviewed: 26 April 2016 Published: 21 September 2016

DOI: 10.5772/63959

From the Edited Volume

Updates in Volcanology - From Volcano Modelling to Volcano Geology

Edited by Karoly Nemeth

Chapter metrics overview

2,725 Chapter Downloads

View Full Metrics


Magma transport and storage beneath active volcanoes occurs in the so-called volcanic plumbing system (VPS), a network of different magmatic sheet intrusions and magma reservoirs. The complex physical and chemical processes, which occur in the volcanic plumbing system, are key parameters that control the occurrence of an eruption, as well as type and size of the eruption. It is therefore imperative to assess plumbing system processes and their dynamics. Traditionally, plumbing system research is done as a part of various scientific disciplines, each with its own research questions, methods, and terms. As a consequence, there is often little overlap and communication between the disciplines. In this chapter, we give an overview of the history of plumbing system research and outline the state of the art of the main scientific disciplines involved. We summarise the potential and limitations of each discipline and then discuss three key components to foster multidisciplinary research—namely communication, information, and education—which are essential to promote a better understanding of the complexity of volcanic plumbing systems.


  • volcanic plumbing system
  • magma transport
  • magma storage
  • multidisciplinary research
  • volcanology

1. Introduction

Volcanic plumbing systems (VPS) form a plexus of magma channels and reservoirs that are governed by complex interactions of chemical and mechanical processes that control how magmas are emplaced and how they propagate through the Earth’s crust to an eventual eruption (Figure 1; [1]). Volcanic plumbing systems thus set the stage for volcanic eruptions and govern the style and magnitude of eruptive activity including dramatic volcano–tectonic phenomena, such as caldera and sector collapses.

Figure 1.

Schematic sketch of different components of volcanic plumbing systems, highlighting the complexity of the different types of magma channels and reservoirs. The inset illustrates the processes that may occur within magma bodies.

Traditionally, however, the study of the plumbing system components, such as dykes, sills, and larger magma bodies (Figure 1), as well as their dynamics is strongly method-based, for example, focussing exclusively on the composition of plutonic bodies or the seismicity of magma ascent. To date, relatively few bridges between the distinct disciplines exist. In this chapter, we will give a short overview of the historical development of the main concepts on volcanic plumbing systems and the diversification of research disciplines that study plumbing systems. We then proceed to outline the aims and approaches, as well as the potential and limitations, of the main disciplines, namely field and structural analysis, igneous petrology and geochemistry, geophysics, geodesy, and modelling. We will then discuss challenges and opportunities of combining approaches from different disciplines to overcome some of the limitations and contribute to a better understanding of the complex and dynamic processes within active volcanoes.


2. History of the study of volcanic plumbing systems

While naturalists of the eighteenth century debated the origin of basalt, James Hutton was the first to propose that dykes and sills formed by the solidification of “subterraneous lava” and that even granite was once in a molten state [2]. Hutton’s views were based on meticulous field observations of magmatic intrusions exposed at the Earth’s surface. Performing one of the first systematic petrological experiments, Sir James Hall [3] delivered supporting evidence by melting and cooling volcanic glass and basalt to produce a variety of crystallinity in the experimental products.

By the time the first petrological classifications of igneous rocks had been established and microscopy and experimental petrology had started to emerge as distinct approaches in the second half of the nineteenth century, field mapping had become more focused on the shapes and structures of igneous intrusions [4]. Studies of igneous rock bodies from that time added laccoliths [5], batholiths [6], and other intrusion geometries to the previously described dykes, veins, and sills. Based on his observation of deformation in the rocks surrounding igneous intrusions in the Henry Mountains in Utah, Gilbert [5] started to study the mechanics of magma emplacement to unravel why magma at some occasions gets trapped in the crust or erupts at other occasions. In the following decades, researchers such as Daly [7] found evidence for different types of magma emplacement mechanisms. These mechanisms were gradually linked to specific types of intrusive bodies, a classification of which was established in the early twentieth century [8]. By this time, field observations in the eroded volcanic complexes of the British Isles had identified ring dykes and cone sheets as specific intrusion types [9, 10].

Petrologists of the early twentieth century focused in turn on fractional crystallisation and magma differentiation in sills and layered intrusions [1114]. Hence, from the first half of the twentieth century, field-based research concerning volcanic plumbing systems was conducted in parallel as part of the disciplines of igneous petrology, structural geology, and volcanic geology, as well as by distinct communities studying granitic plutons, fossil subvolcanic complexes, and extrusive rock suites in volcanic areas.

The second half of the twentieth century saw major advances in the understanding of volcanic deposits and eruptions, as well as how these are linked to the underlying volcanic plumbing system. One of the most influential researchers of this time is George P. L. Walker who, together with his students, and based on systematic observations and quantitative mapping of extinct volcanic complexes in the British Isles and Eastern Iceland, fundamentally advanced our understanding of magma transport in dykes and the architecture of shallow magma plumbing systems in central volcanoes [15, 16].

Apart from the subdiscipline of experimental petrology, even other types of modelling have been employed to study volcanic plumbing systems. Analogue modelling (or laboratory modelling) of intrusive and volcanic processes became an increasingly important method (see Section 3.5), particularly after Hubbert [17] and later Ramberg [18, 19] had introduced the principle of scaling to ensure similarity between geological and model systems [20]. Numerical modelling of plumbing systems, on the other hand, is rooted in the analytical solutions of mostly fluid-dynamic and stress-field problems within classical structural geology. Development of computer software that is often originally designed for materials science simulations and the ever increasing computational capacity have led to a boost in numerical modelling studies of—among many other Earth science topics—volcanic plumbing systems.

The second half of the twentieth century was also the beginning for geophysical studies of volcanic plumbing systems. While gravimetric methods are nowadays routinely applied to study granitic plutons and subvolcanic complexes [21, 22], geophysical monitoring techniques have become a standard in active volcanic areas (see Section 3.3). In addition, since the launch of satellites and satellite-based mapping of the Earth’s surface in the 1960s, geodesy has become an increasingly popular method to infer processes related to magma transport in active volcanoes (see Section 3.4).

The history of plumbing system research reflects the methodological and conceptual development and diversification of the Earth sciences in general. Since the days of the early naturalists who described natural phenomena as a whole, our understanding of volcanic plumbing systems has deepened considerably. However, at the same time, as we started to dig deeper into specific aspects, we mostly lost our view of the big picture. Currently, volcanic plumbing system research is carried out in the disciplines of igneous petrology, structural geology, volcanology, geophysics, and geodesy, employing methods such as field mapping, major and trace element analysis, seismology, analogue modelling, and interferogrametry of radar measurements from satellites. Section 3 gives an overview of the state of the art of the most important approaches.


3. Methods commonly used to study volcanic plumbing systems

3.1. Field geology and structural analysis of volcanic plumbing systems

Field based studies of volcanic plumbing systems have initially aimed at a qualitative description of intrusive phenomena and focused on a classification of intrusion lithologies and morphologies to understand how magma is transported, stored, and evolves in the crust [23]. With time, field work has become increasingly quantitative, producing a detailed record of the compositions and structures associated with magmatic intrusions [24]. Analyses of the composition, absolute and relative ages, and dimensions of the components of magmatic plumbing systems have produced a more and more systematic view of the emplacement and evolution of plumbing systems [25, 26]. Besides the recording of variations in lithology and emplacement-related structures, mapping of magmatic intrusions often includes the study of magmatic fabrics recorded within the igneous rocks [27], such as preferred orientation of phenocrysts and of magnetic minerals using their so-called anisotropy of magnetic susceptibility (AMS; [28, 30]). Recently, classical methods, such as field mapping with paper maps and compass, have become complemented by modern digital mapping techniques using global positioning system (GPS) and smart phones (Figure 2). At the same time, the possibilities to analyse structural data collected in the field become more and more sophisticated. In addition to thorough statistical analysis and stereographic projection [31, 32], three-dimensional (3D) structural modelling is used to visualise, reconstruct, and interpret structural field data [3335].

Figure 2.

Examples of field and structural studies of volcanic plumbing systems. (a) Geological map of the Isle of Rum draped over a Digital Elevation Model. Map modified from Scottish Natural Heritage 1:20,000 map; © SNH. (b) High-precision GPS mapping of the roof of a granitic pluton (cf. [40, 41]). (c) Result of 3D reconstruction of a granitic pluton based on GPS mapping. (d) Measuring fractures in a granitic intrusion using a smart phone and Field Move Clino app by Midland Valley Ltd. (e) Outcrop of basaltic inclined sheets and corresponding projection of poles to orientation planes and density of 321 such sheet intrusions in an equal area, lower hemisphere plot (cf. [42]). (f) Probability Density Function (PDF) of the thickness of magmatic sheets, such as shown in (e) compared to a selection of statistical distributions. (g) 3D structural model of the plumbing system of the Ardnamurchan central complex, Scotland, produced using Move by Midland Valley Ltd (modified from [34]).

Igneous plumbing systems in outcrops regularly extend over several kilometres. Such wide extents represent a substantial challenge for (1) having a correct overview understanding of the structure of the exposed VPS and (2) completing structural field surveying in a manageable time. In addition, numerous field areas of interest are hardly reachable for direct observations, such as very steep mountains and crevasses. Recent technologies, including LiDAR scanners, drones, and robots, have recently started overcoming these challenges. The one hand, LiDAR scanners have been used to produce high-precision and high-resolution textured virtual outcrop models [36]. Such digital models allow detailed, quantitative 2D and 3D fracture mapping of extensive, mostly subvertical outcrops [37]. On the other hand, drone surveys combined with photogrammetric tools are very helpful to produce virtual outcrop models and orthorectified images of extensive subhorizontal or gently dipping outcrops, and provide new observational perspectives for extensive structural mapping [38]. Finally, robots equipped with monitoring tools can explore the Earth's interior that is inaccessible for humans [39]. The data produced by these modern tools allow new possibilities for post-field digital mapping of extensive areas, thus shortening the field campaigns.

Therefore, field studies of exhumed plumbing systems form the foundation of our conceptual understanding of the individual components of volcanic plumbing systems, their morphologies, sizes, and emplacement mechanisms, as well as of characteristic structures in the host rock associated with the evolution of magmatic intrusions. Field work in eroded volcanic areas is therefore of critical importance, because it provides us with fundamental information that may be used to benchmark numerical and laboratory models (see Section 3.5) and to interpret geophysical and geodetic data (see Sections 3.3 and 3.4).

The main limitation of the field-based approach to study volcanic plumbing systems is that fossil and eroded volcanic plumbing systems represent a snapshot of the final state of the magmatic system only, while active volcanoes do not permit a detailed look inside their VPS. An outcrop in an eroded volcano is essentially a snapshot of the sum of all superimposed processes a suite of rock has experienced. It is thus not always straightforward to extract accurate age relationship between individual units or even distinct boundaries at times, as well as to deduce what dynamic processes were contributing to the final picture.

3.2. Petrological and geochemical studies of volcanic plumbing systems

Igneous petrology and geochemistry are among the classic approaches used to characterise volcanic plumbing systems (see also Section 2) and aim to describe the conditions, time scales, and characteristics of the chemical evolution of magma. This characterisation is generally based on the mineralogy and textures of igneous rocks, as well as their major and trace element composition (e.g. Figure 3). In order to study the minerals, textures, and compositions of rock samples collected in the field, igneous petrologists employ a wide range of analytical techniques, which have been developed simultaneously with, and strongly facilitated, an increasing understanding of the chemical evolution of minerals and melts. These analytical techniques have emerged in a rapid succession and have become increasingly precise since igneous petrology was established as a discipline and the first microscopes were built in the second half of the nineteenth century (see Section 2). Today, rock textures and the type and associations of minerals in thin sections of rock samples are studied using a variety of microscopes. The bulk, or whole-rock, composition of an igneous rock that mostly comprises the quantity of major elements or their oxides can be derived through, for example, exposing a powdered sample of a rock to X-rays (X-ray fluorescence (XRF)). Other analytical techniques, such as mass spectrometry, can quantify the isotopic composition of the major and trace elements of a sample.

Data on the mineralogy and composition of igneous rocks are probes into the chemical evolution of volcanic plumbing systems and can be used in many ways. In many cases, the whole-rock composition is characteristic of the geodynamic setting and origin of a magma and can also be used to discriminate processes in the plumbing system, such as fractional crystallisation (e.g. [43, 44]). The concentration of trace elements can be used in geochemical modelling to quantify processes in the volcanic plumbing system, such as assimilation of country rocks into the magma (e.g. [45]). The decay of radiogenic isotopes is used as a standard tool to determine the absolute age of a rock sample, which has led to detailed insights into the time scales of magma emplacement in volcanic plumbing systems (e.g. [46]).

At the scale of individual crystals, crystal size distributions, chemical zoning, and textures can serve as records of the chemical and thermal evolution of their host magmas (e.g. [4749]). Crystal growth rates and chemical diffusion across crystal zones can furthermore be used to quantify the time scales of, for example, magma storage or replenishment (e.g. [50]). Moreover, pressure and/or temperature dependent mineral compositions can be used as so-called geobarometers and/or geothermometers that reveal the depth and conditions at which certain minerals grew, which usually corresponds to magma reservoir depths (e.g. [51, 52]).

In order to quantify the chemical evolution of minerals and magmas in situ, petrological experiments are used to simulate, for example, the influence of pressure and/or temperature on mineral compositions [53] and the reaction of magma with crustal rocks [54, 55].

Petrological and geochemical studies of volcanic plumbing systems are thus the foundation of our understanding of the chemical and thermal processes during magma storage and offer insight into the time scales of magma transport and evolution. Estimates of the depth of magma storage provide valuable constraints on the interpretation of geodetic and geophysical monitoring data of active volcanoes (see Sections 3.3 and 3.4). Furthermore, a characterisation of the types of magma erupted from a volcano and an understanding of the processes of magma evolution in the plumbing system allow evaluating the probable type of eruption in the future.

The main limitation of the petrological and geochemical approaches is that magma samples from active volcanic plumbing systems are generally not available. Although many of the processes of magma evolution can be constrained based on petrological experiments, the major and trace element compositions of rocks and minerals can in many cases not be attributed to any specific process. The contribution of each individual process is often difficult to quantify. Moreover, the insights derived from petrological and geochemical approaches apply to the geological past of the plumbing system only. It is therefore not always straightforward to conclude on the present state.

Figure 3.

Examples of petrological and geochemical studies of volcanic plumbing systems. (a) High-precision chemical analysis and element mapping using a field emission gun electron probe microanalyser. (b) Thin section of a zoned plagioclase crystal (zonation indicated by dashed line) in basalt of the 2015 Holhraun eruption. Crossed polars. Crystal ca. 1 mm across. Crystal textures can be used to reconstruct processes in the volcanic plumbing system [57]. (c) Scanning electron microscope (SEM) image of vesiculated xenolith erupted offshore El Hierro [56]. (d) False-colour SEM image of a zoned clinopyroxene crystal in lava erupted at Holuhraun, Iceland. Zones and inclusions can be used to reconstruct processes in the volcanic plumbing system [58]. (e) Results of geobarometric modelling of plagioclase and clinopyroxene in tephras from Katla volcano, Iceland [58]. (f) Mineralogy of pyroxene crystals in Katla tephras [58]. (g) Analysing trace-element contents in igneous rocks. Image courtesy of Christophe Galerne. (h) Zr versus Nb plot of igneous rocks and xenoliths from the Canary Islands, Spain. Trace element compositions can be used to understand the origin of, and relationships between, rock groups [56].

3.3. Geophysical studies of volcanic plumbing systems

Geophysical studies of volcanic plumbing systems employ a variety of methods (Figure 4) that detect and quantify either the physical properties of different geomaterials, such as magma versus solid rock or igneous versus sedimentary rocks, or the effects of active physical processes, such as seismicity caused by the movement of magma through the crust. Geophysical methods are therefore applied to study both active volcanoes and extinct subsurface or eroded plumbing systems.

Figure 4.

Examples of geophysical studies of volcanic plumbing systems. (a) Seismicity related to the emplacement of a dyke preceding caldera collapse in Miyakejima volcano, Japan, in 2010 [62]; see references therein for source data). (b) and (c) Linear gravity anomalies (marked by black dots) interpreted to be related to dykes on the Moon [59]. (d) Above: Vertical section of a 3D seismic dataset showing the signal related to a sill. Below: 3D seismic image of the same saucer-shaped sill [63]. (e) Volcanic plumbing system of the Toba caldera inferred from seismic tomography [64].

In order to characterise physical properties of different materials in the crust, a wide range of geophysical methods, such as gravimetry, magnetometry, and electrics, are commonly applied to map volcanic plumbing systems and associated hydrothermal systems. These methods detect anomalies in the Earth’s gravimetric, magnetic, and electrical resistivity fields. The anomalies are caused by the presence of rock types with properties contrasting to those of the surrounding country rocks, for example, higher or lower density, resistivity etc. Depending on the type and size of target, for example, mapping a granitic intrusion, and the expected depth range of the target, appropriate methods, and acquisition techniques are chosen. Instruments may be mounted on airplanes and satellites or installed in arrays at the Earth’s surface. Acquired geophysical data can then be used to map the physical properties of the Earth’s crust. Such maps allow locating the extent and map the outline of magmatic intrusions, which is extremely useful for large-scale reconnaissance and mapping in poorly exposed or inaccessible areas, such as moons or other planets (e.g. [59]). The spatial resolution of the data mainly depends on the contrast in physical properties between the target and the country rock. The acquired data on geophysical anomalies can also be used in inversion models that infer the subsurface distribution of a particular physical property by reproducing the pattern and magnitude of the anomaly making simple assumptions about the geology, for example, homogeneity and isotropy of the involved rocks. For instance, the results of inversion models can be applied to determine the subsurface shape and volume of a granitic intrusion (e.g. [60, 61]). The main limitation of geophysical mapping at depth using inversion modelling is that the data are acquired at the Earth’s surface (more or less in 2D) is then used to interpret the 3D distribution of physical properties. Moreover, inversion modelling that often produces as very similar solutions in spite of different input parameters. Hence, the number of fitting solutions may be infinite, and thus, model interpretations are often non-unique.

Geophysical methods are also frequently applied to study active volcanoes in order to locate magma storage levels and to infer processes related to magma movement. For instance, microgravimetry monitoring can be applied to detect magma flow into the plumbing system by monitoring changes in the local gravimetric field, a method often combined with geodetic monitoring [65]. In the brittle crust, volcano-tectonic earthquakes that are interpreted to be related to magma movement are studied in the subdiscipline of volcano seismology [66]. Using a network of seismometers, seismicity in active volcanoes can be monitored, for example, to trace propagating sheet intrusions, which gives valuable insight into dyke emplacement mechanisms [67, 68]. The properties of seismic waves (frequencies, wave forms, etc.) can be interpreted to derive information about processes such as the rise of magma in a conduit or the intrusion of a cryptodome [60]. Travel times of seismic waves through the crust beneath a volcano can be used as input parameters for inversion models, a technique called seismic tomography. The seismic velocity structure of a volcano can be used to map subsurface areas with an increased percentage of melt, potentially corresponding to magma reservoirs [69, 51]. However, the detection limit of seismic tomography for magma bodies is on the order of several hundred metres, depending on the seismic array and acquisition conditions on the one hand and the structure of the individual volcano on the other. Therefore, most parts of the volcanic plumbing system, such as average sized dykes and sills, may be undetectable using seismic tomography.

During the last decade, seismic reflection and refraction have been extensively implemented to image volcanic plumbing systems in sedimentary basins [70, 71]. While all geophysical methods described earlier employ passive measurements of natural rock properties, seismic reflection and refraction methods represent active geophysical surveying: acoustic waves are sent into the subsurface and reflected and refracted at interfaces between rocks of contrasting impedance. Thus, the resulting seismic images are a direct echography of subsurface structures and not the result of data inversion. Initially designed for studying the structure of sedimentary basins and hydrocarbon exploration, seismic data proved essential to document the presence of voluminous sill and laccolith complexes in numerous basins worldwide [73, 74]. On seismic images, igneous intrusions are very prominent because of the strong impedance contrast between igneous and sedimentary rocks. The main advantages of seismic data are that (1) they can image the shapes of entire intrusions in three dimension [74, 75], which can potentially constrain magma flow directions [64, 76], and (2) it is possible to constrain the intrusion-scale deformation induced by magma emplacement [77, 78, 83] and intrusion-fault interactions [79]. The main limitations of seismic data are (1) the limited spatial resolution (about 20 to 40 m), (2) the limited possibility to image subvertical features, such as dykes, (3) seismic artefacts that produce interpretable features that do not exist, and (4) the availability of the seismic data, as they are often kept confidential by oil companies.

Geophysical methods are thus powerful tools to study physical properties and signals produced by processes of volcanic plumbing systems. Despite their limitations, these methods are often the only way to derive information about the location and properties of plumbing systems at depth.

3.4. Geodetic studies of volcanic plumbing systems

Volcano observatories monitor active volcanoes with a combination of techniques, including seismic networks, temperature, and gas monitoring, as well as geodetic techniques (; Figure 5). Geodetic monitoring of the deformation of the Earth’s surface relies mainly on ground-based tiltmeters, GPS networks, and satellite-based Interferometry of Synthetic Aperture Radar (InSAR). A network of GPS stations has to be deployed and maintained in the field and typically delivers daily data for each station. Continuous monitoring can provide surface deformation data with an accuracy of about 1 mm/yr, but the data quality strongly depends on the station network density [80]. On the other hand, InSAR records surface deformation in the line of sight of the satellite at 10–40 days intervals. So, InSAR can reach accuracies of less than 1 cm for a stack of interferograms. Apart from disturbances due to vegetation and glaciation, InSAR allows geodetic monitoring of volcanoes worldwide, irrespective of accessibility [80]. Thanks to InSAR monitoring data, we now know of more than 140 currently deforming volcanoes (

The corrected and processed surface deformation data are interpreted using geodetic volcano-deformation modelling, which reproduces the pattern and magnitude of the surface signal with analytical models based on fluid and solid mechanics [84]. The results of geodetic models have been used to interpret volcano deformation in many volcanoes worldwide [8587] and have led to a better understanding of magma movements. Even though surface deformation in volcanoes is not always related to an imminent eruption, geodetic monitoring, and modelling have become standard tools for eruption forecasting [8890].

Geodetic models are analytical or numerical solutions that comprise three components: A. a deformation source representing the VPS, B. a static process in the deformation source, and C. the model crust bounded by the free surface and characterised by a rheological law. The combination of A, B, and C produces deformation of the model surface, which is compared to the measured surface deformation. As the solution of geodetic models with highly different input parameters can produce very similar results [91, 92], the geodetic modeller is left with an infinite number of possible, non-unique solutions. Non-linear inversion is thus applied to select likely source parameters, to estimate uncertainties, and to compare the goodness-of-fit of different models [93, 94]. When comparing the fit of models, the simplest one, which reproduces the recorded data best, is usually preferred, even though it may be geologically implausible for the studied volcano.

Figure 5.

Examples of geodetic studies of volcanic plumbing systems. (a) GPS data (black arrows) and results of Mogi models (red arrows) of ground deformation of Santorini. The Mogi source (red dot) is given with confidence intervals in per cent [81]. (b) Post-eruptive deformation of Fernandina volcano, Galapagos islands, measured with InSAR [82]. (c) Structure of the volcanic plumbing system of Nyamulagira volcano, Congo, inferred from InSAR and stress modelling [83].

Since it is impossible to quantify the uncertainties of geodetic model results on active geological processes that occur in the subsurface, it is important to constrain the models by using model input parameters provided by depth and location of earthquakes related to moving magma [86], earthquake “shadows” around magma reservoirs (see Section 3.3; [69]), or by volcanic eruption volumes and volcanic plume heights [95]. However, these methods miss small- and medium-sized magma chambers, as well as the shape of the plumbing system. Furthermore, geobarometry (see Section 3.2) may indicate the depth and duration of long-term magma storage, but cannot resolve levels of short-time storage [51]. More importantly, however, the lack of collaborative work between different disciplines that study volcanoes has so far prevented a more realistic interpretation of surface deformation in volcanoes. So in practice, deformation source geometries are often arbitrarily chosen without geological validation.

3.5. Laboratory, theoretical, and numerical modelling of volcanic plumbing systems

Since active volcanic plumbing systems are located within the crust and therefore not accessible for direct observations and even the shallowest parts, that is, volcanic vents, are extremely hazardous to approach, our understanding of the dynamic processes at work in VPS is based on either direct observations of extinct plumbing systems (field geology; see Section 3.1) or indirect observations through geophysics, geodetic surveys, and geochemistry and petrology of volcanic products (see Sections 3.2 to 3.4). However, these approaches are insufficient to identify and quantify the physical parameters and laws that dominantly control the dynamics of VPS.

Figure 6.

Examples of modelling studies of volcanic plumbing systems. (a) Experimental set-up (left) and characteristic laboratory modelling results (right) of ground deformation induced by shallow magma emplacement [96]. Results are differential DEMs (top row) and divergence maps (bottom row) calculated from horizontal displacement maps at two times of an experiment; positive divergence means dilation (tensile cracks), whereas negative divergence means compression (reverse faulting). (b) Numerical results displaying the evolution of density cross sections (left column) and deviatoric stress cross sections (right column) calculated from 2D finite volume simulations of explosive venting [97]. (c) Set-up and boundary conditions of a theoretical model of flat-lying intrusion of radius a, emplaced at the base of an elastic plate of thickness h and weight q0 [98]. The magma has a heterogeneous pressure P(r).

To overcome these limitations, modelling techniques, such as laboratory, theoretical, and numerical modelling are employed to study processes in volcanic plumbing systems (Figure 6). By definition, modelling implies simulating processes in a model system to understand the dynamics of the same processes in nature. However, as the model system is a crude simplification of a part of nature, we cannot expect the model to represent the full complexity of natural processes. The assumptions and simplifications we make when setting up the model, will have a large impact on the model outcome. It is therefore important to keep in mind that models are models, and not nature, and that model results are only as good as the input data. Consequently, modelling of volcanic plumbing systems should always use input data from plumbing systems in nature, for example, field data, and should always be validated against observations in nature.

3.5.1. Laboratory modelling

The aim of laboratory modelling is to simulate processes in volcanic plumbing systems at manageable scales, such as laboratory lengths (sand boxes of a few metres length maximum) and time scales (minutes to days) (Figure 4ac; [20]). In the nineteenth and early twentieth centuries, laboratory models were phenomenological only and used as source of inspiration and/or proof of concepts to demonstrate the existence of a phenomenon [99]. The main limitations of these models were their applicability to geological systems. Laboratory models entered a new era after Hubbert [17] introduced the scaling theory that established scale relationships between laboratory-scale model and geological-scale system. However, these scale relationships address the similarity between the two systems only, that is, they cannot be used to understand the physics behind the modelled processes. Hence, the development of the dimensional analysis concept within the field of physics was a breakthrough for laboratory modelling. Dimensional analysis identifies fundamental physical laws, so-called scaling laws, that are simple relationships between dimensionless parameters, that is, their validity is scale independent ([100] and references therein). Combining dimensional analysis and similarity principles to laboratory models allow the modeller to unravel the fundamental physical laws and apply them to geological-scale volcanic plumbing systems [20].

Successful implementation of laboratory models requires the use of model materials with relevant properties. Classic rock analogues are gels [101, 102] and loose sand [103, 104] to simulate elastic rocks and cohesion-less Coulomb (i.e. frictional) rocks, respectively, which are two end-member behaviours of natural rocks. Recently, new materials of more complex rheology have been successfully used [105, 106]. Cohesive granular materials, some of which have variable cohesion, account for the complex elastoplastic properties of the brittle upper crust, and can simulate both mode I (tensile) and mode II (shear) fracturing [107]. Gels offer the possibility to study the viscoelastic properties of the Earth’s crust [108110]. On the other hand, many different fluids of diverse rheology and viscosity have been used as magma analogues [20]. Most noticeable is the implementation of two-phase fluids to model the dynamics of magma-bubble suspensions [111, 112]. These materials offer the possibility to model the natural complexity of geological systems, and hence, the physics of the complex processes in volcanic plumbing systems is studied.

Another recent enhancement of laboratory models has been the implementation of various monitoring methods that allow for more quantitative data acquisition and analysis [113]. The most commonly used methods are Particle Image Velocimetry (PIV) and Digital Image Correlation (DIC) [62, 114, 115] and stereo-photogrammetry [96] to monitor deformation of surfaces, as well as X-ray scanners to monitor in situ model interiors [116, 117]. The resulting quantitative data have become essential for (i) constraining the physical laws governing the modelled processes and (ii) integrating laboratory results with quantitative geological, geophysical and geodetic data.

The main limitations of laboratory models to be considered include the scale gap between the laboratory and geological systems, which is often critical and hard to fully constrain. The full characterization of model materials is not straightforward, even though it is essential, and requires a solid fundamental mechanics background. For example, cohesive and cohesion-less granular materials exhibit elastic properties that are difficult to measure, though they are crucial for the dynamics of the brittle crust [118, 119]. Another challenge is connected to recording the modelled processes, which often requires advanced techniques, such as laser or X-ray scanning, which are very costly and heavy to handle. Moreover, stress fields in granular material experiments cannot be measured, although stress distribution is a key factor on the dynamics and evolution of VPS. Finally, first-order assumptions need to be made due to technical limitations of model systems, neglecting for instance thermal effects [120] and chemical evolution of magmas.

3.5.2. Theoretical modelling

In order to assess first-order scaling laws governing a process in volcanic plumbing systems, theoretical (or analytical) models are employed. Theoretical models solve analytically or semianalytically the mathematical equations governing the studied processes (Figure 4d). During the mathematical derivations, the first-order scaling parameters commonly appear spontaneously [98, 121, 122], such that the first-order effect of the key physical parameters is obvious, and it is possible to identify whether distinct physical behaviours can be expected. Pioneering theoretical models of hydraulic fractures [123, 124], for example, are the foundation for most subsequent models of dyke or laccolith emplacement.

Although theoretical models are powerful, they require an advanced level of mathematical skill, which is often beyond the training level of most geoscientists. Moreover, theoretical models exhibit limitations, such that the solution of the equations requires a large number of simplifying assumptions. As a consequence, many geoscientists are not aware of the assumptions behind the mathematical models and extrapolate the model results beyond their domain of validity. Furthermore, geoscientists often accept the assumptions underlying the theoretical models as rules, such that the assumptions are rarely questioned, in spite of contradicting evidence from nature. For example, it is commonly assumed that dykes and sills are hydraulic fractures emplaced in purely elastic host rocks, although clear field, geophysical, and laboratory evidence of first-order inelastic deformation accommodating their emplacement are frequently observed.

3.5.3. Numerical modelling

Numerical modelling uses mathematical equations to simulate a simple process or a combination of several processes (Figure 4e). Using numerical codes and software, which have often been developed for materials science applications, the studied system is subdivided into subsystems, so-called elements, and the equations solved for each element. Hence, numerical models can overcome some of the main limitations of theoretical models by accounting for boundaries of complex shapes [125], complex heterogeneities, etc. The models can take into account material properties, temperatures and pressures and calculate stresses and strains, as well as changes in temperature, pressure, and material properties. Furthermore, numerical models can account for static [86, 126], quasistatic, and transient processes [127, 128]. Therefore, in a transient model, it is possible to calculate at each time step, the stress field [129], flow field [130], and/or temperature field [131]. This represents a major advantage with respect to laboratory models, in which the material properties and processes within the model are challenging to measure and monitor. Another advantage is that the scales of numerical models can be directly set as the scale of geological systems. As computational power is the main constraint on how complex a numerical model may be, it is easy to systematically vary model parameters, such as the material properties, a major advantage compared to laboratory modelling.

There are two main types of numerical models with fundamentally different approaches, continuum models and discrete models. Continuum models solve the fundamental equations of continuum mechanics, such as Stokes and Navier-Stokes equations for fluid flow, Hook’s law for elastic deformation, the heat equation calculating heat diffusion and/or advection, and Darcy’s law for porous flow. Various solving methods exist, such as finite difference [132], finite element (e.g. [128]), and finite volume [97]. In volcanic plumbing system research, continuum models are commonly used to model the elastic deformation induced by magmatic intrusions [42, 86, 133], magma flow within intrusions [125], thermal impacts of intrusions on their host rock [131, 134], thermal convection within magma reservoirs [130], and recently magma emplacement [128, 132]. The main limitation of continuum models is that they cannot simulate the formation of new discontinuities, such as fractures or dykes, which is why the geometry of discontinuities has to be prescribed [42].

Discrete models calculate the behaviour of a pack of particles that interact with each other. Some main methods are discrete element models (DEM) and Lattice Boltzman (LBM). The particle interactions are represented by bounds with elastic, friction, or fluid-like properties and strengths [135]. Discrete models are very valuable to simulate processes that create discontinuities, such as fractures, and they appear very suitable to simulate caldera collapse [127], dyke, sill, and laccolith emplacement [129, 136]. The main limitation of discrete models is to link the particle-scale interactions (e.g. bound stiffness, bound strength) with the bulk properties of the particle packing (Young’s modulus, cohesion, internal friction), requiring heavy systematic calibration before running the models and interpreting the results.


4. Challenge: integrating the methods

Since our training and level of experience fundamentally influence how we interpret geological problems [137], different studies of volcanic plumbing systems that employ methods of different disciplines often lead to contradictory conclusions. For instance, shallow magma storage beneath Katla volcano has been suggested based on seismic and geodetic data [138, 139], whereas petrological studies so far found evidence for lower- to mid-crustal storage only [140]. Such contrasting results may lead to fundamentally different interpretations when it comes to volcanic risk and hazard assessment and may therefore have unforeseeable consequences (cf. [58]). However, these contrasting results may be due to the limitations of individual methods that commonly assess a fraction of the complex magmatic dynamics only, a discrepancy that is hardly possible to assess for decision makers.

The disciplinary boundaries in plumbing system research in part reflect the historical development of the Earth sciences in general (see Section 2). During the twentieth century, each of the disciplines involved in the study of magma transport and storage has become more specialised, and new disciplines, such as volcano geodesy, have emerged. Methods have become more and more sophisticated, revealing the complexity of individual, and the interplay between, physical and chemical processes at scales ranging from the size of the crystal lattice to the thickness of the lithosphere. While specialisation and methodological progress naturally continue, the true challenge of the twenty-first century is to overcome the boundaries of our disciplines in order to truly assess the complexity of volcanic plumbing system dynamics.

There are three key components required to foster multidisciplinary research in general, and VPS research in particular: communication, information, and education. The methodological specialisation of each discipline has been naturally accompanied by the development and evolution of specific terms, and in most cases even jargon, that are used within individual disciplines, but might not be accessible to researchers outside of the particular discipline. For instance, while structural geologists may refer to a magma body of a certain size as “pluton”, the same intrusion may be a “magma reservoir” to igneous petrologists, a “shadow in s-wave attenuation” to geophysicists, and a “deformation source” to geodesists. On the other hand, different disciplines may use the same term but with different connotation. In order to overcome disciplinary boundaries, we need to be able to understand each other, which is preferentially achieved by developing a common language and avoiding jargon. In addition, it is key to learn to communicate one’s research to a broad audience, for example, by explaining the essential message of a certain diagram or equation in terms of the research question or by summarising research results in the end of a publication in a way that is generally understandable.

Field studies Petrology Geophysics Geodesy Modelling
Key questions  How is magma transported, stored & erupted? Is there an interaction with the host rock? How does magma form and evolve? When did the magma form? What triggers eruptions? Where does the signal come from? What does it mean? Is the signal magmatic? What does the deformation signal mean? Is the signal magmatic? How does the magmatic plumbing system work? Are there general laws we can apply?
Contribution from [below] towards [above]
Field studies Relative ages; structure; context
& correlation; mechanical
Lithology; types of structures; resolution Shape &
evolution of sources; complexity; response of rocks to deformation
Interpretation of kinematics; structural relationships; geometry; mechanics; properties
Petrology Absolute ages; origin of magma; source evolution; eruption types; chemostratigraphy Composition of samples; depths
& rates; thermal evolution; rock properties
Depths of
storage; time scales; sources
Time scales; composition; feedback
Geophysics Depths & sizes; 4D imaging at depth Rates; depths & sizes; 3D to 4D imaging at depth 3D to 4D
imaging at
depth; rock properties & geometries; rates
Rock properties; rates of processes
Geodesy 3D crustal structure; depths & sizes; rates Rates; depths & sizes; temperature; density Rates Response signal of rocks due to processes
Modelling Feedback on
mechanisms &
of structures
Conceptual context; quantify conditions of magma evolution Mechanisms & processes; state
of stress
Mechanisms & processes

Table 1.

Overview of the main disciplines involved in studying volcanic plumbing systems, their main foci and the contribution each one of them can make to the other disciplines.

Another requirement to overcome disciplinary boundaries is to inform researchers from other disciplines about the capabilities and limitations of our own approaches, a process that requires clear communication. At the same time, we need to learn how other approaches can be used to complement our own methods. Table 1 provides an overview of how approaches from each of the disciplines described earlier (Section 3) can complement research in other disciplines. The key questions and strengths of any of the disciplines can provide complementary data for research in other disciplines. For instance, quantitative constraints on the dimensions of sheet intrusions derived from field work in eroded volcanic systems [32, 141] can be used as input parameters for geodetic modelling of volcano deformation. Geobarometric calculations of the crystallisation depth of crystals in eruptive products of a volcano can be combined with seismic tomography to constrain the depth of magma storage [51]. 3D seismic imaging of intrusions in sedimentary basins can complement field studies on the shape and emplacement of sills [79, 129]. Geodetic monitoring data may help to quantify the amount of magma drained from the volcanic plumbing system during dyke intrusion [142]. Modelling of the emplacement of sheet intrusions can reveal the parameters that control what type of intrusion forms and why different sheet intrusion types are found in the same volcanic plumbing system [143]. These examples highlight only some of the benefits of combining different approaches.

Finally, we need to educate future generations of researchers to be open-minded and aware of the potential of other methods and to be able to communicate their research to a broad audience and collaborate beyond disciplinary boundaries. Multidisciplinary conferences, workshops, and research meetings should become common practice, as should multidisciplinary research projects and teams.

It is up to each one of us to try and understand the potential of other disciplines’ methods, to communicate our research so that we can be understood beyond our own disciplines, and to educate our students to be open-minded and have a broad overview of the approaches in volcanic plumbing system research. Although this development may at first glance look like an about-face from the trend of specialisation in our disciplines, in reality, it implies the chance for a more thorough characterisation of the complexity of volcanic plumbing systems.



The authors are grateful for fruitful discussions with the participants, and in particular our coorganiser Valentin Troll, of the 2013 workshop “Mechanics of magma emplacement and volcanotectonics”, which was funded by an ESF grant within the MeMoVolc scheme. Furthermore, we would like to thank Maria de los Angeles García Juanatey for discussions on geophysical methods and all our colleagues who provided photographs of plumbing system researchers at work through ResearchGate and other channels. The authors also acknowledge funding of their research from the Swedish Research Council (VR), the Swedish Royal Academy of Sciences (KVA), the Centre for Natural Disaster Science (CNDS), the Department of Earth Sciences, Uppsala University and the University in Oslo.


  1. 1. Jerram DA, Bryan SE. Plumbing systems of shallow level intrusive complexes. In: Breitkreuz C, Rocchi E, editors. Physical Geology of Shallow Magmatic Systems. Advances in Volcanology. Springer. Edinburgh 2015.
  2. 2. Hutton J. Abstract of a Dissertation read in the Royal Society of Edinburgh, upon the Seventh of March, and Fourth of April, MDCCLXXXV, concerning the System of the Earth, its Duration, and Stability, Royal Society of Edinburgh, Edinburgh 1785.
  3. 3. Hall J. Experiments on whinstone and lava. Journal of Natural Philosophy, Chemistry, and the Arts. 1800;4:8–18, 56–65.
  4. 4. Young DA. Mind over Magma: The Story of Igneous Petrology. Princeton University Press; Princeton University Press, Princeton. 2003. 686 p.
  5. 5. Gilbert GK. Report on the Geology of the Henry Mountains. Washington, DC., Government Printing Office; 1877.
  6. 6. Suess E. The Face of the Earth Germany. Freytag, Leipzig, Germany. 1885.
  7. 7. Daly RA. The mechanics of igneous intrusions. American Journal of Science. 1903;15:269–298.
  8. 8. Daly RA. Igneous Rocks and Their Origin. McGraw Hill; New York 1914.
  9. 9. Clough CT, Maufe HB, Bailey EB. The cauldron subsidence of Glen Coe, and associated igneous phenomena. Quarterly Journal of the Geological Society of London. 1909;65:611–678.
  10. 10. Bailey EB, Clough CT, Wright WB, Richey JE, Wilson GV. Tertiary and Post-Tertiary Geology of Mull, Loch Aline, and Oban. Geological Survey of Scotland Memoir. Edinburgh 1924.
  11. 11. Bowen NL. Crystallisation-differentiation of silicate liquids. American Journal of Science. 1915;39:175–191.
  12. 12. Grout FF. A type of igneous differentiation. Journal of Geology. 1918;26:626–658.
  13. 13. Hall AL. The Bushveld Igneous Complex of the Central Transvaal. Geological Survey of South Africa Memoir. Pretoria Government Printer, Pretoria 1932;28.
  14. 14. Wager LR, Deer WA. Geological investigations in East Greenland, part III. The petrology of the Skaergaard Intrusion, Kangerdluqssuad, East Greenland. Meddelelser om Grønland. 1939;105:1–352.
  15. 15. Walker GPL. Zeolite zones and dyke distribution in relation to the structure of eastern Iceland. Journal of Geology. 1960;68:515–528.
  16. 16. Sparks RSJ. The legacy of George Walker to volcanology. In: Thordarson T, Self S, Larsen G, Rowland SK, Höskuldsson A, editors. Studies in Volcanology. The Legacy of George Walker. Special Publications of IAVCEI. London 2009;2:pp. 1–15.
  17. 17. Hubbert MK. Theory of scale models as applied to the study of geologic structures. Geological Society of America Bulletin. 1937;48:1459–1520.
  18. 18. Ramberg H, editor. Model studies in relation to intrusion of plutonic bodies, Mechanism of igneous intrusion. Geological Journal Special Issue. 1970;2.
  19. 19. Ramberg H. Gravity, Deformation and the Earth’s Crust. Academic Press; Cambridge, Massachusetts 1981.
  20. 20. Galland O, Holohan EP, van Wyk de Vries B, Burchardt S. Laboratory Modelling of Volcano Plumbing Systems: A Review. In: Breitkreuz C, Rocchi S, editors. Physical Geology of Shallow Magmatic Systems – Dykes, Sills and Laccoliths. Springer; Advances in Volcanology. Berlin 2015;pp. 1–68.
  21. 21. Bott MHP, Tuson J. Deep structure beneath the Tertiary volcanic regions of Skye, Mull and Ardnamurchan, North-west Scotland. Nature Physical Science. 1973;242:114–116.
  22. 22. Ameglio L, Vigneresse J-L. Geophysical imaging of the shape of granitic intrusions at depth: a review. In: Castro A, Fernández C, Vigneresse J-L, editors. Understanding Granites: Integrating New and Classical Techniques. Geological Society, London, Special Publications. 1999;168:pp. 39–54.
  23. 23. Walker GPL. The Breiddalur central volcano, eastern Iceland. Quarterly Journal of the Geological Society. 1963;119:29–63.
  24. 24. Mahan KH, Bartley JM, Coleman DS, Glazner AF, Carl BS. Sheeted intrusion of the synkinematic McDoogle pluton, Sierra Nevada, California. Geological Society of America Bulletin. 2003;115:1570–1582.
  25. 25. Emeleus CH, Cheadle MJ, Hunter RH, Upton BGJ, Wadsworth WJ. The rum layered suite. Developments in Petrology. 1996;15:403–439.
  26. 26. Cruden AR, McCaffrey KJW. Growth of plutons by floor subsidence: implications for rates of emplacement, intrusion spacing and melt-extraction mechanisms. Physics and Chemistry of the Earth, Part A: Solid Earth and Geodesy. 2001;26:303–315.
  27. 27. Paterson SR, Fowler TK, Schmidt KL, Yoshinobu AS, Yuan ES, Miller RB. Interpreting magmatic fabric patterns in plutons. Lithos. 1998;44:53–82.
  28. 28. O’Driscoll B, Troll VR, Reavy RJ, Turner P. The Great Eucrite intrusion of Ardnamurchan, Scotland: reevaluating the ring-dike concept. Geology. 2006;34:189–192.
  29. 29. Stevenson CT, Owens WH, Hutton DH, Hood DN, Meighan IG. Laccolithic, as opposed to cauldron subsidence, emplacement of the Eastern Mourne pluton, N. Ireland: evidence from anisotropy of magnetic susceptibility. Journal of the Geological Society. 2007;164(1):99–110.
  30. 30. Roni E, Westerman DS, Dini A, Stevenson C, Rocchi S. Feeding and growth of a dyke–laccolith system (Elba Island, Italy) from AMS and mineral fabric data. Journal of the Geological Society. 2014;171:413–424.
  31. 31. Tibaldi A, Pasquarè AF, Rust D. New insights into the cone sheet structure of the Cuillin Complex, Isle of Skye, Scotland. Journal of the Geological Society. 2011;168:689–704.
  32. 32. Krumbholz M, Hieronymus C, Burchardt S, Troll VR, Tanner D, Friese N. Weibull distributed dyke thickness reflects probalistic character of host-rock strength. Nature Communications. 2014;5:3272.
  33. 33. Burchardt S, Tanner DC, Troll VR, Krumbholz M, Gustafsson LE. Three-dimensional geometry of concentric intrusive sheet swarms in the Geitafell and the Dyrfjöll Volcanoes, Eastern Iceland. Geochemistry Geophysics Geosystems. 2011;12:Q0AB09.
  34. 34. Burchardt S, Troll VR, Mathieu L, Emeleus HC, Donaldson CH. Ardnamurchan 3D cone-sheet architecture explained by a single elongate magma chamber. Scientific Reports. 2013;3:2891.
  35. 35. Mathieu L, Burchardt S, Troll VR, Krumbholz M, Delcamp A. Geological constraints on the dynamic emplacement f cone-sheets – the Ardnamurchan cone-sheet swarm, NW Scotland. Journal of Structural Geology. 2015;80:133–141.
  36. 36. Buckley SJ, Enge HD, Carlsson C, Howell JA. Terrestrial laser scanning for use in virtual outcrop geology. The Photogrammetric Record. 2010;25:225–239. DOI: 10.1111/j.1477-9730.2010.00585.x.
  37. 37. Senger K, Buckley SJ, Chevallier L, Fagreng Å, Galland O, Kurz TKO, Planke S, Tveranger J. Fracturing of doleritic intrusions and associated contact zones: insights from the Eastern Cape, South Africa. Journal of African Earth Sciences. 2015;102:70–85.
  38. 38. Townsend M, Pollard DD, Johnson K, Culha C. Jointing around magmatic dikes as a precursor to the development of volcanic plugs. Bulletin of Volcanology. 2015;77:1–13. DOI: 10.1007/s00445-015-0978-z.
  39. 39. Parcheta CE, Pavlov CA, Wiltsie N, Carpenter KC, Nash J, Parness A, Mitchell KL. A robotic approach to mapping post-eruptive volcanic fissure conduits. Journal of Volcanology and Geothermal Research. DOI:
  40. 40. Burchardt S, Tanner DC, Krumbholz M. Mode of emplacement of the Slaufrudalur Pluton, Southeast Iceland inferred from three-dimensional GPS mapping and model building. Tectonophysics. 2010;480(1):232–240. DOI: 10.1016/j.tecto.2009.10.010.
  41. 41. Burchardt S, Tanner D, Krumbholz M. The Slaufrudalur pluton, southeast Iceland—an example of shallow magma emplacement by coupled cauldron subsidence and magmatic stoping. Geological Society of America Bulletin. 2012;124(1–2), 213–227.
  42. 42. Burchardt S. New insights in the mechanics of sill emplacement provided by field observations of the Njardvik Sill, Northeast Iceland. Journal of Volcanology and Geothermal Research. 2008;173:280–288.
  43. 43. Langmuir CH. Geochemical consequences of in situ crystallization. Nature. 1989;340:199–205.
  44. 44. Rickwood PC. Boundary lines within petrologic diagrams which use oxides of major and minor elements. Lithos. 1989;22:247–263.
  45. 45. Spera FJ, Bohrson WA. Energy-constrained open-system magmatic processes I: general model and energy-constrained assimilation and fractional crystallization (EC-AFC) formulation. Journal of Petrology. 2001;42(5):999–1018.
  46. 46. Glazner AF, Bartley JM, Coleman DS, Gray W, Taylor RZ. Are plutons assembled over millions of years by amalgamation from small magma chambers? GSA Today. 2004;14:4–12.
  47. 47. Ginibre C, Wörner G, Kronz A. Crystal zoning as an archive for magma evolution. Elements. 2007;3:261–266.
  48. 48. Mock A, Jerram DA. Crystal size distributions (CSD) in three dimensions: insights from the 3D reconstruction of a highly porphyritic rhyolite. Journal of Petrology. 2005;46:1525–1541.
  49. 49. Cashman K, Blundy J. Petrological cannibalism: the chemical and textural consequences of incremental magma body growth. Contributions to Mineralogy and Petrology. 2013;166:703–729.
  50. 50. Saunders K, Buse B, Kilburn MR, Kearns S, Blundy J. Nanoscale characterisation of crystal zoning. Chemical Geology. 2014;364:20–32.
  51. 51. Dahren B, Troll VR, Andersson UB, Chadwick JP, Gardner MF, Jaxybulatov K, Koulakov I. Magma plumbing beneath Anak Krakatau volcano, Indonesia: evidence for multiple magma storage regions. Contributions to Mineralogy and Petrology. 2012;163:631–651.
  52. 52. Chadwick JP, Troll VR, Waight TE, van der Zwan FM, Schwarzkopf LM. Petrology and geochemistry of igneous inclusions in recent Merapi deposits: a window into the sub-volcanic plumbing system. Contributions to Mineralogy and Petrology. 2013;165:259–282.
  53. 53. Muir DD, Blundy JD, Rust AC, Hickey J. Experimental constraints on dacite pre-eruptive magma storage conditions beneath Uturuncu volcano. Journal of Petrology. 2014;55:749–767.
  54. 54. Deegan FM, Troll VR, Freda C, Misiti V, Chadwick JP, McLeod CL, Davidson JP. Magma–carbonate interaction processes and associated CO2 release at Merapi Volcano, Indonesia: insights from experimental petrology. Journal of Petrology. 2010;51:1027–1051.
  55. 55. Jolis EM, Freda C, Troll VR, Deegan FM, Blythe LS, McLeod CL, Davidson JP. Experimental simulation of magma-carbonate interaction beneath Mt. Vesuvius, Italy. Contributions to Mineralogy and Petrology. 2013;166:1335–1353.
  56. 56. Troll VR, Klügel A, Longpré M-A, Burchardt S, Deegan FM, Carracedo JC, Wiesmaier S, Kueppers U, Dahren B, Blythe LS, Hansteen T, Freda C, Budd D, Jolis EM, Jonsson E, Meade F, Berg S, Mancini L, Polacci M. Floating stones off El Hierro, Canary Islands: xenoliths of pre-island sedimentary origin in the early products of the October 2011 eruption. Solid Earth. 2012;3:97–100.
  57. 57. Geiger H, Mattsson T, Deegan FM, Troll VR, Burchardt S, Gudmundsson O, Tryggvason A, Krumbholz M, Harris C. Magma plumbing for the 2014–2015 Holuhraun eruption Iceland. Geochemistry Geophysics Geosystems DOI: 10.1002/2016GC006317.
  58. 58. Budd DA, Troll VR, Dahren B, Burchardt S. 2016. Persistent multitiered plumbing beneath Katla volcano, Iceland. Geochemistry Geophysics Geosystems. DOI: 10.1002/2015GC006118.
  59. 59. Andrews-Hanna JC, Asmar SW, Head JW, Kiefer, WS, Konopliv AS, Lemoine FG et al. Ancient igneous intrusions and early expansion of the Moon revealed by GRAIL gravity gradiometry. Science. 2013;339(6120):675–678.
  60. 60. Améglio L, Vigneresse J. L. Bouchez, D. H. W. Hutton, W. E. Stephens. Granite pluton geometry and emplacement mode inferred from combined fabric and gravity data. In: Granite: From Segregation of Melt to Emplacement Fabrics. Springer. Berlin 1997; pp. 199–214.
  61. 61. Drenth BJ, Keller GR, Thompson RA. Geophysical study of the San Juan Mountains batholith complex, southwestern Colorado. Geosphere. 2012;8:669–684.
  62. 62. Burchardt S, Walter TR. Propagation, linkage, and interaction of caldera ring-faults: comparison between analogue experiments and caldera collapse at Miyakejima, Japan, in 2001. Bulletin of Volcanology. 2010;72:297–308. DOI: 10.1007/s00445-009-0321-7.
  63. 63. Hansen DM, Cartwright JA. Saucer-shaped sill with lobate morphology revealed by 3D seismic data: implications for resolving a shallow-level sill emplacement mechanism. Journal of the Geological Society. 2006;163(3):509–523.
  64. 64. Jaxybulatov K, Shapiro NM, Koulakov I, Mordret A, Landès M, Sens-Schönfelder C. A large magmatic sill complex beneath the Toba caldera. Science. 2014;346(6209):617–619.
  65. 65. Battaglia M, Gottsmann J, Carbone D, Fernández J. 4D volcano gravimetry. Geophysics. 2008;73:WA3–WA18.
  66. 66. Zobin VM. Introduction to Volcanic Seismology. Elsevier; Amsterdam 2012.
  67. 67. Sigmundsson F, Hooper A, Hreinsdóttir S, Vogfjörd KS, Ófeigsson BG, Heimisson ER., et al. Segmented lateral dyke growth in a rifting event at Bardarbunga volcanic system, Iceland. Nature. 2015;517:191–195.
  68. 68. White RS, Drew J, Martens HR, Key J, Soosalu H, Jakobsdóttir SS. Dynamics of dyke intrusion in the mid-crust of Iceland. Earth and Planetary Science Letters. 2011;304:300–312.
  69. 69. Einarsson P. S-wave shadows in the Krafla caldera in NE-Iceland, evidence for a magma chamber in the crust. Bulletin of Volcanology. 1978;41:187–195.
  70. 70. Thomson K. Volcanic features of the North Rockall Trough: application of visualisation techniques on 3D seismic reflection data. Bulletin of Volcanology. 2004;67:116–128.
  71. 71. Planke S, Rasmussen T, Rey SS, Myklebust R. Seismic characteristics and distribution of volcanic intrusions and hydrothermal vent complexes in the Vøring and Møre basins. In: Doré AG, Vining BA, editors. Proceedings of the 6th Petroleum Geology Conference of the Geological Society, London. 2005.
  72. 72. Polteau S, Mazzini A, Galland O, Planke S, Malthe-Sørenssen A. Saucer-shaped intrusions: occurrences, emplacement and implications. Earth and Planetary Science Letters. 2008;266:195–204.
  73. 73. Jackson CAL, Schofield N, Golenkov B. Geometry and controls on the development of igneous sill-related forced-folds: a 2D seismic reflection case study from offshore southern Australia. Geological Society of America Bulletin. 2013;125:1874–1890.
  74. 74. Hansen DM, Cartwright JA, Thomas D. 3D seismic analysis of the geometry of igneous sills and sill junction relationships. In: Davies RJ, Cartwright J, Stewardt SA, Lappin M, Underhill JR, editors. 3D Seismic Technology: Application to the Exploration of Sedimentary Basins. Geological Society of London Memoir, London. 2004.
  75. 75. Rateau R, Schofield N, Smith M. The potential role of igneous intrusions on hydrocarbon migration, West of Shetland. Petroleum Geoscience. 2013;19:259–272. DOI: 10.1144/petgeo2012-035.
  76. 76. Schofield N, Heaton L, Holford SP, Archer, SG, Jackson CAL, Jolley DW. Seismic imaging of “broken bridges”: linking seismic to outcrop-scale investigations of intrusive magma lobes. Journal of the Geological Society. 2012;169:421–426. DOI: 10.1144/0016-76492011-150.
  77. 77. Trude J, Cartwright J, Davies RJ, Smallwood J. New technique for dating igneous sills. Geology. 2003;31:813–816.
  78. 78. Hansen DM, Cartwright JA. The three-dimensional geometry and growth of forced folds above saucer-shaped igneous sills. Journal of Structural Geology. 2006b;28:1520–1535.
  79. 79. Jackson CAL, Schofield N. The influence of normal fault geometry on igneous sill emplacement and morphology. Geology. 2013;41:407–410. DOI: 10.1130/g33824.1.
  80. 80. Dzurisin D. Volcano deformation: new geodetic monitoring techniques. Springer; Berlin 2007.
  81. 81. Lagios E, Sakkas V, Novali F, Bellotti F, Ferretti A, Vlachou K, Dietrich V. SqueeSAR and GPS ground deformation monitoring of Santorini Volcano (1992–2012): tectonic implications. Tectonophysics. 2013;594:38–59.
  82. 82. Chadwick Jr, WW, Jónsson S, Geist DJ, Poland M, Johnson DJ, Batt S, Harpp KS, Ruiz A. The May 2005 eruption of Fernandina volcano, Galápagos: the first circumferential dike intrusion observed by GPS and InSAR. Bulletin of Volcanology. 2011;73(6):679–697.
  83. 83. Wauthier C, Cayol V, Smets B, d’Oreye N, Kervyn F. Magma pathways and their interactions inferred from InSAR and stress modeling at Nyamulagira Volcano, DR Congo. Remote Sensing. 2015;7(11):15179–15202.
  84. 84. Lisowski M. Analytical volcano deformation source models. In: Dzurisin D, editor. Volcano Deformation: New Geodetic Monitoring Techniques. Springer. Berlin 2007.
  85. 85. Pritchard ME, Simons M. A satellite geodetic survey of large-scale deformation of volcanic centres in the central Andes. Nature. 2002;418:167–171.
  86. 86. Sigmundsson F, Hreinsdóttir S, Hooper A, Arnadottir, T, Pedersen R, Roberts MJ, Oskarsson N, Auriac A, Decriem J, Einarsson P, Geirsson H, Hensch M, Ofeigsson BG, Sturkell E, Sveinbjornsson H, Feigl KL. Intrusion triggering of the 2010 Eyjafjallajökull explosive eruption. Nature. 2010;468:426–430. DOI: 10.1038/nature09558.
  87. 87. Sturkell E et al. New insights into volcanic activity from strain and other deformation data for the Hekla 2000 eruption. Journal of Volcanology and Geothermal Research. 2013;256:78–86.
  88. 88. Klein FW. Eruption forecasting at Kilauea volcano, Hawaii. Journal of Geophysical Research: Solid Earth. 1984;89:3059–3073.
  89. 89. Sparks RSJ. Forecasting volcanic eruptions. Earth and Planetary Science Letters. 2003;210:1–15.
  90. 90. Cervelli PF, Fournier T, Freymueller J, Power JA. Ground deformation associated with the precursory unrest and early phases of the January 2006 eruption of Augustine Volcano, Alaska. Geophysical Research Letters. 2006;33:L18304.
  91. 91. Dieterich JH, Decker RW. Finite element modeling of surface deformation associated with volcanism. Journal of Geophysical Research. 1975;80:4094–4102.
  92. 92. Fialko Y, Khazan Y, Simons M. Deformation due to a pressurized horizontal circular crack in an elastic half-space, with applications to volcano geodesy. Geophysical Journal International. 2001;146:181–190.
  93. 93. Murray MH, Marshall GA, Lisowski M, Stein RS. The 1992 M=7 Cape Mendocino, California earthquake: Coseismic deformation at the south end of the Cascadia megathrust. Journal of Geophysical Research. 1996;101:17707–17725.
  94. 94. Cervelli PF, Murray MH, Segall P, Aoki Y, Kato T. Estimating source parameters from deformation data, with an application to the March 1997 earthquake swarm off the Izu Peninsula, Japan. Journal of Geophysical Research Solid Earth. 2001;106:11217–11237.
  95. 95. Hreinsdóttir S, Sigmundsson F, Roberts MJ, Björnsson H, Grapenthin R, Arason P. Volcanic plume height correlated with magma-pressure change at Grimsvotn Volcano, Iceland. Nature Geoscience. 2014;7:214–218.
  96. 96. Galland O, Bertelsen HS, Guldstrand F, Girod L, Johannessen RF, Bjugger F, Burchardt S, Mair K. Application of open-source photogrammetric software MicMac for monitoring surface deformation in laboratory models. Journal of Geophysical Research. 2016. DOI: 10.1002/2015JB012564
  97. 97. Galland O, Gisler GR, Haug ØT. Morphology and dynamics of explosive vents through cohesive rock formations. Journal of Geophysical Research. 2014:119:LB011050, doi: 10.1002/2014JB011050.
  98. 98. Galland O, Scheibert J. Analytical model of surface uplift above axisymmetric flat-lying magma intrusions: implications for sill emplacement and geodesy. Journal of Volcanology and Geothermal Research. 2013;253:114–130. DOI:
  99. 99. Daubrée A. Recherches expérimentales sur le rôle possible des gaz à hautes températures doués de très fortes pressions et animés d’un mouvement fort rapide dans divers phénomènes géologiques. Bulletin de la Société géologique de France. 1891;19:313–354.
  100. 100. Barenblatt GI. Scaling. Cambridge University Press; Cambridge 2003.
  101. 101. Hubbert MK, Willis DG. Mechanics of hydraulic fracturing. AIME Transactions. 1957;210:153–168.
  102. 102. Takada A. Experimental study on propagation of liquid-filled crack in gelatin: shape and velocity in hydrostatic stress condition. Journal of Geophysical Research. 1990;95:8471–8481.
  103. 103. Román-Berdiel T, Gapais D, Brun J-P. Analogue models of laccolith formation. Journal of Structural Geology. 1995;17:1337–1346.
  104. 104. Corti G, Bonini M, Innocenti F, Manetti P, Mulugeta G. Centrifuge models simulating magma emplacement during oblique rifting. Journal of Geodynamics. 2001;31:557–576.
  105. 105. Galland O, Cobbold PR, Hallot E, de Bremond d’Ars J, Delavaud G. Use of vegetable oil and silica powder for scale modelling of magmatic intrusion in a deforming brittle crust. Earth and Planetary Science Letters. 2006;243:786–804.
  106. 106. Galland O. Experimental modelling of ground deformation associated with shallow magma intrusions. Earth and Planetary Science Letters. 2012;317–318:145–156. DOI: 10.1016/j.epsl.2011.10.017.
  107. 107. Abdelamak MM, Bulois C, Mourgues R, Galland O, Legland JB, Gruber C. Description of new dry granular materials of variable cohesion and friction coefficient: implications for laboratory modelling of the brittle crust. Tectonophysics. 2016. in press:
  108. 108. Hirata T. Fracturing due to fluid intrusion into viscoelastic materials. Physical Review E. 1998;57:1772–1779.
  109. 109. Sumita I, Ota Y. Experiments on buoyancy-driven crack around the brittle–ductile transition. Earth and Planetary Science Letters. 2011;304:337–346. DOI:
  110. 110. Di Giuseppe E, Corbi F, Funiciello F, Massmeyer A, Santimano TN, Rosenau M, Davaille A. Characterization of Carbopol® hydrogel rheology for experimental tectonics and geodynamics. Tectonophysics. 2015;642:29–45. DOI: tp://
  111. 111. Vergniolle S, Jaupart C. Separated two-phase flow and basaltic eruptions. Journal of Geophysical Research: Solid Earth. 1986;91:12842–12860. DOI: 10.1029/JB091iB12p12842.
  112. 112. Pioli L, Bonadonna C, Azzopardi BJ, Phillips JC, Ripepe M. Experimental constraints on the outgassing dynamics of basaltic magmas. Journal of Geophysical Research: Solid Earth. 2012;117. DOI: 10.1029/2011JB008392.
  113. 113. Leever KA, Galland O, Acocella V. The science behind laboratory-scale models of the earth. Eos, Transactions American Geophysical Union. 2014;95:30. DOI: 10.1002/2014eo030008.
  114. 114. Leever KA, Gabrielsen RH, Sokoutis D, Willingshofer E. The effect of convergence angle on the kinematic evolution of strain partitioning in transpressional brittle wedges: insight from analog modeling and high-resolution digital image analysis. Tectonics. 2011;30:TC2013. DOI: 10.1029/2010tc002823.
  115. 115. Abdelmalak MM, Mourgues R, Galland O, Bureau D. Fracture mode analysis and related surface deformation during dyke intrusion: Results from 2D experimental modelling. Earth and Planetary Science Letters. 2012;359–360:93–105. DOI: 10.1016/j.epsl.2012.10.008
  116. 116. Adam J, Klinkmüller M, Schreurs G, Wieneke B. Quantitative 3D strain analysis in analogue experiments simulating tectonic deformation: integration of X-ray computed tomography and digital volume correlation techniques. Journal of Structural Geology. 2013;55:127–149. DOI:
  117. 117. Poppe S, Holohan EP, Pauwels E, Cnudde V, Kervyn M. Sinkholes, pit craters, and small calderas: analog models of depletion-induced collapse analyzed by computed X-ray microtomography. Geological Society of America Bulletin. 2015;127:281–296, DOI: 10.1130/b30989.1.
  118. 118. Lohrmann J, Kukowski N, Adam J, Onken O. The impact of analogue materials properties on the geometry, kinematics, and dynamics of convergent sand wedges. Journal of Structural Geology. 2003;25:1691–1711.
  119. 119. Rosenau M, Nerlich R, Brune S, Oncken O. Experimental insights into the scaling and variability of local tsunamis triggered by giant subduction megathrust earthquakes. Journal of Geophysical Research. 2010;115:B09314. DOI: 10.1029/2009jb007100.
  120. 120. Taisne B, Tait S. Effect of solidification on a propagating dike. Journal of Geophysical Research: Solid Earth. 2011;116:B01206, DOI: 10.1029/2009jb007058.
  121. 121. Pollard DD. Derivation and evaluation of a mechanical model for sheet intrusions. Tectonophysics. 1973;19:233–269, DOI:
  122. 122. Sparks RSJ, Huppert HE, Turner JS, Sakuyama M, O’Hara MJ. The fluid dynamics of evolving magma chambers [and discussion]. Philosophical Transactions of the Royal Society of London, Series A. 1984;310:511–534.
  123. 123. Green AE. On Boussinesq’s problem and penny-shaped cracks. Mathematical Proceedings of the Cambridge Philosophical Society. 1949;45:251–257.
  124. 124. Weertman J. Theory of water-filled crevasses in glaciers applied to vertical magma transport beneath oceanic ridges. Journal of Geophysical Research. 1971;76:1171–1183.
  125. 125. Petford N, Mirhadizadeh S. Eddy flow and associated particle dynamics during magma intrusion: the Basement Sill, Antarctica. EGU General Assembly Conference Abstracts. 2014;4182.
  126. 126. Chestler SR, Grosfils EB. Using numerical modeling to explore the origin of intrusion patterns on Fernandina volcano, Galápagos Islands, Ecuador. Geophysical Research Letters. 2013;40:50833. DOI: 10.1002/grl.50833.
  127. 127. Holohan EP, Schöpfer MPJ, Walsh JJ. Mechanical and geometric controls on the structural evolution of pit crater and caldera subsidence. Journal of Geophysical Research. 2011;116:B07202. DOI: 10.1029/2010jb008032.
  128. 128. Keller T, May DA, Kaus BJP. Numerical modelling of magma dynamics coupled to tectonic deformation of lithosphere and crust. Geophysical Journal International. 2013;195:1406–1442, DOI: 10.1093/gji/ggt306.
  129. 129. Malthe-Sørenssen A, Planke S, Svensen H, Jamtveit B. Formation of saucer-shaped sills. In: Breitkreuz C, Petford N, editors. Physical geology of high-level magmatic systems. Geological Society of London. Special Publication. London 2004;pp. 215–227.
  130. 130. Longo A, Vassalli M, Papale P, Barsanti M. Numerical simulation of convection and mixing in magma chambers replenished with CO2-rich magma. Geophysical Research Letters. 2006;33:027760. DOI: 10.1029/2006GL027760.
  131. 131. Aarnes I, Svensen H, Polteau S, Planke S. Contact metamorphic devolatilization of shales in the Karoo Basin, South Africa, and the effects of multiple sill intrusions. Chemical Geology. 2011;281:181–194. DOI: 10.1016/j.chemgeo.2010.12.007
  132. 132. Gerya TV, Burg J-P. Intrusion of ultramafic magmatic bodies into the continental crust: numerical simulation. Physics of the Earth and Planetary Interiors. 2007;160:124–142.
  133. 133. Amelung F, Jonsson S, Zebker H, Segall P. Widespread uplift and “trapdoor” faulting on Galápagos volcanoes observed with radar interferometry. Nature. 2000;407:993–996. DOI: 10.1038/35039604.
  134. 134. Aarnes I, Podladchikov YY, Svensen,H. Devolatilization-induced pressure build-up: implications for reaction front movement and breccia pipe formation. Geofluids. 2012;12:265–279. DOI: 10.1111/j.1468-8123.2012.00368.x.
  135. 135. Weatherley D, Hancock W, Abe S, Boros V. ESyS-Particle Tutorial and User’s Guide Version 2.2. 2. 2013.
  136. 136. Zhao C, Hobbs BE, Ord A, Peng S. Particle simulation of spontaneous crack generation associated with the laccolithic type of magma intrusion processes. International Journal for Numerical Methods in Engineering. 2008;75:1172–1193. DOI: 10.1002/nme.2287.
  137. 137. Bond CE, Gibbs AD, Shipton ZK, Jones S. What do you think this is? “Conceptual uncertainty” in geoscience interpretation. GSA Today. 2007;17:4–10.
  138. 138. Gudmundsson O, Brandsdottir B, Menke W, Sigvaldason GE. The crustal magma chamber of the Katla volcano in south Iceland revealed by 2-D seismic undershooting. Geophysical Journal International. 1994;119:277–296.
  139. 139. Sturkell E, Sigmarsson F. Recent unrest and magma movements at Eyjafjallajökull and Katla volcanoes, Iceland. Journal of Geophysical Research. 2003;108:2369.
  140. 140. Oladottir BA, Sigmarsson O, Larsen G, Thordarson T. Katla volcano, Iceland: magma composition, dynamics and ruption frequency as recorded by Holocene tephra layers. Bulletin of Volcanology. 2008;70:475–493.
  141. 141. Kavanagh JL, Sparks RSJ. Insights of dyke emplacement mechanics from detailed 3D dyke thickness datasets. Journal of the Geological Society. 2011;168:965–978.
  142. 142. Riel B, Milillo P, Simons M, Lundgren P, Kanamori H, Samsonov S. The collapse of Bárðarbunga caldera, Iceland. Geophysical Journal International. 2015;202:446–453.
  143. 143. Galland O, Burchardt S, Hallot E, Mourgues R. Dynamics of dikes versus cone sheets in volcanic systems. Journal of Geophysical Research – Solid Earth. 2014;119:6178–6192.

Written By

Steffi Burchardt and Olivier Galland

Reviewed: 26 April 2016 Published: 21 September 2016