The changing morphology of an open lava channel on Mt. Etna

An open channel lava flow on Mt. Etna (Sicily) was observed during May 30–31, 2001. Data collected using a forward looking infrared (FLIR) thermal camera and a Minolta-Land Cyclops 300 thermal infrared thermometer showed that the bulk volume flux of lava flowing in the channel varied greatly over time. Cyclic changes in the channel's volumetric flow rate occurred over several hours, with cycle durations of 113–190 min, and discharges peaking at 0.7 m³ s⁻¹ and waning to 0.1 m³ s⁻¹. Each cycle was characterized by a relatively short, high-volume flux phase during which a pulse of lava, with a well-defined flow front, would propagate down-channel, followed by a period of waning flow during which volume flux lowered. Pulses involved lava moving at relatively high velocities (up to 0.29 m s⁻¹) and were related to some change in the flow conditions occurring up-channel, possibly at the vent. They implied either a change in the dense rock effusion rate at the source vent and/or cyclic-variation in the vesicle content of the lava changing its bulk volume flux. Pulses would generally overspill the channel to emplace pāhoehoe overflows. During periods of waning flow, velocities fell to 0.05 m s^–1. Blockages forming during such phases caused lava to back up. Occasionally backup resulted in overflows of slow moving ‘a‘ā that would advance a few tens of meters down the levee flank. Compound levees were thus a symptom of unsteady flow, where overflow levees were emplaced as relatively fast moving pāhoehoe sheets during pulses, and as slow-moving ‘a‘ā units during backup. Small, localized fluctuations in channel volume flux also occurred on timescales of minutes. Volumes of lava backed up behind blockages that formed at constrictions in the channel. Blockage collapse and/or enhanced flow under/around the blockage would then feed short-lived, wave-like, down-channel surges. Real fluctuations in channel volume flux, due to pulses and surges, can lead to significant errors in effusion rate calculations. Editorial responsibility: A. Woods     

Construction dynamics of a lava channel

We use a kinematic GPS and laser range finder survey of a 200 m-long section of the Muliwai a Pele lava channel (Mauna Ulu, Kilauea) to examine the construction processes and flow dynamics responsible for the channel–levee structure. The levees comprise three packages. The basal package comprises an 80–150 m wide ′a′a flow in which a ∼2 m deep and ∼11 m wide channel became centred. This is capped by a second package of thin (<45 cm thick) sheets of pahoehoe extending no more than 50 m from the channel. The upper-most package comprises localised ′a′a overflows. The channel itself contains two blockages located 130 m apart and composed of levee chunks veneered with overflow lava. The channel was emplaced over 50 h, spanning 30 May–2 June, 1974, with the flow front arriving at our section (4.4 km from the vent) 8 h after the eruption began. The basal ′a′a flow thickness yields effusion rates of 35 m³ s⁻¹ for the opening phase, with the initial flow advancing across the mapped section at ∼10 m/min. Short-lived overflows of fluid pahoehoe then built the levee cap, increasing the apparent channel depth to 4.8 m. There were at least six pulses at 90–420 m³ s⁻¹, causing overflow of limited extent lasting no more than 5 min. Brim-full flow conditions were thus extremely short-lived. During a dominant period of below-bank flow, flow depth was ∼2 m with an effusion rate of ∼35 m³ s⁻¹, consistent with the mean output rate (obtained from the total flow bulk volume) of 23–54 m³ s⁻¹. During pulses, levee chunks were plucked and floated down channel to form blockages. In a final low effusion rate phase, lava ponded behind the lower blockage to form a syn-channel pond that fed ′a′a overflow. After the end of the eruption the roofed-over pond continued to drain through the lower blockage, causing the roof to founder. Drainage emplaced inflated flows on the channel floor below the lower blockage for a further ∼10 h. The complex processes involved in levee–channel construction of this short-lived case show that care must be taken when using channel dimensions to infer flow dynamics. In our case, the full channel depth is not exposed. Instead the channel floor morphology reflects late stage pond filling and drainage rather than true channel-contained flow. Components of the compound levee relate to different flow regimes operating at different times during the eruption and associated with different effusion rates, flow dynamics and time scales. For example, although high effusion rate, brim-full flow was maintained for a small fraction of the channel lifetime, it emplaced a pile of pahoehoe overflow units that account for 60% of the total levee height. We show how time-varying volume flux is an important parameter in controlling channel construction dynamics. Because the complex history of lava delivery to a channel system is recorded by the final channel morphology, time-varying flow dynamics can be determined from the channel morphology. Developing methods for quantifying detailed flux histories for effusive events from the evidence in outcrop is therefore highly valuable. We here achieve this by using high-resolution spatial data for a channel system at Kilauea. This study not only indicates those physical and dynamic characteristics that are typical for basaltic lava flows on Hawaiian volcanoes, but also a methodology that can be widely applied to effusive basaltic eruptions.     

