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.     

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     

A laboratory analog study of the surface morphology of lava flows extruded from point and line sources

In laboratory experiments designed to model lava flow processes, liquid polyethylene glycol wax is forced through either a small hole or a long narrow slit onto the base of a tank of cold water, where it spreads laterally while cooling and solidifying at its surface. We observe the surface structure of the flow, and its dependence on the flow rate, thermal conditions and basal roughness. In each case, solidification of a crust during spreading gives rise to a number of different surface morphologies, each of which forms under a restricted range of conditions. The dominant morphologies, referred to as “pillows”, “rifts”, transverse folds, and marginal levees, correspond to features observed on natural lava flows. Results for radial spreading over a rough base and for spreading in two directions from the line source over both smooth and rough bases complement those reported earlier for radial spreading of solidifying wax from a point source on a smooth base. Together they indicate a robust dependence of morphology on the distance from the vent at which solid crust begins to form. This distance is, in turn, determined by the extrusion rate, the rheology of the liquid wax (or magma), the reduced gravity, the magnitude of the surface heat flux, and the amount of cooling required to solidify the flow surface. The results also indicate factors influencing the distribution of crust and its deformation, and may provide a means by which observations of surface morphology can be used to place constraints on the emplacement conditions of lava flows.     

Effects of eruption history and cooling rate on lava dome growth

To better understand the factors controlling the shapes of lava domes, laboratory simulations, measurements from active and prehistoric flows and dimensional analysis were used to explore how effusion history and cooling rate affect the final geometry of a dome. Fifty experiments were conducted in which a fixed volume of polyethylene glycol wax was injected into a tank of cold sucrose solution, either as one continuous event or as a series of shorter pulses separated by repose periods. When the wax cooling rates exceeded a critical minimum value, the dome aspect ratios (height/diameter) increased steadily with erupted volume over the course of a single experiment and the rate at which height increased with volume depended linearly on the time-averaged effusion rate. Thus the average effusion rate could be estimated from observations of how the dome shape changed with time. Our experimental results and dimensional analyses were compared with several groups of natural lava flows: the recently emplaced Mount St Helens and Soufrière domes, which had been carefully monitored while active; three sets of prehistoric rhyolite domes that varied in eruptive style and shape; and two sets of Holocene domes with similar shapes, but different compositions. Geometric measurements suggest that dome morphology can be directly correlated with effusion rate for domes of similar composition from the same locality, and that shape alone can be related to a dimensionless number comparing effusion rate and cooling rate. Extrapolation to the venusian pancake domes suggests that they formed from relatively viscous lavas extruded either episodically or at average effusion rates low enough to allow solidified surface crust to exert a dominating influence on the final morphology.     

Factors controlling the lengths of channel-fed lava flows

Factors which control lava flow length are still not fully understood. The assumption that flow length as mainly influenced by viscosity was contested by Walker (1973) who proposed that the length of a lava flow was dependent on the mean effusion rate, and by Malin (1980) who concluded that flow length was dependent on erupted volume. Our reanalysis of Malin's data shows that, if short duration and tube-fed flows are eliminated, Malin's Hawaiian flow data are consistent with Walker's assertion. However, the length of a flow can vary, for a given effusion rate, by a factor of 7, and by up to 10 for a given volume. Factors other than effusion rate and volume are therefore clearly important in controlling the lengths of lava flows. We establish the relative importance of the other factors by performing a multivariate analysis of data for recent Hawaiian lava flows. In addition to generating empirical equations relating flow length to other variables, we have developed a non-isothermal Bingham flow model. This computes the channel and levee width of a flow and hence permits the advance rates of flows and their maximum cooling-limited lengths for different gradients and effusion rates to be calculated. Changing rheological properties are taken into account using the ratio of yield strength to viscosity; available field measurements show that this varies systematically from the vent to the front of a lava flow. The model gives reasonable agreement with data from the 1983–1986 Pu'u Oo eruptions and the 1984 eruption of Mauna Loa. The method has also been applied to andesitic and rhyolitic lava flows. It predicts that, while the more silicic lava flows advance at generally slower rates than basaltic flows, their maximum flow lengths, for a given effusion rate, will be greater than for basaltic lava flows.     

