Morpho-Structure of the 1982 Lava Flow Field at Mount Cameroon Volcano, West-Central Africa

M. N. Wantim, M. Kervyn, G. G. J. Ernst, M.-A. Del Marmol, C. E. Suh, P. Jacobs

Department of Geography and Earth System Science, Vrije Universiteit Brussel, Brussels, Belgium.
Department of Geology and Soil Science, Ghent University, Ghent, Belgium.
Department of Geology, University of Buea, South West Region, Buea, Cameroon.

DOI: 10.4236/ijg.2013.43052                

Abstract

Basaltic eruptions have been observed to produce structurally complex, compound 'a'ā lava flow fields but their morphometry has only rarely been systematically documented. We document the morphology and structures that developed during the emplacement of the 1982 basaltic lava flow field at Mount Cameroon (MC) volcano over a period of one month. Topographic cross-sections (13 in total) were made from the main vent (~2700 m above sea level (a.s.l)) down to a distance of 5.5 km on the cooled lava surface. Details obtained from these cross-sections include: channel width and depth, levee slope, lava surface morphology and structures. These details enabled us to describe the physical characteristics of the 1982 lava flow field. The inclined (12° - 19°) underlying slopes on which this flow field was emplaced resulted in a characteristic channelized basaltic 'a'ā flow field morphology. This includes a proximal zone characterised by reduced flow width and depth with no subsidiary channels. Slab-crusted lava dominates the proximal channel distinctively bent into convex upward shapes. 7 secondary vents were observed for the first time ~2.5 km from the main vent, with heights of 3 - 15 m. This is a very significant observation since it points to the fact that the flow field emplacement may have been a product of 2 eruption sites as observed at other historical MC lava flow fields. This supposition was ruled out by further evidence obtained from other surface features within the flow field. The presence of these secondary vents still has an important bearing in lava flow hazard assessment. Field observations also revealed the presence of tumulus. This is a novel feature for MC lava flow fields. It displayed a close similarity to those observed at other basaltic volcanoes occurring in association with clinker 'a'ā lava, lava tubes, squeeze-ups and pressure ridges. Channels are well-defined, bounded by levees. Accretional and overflow levees dominate in this flow field. This lava flow-field attained a final length of 7.5 km, an area of 2.6 × 10⁶ m² and volume of 1.3 × 10⁷ m³. The presence of tumulus indicates internal inflation together with structures such as pressure ridges and squeeze-ups which are also attributed to compressive forces. Our observations suggest that real-time monitoring of compound lava flow fields evolution at MC may reveal the emplacement mechanisms of complex structures such as the secondary vents (~2180 - 2011 m a.s.l.) observed within the flow field. In addition, documenting the occurrence, morphology and link between lava tubes, tumulus and squeeze-ups may allow us to determine the risk of reactivation of a stalled flow front. This will thereby enhance the ability to track and assess hazards posed by lava flow emplacement from MC-like volcanoes.

Keywords

1982 Eruption; Channels; Levees; Secondary Vents; Slope; Tumulus

1. Introduction

The advance of lava flows produced by volcanic eruptions has been studied through field observations [1-12]; remote sensing techniques [13-17], as well as through analytical or numerical modelling [18-23]. The above cited methods have greatly improved understanding of the emplacement dynamics of lava flows as they provide clues about key processes occurring during eruptions even for flows not witnessed. However, the physical volcanology of most lava flow fields in developing countries has received limited systematic attention.

Lava morphological data are significant to understand lava emplacement mechanism and anticipate impacts from effusive eruptions. Lava flows can cover long distances in 3 - 6 days, damaging properties and threatening lives [2,24,25]. The complex morphologies characteristic of long duration basaltic 'a'ā lava flow fields reflect the importance of processes such as inflation, formation of lava tubes and secondary vents [17]. The identification of secondary vent sites of lava flow fields is critical for lava flow hazard assessment [17]. They have been observed to allow flows to lengthen significantly over their forecasted cooling-limited lengths [26]. In contrast, overflow events or new lobes development favour flow field widening [27] instead of flow lengthening.

Advances in lava flow modelling in the past 25 years show that it is increasingly possible to anticipate the final lengths and spatial spread of single flow units [18,28]. However, during basaltic eruptions, complexities like new flow units are emplaced alongside and on top of earlier units resulting in compound flow fields [5,6]; channels crust over and lava tubes, secondary vents, squeeze-ups and tumulus develop. The development of these features hinders accurate modelling of long-lived (>1 week) basaltic 'a'ā lava flow fields [17] which is not accounted for in these models. To improve understanding on the development of compound lava flow fields and the dynamics of any volcanic system, it is recommended to quantify lava morpho-structures in detail both during and after emplacement [6,29].

Studies of terrestrial lavas suggest that the overall development of flow fields is systematic and that a general, normalized relationship can be established linking the final dimensions of a flow field to underlying slope, eruption duration, discharge rate, gravitational acceleration, lava density and rheology [27,30,31]. Linking qualitative and quantitative measurements of lava flow surface morphology with historical observations of eruptions is an important, but yet underexploited route, to constrain emplacement mechanism of basaltic lavas.

Basaltic lava flow fields often demonstrate compound morphology. That is they are comprised of several flow units and lobes with some superposed on each other [5, 17,32]). A flow lobe here refers to an individual package of lava surrounded by a chilled crust [1]. Compound here refers to a lava flow field made up of two or more flow lobes of any geometry or size [1,33]. The development of compound basaltic 'a'ā lava flow fields is a common phenomena at Mount Cameroon (MC) volcano as exemplified in the 1959, 1982, and 1999 lava flow fields [2,8,34]. Most published data existing at MC (Figure 1) are limited to its seismicity, eruption dynamics, petrology and geochemistry of specific eruptions that occurred in the 20^th and 21^st centuries [2,24,25,34-40]. Recently, a few studies started focusing on lava flow morphology [2,8,11] and hazards [35,41-43] at MC.

Eruptions that produced basaltic lava flow fields from MC have been associated with significant impacts over the past years. The 1922, 1959 and 1999 eruptions posed major threats to agro-allied complexes and road infrastructure around the SW [24,25] and NE [2] flanks (Figure 1). Common phenomena observed for most historical eruptions at MC are that they are all fissure eruptions with venting typically at more than one site (1922, 1959, 1999, 2000; Figure 2; [2,8,11,24,25,40]). Upper sites are dominated by more explosive activity building large pyroclastic cones, whereas lower sites emit the largest volume of lava. The principal lava type produced ranges from pahoehoe (ultra-proximal) to 'a'ā (dominant flow morphology) and blocky lava observed at the distal end [11]. Petrographically, all these lavas are porphyritic [2,25,39] with the exception of the 1982 lava that is nearly aphyric [2,25,34].

