Replenishment rates of crustal magma and their bearing on potential sources of thermal energy

A general model of magma intrusion into the crust is developed which is based on a viscous-dissipation, forced-convection flow process driven by gravitational-buoyancy forces. Although some of the points in this general model have been studied before, it is possible with the present model to go further and calculate magma volumetric intrusion rates from fundamental properties and parameters. Equations for forced convection in a conduit with viscous dissipation are combined with results for the temperature dependence of magma viscosity. The volumetric intrusion rate is shown to be not a function of viscosity as might be expected, but rather a function primarily of the rate of change of viscosity with temperature. The model predictions for intrusion rate correlate well with field results for several sites where data exist for both intrusion or extrusion rate and for the temperature-dependent behavior of magma viscosity. The model predicts magma chamber replenishment rates equivalent to thermal energy rates on the order 10 GW (gigawatts) for a single active magma site. Assuming active magma sites on a 50-km spacing along volcanic lineaments leads to an estimate of a renewable magma intrusion rate into the crust of the western U.S. on the order of 2 TW (terawatts).     

Nonequilibrium flow in volcanic conduits and application to the eruptions of Mt. St. Helens on May 18, 1980, and Vesuvius in AD 79

A steady-state, one-dimensional, and nonhomogeneous two-phase flow model was developed for the prediction of local flow properties in volcanic conduits. The model incorporates the effects of relative velocity between the phases and for the variable magma viscosity. The resulting set of nonlinear differential equations was solved by a stiff numerical solver and the results were verified with the results of basaltic fissure eruptions obtained by a homogeneous two-phase flow model, before applying the model to the eruptions of Mt. St. Helens and Vesuvius volcanoes. This verification, and a study of the sensitivity of several modeling parameters, proved effective in establishing the confidence in the predicted nonequilibrium results of flow distribution in the conduits when the mass flow rate is critical or maximum. The application of the model to the plinian eruptions of Mt. St. Helens on May 18, 1980, and Vesuvius in AD 79, demonstrates the sensitivity of the magma discharge rate and distributions of pressure, volumetric fraction, and velocities of phases, on the hydrous magma viscosity feeding the volcanic conduits. Larger magma viscosities produce smaller mass discharge rates (or greater conduit diameters), smaller exit pressures, larger disequilibrium between the phases, and larger difference between the local lithostatic and fluid pressures in the conduit. This large pressure difference occurs when magma fragments and may cause a rupture of the conduit wall rocks, producing a closure of the conduit and cessation of the volcanic eruption, or water pouring into the conduit from underground aquifers leading to phreatomagmatic explosions. The motion of the magma fragmentation zone along a conduit during an eruption can be caused by the varying viscosity of magma feeding the volcanic conduit and may cause intermittent phreatomagmatic explosions during the plinian phases as different underground aquifers are activated at different depths. The variation of magma viscosity during the eruptions of Mt. St. Helens in 1980 and Vesuvius in AD 79 is normally associated with the tapping of magmas from different depths of the magma chambers. This variation of viscosity, which can include different crystal and dissolved water contents, can also produce conduit wall erosion, the onset and collapse of volcanic columns above the vent, and the onset and cessation of pyroclastic flows and surges.     

Eruptive style of the young high-Mg basaltic-andesite Pelagatos scoria cone,southeast of México City

