This study has the objective to better constrain the role played by thermal erosion by turbulent lava in the formation of large channels on Mars and the Moon. On Mars, a rigorous one-dimensional model was used to test whether lava might have excavated the Athabasca Valles outflow channel. Calculated erosion depths are much lower than the measured depths of the channel, and suggest a limited role played by thermal erosion in excavating it. On the Moon, the investigation focused on the outer and inner sinuous rilles of Vallis Schröteri. At this site, erosional features cannot be explained by one- and two-dimensional models. The first 3-D model of thermal erosion by turbulent lava on the Moon was created to relate the spatial distribution of erosion rates over the bed and banks of a channel with changes in fluid- and thermodynamic parameters. The turbulence model chosen for each steady-state simulation is the Shear Stress Transport (SST) k-ω model and OpenFOAM is the Computational Fluid Dynamics software used. At the 150-km-long, 4-km-wide, and up-to 700-m-deep outer rille, I aimed to determine maximum erosion rates at/near the lava source and rille segments 1-km-long and 4-km-wide were chosen for the simulations. By adopting the obtained maximum erosion rates of 1 m/day, lava might have taken ~2 years to excavate the 700-m-deep depression. These fast erosion rates were unlikely maintained downstream of the lava source unless lava flowed in a tube. Besides, observational evidence suggests that tectonics and constructional processes likely contributed to rille development. On these grounds, thermal or thermo-mechanical erosion might have contributed to rille formation at a later stage. At the Vallis Schröteri inner rille, 1-km-long and 160-m-wide meandering channels were chosen. In one scenario, lava loses heat by radiation, in the other flows in a tube. Using the calculated (and conservative) erosion rate of 50 cm/day, it would have taken ~6 months for the 90-m deep inner rille to be excavated. A mechanism of secondary flow circulation analogous to that found in meandering rivers potentially explains meander generation. At each bend, downstream and cross-stream velocity variations lead to local temperature/ erosion enhancements.

Lava flow emplacement in the laboratory and on the surface of Mars was investigated. In the laboratory, the effects of unsteady effusion rates at the vent on four modes of emplacement common to lava flow propagation: resurfacing, marginal breakouts, inflation, and lava tubes was addressed. A total of 222 experiments were conducted using a programmable pump to inject dyed PEG wax into a chilled bath (~ 0° C) in tanks with a roughened base at slopes of 0, 7, 16, and 29°. The experiments were divided into four conditions, which featured increasing or decreasing eruption rates for either 10 or 50 s. The primary controls on modes of emplacement were crust formation, variability in the eruption rate, and duration of the pulsatory flow rate. Resurfacing – although a relatively minor process – is inhibited by an extensive, coherent crust. Inflation requires a competent, flexible crust. Tube formation requires a crust and intermediate to low effusion rates. On Mars, laboratory analogue experiments combined with models that use flow dimensions to estimate emplacement conditions and using high resolution image data and digital terrain models (e.g. THEMIS IR, CTX, HRSC), the eruption rates, viscosities, and yield strengths of 40 lava flows in the Tharsis Volcanic Province have been quantified. These lava flows have lengths, mean widths, and mean thicknesses of 15 – 314 km, 0.5 – 29 km, and 11 – 91 m, respectively. Flow volumes range from ~1 – 430 km3. Based on laboratory experiments, the 40 observed lava flows were erupted at 0.2 – 6.5x103 m3/s, while the Graetz number and Jeffrey’s equation when applied to 34 of 40 lava flows indicates eruption rates and viscosities of 300 – ~3.5 x 104 m3/s and ~105 – 108 Pa s, respectively. Another model which accounts for mass loss to levee formation was applied to a subset of flows, n = 13, and suggests eruption rates and viscosities of ~30 – ~1.2 x 103 m3/s and 4.5 x 106 – ~3 x 107 Pa s, respectively. Emplacement times range from days to centuries indicating the necessity for long-term subsurface conduits capable of delivering enormous volumes of lava to the surface.

