Effusive eruption of viscous silicic magma triggered and driven by recharge: a case study of the Cerro Chascon-Runtu Jarita Dome Complex in Southwest Bolivia

 The Cerro Chascon-Runtu Jarita Complex is a group of ten Late Pleistocene (∼85 ka) lava domes located in the Andean Central Volcanic Zone of Bolivia. These domes display considerable macroscopic and microscopic evidence of magma mixing. Two groups of domes are defined chemically and geographically. A northern group, the Chascon, consists of four lava bodies of dominantly rhyodacite composition. These bodies contain 43–48% phenocrysts of plagioclase, quartz, sanidine, biotite, and amphibole in a microlite-poor, rhyolitic glass. Rare mafic enclaves and selvages are present. Mineral equilibria yield temperatures from 640 to 750  °C and log ƒO₂ of –16. Geochemical data indicate that the pre-eruption magma chamber was zoned from a dominant volume of 68% to minor amounts of 76% SiO₂. This zonation is best explained by fractional crystallization and some mixing between rhyodacite and more evolved compositions. The mafic enclaves represent magma that intruded but did not chemically interact much with the evolved magmas. A southern group, the Runtu Jarita, is a linear chain of six small domes (<1 km³ total volume) that probably is the surface expression of a dike. The five most northerly domes are composites of dacitic and rhyolitic compositions. The southernmost dome is dominantly rhyolite with rare mafic enclaves. The composite domes have lower flanks of porphyritic dacite with ∼35 vol.% phenocrysts of plagioclase, orthopyroxene, and hornblende in a microlite-rich, rhyodacitic glass. Sieve-textured plagioclase, mixed populations of disequilibrium plagioclase compositions, xenocrystic quartz, and sanidine with ternary composition reaction rims indicate that the dacite is a hybrid. The central cores of the composite domes are rhyolitic and contain up to 48 vol.% phenocrysts of plagioclase, quartz, sanidine, biotite, and amphibole. This is separated from the dacitic flanks by a banded zone of mingled lava. Macroscopic, microscopic, and petrologic evidence suggest scavenging of phenocrysts from the silicic lava. Mineral equilibria yield temperatures of 625–727  °C and log ƒO₂ of –16 for the rhyolite and 926–1000  °C and log ƒO₂ of –9.5 for the dacite. The rhyolite is zoned from 73 to 76% SiO₂, and fractionation within the rhyolite composition produced this variation. Most of the 63–73% SiO₂ compositional range of the lava in this group is the result of mixing between the hybrid dacite and the rhyolite. Eruption of both groups of lavas apparently was triggered by mafic recharge. A paucity of explosive activity suggests that volatile and thermal exchanges between reservoir and recharge magmas were less important than volume increase and the lubricating effects of recharge by mafic magmas. For the Runtu Jarita group, the eruption is best explained by intrusion of a dike of dacite into a chamber of crystal-rich rhyolite close to its solidus. The rhyolite was encapsulated and transported to the surface by the less-viscous dacite magma, which also acted as a lubricant. Simultaneous effusion of the lavas produced the composite domes, and their zonation reflects the subsurface zonation. The role of recharge by hotter, more fluid mafic magma appears to be critical to the eruption of some highly viscous silicic magmas. Received: 23 August 1998 / Accepted: 10 March 1999     

Subaqueous, basaltic lava dome and carapace breccia on King George Island, South Shetland Islands, Antarctica

 On King George Island during latest Oligocene/earliest Miocene time, submarine eruptions resulted in the emplacement of a small (ca. 500 m estimated original diameter) basalt lava dome at Low Head. The dome contains a central mass of columnar rock enveloped by fractured basalt and basalt breccia. The breccia is crystalline and is a joint-block deposit (lithic orthobreccia) interpreted as an unusually thick dome carapace breccia cogenetic with the columnar rock. It was formed in situ by a combination of intense dilation, fracturing and shattering caused by natural hydrofracturing during initial dome effusion and subsequent endogenous emplacement of further basalt melt, now preserved as the columnar rock. Muddy matrix with dispersed hyaloclastite and microfossils fills fractures and diffuse patches in part of the fractured basalt and breccia lithofacies. The sparse glass-rich clasts formed by cooling-contraction granulation during interaction between chilled basalt crust and surrounding water. Together with muddy sediment, they were injected into the dome by hydrofracturing, local steam fluidisation and likely explosive bulk interaction. The basalt lava was highly crystallised and degassed prior to extrusion. Together with a low effusion temperature and rapid convective heat loss in a submarine setting, these properties significantly affected the magma rheology (increased the viscosity and shear strength) and influenced the final dome-like form of the extrusion. Conversely, high heat retention was favoured by the degassed state of the magma (minimal undercooling), a thick breccia carapace and viscous shear heating, which helped to sustain magmatic (eruption) temperatures and enhanced the mobility of the flow. Received: 1 August 1996 / Accepted: 15 September 1997     

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.     

