B. van Wyk de Vries

Sub-surface structures and collapse mechanisms of summit pit craters

Abstract Summit pit craters are found in many types of volcanoes and are generally thought to be the product of collapse into an underpressured reservoir caused by magma withdrawal. We investigate the mechanisms and structures associated with summit pit crater formation by scaled analogue experiments and make comparisons with natural examples. Models use a sand plaster mixture as analogue rock over a cylinder of silicone simulating an underpressured magma reservoir. Experiments are carried out using different roof aspect ratios (roof thickness/roof width) of 0.2–2. They reveal two basic collapse mechanisms, dependant on the roof aspect ratio. One occurs at low aspect ratios (≤1), as illustrated by aspect ratios of 0.2 and 1. Outward dipping reverse faults initiated at the silicone margins propagates through the entire roof thickness and cause subsidence of a coherent block. Collapse along the reverse faults is accommodated by marginal flexure of the block and tension fractures at the surface (aspect ratio of 0.2) or by the creation of inward dipping normal faults delimiting a terrace (aspect ratio of 1). At an aspect ratio of 1, overhanging pit walls are the surface expressions of the reverse faults. Experiments at high aspect ratio (>1.2) reveal a second mechanism. In this case, collapse occurs by stopping, which propagates upwards by a complex pattern of both reverse faults and tension fractures. The initial underground collapse is restricted to a zone above the reservoir and creates a cavity with a stable roof above it. An intermediate mechanism occurs at aspect ratios of 1.1–1.2. In this case, stopping leads to the formation of a cavity with a thin and unstable roof, which collapses suddenly. The newly formed depression then exhibits overhanging walls. Surface morphology and structure of natural examples, such as the summit pit craters at Masaya Volcano, Nicaragua, have many of the features created in the models, indicating that the internal structural geometry of experiments can be applied to real examples. In particular, the surface area and depth of the underpressured reservoir can be roughly estimated. We present a morphological analysis of summit pit craters at volcanoes such as Kilimanjaro (Tanzania), San Cristobal, Telica and Masaya (Nicaragua), and Ubinas (Peru), and indicate a likely type of subsidence and possible position of the former magma reservoir responsible for collapse in each case.

Hummocks: how they form and how they evolve in rockslide-debris avalanches

Hummocks are topographic features of large landslides and rockslide-debris avalanches common in volcanic settings. We use scaled analog models to study hummock formation and explore their importance in understanding landslide kinematics and dynamics. The models are designed to replicate large-scale volcanic collapses but are relevant also to non-volcanic settings. We characterize hummocks in terms of their evolution, spatial distribution, and internal structure from slide initiation to final arrest. Hummocks initially form by extensional faulting as a landslide begins to move. During motion, individual large blocks develop and spread, creating an initial distribution, with small hummocks at the landslide front and larger ones at the back. As the mass spreads, hummocks can get wider but may decrease in height, break up, or merge to form bigger and long anticlinal hummocks when confined. Hummock size depends on their position in the initial mass, modified by subsequent breakup or coalescence. A hummock has normal faults that flatten into low-angle detachments and merge with a basal shear zone. In areas of transverse movement within a landslide, elongate hummocks develop between strike–slip flower structures. All the model structures are consistent with field observations and suggest a general brittle-slide emplacement for most landslide avalanches. Absence of hummocks and fault-like features in the deposit may imply a more fluidal flow of emplacement or very low cohesion of lithologies. Hummocks can be used as kinematic indicators to indicate landslide evolution and reconstruct initial failures and provide a framework with which to study emplacement dynamics.

A sagging-spreading continuum of large volcano structure

Gravitational deformation strongly infl uences the structure and eruptive behavior of large volcanoes. Using scaled analog models, we characterize a range of structural architectures produced by volcano sagging and volcano spreading. These arise from the interplay of variable basement rigidity and volcano-basement (de-)coupling. From comparison to volcanoes on Earth (La Reunion and Hawaii) and Mars (Elysium and Olympus Montes), the models highlight a structural continuum in which large volcanoes throughout the Solar System lie.

N-view reconstruction: a new method for morphological modelling and deformation measurement in volcanology

Abstract We present a method of reconstructing volcanic morphology using multiple digital views ( N -view), captured at different angles around an object. This approach uses recent advances in the field of Computer Vision to provide accurate 3-D measurements of volcanic surfaces. Videogrammetry (digital image reconstruction) is used, as it is best adapted to numerical processing. The method is tested and now used in the laboratory on analogue volcanic cones. The method begins with calibrating the camera and finding image positions, using an accurate N -view calibration method. This is done by estimating internal and external parameters of the camera using several views of a specially constructed calibration target. The N -view reconstruction of the real object is then done by iteratively deforming an initial theoretical model of the surface. Laboratory tests show that reconstruction accuracy is about 10 −4 m for a distance between the object and the camera of 0.5 m, and is potentially several orders of magnitude higher for surfaces of finer texture and using higher precision sensors. This is easily high enough to be useful for the accuracy required for morphological studies. It is also sufficient for monitoring most types of volcano deformation. The technique has the potential to detect morphology changes of the order of mm. Use of the method in the field requires a slightly different approach from that in the laboratory: textures and lighting are more variable, and the sensor and ground control point location and model calibration must be done differently. We provide case studies from laboratory tests and qualitative image analysis for two field cases: Piton de la Fournaise (Indian Ocean) and Santaguito (Guatemala). These illustrate the technique’s potential and explore problems of field application. Using current sensors, the method has the potential to provide sufficient precision for fine (mm–cm) scale reconstruction, and will represent a valuable, simple and flexible tool when compared with classical stereophotogrammetry techniques.