In October 1982, an eruption that produced seismic swarms widely felt around Buea (SE flank; Figure 2), occurred on the SW flank of MC. The physical characteristics, composition and evolution of the 1982 lava flow field have been described by [34]. Other aspects including its petrology and geochemistry have been analysed by [2,25,39]. So far, no systematic field survey has been carried out on the 1982 flow field to document its morphology or quantify down channel geometry. Here, a systematic field study was carried out on the 1982 flow field with the primary objective of obtaining geometric parameters (width and depth) of the stable channel zone that fed the flow from source to the distal transitional-dispersed flow zone.

The goal of this paper is to describe the morphology and structures observed within the stable and transitional channel zones of the 1982 lava flow field as defined by [44]. Based on their observations of basaltic 'a'ā flow fields at Mauna Loa, Hawaii, [44] divided lava flow fields into four zones from vent to toe: 1) stable channel; 2) transitional channel; 3) dispersed flow zone and 4) flow toe. We present observations and field measurements made after the emplacement of this lava flow field, covering the proximal part fed principally from the source scoria cone (emplaced on characteristic slopes of 14˚ - 18˚) and the transitional zone (12˚ - 19˚).

The 1982-lava flow field offers an excellent opportunity to examine the large-scale structural evolution of a compound flow field as well as the complex surface features that are ubiquitous in most 'a'ā flows. Such intermediate-long lived eruptions emplaced on inclined slopes offer good opportunities to investigate the effect of slope on 'a'ā lava flow fields in terms of surface morphology and inferred emplacement. So far few studies have described long-lived basaltic lava flow field morphologies and emplacement processes on inclined [11,31] and extremely steep slopes [45]. Numerous studies however exist on the morphology of long-lived basaltic lava flow fields emplaced on flat to gently sloping surfaces [27], [46]. From this analysis, we document in detail the latestage development of this flow field which includes the formation of the secondary vents (short-lived ephemeral vents), lava tubes and tumulus and consider implications for lava flow hazard assessment. We also consider the role of effusion rate fluctuations for channel geometry fluctuation. Field measurements were used to constrain flow rheology. Such constraints are essential for a follow-up effort where lava flow hazard will be quantified.

2. The Geology of Mount Cameroon and the 1982 Eruption Chronology

2.1. Geological Context of Mount Cameroon (MC)

MC is the largest and most active of the continental volcanoes of the Cameroon Volcanic Line (CVL; Figure 1). The CVL is a major tectonic feature in West-Central Africa that runs SW-NE following a major left-lateral fault system that extends for more than 2000 km, from Pagalu Island into West-Central Africa (Figure 1). The origin, geology, structure and petrology of the CVL have been discussed by several authors: [35,47-51]. The CVL’s origin is still subjected to debate. The most widely accepted structural explanation for the origin of the CVL is that it is a product of Cretaceous reactivation of PanAfrican strike-slip faults trending N70E [38,48,52].