The morphology and evolution of the Stromboli 2002–2003 lava flow field: an example of a basaltic flow field emplaced on a steep slope

The use of a hand-held thermal camera during the 2002–2003 Stromboli effusive eruption proved essential in tracking the development of flow field structures and in measuring related eruption parameters, such as the number of active vents and flow lengths. The steep underlying slope on which the flow field was emplaced resulted in a characteristic flow field morphology. This comprised a proximal shield, where flow stacking and inflation caused piling up of lava on the relatively flat ground of the vent zone, that fed a medial–distal lava flow field. This zone was characterized by the formation of lava tubes and tumuli forming a complex network of tumuli and flows linked by tubes. Most of the flow field was emplaced on extremely steep slopes and this had two effects. It caused flows to slide, as well as flow, and flow fronts to fail frequently, persistent flow front crumbling resulted in the production of an extensive debris field. Channel-fed flows were also characterized by development of excavated debris levees in this zone (Calvari et al. 2005). Collapse of lava flow fronts and inflation of the upper proximal lava shield made volume calculation very difficult. Comparison of the final field volume with that expecta by integrating the lava effusion rates through time suggests a loss of ~70% erupted lava by flow front crumbling and accumulation as debris flows below sea level. Derived relationships between effusion rate, flow length, and number of active vents showed systematic and correlated variations with time where spreading of volume between numerous flows caused an otherwise good correlation between effusion rate, flow length to break down. Observations collected during this eruption are useful in helping to understand lava flow processes on steep slopes, as well as in interpreting old lava–debris sequences found in other steep-sided volcanoes subject to effusive activity.     

Cooling and crystallization of lava in open channels, and the transition of Pāhoehoe Lava to 'A'ā

 Samples collected from a lava channel active at Kīlauea Volcano during May 1997 are used to constrain rates of lava cooling and crystallization during early stages of flow. Lava erupted at near-liquidus temperatures (∼1150  °C) cooled and crystallized rapidly in upper parts of the channel. Glass geothermometry indicates cooling by 12–14  °C over the first 2 km of transport. At flow velocities of 1–2 m/s, this translates to cooling rates of 22–50  °C/h. Cooling rates this high can be explained by radiative cooling of a well-stirred flow, consistent with observations of non-steady flow in proximal regions of the channel. Crystallization of plagioclase and pyroxene microlites occurred in response to cooling, with crystallization rates of 20–50% per hour. Crystallization proceeded primarily by nucleation of new crystals, and nucleation rates of ∼10⁴/cm³s are similar to those measured in the 1984 open channel flow from Mauna Loa Volcano. There is no evidence for the large nucleation delays commonly assumed for plagioclase crystallization in basaltic melts, possibly a reflection of enhanced nucleation due to stirring of the flow. The transition of the flow surface morphology from pāhoehoe to 'a'ā occurred at a distance of 1.9 km from the vent. At this point, the flow was thermally stratified, with an interior temperature of ∼1137  °C and crystallinity of ∼15%, and a flow surface temperature of ∼1100  °C and crystallinity of ∼45%. 'A'ā formation initiated along channel margins, where crust was continuously disrupted, and involved tearing and clotting of the flow surface. Both observations suggest that the transition involved crossing of a rheological threshold. We suggest this threshold to be the development of a lava yield strength sufficient to prevent viscous flow of lava at the channel margin. We use this concept to propose that 'a'ā formation in open channels requires both sufficiently high strain rates for continued disruption of surface crusts and sufficient groundmass crystallinity to generate a yield strength equivalent to the imposed stress. In Hawai'i, where lava is typically microlite poor on eruption, these combined requirements help to explain two common observations on 'a'ā formation: (a) 'a'ā flow fields are generated when effusion rates are high (thus promoting crustal disruption); and (b) under most eruption conditions, lava issues from the vent as pāhoehoe and changes to 'a'ā only after flowing some distance, thus permitting sufficient crystallization. Received: 3 September 1998 / Accepted: 12 April 1999     