The eruption of the Pelagatos scoria cone in the Sierra Chichinautzin monogenetic field near the southern suburbs of Mexico City occurred less than 14,000 years ago. The eruption initiated at a fissure with an effusive phase that formed a 7-km-long lava flow, and continued with a phase of alternating and/or simultaneous explosive and effusive activity that built a 50-m-high scoria cone on the western end of the fissure and formed a compound lava flow-field near the vent. The eruption ended with the emplacement of a short lava flow that breached the cone and was accompanied by weak explosions at the crater. Products consist of a microlite-rich high-Mg basaltic andesite. Samples were analyzed to determine the magma’s initial properties as well as the effects of degassing-induced crystallization on eruptive style. Although distal ash fallout deposits from this eruption are not preserved, a recent quarry exposes a large section of the scoria cone. Detailed study of exposed layers allows us to elucidate the mode of cone-building activity. Petrological and textural data, combined with models calibrated by experimental work and melt-inclusion analyses of similar magmas elsewhere, indicate that the magma was initially hot (>1,200°C), gas-rich (up to 5 wt.% H₂O), crystal-poor (~10 vol.% Fo₉₀ olivine phenocrysts) and thus poorly viscous (40–80 Pa s). During the early phase, low magma ascent velocity at the fissure vent allowed low-viscosity magma to degas and crystallize during ascent, producing lava flows with elevated crystal contents at T < 1,100°C, and blocky surfaces. Later, the closure of the fissure by cooling dikes focused the magma flow at a narrow section of the fissure. This led to an increased magma ascent velocity. Rapid and shallow degassing (<3 km deep) triggered ~40 vol.% microlite crystallization. Limited times for gas-escape and higher magma viscosity (6 × 10⁵–4 × 10⁶ Pa s) drove strong explosions of highly (60–80 vol.%) and finely vesicular magma. Coarse clasts broke on landing, which implies brittle behavior due to complete solidification. This requires sufficient time to cool and in turn implies ejection heights of over 1 km, which is much higher than “normal” Strombolian activity. Hence, magma viscosity significantly impacts eruption style at monogenetic volcanoes because it affects the kinetics of shallow degassing. The long-lasting eruptions of Jorullo and Paricutin, which produced similar magmas in western México, were more explosive. This can be related to higher magma fluxes and total erupted volumes. Implications of this study are important because basaltic andesites are commonly erupted to form monogenetic scoria cones of the Trans-Mexican Volcanic Belt.     

Experimental study of dense pyroclastic density currents using sustained,gas-fluidized granular flows

Submarine lava flow morphology is commonly used to estimate relative flow velocity, but the effects of crystallinity and viscosity are rarely considered. We use digital petrography and quantitative textural analysis techniques to determine the crystallinity of submarine basaltic lava flows, using a set of samples from previously mapped lava flow fields at the hotspot-affected Galápagos Spreading Center. Crystallinity measurements were incorporated into predictive models of suspension rheology to characterize lava flow consistency and rheology. Petrologic data were integrated to estimate bulk lava viscosity. We compared the crystallinity and viscosity of each sample with its flow morphology to determine their respective roles in submarine lava emplacement dynamics. We find no correlation between crystallinity, bulk viscosity, and lava morphology, implying that flow advance rate is the primary control on submarine lava morphology. However, we show systematic variations in crystal size and shape distribution among pillows, lobates, and sheets, suggesting that these parameters are important indicators of eruption processes. Finally, we compared the characteristics of lavas from two different sampling sites with contrasting long-term magma supply rates. Differences between lavas from each study site illustrate the significant effect of magma supply on the physical properties of the oceanic upper crust.     

Modelling the dynamics and thermodynamics of volcanic degassing

 The rates of passive degassing from volcanoes are investigated by modelling the convective overturn of dense degassed and less dense gas-rich magmas in a vertical conduit linking a shallow degassing zone with a deep magma chamber. Laboratory experiments are used to constrain our theoretical model of the overturn rate and to elaborate on the model of this process presented by Kazahaya et al. (1994). We also introduce the effects of a CO₂–saturated deep chamber and adiabatic cooling of ascending magma. We find that overturn occurs by concentric flow of the magmas along the conduit, although the details of the flow depend on the magmas' viscosity ratio. Where convective overturn limits the supply of gas-rich magma, then the gas emission rate is proportional to the flow rate of the overturning magmas (proportional to the density difference driving convection, the conduit radius to the fourth power, and inversely proportional to the degassed magma viscosity) and the mass fraction of water that is degassed. Efficient degassing enhances the density difference but increases the magma viscosity, and this dampens convection. Two degassing volcanoes were modelled. At Stromboli, assuming a 2 km deep, 30% crystalline basaltic chamber, containing 0.5 wt.% dissolved water, the ∼700 kg s^–1 magmatic water flux can be modelled with a 4–10 m radius conduit, degassing 20–100% of the available water and all of the 1 to 4 vol.% CO₂ chamber gas. At Mount St. Helens in June 1980, assuming a 7 km deep, 39% crystalline dacitic chamber, containing 4.6 wt.% dissolved water, the ∼500 kg s^–1 magmatic water flux can be modelled with a 22–60 m radius conduit, degassing ∼2–90% of the available water and all of the 0.1 to 3 vol.% CO₂ chamber gas. The range of these results is consistent with previous models and observations. Convection driven by degassing provides a plausible mechanism for transferring volatiles from deep magma chambers to the atmosphere, and it can explain the gas fluxes measured at many persistently active volcanoes. Received: 26 September 1997 / Accepted: 11 July 1998     

