Thesis (Ph. D.)--University of Oregon, 2006. / Typescript. Includes vita and abstract. Includes bibliographical references (leaves 162-173). Also available for download via the World Wide Web; free to University of Oregon users.

The thermal evolution of pyroclastic density currents (PDCs) is the result of entrainment of ambient air, particle concentration, and initial eruptive temperature, which all impact PDC dynamics and their hazards, such as runout distance. The associated hazards and opaqueness of PDCs make it impossible for in-situ entrainment efficiencies or concentration measurements that would provide critical information on the thermal evolution and physical processes of PDCs. The thermal evolution of explosive eruptive events such as volcanic plumes and pyroclastic density currents (PDCs) is reflected in the textures of the material they deposit. A multiscale model is developed to evaluate how the rinds of breadcrust bombs can be used as a unique thermometer to examine the thermal evolution of PDCs. The multiscale, integrated model examines how bubble growth, pyroclast cooling, and dynamics of PDC and projectile pyroclasts form unique pyroclast morphology. Rind development is examined as a function of transport regime (PDC and projectile), transport properties (initial current temperature and current density), and pyroclast properties (initial water content and radius). The model reveals that: 1) rinds of projectile pyroclasts are in general thicker and less vesicular than those of PDC pyroclasts; 2) as the initial current temperature decreases due to initial air entrainment, the rinds on PDC pyroclasts progressively increase in thickness; and 3) rind thickness increases with decreasing water concentration and decreasing clast radius. Therefore, the modeled pyroclast’s morphology is dependent not only on initial water concentration but also on the cooling rate, which is determined by the transport regime. The developed secondary thermal proxy is then applied to the 2006 PDCs from the Tungurahua eruption to constrain the entrainment efficiency and thermal evolution of PDCs. A three-dimensional multiphase Eulerian-Eulerian-Lagrangian (EEL) model is coupled to topography and field data such as paleomagnetic data and rind thicknesses of collected pyroclasts to study the entrainment efficiency and thus the thermal history of PDCs at Tungurahua volcano, Ecuador. The modeled results that are constrained with observations and thermal proxies demonstrate that 1) efficient entrainment of air to the upper portion of the current allows for rapid cooling, 2) the channelized pyroclastic density currents may have developed a stable bed load region that was inefficient at cooling and 3) the PDCs had temperatures of 600-800K in the bed load region but the upper portion of the currents cooled down to ambient temperatures. The results have shown that PDCs can be heterogeneous in particle concentration, temperature, and dynamics and match observations of PDCs down a volcano and the thermal proxies. Lastly, the entrainment efficiencies of PDCs increases with increasing PDC temperature and entrainment varies spatially and temporally. Therefore, the assumption of a well-mixed current with a single entrainment coefficient cannot fully solve the thermal evolution and dynamics of the PDC.

Mafic magma is the most common magma erupted at the surface of the earth. It is generated from partial melting of the mantle, which has been subdivided into end-members based on unique geochemical signatures. One reason these end members, or heterogeneities, exist is subduction of lithospheric plates back into the mantle. The amount of elements, such as Cl and K, removed during subduction and recycled into the deep mantle, is poorly constrained. Additionally, the amount of volatiles, such as Cl, that are recycled into the deep mantle will strongly affect the behavior of the system. I have looked at Cl and K in HIMU source melts to see how it varies. Cl/Nb and K/Nb suggest that elevated Cl/K ratios are the result of depletion of K rather than increased Cl recycled into the deep mantle. After the mantle has partially melted and mafic melt has migrated to the surface, it usually erupts effusively or with low explosivity because of its low viscosity, but it is possible for larger eruptions to occur. These larger, Plinian eruptions, are not well understood in mafic systems. It is generally thought that basalt has a viscosity that is too low to allow for such an eruption to occur. Plinian eruptions require fragmentation to occur, which means the melt must undergo brittle failure. This may occur if the melt ascends rapidly enough to allow pressure to build in bubbles without the bubbles expanding. To test this, I have done decompression experiments to try to bracket the ascent rate for two Plinian eruptions. One eruption has a fast ascent, faster than those seen in more silicic melts, whereas the other eruption is unable to be reproduced in the lab, however it began with a increased viscosity in the partly crystallized magma. After fragmentation and eruption, it is generally thought that tephra do not continue to crystallize. We have found that crystallinity increases from rim to core in two basaltic pumice. Textural data along with a cooling model has allowed us to estimate growth rates in a natural system, which are similar to experimental data. / text