Viscous Newtonian laminar flow in a rectangular channel: application to Etna lava flows

 We introduce a 3D model for near-vent channelized lava flows. We assume the lava to be an isothermal Newtonian liquid flowing in a rectangular channel down a constant slope. The flow velocity is calculated with an analytical steady-state solution of the Navier-Stokes equation. The surface velocity and the flow rate are calculated as functions of the flow thickness for different flow widths, and the results are compared with those of a 2D model. For typical Etna lava flow parameters, the influence of levees on the flow dynamics is significant when the flow width is less than 25 m. The model predicts the volume flow rate corresponding to the surface velocity, taking into account that both depend on flow thickness. The effusion rate is a critical parameter to evaluate lava flow hazard. We propose a model to calculate the effusion rate given the lava flow width, the topograhic slope, the lava density, the surface flow velocity, and either the lava viscosity or the flow thickness. Received: 20 January 1998 / Accepted: 8 January 1999     

A laboratory model of surface crust formation and disruption on lava flows through non-uniform channels

Crust formation on basaltic lava flows dictates conditions of both flow cooling and emplacement. For this reason, flow histories are dramatically different depending on whether lava is transported through enclosed lava tubes or through open channels. Recent analog experiments in straight uniform channels (Griffiths et al. J Fluid Mech 496:33–62, 2003) have demonstrated that tube flow, dictated by a stationary surface crust, can be distinguished from a mobile crust regime, where a central solid crust is separated from channel walls by crust-free shear zones, by a simple dimensionless parameter ϑ, such that ϑ<25 produces tube flow and ϑ>25 describes the mobile crust regime. ϑ combines a previously determined parameter ψ, which describes the balance between the formation rate of surface solid and the shear strain that disrupts the solid crust, with the effects of thermal convection (described by the Rayleigh number Ra).Here we explore ways in which ϑ can be used to describe the behavior of basaltic lava channels. To do this we have extended the experimental approach to examine the effects of channel irregularities (expansions, contractions, sinuosity, and bottom roughness) on crust formation and disruption. We find that such changes affect local flow behavior and can thus change channel values of ϑ. For example, gradual widening of a channel results in a decrease in flow velocity that causes a decrease in ϑ and may allow a down-flow transition from the mobile crust to the tube regime. In contrast, narrowing of the channel causes an increase in flow velocity (increasing ϑ), thus inhibiting tube formation.We also quantify the fraction of surface covered by crust in the mobile crust regime. In shallow channels, variations in crust width (d _c) with channel width (W) are predicted to follow d _c∼W ^5/3. Analysis of channelized lava flows in Hawaii shows crustal coverage consistent with this theoretical result along gradually widening or narrowing channel reaches. An additional control on crustal coverage in both laboratory and basaltic flows is disruption of surface crust because of flow acceleration through constrictions, around bends, and over breaks in slope. Crustal breakage increases local rates of cooling and may cause local blockage of the channel, if crusts rotate and jam in narrow channel reaches. Together these observations illustrate the importance of both flow conditions and channel geometry on surface crust development and thus, by extension, on rates and mechanisms of flow cooling. Moreover, we note that this type of analysis could be easily extended through combined use of FLIR and LiDAR imaging to measure crustal coverage and channel geometry directly.Editorial responsibility: A. Harris     

The pit-craters and pit-crater-filling lavas of Masaya volcano

Lava flowing into a pit crater will become entrapped to form an inactive lava lake. At Masaya volcano (Nicaragua) pit filling lavas are exposed in the walls of Nindiri, Santiago and San Pedro pits. Mapping of these lavas shows that fill can involve emplacement of both ’a’a and pahoehoe, with single fill units ranging in thickness from 2 to 22 m. Thick units with columnar joints were emplaced as simple inactive lava lakes during high effusion rate episodes. Sequences of thinner units, which can form pit floor shields or compound lakes, were emplaced at lower effusion rates. Lava withdrawal caused unsupported sections of three 20-m-thick units to subside, resulting in unit flexure and faulting, and viscous peeling features reveal that subsidence occurred while at least one unit was still partially molten. Where withdrawal has not occurred, fill sequences are flat lying and symmetrically distributed around the feeder structures (cinder cones and dykes). The filled Nindiri pit holds 5 × 10⁷ m³ of lava in a 215-m-thick sequence. Partial fill of Santiago pit with 1 × 10⁷ m³ of lava has filled the pit with a 110-m-thick lava sequence, of which ∼50% has been consumed by formation of a secondary pit. Altogether, 6.4 × 10⁷ m³ of lava was erupted into Nindiri and Santiago during 1525–1965, with 94% of this volume remaining pit-contained; the remainder forms a north flank lava flow field. Pit development and filling is a dynamic and ephemeral process, having short-lived effects on volcano morphology, where pits develop and fill over hours-to-centuries. However, pits play an important role in shaping an edifice, representing lava sinks and controlling whether lavas are trapped or able to spread onto the flanks.     