Results of repeated leveling surveys at Newberry Volcano, Oregon, and near Lassen Peak Volcano, California

 Personnel from the U.S. Geological Survey's Cascades Volcano Observatory conducted first-order, class-II leveling surveys near Lassen Peak, California, in 1991 and at Newberry Volcano, Oregon, in 1985, 1986, and 1994. Near Lassen Peak no significant vertical displacements had occurred along either of two traverses, 33 and 44 km long, since second-order surveys in 1932 and 1934. At Newberry, however, the 1994 survey suggests that the volcano's summit area had risen as much as 97±22 mm with respect to a third-order survey in 1931. The 1931 and 1994 surveys measured a 37-km-long, east–west traverse across the entire volcano. The 1985 and 1986 surveys, on the other hand, measured only a 9-km-long traverse across the summit caldera with only one benchmark in common with the 1931 survey. Comparison of the 1985, 1986, and 1994 surveys revealed no significant differential displacements inside the caldera. A possible mechanism for uplift during 1931–1994 is injection of approximately 0.06 km³ of magma at a depth of approximately 10 km beneath the volcano's summit. The average magma supply rate of approximately 1×10^–3 km³/year would be generally consistent with the volcano's growth rate averaged over its 600,000-year history (0.7–1.7×10^–3 km³/year). Received: 10 September 1998 / Accepted: 4 December 1998     

Morphological and textural diversity of the Steens Basalt lava flows,Southeastern Oregon,USA: implications for emplacement style and nature of eruptive episodes

This study focuses on Middle Miocene tholeiitic flood basalt lava flows from the Oregon Plateau, northwestern USA (Steens Basalt), and is the first to comprehensively document and evaluate their morphology. Field observations of flows from several sections within and proximal to the main exposures at Steens Mountain have been supplemented with textural and geochemical data, and are used to offer preliminary insights into their emplacement. Compound pahoehoe flows of variable thickness appear to be common throughout the study area, laterally and vertically. These tend to be plagioclase phyric and the morphology and disposition of constituent flow lobes are quite similar to those from other provinces such as Hawaii and the Snake River Plain. Classic a’a flows with brecciated upper and basal crusts are not abundant, but by no means uncommon. Flows with characters different from typical pahoehoe and a’a are also common. Such flows display a range in morphology; flows with preserved upper crusts but brecciated basal crusts, as well as those displaying well-developed flow-top breccias and preserved basal crusts (rubbly pahoehoe) are observed. The Steens Basalt appears to display greater morphological and textural diversity at the outcrop scale than that described for some other flood basalt provinces. The abundant compound pahoehoe flows (often rich in plagioclase phenocrysts) were likely emplaced during slow but sustained eruptive episodes; their constituent lobes show clear evidence for endogenous growth. The relatively aphyric flows with brecciated surfaces (including a’a) hint at higher strain rates and/or higher viscosity, probably caused by higher effusion rates. A couple of sections are characterized by compositionally similar, but morphologically different flows that were possibly part of the same eruption. While differences in pre-eruptive topography could explain this, it is also possible that certain physical parameters changed substantially and abruptly during eruption and that such changes were accompanied by differentiation processes within the plumbing system. It is possible that such observations indicate temporal fluctuations within complex magmatic and eruptive systems, and deserve closer scrutiny.     