Volcanoes can experience multiple eruption styles throughout their eruptive histories. Among the most complex and most common are eruptions of intermediate explosivity, such as Vulcanian and sub-Plinian eruptions. Vulcanian eruptions are characterized by small-scale, short-lived, ash-rich eruptions initiated by the failure of a dense magma plug or overlying dome that had sealed an overpressured conduit. Sub-Plinian eruptions are characterized by sustained columns that reach tens of kilometers in height.

Multiple eruption styles can be observed in a single eruptive sequence. In recent decades, transitions in eruption style during well-documented eruptions have been described in detail, with some workers proposing explanatory mechanisms for the transitions. These proposed mechanisms may be broadly classified into processes at depth, processes in the conduit, or some combination of both.

The present study is focused on the Pietre Cotte sequence because it may have encompassed up to three different eruptive cycles, each representing different degrees of explosivity. The first deposits are composed of repeated layers of fine ash and lapilli composed of latite and rhyolite endmembers, efficiently mixed at sub-cm scales. The thin layers and bubble/crystal textures indicate that the magma underwent numerous decompressions and open-system degassing, and that the eruptions waned with time. The second phase of the sequence appears to have been initiated by cm-scale mixing between a volatile-rich, mafic magma from deeper in the system and a shallow silicic body. Textures indicate that the magma ascended rapidly and experienced little to no open-system degassing. The final phase of the sequence again produced repeated layers of fine ash and lapilli, of uniform trachyte composition, and waned with time. The first and last phases were likely produced in Vulcanian eruptions, while the pumice-rich layers were likely produced in Vulcanian to sub-Plinian eruptions.

In summary, the Pietre Cotte sequence is characterized by up to three magma recharge events in ~200 years. The differences in eruptive style appear to have been controlled by variations in the volatile content of the recharge magma, as well as the efficiency and scale of magma mixing and resulting overpressures.

The seasonal deposition of CO2 on the polar caps is one of the most dynamic processes on Mars and is a dominant driver of the global climate. Remote sensing temperature and albedo data were used to estimate the subliming mass of CO2 ice on south polar gullies near Sisyphi Cavi. Results showed that column mass abundances range from 400 - 1000 kg.m2 in an area less than 60 km2 in late winter. Complete sublimation of the seasonal caps may occur later than estimated by large-scale studies and is geographically dependent. Seasonal ice depth estimates suggested variations of up to 1.5 m in depth or 75% in porosity at any one time. Interannual variations in these data appeared to correlate with dust activity in the southern hemisphere. Correlation coefficients were used to investigate the relationship between frost-free surface properties and the evolution of the seasonal ice in this region. Ice on high thermal inertia units was found to disappear before any other ice, likely caused by inhibited deposition during fall. Seasonal ice springtime albedo appeared to be predominantly controlled by orientation, with north-facing slopes undergoing brightening initially in spring, then subliming before south-facing slopes. Overall, the state of seasonal ice is far more complex than globally and regionally averaged studies can identify.

The discovery of cryovolcanic features on Charon and the presence of ammonia hydrates on the surfaces of other medium-sized Kuiper Belt Objects suggests that cryovolcanism may be important to their evolution. A two-dimensional, center-point finite difference, thermal hydraulic model was developed to explore the behavior of cryovolcanic conduits on midsized KBOs. Conduits on a Charon-surrogate were shown to maintain flow through over 200 km of crust and mantle down to radii of R = 0.20 m. Radii higher than this became turbulent due to high viscous dissipation and low thermal conductivity. This model was adapted to explore the emplacement of Kubrik Mons. Steady state flow was achieved with a conduit of radius R = 0.02 m for a source chamber at 2.3 km depth. Effusion rates computed from this estimated a 122 - 163 Myr upper limit formation timescale.