The 1975 sub-terminal lavas, mount etna: a case history of the formation of a compound lava field

The 1975 sub-terminal activity was characterised by low effusion rates (0.3–0.5 m³ s⁻¹) and the formation of a compound lava field composed of many thousands of flow units. Several boccas were active simultaneously and effusion rates from individual boccas varied from about 10⁻⁴ to 0.25 m³s⁻¹. The morphology of lava flows was determined by effusion rate (E): aa flows with well-developed channels and levees formed when E > 2 × 10⁻³ m³ s⁻¹, small pahoehoe flows formed when 2 × 10⁻³ m³ s⁻¹ >E > 5 > 10⁻⁴ m³ s⁻¹ and pahoehoe toes formed when E < 5 × 10⁻⁴ m³ s⁻¹. There was very little variation with time in the effusion temperature, composition or phenocryst content of the lava.New boccas were commonly formed at the fronts of mature lava flows which had either ceased to flow or were moving slowly. These secondary boccas developed when fluid lava in the interior of mature aa flows either found a weakness in the flow front or was exposed by avalanching of the moving flow front. The resulting release of fluid lava was accompanied by either partial drainage of the mature flow or by the formation of a lava tube in the parent flow. The temperature of the lava forming the new bocca decreased with increasing distance from the source bocca (0.035°C m⁻¹). It is demonstrated from the rate of temperature decrease and from theoretical considerations that many of the Etna lavas still contained a substantial proportion of uncooled material in their interior as they came to rest. The formation of secondary boccas is postulated to be one reason why direct measurements of effusion rates tend, in general, to overestimate the total effusion rates of sub-terminal Etna lava fields.     

Rootless shield and perched lava pond collapses at Kīlauea Volcano,Hawai'i

Effusion rate is a primary measurement used to judge the expected advance rate, length, and hazard potential of lava flows. At basaltic volcanoes, the rapid draining of lava stored in rootless shields and perched ponds can produce lava flows with much higher local effusion rates and advance velocities than would be expected based on the effusion rate at the vent. For several months in 2007–2008, lava stored in a series of perched ponds and rootless shields on Kīlauea Volcano, Hawai'i, was released episodically to produce fast-moving 'a'ā lava flows. Several of these lava flows approached Royal Gardens subdivision and threatened the safety of remaining residents. Using time-lapse image measurements, we show that the initial time-averaged discharge rate for one collapse-triggered lava flow was approximately eight times greater than the effusion rate at the vent. Though short-lived, the collapse-triggered 'a'ā lava flows had average advance rates approximately 45 times greater than that of the pāhoehoe flow field from which they were sourced. The high advance rates of the collapse-triggered lava flows demonstrates that recognition of lava accumulating in ponds and shields, which may be stored in a cryptic manner, is vital for accurately assessing short-term hazards at basaltic volcanoes.     

Development of the 1990 Kalapana Flow Field,Kilauea Volcano,Hawaii

The 1990 Kalapana flow field is a complex patchwork of tube-fed pahoehoe flows erupted from the Kupaianaha vent at a low effusion rate (approximately 3.5 m³/s). These flows accumulated over an 11-month period on the coastal plain of Kilauea Volcano, where the pre-eruption slope angle was less than 2°. the composite field thickened by the addition of new flows to its surface, as well as by inflation of these flows and flows emplaced earlier. Two major flow types were identified during the development of the flow field: large primary flows and smaller breakouts that extruded from inflated primary flows. Primary flows advanced more quickly and covered new land at a much higher rate than breakouts. The cumulative area covered by breakouts exceeded that of primary flows, although breakouts frequently covered areas already buried by recent flows. Lava tubes established within primary flows were longer-lived than those formed within breakouts and were often reoccupied by lava after a brief hiatus in supply; tubes within breakouts were never reoccupied once the supply was interrupted. During intervals of steady supply from the vent, the daily areal coverage by lava in Kalapana was constant, whereas the forward advance of the flows was sporadic. This implies that planimetric area, rather than flow length, provides the best indicator of effusion rate for pahoehoe flow fields that form on lowangle slopes.     

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     

