Fracture patterns at lava–ice contacts on Kokostick Butte, OR, and Mazama Ridge, Mount Rainier, WA: Implications for flow emplacement and cooling histories

Cooling lava commonly develop polygonal joints that form equant hexagonal columns. Such fractures are formed by thermal contraction resulting in an isotropic tensional stress regime. However, certain linear cooling fracture patterns observed at some lava–ice contacts do not appear to fit the model for formation of cooling fractures and columns because of their preferred orientations. These fracture types include sheet-like (ladder-like rectangular fracture pattern), intermediate (pseudo-aligned individual column-bounding fractures), and pseudopillow (straight to arcuate fractures with perpendicular secondary fractures caused by water infiltration) fractures that form the edges of multiple columns along a single linear fracture. Despite the relatively common occurrence of these types of fractures at lava–ice contacts, their significance and mode of formation have not been fully explored. This study investigates the stress regimes responsible for producing these unique fractures and their significance for interpreting cooling histories at lava–ice contacts.Data was collected at Kokostick Butte dacite flow at South Sister, OR, and Mazama Ridge andesite flow at Mount Rainier, WA. Both of these lava flows have been interpreted as being emplaced into contact with ice and linear fracture types have been observed on their ice-contacted margins. Two different mechanisms are proposed for the formation of linear fracture networks. One possible mechanism for the formation of linear fracture patterns is marginal bulging. Melting of confining ice walls will create voids into which flowing lava can deform resulting in margin-parallel tension causing margin-perpendicular fractures. If viewed from the ice-wall, these fractures would be steeply dipping, linear fractures. Another possible mechanism for the formation of linear fracture types is gravitational settling. Pure shear during compression and settling can result in a tensional environment with similar consequences as marginal inflation. In addition to this, horizontally propagating cooling fractures will be directly influenced by viscous strain caused by the settling of the flow. This would cause preferential opening of fractures horizontally, resulting in vertically oriented fractures.It is important to note that the proposed model for the formation of linear fractures is dependent on contact with and confinement by glacial ice. The influence of flow or movement on cooling fracture patterns has not been extensively discussed in previous modeling of cooling fractures. Rapid cooling of lava by the interaction with water and ice will increase the ability to the capture and preserve perturbations in the stress regime.     

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.     

The effect of crystallization on the rheology and dynamics of lava flows

The dynamics of a lava flow is studied by a two-dimensional model describing a viscous fluid with Bingham rheology, flowing down a slope. The temperature in the flow is calculated assuming that heat is transferred through the plug by conduction and is lost by radiation to the atmosphere at the top of the flow. Taken into account is that the increasing crystallization takes place in the flow as a consequence of cooling. The lava viscosity and yield stress are expressed as a function of crystallization degree as well as of temperature: in particular it is assumed that yield stress reaches a maximum value above the solidus temperature, according to experimental data. Dynamical variables, such as velocity and thickness of the flow, are calculated for different values of the maximum crystallization degree and the flow rate. The model shows how the lava flow dynamics is affected by cooling and crystallization. The cooling of the flow is controlled by the increase of yield stress, which produces a thicker plug and makes the heat loss slower. The increasing crystallization has two opposing effects on viscosity: it produces an increase of viscosity, but at the same time produces an increase of yield stress and hence reduces the heat loss and keeps the internal temperature high. As a consequence, lava flows are significantly affected by the dependence of yield stress on temperature and scarcely by the maximum crystallization degree.     

Rheology of arc dacite lavas: experimental determination at low strain rates