MC is a steep volcanic shield covered by successions of lava flows [11,53] and subsidiary scoriae deposits (Figure 2). It has a flat summit plateau, a rift zone defined by a linear cluster of eruptive vents, a deep valley (elephant valley) in the N flank, topographic steps at the base [53] together with numerous faults and fissures mostly trending N40E [35]. It is composed entirely of moderately alkaline basic lavas (alkaline basalts, hawaiites, picrites and mugearites [35,48,49]. Field observations indicate that lava flows are massive, porphyritic or vesicular. Moderately explosive Strombolian activity is common at high elevation vents [11,24,25], but lava flows are the most common products released from MC eruptions. Over historical times, most lava flow activity has

Figure 1. Sketch map of Cameroon showing Mount Cameroon and other volcanoes (both oceanic and continental sectors of the line) along the Cameroon Volcanic Line (CVL); Inset map of Mount Cameroon showing historical lava flows for the last 110 years. The reported ages for these volcanoes were obtained from Marzoli et al. (2000).

Figure 2. Map of the geology of Mount Cameroon and environs compiled from past dissertations. Pyroclastic cones and the historical lava flows indicated by their years with the numbers in bracket (1, 2, 3) representing the different eruption sites were mapped from Digital Elevation Models (DEMs; 30 and 90 m) and other satellite imagery (multispectral Landsat TM, ETM+ and ASTER scenes). The rectangle shows where the 1982 flow was buried by the lower 1999 vents and lava.

been confined to the summit, SW and NE flanks of MC (Figure 2). This has effectively protected the SE and NW flanks from lava inundation.

Historical lava flows have flow lengths from ~850 m (upper 2000) to 11.5 km (1922, lower 1999; Figure 2; [11,41]). They cover an estimated surface area of 26.5 km² in the last 110 years at MC (Figure 2). Observations of effusive events at the MC summit are extremely minor, characterised by flows with limited extents that stay within the confines of the crater (e.g. 2000 site 1 flow; Figure 2). In addition, vent types differ between minor and major effusive events. Minor effusive events are fed by persistently active strombolian vents within the crater that host lava lakes which release the lava by breaching of the cones. In contrast, most vents on MC that feed major effusive events result from secondary vents/ephemeral boccas breaching or tapping from the central magma column somewhere below the level of the summit vents [11].

The next most abundant products observed at MC are pyroclastic cones (~340 cones have been mapped) aligned along a NE-SW trend (Figure 2), corresponding to the CVL alignment (Figure 1). They have the highest spatial density on the upper SW flank of MC (Figure 2). These cones result mostly from the observed moderately explosive activity and range from small spatter mounds to larger spatter-scoria cones. Other observed formations include lahars (volcanic mudflows; Figure 2) graded as a possible ancient hazard around MC [42]. Etindites (nephelinite lavas) are restricted to Mount Etinde (small volcano probably of Miocene age on the lower SW flank of MC; Figure 1; [35,54]). Lastly, are underlying sediments upon which all the above rock types lie (Figure 2; [35]).

2.2. Evolution of the 1982 Lava Flow Field

On 16 October 1982, new effusive activity started at MC and continued until 23 November 1982. This eruption occurred after 23 years of quiescence and emplaced a ~7.5 km-long basaltic 'a'ā lava flow field (Figure 3). A detailed account of the chronology/sequence of this eruption has been given by [47]. Thus, we summarise only elements of the chronology pertinent to the analysis of the lava flow field evolution presented here. Lava, ash, gases and tephra were emitted from a SW-NE trending fissure at an elevation of ~2700 m above sea level (a.s.l.) on the SW

Figure 3. (Colour online) 1982 lava flow field showing the location of the main scoria cone (represented by a square) and the 7 secondary vents (indicated by circles) that contributed to the flow width, thickness and lava volume emplaced for this eruption. (a) Photograph of active cone taken two days (7/11/1982) before the end of the eruption by G. J. Fitton; (b) emplaced scoria cone taken 28 years after (c) firsthand view of the lava channel captured by G. J. Fitton, 24 days after the onset of the eruption (9/11/1982); and (d) landslide debris as observed on the 9/11/1982 by G. J. Fitton.

flank of MC. The observed fissure extended over a distance of ~1 km down slope.

Moderately explosive activity at the upper end of the fissure led to the formation of a 25 m high scoria cone by November 7 (Figures 3(a) and (b)). A lava channel with an initial width of 3 m was observed in the first few days of eruption (Figure 3(c)). Lava flowed in this channel at a velocity of 5.2 m∙s⁻¹ and initial magma discharge rate of 10 - 30 m³∙s⁻¹ which dropped to ~0.1 m³∙s⁻¹ in the last 24 hours of eruption. Approximately 900 m down slope from the cone (8 November), measurements of an active lava lobe showed that it was 3.5 m wide, 1.25 m thick and advanced with a mean velocity of 0.2 m∙s⁻¹ on a slope of 18˚. Temperature readings at this point in the flow using a thermocouple produced values of 1045˚C ± 5˚C - 1070˚C ± 5˚C. This led to estimation of initial eruption temperature at 1160˚C and to average cooling rates of 0.1˚C - 0.3˚C∙s⁻¹. This lava lobe entered dense forest at 2400 m a.s.l., 1.5 km from the cone. Flow front measurements from the main flow lobe (Figure 3) at 1000 m a.s.l. gave values of ~200 m wide and 20 m high.

Between 30 October and 4 November, a landslide occurred just below the main cone (Figure 3(d)) that initiated a debris flow composed of a mixture of older volcanic material and the only just emplaced 1982 lava flow. This debris flow modified the upper part of the 1982 lava (Figure 3(d)) acting as a blockage zone and stopped 1.5 km from the vent. It produced an elongated depression with a length, width and depth of 500 m, 200 m and 150 m respectively. [47] estimated the total lava volume extruded at ~10⁷ m³. Other estimated parameters for this eruption include maximum pre-eruptive water content for the magma at 0.6 wt%, viscosity (10⁴ Pa∙s), yield strength (10⁴ Pa) as well as petrographic and geochemical data.

The emplacement of this lava flow field on moderately inclined slopes (12˚ - 19˚) resulted in distinctive flow field morphologies. The development of a compound lava flow field characterised by several lobes (Figure 3), present an excellent opportunity to examine the morphology and evolution of an intermediate/long-lived 'a'ā lava flow field emplaced on inclined surfaces. Unfortunately, during the 1982 eruption, virtually no information (e.g. GPS and ground control points for observed features), was collected from which the flow field morphology and evolution could be examined. Here, fieldwork was carried out to map out observed features and structures on this flow field and describe the surface characteristics of the flow. This data was then used to present details of the morphology and evolution of the flow field, and relate this evolution to quantitative parameters such as slope, effusion rate, number of active vents, flow length, lava tubes and number of active flow lobes. In doing this, we define the characteristics of the 1982 basaltic lava flow field inferred from its emplacement dynamics. From field observations and the narrative given by [47], the lava flow was extremely unstable and suffered from almost constant flow front collapses to feed other flows extending down-slope from stalled flow fronts.

3. Materials and Methods

Field observations and measurements (channel width, depth and levee angle) were made at the proximal, medial and dispersed flow portions of the flow field. Details of the surface morphology for the different lava types and structures observed in this flow field were obtained across transversal profiles made within the stable and transitional flow zones. 13 profiles (Figure 4) were measured using a 30 m tape, an abney level and a compass clinometer. They enabled to describe the physical characteristics of the 1982 lava flow field. Four representative profiles are illustrated here to show changes in channel morphology down-flow (Figures 4(a)-(d)).

Morphology of the flow field at a distance (length) above 5.5 km (that comprised the dispersed and flow front zones) could not be observed because this portion of the flow is buried by the lower 1999 vents and lava (Figures 1 and 2). Several bifurcations, islands and flow lobes are observed. Our focus is on the WSW branch of the flow field (Figure 3).

Lava thickness was derived from trigonometry. These values were substituted in equations from [55,56] and into Jeffrey’s equation to estimate yield strength (levees)

Figure 4. (Colour online) 1982 lava flow map showing profile locations indicated as points that represent the start and end GPS positions of some of the profiles. The illustrated profiles (P1-P4) were made (a) 1 m from the vent immediately after the observed zone of blockage in Figure 5(a); (b) 200 m from the vent showing a characteristic trend of clinker 'a'ā morphology that covers the entire section; (c) 1 km from the vent, dominated by blocky 'a'ā lava morphology and (d) 5 km from the vent showing a dominant clinker 'a'ā morphology (larger blocks) with accreted channel margins. P1 represent the entire flow field width, while P2-P4 represents just the stable channel width. Question marks and arrows are indicative of side channels and the continuation of the flow field.

and mean channel velocity (Equations (1) and (2)).

(1)

where t is yield strength, r is the dense rock equivalent density, g is acceleration due to gravity which is equal to 9.8 m∙s⁻², t is levee thickness and a is the gradient of slope (slope here represents pre-eruptive down flow slope).

(2)

where V is mean velocity, a is the pre-eruptive down flow slope (obtained from the field and the 30 m Digital Elevation Model (DEM) for MC) and B is the shape constant considered to be equal to 3 (value is representative of wide channels).

From the estimated velocity (V), channel depth (d) and width (W), effusion rate for lava flowing in these channels was derived (see [57]). A uniform density of 2700 kg∙m⁻³ was assumed in estimating the above rheological parameters. This is an average density for this lava flow estimated based on previously published whole-rock geochemical compositional data from [25,39,47]) following [58] approach. Lava volume for this eruption was reestimated by the planimetric approach proposed by [59]. This involves the measurement of the area covered by the lava flow field (obtained from GIS estimate) multiplied by an estimated mean lava thickness for the flow field.

For the purpose of reconstructing the stable channel geometry, profiles were made at intervals of 100/200 m (for the first one kilometre) and at 1 km interval from then onwards (Table 1) starting with a profile within a few metres from the vent (Figure 4(a)). Other observations from this flow were obtained between profiles.

4. Results

4.1. Lava Flow Field Morphology

The 1982 compound 'a'ā lava flow field (Figures 4(a)-(d)) has 6 main and uncountable subsidiary lava flow lobes which branch out in several directions, with some of them meeting up again towards the flow front (Figure 3). The principal flow morphology observed from near-vent to 5.5 km distance for this lava flow field is clinker 'a'ā lava followed closely by blocky 'a'ā lava. Upper flow proximal lava in the first 17 m from the vent is slabcrusted lava within the stable channel according to the terminology given by [60]. Most morphological features indicative of lava flow emplacement mechanism and late stage deformation were recorded.

4.1.1. Large-Scale Flow Field Development and Structures