Includes appendix on CD-ROM (p. 157). / Thesis (Ph. D.)--University of Hawaii at Manoa, 2004. / Includes bibliographical references (leaves 158-175). / Also available by subscription via World Wide Web / xiv, 175 leaves, bound ill. (some col., one folded), maps (some col.) 29 cm. +

Thesis (Ph. D.)--Oregon State University, 2006. / Includes map in pocket. Printout. Includes bibliographical references (leaves 170-185). Also available on the World Wide Web.

Magnetic flux emergence is a fundamental process in the Sun, during which magnetic fields emerge from the solar interior to the surface, to build up active regions and give onset to spectacular dynamic phenomena, such as eruptions and jets. In this thesis, we performed 3D, resistive MHD simulations to study the emergence and the associated magnetic activity of a quadrupolar region in the Sun. Our aim behind the setup of this initial condition (i.e. a quadrupolar region) was to study a magnetic field configuration, which has not been studied in detail before, although it has been repeatedly observed in the Sun and it has been shown that it can host intense magnetic activity (e.g. in the form of jets, flares and eruptions). The results of our experiments showed that the internal dynamics of such regions leads to the onset of eruptions in the form of twisted magnetic flux tubes (flux ropes). These eruptions are recurrent but they cannot escape the outermost field of the emerging flux (envelope field). They remain confined within the envelope field, as the downward tension of the outermost field lines overwhelms the upward Lorentz force of the erupting field. When we add an ambient magnetic field in the solar atmosphere, external reconnection between the emerging and the ambient field triggers the emission of (standard) reconnection jets. The external reconnection also releases the tension of the ambient field lines and, thus, the eruptions move in an ejective way towards the outer space. Namely, the confined eruptions become ejective eruptions, which escape from the numerical domain. These ejective eruptions drive a newly observed class of jets, the so called "blowout" jets. Our experiments reproduce some of the main observed characteristics of the "blowout" jets. We showed that "blowout" jets emit hot and cool plasma into the outer solar atmosphere simultaneously, and they undergo untwisting motion due to the relaxation of twist during their ejection. We found that the untwisting motion of the "blowout" jets is associated with the propagation of torsional Alfvén waves. Finally, we performed a parametric study to explore the effect of the ambient field strength on the onset and dynamics of the eruptive events. We found that one of the main effects is that the stronger ambient field suppresses the vertical expansion of the magnetic envelope of the quadrupolar region due to the higher magnetic pressure above it. This result has an effect on the emission of jets, which are emitted due to reconnection between the two fields. When the ambient field is relatively weak, it is pushed away from the strong emerging field and reconnection between them is not so persistent. On the other hand, when the ambient field is relatively strong, we find that more jets are ejected due to more efficient and more frequent reconnection between the two flux systems. As a consequence, we find that more mass and flux is being transferred into the solar corona by the reconnection jets. Also, we find that there are more eruptions when the ambient field is stronger. The study of the total energy flux carried by the jets showed that it is sufficient to provide the energy required to accelerate the high speed solar wind. This indicates that the "blowout" jets may play an important role in driving the solar wind.

