Explosive Volcanic Eruptions from Linear Vents on Earth, Venus and Mars: Comparisons with Circular Vent Eruptions

Conditions required to support buoyant convective plumes are investigated for explosive volcanic eruptions from circular and linear vents on Earth, Venus, and Mars. Vent geometry (linear versus circular) plays a significant role in the ability of an explosive eruption to sustain a buoyant plume. On Earth, linear and circular vent eruptions are both capable of driving buoyant plumes to equivalent maximum rise heights, however, linear vent plumes are more sensitive to vent size. For analogous mass eruption rates, linear vent plumes surpass circular vent plumes in entrainment efficiency approximately when L(sub o) > 3r(sub o) owing to the larger entrainment area relative to the control volume. Relative to circular vents, linear vents on Venus favor column collapse and the formation of pyroclastic flows because the range of conditions required to establish and sustain buoyancy is narrow. When buoyancy can be sustained, however, maximum plume heights exceed those from circular vents. For current atmospheric conditions on Mars, linear vent eruptions are capable of injecting volcanic material slightly higher than analogous circular vent eruptions. However, both geometries are more likely to produce pyroclastic fountains, as opposed to convective plumes, owing to the low density atmosphere. Due to the atmospheric density profile and water content on Earth, explosive eruptions enjoy favorable conditions for producing sustained buoyant columns, while pyroclastic flows would be relatively more prevalent on Venus and Mars. These results have implications for the injection and dispersal of particulates into the planetary atmosphere and the ability to interpret the geologic record of planetary volcanism

Simulation of Cooling and Pressure Effects on Inflated Pahoehoe Lava Flows

Pahoehoe lobes are often emplaced by the advance of discrete toes accompanied by inflation of the lobe surface. Many random effects complicate modeling lobe emplacement, such as the location and orientation of toe breakouts, their dimensions, mechanical strength of the crust, micro-topography and a host of other factors. Models that treat the movement of lava parcels as a random walk have explained some of the overall features of emplacement. However, cooling of the surface and internal pressurization of the fluid interior has not been modeled. This work reports lobe simulations that explicitly incorporate 1) cooling of surface lava parcels, 2) the propensity of breakouts to occur at warmer margins that are mechanically weaker than cooler ones, and 3) the influence of internal pressurization associated with inflation. The surface temperature is interpreted as a surrogate for the mechanic strength of the crust at each location and is used to determine the probability of a lava parcel transfer from that location. When only surface temperature is considered, the morphology and dimensions of simulated lobes are indistinguishable from equiprobable simulations. However, inflation within a lobe transmits pressure to all connected fluid locations with the warmer margins being most susceptible to breakouts and expansion. Simulations accounting for internal pressurization feature morphologies and dimensions that are dramatically different from the equiprobable and temperature-dependent models. Even on flat subsurfaces the pressure-dependent model produces elongate lobes with distinct directionality. Observables such as topographic profiles, aspect ratios, and maximum extents should be readily distinguishable in the field

Topographic and Stochastic Influences on Pahoehoe Lava Lobe Emplacement

A detailed understanding of phoehoe emplacement is necessary for developing accurate models of flow field development, assessing hazards, and interpreting the significance of lava morphology on Earth and other planetary surfaces. Active pahoehoe lobes on Kilauea Volcano, Hawaii, were examined on 21-26 February 2006 using oblique time-series stereo-photogrammetry and differential global positioning system (DGPS) measurements. During this time, the local discharge rate for peripheral lava lobes was generally constant at 0.0061 +/- 0.0019 m3/s, but the areal coverage rate of the lobes exhibited a periodic increase every 4.13 +/- 0.64 minutes. This periodicity is attributed to the time required for the pressure within the liquid lava core to exceed the cooling induced strength of its margins. The pahoehoe flow advanced through a series of down slope and cross-slope breakouts, which began as approximately 0.2 m-thick units (i.e., toes) that coalesced and inflated to become approximately meter-thick lobes. The lobes were thickest above the lowest points of the initial topography and above shallow to reverse facing slopes, defined relative to the local flow direction. The flow path was typically controlled by high-standing topography, with the zone directly adjacent to the final lobe margin having an average relief that was a few centimeters higher than the lava inundated region. This suggests that toe-scale topography can, at least temporarily, exert strong controls on pahoehoe flow paths by impeding the lateral spreading of the lobe. Observed cycles of enhanced areal spreading and inflated lobe morphology are also explored using a model that considers the statistical likelihood of sequential breakouts from active flow margins and the effects of topographic barriers