Previous workers hypothesized that lunar Localized Pyroclastic Deposits (LPDs) represent products of vulcanian-style eruptions, since some have low proportions of juvenile material. The objective of the first study is to determine how juvenile composition, calculated using deposit and vent volumes, varies among LPDs. I used Lunar Reconnaissance Orbiter Camera Narrow Angle Camera (LROC NAC) digital terrain models (DTMs) to generate models of pre-eruption surfaces for 23 LPDs and subtracted them from the NAC DTMs to calculate deposit and vent volumes. Results show that LPDs have a wide range of juvenile compositions and thinning profiles, and that there is a positive relationship between juvenile material proportion and deposit size. These findings indicate there is greater diversity among LPDs than previously understood, and that a simple vulcanian eruption model may only apply to the smallest deposits.

There is consensus that martian outflow channels were formed by catastrophic flooding events, yet many of these channels exhibit lava flow features issuing from the same source as the eroded channels, leading some authors to suggest that lava may have served as their sole agent of erosion. This debate is addressed in two studies that use Context Camera images for photogeologic analysis, geomorphic mapping, and cratering statistics: (1) A study of Mangala Valles showing that it underwent at least two episodes of fluvial activity and at least three episodes of volcanic activity during the Late Amazonian, consistent with alternating episodes of flooding and volcanism. (2) A study of Maja Valles finds that it is thinly draped in lava flows sourced from Lunae Planum to the west, rendering it analogous to the lava-coated Elysium outflow systems. However, the source of eroded channels in Maja Valles is not the source of the its lava flows, which instead issue from south Lunae Planum. The failure of these lava flows to generate any major channels along their path suggests that the channels of Maja Valles are not lava-eroded.

Finally, I describe a method of locating sharp edges in out-of-focus images for application to automated trajectory control systems that use images from fixed-focus cameras to determine proximity to a target.

Many planetary science missions study thermophysical properties of surfaces using infrared spectrometers and infrared cameras. Thermal inertia is a frequently derived thermophysical property that quantifies the ability for heat to exchange through planetary surfaces.

To conceptualize thermal inertia, the diffusion equation analogies are extended using a general effusivity term: the square root of a product of conductivity and capacity terms. A hypothetical thermal inductance was investigated for diurnal planetary heating. The hyperbolic heat diffusion equation was solved to derive an augmented thermal inertia. The hypothetical thermal inductance was modeled with negligible effect on Mars.

Extending spectral performance of infrared cameras was desired for colder bodies in the outer solar system where peak infrared emission is at longer wavelengths. The far-infrared response of an infrared microbolometer array with a retrofitted diamond window was determined using an OSIRIS-REx—OTES interferometer. An instrument response function of the diamond interferometer-microbolometer system shows extended peak performance from 15 µm out to 20 µm and 40% performance to at least 30 µm. The results are folded into E-THEMIS for the NASA flagship mission: Europa Clipper.

Infrared camera systems are desired for the expanding smallsat community that can inherit risk and relax performance requirements. The Thermal-camera for Exploration, Science, and Imaging Spacecraft (THESIS) was developed for the Prox-1 microsat mission. THESIS, incorporating 2001 Mars Odyssey—THEMIS experience, consists of an infrared camera, a visible camera, and an instrument computer. THESIS was planned to provide images for demonstrating autonomous proximity operations between two spacecraft, verifying deployment of the Planetary Society’s LightSail-B, and conducting remote sensing of Earth. Prox-1—THESIS was selected as the finalist for the competed University Nanosatellite Program-7 and was awarded a launch on the maiden commercial SpaceX Falcon Heavy. THESIS captures 8-12 µm IR images with 100 mm optics and RGB color images with 25 mm optics. The instrument computer was capable of instrument commanding, automatic data processing, image storage, and telemetry recording. The completed THESIS has a mass of 2.04 kg, a combined volume of 3U, and uses 7W of power. THESIS was designed, fabricated, integrated, and tested in ASU’s 100K clean lab.