Morphological complexities and hazards during the emplacement of channel-fed `a`ā lava flow fields: A study of the 2001 lower flow field on Etna

Long-lived basaltic eruptions often produce structurally complex, compound `a`ā flow fields. Here we reconstruct the development of a compound flow field emplaced during the 2001 eruption of Mt. Etna (Italy). Following an initial phase of cooling-limited advance, the reactivation of stationary flows by superposition of new units caused significant channel drainage. Later, blockages in the channel and effusion rate variations resulted in breaching events that produced two new major flow branches. We also examined small-scale, late-stage ‘squeeze-up’ extrusions that were widespread in the flow field. We classified these as ‘flows’, ‘tumuli’ or ‘spines’ on the basis of their morphology, which depended on the rheology, extrusion rate and cooling history of the lava. Squeeze-up flows were produced when the lava was fluid enough to drain away from the source bocca, but fragmented to produce blade-like features that differed markedly from `a`ā clinker. As activity waned, increased cooling and degassing led to lava arriving at boccas with a higher yield strength. In many cases this was unable to flow after extrusion, and laterally extensive, near-vertical sheets of lava developed. These are considered to be exogenous forms of tumuli. In the highest yield strength cases, near-solid lava was extruded from the flow core as a result of ramping, forming spines. The morphology and location of the squeeze-ups provides insight into the flow rheology at the time of their formation. Because they represent the final stages of activity of the flow, they may also help to refine estimates of the most advanced rheological states in which lava can be considered to flow. Our observations suggest that real-time monitoring of compound flow field evolution may allow complex processes such as channel breaching and bocca formation to be forecast. In addition, documenting the occurrence and morphology of squeeze-ups may allow us to determine whether there is any risk of a stalled flow front being reactivated. This will therefore enhance our ability to track and assess hazard posed by lava flow emplacement.     

Effects of regional slope on viscous flows: a preliminary study of lava terrace emplacement at submarine volcanic rift zones

To understand how large submarine lava terraces form and why they are not commonly observed on land, we developed an isoviscous gravity flow model on an inclined surface to simulate the evolution and emplacement of lava flows under submarine conditions. By solving this preliminary model using a finite difference method, we are able to quantify how lava viscosity, pre-existing topographic slope, effusion rate, and lava volume affect meso-scale lava morphology. Our simulations show that, in general, high lava viscosity, gentle regional slope, and low effusion rate favor the formation of large terraces, but environmental conditions also play an important role. A gravity flow spreads more slowly underwater than subaerially. We also conclude that for low viscosity basaltic lava, the cooling of the lava body is one of the most critical factors that affect its shape. This study shows that the isoviscous model, though oversimplified, provides a quantitative tool to relate lava morphology to eruption characteristics. To gain a better understanding of the controls on submarine lava terrace formation, future models must take into account the temporal and spatial variation of lava viscosity, especially the effects of a brittle outer shell.     

Tumulus development on lava flows: insights from observations of active tumuli and analysis of formation models

Here, we use observations of active flows along with detailed morphometric field measurements of more than 70 tumuli on flows at Mount Etna (Italy), Kilauea, and Hualalai (US) volcanoes to constrain a previously published model that estimates the pressure needed to form tumuli. In an attempt to discover the nature and magnitude of pressure variations within active lava flow interiors, we then consider how tumuli differ from idealized circular plates. We incorporate observations of active tumuli and find that they may grow asymmetrically yet produce a symmetrical tumulus and can form where the flow path significantly changes direction. Bending models of clamped edges provide the most reasonable head estimates for the tumuli in our study. Tumulus formation requires the proper combination of cooling and effusion rate. If cooling is too extensive and effusion rate too low, the crust will provide too much resistance to bending. If cooling is too limited and effusion rates too high, crusts will not develop or have insufficient strength to resist fracture and subsequent breakouts. We do not find it surprising that tumuli are rarely found over well-established lava tubes that typically have rigid, walls/overlying crusts that exceed 2 m in thickness and provide too much resistance to bending. Silicic flows lack tumuli because the viscosity gradients within the flow are insufficient to concentrate stress in a localized area.     