The Influence of Slope Breaks on Lava Flow Surface Disruption

Changes in the underlying slope of a lava flow impart a significant fraction of rotational energy beyond the slope break. The eddies, circulation and vortices caused by this rotational energy can disrupt the flow surface, having a significant impact on heat loss and thus the distance the flow can travel. A basic mechanics model is used to compute the rotational energy caused by a slope change. The gain in rotational energy is deposited into an eddy of radius R whose energy is dissipated as it travels downstream. A model of eddy friction with the ambient lava is used to compute the time-rate of energy dissipation. The key parameter of the dissipation rate is shown to be rho R(sup 2/)mu, where is the lava density and mu is the viscosity, which can vary by orders of magnitude for different flows. The potential spatial disruption of the lava flow surface is investigated by introducing steady-state models for the main flow beyond the steepening slope break. One model applies to slow-moving flows with both gravity and pressure as the driving forces. The other model applies to fast-moving, low-viscosity, turbulent flows. These models provide the flow velocity that establishes the downstream transport distance of disrupting eddies before they dissipate. The potential influence of slope breaks is discussed in connection with field studies of lava flows from the 1801 Hualalai and 1823 Keaiwa Kilauea, Hawaii, and 2004 Etna eruptions

Quantitative Studies in Planetary Volcanism

Scientific research was conducted on volcanic processes on Mars, Venus, Io, the moon, and the Earth. The achievements led to scientific advances in the understanding of volcanic plumes, lava flow emplacements, coronae, and regoliths on the solid surfaces. This research led to multiple publications on each of the main topics of the proposal. Research was also presented at the annual Lunar and Planetary Science Conference at Houston. Typically, this grant contributed to 3-4 presentations each year. This grant demonstrated, numerous times, the usefulness of NASA mission data for advancing the understanding of volcanic processes on other planetary surfaces and the Earth

Topographic Effects on Geologic Mass Movements

This report describes research directed toward understanding the response of volcanic lahars and lava flows to changes in the topography along the path of the flow. We have used a variety of steady-state and time-dependent models of lahars and lava flows to calculate the changes in flow dynamics due to variable topography. These models are based on first-order partial differential equations for the local conservation of volume. A global volume conservation requirement is also imposed to determine the extent of the flow as a function of time and the advance rate. Simulated DEMs have been used in this report

Constraints on Determining the Eruption Style and Composition of Terrestrial Lavas from Space

The surface temperatures of active lavas relate to cooling rates, chemistry, and eruption style. We analyzed 61 hyperspectral satellite images acquired by the National Aeronautics and Space Administration s Earth Observing-1 (EO-1) Hyperion imaging spectrometer to document the surface temperature distributions of active lavas erupted at 13 volcanoes. Images were selected to encompass the range of common lava eruption styles, specifically, lava fountains, flows, lakes, and domes. Our results reveal temperature distributions for terrestrial lavas that correlate with composition (i.e., a statistically significant difference in the highest temperatures retrieved for mafic lavas and intermediate and felsic lavas) and eruption style. Maximum temperatures observed for mafi c lavas are approx.200 C higher than for intermediate and felsic lavas. All eruption styles exhibit a low-temperature mode at approx.300 C; lava fountains and 'a' a flows also exhibit a higher-temperature mode at approx.700 C. The observed differences between the temperatures are consistent with the contrasting rates at which the lava surfaces are thermally renewed. Eruption styles that allow persistent and pervasive thermal renewal of the lava surface (e.g., fractured crusts on channel-fed 'a' a flows) exhibit a bimodal temperature distribution; eruption styles that do not (e.g., the continuous skin of pahoehoe lavas) exhibit a single mode. We conclude that insights into composition and eruption style can only be gained remotely by analyzing a large spatio-temporal sample of data. This has implications for determining composition and eruption style at the Jovian moon Io, for which no in situ validation is available