Andesitic–dacitic volcanoes exhibit a large variety of eruption styles, including explosive eruptions, endogenous and exogenous dome growth, and kilometer-long lava flows. The rheology of these lavas can be investigated through field observations of flow and dome morphology, but this approach integrates the properties of lava over a wide range of temperatures. Another approach is through laboratory experiments; however, previous studies have used higher shear stresses and strain rates than are appropriate to lava flows. We measured the apparent viscosity of several lavas from Santiaguito and Bezymianny volcanoes by uniaxial compression, between 1,109 and 1,315?K, at low shear stress (0.085 to 0.42?MPa), low strain rate (between 1.1?×?10^?8 and 1.9?×?10^?5?s^?1), and up to 43.7 % total deformation. The results show a strong variability of the apparent viscosity between different samples, which can be ascribed to differences in initial porosity and crystallinity. Deformation occurs primarily by compaction, with some cracking and/or vesicle coalescence. Our experiments yield apparent viscosities more than 1 order of magnitude lower than predicted by models based on experiments at higher strain rates. At lava flow conditions, no evidence of a yield strength is observed, and the apparent viscosity is best approached by a strain rate- and temperature-dependent power law equation. The best fit for Santiaguito lava, for temperatures between 1,164 and 1,226?K and strain rates lower than 1.8?×?10^?4?s^?1, is $ \log {\eta_{\text{app}}} = - 0.738 + 9.24 \times {10^3}{/}T(K) - 0.654 \cdot \log \dot{\varepsilon } $ where η _app is apparent viscosity and $ \dot{\varepsilon } $ is strain rate. This equation also reproduced 45 data for a sample from Bezymianny with a root mean square deviation of 0.19 log unit Pa?s. Applying the rheological model to lava flow conditions at Santiaguito yields calculated apparent viscosities that are in reasonable agreement with field observations and suggests that internal shear heating may be significant ongoing heat source within these flows, enabling highly viscous lava to travel long distances.     

Causes and consequences of pressurisation in lava dome eruptions

High total and fluid pressures develop in the interior of high-viscosity lava domes and in the uppermost parts of the feeding conduit system as a consequence of degassing. Two effects are recognised and are modelled quantitatively. First, large increases in magma viscosity result from degassing during magma ascent. Strong vertical gradients in viscosity result and large excess pressures and pressure gradients develop at the top of the conduit and in the dome. Calculations of conduit flow show that almost all the excess pressure drop from the chamber in an andesitic dome eruption occurs during the last several hundred metres of ascent. Second, microlites grow in the melt phase as a consequence of undercooling caused by gas loss. Rapid microlite growth can cause large excess fluid pressures to develop at shallow levels. Theoretically closed-system microlite crystallization can increase local pressure by a few tens of MPa, although build up of pressure will be countered by gas loss through permeable flow and expansion by viscous flow. Microlite crystallization is most effective in causing excess gas pressures at depths of a few hundred metres in the uppermost parts of the conduit and dome interior. Some of the major phenomena of lava dome eruptions can be attributed to these pressurisation effects, including spurts of growth, cycles of dome growth and subsidence, sudden onset of violent explosive activity and disintegration of lava during formation of pyroclastic flows. The characteristic shallow-level, long-period and hybrid seismicity, characteristic of dome eruptions, is attributed to the excess fluid pressures, which are maintained close to the fracture strength of the dome and wallrock, resulting in fluid movement during formation of tensile and shear fractures within the dome and upper conduit.     

Fault textures in volcanic conduits: evidence for seismic trigger mechanisms during silicic eruptions