The emplacement of pahoehoe toes: field observations and comparison to laboratory simulations

We observed active pahoehoe lobes erupted on Kilauea during May-June 1996, and found a range of emplacement styles associated with variations in local effusion rate, flow velocity, and strain rate. These emplacement styles were documented and quantified for comparison with earlier laboratory experiments.At the lowest effusion rates, velocities, and strain rates, smooth-surfaced lobes were emplaced via swelling, where new crust formed along an incandescent lip at the front of the lobe and the rest of the lobe was covered with a dark crust. At higher effusion rates, strain rates and velocities, lobes were emplaced through tearing or cracking. Tearing was characterized by ripping of the ductile crust near the initial breakout point, and most of the lobe surface was incandescent during its emplacement. This mechanism was observed to generate both smooth-surfaced lobes, and, when the lava encountered an obstacle, folded lobes. Cracking lobes were similar to those emplaced via tearing, but involved breaking of a thicker, brittle crust at the initial breakout of the lobe and therefore required somewhat higher flow rates than did tearing. Cracking lobes typically formed ropy folds in the center of the lobe, and smooth margins. At the highest effusion rates, strain rates, and flow velocities, the lava formed open channels with distinct levees.The final lobe morphologies were compared to results from laboratory simulations, which were designed to infer effusion rate from final flow morphology, to quantitatively test the laboratory results on the scale of individual natural pahoehoe lobes. There is general agreement between results from laboratory simulations and natural lavas on the scale of individual pahoehoe lobes, but there are disparities between laboratory flows and lava flows on the scale of an entire pahoehoe lava flow field.Editorial responsibility: A. Woods     

Slabby pahoehoe from the western Deccan Volcanic Province: evidence for incipient pahoehoe–aa transitions

The pahoehoe–aa transition for a flow exposed near Bodshil village from the western part of the Deccan Volcanic Province (DVP) is reported for the first time. The 1-km-long Bodshil flow issued as a small sheet from a pre-existing lobe. Near the source, the crust is characterised by numerous squeeze-ups. A number of gaping fractures, parallel to sub-parallel to the flow direction, are exposed on the surface in the medial portion of the flow. About 800 m away, the flow completely transforms to slabby pahoehoe. The terminal portion of the flow is characterised by concentrations of slabs, blocks and lava balls. The size and concentrations of the slabs and lava balls appear to increase along the length of the flow. Petrographic studies reveal a dominant hypohyaline texture. The flow core is coarse and is characterised by plagioclase set in a glassy matrix. The presence of clinopyroxene in addition to plagioclase and glass distinguishes the crust and interslab crust from the core. On the basis of mineralogy, a temperature range of 1146±15°C to 1169±15°C is inferred for the Bodshil flow. Increased vesicle deformation across the transition is discernible and an average D-value of <0.4 indicates moderate strain rates during emplacement. In light of the morphology and petrography, the cooling history and the mode of emplacement of the Bodshil flow is discussed. The flow originated as a small toe at the leading edge of a pahoehoe flow, and grew into a sheet by the mechanism of inflation. Continuous inflation caused the brittle crust to uplift and produce a network of inflation clefts that were subsequently occupied by squeeze-ups. Temporary stagnation of the flow due to cessation of lava supply or storage allowed the crust to grow and thicken. Renewed movement of the stored and cooled lava to the flow front at a fairly high volumetric rate was responsible for the initial disruption of the crust. High rates of crustal disruption induced higher rates of degassing and cooling, which resulted in rapid crystallisation of the fluid core. Increase in crystallinity lead to the onset of yield strength, and it is envisaged that at least the terminal parts of the flow behaved as a Bingham fluid. The Bodshil flow is unique to the DVP because it is the first to record slabby pahoehoe and provide evidence for the incipient transformation of basaltic lava from pahoehoe to aa.     

