Abstract

Analogue models are commonly used to gain insights into large-scale volcano-tectonic processes. Documenting model surface topography and the three-dimensional (3D) aspect of deformation structures remains the greatest challenge in understanding the simulated processes. Here we present the results of volcano analogue models imaged with an X-ray computerized micro-tomography (μCT) system developed at the Ghent University Centre for Tomography (UGCT). Experiments simulate volcano deformation due to gravitational loading over a ductile layer, a process affecting many natural volcanoes built over a sedimentary substratum. Results show that μCT is able to provide a 3D reconstruction of the model topography with unprecedented resolution. Virtual cross sections through reconstructed models enable us to map the main structures at depth and to document the deformation of the brittle-ductile interface due to contrasting X-ray attenuation. Results for lateral spreading and vertical sagging into thin and thick ductile layers, respectively, are illustrated for circular cones and elongated ridges. Results highlight structural patterns not seen in previous models, such as: 1) the 3D form of a polygonal brecciated zone at the center of spreading cones; 2) the complete lack of such a zone in sagging cones; and 3) relay structures between graben-bounding faults in spreading cones. In addition, detailed imaging of tension gashes and of the flexure surface below sagging cones enables the 3D strain distribution to be explored. Experiments with non-cohesive and low cohesion granular materials present striking differences in surface topography and fault characteristics. Despite limitations associated with the scan duration, μCT reconstruction of analogue models appears a powerful tool for better understanding the complex 3D deformation associated with volcano-tectonic processes.

INTRODUCTION

Analogue modeling has been used for over a century to simulate large scale geological processes. Laboratory models enable researchers to understand the mechanisms controlling the kinematics of large scale deformation processes and the geometry of the resulting structures. In volcanology in particular, physical models enable the investigation of processes such as volcano-scale deformation under gravitational loading (e.g., Merle and Borgia, 1996), volcano flank collapse (e.g., Donnadieu and Merle, 1998; Acocella, 2005; Andrade and van Wyk de Vries 2010), volcanic debris avalanche (Shea and de Vries, 2008), interactions of volcanic and tectonic structures (e.g., Lagmay et al., 2000; Merle et al., 2001), as well as intrusive processes (e.g., Mathieu et al., 2008; Kervyn et al., 2009) and caldera formation (e.g., Marti et al., 1994; Roche et al., 2000) on a time frame and at a scale unobtainable by other means. The main objective of analogue models is to create a quantitatively scaled analogue of a natural process in order to analyze the relationships between the imposed boundary conditions and the resulting structures.

This study focuses on analogue models simulating volcano deformation due to gravitational loading over a ductile substratum. This process was first proposed by van Bemmelen (1949) and followed by Borgia et al. (1992) to explain structures observed at natural volcanoes. Later, factors controlling the type of deformation were studied with analogue experiments, using silicone to simulate the viscous material deforming under the volcano load (e.g., Merle and Borgia, 1996; Delcamp et al., 2008; Le Corvec and Walter, 2009) and numerical models (e.g., van Wyk de Vries and Matela, 1998; Morgan and McGovern, 2005). Gravitational loading over a ductile layer is associated with radial spreading of the volcanic cone due to lateral flow of the underlying viscous material, such as clayey sediments. If the viscous substratum is thicker, such as asthenospheric mantle material, thick sedimentary basins, or the crusts over long time spans, the cone tends to sag vertically into the substratum (see van Wyk de Vries and Matela, 1998). Evidence for deformation due to gravitational loading has now been recognized at many volcanoes, from strato-volcanoes (e.g., Borgia and van Wyk de Vries, 2003) to large terrestrial and planetary shields (e.g., Morgan and McGovern, 2005; Byrne et al., 2009; Le Corvec and Walter, 2009).

Detailed documentation of the surface topography, vertical and horizontal surface deformations, and structures accommodating this deformation within the analogue model is essential for exploring the ongoing deformation processes and for enabling comparison with natural volcanoes. Previously, analogue model deformations have been analyzed based on plan view pictures acquired at regular intervals during experiments. Analysis is often restricted to the extraction of horizontal deformation, based on displacement of discrete point markers (e.g., Kervyn et al., 2009). More advanced techniques for topographic reconstruction of analogue models, based on multiple images (Cecchi et al., 2003) or stereoscopic image acquisition (Donnadieu et al., 2003; Delcamp et al., 2008), enabled the characterization of vertical and horizontal deformation, as well as the resulting model topography with mm-scale accuracy when sufficient surface texture is available.

Dip and kinematics of structures accommodating the deformation are either extrapolated from surface traces of these structures and the observed surface deformation, or directly observed from serial cross-sections. Cross-sections require that the model be physically cut, after having been impregnated with water to consolidate the model. Structures are then analyzed based on relative movement of colored layers originally built into the model. This procedure has certain limitations: tension gashes are closed by impregnation and compaction, structures with limited surface displacement might be undetected, and the model is partially or totally destroyed, limiting the number of orientations analyzed. Cross-sections are further complicated by the presence of viscous material in the model that cannot be consolidated (e.g., silicone, golden syrup). Experiments were also carried out in two-dimensional (2D) configurations against a transparent glass wall, enabling detailed characterization of the development of structures at depth during deformation using advanced image correlation techniques (e.g., Burchardt and Walter, 2009; Le Corvec and Walter, 2009). The 2D configuration and the potential boundary effects, however, limit the capability to extrapolate these observations to 3D systems.