Explosive mafic (basaltic) volcanism is not easily explained by current eruption models, which predict low energy eruptions from low viscosity magma due to decoupling of volatiles (gases). Sunset Crater volcano provides an example of an alkali basalt magma that produced a highly explosive sub-Plinian eruption. I investigate the possible role of magmatic volatiles in the Sunset Crater eruption through study of natural samples of trapped volatiles (melt inclusions) and experiments on mixed-volatile (H2O-CO2) solubility in alkali-rich mafic magmas.

I conducted volatile-saturated experiments in six mafic magma compositions at pressures between 400 MPa and 600 MPa to investigate the influence of alkali elements (sodium and potassium) on volatile solubility. The experiments show that existing volatile solubility models do not accurately describe CO2 solubility at mid-crustal depths. I calculate thermodynamic fits for solubility in each composition and calibrate a general thermodynamic model for application to other mafic magmas. The model shows that the relative percent abundances of sodium, calcium, and potassium have the greatest influence on CO2 solubility in mafic magmas.

I analyzed olivine-hosted melt inclusions (MIs) from Sunset Crater to investigate pre-eruptive volatiles. I compared the early fissure activity to the sub-Plinian eruptive phases. The MIs are similar in major element and volatile composition suggesting a relatively homogeneous magma. The H2O content is relatively low (~1.2 wt%), whereas the dissolved CO2 content is high (~2300 ppm). I explored rehomogenization and Raman spectroscopy to quantify CO2 abundance in MI vapor bubbles. Calculations of post-entrapment bubble growth suggest that some MI bubbles contain excess CO2. This implies that the magma was volatile-saturated and MIs trapped exsolved vapor during their formation. The total volatile contents of MIs, including bubble contents but excluding excess vapor, indicate pre-eruptive magma storage from 10 km to 18 km depth.

The high CO2 abundance found in Sunset Crater MIs allowed the magma to reach volatile-saturation at mid-crustal depths and generate overpressure, driving rapid ascent to produce the explosive eruption. The similarities in MIs and volatiles between the fissure eruption and the sub-Plinian phases indicate that shallow-level processes also likely influenced the final eruptive behavior.

Silicic volcanoes produce many styles of activity over a range of timescales. Eruptions vary from slow effusion of viscous lava over many years to violent explosions lasting several hours. Hazards from these eruptions can be far-reaching and persistent, and are compounded by the dense populations often surrounding active volcanoes. I apply and develop satellite and ground-based remote sensing techniques to document eruptions at Merapi and Sinabung Volcanoes in Indonesia. I use numerical models of volcanic activity in combination with my observational data to describe the processes driving different eruption styles, including lava dome growth and collapse, lava flow emplacement, and transitions between effusive and explosive activity.

Both effusive and explosive eruptions have occurred recently at Merapi volcano. I use satellite thermal images to identify variations during the 2006 effusive eruption and a numerical model of magma ascent to explain the mechanisms that controlled those variations. I show that a nearby tectonic earthquake may have triggered the peak phase of the eruption by increasing the overpressure and bubble content of the magma and that the frequency of pyroclastic flows is correlated with eruption rate. In 2010, Merapi erupted explosively but also shifted between rapid dome-building and explosive phases. I explain these variations by the heterogeneous addition of CO2 to the melt from bedrock under conditions favorable to transitions between effusive and explosive styles.

At Sinabung, I use photogrammetry and satellite images to describe the emplacement of a viscous lava flow. I calculate the flow volume (0.1 km3) and average effusion rate (4.4 m3 s-1) and identify active regions of collapse and advance. Advance rate was controlled by the effusion rate and the flow’s yield strength. Pyroclastic flow activity was initially correlated to the decreasing flow advance rate, but was later affected by the underlying topography as the flow inflated and collapsed near the vent, leading to renewed pyroclastic flow activity.

This work describes previously poorly understood mechanisms of silicic lava emplacement, including multiple causes of pyroclastic flows, and improves the understanding, monitoring capability, and hazard assessment of silicic volcanic eruptions.