Le 26 octobre 2010, le Merapi, volcan le plus peuplé et le plus actif de l'île de Java en Indonésie, est entré en éruption. Depuis plusieurs mois, l'activité du volcan avait crû de manière exponentielle et les scientifiques s'attendaient à une éruption exceptionnelle. Leurs pronostics furent confirmés puisqu'un blast parcourut une distance de 8 km, rasa en partie deux villages évacués, et tua mbah Marijan, le « gardien des clés du volcan » et son entourage, ce qui perturba l'équilibre culturel de la région. Dans les jours suivants, l'activité du volcan continua de s'intensifier et des évacuations massives furent organisées. L'éruption paroxysmale du Merapi eut lieu dans la nuit du 4 au 5 novembre ; alors que plusieurs dizaines de milliers de villageois étaient en train de fuir le volcan, de puissantes coulées pyroclastiques dévastèrent les flancs du volcan, tuant et détruisant tout sur leur passage sur une distance maximale de 17 km. Cette nuit-là, plus de 250 personnes furent tuées et des milliers d'autres blessées. Le secteur agricole subit de lourdes pertes puisque de nombreuses récoltes et un grand nombre de têtes de bétail furent perdus. L'activité touristique de la région connut une réaction en chaîne d'impacts, générée en grande partie par la fermeture de l'aéroport de Yogyakarta. La thèse propose de faire une étude transversale et systémique de ces impacts en analysant comment gouvernement et populations se remettent de cette éruption et s'apprêtent à faire face à une future crise causée par le Merapi. / On October 26th 2010, the Merapi volcano, the most populated and active volcano on the Java Island in Indonesia, began erupting. During several months, the volcano's activity grew exponentially and scientists were expecting an exceptional eruption. Their predictions got confirmed as a blast ran through a distance of 8 km, partially destructed two villages (that got evacuated) and killed mbah Marijan, the « Volcano's keys guardian » and his entourage. As a consequence, the cultural balance of the region got disrupted for a while. In the following days, the volcano's activity kept increasing and mass evacuations were organized. The paroxysmal eruption of Merapi occurred in the night of November 4th to 5th 2010, while tens of thousands of villagers were running away from the volcano. Powerful pyroclastic flows ran down its flanks, killing and destroying everything on their way over a maximum distance of 17 km. During that night, over 250 people were killed and thousands were injured. The agricultural sector suffered from heavy losses as many crops and livestock were lost. The tourism industry in the region experienced a chain reaction of impacts generated by the closure of Yogyakarta airport. The thesis offers to make a transversal and systemic study of these impacts by analyzing how the government and people recovered from the eruption and how they are to face a future crisis caused by the Merapi.