Thermal structure and heat loss at the summit crater of an active lava dome

Forward-Looking Infrared (FLIR) nighttime thermal images were used to extract the thermal and morphological properties for the surface of a blocky-to-rubbley lava mass active within the summit crater of the Caliente vent at Santiaguito lava dome (Guatemala). Thermally the crater was characterized by three concentric regions: a hot outer annulus of loose fine material at 150–400°C, an inner cold annulus of blocky lava at 40–80°C, and a warm central core at 100–200°C comprising younger, hotter lava. Intermittent explosions resulted in thermal renewal of some surfaces, mostly across the outer annulus where loose, fine, fill material was ejected to expose hotter, underlying, material. Surface heat flux densities (radiative + free convection) were dominated by losses from the outer annulus (0.3–1.5 × 10⁴  s⁻¹m⁻²), followed by the hot central core (0.1–0.4 × 10⁴ J s⁻¹m⁻²) and cold annulus (0.04–0.1 × 10⁴ J s⁻¹m⁻²). Overall surface power output was also dominated by the outer annulus region (31–176 MJ s⁻¹), but the cold annulus contributed equal power (2.41–7.07 MJ s⁻¹) as the hot central core (2.68–6.92 MJ s⁻¹) due to its greater area. Cooled surfaces (i.e. the upper thermal boundary layer separating surface temperatures from underlying material at magmatic temperatures) across the central core and cold annulus had estimated thicknesses, based on simple conductive model, of 0.3–2.2 and 1.5–4.3 m. The stability of the thermal structure through time and between explosions indicates that it is linked to a deeper structural control likely comprising a central massive plug, feeding lava flow from the SW rim of the crater, surrounded by an arcuate, marginal fracture zone through which heat and mass can preferentially flow.     

Three‐dimensional flow structure and patterns of bed shear stress in an evolving compound meander bend

Compound meander bends with multiple lobes of maximum curvature are common in actively evolving lowland rivers. Interaction among spatial patterns of mean flow, turbulence, bed morphology, bank failures and channel migration in compound bends is poorly understood. In this paper, acoustic Doppler current profiler (ADCP) measurements of the three‐dimensional (3D) flow velocities in a compound bend are examined to evaluate the influence of channel curvature and hydrologic variability on the structure of flow within the bend. Flow structure at various flow stages is related to changes in bed morphology over the study timeframe. Increases in local curvature within the upstream lobe of the bend reduce outer bank velocities at morphologically significant flows, creating a region that protects the bank from high momentum flow and high bed shear stresses. The dimensionless radius of curvature in the upstream lobe is one‐third less than that of the downstream lobe, with average bank erosion rates less than half of the erosion rates for the downstream lobe. Higher bank erosion rates within the downstream lobe correspond to the shift in a core of high velocity and bed shear stresses toward the outer bank as flow moves through the two lobes. These erosion patterns provide a mechanism for continued migration of the downstream lobe in the near future. Bed material size distributions within the bend correspond to spatial patterns of bed shear stress magnitudes, indicating that bed material sorting within the bend is governed by bed shear stress. Results suggest that patterns of flow, sediment entrainment, and planform evolution in compound meander bends are more complex than in simple meander bends. Moreover, interactions among local influences on the flow, such as woody debris, local topographic steering, and locally high curvature, tend to cause compound bends to evolve toward increasing planform complexity over time rather than stable configurations. Copyright © 2016 John Wiley & Sons, Ltd.     