It is proposed that fault textures in two dissected rhyolitic conduits in Iceland preserve evidence for shallow seismogenic faulting within rising magma during the emplacement of highly viscous lava flows. Detailed field and petrographic analysis of such textures may shed light on the origin of long-period and hybrid volcanic earthquakes at active volcanoes. There is evidence at each conduit investigated for multiple seismogenic cycles, each of which involved four distinct evolutionary phases. In phase 1, shear fracture of unrelaxed magma was triggered by shear stress accumulation during viscous flow, forming the angular fracture networks that initiated faulting cycles. Transient pressure gradients were generated as the fractures opened, which led to fluidisation and clastic deposition of fine-grained particles that were derived from the fracture walls by abrasion. Fracture networks then progressively coalesced and rotated during subsequent slip (phase 2), developing into cataclasite zones with evidence for multiple localised slip events, fluidisation and grain size reduction. Phase 2 textures closely resemble those formed on seismogenic tectonic faults characterised by friction-controlled stick-slip behaviour. Increasing cohesion of cataclasites then led to aseismic, distributed ductile deformation (phase 3) and generated deformed cataclasite zones, which are enriched in metallic oxide microlites and resemble glassy pseudotachylite. Continued annealing and deformation eventually erased all structures in the cataclasite and formed microlite-rich flow bands in obsidian (phase 4). Overall, the mixed brittle–ductile textures formed in the magma appear similar to those formed in lower crustal rocks close to the brittle–ductile transition, with the rheological response mediated by strain-rate variations and frictional heating. Fault processes in highly viscous magma are compared with those elsewhere in the crust, and this comparison is used to appraise existing models of volcano seismic activity. Based on the textures observed, it is suggested that patterns of long-period and hybrid earthquakes at silicic lava domes reflect friction-controlled stick-slip movement and eventual healing of fault zones in magma, which are an accelerated and smaller-scale analogue of tectonic faults.Editorial responsibility: J. Stix     

Flow and fracture maps for basaltic rock deformation at high temperatures

Understanding how the strength of basaltic rock varies with the extrinsic conditions of stress state, pressure and temperature, and the intrinsic rock physical properties is fundamental to understanding the dynamics of volcanic systems. In particular it is essential to understand how rock strength at high temperatures is limited by fracture. We have collated and analysed laboratory data for basaltic rocks from over 500 rock deformation experiments and plotted these on principal stress failure maps. We have fitted an empirical flow law (Norton’s law) and a theoretical fracture criterion to these data. The principal stress failure map is a graphical representation of ductile and brittle experimental data together with flow and fracture envelopes under varying strain rate, temperature and pressure. We have used these maps to re-interpret the ductile–brittle transition in basaltic rocks at high temperatures and show, conceptually, how these failure maps can be applied to volcanic systems, using lava flows as an example.     

Collapse mechanism of the Paleogene Sakurae cauldron, SW Japan

The Paleogene Sakurae cauldron of SW Japan is characterized by a nested structure with a polygonal outline (21×13 km²) including a circular collapsed part (5 km in diameter). Total thickness of the caldera infill amounts to 2,000 m. The lower member of the infill consists mainly of felsic crystal tuff and lesser intercalated andesitic lava flows, whereas the upper member is composed of high-grade ignimbrite capped with a large rhyolitic lava dome. These members represent the first and second stage eruptions, respectively. Faults bounding the cauldron rim comprise intersecting radial and concentric faults, producing the polygonal outline of this cauldron. The primary collapse of this cauldron thus occurred as a polygonal caldera basin where products of the first stage eruption accumulated. In contrast, the inner collapse part is defined by a ring fracture system. This sector subsided concurrently with accumulation of the high-grade ignimbrite of the second stage eruption. This inner circular collapse thus represents syn-eruptional subsidence concurrent with the climactic eruption. Magma drainage during the first stage probably induced outward-dipping ring fractures in the chamber roof. Opening of the ring fractures following subsidence of the central bell-jar block caused rapid evacuation of magma as voluminous pumice flows, even though magma pressure may have decreased to some degree.     

Formation of clastogenic lava flows during fissure eruption and scoria cone collapse: the 1986 eruption of Izu-Oshima Volcano, eastern Japan