New formulae for estimating lava flow volumes at Mt. Etna Volcano, Sicily

 A new data set of Etna lava flows erupted since 1868 has been compiled from eight topographic maps of the volcano published at intervals since then. Volumes of 59 flows or groups of flows were measured from topographic difference maps. Most of these volumes are likely to be considerably more accurate than those published previously. We cut the number of flow volumes down to 25 by selecting those examples for which the volume of an individual eruption could be derived with the highest accuracy. This refined data set was searched for high correlations between flow volume and more directly measurable parameters. Only two parameters showed a correlation coefficient of 70% or greater: planimetric flow area A (70%) and duration of the eruption D (79%). If only short duration (<18 days) flows were used, flow length cubed, L³, had a correlation coefficient of 98%. Using combinations of measured parameters, much more significant correlations with volume were found. Dh had a correlation coefficient of 90% (h is the hydrostatic head of magma above the vent), and  , 92% (where W is mean width and E is the degree of topographic enclosure), and a combination of the two , 97%. These latter formulae were used to derive volumes of all eruptions back to 1868 to compare with those from the complete data set. Values determined from the formulae were, on average, lower by 16% (Dh), 7% (, and 19% . Received: 30 November 1998 / Accepted: 20 June 1999     

Pahoehoe toe dimensions, morphology, and branching relationships at Mauna Ulu, Kilauea Volcano, Hawai'i

 Pahoehoe toe dimensions, morphology, and branching relationships were analyzed in flows emplaced during 1972 at Mauna Ulu, a satellitic shield on the east rift zone of Kilauea Volcano, Hawai'i. In order to characterize regions within flow fields dominated by networks of pahoehoe toes, measurements of toe length, width, thickness, and orientation were completed for 445 toes at 13 sites. Variations in site characteristics, including slope, substrate, and position in the flow field allow an evaluation of the effects of such parameters on toe dimensions. Toe surface morphology (ropy or smooth), local flow lobe position (interior or margin), and connective relationships between toes were documented in the form of detailed toe maps. These maps show the number of branches connecting a given toe to other toes in its local pahoehoe network and illustrate branching patterns. Statistical analyses of toe dimensions and comparisons of pahoehoe toe study sites and sub-populations combined with field observations, evaluation of toe maps, and qualitative examination of toe dimension size distributions show the following: (a) Although there are significant variations at a given site, toes typically have mean lengths (101 cm) greater than mean widths (74 cm) and mean widths greater than mean thicknesses (19 cm); sites that have mean widths greater than mean lengths are those with lower slopes. (b) Where significant site-to-site variations in mean values of a given toe dimension were apparent, these differences could not be directly related to site characteristics. (c) Ropy toes have significantly larger mean values of length, width, and number of branches than smooth toes, and toes with three or more branches have greater lengths, widths, and thicknesses than toes with two or fewer branches, suggesting concentration of flow in these toe types. (d) The skewness of all size distributions of toe length and width to larger values suggests that toes are transitional to larger sheets and channels, consistent with field observations; and (e) Two distinct types of branching patterns (called monolayer and centrally ridged) were observed in preserved pahoehoe flow lobes. The significant variability in measured toe dimensions at Mauna Ulu suggests that toe dimensions are influenced by numerous locally defined, random factors, and that an approach based on stochastic methods can be used to model pahoehoe flow emplacement. Received: 19 January 1998 / Accepted: 24 March 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     

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     

The dynamics of a channel-fed lava flow on Pico Partido volcano, Lanzarote

A short length of channel on Pico Partido volcano, Lanzarote, provides us the opportunity to examine the dynamics of lava flowing in a channel that extends over a sudden break in slope. The 1–2-m-wide, 0.5–2-m-deep channel was built during the 1730–1736 eruptions on Lanzarote and exhibits a sinuous, well-formed channel over a steep (11° slope) 100-m-long proximal section. Over-flow units comprising smooth pahoehoe sheet flow, as well as evidence on the inner channel walls for multiple (at least 11) flow levels, attest to unsteady flow in the channel. In addition, superelevation is apparent at each of the six bends along the proximal channel section. Superelevation results from banking of the lava as it moves around the bend thus causing preferential construction of the outer bank. As a result, the channel profile at each bend is asymmetric with an outer bank that is higher than the inner bank. Analysis of superelevation indicates flow velocities of ~8 m s^–1. Our analysis of the superelevation features is based on an inertia-gravity balance, which we show is appropriate, even though the down-channel flow is in laminar flow. We use a viscosity-gravity balance model, together with the velocities calculated from superelevation, to obtain viscosities in the range 25–60 Pa s (assuming that the lava behaved as a Newtonian liquid). Estimated volume fluxes are in the range 7–12 m³ s^–1. An apparent down-flow increase in derived volume flux may have resulted from variable supply or bulking up of the flow due to vesiculation. Where the channel moves over a sharp break in slope and onto slopes of ~6°, the channel becomes less well defined and widens considerably. At the break of slope, an elongate ridge extends across the channel. We speculate that this ridge was formed as a result of a reduction in velocity immediately below the break of slope to allow deposition of entrained material or accretion of lava to the channel bed as a result of a change in flow regime or depth.     

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.     

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.     

Kīlauea summit overflows: their ages and distribution in the Puna District, Hawai'i

 The tube-fed pāhoehoe lava flows covering much of the northeast flank of Kīlauea Volcano are named the 'Ailā'au flows. Their eruption age, based on published and six new radiocarbon dates, is approximately AD 1445. The flows have distinctive paleomagnetic directions with steep inclinations (40°–50°) and easterly declinations (0°–10°E). The lava was transported ∼40 km from the vent to the coast in long, large-diameter lava tubes; the longest tube (Kazumura Cave) reaches from near the summit to within several kilometers of the coast near Kaloli Point. The estimated volume of the 'Ailā'au flow field is 5.2±0.8 km³, and the eruption that formed it probably lasted for approximately 50 years. Summit overflows from Kīlauea may have been nearly continuous between approximately AD 1290 and 1470, during which time a series of shields formed at and around the summit. The 'Ailā'au shield was either the youngest or the next to youngest in this series of shields. Site-mean paleomagnetic directions for lava flows underlying the 'Ailā'au flows form only six groups. These older pāhoehoe flows range in age from 2750 to <18,000 BP, and the region was inundated by lava flows only three times in the past 5000 years. The known intervals between eruptive events average ∼1600 years and range from ∼1250 years to >2200 years. Lava flows from most of these summit eruptions also reached the coast, but none appears as extensive as the 'Ailā'au flow field. The chemistry of the melts erupted during each of these summit overflow events is remarkably similar, averaging approximately 6.3 wt.% MgO near the coast and 6.8 wt.% MgO near the summit. The present-day caldera probably formed more recently than the eruption that formed the 'Ailā'au flows (estimated termination ca. AD 1470). The earliest explosive eruptions that formed the Keanakāko'i Ash, which is stratigraphically above the 'Ailā'au flows, cannot be older than this age. Received: 10 October 1998 / Accepted: 12 May 1999     

降温过程对粗砂土中甲烷水合物形成的影响

多孔介质中甲烷水合物的形成和分解研究对于了解自然界多年冻土中的天然气水合物具有重要意义.本文通过3组10次甲烷水合物形成实验,研究了降温过程对粗砂土中甲烷水合物形成的影响.实验结果表明,降温过程对粗砂土中甲烷水合物的形成过程有较大的影响,在反应釜内温度降到0 ℃之前约90%的甲烷水合物已经形成,且降温速率越慢,水转化为甲烷水合物的转化率越高.0 ℃以下冻结过程对甲烷水合物的形成基本没有影响,但在冻结土体开始融化时发现有水合物二次生成的现象.     

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.     

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.     

Time-dependent analytic solutions of quasi-steady shocks with cooling

The author presents time-dependent analytical solutions of quasi-steady shocks with cooling, where quasi-steady shocks are objects composed of truncated steady-state models of shocks at any intermediate time. These solutions are compared to simulations with a hydrodynamical code and finally quasi-steady shocks are discussed as approximations to time-dependent shocks. Large departure of both the adiabatic and steady-state approximations from the quasi-steady solution emphasise the importance of the cooling history in determining the trajectory of a shock.     

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.     

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.