Digital photogrammetry as a tool in analogue modelling: applications to volcano instability

Abstract Three techniques of digital photogrammetry have been applied successfully to laboratory analogue models to study surface displacements caused by various volcano deformation types. Firstly, side-perspective videos are used to differentiate profile displacements between cryptodome intrusion models and models deforming by ductile inner-core viscous flow. Both models show similar morphologic features including a bulged flank and an asymmetric upper graben. However, differences in displacement trajectories of the bulge crest reflect upward intrusion push contrasting with essentially downward displacement vectors of weak core models. The other two techniques use vertical views correlated automatically either as time-sequence monoscopic views or as coeval stereoscopic pairs. This exploits to a maximum the method’s potential by imaging surface displacements over the whole model. Successive monoscopic photograms, because they suffer only moderate numerical processing for topographic effect removal, can detect very small displacements occurring early in deformation processes. As illustrated by analysis of intrusion models, the monoscopic method allows prediction of fault locations and main displacement locations. It can also anticipate the principal strain directions, and separate different deformation stages. On the other hand, the stereo-photogrammetry technique, although more complicated, provides topography and volume changes, as well as pictures of surface displacements in three dimensions. Results are presented for the spreading of volcano models on a ductile substratum and viscous cored cones. We have found digital photogrammetry to be a useful tool for analogue modelling, because it provides quantitative data on surface displacements, including movement invisible to the eye, as well as topographic changes. It is a good method for investigating and comparing different deformation mechanisms. It is especially useful for interpretation of displacement patterns obtained from monitoring of natural active volcanoes. In fact, results of the methods used in the laboratory can be directly compared with field data from geodetic or photogrammetric surveys, as at Mount St. Helens in 1980.

Craters of elevation revisited: forced-folds, bulging and uplift of volcanoes

The eighteenth/nineteenth century ‘craters of elevation’ theory required magma to uplift strata, doming the surface and creating a central down-fallen ‘crater’ or graben. Exponents of craters of elevation attempted to apply it to explain the origin of all volcanoes, and rapidly the contemporary competing ‘craters of eruption’ theory replaced it as the paradigm for volcano construction. Several historic examples have shown that intrusions can cause uplift, termed bulges and can create features like those proposed for craters of elevation (e.g. at Usu 1944, Bezymianny 1955 and Mt. St. Helens 1980). Work on sedimentary basins that have had igneous activity has shown that intrusions create ‘forced folds’ that uplift and deform strata in a similar way to that originally proposed for craters of elevation. In view of the above, we investigate large-scale intrusion-related topographic changes at two sites where the craters of elevation theory was developed: the monogenetic volcanoes of the Chaîne des Puys, France and the Teide stratovolcano, Tenerife. We combine observations of such features with examples of forced folding to integrate the two fields of research. Our observations in the Chaîne des Puys show that: (1) the Petit Puy de Dôme has a bulge of up to 150-m uplift. The uplift has a central depressed area (a graben), a dense network of normal faults, basal thrusts and an aborted landslide. (2) The Grosmanaux volcano is a forced fold created by uplift of a previously flat-lying area, and has dense faulting and a graben on the resultant topographic bulge. It was the site also of a major vulcanian eruption from the associated Kilian crater. (3) The Gouttes volcano was uplifted by an intrusion like the Petit Puy de Dôme, but then collapsed to generate a landslide and lateral blast. (4) Excavation in the Lemptégy Volcano exposes intra-eruption intrusions with associated uplift, providing examples in cross-section of the internal deformation likely to be found inside other Chaîne des Puys uplifted bulges. On Teide, a bulge near the summit shows similar structures and surface tilting as seen on the Petit Puy de Dôme and this bulging may have formed during the eruption of the Lavas Negras, the most recent activity on the summit area. Fault scarps on Teide also expose small cryptodomes, like those seen at Lemptégy. These examples, coupled with field studies on eroded intrusions, data on forced folds in basins and analogue models, show how large-scale topographic remodelling and structural change can be created by intrusions. These can rapidly and significantly change the volcanic edifice. A crater of elevation bulge, or forced fold that is stabilised by the cooling of the intrusion, will remain an important structural element in a volcano. This process starts even at the small scale of monogenetic volcanoes, and could occur through the lifetime of any growing stratovolcano. Such activity may be commonplace, but may be masked by concomitant eruption or removed by subsequent collapse. Monitoring and hazard strategies should be ready to deal with such large-scale events that will seriously modify the eruptive activity and stability of a volcano within days or weeks.