The character of long-term eruptions: inferences from episodes 50–53 of the Pu'u 'Ō'ō-Kūpaianaha eruption of Kīlauea Volcano

 The Pu'u 'Ō'ō-Kūpaianaha eruption on the east rift zone of Kīlauea began in January 1983. The first 9 years of the eruption were divided between the Pu'u 'Ō'ō (1983–1986) and Kūpaianaha (1986–1992) vents, each characterized by regular, predictable patterns of activity that endured for years. In 1990 a series of pauses in the activity disturbed the equilibrium of the eruption, and in 1991, the output from Kūpaianaha steadily declined and a short-lived fissure eruption broke out between Kūpaianaha and Pu'u 'Ō'ō. In February 1992 the Kūpaianaha vent died, and, 10 days later, eruptive episode 50 began as a fissure opened on the uprift flank of the Pu'u 'Ō'ō cone. For the next year, the eruption was marked by instability as more vents opened on the flank of the cone and the activity was repeatedly interrupted by brief pauses in magma supply to the vents. Episodes 50–53 constructed a lava shield 60 m high and 1.3 km in diameter against the steep slope of the Pu'u 'Ō'ō cone. By 1993 the shield was pockmarked by collapse pits as vents and lava tubes downcut as much as 29 m through the thick deposit of scoria and spatter that veneered the cone. As the vents progressively lowered, the level of the Pu'u 'Ō'ō pond also dropped, demonstrating the hydraulic connection between the two. The downcutting helped to undermine the prominent Pu'u 'Ō'ō cone, which has diminished in size both by collapse, as a large pit crater formed over the conduit, and by burial of its flanks. Intervals of eruptive instability, such as that of 1991–1993, accelerate lateral expansion of the subaerial flow field both by producing widely spaced vents and by promoting surface flow activity as lava tubes collapse and become blocked during pauses. Received: 1 July 1997 / Accepted: 23 October 1997     

Thick lava flows of Karisimbi Volcano, Rwanda: insights from SIR-C interferometric topography

 We use a digital elevation model (DEM) derived from interferometrically processed SIR-C radar data to estimate the thickness of massive trachyte lava flows on the east flank of Karisimbi Volcano, Rwanda. The flows are as long as 12 km and average 40–60 m (up to >140 m) in thickness. By calculating and subtracting a reference surface from the DEM, we derived a map of flow thickness, which we used to calculate the volume (up to 1 km³ for an individual flow, and 1.8 km³ for all the identified flows) and yield strength of several flows (23–124 kPa). Using the DEM we estimated apparent viscosity based on the spacing of large folds (1.2×10¹² to 5.5×10¹² Pa s for surface viscosity, and 7.5×10¹⁰ to 5.2×10¹¹ Pa s for interior viscosity, for a strain interval of 24 h). We use shaded-relief images of the DEM to map basic flow structures such as channels, shear zones, and surface folds, as well as flow boundaries. The flow thickness map also proves invaluable in mapping flows where flow boundaries are indistinct and poorly expressed in the radar backscatter and shaded-relief images. Received: 6 September 1997 / Accepted: 15 May 1998     

Viscosity of microlite-bearing rhyolitic obsidians: an experimental study

 To investigate the influence of microlites on lava flow rheology, the viscosity of natural microlite-bearing rhyolitic obsidians of calc-alkaline and peralkaline compositions containing 0.1–0.4 wt.% water was measured at volcanologically relevant temperatures (650–950  °C), stresses (10³–10⁵ Pa) and strain rates (10^–5 to 10^–7 s^–1). The glass transition temperatures (T _g ) were determined from scanning calorimetric measurements on the melts for a range of cooling/heating rates. Based on the equivalence of enthalpic (calorimetric) and shear (viscosity) relaxation, we calculated the viscosity of the melt in crystal-bearing samples from the T _g data. The difference between the calculated viscosity of the melt phase and the measured viscosity for the crystal-bearing samples is interpreted to be the physical effect of microlites on the measured viscosity. The effect of <5 vol.% rod-like microlites on the melt rheology is negligible. Microlite-rich and microlite-poor samples from the same lava flow and with identical bulk chemistry show a difference of 0.6 log₁₀ units viscosity (Pa s), interpreted to be due to differences in melt chemistry caused by the presence of microlites. The only major differences between measured and calculated viscosities were for two samples: a calc-alkaline rhyolite with 1 vol.% branching crystals, and a peralkaline rhyolite containing crystal-rich bands with >45 vol.% crystals. For both of these samples a connectivity factor is apparent, with, for the latter, a close packing framework of crystals which is interpreted to influence the apparent viscosity. Received: 14 March 1996 / Accepted: 30 May 1996     

Lava flow hazard at Nyiragongo volcano,D.R.C.