Their conical shape causes volcanoes to deform in all directions under the effect of gravity, resulting in a complex network of structures with changing orientations that link and interact with each other. Computerized X-ray tomography (CT), as a non-destructive imaging system, has the capability to enable the 3D reconstruction of analogue models, permitting virtual cross-sections and the documentation of the 3D complexity of structures. This technique has been used to document analogue models simulating 2D tectonic deformation in compressive (e.g., Colletta et al., 1991; Schreurs et al., 2001), transcurrent (e.g., Schreurs, 1994; Ueta et al., 2000) and extensional (e.g., Mandl, 1988; Schreurs et al., 2003) regimes. These applications used a medical X-ray scanner, providing parallel X-ray cross-sections through the model, at intervals of a millimeter to several centimeters. These studies highlighted the capabilities of CT imaging for analysis of the 3D geometry and the kinematic evolution of internal structures. The structure identification is based on the fact that fault planes forming in granular material correspond to dilatant zones along which the normal arrangement of grains is disturbed, resulting in lower density and thus lower X-ray attenuation, irrespective of the type of fault (i.e., normal, strike-slip, or reverse).

The objective of this contribution is to test and highlight the added value and limitations of computerized X-ray micro-tomography (μCT) to image volcanic analogue models. Volcano gravitational deformation is used as an example of complex 3D topography and deformation. Preliminary results from end-member experiments with contrasted deformation styles, i.e., volcano spreading or sagging, and with contrasted volcano shape, i.e., circular cone or elongated ridges, are presented and analyzed. Specifically, the 3D results are used to highlight the impact of gravitational deformation on the internal volcano structure, including the formation of a brecciated central zone in spreading volcanoes and the lack of it in sagging volcanoes. The documentation of interactions at depth between faults observed at the surface and relationships between the surface and the silicone deformation are other original contributions of this study.

METHODOLOGY

Experimental Material, Setup, and Scaling

The experiment setup and materials used were similar to previous analogue modeling of gravitational deformation of volcanic cones (e.g., Merle and Borgia, 1996; Walter et al., 2006; Delcamp et al., 2008). Analogue models were built in cylindrical plastic containers, 15 cm in diameter and 5 cm in height, open in their upper part. The small size of the setup was dictated by the relationship between setup size and the resolution of the μCT images. The ductile sediment layer was simulated using a 4.5 or 11 mm thick layer of Polydimethylsiloxane (PDMS; Dow Corning SGM 36), a Newtonian fluid of 4 × 104 Pa s viscosity. A horizontal layer of granular material (i.e., analogue for the brittle crust), half the thickness of the ductile layer, was sifted over the PDMS layer, before emplacing a 2–2.5 cm high circular or elongated cone in the center of the model.

In order to be comparable with natural volcanoes, analogue models need to be geometrically, kinematically, and dynamically scaled. For scaling we used the approach developed in previous experimental work on volcano gravitational spreading (Merle and Borgia, 1996; Borgia et al., 2000; Walter, 2003; Walter et al., 2006; Delcamp et al., 2008). Table 1 shows that geometrical parameters are scaled with a model/nature ratio of 10−5 to 5 × 10−6. One centimeter in the model represents 1–2 km in nature. The dynamic scaling is achieved if the gravitational stress ratio σ* = ρ* × g* × h*, where ρ*, g* and h* are the model/nature ratios for the density, gravitational acceleration, and height of the volcanic cone, respectively. This calculation yields a stress ratio of ~2.5–5 × 10−6. Hence the analogue volcano should be ~2.5–5 × 106 times mechanically weaker than a real volcano. Using very low cohesion granular material enables simulation of volcanoes that have a bulk cohesion of 106 – 107 Pa, the approximate cohesion of fresh non-fractured rock.

Natural white quartz sand sieved to 125–180 μm grain size range was used as granular material to simulate the brittle crust and the volcano. The sand was sifted to build up the brittle crust and the cone in order to obtain a homogeneous, densely packed sand as suggested by Lohrmann et al. (2003) and Panien et al. (2006). To test the effect of cohesion in the model, experiments were also run with a 10:1 mix of the same sand with plaster (<50 μm). This plaster has a higher coefficient of internal friction (0.78 instead of 0.62) and a low cohesion, 70 Pa, in contrast to the quasi-cohesionless pure sand.

The kinematic scaling, i.e., the speed at which the analogue model deforms relative to natural cases, is poorly constrained. Deformation velocities are entirely constrained by the balance between the brittle crust and volcano load (and thus its density and height) and the viscosity of the underlying material (Delcamp et al., 2008). Natural viscosity values are ill-constrained and range over several orders of magnitude. Therefore, the deformation velocity of analogue models cannot be directly compared with natural cases, but our approach focuses here on the structures associated with the deformation, which are controlled by the system geometry rather than the deformation velocity (Delcamp et al., 2008; Le Corvec and Walter, 2009).