1) Scoria Cone The development of the main cone thought to have fed the greatest volume of the lava flow field was partially documented by [47]. This cone is visible both on Landsat TM, ETM+ scenes, and in the field (Figure 5(a)) emplaced on slopes of 16˚ - 18˚ and breached in the downflow direction. Measurements (GIS) made from the already emplaced 25 m high cone gave ~350 m base diameter and ~85 m crater diameter, occupying an area of ~1.2 × 10⁵ m² and volume of ~2.5 × 10⁶ m³ (~20 m thick) respectively. Lava-fountaining up to 300 - 400 m high from this cone ([47]; Figure 3(a)) in its explosive phase, led to emplacement of a tephra layer observed for a distance of up to ~200 m away in the NE direction (along the route used to get to the vent). This layer was sparsely dotted by large bombs (~1 m) made up mostly of lapilli and sand-sized particles close to vent that became finer away from it (Figure 5(a)). The scoria and lapilli range in size from 2 to 65 mm in diameter.

The blockage (Figure 5(a)) observed immediately after the cone breached zone was produced ~14 days after eruption onset [47]. It resulted from a landslide below the cone that produced a debris flow of old volcanic material (ash and scoria) mixed with disrupted remains of the 1982 lava [47]. This blockage came after the establishment of the stable channel (Figure 5(a)) observed immediately after this zone [47]. The observed blockage at the time of eruption prevented the debris flow from taking the direction of the channel. The debris flow deposit was still preserved in the field at the time of the study in the upper sector of the flow field (Figure 5(a)). This debris flow ended with a sharply defined lobate front as described by [47]. It obstructed the flow of lava in the established channel (Figure 5(a)) causing diversion of the flow that led to the production of a secondary channel to the West of the vent which merged a few metres downflow with the primary vent.

2) Secondary Vents Field investigations revealed a series of 7 late-forming secondary and/or ephemeral (short-lived) vents (Figures 5(b) and (c)). These vents are observed ~2.5 km away from the main vent (Figures 5(b) and (c)) at ~2180 m a.s.l, emplaced on slopes of 15˚ - 19˚, all aligned in the E-W direction within the flow field. These small vent constructs are 3 to 15 m high (Figure 5(b)), are most often breached in the SE direction (Figure 5(c)), and are spaced at intervals of 2 to 10 m of each other. Threequarters of them produced lava with clinker 'a'ā morphology that flowed in distinct channels within the flow field which disappeared ~20 m down flow. Most of their walls have already been colonised by mosses. Close inspection revealed that some of these walls are plastered with spatter and clinker (Figures 5(c) and (d)).

3) Levees Few initial levees were preserved because of variations in flow level and blockages within channels that fa-

Table 1. Geometrical parameters of the stable channel zone of the 1982 lava flow field.

voured over flow events and formation of overflow levees. These events modify levees (Figure 6(a)) and cause flow field widening. Initial levees are characterised by outer slopes of 26˚ to 66˚. These levees are 0.3 to 2 m thick (t) with corresponding levee widths of 0.4 to 3 m (Figures 4(a)-(d), Profiles 1 - 4). Yield strength of initial lava estimated using outer initial levee thickness and angle of rest [55] fluctuates between 10³ to 10⁴ Pa across the different profiles. The main lava types that characterised these levees are loose clinker and blocky 'a'ā lava

Figure 5. (Colour online) Photographs showing (a) scoria cone with breached zone, tephra, debris flow and the start of channel and (b) and (c). Alignment of secondary vents within channel showing a breached zone.

(Figures 4(a)-(d), Profiles 1 - 4). The inner walls of these levees show horizontal layering and shearing features (Figure 6(b)).

Overflow levees could be distinguished from initial levees because they showed layering of different lava levels believed to represent different flow episodes. Overflow material is basically clinker and blocky 'a'ā lava (Figure 6(a)) and to a limited extent rubble 'a'ā lava (Figure 4(c), profile 3). Typical outer slopes for overflow levees are 30˚ to 50˚ with levee thickness (t) values of 0.4 to 2.5 m and levee widths of 0.3 to 2.6 m.

Accretionary levees (Figures 6(b)-(e)) are the dominant levee type. Accreted/agglutinated lava surfaces were observed both on the outer and inner walls of channels. They are characterised by brecciated surfaces (Figure 6(c)) showing shearing features, tension fractures and cooling cracks (Figure 6(d)). They are 0.5 to 8 m thick. Virtually all of the internal accretionary levees observed are characterised by sub-vertical to vertical inner levee walls (Figures 4(a)-(d), Profiles 1 - 4) with a massive interior texture (Figure 6(e)). Yield strength estimated for this levee type is stable at 10⁴ Pa. Since the flow front could not be observed at the time of this study, no flow front measurements were made. However, [47] gave a width of ~200 m with a thickness of 20 m for the main flow lobe.

4) Lava Channels In this flow field, as for other 'a'ā dominated lava flow fields at MC (e.g. 1959, 1999) channels were the dominant preserved transport pathway. Lava in proximal channels (first 100 m) flowed with initial estimated velocity and effusion rate of 4 m∙s⁻¹ and 3.3 m³∙s⁻¹ respectively. Channels are 0.85 - 5 m wide, dominated by slab-crusted lava with clinker 'a'ā margins (Figure 5(a); Table 1) observed up to a distance of 205 m (Figure 7(a)) away from vent. However, these crusted surfaces are not continuous for the whole length of the 205 m. They broke up into slabs after covering distances of 10 to 17 m (Figure 7(b)). These distances correspond to an increase in slope from 16˚ (at-vent) to 18˚ (300 m). These broken slabs probably acted as temporary blockage zones within these channels (Figure 7(c)) while this flow was active. Broken slab material usually piles up forming heaps up to 13 m long at some locations (Figure 7(c)).

Within a few metres (~20 m) from vent, secondary channels were observed at the sides of the primary or main channel (Figure 7(c)). These channels are 2 to 3 m wide, 100 m from vent Figure 7(c)), and 30 m wide ~400 m from vent (Table 1). Effusion rate increased to ~29 m³∙s⁻¹ (100 m distance) and dropped to 13 m³∙s⁻¹ 200 m away from vent. This drop in effusion rate led to the production of channels (~200 m from vent) characterised by a reduced width (1.5 m) and an increased depth (5.4 m) with a slab-crusted surface that continued for a distance of 10 m before breaking up into slabs (Figures 7(a) and (d)). Velocity estimates at these distances had dropped to 2 and 1.4 m∙s⁻¹. Both effusion rate and velocity (12 m∙s⁻¹) increased from these distances up to 1 km distance. This favoured the formation of wider channels (6 - 21 m) and clinker and blocky 'a'ā lava emplaced on slopes of 14˚ - 17˚. The widest channel (113.5 m) in this flow field was observed ~2 km from vent