The 2002 eruption of Nyiragongo volcano constitutes the most outstanding case ever of lava flow in a big town. It also represents one of the very rare cases of direct casualties from lava flows, which had high velocities of up to tens of kilometer per hour. As in the 1977 eruption, which is the only other eccentric eruption of the volcano in more than 100 years, lava flows were emitted from several vents along a N–S system of fractures extending for more than 10 km, from which they propagated mostly towards Lake Kivu and Goma, a town of about 500,000 inhabitants. We assessed the lava flow hazard on the entire volcano and in the towns of Goma (D.R.C.) and Gisenyi (Rwanda) through numerical simulations of probable lava flow paths. Lava flow paths are computed based on the steepest descent principle, modified by stochastically perturbing the topography to take into account the capability of lava flows to override topographic obstacles, fill topographic depressions, and spread over the topography. Code calibration and the definition of the expected lava flow length and vent opening probability distributions were done based on the 1977 and 2002 eruptions. The final lava flow hazard map shows that the eastern sector of Goma devastated in 2002 represents the area of highest hazard on the flanks of the volcano. The second highest hazard sector in Goma is the area of propagation of the western lava flow in 2002. The town of Gisenyi is subject to moderate to high hazard due to its proximity to the alignment of fractures active in 1977 and 2002. In a companion paper (Chirico et al., Bull Volcanol, in this issue, 2008) we use numerical simulations to investigate the possibility of reducing lava flow hazard through the construction of protective barriers, and formulate a proposal for the future development of the town of Goma.     

Voluminous lava flows at Oldoinyo Lengai in 2006: chronology of events and insights into the shallow magmatic system

The largest natrocarbonatite lava flow eruption ever documented at Oldoinyo Lengai, NW Tanzania, occurred from March 25 to April 5, 2006, in two main phases. It was associated with hornito collapse, rapid extrusion of lava covering a third of the crater and emplacement of a 3-km long compound rubbly pahoehoe to blocky aa-like flow on the W flank. The eruption was followed by rapid enlargement of a pit crater. The erupted natrocarbonatite lava has high silica content (3% SiO₂). The eruption chronology is reconstructed from eyewitness and news media reports and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data, which provide the most reliable evidence to constrain the eruption’s onset and variations in activity. The eruption products were mapped in the field and the total erupted lava volume estimated at 9.2 ± 3.0 × 10⁵ m³. The event chronology and field evidence are consistent with vent construct instability causing magma mixing and rapid extrusion from shallow reservoirs. It provides new insights into and highlights the evolution of the shallow magmatic system at this unique natrocarbonatite volcano.     

Emplacement of subaerial pahoehoe lava sheet flows into water: 1990 Kūpaianaha flow of Kilauea volcano at Kaimū Bay, Hawai`i

Episode 48 of the ongoing eruption of Kilauea, Hawai`i, began in July 1986 and continuously extruded lava for the next 5.5 years from a low shield, Kūpaianaha. The flows in March 1990 headed for Kalapana and inundated the entire town under 15–25 m of lava by the end of August. As the flows advanced eastward, they entered into Kaimū Bay, replacing it with a plain of lava that extends 300 m beyond the original shoreline. The focus of our study is the period from August 1 to October 31, 1990, when the lava buried almost 406,820 m² of the 5-m deep bay. When lava encountered the sea, it flowed along the shoreline as a narrow primary lobe up to 400 m long and 100 m wide, which in turn inflated to a thickness of 5–6 m. The flow direction of the primary lobes was controlled by the submerged delta below the lavas and by damming up lavas fed at low extrusion rates. Breakout flows through circumferential and axial inflation cracks on the inflating primary lobes formed smaller secondary lobes, burying the lows between the primary lobes and hiding their original outlines. Inflated flow lobes eventually ruptured at proximal and/or distal ends as well as mid-points between the two ends, feeding new primary lobes which were emplaced along and on the shore side of the previously inflated lobes. The flow lobes mapped with the aid of aerial photographs were correlated with daily observations of the growing flow field, and 30 primary flow lobes were dated. Excluding the two repose periods that intervened while the bay was filled, enlargement of the flow field took place at a rate of 2,440–22,640 square meters per day in the bay. Lobe thickness was estimated to be up to 11 m on the basis of cross sections of selected lobes measured using optical measurement tools, measuring tape and hand level. The total flow-lobe volume added in the bay during August 1–October 31 was approximately 3.95 million m³, giving an average supply rate of 0.86 m³/s.     

Lava penetrating water: the different behaviours of pāhoehoe and ‘a‘ā at the Nesjahraun, Tingvellir, Iceland