The 1986 eruption of B fissure at Izu-Oshima Volcano, Japan, produced, among other products, one andesite and two basaltic andesite lava flows. Locally the three flows resemble vent-effused holocrystalline blocky or aa lava; however, remnant clast outlines can be identified at most localities, indicating that the flows were spatter fed or clastogenic. The basaltic andesite flows are interpreted to have formed by two main processes: (a) reconstitution of fountain-generated spatter around vent areas by syn-depositional agglutination and coalescence, followed by extensional non-particulate flow, and (b) syn-eruptive collapse of a rapidly built spatter and scoria cone by rotational slip and extensional sliding. These processes produced two morphologically distinct lobes in both flows by: (a) earlier non-particulate flow of agglutinate and coalesced spatter, which formed a thin lobe of rubbly aa lava (ca. 5 m thick) with characteristic open extension cracks revealing a homogeneous, holocrystalline interior, and (b) later scoria-cone collapse, which created a larger lobe of irregular thickness (<20 m) made of large detached blocks of scoria cone interpreted to have been rafted along on a flow of coalesced spatter. The source regions of these lava flows are characterized by horseshoe-shaped scarps (<30 m high), with meso-blocks (ca. 30 m in diameter) of bedded scoria at the base. One lava flow has a secondary lateral collapse zone with lower (ca. 7 m) scarps. Backward-tilted meso-blocks are interpreted to be the product of rotational slip, and forward-tilted blocks the result of simple toppling. Squeeze-ups of coalesced spatter along the leading edge of the meso-blocks indicate that coalescence occurred in the basal part of the scoria cone. This low-viscosity, coalesced spatter acted as a lubricating layer along which basal failure of the scoria cone occurred. Rotational sliding gave way to extensional translational sliding as the slide mass spread out onto the present caldera floor. Squeeze-ups concentrated at the distal margin indicate that the extensional regime changed to one of compression, probably as a result of cooling of the flow front. Sliding material piled up behind the slowing flow front, and coalesced spatter was squeezed up from the interior of the flow through fractures and between rafted blocks. The andesite flow, although morphologically similar to the other two flows, has a slightly different chemical composition which corresponds to the earliest stage of the eruption. It is a much smaller lava flow emitted from the base of the scoria cone 2 days after the eruption had ceased. This lava is interpreted to have been formed by post-depositional coalescence of spatter under the influence of the in-situ cooling rate and load pressure of the deposit. Extrusion occurred through the lower part of the scoria cone, and subsequent non-particulate flow of coalesced material produced a blocky and aa lava flow. The mechanisms of formation of the lava flows described may be more common during explosive eruptions of mafic magma than previously envisaged. Received: 30 May 1997 / Accepted: 19 May 1998     

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.     

Mt. Niragongo: renewed activity of the lava lake

On the 21st of June, 1982, Mt. Niragongo ended a period of dormancy that had begun on January 11, 1977, and fresh lava began to flow into the 800-m-deep crater. On October 3, a huge lava lake, wider and deeper than any previously observed (500 m across and close to 400 m deep) rose to within 440 m of the crater rim. The observed activity consisted of a large, central upwelling fountain of very fluid lava from which concentric lava waves expanded radially; numerous small, relatively viscous lava flows creeped over the surrounding thin solidified crust, that covered about 95% of the lake area. These observed features seem to characterize the upper part of a large convective system. The persistence of such an extraordinarily large steady-state lava lake may be due to the equally exceptional fluidity of the magma rising at the intersection of four different tectonic trends of fractures in the subvolcanic basement.     

Experiments on internal strain in lava dome cross sections

Simple experiments have been conducted to study the strain evolution in lava dome cross sections. A viscous fluid is injected vertically from a reservoir into a feeding conduit. Silicone putty is used as analogue magma. Two-dimensional experiments allow the assessment of the internal strain within the dome. Particle paths are symmetrical on either side of a central line passing through the feeding conduit and display parabolic trajectories. The highest strain zone is located above the extrusion zone. In cross sections, stretch trajectories show a remarkable concentric pattern, wrapping around the extrusion zone of the analogue magma. To the lateral margins, a triple junction of stretch trajectories defines an isotropic point in the strain field. In the main central part of the dome, an intermediate zone of reversed sense of shearing is caused by a change in the sign of the velocity gradient with respect to that in the upper and lower zones. Knowledge of this evolving strain pattern can provide a better understanding of the evolution of natural domes. Also, it can help to unravel the kinematic history of ancient domes partly removed by erosion.     