Volatile-rich magma injection into the feeding system during the 2001 eruption of Mt. Etna (Italy): its role on explosive activity and change in rheology of lavas

The explosive behavior and the rheology of lavas in basaltic volcanoes, usually driven by differentiation, can also be significantly affected by the kinetics of magma degassing in the upper portion of the feeding system. The complex eruption of 2001 at Mt. Etna, Italy, was marked by two crucial phenomena that occurred at the Laghetto vent on the southern flank of the volcano: 1) intense explosive activity and 2) at the end of the eruption, emission of a lava flow with higher viscosity than flows previously emitted from the same vent. Here, we investigate the hypothesis that these events were driven by the injection of volatile-rich magma into the feeding system. The input and mixing of this magma into a reservoir containing more evolved magma had the twofold effect of increasing 1) the overall concentration of volatiles and 2) their exsolution with consequent efficient vesiculation and degassing. This led to an explosive stage of the eruption, which produced a ~75-m-high cinder cone. Efficient volatile loss and the consequent increase of the liquidus temperature brought about the nucleation of Fe-oxides and other anhydrous crystalline phases, which significantly increased the magma viscosity in the upper part of the conduit, leading to the emission of a high viscosity lava flow that ended the eruption. The 2001 eruption has offered the opportunity to investigate the important role that input of volatile-rich magma may exert in controlling not only the geochemical features of erupted lavas but also the eruption dynamics. These results present a new idea for interpreting similar eruptions in other basaltic volcanoes and explaining eruptions with uncommonly high explosivity when only basic magmas are involved.     

Magma dynamics, crystallization, and chemical differentiation of the 1959 Kilauea Iki lava lake, Hawaii, revisited

Using constraints from an extensive database of geological and geochemical observations along with results from fluid mechanical studies of convection in magma chambers, we identify the main physical processes at work during the solidification of the 1959 Kilauea Iki lava lakes. In turn, we investigate their quantitative influence on the crystallization and chemical differentiation of the magma, and on the development of the internal structure of the lava lake. In contrast to previous studies, vigorous stirring in the magma, driven predominately by the descent of dense crystal-laden thermal plumes from the roof solidification front and the ascent of buoyant compositional plumes due to the in situ growth of olivine crystals at the floor, is predicted to have been an inevitable consequence of very strong cooling at the roof and floor. The flow is expected to have caused extensive but imperfect mixing over most of the cooling history of the magma, producing minor compositional stratification at the roof and thermal stratification at the floor. The efficient stirring of the large roof cooling is expected to have resulted in significant internal nucleation of olivine crystals, which ultimately settled to the floor. Additional forcing due to either crystal sedimentation or the ascent of gas bubbles is not expected to have increased significantly the amount of mixing. In addition to convection in the magma, circulation driven by the convection of buoyant interstitial melt in highly permeable crystal-melt mushes forming the roof and the floor of the lava lake is envisaged to have produced a net upward flow of evolved magma from the floor during solidification. In the floor zone, mush convection may have caused the formation of axisymmetric chimneys through which evolved magma drained from deep within the floor into the overlying magma and potentially the roof. We hypothesize that the highly evolved, pipe-like ‘vertical olivine-rich bodies’ (VORBs) [Bull. Volcanol. 43 (1980) 675] observed in the floor zone, of the lake are fossil chimneys. In the roof zone, buoyant residual liquid both produced at the roof solidification front and gained from the floor as a result of incomplete convective mixing is envisaged to have percolated or ‘leaked‘ into the overlying highly-permeable cumulate, displacing less buoyant interstitial melt downward. The results from Rayleigh fractionation-type models formulated using boundary conditions based on a quantitative understanding of the convection in the magma indicate that most of the incompatible element variation over the height of the lake can be explained as a consequence of a combination of crystal settling and the extensive but imperfect convective mixing of buoyant residual liquid released from the floor solidification front. The remaining chemical variation is understood in terms of the additional influences of mush convection in the roof and floor on the vertical distribution of incompatible elements. Although cooling was concentrated at the roof of the lake, the floor zone is found to be thicker than the roof zone, implying that it grew more quickly. The large growth rate of the floor is explained as a consequence of a combination of the substantial sedimentation of olivine crystals and more rapid in situ crystallization due to both a higher liquidus temperature and enhanced cooling resulting from imperfect thermal and chemical mixing.     

A preliminary study of the effects of some flow parameters in the generation of poloidal and toroidal energies within a 3-D spherical thermal-convective system with variable viscosity