The Nesjahraun is a basaltic lava flow erupted from a subaerial fissure, extending NE along the Tingvellir graben from the Hengill central volcano that produced pāhoehoe lava followed by ‘a‘ā. The Nesjahraun entered Iceland’s largest lake, Tingvallavatn, along its southern shore during both phases of the eruption and exemplifies lava flowing into water in a lacustrine environment in the absence of powerful wave action. This study combines airborne light detection and ranging, sidescan sonar and Chirp seismic data with field observations to investigate the behaviour of the lava as it entered the water. Pāhoehoe sheet lava was formed during the early stages of the eruption. Along the shoreline, stacks of thin (5–20 cm thick), vesicular, flows rest upon and surround low (<5 m) piles of coarse, unconsolidated, variably oxidised spatter. Clefts within the lava run inland from the lake. These are 2–5 m wide, >2 m deep, ∼50 m long, spaced ∼50 m apart and have sub-horizontal striations on the walls. They likely represent channels or collapsed tubes along which lava was delivered into the water. A circular rootless cone, Eldborg, formed when water infiltrated a lava tube. Offshore from the pāhoehoe lavas, the gradient of the flow surface steepens, suggesting a change in flow regime and the development of a talus ramp. Later, the flow was focused into a channel of ‘a‘ā lava, ∼200–350 m wide. This split into individual flow lobes 20–50 m wide along the shore. ‘A‘ā clinker is exposed on the water’s edge, as well as glassy sand and gravel, which has been locally intruded by small (<1 m), irregularly shaped, lava bodies. The cores of the flow lobes contain coherent, but hackly fractured lava. Mounds consisting predominantly of scoria lapilli and the large paired half-cone of Grámelur were formed in phreatomagmatic explosions. The ‘a‘ā flow can be identified underwater over 1 km offshore, and the sidescan data suggest that the flow lobes remained coherent flowing down a gradient of <10°. The Nesjahraun demonstrates that, even in the absence of ocean waves, phreatomagmatic explosions are ubiquitous and that pāhoehoe flows are much more likely to break up on entering the water than ‘a‘ā flows, which, with a higher flux and shallow underlying surface gradient, can penetrate water and remain coherent over distances of at least 1 km.     

Formation of submarine flat-topped volcanic cones in Hawai'i

 High-resolution bathymetric mapping has shown that submarine flat-topped volcanic cones, morphologically similar to ones on the deep sea floor and near mid-ocean ridges, are common on or near submarine rift zones of Kilauea, Kohala (or Mauna Kea), Mahukona, and Haleakala volcanoes. Four flat-topped cones on Kohala were explored and sampled with the Pisces V submersible in October 1998. Samples show that flat-topped cones on rift zones are constructed of tholeiitic basalt erupted during the shield stage. Similarly shaped flat-topped cones on the northwest submarine flank of Ni'ihau are apparently formed of alkalic basalt erupted during the rejuvenated stage. Submarine postshield-stage eruptions on Hilo Ridge, Mahukona, Hana Ridge, and offshore Ni'ihau form pointed cones of alkalic basalt and hawaiite. The shield stage flat-topped cones have steep (∼25°) sides, remarkably flat horizontal tops, basal diameters of 1–3 km, and heights <300 m. The flat tops commonly have either a low mound or a deep crater in the center. The rejuvenated-stage flat-topped cones have the same shape with steep sides and flat horizontal tops, but are much larger with basal diameters up to 5.5 km and heights commonly greater than 200 m. The flat tops have a central low mound, shallow crater, or levees that surrounded lava ponds as large as 1 km across. Most of the rejuvenated-stage flat-topped cones formed on slopes <10° and formed adjacent semicircular steps down the flank of Ni'ihau, rather than circular structures. All the flat-topped cones appear to be monogenetic and formed during steady effusive eruptions lasting years to decades. These, and other submarine volcanic cones of similar size and shape, apparently form as continuously overflowing submarine lava ponds. A lava pond surrounded by a levee forms above a sea-floor vent. As lava continues to flow into the pond, the lava flow surface rises and overflows the lowest point on the levee, forming elongate pillow lava flows that simultaneously build the rim outward and upward, but also dam and fill in the low point on the rim. The process repeats at the new lowest point, forming a circular structure with a flat horizontal top and steep pillowed margins. There is a delicate balance between lava (heat) supply to the pond and cooling and thickening of the floating crust. Factors that facilitate construction of such landforms include effusive eruption of lava with low volatile contents, moderate to high confining pressure at moderate to great ocean depth, long-lived steady eruption (years to decades), moderate effusion rates (probably ca. 0.1 km³/year), and low, but not necessarily flat, slopes. With higher effusion rates, sheet flows flood the slope. With lower effusion rates, pillow mounds form. Hawaiian shield-stage eruptions begin as fissure eruptions. If the eruption is too brief, it will not consolidate activity at a point, and fissure-fed flows will form a pond with irregular levees. The pond will solidify between eruptive pulses if the eruption is not steady. Lava that is too volatile rich or that is erupted in too shallow water will produce fragmental and highly vesicular lava that will accumulate to form steep pointed cones, as occurs during the post-shield stage. The steady effusion of lava on land constructs lava shields, which are probably the subaerial analogs to submarine flat-topped cones but formed under different cooling conditions. Received: 30 September 1999 / Accepted: 9 March 2000     