Emplacement conditions of the c. 1,600-year bp Collier Cone lava flow, Oregon: a LiDAR investigation

A long-standing question in lava flow studies has been how to infer emplacement conditions from information preserved in solidified flows. From a hazards perspective, volumetric flux (effusion rate) is the parameter of most interest for open-channel lava flows, as the effusion rate is important for estimating the final flow length, the rate of flow advance, and the eruption duration. The relationship between effusion rate, flow length, and flow advance rate is fairly well constrained for basaltic lava flows, where there are abundant recent examples for calibration. Less is known about flows of intermediate compositions (basaltic andesite to andesite), which are less frequent and where field measurements are limited by the large block sizes and the topographic relief of the flows. Here, we demonstrate ways in which high-resolution digital topography obtained using Light Detection and Ranging (LiDAR) systems can provide access to terrains where field measurements are difficult or impossible to collect. We map blocky lava flow units using LiDAR-generated bare earth digital terrain models (DTMs) of the Collier Cone lava flow in the central Oregon Cascades. We also develop methods using geographic information systems to extract and quantify morphologic features such as channel width, flow width, flow thickness, and slope. Morphometric data are then analyzed to estimate both effusion rates and emplacement times for the lava flow field. Our data indicate that most of the flow outline (which comprises the earliest, and most voluminous, flow unit) can be well explained by an average volumetric flux ～14–18?m³/s; channel data suggest an average flux ～3?m³/s for a later, channel-filling, flow unit. When combined with estimates of flow volume, these data suggest that the Collier Cone lava flow was most likely emplaced over a time scale of several months. This example illustrates ways in which high-resolution DTMs can be used to extract and analyze morphologic measurements and how these measurements can be analyzed to estimate emplacement conditions for inaccessible, heavily vegetated, or extraterrestrial lava flows.     

Downslope flow models of a Bingham liquid: Implications for lava flows

It is widely recognized that lavas behave as Bingham liquids, which are characterized by a yield stress σ and a plastic viscosity η. We consider two models describing downslope flows of a Bingham liquid with different aspect ratios A (= flow height/flow width): model 1 with A 1 and model 2 with A ≈ 1. Sufficiently uphill with respect to the front, such flows can be considered as laminar and locally isothermal. For both models, we obtain analytically the steady-state solution of the Navier-Stokes equations and the constitutive equation for a Bingham liquid. We study the flow height and velocity as functions of flow rate, rheological parameters and ground slope. It is found that such flows remain in the Newtonian regime at low yield stresses (σ 10³dyne/cm²), but the transition to the Bingham regime also depends on flow rate and occurs at higher values of σ for higher flow rates: for instance, a high aspect ratio flow (model 2) is still very close to the Newtonian regime at σ = 10⁴ dyne/cm², if the flow rate is greater than 10⁵ g/s. In the Bingham regime, flow heights are generally greater and flow velocities are smaller than in the Newtonian regime; moreover, flow heights are independent of flow rate, so that a change in flow rate results exclusively in a velocity change. After assuming a specific temperature dependence of σ and η between the solidus and the liquidus temperatures of an ideal Bingham liquid (1000°C and 1200 °C respectively), flow heights and velocities are examined as functions of temperature along the flow. Several effects observed in lava flows are predicted by these models and allow a more quantitative insight into the behaviour of lava flows.     

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     

Thermal–mechanical modeling of cooling history and fracture development in inflationary basalt lava flows