Figure 6. (Colour online) Levees observed at the 1982 lava flow field showing (a) initial/overflow levee over-underlain by clinker 'a'ā; (b) accretionary levee wall showing horizontal layering (broken lines) and fracturing; (c) outer accreted levee wall with highly brecciated surface; (d) levee wall with cooling cracks and tension fractures and (e). massive levee interior with arrow showing lava flow direction within channel.

hosting sealed lava tube surfaces, characterised by reduced velocity (4.5 m∙s⁻¹) and high flow rates (68 m³∙s⁻¹). Channel material is clinker 'a'ā emplaced on a slope of 15˚.

Lava channels were still considerable wide (10 - 22 m) from this distance but dropped to 6 m, 5.5 km from vent. At this distance, velocity and flow rate had dropped to 0.3 m∙s⁻¹ and 13 m³∙s⁻¹ respectively. The lava type here is clinker and blocky 'a'ā lava emplaced on slopes of 12˚ - 14˚. From the estimates made, an average velocity and flow rate of 4.5 m∙s⁻¹ and 26 m³∙s⁻¹ were produced for this flow. Taking an average lava thickness of 5 m (±2 m; field estimates), an area of 2.6 × 10⁶ m² (GIS estimate), this eruption produced lava with a volume of ~1.3 × 10⁷ m³.

5) Lava Tubes Lava tubes were first observed in channels at a distance 600 m from vent (Figure 7(b)). They could be distinguished based on their elongated forms, continuous surface morphology, raised margins and extensive lengths relative to that of the surrounding rocks (Figures 7(b) and (c)). These tubes are characterised by a spinose surface with crusts just a few cm thick (Figure 7(b)). At 600 m distance, an observed lava tube is 10 m long and ~2 m wide (Figure 7(b)). These lava tubes were abundant within the primary channel (~2 km from vent) where the channel is widest. From their surface appearance, 5 could be distinguished with extensive lengths (~10 - 20 m) that became buried down-slope by channel material. This is 500 m (up-slope) from the point where

Figure 7. (Colour online) Channel morphology across the 1982 lava flow field showing (a) intermediate stable channel with reduced width and increased depth with slab-crusted floor; (b) lava-tube/slab-crusted surface with a broken surface (encircled); (c) blockage within channel characterised by tilted pahoehoe-slabs hosting side channels on both sides of the stable channel zone bounded by levees and (d) continuation of intermediate stable channel zone observed further down-stream after slabs had broken up in (b) above. Arrows within channels are indicative of flow direction.

secondary vents were observed.

4.2. Down-Flow Channel Geometry and Small Scale-Structures

Here, on the basis of channel geometry fluctuation, we split the stable channel zone of the flow into three sections: 1) an upper section characterised by wider than deep channels; 2) an intermediate section governed by deeper than wide channels and 3) a lower section where channels resume their normal trend as in the upper section being wider than deep. This zone is followed by the transitional channel zone observed ~4 km down-flow.

4.2.1. Stable Channel Zone