Shallow earthquakes in the upper part of the overriding plate of subduction zones can be devastating due to their proximity to population centers despite the smaller rupture extents than commonly occur on subduction megathrusts that produce the largest earthquakes. Damaging effects can be greater in volcanic arcs like Java because ground shaking is amplified by surficial deposits of uncompacted volcaniclastic sediments. Identifying the upper-plate structures and their potential hazards is key for minimizing the dangers they pose. In particular, the knowledge of the regional stress field and deformation pattern in this region will help us to better understand how subduction and collision affects deformation in this part of the overriding plate. The majority of the upper plate deformation studies have been focused on the deformation in the main thrusts of the fore-arc region. Study of deformation within volcanic arc is limited despite the associated earthquake hazards. In this study, I use maps of active upper-plate structures, earthquake moment tensor data and stress orientation deduced from volcano morphology analysis to characterize the strain field of Java arc. In addition, I use sandbox analog modeling to evaluate the mechanical factors that may be important in controlling deformation. My field- and remotely-based mapping of active faults and folds, supplemented by results from my paleoseismic studies and physical models of the system, suggest that Java’s deformation is distributed over broad areas along small-scale structures. Java is segmented into three main zones based on their distinctive structural patterns and stress orientation. East Java is characterized by NW-SE normal and strike-slip faults, Central Java has E-W folds and thrust faults, and NE-SW strike-slip faults dominate West Java. The sandbox analog models indicate that the strain in response to collision is partitioned into thrusting and strike-slip faulting, with the dominance of margin-normal thrust faulting. My models test the effects of convergence obliquity, geometry, preexisting weaknesses, asperities, and lateral strength contrast. The result suggest that slight variations in convergence obliquity do not affect the deformation pattern significantly, while the margin shape, lateral strength contrast, and perturbation of deformation from asperities each have a greater impact on deformation.

Olympus Mons is the largest volcano on Mars. Previous studies have focused on large scale features on Olympus Mons, such as the basal escarpment, summit caldera complex and aureole deposits. My objective was to identify and characterize previously unrecognized and unmapped small scale features to understand the volcanotectonic evolution of this enormous volcano. For this study I investigated flank vents and arcuate graben. Flank vents are a common feature on composite volcanoes on Earth. They provide information on the volatile content of magmas, the propagation of magma in the subsurface and the tectonic stresses acting on the volcano. Graben are found at a variety of scales in close proximity to Martian volcanoes. They can indicate flexure of the lithosphere in response to the load of the volcano or gravitation spreading of the edifice. Using Context Camera (CTX), High Resolution Imaging Science Experiment (HiRISE), Thermal Emission Imaging System (THEMIS), High Resolution Stereo Camera Digital Terrain Model (HRSC DTM) and Mars Orbiter Laser Altimeter (MOLA) data, I have identified and characterized the morphology and distribution of 60 flank vents and 84 arcuate graben on Olympus Mons. Based on the observed vent morphologies, I conclude that effusive eruptions have dominated on Olympus Mons in the Late Amazonian, with flank vents playing a limited role. The spatial distribution of flank vents suggests shallow source depths and radial dike propagation. Arcuate graben, not previously observed in lower resolution datasets, occur on the lower flanks of Olympus Mons and indicate a recent extensional state of stress. Based on spatial and superposition relationships, I have constructed a developmental sequence for the construction of Olympus Mons: 1) Construction of the shield via effusive lava flows.; 2) Formation of the near summit thrust faults (flank terraces); 3) Flank failure leading to scarp formation and aureole deposition; 4) Late Amazonian effusive resurfacing and formation of flank vents; 5) Subsidence of the caldera, waning volcanism and graben formation. This volcanotectonic evolution closely resembles that proposed on Ascraeus Mons. Extensional tectonism may continue to affect the lower flanks of Olympus Mons today.