Le soleil est une étoile active, et les éruptions solaires sont une des manifestations de cette activité. Il est admis que l'énergie disponible pour les éruptions solaires a une origine magnétique, et est transmise au milieu lors de phénomènes de reconnexion magnétique dans la couronne. Une partie de cette énergie permet d'accélérer les particules du milieu (électrons et ions). Cependant, les détails concernant les conditions dans lesquelles les particules sont accélérées et se propagent des régions d'accélération aux sites d'interaction lors des éruptions solaires ne sont pas encore tous compris.Plusieurs modèles d'accélération de particules ont été développés dans le cadre de l'étude des éruptions solaires. Dans certains modèles, les particules sont accélérées par un champ électrique généré au niveau de couches de courants électriques, qui peuvent être fragmentées, et qui sont préférentiellement localisées au niveau de surfaces quasi-séparatrices. Afin d'étudier le lien entre l'accélération de particules et le champ électrique direct produit au niveau de couches de courants, nous avons recherché s'il y avait des corrélations entre les sites d'émission des particules énergétiques et les courants électriques mesurés au niveau de la photosphère. Les observations X (dur) représentent les diagnostics les plus directs des électrons énergétiques produits pendant les éruptions solaires (rayonnement de freinage des électrons dans l'atmosphère solaire) et nous avons donc utilisé les observations X du satellite RHESSI (Reuven Ramaty High Energy Solar Spectrometric Imager) afin de produire des images et des spectres du rayonnement X dur des électrons énergétiques. Afin de caractériser les courants électriques dans la région éruptive, nous avons utilisé les données spectropolarimétriques de l'instrument HMI (Helioseismic and Magnetic Imager) du satellite SDO (Solar Dynamic Observatory) et nous avons calculé les densités de courants verticales photosphériques à partir du champ magnétique vectoriel reconstruit. Une corrélation entre les émissions X coronales (dues aux particules énergétiques proches du site d'accélération) et les rubans de forte densité de courants photosphériques (traces des couches de courants coronales) a été mise en évidence pour les cinq éruptions de classe X étudiées. De plus, grâce à la cadenc / The Sun is an active star and one manifestation of its activity is the production of solar flares. It is currently admitted that solar flares are caused by the release of magnetic energy during the process of magnetic reconnection in the solar upper atmosphere, the solar corona. During these flares, a large fraction of the magnetic energy is transferred to the acceleration of particles (electrons and ions). However, the details of particle acceleration during flares are still not completely understood.Several scenarios and models have been developed to explain particle acceleration. In some of them, electric fields, produced at the location of current sheets, which can be fragmented or collapsing, and which are preferentially located on quasi-separatrix layers (QSLs), are accelerating particles. To investigate a possible link between energetic particles and direct electric fields produced at current sheet locations, we looked for a correlation between X-ray emission from energetic electrons and electric currents which can be measured at the photospheric level. We used the Reuven Ramaty High Energy Solar Spectrometric Imager (RHESSI) data to produce spectra and images of the X-ray emissions during GOES X-class flares, and spectropolarimetric data from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) to calculate the vertical current densities from the reconstructed 3D vector magnetic field. A correlation between the coronal X-ray emissions (tracing the energetic electrons near the acceleration site) and the strong current ribbons at the photospheric level (tracing the coronal current sheet) was found in the five studied X-class flares. Moreover, thanks to the 12-minute time cadence of SDO/HMI, we could study for the first time the time evolution of electric currents: in several flares, a change in the current intensity, occurring during the flare peak, was found to be spatially correlated with X-ray emission sites. These observations enlighten a common evolution of both electric currents and X-ray emissions during the

Volcanic eruptions pose hazards to life and insfrastructure, and contribute to the resurfacing of earth and other planetary bodies. Lavas and magmas are multi-phase suspensions of silicate melts (liquids), solid crystals, and vapor bubbles, and solidify into glass and rock upon cooling. The interactions between phases place important controls on the dynamics and timescales of magma and lava transport and emplacement. The purpose of this thesis is to explore the role of multiphase interactions in controlling eruption dynamics and inform conceptual and numerical models for hazard prediction. In Chapters 1 and 2, centimeter to meter scale analog experiments are used to explore the multi-phase rheological properties and flow behaviors of bubble- and particle-bearing suspensions. Optical imaging of dam-break experiments presented in Chapter 1 expand existing experimental parameter ranges for lava analogs to higher bubble concentrations than existing datasets (up to 82% by volume bubbles and 37% by volume particles). I develop a constitutive relationship for threephase relative viscosity, and demonstrate that at the low strain-rate conditions relevant to many natural lava flows, accounting for the rheological effect of bubbles can result in the prediction of slower runout speeds. Chapter 2 expands upon the work of Chapter 1 using different analog materials observed using nuclear magnetic resonance imaging (MRI) phase-contrast velocimetry (PCV) to measure velocity in the flow interior of three-phase dam-break experiments. I find that for high-aspect ratio particles (sesame seeds), phase segregation into shear bands readily occurs, even at low particle fraction (20%) and results in strain localization. I suggest that the presence of shear bands can lead to faster flow runout than predicted using assumptions of bulk rheology. Chapter 3 analyzes thermal infrared (IR) time-lapse photography and videography of Hawaiian to Strombolian explosive activity during the 2021 eruption of Cumbre Vieja volcano, La Palma, Canary Islands, Spain. Images are analyzed to find time series of apparent plume radius, velocity, and apparent volume flux of high-temperature gas and lava. I compare with other measures of eruptive activity, including remote observations of plume height, SO₂ flux, effusive flux, tremor, and events at the volcano edifice including edifice collapses and the opening of new vents. I find correlations between tremor and explosive flux, but no correlation with SO2 flux or effusive flux, which I interpret as evidence of bubble segregation, highlighting the role of phase segregation and temporal variability in material properties in natural systems. Finally, in Chapter 4, I develop a novel finite element model to explore the interaction between a viscous flow with a solidified crust, and the effect of these interactions on lava flow and lava dome emplacement. I develop a model that couples a temperature-dependent viscous interior with an elastic shell flowing into air, water, or dense atmospheres. The model expands upon existing numerical simulations used in volcanology to have direct applications to lava flows and domes on the sea floor, which accounts for a large portion of the volcanism on Earth, and volcanism on other planetary bodies. Additionally, the formation of levees or solidified flow fronts that fracture and lead to a restart of flow. These lava flow breakouts pose a significant hazard, but there are currently no volcanological community codes capable of using a physics-based approach to predict the timing or location of breakouts. The model in Chapter 4 is the first to allow for assessment of the likelihood of failure at the scale of a flow lobe. Chapter 4 describes the model formulation and verification, and validation against centimeter-scale molten basalt experiments. The dissertation as a whole integrates work using a variety of methods including analog experiments, observations of natural eruptions, and numerical simulations to contribute to our understanding of the effects of multi-phase interactions on volcanic eruptions.

学位プログラム名: 京都大学大学院思修館 / 京都大学 / 新制・課程博士 / 博士(総合学術) / 甲第23343号 / 総総博第16号 / 新制||総総||3(附属図書館) / 京都大学大学院総合生存学館総合生存学専攻 / (主査)教授 山敷 庸亮, 教授 寶 馨, 准教授 浅井 歩 / 学位規則第4条第1項該当 / Doctor of Philosophy / Kyoto University / DFAM