Morphology and geometry of the distal ramparts of Martian impact craters

From the proceedings of the Workshop on the Role of Volatiles and Atmospheres on Martian Impact Craters held on July 11-14, 2005, at the Johns Hopkins University Applied Physics Laboratory.We used Mars Orbiter Laser Altimeter (MOLA), Thermal Emission Imaging System visible light (THEMIS VIS), and Mars Orbiter Camera (MOC) data to identify and characterize the morphology and geometry of the distal ramparts surrounding Martian craters. Such information is valuable for investigating the ejecta emplacement process, as well as searching for spatial variations in ejecta characteristics that may be due to target material properties and/or latitude, altitude, or temporal variations in the climate. We find no systematic trend in rampart height that would indicate regional variations in target properties for 54 ramparts at 37 different craters 5.7-35.9 km in diameter between 52.3 degrees S to 47.6 degrees N. Rampart heights for multi-lobe and single-lobe ejecta are each normally distributed with a common standard deviation, but statistically distinct mean values. Ramparts range in height from 20-180 m, are not symmetric, are typically steeper on their distal sides, and may be as much as ~4 km wide. The ejecta blanket proximal to parent crater from the rampart may be very thin (<5 m). A detailed analysis of two craters, Toconao crater (21 degrees S, 285 degrees E) (28 measurements), and an unnamed crater within Chryse Planitia (28.4 degrees N, 319.6 degrees E) (20 measurements), reveals that ejecta runout distance increases with an increase in height between the crater rim and the rampart, but that rampart height is not correlated with ejecta runout distance or the thickness of the ejecta blanket.The Meteoritics & Planetary Science archives are made available by the Meteoritical Society and the University of Arizona Libraries. Contact [email protected] for further information.Migrated from OJS platform February 202

Transport of atmospheric water vapor by volcanic eruption columns

Contrary to assumptions often made in the literature, explosive volcanic eruptions are capable of transporting significant amounts of water into the stratosphere. In addition to the magmatic water component, atmospheric water vapor is entrained by the column at lower levels. A theoretical model for the conservation of mass, momentum, and thermal energy of four separate components (dry air, water vapor, liquid condensates, and solid particles) is used to determine the extent of atmospheric water redistribution. We examine the effects of water vapor condensation on dynamical characteristics and ambient water vapor transport. A simple technique is presented for deriving canonical forms for the complex system of ordinary differential equations governing the column components. Solutions of this model are presented that show the influence of different volcanic boundary conditions and a range of ambient water vapor distributions on transport of the buoyant column. We show that the water component (vapor + liquid) of small eruption columns rising through a wet atmosphere is dominated by entrained water, whereas larger columns are dominated by the magmatic water. This is due, in part, to the proportionately smaller entrainment surface area in relation to the control volume for the larger columns. We also show that a maintained column with an initial mass flux of 2.7 × 108 Kg s-1 erupted into a wet atmosphere would inject 96 Mt of water vapor into the stratosphere over 24 hours, comparable to the annual input from methane oxidation or 100 midlatitude thunderstorms. This increase may accelerate the conversion of simultaneously erupted volcanic SO2 into sulfuric acid

Corrigendum to "New Approaches to Inferences for Steep-Sided Domes on Venus" [J. Volcanol. Geotherm. Res. 319 (2016) 93-105]

A typographical error contained in Quick et al. (2016) indicates the incorrect units for the value of the combined quantity (roh(exp. 3)o) that is the basis of Figs. 5, 6, and 7, and Tables 2 and 3. Using the values of ro and ho provided in Table 2, it can easily be shown that the combined quantity is correctly stated as roh(exp. 3)o =0.617 km(exp. 4). As correctly stated in Quick et al. (2016), the combined quantity of (roh(exp. 3)o) determines the family of curves shown in Fig. 5. The derivation of this relationship is shown below for completeness. Note that all results as reported in Quick et al. (2016) remained unchanged