Thermal–mechanical analyses of isotherms in low-volume basalt flows having a range of aspect ratios agree with inferred isotherm patterns deduced from cooling fracture patterns in field examples on the eastern Snake River Plain, Idaho, and highlight the caveats of analytical models of sheet flow cooling when considering low-volume flows. Our field observations show that low-volume lava flows have low aspect ratios (width divided by thickness), typically < 5. Four fracture types typically develop: column-bounding, column-normal, entablature (all of which are cooling fractures), and inflation fractures. Cooling fractures provide a proxy for isotherms during cooling and produce patterns that are strongly influenced by flow aspect ratio. Inflation fractures are induced by lava pressure-driven inflationary events and introduce a thermal perturbation to the flow interior that is clearly evidenced by fracture patterns around them. Inflation fracture growth occurs incrementally due to blunting of the lower tip within viscoelastic basalt, allowing the inflation fracture to pivot open. The final stage of growth involves propagation beyond the blunted tip towards the stress concentration at the tapered tip of a lava core, resulting in penetration of the core that causes quenching of the lava and the formation of a densely fractured entablature. We present numerical models that include the effects of inflation fractures on lava cooling and which support field-based inferences that inflation fractures depress the isotherms in the vicinity of the fracture, cause a subdivision of the lava core, control the location of the final portion of the lava flow to solidify, and cause significant changes in the local cooling fracture orientations. In addition to perturbing isotherms, inflation fractures cause a lava flow to completely solidify in a shorter amount of time than an identically shaped flow that does not contain an inflation fracture.     

Seismicity of an andesitic volcano during block-lava effusion: Volcán de Colima, México, November 1998–January 1999

The block-lava effusion at Volcán de Colima, México began on November 20, 1998, after 12 months of seismic activity, and ended about 80 days later. Three types of seismic events were observed during the lava effusion. Volcano—tectonic earthquakes occurred mainly at the very beginning and after the termination of lava effusion. Explosion earthquakes occurred frequently during the period of the maximum rate in lava effusion. The remainder of the seismic signals were associated with pyroclastic flows and rockfalls from the lava dome. These latter signals increased sharply in number at the onset of lava effusion. The rate of occurrence remained high when the lava discharge rate decreased but gradually decreased after the termination of lava effusion. Maximum daily durations of seismic signals are proportional to the daily volumetric output of lava, indicating the dependence of the number of pyroclastic flows on the rate of lava output. A log-log plot of seismic signal duration vs. number of events with this duration displays a linear relationship. The short-period seismic signals can be divided into three categories based on duration: short events with durations less than 100 s; intermediate events with durations between 100 and 250 s; and long events with durations longer than 250 s. We infer that long events correspond to pyroclastic flows with mean deposit volume 2×10⁵ m³, and intermediate events represent pyroclastic flows with mean deposit volume 1×10³ m³.Editorial responsibility: J McPhie     

Persistent activity and violent strombolian eruptions at Vesuvius between 1631 and 1944

During the period 1631–1944, Vesuvius was in persistent activity with alternating mild strombolian explosions, quiet effusive eruptions, and violent strombolian eruptions. The major difference between the predominant style of activity and the violent strombolian stages is the effusion rate. The lava effusion rate during major eruptions was in the range 20–100 m³/s, higher than during mild activity and quiet effusion (0.1–1 m³/s). The products erupted during the mild activity and major paroxysms have different degree of crystallization. Highly porphyritic lava flows are slowly erupted during years-long period of mild activity. This activity is fed by a magma accumulating at shallow depth within the volcanic edifice. Conversely, during the major paroxysms, a fast lava flow precedes the eruption of a volatile-rich, crystal-poor magma. We show that the more energetic eruptions are fed by episodic, multiple arrival of discrete batches of magma rising faster and not degassing during the ascent. The rapidly ascending magma pushes up the liquid residing in the shallow reservoir and eventually reaches the surface with its full complement of volatiles, producing kilometer-high lava fountains. Rapid drainage of the shallow reservoir occasionally caused small caldera collapses. The major eruptions act to unplug the upper part of the feeding system, erupting the cooling and crystallizing magma. This pattern of activity lasted for 313 y, but with a progressive decrease in the number of more energetic eruptions. As a consequence, a cooling plug blocked the volcano until it eventually prevented the eruption of new magma. The yearly probability of having at least one violent strombolian eruption has decreased from 0.12 to 0.10 from 1944 to 2007, but episodic seismic crises since 1979 may be indicative of new episodic intrusions of magma batches.