Table 1 and Figures 8(a) and (b) give details on the stable channel geometry at a distance starting from 1 m from the vent, progressively down-stream. The widths of breakaway/side channels are also given at the corresponding

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	S. R. Passey and B. R. Bell, “Morphologies and Emplacement Mechanisms of the Lava Flows of the Faroe Islands Basalt Group, Faroe Islands, NE Atlantic Ocean,” Bulletin of Volcanology, Vol. 70, No. 2, 2007, pp. 139-156. doi:10.1007/s00445-007-0125-6
[2]	M. S. Njome, C. E. Suh, R. J. Sparks, S. N. Ayonghe and J. G. Fitton, “The Mount Cameroon 1959 Compound Lava Flow Field: Morphology, Petrology and Geochemistry,” Swiss Journal of Geosciences, Vol. 101, No. 1, 2008, pp. 85-98. doi:10.1007/s00015-007-1245-x
[3]	H. C. Sheth, J. S. Ray, R. Bhutani, A. Kumar and R. S. Smitha, “Volcanology and Eruptive Styles of Barren Island: An Active Mafic Stratovolcano in the Andaman Sea, NE Indian Ocean,” Bulletin of Volcanology, Vol. 71, No. 9, 2009, pp. 1021-1039. doi:10.1007/s00445-009-0280-z
[4]	H. C. Sheth, J. S. Ray, A. Kumar, R. Bhutani and N. Awasthi, “Toothpaste Lava from the Barren Island Volcano (Andaman Sea),” Journal of Volcanology and Geothermal Research, Vol. 202, No. 1-2, 2011, pp. 73-82. doi:10.1016/j.jvolgeores.2011.01.006
[5]	L. J. Applegarth, H. Pinkerton, M. R. James and S. Calvari, “Lava Flow Superposition: The Reactivation of Flow Units in Compound 'a'ā Flows,” Journal of Volcanology and Geothermal Research, Vol. 194, No. 4, 2010, pp. 100-106. doi:10.1016/j.jvolgeores.2010.05.001
[6]	L. J. Applegarth, H. Pinkerton, M. R. James and S. Calvari, “Morphological Complexities and Hazards during the Emplacement of Channel-Fed 'a'ā Lava Flow Fields: A Study of the 2001 Lower Flow Field on Etna,” Bulletin of Volcanology, Vol. 72, No. 6, 2010, pp. 641-656. doi:10.1007/s00445-010-0351-1
[7]	H. El Hachimi, N. Youbi, J. Madeira, L. Martins, A. Marzoli, H. Bertrand, G. Bellieni, S. Callegaro, J. Mata, J. M. Munhá, F. Medina, A. Mahmoudi, K. M. Bensalah, B.M. Abbou, “Morphologies and Emplacement Mechanisms of the Lava Flows of the Central Atlantic Magmatic Province (CAMP) of Morocco,” Proceedings of the II Central & North Atlantic Conjugate Margins Conference, Morocco, 29 September-1 October 2010, pp. 96-100.
[8]	C. E. Suh, S. A. Stansfield, R. S. J Sparks, M. S. Njome, M. N. Wantim and G. G. J. Ernst, “Morphology and Structure of the 1999 Lava Flows at Mount Cameroon Volcano (West Africa) and their Bearing on the Emplacement Dynamics of Volume-Limited Flows,” Geological Magazine, 2010, pp. 1-13.
[9]	D. Takagi and H. E. Huppert, “Initial Advance of Long Lava Flows in Open Channels,” Journal of Volcanology and Geothermal Research, Vol. 195, No. 2-4, 2010, pp. 121-126. doi:10.1016/j.jvolgeores.2010.06.011
[10]	R. J. Brown, S. Blake, N. R. Bondre, V. M. Phadnis and S. Self, “'a'ā Lava Flows in the Deccan Volcanic Province, India, and their Significance for the Nature of Continental Flood Basalt Eruptions,” Bulletin of Volcanology, Vol. 73, No. 6, 2011, pp. 737-752. doi:10.1007/s00445-011-0450-7
[11]	M. N. Wantim, C. E. Suh, G. G. J. Ernst, M. Kervyn and P. Jacobs, “Characteristics of the 2000 Fissure Eruption and Lava Flow Fields at Mount Cameroon Volcano, West Africa: A Combined Field Mapping and Remote Sensing Approach,” Geological Journal, Vol. 46, No. 4, 2011, pp. 344-363. doi:10.1002/gj.1277
[12]	S. W. Anderson, S. E. Smrekar and E. R. Stofan, “Tumulus Development on Lava Flows: Insights from Observations of Active Tumuli and Analysis of Formation Models,” Bulletin of Volcanology, Vol. 74, No. 4, 2012, pp. 931-946. doi:10.1007/s00445-012-0576-2
[13]	J. E. Bailey, A. J. L. Harris, J. Dehn, S. Calvari and S. K. Rowland, “The Changing Morphology of an Open Channel on Mount Etna,” Bulletin of Volcanology, Vol. 68, No. 6, 2006, pp. 497-515. doi:10.1007/s00445-005-0025-6
[14]	A. J. L. Harris, L. P. Flynn, O. Matias, W. I. Rose and J. Cornejo, “The Evolution of An Active Silicic Lava Flow Field: An ETM+ Perspective,” Journal of Volcanology and Geothermal Research, Vol. 135, No. 1-2, 2004, pp. 147-168. doi:10.1016/j.jvolgeores.2003.12.011
[15]	A. J. L. Harris, M. Favalli, M. Francesco and C. W. Hamilton, “Construction Dynamics of A Lava Channel,” Bulletin of Volcanology, Vol. 71, No. 4, 2009, pp. 459-474. doi:10.1007/s00445-008-0238-6
[16]	M. Favalli, A. J. L. Harris, A. Fornaciai, M. T. Pareschi and F. Mazzarini, “The Distal Segment of Etna’s 2001 Basaltic Lava Flow,” Bulletin of Volcanology, Vol. 72, No. 4, 2010, pp. 119-127. doi:10.1007/s00445-009-0300-z
[17]	M. R. James, L. J. Applegarth and H. Pinkerton, “Lava Channel Roofing, Overflows, Breaches and Switching: Insights from the 2008-2009 Eruption of Mt. Etna,” Bulletin of Volcanology, Vol. 74, No. 1, 2012, pp. 107-117. doi:10.1007/s00445-011-0513-9
[18]	A. J. L. Harris and S. K. Rowland, “FLOWGO: A Kinematic Thermo-Rheological Model for Lava Flowing in A Channel,” Bulletin of Volcanology, Vol. 63, No. 1, 2001, pp. 20-44. doi:10.1007/s004450000120
[19]	G. M. Crisci, S. Di Gregorio, R. Rongo, M. Scarpelli, W. Spataro and S. Calvari, “Revisiting the 1669 Etnean Eruptive Crisis Using Cellular Automata Model and Implications for Volcanic Hazard in the Catania Area,” Journal of Volcanology and Geothermal Research, Vol. 123, No. 1-2, 2003, pp. 211-230. doi:10.1016/S0377-0273(03)00037-4
[20]	G. M. Crisci, R. Rongo, S. Di Gregorio and W. Spataro, “The Simulation Model Sciara: The 1991 and 2001 Lava Flows at Mount Etna,” Journal of Volcanology and Geothermal Research, Vol. 132, No. 2-3, 2004, pp. 253-267. doi:10.1016/S0377-0273(03)00349-4
[21]	G. M. Crisci, G. Iovine, S. Di Gregorio and V. Lupiano, “Lava-flow Hazard on the SE Flank of Mt. Etna (Southern Italy),” Journal of Volcanology and Geothermal Research, Vol. 177, No. 4, 2008, pp. 778-796. doi:10.1016/j.jvolgeores.2008.01.041
[22]	A. Ciraudo, C. Del Negro, A. Herault and A. Vicari, “Advances in Modelling Methods for Lava Flows Simulation,” Communications to SIMAI Congress, Vol. 2, No. 4, 2007, pp. 1827-9015.
[23]	A. Vicari, H. Alexis, C. Del Negro, M. Coltelli, M. Marsella and C. Proietti, “Modelling of the 2001 Lava Flow at Etna Volcano by a Cellular Automata Approach,” Environmental Modelling and Software, Vol. 22, No. 4, 2007, pp. 1465-1471. doi:10.1016/j.envsoft.2006.10.005
[24]	Ruxton, “Report on Volcanic Eruptions on the Cameroons Mountain,” 1922.
[25]	C. E. Suh, R. S. J. Sparks, J. G. Fitton, S. N. Ayonghe, C. Annen, R. Nana and A. Luckman, “The 1999 and 2000 Eruptions of Mount Cameroon: Eruption Behaviour and Petrochemistry of Lava,” Bulletin of Volcanology, Vol. 65, No. 4, 2003, pp. 267-281. doi:10.1007/s00445-002-0257-7
[26]	S. Calvari and H. Pinkerton, “Formation of Lava Tubes and Extensive Flow Field during the 1991-93 Eruption of Mount Etna,” Journal of Geophysical Research, Vol. 103, No. B11, 1998, pp. 27291-27302. doi:10.1029/97JB03388
[27]	C. R. J. Kilburn and R. M. C. Lopes, “The Growth of 'a'ā Lava Flow Fields on Mount Etna, Sicily,” Journal of Geophysical Research, Vol. 93, No. B12, 1988, pp. 759-772. doi:10.1029/JB093iB12p14759
[28]	A. Felpeto, V. Arana, R. Ortiz, M. Astiz and A. Garcia, “Assessment and Modelling of Lava Flow Hazard on Lanzarote (Canary Islands),” Natural Hazards, Vol. 23, No. 2-3, 2001, pp. 247-257. doi:10.1023/A:1011112330766
[29]	R. A. Duraiswami, N. R. Bondre, G. Dole, V. M. Phadnis and V. S. Kale, “Tumuli and Associated Features from the Western Deccan Volcanic Province, India,” Bulletin of Volcanology, Vol. 63, No. 7, 2001, pp. 435-442. doi:10.1007/s004450100160
[30]	C. R. J. Kilburn and R. M. C. Lopes, “General Patterns of Flow Field Growth: 'a'ā and Blocky Lavas,” Journal of Geophysical Research, Vol. 96, No. B12, 1991, pp. 721-732. doi:10.1029/91JB01924
[31]	J. E. Guest, A. M. Duncan, E. R. Stofan and S. W. Anderson, “Effect of Slope on Development of Pahoehoe Flow Fields: Evidence from Mount Etna,” Journal of Volcanology and Geothermal Research, Vol. 219, No. 11, 2012, pp. 52-62. doi:10.1016/j.jvolgeores.2012.01.006
[32]	C. R. J. Kilburn, “Lava Flows and Flow Fields,” In: H. Sigurdsson, Ed., Encyclopaedia of Volcanoes, Academic Press, San Diego, 2000, pp. 291-305.
[33]	G. P. L. Walker, “Structure and Origin by Injection of Lava under Surface Crust, of Tumuli, ‘Lava Rises’, ‘Lava-Rise Pits’, and ‘Lava-Inflation Clefts’ in Hawaii,” Bulletin of Volcanology, Vol. 53, No. 7, 1991, pp. 546-558. doi:10.1007/BF00298155
[34]	J. G. Fitton, C. R. J Kilburn, M. F. Thirlwall and D. J. Hughes, “1982 Eruption of Mount Cameroon, West Africa,” Nature, Vol. 306, No. 5941, 1983, pp. 327-332. doi:10.1038/306327a0
[35]	B. Deruelle, J. N’ni and R. Kambou, “Mount Cameroon: An Active Volcano of the Cameroon Line,” Journal of African Earth Sciences, Vol. 6, No. 2, 1987, pp. 197-214.
[36]	B. Ateba and N. Ntepe, “Post-Eruptive Activity of Mount Cameroon (Cameroon) West Africa: A Statistical Analysis,” Journal of Volcanology and Geothermal Research, Vol. 79, No. 1-2, 1997, pp. 25-45. doi:10.1016/S0377-0273(97)00022-X
[37]	R. U. Ubangoh, B. Ateba, S. N. Ayonghe and G. E. Ekodeck, “Earthquake Swarms of Mount Cameroon, West Africa,” Journal of African Earth Sciences, Vol. 24, No. 4, 1997, pp. 413-424. doi:10.1016/S0899-5362(97)00072-9
[38]	C. E. Suh, S. N. Ayonghe and E. S. Njumbe, “Neotectonic Earth Movement Related to the 1999 Eruption of Cameroon Mountain, West Africa,” Episodes, Vol. 24, No. 1, 2001, pp. 9-12.
[39]	C. E. Suh, J. F. Luhr and M. S. Njome, “Olivine-Hosted Glass Inclusions from Scoriae Erupted in 1954-2000 at Mount Cameroon Volcano, West Africa,” Journal of Volcanology and Geothermal Research, Vol. 169, No. 1-2, 2008, pp. 1-33. doi:10.1016/j.jvolgeores.2007.07.004
[40]	B. Ateba, C. Dorbath, L. Dorbath, N. Ntepe, M. Frogneux, F. T. Aka, J. V. Hell, J. C. Delmond and D. Manguelle, “Eruptive and Earthquake Activities Related to the 2000 Eruption of Mount Cameroon Volcano (West Africa),” Journal of Volcanology and Geothermal Research, Vol. 179, No. 3-4, 2009, pp. 206-216. doi:10.1016/j.jvolgeores.2008.11.021
[41]	K. Bonne, M. Kervyn, L. Cascones, S. Njome, E. V. Ranst, E. Suh, S. Ayonghe, P. Jacobs and G. Ernst, “A New Approach to Assess Long-Term Lava Flow Hazard and Risk Using GIS and Low-Cost Remote Sensing: the Case of Mount Cameroon, West Africa,” International Journal of Remote Sensing, Vol. 29, No. 22, 2008, pp. 6537-6562. doi:10.1080/01431160802167873
[42]	P. Thierry, L. Stieltjes, E. Kouokam, P. Nguéya and P. M. Salley, “Multi-Hazard Risk Mapping and Assessment on An Active Volcano: The GRINP Project at Mount Cameroon,” Natural Hazards, Vol. 45, No. 3, 2008, pp. 429-456. doi:10.1007/s11069-007-9177-3
[43]	M. Favalli, S. Tarquini, P. Papale, A. Fornaciai and E. Boschi, “Lava Flow Hazard and Risk at Mt. Cameroon Volcano,” Bulletin of Volcanology, Vol. 74, No. 2, 2012, pp. 423-439. doi:10.1007/s00445-011-0540-6
[44]	P. Lipman and N. G. Banks, “'a'ā Flow Dynamics, Mauna Loa.” In: R. W. Decker, T. L. Wright and P. H. Stauffer, Eds., Volcanism in Hawaii, US Geological Survey Professional Paper, Vol. 1350, No. 2, 1987, pp. 1527-1567.
[45]	L. Lodato, L. Spampinato, A. J. L. Harris, S. Calvari, J. Dehn and M. Patrick, “The Morphology and Evolution of the Stromboli 2002-2003 Lava Flow Field: An Example of A Basaltic Flow Field Emplaced on A Steep Slope,” Bulletin of Volcanology, Vol. 69, No. 6, 2007, pp. 661-697. doi:10.1007/s00445-006-0101-6
[46]	K. Hon, J. Kauahikaua, R. Denlinger and K. Mackay, “Emplacement and Inflation of Pahoehoe Sheet Flows: Observations and Measurements of Active Lava Flows on Kilauea Volcano, Hawaii,” Geological Society of America Bulletin, Vol. 106, No. 3, 1994, pp. 351-370. doi:10.1130/0016-7606(1994)106<0351:EAIOPS>2.3.CO;2
[47]	J. G. Fitton and H. M. Dunlop, “The Cameroon Line, West Africa, and Its Bearing on the Origin of Oceanic and Continental Alkali Basalt,” Earth and Planetary Science Letters, Vol. 72, No. 1, 1985, pp. 23-28. doi:10.1016/0012-821X(85)90114-1
[48]	C. Moreau, J-M. Regnoult, B. Deruelle and B. Robineau, “A New Tectonic Model for the Cameroon Line, Central Africa,” Tectonophysics, Vol. 139, 1987, pp. 317-334.
[49]	A. N. Halliday, A. P. Dickin, A. E. Fallick and J. G. Fitton, “Mantle Dynamics: A Nd, Sr, Pb and O isotopic study of the Cameroon Line Volcanic Chain,” Journal of Petrology, Vol. 29, No. 1, 1988, pp. 181-211. doi:10.1093/petrology/29.1.181
[50]	C. T. Tabod, J. D. Fairhead, G. W. Stuart, B. Ateba and N. Ntepe, “Seismicity of the Cameroon Volcanic Line, 1982-1990,” Tectonophysics, Vol. 212, No. 3-4, 1992, pp. 303-320. doi:10.1016/0040-1951(92)90297-J
[51]	F. T. Aka, K. Nagao, M. Kusakabe, H. Sumino, G. Tanyileke, B. Ateba and J. Hell, “Symmetrical Helium Isotope Distribution on the Cameroon Volcanic Line, West Africa,” Chemical Geology, Vol. 203, No. 3-4, 2004, pp. 205-223. doi:10.1016/j.chemgeo.2003.10.003
[52]	B. Deruelle, I. Ngounouno and D. Demaiffe, “The ‘Cameroon Hot Line’ (CHL): A Unique Example of Active Alkaline Intraplate Structure in both Oceanic and Continental Lithospheres,” Comptes Rendus Geoscience, Vol. 339, No. 9, 2007, pp. 589-600. doi:10.1016/j.crte.2007.07.007
[53]	L. Mathieu, M. Kervyn and G. G. J. Ernst, “Field Evidence for Flank Instability, Basal Spreading and Volcano-Tectonic Interactions at Mt. Cameroon, West Africa,” Bulletin of Volcanology, Vol. 73, No. 7, 2011, pp. 851-867. doi:10.1007/s00445-011-0458-z
[54]	C. Nkoumbou, B. Deruelle and D. Velde, “Petrology of Mt. Etinde Nephelinite Series,” Journal of Petrology, Vol. 36, No. 2, 1995, pp. 373-393. doi:10.1093/petrology/36.2.373
[55]	G. Hulme, “The Interpretation of Lava Flow Morphology,” Geophysical Journal Royal Astronomical Society, Vol. 39, No. 2, 1974, pp. 361-383. doi:10.1111/j.1365-246X.1974.tb05460.x
[56]	H. J. Moore, “Preliminary Estimates of the Rheological Properties of 1984 Mauna Loa Lava.” In: R. W. Decker, T. L. Wright, P. H. Stauffer, Eds., Volcanism in Hawaii, United States Government Printing Office, Washington DC, 1987, pp. 1569-1588.
[57]	A. J. L. Harris, A. L. Butterworth, R. W. Carlton, I. Downey, P. Miller, P. Navarro and D. A. Rothery, “Low-Cost Volcano Surveillance from Space: Case Studies from Etna, Rafla, Cerro Negro, Fogo, Lascar and Erebus,” Bulletin of Volcanology, Vol. 59, No. 1, 1997, pp. 49-64. doi:10.1007/s004450050174
[58]	Y. Bottinga and D. F. Weill, “Densities of Liquid Silicate Systems Calculated from Partial Molar Volumes of Oxide Components,” American Journal of Science, Vol. 269, No. 2, 1970, pp. 169-182. doi:10.2475/ajs.269.2.169
[59]	N. F. Stevens, G. Wadge and J. B. Murray, “Lava Flow Volume and Morphology from Digitised Contour Maps: A Case Study at Mount Etna, Sicily,” Geomorphology, Vol. 28, No. 3-4, 1999, pp. 251-261. doi:10.1016/S0169-555X(98)00115-9
[60]	M. Polacci and P. Papale, “The Evolution of Lava Flows from Ephemeral Vents at Mount Etna: Insights from Vesicle Distribution and Morphological Studies,” Journal of Volcanology and Geothermal Research, Vol. 76, No. 1-2, 1997, pp. 1-17. doi:10.1016/S0377-0273(96)00070-4
[61]	J. A. Naranjo, R. S. J. Sparks, M. V. Stasiuk, H. Moreno and G. J. Ablay, “Morphological, Structural and Textural Variations in the 1988-1990 Andesite Lava of Lonquimay Volcano, Chile,” Geological Magazine, Vol. 129, No. 6, 1992, pp. 657-678. doi:10.1017/S0016756800008426
[62]	M. J. Kennish and R. A. Lutz, “Morphology and Distribution of Lava Flows on Mid-Ocean Ridges: A Review,” Earth Science Reviews, Vol. 43, No. 3-4, 1998, pp. 63-90. doi:10.1016/S0012-8252(98)00006-3
[63]	S. Calvari and H. Pinkerton, “Lava Tube Morphology on Etna and Evidence for Lava Flow Emplacement Mechanisms,” Journal of Volcanology and Geothermal Research, Vol. 90, No. 3-4, 1999, pp. 263-280. doi:10.1016/S0377-0273(99)00024-4
[64]	M. Dragoni and A. Piombo, “Thermoelastic Deformation Associated with A Lava Tube,” Bulletin of Volcanology, Vol. 71, No. 4, 2009, pp. 409-418. doi:10.1007/s00445-008-0248-4
[65]	R. S. J. Sparks, H. Pinkerton and G. Hulme, “Classification and Formation of Lava Levees on Mount Etna, Sicily,” Geology, Vol. 4, No. 5, 1976, pp. 269-271. doi:10.1130/0091-7613(1976)4<269:CAFOLL>2.0.CO;2
[66]	F. Quareni, A. Tallarico and M. Dragoni, “Modelling of the Steady-State Temperature Field in Lava Flow Levees,” Journal of Volcanology and Geothermal Research, Vol. 132, No. 2-3, 2004, pp. 241-251. doi:10.1016/S0377-0273(03)00348-2
[67]	G. A. Macdonald, “Pahoehoe, 'a'ā and Block Lava,” American Journal of Science, Vol. 251, No. 3, 1953, pp. 169-191. doi:10.2475/ajs.251.3.169
[68]	C. R. J. Kilburn, “Fracturing as a Quantitative Indicator of Lava Flow Dynamics,” Journal of Volcanology and Geothermal Research, Vol. 132, No. 2-3, 2004, pp. 209-224. doi:10.1016/S0377-0273(03)00346-9
[69]	J. E. Guest and E. R. Stofan, “The Significance of Slab-Crusted Lava Flows for Understanding Controls on Flow Emplacement at Mount Etna, Sicily,” Journal of Volcanology and Geothermal Research, Vol. 142, No. 3-4, 2006, pp. 193-205. doi:10.1016/j.jvolgeores.2004.09.003
[70]	P. Francis and C. Oppenheimer, “Volcanoes,” 2nd Edition, Oxford University Press, Oxford, 2004.
[71]	M. J. Rossi and A. Gudmundsson, “The Morphology and Formation of Flow-Lobe Tumuli on Icelandic Shield Volcanoes,” Journal of Volcanology and Geothermal Research, Vol. 72, No. 3-4, 1996, pp. 291-308. doi:10.1016/0377-0273(96)00014-5
[72]	A. M. Duncan, J. E. Guest, E. R. Stofan, S. W. Anderson, H. Pinkerton and S. Calvari, “Development of Tumuli in the Medial Portion of the 1983 'a'ā Flow field, Mount Etna, Sicily,” Journal of Volcanology and Geothermal Research, Vol. 132, No. 2-3, 2004, pp. 173-187. doi:10.1016/S0377-0273(03)00344-5
[73]	R. Greeley, “The Role of Lava Tubes in Hawaiian Volcanoes.” In: R. W. Decker, T. L. Wright, P. H Stauffer, Eds., Volcanism in Hawaii, United States Government Printing Office, Washington, 1993, pp. 1589-1602.