Dykes and structures of the NE rift of Tenerife, Canary Islands: a record of stabilisation and destabilisation of ocean island rift zones

Many oceanic island rift zones are associated with lateral sector collapses, and several models have been proposed to explain this link. The North–East Rift Zone (NERZ) of Tenerife Island, Spain offers an opportunity to explore this relationship, as three successive collapses are located on both sides of the rift. We have carried out a systematic and detailed mapping campaign on the rift zone, including analysis of about 400 dykes. We recorded dyke morphology, thickness, composition, internal textural features and orientation to provide a catalogue of the characteristics of rift zone dykes. Dykes were intruded along the rift, but also radiate from several nodes along the rift and form en échelon sets along the walls of collapse scars. A striking characteristic of the dykes along the collapse scars is that they dip away from rift or embayment axes and are oblique to the collapse walls. This dyke pattern is consistent with the lateral spreading of the sectors long before the collapse events. The slump sides would create the necessary strike-slip movement to promote en échelon dyke patterns. The spreading flank would probably involve a basal decollement. Lateral flank spreading could have been generated by the intense intrusive activity along the rift but sectorial spreading in turn focused intrusive activity and allowed the development of deep intra-volcanic intrusive complexes. With continued magma supply, spreading caused temporary stabilisation of the rift by reducing slopes and relaxing stress. However, as magmatic intrusion persisted, a critical point was reached, beyond which further intrusion led to large-scale flank failure and sector collapse. During the early stages of growth, the rift could have been influenced by regional stress/strain fields and by pre-existing oceanic structures, but its later and mature development probably depended largely on the local volcanic and magmatic stress/strain fields that are effectively controlled by the rift zone growth, the intrusive complex development, the flank creep, the speed of flank deformation and the associated changes in topography. Using different approaches, a similar rift evolution has been proposed in volcanic oceanic islands elsewhere, showing that this model likely reflects a general and widespread process. This study, however, shows that the idea that dykes orient simply parallel to the rift or to the collapse scar walls is too simple; instead, a dynamic interplay between external factors (e.g. collapse, erosion) and internal forces (e.g. intrusions) is envisaged. This model thus provides a geological framework to understand the evolution of the NERZ and may help to predict developments in similar oceanic volcanoes elsewhere.

Remote sensing of the 1998 mudflow at Casita volcano, Nicaragua

A devastating lahar (volcanic mudflow) occurred at Casita volcano (Nicaragua) in 1998, triggered by excessive precipitation associated with Hurricane Mitch. We investigate here the morphology and drainage structure of the flow deposition area, primarily using satellite optical and radar imagery. Because the lahar destroyed several towns and villages, killing over 2500 people, we also assess the utility of images available at the time of the event for disaster management. We find that SPOT multispectral data are most suited to delineate and characterize the flow field, but show limitations for damage assessment, and problems with cloud contamination. ERS Synthetic Aperture Radar (SAR) imagery largely failed to detect the lahar deposits, with RADARSAT performing slightly better. The relatively coarse-grained deposits, together with a dense cover of wooden debris, made the flow nearly imperceptible to the short wavelength C-Band radar, especially at ERSs steep incident angle. Because of spatial, spectral and radiometric limitations in all types of imagery used, the synergistic potential of optical and radar, as well as high- and low-resolution optical data was explored. The best synergy was found not in merged imagery, but in the incorporation of auxiliary information, such as elevation, map and GIS data. Given the cloud problem, forthcoming radar satellites with variable polarization and incident angles are expected to provide better results in comparable, future situations.

3D imaging of volcano gravitational deformation by computerized X-ray micro-tomography

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.

Directional flank spreading at Mount Cameroon volcano: Evidence from analogue modeling

Mount Cameroon is characterized by an elongated summit plateau, steep flanks, and topographic terraces around its base. Although some of these features can be accounted for by intrusion-induced deformation, we here focus on the contribution of edifice-scale gravitational spreading in the structure of Mount Cameroon. We review the existing geological and geophysical data and morphostructural features of Mount Cameroon and surrounding sedimentary basins. Volcanic ridge gravitational spreading is then simulated by scaled analogue models on which fault formation is recorded using digital image correlation. Three sets of models are presented (i) models recorded in cross section (Type I), (ii) models recorded from above with a uniform (Type IIa), and (iii) nonuniform ductile layer (Type IIb). Type I models illustrate the formation of faults accommodating summit subsidence and lower flank spreading. Type IIa models favor displacement perpendicular to the long axis, with formation of a summit graben and basal folds, but fail to reproduce the steep flanks. Type IIb models investigate the effect of spatial variations in sediment thickness and/or properties consistent with geological evidence. Directional spreading of the volcanos central part perpendicular to the long axis is accounted for by a sediment layer with restricted lateral extent and increasing thickness away from the volcano axis. The later model closely reproduces key features observed at Mount Cameroon: steep upper flanks are accounted for by enhanced lateral spreading of the lower flanks relative to the summit. The relevance of these findings for understanding flank instabilities at large oceanic volcanoes is finally highlighted.