Simulations of the 2004 lava flow at Etna volcano using the magflow cellular automata model

Since the mechanical properties of lava change over time, lava flows represent a challenge for physically based modeling. This change is ruled by a temperature field which needs to be modeled. MAGFLOW Cellular Automata (CA) model was developed for physically based simulations of lava flows in near real-time. We introduced an algorithm based on the Monte Carlo approach to solve the anisotropic problem. As transition rule of CA, a steady-state solution of Navier-Stokes equations was adopted in the case of isothermal laminar pressure-driven Bingham fluid. For the cooling mechanism, we consider only the radiative heat loss from the surface of the flow and the change of the temperature due to mixture of lavas between cells with different temperatures. The model was applied to reproduce a real lava flow that occurred during the 2004–2005 Etna eruption. The simulations were computed using three different empirical relationships between viscosity and temperature.     

Lava channel roofing,overflows, breaches and switching: insights from the 2008–2009 eruption of Mt. Etna

During long-lived basaltic eruptions, overflows from lava channels and breaching of channel levées are important processes in the development of extensive 'a'ā lava flow-fields. Short-lived breaches result in inundation of areas adjacent to the main channel. However, if a breach remains open, lava supply to the original flow front is significantly reduced, and flow-field widening is favoured over lengthening. The development of channel breaches and overflows can therefore exert strong control over the overall flow-field development, but the processes that determine their location and frequency are currently poorly understood. During the final month of the 2008–2009 eruption of Mt. Etna, Sicily, a remote time-lapse camera was deployed to monitor events in a proximal region of a small ephemeral lava flow. For over a period of ~10 h, the flow underwent changes in surface elevation and velocity, repeated overflows of varying vigour and the construction of a channel roof (a required prelude to lava tube formation). Quantitative interpretation of the image sequence was facilitated by a 3D model of the scene constructed using structure-from-motion computer vision techniques. As surface activity waned during the roofing process, overflow sites retreated up the flow towards the vent, and eventually, a new flow was initiated. Our observations and measurements indicate that flow surface stagnation and flow inflation propagated up-flow at an effective rate of ~6 m h⁻¹, and that these processes, rather than effusion rate variations, were ultimately responsible for the most vigorous overflow events. We discuss evidence for similar controls during levée breaching and channel switching events on much larger flows on Etna, such as during the 2001 eruption.     

Widespread inflation and drainage of a pāhoehoe flow field: the Nesjahraun, Tingvellir, Iceland

This study describes the emplacement of the Nesjahraun, a basaltic lava flow that entered the lake Tingvallavatn, SW Iceland. High-resolution remotely sensed data were combined with fieldwork to map the flow field. Onshore, the Nesjahraun exhibits a variety of textures related to the widespread inflation and collapse of a pāhoehoe flow field. Its emplacement is interpreted as follows: Initially, the eruption produced sheet pāhoehoe. In the central part of the flow field, the lava has a platy-ridged surface, which is similar to some other lava flows in Iceland and on Mars. Here, the texture is interpreted to have formed by unsteady inflation of the brittle crust of stationary sheet pāhoehoe, causing it to break into separate plates. The ridges of broken pāhoehoe slabs formed as the plates of crust moved vertically past each other in a process similar to the formation of shatter rings. Upstream, fresh lava overflowed repeatedly from channels and tubes, covering the surface with shelly pāhoehoe. Formation of a 250-m-wide open channel through the flow field allowed the inflated central part of the flow to drain rapidly. This phase produced ‘a‘ā lava, which eroded the channel walls, carrying broken pāhoehoe slabs, lava balls and detaching large (>200 m long) rafts of compound shelly pāhoehoe lava. Much of this channelized lava flowed into the lake, leaving a network of drained channels and tubes in the upstream part of the flow. As in other locations, the platy-ridged texture is associated with a low underlying slope and high eruption rate. Here, its formation was possibly enhanced by lateral confinement, hindered entry into the lake and an elevated vent location. We suggest that formation of this type of platy-ridged lava, where the plates are smooth and the ridges are slabs of broken pāhoehoe, can occur without significant horizontal transport, as the surface crust is broken into plates in situ. This reconstruction of the emplacement of the Nesjahraun also demonstrates that high-resolution aerial survey data are extremely useful in the mapping and measurement of lithofacies distributions in large flow fields, but that fieldwork is still necessary to obtain the detailed textural information necessary to interpret them.