The initial cooling of pahoehoe flow lobes

 In this paper we describe a new thermal model for the initial cooling of pahoehoe lava flows. The accurate modeling of this initial cooling is important for understanding the formation of the distinctive surface textures on pahoehoe lava flows as well as being the first step in modeling such key pahoehoe emplacement processes as lava flow inflation and lava tube formation. This model is constructed from the physical phenomena observed to control the initial cooling of pahoehoe flows and is not an empirical fit to field data. We find that the only significant processes are (a) heat loss by thermal radiation, (b) heat loss by atmospheric convection, (c) heat transport within the flow by conduction with temperature and porosity-dependent thermal properties, and (d) the release of latent heat during crystallization. The numerical model is better able to reproduce field measurements made in Hawai'i between 1989 and 1993 than other published thermal models. By adjusting one parameter at a time, the effect of each of the input parameters on the cooling rate was determined. We show that: (a) the surfaces of porous flows cool more quickly than the surfaces of dense flows, (b) the surface cooling is very sensitive to the efficiency of atmospheric convective cooling, and (c) changes in the glass forming tendency of the lava may have observable petrographic and thermal signatures. These model results provide a quantitative explanation for the recently observed relationship between the surface cooling rate of pahoehoe lobes and the porosity of those lobes (Jones 1992, 1993). The predicted sensitivity of cooling to atmospheric convection suggests a simple field experiment for verification, and the model provides a tool to begin studies of the dynamic crystallization of real lavas. Future versions of the model can also be made applicable to extraterrestrial, submarine, silicic, and pyroclastic flows. Received: 26 November 1994 / Accepted: 1 December 1995     

Propagation pathways and fluid transport of hydrofractures in jointed and layered rocks in geothermal fields

Palaeofluid-transporting systems, observed as networks of mineral-filled veins in deeply eroded parts of extinct geothermal fields, indicate that hydrofractures commonly supply fluids to geothermal fields. Here we examine well-exposed vein networks that occur at crustal depths of around 1.5 km below the initial surface of the Tertiary lava pile in North Iceland. The veins are located in the damage zone of a major fault zone that dissects basaltic lava flows, the most common host rocks of geothermal fields in Iceland. The lava flows contain numerous weaknesses, particularly columnar (cooling) joints and contacts. For hydrofractures to supply fluids to geothermal fields, the fractures must be able to propagate, and transport fluids, to the surface. We explore hydrofracture pathway formation using boundary-element models of hydrofractures with fluid overpressure varying linearly from 10 MPa at the fracture centre to 0 MPa at the fracture tip (or the fluid front). The hydrofractures propagate through a vertically jointed and horizontally layered pile of lava flows with a general rock-matrix Young’s modulus of 1×10¹⁰ Pa and a Poisson’s ratio of 0.25. The joints and contacts between layers are modelled as internal springs, each with a stiffness (‘strength’) of 6 MPa/m. The location and sizes of discontinuities, as well as the location of the hydrofracture tip, vary between the models. The results indicate that tensile stresses generated at the tip of an overpressured hydrofracture can open up horizontal and vertical discontinuities out to a considerable distance from the tip, and that these discontinuities eventually link up to form the hydrofracture pathway. Analytical models indicate that for a hot spring of a given yield associated with a fault, the dimensions of the fluid-transporting part of the fault are likely to be similar for a typical normal fault and a strike-slip fault. Also, a hot spring of yield 180 l/s (the maximum in the low-temperature fields of Iceland) can be supplied through a hydrofracture of aperture 3 mm and trace length 1.2 m. These dimensions are very similar to those of typical veins in the studied networks. Buoyancy, rather than excess pressure in the fluid source, appears to be the primary driving force of hydrofractures in the geothermal fields of Iceland.     

A laboratory investigation into the effects of slope on lava flow morphology