The effects of variable viscosity on flow dynamics within spherical shells are investigated using a finite-element thermal convection model, and preliminary result for cases with relatively low Rayleigh numbers and small viscosity contrasts are reported. These results demonstrate some general effects of viscosity variation on mantle dynamics, and, in particular, the generation of toroidal energy. Since lateral viscosity variations are necessary in the generation of toroidal motion in a thermally driven convective system, it is not surprising our results show that flows with greater viscosity contrasts produce greater amounts of toroidal energy. Our preliminary study further shows that solutions become more time-dependent as viscosity contrasts increase. Increasing the Rayleigh number is also found to increase the magnitude of toroidal energy. Internal heating, on the other hand, appears to lead to less toroidal energy compared wth bottom heating because it tends to produce a thermally more uniform interior and thus smaller viscosity variations.     

Some thermal problems associated with magma migration

The mechanisms by which magma migrates from the point in the earth's interior where melting occurs to the earth's surface are poorly understood. In this paper several aspects of this problem are examined. Magma can migrate upward due to its differential buoyancy on the scale of crystalline grains or as large dispairs. Magma transport is an effective means of heat transport. Magma transport at a rate of 0.15 cm/yr is equivalent to a heat flow of 10^-6 cal/cm² s. If magma encounters country rock witha lower melting point the original magma is likely to solidify while melting the country rock. This would be an effective mechanism of purging silicic rocks and incomparable elements from the lower crust. Under some circumstances magma must penetrate up to 100 km or more of cold lithospheric rock. In order for this magma to reach the surface without solidification a heated path must be provided. The heating of this path requires the solidification of some magma. It is estimated that magma penetrates the lithosphere in about 5000 years and that the crack is lined by several hundred meters of frozen basalt.     

An experimental study on the formation of composite intrusions from zoned magma chambers

We present a model which accounts for the common, but paradoxical arrangement of composite intrusions (i.e. silicic core and mafic margins) on the basis of analogue experiments using gelatin and aqueous solutions. The present model involves simultaneous flow-out of the upper and lower magmas from a longitudinal crack along the chamber wall. Experimental results suggest that the mafic magma from the lower layer leaks from the side-wall of the chamber and travels faster than the silicic magma because of its lower viscosity, so that the mafic magma reaches the tip of the crack first. Once the mafic magma reaches the crack tip, then the rate of dyke propagation becomes determined by the viscosity of the less viscous mafic magma, and so it can advance rapidly. The viscous silicic magma can flow efficiently into the center of the dyke, being lubricated by the mafic magma margins. This model accounts for the common arrangement of composite intrusions and gives an efficient mechanism of flow of viscous silicic magmas.     

Effects of planetary thermal structure on the ascent and cooling of magma on Venus

Magellan radar images of the surface of Venus show a spatially broad distribution of volcanic features. Models of magmatic ascent processes to planetary surfaces indicate that the thermal structure of the interior significantly influences the rate of magmatic cooling and thus the amount of magma that can be transported to the surface before solidification. In order to understand which aspects of planetary thermal structure have the greatest influence on the cooling of buoyantly ascending magma, we have constructed magma cooling profiles for a plutonic ascent mechanism, and evaluated the profiles for variations in the surface and mantle temperature, surface temperature gradient, and thermal gradient curvature. Results show that, for a wide variety of thermal conditions, smaller and slower magma bodies are capable of reaching the surface on Venus compared to Earth, primarily due to the higher surface temperature of Venus. Little to no effect on the cooling and transport of magma are found to result from elevated mantle temperatures, elevation-dependent surface temperature variations, or details of the thermal gradient curvature. The enhanced tendency of magma to reach the surface on Venus may provide at least a partial explanation for the extensive spatial distribution of observed volcanism on the surface.     

地壳中岩浆沿已有断层运移的数值模拟——在长白山天池火山的应用