The physical properties of the analogue materials are given in Table 1 and compared to values expected in natural volcanoes. These experimental materials were chosen because they have been shown to be well suited to CT analysis. PDMS has a lower and well-contrasted X-ray attenuation relative to quartz sand (Colletta et al., 1991; Panien et al., 2006), enabling good discrimination of the viscous layer deformation. Analogue models have to be stable during the ~2 h that the μCT scan lasts, i.e., the remaining displacements have to be smaller than the scan resolution (80 μm). Through vertical imaging of test experiments, it was found that horizontal surface displacements were systematically below this limit 72 h after emplacement for all model types. Four similar models for each set of experimental conditions were built and left to deform for 72 h. Comparable surface deformation structures were observed in each of the repeated models, except for occurrence of asymmetric structures due to model heterogeneities, ensuring that experimental results are repeatable. The most symmetric model for each type was selected to be scanned with the μCT setup.

X-ray Micro-tomography at UGCT

The basic principle of CT imaging is based on the attenuation of X-rays through materials. This attenuation depends mainly on material density, effective atomic number, and thickness. Fault formation in sand-box models is associated with dilatation of the granular packing and thus reduction of the X-ray attenuation, as analyzed and discussed by Panien et al. (2006).

The study was carried out with the Ghent University Centre for Tomography (UGCT) micro-CT scanner (http://www.ugct.ugent.be), which is described by Masschaele et al. (2007). The high-power directional tube head was used for maximum beam intensity at a high voltage of 130 kV and a tube current of 269 μA target current, resulting in 35 W effective target power. The X-ray beam was filtered by 1.5 mm Cu and 3 mm Al to prevent beam hardening effects. A Varian PaxScan 2520 V detector with CsI scintillator was positioned horizontally to obtain an image width of 1880 pixels of 127 × 127 μm2. Due to the flat geometry of the volcano models, only part of the detector height was used. 1500 projection images were recorded; each projection was an average of 4 frames with 500 ms exposure time, resulting in a total scan time of a little less than 2 h. With a source-to-object distance of 516 mm and a source-to-detector distance of 820 mm, a magnification of 1.59 was achieved, which corresponds to a spatial resolution of 80 μm. This resolution is much better than that achieved by medical scanners, but it also has two drawbacks. The duration of the scan is much longer, making a scan of the experiment during the deformation impossible. Because the sample itself is rotating (Fig. 1), the rotation speed must be set sufficiently low to prevent damage from centrifugal effects to the model. It was verified experimentally that a rotation speed of 25° s−1 does not affect the structure. Because influence on scan time is negligible, a safety margin was included by setting the rotation speed to 10° s−1.

All data sets were reconstructed using the software package Octopus (http://www.inct.be). The server-client mode reduced the reconstruction time for around 250 slices of 1880 × 1880 pixels to 5 min for 4 clients. No beam hardening correction was used in the reconstruction. For 3D rendering, the software package VGStudioMax was used.

EXPERIMENT RESULTS

Spreading of a Low Cohesion Circular Cone

The first experiment simulates the deformation of a circular cone over a thin silicone layer (i.e., initial Thd/Hco = 0.2; see Table 1 for symbols and Figures 2–5) using a fine, low-cohesion quartz sand and plaster mix. The cone is originally 7.2 cm in diameter and 2.2 cm high (i.e., initial Hco/Wco = 0.3). After emplacement, the cone rapidly spreads radially due to the loading exerted on the silicone. After 72 h of deformation, the cone is reduced drastically in height to 1 cm above the base and its diameter increased to 7.9 cm (i.e., final Hco/Wco = 0.07).

A Digital Elevation Model (DEM) of the experiment surface is extracted from the μCT scan data. Figure 2A, C show the shaded relief and slope distribution derived from the DEM. The model surface is highly dissected and irregular. The cone is characterized by 8 well-expressed symmetric grabens, 1–3 mm deep, and several immature or asymmetric grabens, bounded incompletely or only on one side by a fault. These grabens closing at the base of the cone form a conjugate “flower pattern” typical of volcano spreading (Merle and Borgia, 1996; Delcamp et al., 2008). Conjugate graben-bounding faults form an angle of 37–43°.

In addition to V-shaped grabens, the surface of the model is characterized by a high density of fractures, not always associated with vertical offset. The central part of the cone is strongly dissected, with a blocky pattern in the summit zone surrounded by a zone with high fracture density (Figs. 2A, 2C, and Fig. 5). Note that such fracture zones have only been reported by Delcamp et al. (2008), and then only for the surface. Slopes >45° characterize the sharp graben-bounding structures and the edge of the summit blocks. The mid-flank slopes are 10–15° on average, whereas lower flanks range between 20° and 38°. Slope angles in the sliding wedges (Fig. 2B, Fig. 5) between grabens (30–37°) are significantly steeper than slopes within grabens (20–30°).

Figure 3 illustrates the 3D reconstruction of the model volume from the μCT data. It illustrates the great opportunity offered by μCT volume reconstruction to provide an infinite number of cross-sections within the 3D volume, enabling us to relate the internal structures with their topographic expression. It can be seen that the cone not only spreads radially, but also sags significantly into the silicone layer, which thins to less than 1 mm below the cone summit.

Figures 4 and 5 illustrate typical cross-sections through the model. These serial cross-sections enable the analysis of the structures that accommodated the cone deformation. Moving from the cone base toward its summit, it is observed, in horizontal cross-sections (Fig. 4A–F), that the outer half of the cone is characterized by undisturbed material except for the graben-bounding structures and some shallow fractures. The graben-bounding faults diverge slightly toward the center and have a curved-outward outline. Structures bounding two adjacent grabens merge into a single structure toward the cone center (Fig. 4A–C). At the cone base level (Fig. 4D), graben-bounding structures penetrate for 50–60% of the cone radius toward the cone center, before significant interaction occurs with adjacent structures. Moving inward into the cone, there is a 0.5–1 cm circular zone with lower density, isolating the central blocks from the outer grabens (Fig. 5). This intermediate zone, characterized at surface by low slopes (Fig. 2C) and general low X-ray attenuation, results from the intersection of the graben-bounding structures.

In vertical cross-sections (Fig. 4G–J), a significant vertical gradient of average X-ray attenuation is visible (Fig. 4G–J). This gradient is associated with variation in the total depth of material crossed by the X-rays. Although this effect limits the capability of identifying structures through the entire cone height, the main graben-bounding structures can be identified. The dip of normal faults bounding grabens is 60–65° on average (Figs. 4G, 4J, and Fig. 5C). The summit zone (i.e., ~30% of the cone radius around the center) is characterized by a complex network of structures delimiting individual compact blocks separated by open fractures and normal faults. The contrast between the dense blocks and the intermediate disturbed zone is less well marked deeper into the cone, as normal faults merge below the summit (Fig. 5B). A high density occurrence of a few millimeter-deep tension gashes and deeper fractures is imaged in the intermediate low density zone and around homogeneous summit blocks. These structures, some as narrow as 200 μm, have a sharp contrast with the surrounding material. Deeper fractures could be the prolongation of graben-bounding faults cutting into the cone center (Figs. 4–5).

Around the cone base, the deformation is associated with formation of concentric, generally asymmetric, anticlines (Fig. 3C). As illustrated by the interpretation of surface topography on Figure 5A, these fold structures display lobate shapes and appear to be related in places to graben-bouding faults through relay ramps. Sectors of cones with less well-expressed basal folding are those associated with a slightly higher initial ratio of brittle to ductile crust thickness (i.e., subtle heterogeneities in the initial model).

μCT also provides a good image of the shape of the sand-silicone contact. Generally thinner below the cone, the silicone is seen to rise into extensional zones defined by graben-structures but also below the cone center and in the circular zone of heavy fracturing below the mid-flanks (Fig. 4F, 5B–C). Silicone upwelling is also associated with an anticline structure around the cone base.

Spreading of a Cohesionless Circular Cone

Figures 2B, 2D, and 6 illustrate the topography and representative sections of an experiment with the same initial conditions as the one presented in the previous section, except for the use of a cohesionless sand for the cone and brittle crust (Table 1). Striking changes in the resulting cone morphology and structures are associated with the absence of material cohesion.

The surface of the model presents a smooth texture with twelve symmetrical V-shaped grabens on the outer flanks. The central half of the cone is mainly flat with no visible surface structure. Within this central zone, steeper slopes (10–25°) are, however, observed to surround the summit area (i.e., 1–1.5 cm around the center). The lower flanks are characterized by slopes close to the angle of repose (25–35°). No significant slope differences are observed within or outside the grabens, nor are the graben-bounding structures characterized by steeper slopes.

Looking at the model topography, the conjugate graben-bounding faults intersect at an angle of 36–42° with smaller opening angles observed on flanks with more closely-spaced grabens. Going up on the slope toward the summit, structures from adjacent grabens intersect with each other before becoming indiscernible. Grabens are ~1.2 cm wide and ~1.8 mm deep in their upper zone. From the cone base they extend ~2 cm inward (~40% of the cone radius).

In horizontal sections, grabens are bounded by sub-radial structures with lower X-ray attenuation than the surrounding material (Fig. 6A–C). These faults are systematically 1.4–1.8 mm wide (or 18–23 voxels). At high levels in the cone, they affect only the outer cone portion. At lower levels, graben-bounding structures are seen to start always at the same absolute distance from the center, affecting up to 60% of the cone radius at the base level. In horizontal sections, conjugate graben-bounding structures slightly diverge from each other in the upper half of the cone when propagating into the cone (Fig. 6B–C). Faults from adjacent grabens intersect each other within the cone isolating horsts, or sliding wedges (Fig. 6C), between successive grabens. At lower levels, the conjugate faults are joined at the cone outer edge but rapidly diverge from each other to become parallel further inward, preserving an internal graben (Fig. 6D–E).

In vertical cross-sections, each graben is seen to be bounded by two linear normal faults converging toward each other at the base of the model. These conjugate normal faults dip at an angle varying from 45 to 70°. Analysis of many successive vertical sections through a single graben shows that these faults are steeper on lower flanks (65–70°) than in the upper zone of the graben (i.e., dip ~45–55°). Cross-sections through the central part of the cone reveal no well-defined structures but illustrate the flat summit topography and the steeper outer flanks (Fig. 6I–J).

In addition to the identification of well-defined structures, μCT imaging also enables the documentation of broader zones with contrasted average X-ray attenuation and thus contrasted material density. As seen in Figure 6B–C, wedges, or horsts, and grabens between normal faults present a bright color, and thus high density, assumed to correspond to the preserved initial cone density. Preservation of the material density, and thus the compaction, suggest that these zones were not affected by fracturing or faulting outside the well-defined fault planes. In contrast, the central part of the cone displays a much lower X-ray attenuation, and thus a lower density. This zone, of constant radius with depth, underwent a significant reduction in compaction during deformation as a result of the complex interactions of the individual graben-bounding faults. As illustrated in Figure 6A and Figure 6I–J, a shallow zone at the summit preserved a higher density, suggesting that most of the decompaction associated with deformation occurred around and below the summit without significantly affecting it.

Finally, μCT provides a good contrast between the granular material and the silicone layer. The contact sand-silicone is even enhanced by a thin layer with higher X-ray attenuation than for the sand, due to infiltration of the silicone between sand grains at this contact layer (Fig. 6E–J). Figure 6G–J shows that cone spreading is associated with significant changes in the silicone layer thickness, thinning significantly below the cone and thickening around its base. This thickening is not uniform, as silicone is seen to rise along radial ridges below the flank grabens, as well as within the low density concentric zone around the summit (Fig. 2B, 6F–H). Directly at the cone base, the silicone forms 2–3 sets of concentric ridges. The first one has an amplitude of 1.5 mm and a wavelength of ~3 mm, resulting in a 0.6–1 cm broad and 1.7 mm high anticline in the brittle crust. This first anticline is asymmetric in places (Fig. 6I) verging toward the cone, with steeper surface slopes on its inward side. Folds developing further away from the cone have lower vertical amplitudes and greater wavelengths.

Sagging of a Circular Cone

Volcanoes can build up over thicker layers of ductile sediments or deform over viscous mantle material (e.g., van Wyk de Vries and Matela, 1998). Here we simulate these conditions by emplacing a volcanic cone of the same shape and size as previous models over a thicker silicone layer (Thd = 1.1 cm instead of 0.45 cm in previous cases). The silicone thickness is about half the cone height (Thd/Hco ~0.5). The brittle crust thickness is kept to half the thickness of the ductile layer. The same experiment was run with non-cohesive and low cohesion granular material.

Figure 7 illustrates the extracted surface topography of the two experiments after 72 h of deformation. Figure 8 provides a 3D view of the reconstructed experiment volume with a vertical cross-section through the cone center and its interpretation (Fig. 8C). Contrary to spreading over thin ductile layers, deformation over a thicker silicone layer causes reduction in the apparent cone radius. Height above its base level was also reduced by 35% and 20%, for the non-cohesive and low cohesion cases, respectively. Starting from the cone summit outward, steep slopes (38–50°) are observed on the upper flanks for both models. Slopes gradually decrease toward the cone base that is characterized by a topographic depression, or moat, whose depth corresponds to a third of the cone height. The surface of the non-cohesive model is smooth and displays no structures in the cone. Greater variability in slope angle is observed for the low-cohesion cones, associated with lobate surface features bordered by curvilinear steep slopes, here called terraces (see Byrne et al., 2009; Fig. 8C). The outer slopes bordering the moat are linear and sharply delineated with slopes between 25° and 45°.

Below the cone center, the silicone is totally depleted, with deformation stopping when the sand reached the base of the experiment setup (Fig. 8). The vertical sagging of the cone was accommodated by lateral flow of the silicone away from the cone. Flexure of the brittle crust caused the formation of the moat and the gradual decrease in the cone's lower slope angle. In the case of the non-cohesive experiment, the silicone layer thickened homogeneously around the cone. The only structure to be observed in the brittle crust is a concentric zone of extension, characterized by lower X-ray attenuation at the outer edge of the moat (Figs. 7A–B and 8A). A similar zone of extension directly outside the moat is observed in the experiment with low cohesion material (Figs. 7C–D and 8B), characterized by tension fractures 0.2–0.5 mm wide and 2–3 mm deep. As silicone thickened non-uniformly, a second set of tension fractures, up to 1 mm wide and cutting through the entire brittle crust (i.e., ~5.5 mm), is observed at greater distance from the cone (Figs. 7C–D and 8B). In places these two zones of extension merge into one (Fig. 8C). The reduced silicone thickness at the experiment border might be caused by the proximity of the setup edge and would not be expected for a ductile substratum with greater lateral extent.

Cross-sections through the reconstructed experiment volumes reveal no evidence for structures accommodating the vertical sagging of the cone into its substratum. A faint, but systematic variation in the average X-ray attenuation across the cone is observed on vertical and horizontal cross-sections. The outer flanks of the cone down to the moat level are characterized by lower X-ray attenuation, with higher attenuation in the center zone (Fig. 8). The width of the low density zone increases downward from the cone summit to the moat level. An additional observation is the variation of X-ray attenuation in the low cohesion experiment, along circular structures in plan view, dipping sub-parallel to the slope in vertical sections (Fig. 8B).

Finally, at the sand-silicone contact, below the lower cone slopes, a concentric zone of silicone upwelling is observed, especially well expressed in the low cohesion experiment. On selected cross-sections, this silicone upwelling is associated with a narrowing-upward structure in the sand with low X-ray attenuation, suggesting that extensional stresses are exerted at the base of the brittle layer in the zone of highest curvature (Fig. 8C).

Gravitational Deformation of an Elongated Ridge

As eruptions often occur along parallel fissures, volcanoes can develop as elongated ridges rather than circular cones. Figure 9 illustrates the cases of spreading and sagging of an non-cohesive elongated sand ridge (i.e., short to long axis ratio of 0.55) over thin and thick ductile substrata, respectively.

To simulate lateral spreading of an elongated ridge, an 11.7 × 6.3 cm wide, 1.8 cm high cone was emplaced over a 4.5 mm thick silicone layer. After deformation, the cone was 12.1 × 7.8 cm wide and 1 cm above its base, with greater spreading in the orientation of its short axis. The spreading of this elongated ridge results in a construct and structures comparable to its circular analogue, with one dominant direction of deformation perpendicular to the elongation axis (Fig. 9A). One main graben bounded on each side by main normal faults forms along the summit ridge. Within this central graben several smaller structures define nested internal grabens. At each extremity of the ridge, this central graben separates into 2 to 3 diverging V-shaped grabens. Graben-bounding normal faults dip at 45–50° in the central part of the main graben, and at 55–60° in the V-shaped grabens on the flanks (Fig. 9G–H). A few poorly expressed grabens are also observed on the long flanks. At depth into the sand ridge, the deformation results in several low density linear structures parallel to the elongation axis, causing an overall decrease in material density along the summit ridge (Fig. 9C). Silicone upwelling along two linear structures at the cone base (Fig. 9D) corresponds to the extensional zones delineated by the two main graben-bounding normal faults and the secondary internal structures (Fig. 9H).

To simulate sagging of a sand ridge, a 11.3 × 6.3 cm wide and 1.7 cm high cone was emplaced over a 11 mm silicone layer. After deformation, the cone was 11.2 × 4.8 and 1 cm above its base, with thus greater apparent basal shortening along the short axis. Vertical cross-sections parallel to the cone short axis present similar characteristics to the experiment of circular cone sagging presented in the previous section (Fig. 9I–J). The slope decreases toward the cone base, and the topographic moat has a depth corresponding to half the remaining cone height. The cone sagged vertically to reach the base of the setup and the silicone is thicker than initially below the flat brittle crust. Concentric zones of extension, ~1 cm wide at the surface and narrowing with depth, defined by their lower X-ray attenuation, are observed adjacent to the moat, where the surface becomes suddenly flatter, but also close to the setup edge probably due to border effects (Fig. 9B, E, I, J). Figure 9F also illustrates that the lower cone flanks are characterized by lower material density than the material in the central axis of the cone. Just above the sand-silicone contact level, a narrow structure of low X-ray attenuation, parallel to the cone long axis, can be observed at the vertical of the moat's deepest point. This structure is interpreted as an extensional structure associated with low amplitude silicone upwelling. Parallel to the long axis, the vertical sagging of the flanks is more gradual and the moat is less deep. No significant variation in material density along the long axis can be observed.

DISCUSSION

3D Gravitational Deformation of Volcanic Cones

Detailed comparison of the experimental results presented here with specific natural volcanoes is outside the scope of this preliminary study. Field evidence for gravitational deformation of natural volcanoes has been provided in numerous previous studies (e.g., Borgia et al., 1992; Borgia, 1994; Merle and Borgia, 1996; van Wyk de Vries et al., 2001; Borgia and van Wyk de Vries, 2003; Morgan and McGovern, 2005; Münn et al., 2006; Walter et al., 2006; Delcamp et al., 2008; Le Corvec and Walter, 2009). μCT imaging of the presented analogue experiments, however, provides new insights into the gravitational deformation of volcanic constructs.

In the case of a spreading cone, key new information is the vertical and horizontal variations of cone density and especially the formation of a central brecciated zone of sub-polygonal shape (Fig. 10). This zone is delimited by the prolongation of the individual graben-bounding faults. Within this central zone, the central part is characterized by undisturbed blocks separated by clear fractures, whereas the outer part of the polygon is more heavily fractured. This variation in fracture density cannot be highlighted from analysis of the model topography nor preserved in physical sections of the model. It is possible to preserve the fracture pattern at depth through the introduction of several thin plaster layers within the cone structure, excavating the cone until reaching each successive layer (Lebas, 2009). Observations of larger spreading cones using such a technique provided a similar fracture network as documented in 3D by μCT (Fig. 11). However, direct excavation is difficult, intrusive, destroys structures, and does not give a clear 3D picture. The scans here provide hitherto inaccessible detail of the fractured core of the cone. This weaker zone, created by the complex interactions of individual graben-bounding structures and the vertical sinking of the volcano summit, will affect how magma intrudes into a volcanic cone, control the structure of hydrothermal systems and ore deposition, and influence the destabilization of the central portion of the edifice.

Delcamp et al. (2008) have shown that graben-bounding structures are transtensional strike-slip faults and not pure normal faults. μCT imaging does not enable the detail constraint of the motion of structures as normal and strike-slip faulting both result in similar material decompaction (Colletta et al., 1991; Ueta et al., 2000). It is, however, worth noting that movements along graben-bounding structures are accommodated at the base by fold or thrust structures (Fig. 10A), and that the two systems of faults are related to each other through relay ramps that enable transfer from transtentional to compressive movement. Delcamp et al. (2008) observed that graben-bounding faults can prolong into the surrounding base as sub-radial transtensional strike-slip faults (Fig. 10B). This situation is characteristic for high values of the ratio of brittle to ductile layers’ thickness and for a gradual slope decrease at the cone base. Comparing idealized spreading models of Figures 10A and B, it can be seen that a continuum exists between these two cases. Relay ramps transform into a compressive structure when they cross at high angles, whereas when they cross each other at a smaller angle, they delineate diamond-shaped sliding wedges.

The sliding wedges that formed between grabens in both models are not deformed and not affected by fracturing, as highlighted by the homogeneous high density in μCT scans. These more preserved zones will be the volumes that will resist erosion at natural volcanoes and will thus form the part of the volcano that will strike out after prolonged erosion periods. Such topographic wedges will provide good evidence of spreading after significant erosion. This is the case, for example, of many eroded strato-volcanoes in the Andes or Kamtchatka (see figure 15 in Delcamp et al., 2008).

In the case of sagging cones, the high resolution DEM of the model enables the identification of subtle terrace features. Such structures have been observed to occur at large natural volcanoes undergoing basal flexure on Earth and Mars (Borgia et al., 1992; Merle and Borgia, 1996; Byrne et al., 2009). These terraces do appear to be related to significant deformation at depth. The density stratification observed in the central cone zone (Fig. 8) is probably to be attributed to grain size segregation through avalanching during cone emplacement and formation of thin, plaster rich layers in the central part of the cone. The decreased X-ray attenuation on the lower flanks of sagging cones is probably due to grain rearrangement associated with the terrace-like thrusts seen on the surface. The scans show that these structures are formed without clear localization along individual structures below the surface. The decompaction is somewhat counter-intuitive as it is in association with an expected contraction in this area (van Wyk de Vries and Matela, 1998), but Colletta et al. (1991) showed that lower X-ray attenuation also characterizes compressive structures. Here, radial and concentric contraction must be accommodated by vertical dilation that appears at the surface as terraces.

μCT provides an image of the model surface as well as the cone base after flexure. This is impossible to do through model section, as the silicone moves rapidly. Comparison of the flexure surface and the surface topography with the original model conditions enables the determination of the strain in the volcano. Extension can, for example, be inferred at the lower boundary of the brittle crust for sagging cones, as illustrated by silicone upwelling and structures of lower X-ray attenuation (Fig. 7C), whereas contraction takes place on the cone's lower flanks on the terraces.

Tension fractures and extensional fractures are exceptionally well imaged by the μCT technique. These features are usually not preserved when physically cutting the model. Extension structures are imaged as tension gashes or narrow lower attenuation zones in the case of low cohesion material, whereas they appear as broader zones of lower attenuation for non-cohesive models (see, for example, extensional structures in the crust around sagging cones, Figures 7–8). These features are essential for indicating directions of extension.

μCT imaging has been proven to provide valuable information about volcano deformation processes through detailed reconstruction of the volcano volume, which enables documentation of the relationships between the surface topography, the structures at depth, and the silicone deformation. The relationships between individual structures; the documentation of the 3D variation in bulk density; and the possibility to infer strain distribution at the surface, at depth, and at the brittle-ductile boundary are other contributions of this technique.

Advantages and Limitations of μCT

First, volcano topography is documented in great detail. The method presented here has the advantage of providing a higher resolution (~80 μm) than that achieved with photogrammetry techniques (200–300 μm; Cecchi et al., 2003; Donnadieu et al., 2003), to be independent of the surface texture and to enable detailed characterization of deep but narrow fractures, usually affected by shadow effects in stereoscopic methods. Characterization of surface topography and slope distribution is essential for comparison with natural volcanoes for which the morphology and surface structures are often the only available evidence for active deformation.

The second added value of μCT is that it enables imaging of internal structures. Horizontal and vertical cross-sections through the models display clear variations in the X-ray attenuation. As discussed by Colletta et al. (1991), X-ray attenuation in granular material depends on the mineral composition of the grains, the grain size, and the bulk density or the compaction of the material. As the material can be assumed to be initially homogeneous, variation in the observed X-ray attenuation after deformation can be attributed to variation in compaction due to rearrangement of grain packing during deformation. In spreading experiments, surface traces of graben-bounding faults always correspond at depth with thin planar structures of lower attenuation.

Striking differences are observed between spreading of a non-cohesive and a low cohesion sand cone. In experiments with sand-plaster mix, the surface topography is dissected, tension fractures are preserved, and deformation structures are sharp at the surface and thin in cross-sections. The X-ray attenuation of faults is well contrasted with the surrounding, undisturbed material. With non-cohesive sand, structures are broader and more diffuse, surface topography is flatter, and individual structures cannot be discriminated in areas of high fault density. For pure quartz sand, the fault width corresponds to ten times the average grain size (~150 μm). This matches with previous observations (Panien et al., 2006 and references therein) and theoretical modeling of shear localization in granular media (Herrmann, 2001). Mixing the sand with a limited fraction of fine cohesive material such as plaster enables the width reduction of these shear zones by a factor of two.

The application of μCT imaging to volcanic analogue models is also constrained by several limitations. First, the spatial resolution of the μCT scan is directly proportional to the size of the model. In order to enable identification of structures smaller than a millimeter, it is necessary to develop analogue models of the order of 10 cm in diameter or smaller, 2–5 times smaller than the ones usually used in analogue models (Merle and Borgia, 1996; Delcamp et al., 2008). To prevent border effects, cone deformation must not approach the model edges too closely. This further reduces the size of the investigated cone. Building homogeneous models with well-controlled crust thickness and constant cone slope becomes a challenge at such a small size.

The signal to noise ratio of μCT data depends partly on the duration of the scan but also on the thickness of the material that the X-ray has to pass through before reaching the X-ray detector. In our experiment, a significant contrast in average X-ray attenuation is observed between the crust zone and the cone, as the material thickness varies by a factor of two. This variation in average attenuation complicates the interpretation of vertical cross-sections in the reconstructed experiments. The lower signal to noise ratio within the crust prevents clear identification of all the structures to their lowest level. This effect also causes higher noise in zones of topographic lows, such as the moat surrounding the sagging cones.

In contrast to medical scanners, which are using a linear detector and higher energy, μCT scans require 30–120 min to image an object; longer scan times enable a higher signal to noise ratio. During the scan, the object is slowly rotated but it needs to be fully stable relative to its rotation axis for accurate 3D reconstruction. This means that the experiments on volcano gravitational deformation could only be imaged after a prolonged deformation period (i.e., 72 h) once continued deformation became negligible. In the simulated gravitational deformation processes, there is no way to stop the process before it reaches a stable situation. The magnitude of deformation can only be reduced by increasing the thickness of the brittle crust relative to the ductile substratum (Merle and Borgia, 1996). This limitation prevents imaging the model at its initial stage and at intermediate stages during deformation, in order to document the vertical and horizontal surface deformation, as with photogrammetry techniques (Donnadieu et al., 2003; Delcamp et al., 2008), or to document the spatio-temporal evolution of deformation structures, as in 2D experiments (Burchardt and Walter, 2009; Le Corvec and Walter, 2009). However, this issue will not limit experiments simulating most other types of volcano-tectonic processes in which deformation (e.g., magma intrusion, magma chamber inflation/deflation, tectonic fault movement) can be implemented incrementally and stopped at different stages. Such experiments would then be imaged at their initial stage and at several stages of deformation, as is done for tectonic analogue models (e.g., Schreurs et al., 2001). Layers of sands with different X-ray attenuation, but comparable friction angle and grain size, or individual denser grains can then be used as markers to document 3D movements at the surface and at depth (e.g., Wilkerson et al., 1992).

CONCLUSIONS

This study aimed at illustrating the capabilities offered by μCT in documenting analogue models of volcano-tectonic processes through a small set of experiments focusing on gravitational deformation of a volcanic construct over a ductile substratum. From the presented observations, we can conclude that:

  1. μCT imaging enables the reconstruction and visualization of the 3D volume of analogue models with a resolution of 80 μm.

  2. Cone spreading is associated with the formation of flank grabens and the formation of anticlines at the cone base, but the scans also show heavy fracturing of the core of the volcanic cone around a relatively preserved summit zone and upwelling of the ductile substratum in the extensional grabens.

  3. A highly accurate DEM of the experiment surface can be extracted from the 3D reconstruction, enabling quantitative documentation of the model topography and the slope distribution.

  4. Virtual cross-sections through the reconstructed model volume with vertical and horizontal orientations make it possible to document variations in X-ray attenuation associated with formation of fractures and faults and with the transition from sand to silicone material.

  5. In comparison with non-cohesive sand, experiments with low cohesion sand display sharper fault structures, preservation of shallow fractures, and formation of fractured blocks in the cone center.

  6. Documentation of analogue models through μCT is constrained by the size of the experiment, the variation of signal to noise ratio with the depth of material crossed by X-rays, and the need to have stable models during the duration of the scan.

Despite the highlighted limitations, μCT has the potential to document the full 3D complexity of analogue models essential for investigation of other volcano-tectonic processes, such as caldera formation, magma emplacement, or volcano interactions with tectonic structures. This technique will prove to be especially powerful for processes whose deformation can be simulated incrementally and repetitively imaged by μCT. Documenting the 3D complexity of deformation structures offers the potential to more closely relate analogue observations with structures accommodating deformation at natural volcanoes as documented by geophysical or geological evidence.

M.K., V.C., and K.F. are supported by the Research Foundation–Flanders (Fonds voor Wetenschappelijk Onderzoek–Vlaanderen).