In an attempt to model the effect of slope on the dynamics of lava flow emplacement, four distinct morphologies were repeatedly produced in a series of laboratory simulations where polyethylene glycol (PEG) was extruded at a constant rate beneath cold sucrose solution onto a uniform slope which could be varied from 1° through 60°. The lowest extrusion rates and slopes, and highest cooling rates, produced flows that rapidly crusted over and advanced through bulbous toes, or pillows (similar to subaerial “toey” pahoehoe flows and to submarine pillowed flows). As extrusion rate and slope increased, and cooling rate decreased, pillowed flows gave way to rifted flows (linear zones of liquid wax separated by plates of solid crust, similar to what is observed on the surface of convecting lava lakes), then to folded flows with surface crusts buckled transversely to the flow direction, and, at the highest extrusion rates and slopes, and lowest cooling rates, to leveed flows, which solidified only at their margins. A dimensionless parameter, Ψ, primarily controlled by effusion rate, cooling rate and flow viscosity, quantifies these flow types. Increasing the underlying slope up to 30° allows the liquid wax to advance further before solidifying, with an effect similar to that of increasing the effusion rate. For example, conditions that produce rifted flows on a 10° slope result in folded flows on a 30° slope. For underlying slopes of 40°, however, this trend reverses, slightly owing to increased gravitational forces relative to the strength of the solid wax. Because of its significant influence on heat advection and the disruption of a solid crust, slope must be incorporated into any quantitative attempt to correlate eruption parameters and lava flow morphologies. These experiments and subsequent scaling incorporate key physical parameters of both an extrusion and its environment, allowing their results to be used to interpret lava flow morphologies on land, on the sea floor, and on other planets.     

Surficial fractures in the Navajo sandstone,south-western USA: the roles of thermal cycles,rainstorms, granular disintegration,and iterative cracking

Deep (> 5 m) sheeting fractures in the Navajo sandstone are evident at numerous sites in southern Utah and derive from tectonic stresses. Strong diurnal thermal cycles are, however, the likely triggers for shallow (< 0.3 m) sheeting fractures. Data from subsurface thermal sensors reveal that large temperature differences between sensors at 2 and 15 cm depth on clear summer afternoons are as great as those that trigger sheeting fractures in exposed California granite. Extensive polygonal patterns in the Navajo sandstone are composed of surface-perpendicular fractures and were produced by contractile stresses. Numerous studies have shown that porewater diminishes the tensile strength of sandstone. Based on our thermal records, we propose that cooling during monsoonal rainstorms triggers polygonal fracturing of temporarily weakened rock. On steep outcrops, polygonal patterns are rectilinear and orthogonal, with T-vertices. Lower-angle slopes host hexagonal patterns (defined by the dominance of Y-vertices). Intermediate patterns with rectangles and hexagons of similar scale are common. We posit that outcropping fractures are advancing downward by iterative steps, and that hexagons on sandstone surfaces (like prismatic columns of basalt) have evolved from ancestral orthogonal polygons of similar scale. In lava flows, fractures elongate intermittently as they follow a steep thermal gradient (the source of stress) as it rapidly moves through the rock mass. In our model, a steep, surficial thermal gradient descends through unfractured sandstone, but at the slow pace of granular disintegration. Through time, as the friable rock on stable slopes erodes, iterative cracking advances into new space. Hexagonal patterns form as new fractures, imperfectly guided by the older ones, propagate in new directions, and vertices drift into a configuration that minimizes the ratio of fracture length to polygon area. © 2020 John Wiley & Sons, Ltd.     

Dynamics of lava flow fronts, Arenal Volcano, Costa Rica