本文利用有限元及有限差分的方法,对壳内岩浆房或岩浆囊中的岩浆在构造应力及由于围岩与岩浆的密度差产生的浮力作用下,沿已有断层向上运移的动力学过程进行了数值模拟.在岩浆囊顶部与上覆岩层接触处,沿着已有微小破裂,岩浆在一定超压力条件下使已有断层张开并继续向上延伸,从而形成岩浆向上运移的通道.研究了岩浆黏度、密度差、模型深度对最小超压力（岩浆运移到地表所需的最小岩浆房超压力）的影响.在10 km深度的地壳中,若岩浆黏度为0.1~10³ Pa·s,当超压力达到17~20 MPa时,岩浆压力可以驱动岩浆运移到地表层;同时,岩浆动力黏度越大,使岩脉运移到地表需要的超压力就越大.当密度差为300~700 kg·m^-3,其变化对超压力的影响比较小.本文亦对比了三维应力条件和二维平面应变条件下不同结果,比较了不同条件下岩浆运移造成的地表垂直位移变化.结合长白山天池火山地区的区域地质环境,对长白山天池火山岩浆运移条件进行了参数试验性计算分析,估算了在给定长白山天池火山模型条件下地下可能存在的岩浆囊的大小,其结果对认识长白山天池火山地区岩浆活动及相关的预测和监控有参考意义.     

Modeling of the ascent of magma during the plinian eruption of Vesuvius in A.D. 79

The ascent of magma during the A.D. 79 eruption of Vesuvius was studied by a steady-state, one-dimensional, and nonequilibrium two-phase flow model. The gas exsolution process was modeled by assuming a chemical equilibrium between the exsolved and dissolved gas, whereas the magma density and viscosity were modeled by accounting for the crystal content in magma. The exsolution, density, and viscosity models consider the effect of different compositions of the white and gray magmas. By specifying the conduit geometry and magma composition, and employing the model to search for the maximum discharge rate of magma which is consistent with the specified geometry and magma composition, the model was then used to establish the two-phase flow parameters along the conduit. It was found that for all considered conditions the magma pressure in the conduit decreases below the lithostatic pressure near the magma fragmentation level, and that in the deep regions of the conduit the white magma pressure is larger and the gray magma pressure is lower than the lithostatic one. The exsolution and fragmentation levels were found to be deeper for the white than for the gray magma, and the changing composition during the eruption causes an increase of the exit pressure and decrease of the exit gas volumetric fraction. The model also predicted a minimum conduit diameter which is consistent with the white and gray magma compositions and mass flow-rates. The predictions of the model were shown to be consistent with column collapses during the gray eruption phase, large presence of carbonate lithics in the gray pumice fall deposit, and magma-water interaction during a late stage of the eruption.     

The influence of geometry on the ascent of magma in open fissures

During steady eruption, the flow conditions (emitted mass flux, exit velocity and exit pressure) depend on the geometry of the conduit in which the eruption occurs. This dependence is examined for the onedimensional, isothermal ascent of a homogeneous basaltic magma with an aqueous volatile phase and newtonian rheology. By fixing the geometry of the feeding fissure, the mass flux flowing in steady conditions can be determined at any depth, as well the magma pressure and vertical velocity. Flow behaviour is analysed for three fissure shapes: constant width, slowly upward narrowing and lenticular. In all the cases examined the magma arrives at the earth's surface with a pressure greater than atmospheric. The results are compared with those obtained when a lithostatic pressure gradient is assumed for the magma column. Some speculations are made, moreover, about the change in eruption style, if conduit geometry varies during a non-steady phase.     

Degassing of stratified magma by compositional convection

A new model is proposed for passive degassing from sub-volcanic magma chambers. The water content in stably stratified shallow magma chamber will be equated to its solubility at the upper boundary by convection. Water from a lower layer high in water content can enrich the contact zone of the upper layer and lead to further convective overturn of this boundary layer. A complete set of equations describing convection with bubble formation and dissolution is reduced to a simplified form by assuming a small bubble content. The development and pattern of flow driven by vesiculation is modeled numerically in a 2D magma chamber for relatively low Raleigh numbers (5×10⁵). Bubbles rising from the magma will collect near the roof in a layer of 8–10 vol% and then escape upward to fumaroles. The Stokes flux of bubbles escaping from an andesitic magma with viscosity 10⁴ P and a top surface of about 500×500 m corresponds with observed total magmatic water fluxes of 35 kg/s. Pressure within the chamber is buffered by elastic (and local visco-elastic) deformations in the solid rocks bounding the chamber to the range between ambient and close to lithostatic values. In a chamber closed to fresh magma inputs, the decrease in volume due to such gentle volatile escape lowers the reference pressure. Bubbles flux from the lower layer induced by variation of the saturation level around stratification boundary may be efficient mechanism for the water transport between layers.     

Volumetric characteristics of lava flows from interferometric radar and multispectral satellite data: the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos

We have used a suite of remotely sensed data, numerical lava flow modeling, and field observations to determine quantitative characteristics of the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos Islands. Flank lava flow areas, volumes, instantaneous effusion rates, and average effusion rates were all determined for these two eruptions, for which only limited syn-eruptive field observations are available. Using data from SPOT, TOPSAR, ERS-1, and ERS-2, we determined that the 1995 Fernandina flow covers a subaerial area of 6.5×10⁶ m² and has a subaerial dense rock equivalent (DRE) volume of 42×10⁶ m³. Field observations, ATSR satellite data, and the FLOWGO numerical model allow us to determine that the effusion rate declined exponentially from a high of ~60–200 m³ s^-1 during the first few hours to <5 m³ s^-1 prior to ceasing after 73 days, with a mean effusion rate of 4–16 m³ s^-1. Integrating the ATSR-derived, exponentially declining effusion rate over the eruption duration produces a total (subaerial + submarine) DRE volume of between 27 and 100×10⁶ m³, the range in values being due to differing assumptions about heat loss characteristics; only values in the higher part of this range are consistent with the independently derived subaerial volume. Using SPOT, TOPSAR, ERS-1, and ERS-2 data, we determine that the 1998 Cerro Azul flow is 16 km long, covers 16 km², and has a DRE volume of 54×10⁶ m³. FLOWGO produces at-vent velocity and effusion rate values of 11 m s^-1 and ~600 m³ s^-1, respectively. The velocity value agrees well with the 12 m s^-1 estimated in the field. The mean effusion rate (total DRE volume/duration) was 7–47 m³ s^-1. Dike dimensions, fissure lengths, and pressure gradients along the conduit based on magma chamber depth estimates of 3–5 km produce mean effusion rates for the two eruptions that range over nearly four orders of magnitude, the range being due to uncertainty in the magma viscosity, dike dimensions, and pressure gradient between magma chamber and vent. Although somewhat consistent with mean effusion rates from other techniques, their wide range makes them less useful. The exponentially declining effusion rates during both eruptions are consistent with release of elastic strain being the driving mechanism of the eruptions. Our results provide independent input parameters for previously published theoretical relationships between magma chamber pressurization and eruption rates that constrain chamber volumes and increases in volume prior to eruption, as well as time constants of exponential decay during the eruption. The results and theoretical relationships combine to indicate that at both volcanoes probably 25–30% of the volumetric increase in the magma chamber erupted as lava onto the surface. In both eruptions the lava flow volumes are less than 1% of the magma chamber volume.     

Conditions for the arrest of a vertical propagating dyke

Magma ascent towards the Earth’s surface occurs through dyke propagation in the vast majority of cases. We investigate two purely mechanical effects unrelated to cooling or solidification that lead to the arrest of propagation, so that no eruption occurs. The first is that the input of magma from the source is not maintained continuously, such that a fixed volume of magma is released. Laboratory experiments show that, in this case, the dyke stops at a finite distance from the source. This behaviour is specific to the fracturing process in 3-D. We derive a relationship for the minimum magma volume required for an eruption as a function of magma buoyancy and source depth. When large magma volumes are available, eruption may also be prevented by a thick low density layer in the upper crust. Numerical studies of dyke propagation show that the dyke continues to rise even though it is negatively buoyant. Magma accumulates in a swollen nose region at the interface between the low density layer and the dense basement. Magma overpressure is largest at this interface and increases with increasing penetration into the upper layer. It may become large enough to induce horizontal fractures in the dyke walls and lateral intrusion of a sill, which prevents eruption. This requires that the thickness of the low density layer exceeds a threshold value that depends on the density contrast between magma and host rock. If the magma volume is smaller than a threshold value, neither sill intrusion nor eruption are possible and magma gets stored in a horizontal blade-shaped dyke straddling the interface. Scaling laws for variations of ascent rate and for the minimum magma volume allow diagnosis of a failed eruption.     

Toothpaste lava: Characteristics and origin of a lava structural type transitional between pahoehoe and aa

Toothpaste lava, an important basalt structural type which illustrates the transition from pahoehoe to aa, is particularly well displayed on the 1960 Kapoho lava of Kilauea Volcano. Its transitional features stem from a viscosity higher than that of pahoehoe and a rate of flow slower than that of aa. Viscosity can be quantified by the limited settling of olivine phenocrysts and rate of flow by field observations related to the low-angle slope on which the lava flowed. Much can be learned about the viscosity, rheologic condition, and flow velocity of lavas long after solidification by analyses of their structural characteristics, and it is possible to make at least a semiquantitative assessment of the numerical values of these parameters.