Arenal Volcano has effused basaltic andesite lava flows nearly continuously since September, 1968. The two different kinds of material in flows, lava and lava debris, have different rheologic properties and dynamic behavior. Flow morphology depends on the relationship between the amount and distribution of the lava and the debris, and to a lesser extent the ground morphology.Two main units characterize the flows: the channel zone and the frontal zone. The channel zone consists of two different units, the levées and the channel proper. A velocity profile in the channel shows a maximum value at the plug where the rate of shear is zero, and a velocity gradient increasing outward until, at the levées, the velocity becomes zero. Cooling produces a marked temperature gradient in the flow, leading to the formation of debris by brittle fracture when a critical value of shear rate to viscosity is reached. When the lava supply ceases, much of this debris and part of the lava is left behind after the flow nucleus drains out, forming a collapsed channel.Processes at the frontal zone include levée formation, debris formation, the change in shape of the front, and the choice of the flow path. These processes are controlled primarily by the rheological properties of the lava.Frontal zone dynamics can be understood by fixing the flow front as the point of reference. The lava flows through the channel into the front where it flows out into the levées, thereby increasing the length of the channel and permitting the front to advance. The front shows a relationship of critical height to the yield strength (τ₀) surface tension, and slope; its continued movement is activated by the pressure of the advancing lava in the channel behind. For an ideal flow (isothermal, homogeneous, and isotropic) the ratio of the section of channel proper to the section of levées is calculated and the distance the front will have moved at any time t_x can be determined once the amount of lava available to the front is known. Assuming that the velocity function of the front {G(t)} during the collapsing stage is proportional to the entrance pressure of the lava at the channel-front boundary, an exponential decrease of velocity through time is predicted, which shows good agreement with actual frontal velocity measurements taken on two flows. Local variations in slope have a secondary effect on frontal velocities.Under conditions of constant volume the frontal zone can be considered as a machine that consumes energy brought in by the lava to perform work (front advancement). While the front will use its potential energy to run the process, the velocity at which it occurs is controlled by the activation energy that enters the system as the kinetic energy of the lava flowing into the front. A relation for the energy contribution due to frontal acceleration is also derived. Finally the entrance pressure, that permits the front to deform, is calculated. Its small value confirms that the lava behaves very much like a Bingham plastic.     

Salt-bearing fumarole deposits in the summit crater of Oldoinyo Lengai, Northern Tanzania: interactions between natrocarbonatite lava and meteoric water

Oldoinyo Lengai in the Northern Tanzania rift is the only active nephelinite–carbonatite stratovolcano. We report the discovery of thermonatrite, aphthitalite, halite and sylvite fumarole deposits on recent natrocarbonatite lava flows erupted in the summit crater during the wet season. These salt deposits occur as delicate, concave fringes or tubes that line the cooling cracks in the lava flows and consist of intergrowths of euhedral crystals. The presence of a dark altered zone, depleted in halides and alkalies, adjacent to cooling cracks and observations of steam fumaroles emanating from the fractures suggest that the salts are formed by sublimation from saturated vapours generated by the extrusion of lavas over meteoric water. The crystallisation sequence recorded in the salts suggests that mixing between meteoric steam and magmatic CO₂ and H₂S occurs at high temperatures resulting in the sublimation of carbonates and sulphates. At lower temperatures the vapours are dominated by meteoric steam and sublimate halides. The high solubility of the fumarole salts within meteoric water and their formation only during the wet season implies that these are ephemeral deposits that are unlikely to be preserved in the geological record.     

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.     

Criteria for recognition of constructional silicic lava flow surfaces

Surface-exposure dating (SED) methods typically rely on the measurement of a geochemical parameter that systematically changes with time. A pivotal task in the calibration of many of these techniques is to demonstrate that lava flow surfaces sampled for dating have not experienced erosion. Although criteria for identification of constructional basaltic lava flow surfaces have been published, no such criteria presently exist for the recognition of constructional silicic flows. Here we present several criteria for identifying constructional silicic lava flow features in the field. First, crease structures are fractures with easily identified, curved, striated walls that are commonly observed on recent and active silicic lava flows. Crease structures form during extrusion, and are resistant to mechanical disintegration because they expose dense material from the flow interior. Second, some crease structures break apart during formation, leaving a deposit of striated blocks on the flow surface. Crease structure blocks are striated on only one side, whereas blocks from internal columnar joints exposed through erosion are striated on two or more sides. Only the striated side of the crease structure block is definitively constructional. Finally, many silicic flow surfaces exhibit expanded or breadcrusted textures. These features consist of a dense, fractured rind, 1 –2cm thick, enclosing highly vesicular material. Breadcrust flow textures appear similar to breadcrust bombs produced during volcanic explosions, so it is imperative to demonstrate that they are part of the lava flow surface. These criteria should enable investigators to positively identify constructional silicic lava flow surfaces when calibrating an SED method.