Silicate Eruption Dynamics on Io: Comparison with Observations
Explore the source record for details and available documents.
Io volcanos silicate magma↗
SEARCH · Search NASA
Results for “Eruption dynamics”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Lava Eruption and Emplacement: Using Clues from Hawaii and Iceland to Probe the Lunar Past
Investigating recent eruptions on Earth is crucial to improving understanding of relationships between eruption dynamics and final lava flow morphologies. In this study, we investigated eruptions in Holuhraun, Iceland, and Kilauea, Hawaii to gain insight into the lava dynamics near the source vent, the initiation of lava channels, and the origin of down-channel features. Insights are applied to Rima Bode on the lunar nearside to deduce the sequence of events that formed this lunar sinuous rille system.These insights are crucial to correctly interpreting whether the volcanic features associated with Rima Bode directly relate to eruption conditions at the vent and, thus, can help us understand those eruption dynamics, or, alternatively, whether the features formed as a result of more localized influences on lava flow dynamics. For example, if the lava channel developed early in the eruption and was linked to pulses in vent activity, its morphology can be analyzed to interpret the flux and duration of the eruption. Conversely, if the lava channel initiated late in the eruption as the result of a catastrophic breaching of lava that had previously pooled within the vent [e.g., 1], then the final channel morphology will not indicate eruption dynamics but rather local dynamics associated with that breach event. Distinguishing between these two scenarios is crucial for correctly interpreting the intensity and duration of volcanic history on the Moon.
sinuous rille↗
Lava Eruption and Emplacement: Using Clues from Hawaii and Iceland to Probe the Lunar Past
Investigating recent eruptions on Earth is crucial to improving understanding of relationships between eruption dynamics and final lava flow morphologies. In this study, we investigated eruptions in Holuhraun, Iceland, and Kilauea, Hawaii to gain insight into the lava dynamics near the source vent, the initiation of lava channels, and the origin of down-channel features. Insights are applied to Rima Bode on the lunar nearside to deduce the sequence of events that formed this lunar sinuous rille system. These insights are crucial to correctly interpreting whether the volcanic features associated with Rima Bode directly relate to eruption conditions at the vent and, thus, can help us understand those eruption dynamics, or, alternatively, whether the features formed as a result of more localized influences on lava flow dynamics. For example, if the lava channel developed early in the eruption and was linked to pulses in vent activity, its morphology can be analyzed to interpret the flux and duration of the eruption. Conversely, if the lava channel initiated late in the eruption as the result of a catastrophic breaching of lava that had previously pooled within the vent [e.g., 1], then the final channel morphology will not indicate eruption dynamics but rather local dynamics associated with that breach event. Distinguishing between these two scenarios is crucial for correctly interpreting the intensity and duration of volcanic history on the Moon.
Needham, Debra Hurwitz↗
Three-dimensional Dynamics of Strongly Twisted Magnetar Magnetospheres: Kinking Flux Tubes and Global Eruptions
Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement of ψ crit ≳ π . In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲10 43 erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks.
79 ASTRONOMY AND ASTROPHYSICS↗
Effects of Optimisation Parameters on Data-Driven Magnetofrictional Modelling of Active Regions
Context . The solar magnetic field plays an essential role in the formation, evolution, and dynamics of large-scale eruptive structures in the corona. The estimation of the coronal magnetic field, the ultimate driver of space weather, particularly in the ‘low’ and ‘middle’ corona, is presently limited due to practical difficulties. Data-driven time-dependent magnetofrictional modelling (TMFM) of active region magnetic fields has been proven to be a useful tool to study the corona. The input to the model is the photospheric electric field that is inverted from a time series of the photospheric magnetic field. Constraining the complete electric field, that is, including the non-inductive component, is critical for capturing the eruption dynamics. We present a detailed study of the effects of optimisation of the non-inductive electric field on the TMFM of AR 12473. Aims . We aim to study the effects of varying the non-inductive electric field on the data-driven coronal simulations, for two alternative parametrisations. By varying parameters controlling the strength of the non-inductive electric field, we wish to explore the changes in flux rope formation and their early evolution and other parameters, for instance, axial flux and magnetic field magnitude. Methods . We used the high temporal and spatial resolution cadence vector magnetograms from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). The non-inductive electric field component in the photosphere is critical for energising and introducing twist to the coronal magnetic field, thereby allowing unstable configurations to be formed. We estimated this component using an approach based on optimising the injection of magnetic energy. Results . Our data show that flux ropes are formed in all of the simulations except for those with the lower values of these optimised parameters. However, the flux rope formation, evolution and eruption time varies depending on the values of the optimisation parameters. The flux rope is formed and has overall similar evolution and properties with a large range of non-inductive electric fields needed to determine the non-inductive electric field component that is critical for energising and introducing twist to the coronal magnetic field. Conclusions . This study shows that irrespective of non-inductive electric field values, flux ropes are formed and erupted, which indicates that data-driven TMFM can be used to estimate flux rope properties early in their evolution without needing to employ a lengthy optimisation process.
A. Kumari↗
Investigating Si Bonding Environments in Mercurian Melts
Mercury’s low oxygen fugacity results in lavas with extraordinarily high sulfur contents (> 1 wt%). The high S lowers melt viscosity, but why this occurs is unclear (Mouser et al., 2021). The detection of Si-S bonds (using NMR) in such highly reduced basaltic glasses, suggest that the decreased viscosity is due to changes in the Si bonding environment (Pommier et al., 2023). To investigate further, we synthesized a series of Mercurian glasses at 1 GPa, with S contents ranging from 0-6 wt% and oxygen fugacities from IW-4 to IW-6. The experimental glasses were analyzed by Si K-edge X-ray Absorption Near Edge Structure (XANES) spectroscopy and 29Si Nuclear Magnetic Resonance (NMR). Ongoing analyses will include Raman and Fourier Transform Infrared (FTIR). The 29Si NMR analyses suggest the presence of Si-S bonds in the melts, in agreement with Pommier et al. (2023). The Si K-edge XANES spectra show a decrease in K-edge position of 0.4 to 0.6 eV from terrestrial basalts (at ~IW+3) to the experimental melts at IW-4 to IW-6. This shift suggests a moderate (<20 mole %) amount of Si has shifted to lower valence states (Boujibar et al, 2019). The melts contain 0-6 wt% S, so if the K-edge shifts are the results of Si-S bonding, it implies that the majority of the S in the melts is bonded to Si. However, S K-edge XANES analyses indicate that much of the S is also bonded with Mg and Ca at these low oxygen fugacities (Anzures et al, 2020). The implication is that non-bridging Si-S-(Mg,Ca) bonds are prevalent in the reduced melts, consistent with the decrease in viscosity of the melts being due to melt depolymerization. Overall, our results suggest that dissolution of S into the melts at low fO2 is accomplished through bonding with Si, and is accompanied by significant structural changes that lower melt viscosity. Lower melt viscosities will affect many aspects of Mercury’s evolution, from magma ocean dynamics, to melt migration in the mantle, to eruption dynamics and degassing.
Emily L Fischer↗
Acoustic noise from volcanoes - Theory and experiment
The paper discusses some theoretical aspects of acoustic investigation of volcanoes and describes a field experiment involving the recording, analysis, and interpretation of acoustic radiation from energetic fumaroles at Volcan Acatenango, Guatemala, during mid-January 1973. Particular attention is given to deriving information about the flow velocity of the erupting medium from acoustics as a means to study eruption dynamics. Theoretical considerations suggest that acoustic power radiated during gaseous volcanic eruptions may be related to gas exit velocity according to appropriate power laws. Eruption acoustics proves useful as a means of quantitative monitoring of volcanic activity.
Woulff, G.↗
Volcano infrasound: progress and future directions
Over the past two decades (2000–2020), volcano infrasound (acoustic waves with frequencies less than 20 Hz propagating in the atmosphere) has evolved from an area of academic research to a useful monitoring tool. As a result, infrasound is routinely used by volcano observatories around the world to detect, locate, and characterize volcanic activity. It is particularly useful in confirming subaerial activity and monitoring remote eruptions, and it has shown promise in forecasting paroxysmal activity at open-vent systems. Fundamental research on volcano infrasound is providing substantial new insights on eruption dynamics and volcanic processes and will continue to do so over the next decade. The increased availability of infrasound sensors will expand observations of varied eruption styles, and the associated increase in data volume will make machine learning workflows more feasible. More sophisticated modeling will be applied to examine infrasound source and propagation effects from local to global distances, leading to improved infrasound-derived estimates of eruption properties. Future work will use infrasound to detect, locate, and characterize moving flows, such as pyroclastic density currents, lahars, rockfalls, lava flows, and avalanches. Infrasound observations will be further integrated with other data streams, such as seismic, ground- and satellite-based thermal and visual imagery, geodetic, lightning, and gas data. The volcano infrasound community should continue efforts to make data and codes accessible and to improve diversity, equity, and inclusion in the field. In summary, the next decade of volcano infrasound research will continue to advance our understanding of complex volcano processes through increased data availability, sensor technologies, enhanced modeling capabilities, and novel data analysis methods that will improve hazard detection and mitigation.
58 GEOSCIENCES↗
Controls on the frequency content of near-source infrasound at open-vent volcanoes: a case study from Volcán Villarrica, Chile
Infrasound signals from open-vent volcanoes can contain information related to that volcano’s eruption dynamics and future activity. Studying a specific volcano’s acoustic signature during periods of relative quiescence is thus important for understanding potential warning signs of impending eruptions. Volcán Villarrica, located in southern Chile, has a long-lived active lava lake that produces continuous infrasound with spectral peaks centered around 1 Hz that vary by + / - ~ 0.2 Hz over day-to-week-long time scales. The infrasound frequency content has been shown to illuminate key volcanic properties such as eruption style and crater shape. Leading up to Villarrica’s most recent paroxysm in 2015, for instance, infrasound spectral changes coincided with a rise in Villarrica’s lava lake level. Quantifying and understanding the regular fluctuation in recorded infrasonic frequencies from Villarrica and other open-vent volcanoes is thus imperative during periods of relative calm. Our experiment entailed a week-long period of crater rim infrasound observations associated with stable, open-vent activity and revealed two independent source processes: spatter bursting events (relatively broadband infrasound) and lava lake-induced tremor (sustained signals peaked at 1 Hz). A comparison of these acoustic signals with results from 3-D finite-difference time-domain wave propagation models shows that while sound speed variability can influence Villarrica’s spectrum, it cannot explain the full extent of the observed frequency excursions. We conclude that source spectrum variability is primarily responsible for the frequency excursions. In conclusion, our work highlights the utility of data collected from a dense distribution of twenty infrasound sensors operating at the summit coupled with numerical modeling of sound radiation.
58 GEOSCIENCES↗
Twenty Years of NASA-EOS Multi-Sensor Satellite Observations at Kīlauea Volcano (2000–2019)
Recent eruptions of Kīlauea volcano, on the Island of Hawai'i, represent an ideal test location for studying volcanic plumes, due to its remote location and minimal anthropogenic pollution. Within this work we exploit the over 20-year data record from NASA EOS satellites to investigate the degree to which space-borne observations can detect shifts in known eruption dynamics at Kīlauea. Through the combined analysis of remotely sensed plume heights, plume particle microphysical properties, lava flow thermal tracking and high resolution SO2 mapping, we provide interpretations of the satellite signals in a volcanic context and link them to ground-based eruption reports for validation. We establish common patterns of plume dispersion and particle evolution and identify when significant shifts in activity occur. We also determine where ambient conditions influence the output retrieved from the space-borne sensors. By comparing the inferences derived from remote-sensing with the extensive record of suborbital observations at Kīlauea, we assess the strengths and limitations of the satellite-based volcanic assessment techniques. The work presented here highlights the capabilities of these techniques, allowing us to more confidently interpret volcanic activity observed from space globally, particularly in cases where ground observations are limited or entirely lacking.
volcano remote sensing↗
Filament Eruptions, Jets, and Space Weather
Previously, from chromospheric H alpha and coronal X‐ray movies of the Sun's polar coronal holes, it was found that nearly all coronal jets (greater than 90%) are one or the other of two roughly equally common different kinds, different in how they erupt: standard jets and blowout jets (Yamauchi et al 2004, Apl, 605, 5ll: Moore et all 2010, Apj, 720, 757). Here, from inspection of SDO/AIA He II 304 A movies of 54 polar x-ray jets observed in Hinode/XRT movies, we report, as Moore et al (2010) anticipated, that (1) most standard x-ray jets (greater than 80%) show no ejected plasma that is cool enough (T is less than or approximately 10(exp 5K) to be seen in the He II 304 A movies; (2) nearly all blownout X-ray jets (greater than 90%) show obvious ejection of such cool plasma; (3) whereas when cool plasma is ejected in standard X-ray jets, it shows no lateral expansion, the cool plasma ejected in blowout X-ray jets shows strong lateral expansion; and (4) in many blowout X-ray jets, the cool plasma ejection displays the erupting-magnetic-rope form of clasic filament eruptions and is thereby seen to be a miniature filament eruption. The XRT movies also showed most blowout X-ray jets to be larger and brighter, and hence to apparently have more energy, than most standard X-ray jets. These observations (1) confirm the dichotomy of coronal jets, (2) agree with the Shibata model for standard jets, and (3) support the conclusion of Moore et al (2010) that in blowout jets the magnetic-arch base of the jet erupts in the manner of the much larger magnetic arcades in which the core field, the field rooted along the arcade's polarity inversion line, is sheared and twisted (sigmoid), often carries a cool‐plasma filament, and erupts to blowout the arcade, producing a CME. From Hinode/SOT Ca II movies of the polar limb, Sterling et al (2010, ApJ, 714, L1) found that chromospheric Type‐II spicules show a dichotomy of eruption dynamics similar to that found here for the cool‐plasma component of coronal X‐ray jets. This favors the idea that Type‐II spicules are miniature counterparts of coronal X‐ray jets. In Moore et al (2011, ApJ, 731, L18), we pointed out that if Type‐II spicules are magnetic eruptions that work like coronal X‐ray jets, they carry an area‐averaged mechanical energy flux of approximately 7x10)(exp 5) erg cm(exp -2) s(exp-1) into the corona in the form of MHD waves and jet outflow, enough to power the heating of the global corona and solar wind. On this basis, from our observations of mini‐filament eruptions in blowout X‐ray jets, we infer that magnetic explosions of the type that have erupting filaments in them are the main engines of both (1) the steady solar wind and (2) the CMEs that produce the most severe space weather by blasting out through the corona and solar wind, making solar energetic particle storms, and bashing the Earth's magnetosphere. We conclude that in focusing on prominences and filament eruptions, Einar had his eye on the main bet for understanding what powers all space weather, both the extreme and the normal.
Moore, Ronald↗
Solar Activity Modeling: From Subgranular Dynamical Scales to the Solar Cycles
The dynamical effects of solar magnetoconvection span a wide range spatial and temporal scales that extend from the interior to the corona and from fast turbulent motions to global magnetic activity. To study the solar activity on short temporal scales (from minutes to hours), we use 3D radiative MHD simulations that allow us to investigate complex turbulent interactions that drive various phenomena, such as plasma eruptions, spontaneous formation of magnetic structures, funnel-like structures and magnetic loops in the corona, and others. In particular, we focus on multi-scale processes of energy exchange across layers of the solar interior and atmosphere, which contribute to coronal heating and eruptive dynamics. For modeling global-scale activity, we use a data assimilation approach that has demonstrated great potential for building reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to the Parker-Kleeorin-Ruzmakin dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of observational limitations on prediction accuracy, and we present the EnKF predictions of the upcoming Solar Cycle (25) based on both the sunspot number series and observed magnetic fields and discuss the uncertainties and potential of the data assimilation approach for modeling and forecasting solar activity.
Kitiashvili, I. N.↗
Solar Activity Modeling: From Subgranular Dynamical Scales to the Solar Cycles
Dynamical effects of solar magnetoconvection span a wide range spatial and temporal scales that extends from the interior to the corona and from fast turbulent motions to the global-Sun magnetic activity. To study the solar activity on short temporal scales (from minutes to hours), we use 3D radiative MHD simulations that allow us to investigate complex turbulent interactions that drive various phenomena, such as plasma eruptions, spontaneous formation of magnetic structures, funnel-like structures and magnetic loops in the corona, and others. In particular, we focus on multi-scale processes of energy exchange across the different layers, which contribute to the corona heating and eruptive dynamics, as well as interlinks between different layers of the solar interior and atmosphere. For modeling the global-scale activity we use the data assimilation approach that has demonstrated great potential for building reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to the Parker-Kleeorin-Ruzmakin dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of observational limitation on the prediction accuracy. We present the EnKF predictions of the upcoming Solar Cycle 25 based on both the sunspot number series and observed magnetic fields, and discuss the uncertainties and potential of the data assimilation approach for modeling and forecasting the solar activity.
Kitiashvili, I. N.↗
The ultralow viscosity of volatile-rich kimberlite magma: Implications for the water content of primitive kimberlite melts
The eruption of deeply sourced kimberlite magma offers the fastest route to bring deep-seated volatiles back to the Earth’s surface. However, the viscosity of kimberlite magma, a factor governing its migration and eruption dynamics within Earth, remains poorly constrained. We conducted synchrotron in situ falling sphere viscometry experiments to examine kimberlite magma with different volatile contents (0 to 5 wt % H 2 O and 2 to 8 wt % CO 2 ) under high pressure-temperature conditions. The results reveal that the viscosity of volatile-rich kimberlite magma is ~1 to 2 orders lower than that of mid-ocean ridge basalt (MORB) and comparable to the ultramobile pure carbonate melt. Using the measured viscosity values, we simulated the ascent and eruption process of kimberlite magma. We found that a minimum content of ~0.5 wt % water in the primitive magma is necessary to allow the ultrafast eruption process of kimberlite, thereby enabling the preservation of diamonds and high-pressure mineral inclusions transported by the magma.
58 GEOSCIENCES↗
The MURAVES Experiment: A Study of the Vesuvius Great Cone with Muon Radiography
The MURAVES experiment aims at the muographic imaging of the internal structure of the summit of Mt. Vesuvius, exploiting muons produced by cosmic rays. Though presently quiescent, the volcano carries a dramatic hazard in its highly populated surroundings. The challenging measurement of the rock density distribution in its summit by muography, in conjunction with data from other geophysical techniques, can help the modeling of possible eruptive dynamics. The MURAVES apparatus consists of an array of three independent and identical muon trackers, with a total sensitive area of 3 square meters. In each tracker, a sequence of 4 XY tracking planes made of plastic scintillators is complemented by a 60 cm thick lead wall inserted between the two downstream planes to improve rejection of background from low-energy muons. The apparatus is currently acquiring data. Preliminary results from the analysis of the first data sample are presented.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Lunar magma transport phenomena
An outline of magma transport theory relevant to the evolution of a possible Lunar Magma Ocean and the origin and transport history of the later phase of mare basaltic volcanism is presented. A simple model is proposed to evaluate the extent of fractionation as magma traverses the cold lunar lithosphere. If Apollo green glasses are primitive and have not undergone significant fractionation en route to the surface, then mean ascent rates of 10 m/s and cracks of widths greater than 40 m are indicated. Lunar tephra and vesiculated basalts suggest that a volatile component plays a role in eruption dynamics. The predominant vapor species appear to be CO CO2, and COS. Near the lunar surface, the vapor fraction expands enormously and vapor internal energy is converted to mixture kinetic energy with the concomitant high-speed ejection of vapor and pyroclasts to form lunary fire fountain deposits such as the Apollo 17 orange and black glasses and Apollo 15 green glass.
Spera, Frank J.↗
Building Reliable Forecasts of Solar Activity
Solar ionizing radiation critically depends on the level of the Sun’s magnetic activity. For robust physics-based forecasts, we employ the procedure of data assimilation, which combines theoretical modeling and observational data such that uncertainties in both the model and the observations are taken into account. Currently we are working in two major directions: 1) development of a new long-term forecast procedure on time-scales of the 11-year solar cycle, using a 2-dimensional mean-field dynamo model and synoptic magnetograms; 2) development of 3-dimensional radiative MHD (Magnetohydrodynamic) simulations to investigate the origin and precursors of local manifestations of magnetic activity, such as the formation of magnetic structures and eruptive dynamics.
Kitiashvili, Irina↗
VESIcal: An Open-source Thermodynamic Model Engine for Mixed Volatile (H2O-CO2) Solubility in Silicate Melts
Modeling the solubility of volatiles in silicate melts is fundamental to the interpretation of volcanic systems and has implications for magma dynamics, eruption style, and material transport between the mantle, crust, and atmosphere. Recent advancements in computational capabilities and access to computing tools has outpaced the functionality and extensibility of previously available modeling platforms. Here we present VESIcal (Volatile Equilibria and Saturation Index calculator), the first comprehensive modeling tool for H2O, CO2, and mixed (H2O-CO2) solubility in silicate melts that: a) allows users access to seven popular models, with easy inter-comparison between models; b) provides universal functionality for all models (e.g., functions for calculating saturation pressures, degassing paths, etc.); c) can process large datasets (1,000’s of samples) automatically; d) can output computed data into an Excel spreadsheet or CSV file for post-modeling analysis; e) integrates plotting capabilities directly within the tool; and f) provides all of this within the framework of a python library, making the tool extensible by the user and allowing any of the model functions to be incorporated into any other code capable of calling python.Here we will provide a demonstration of VESIcal and its capabilities with applications to various volcanic processes affected by volatiles. VESIcal represents the first tool capable of directly comparing multiple solubility models and equations of state. We find that commonly used models predict surprisingly different volatile solubilities, particularly for pure CO2 or mixed CO2-H2O fluids. Even for melt compositions that are well represented in the calibration datasets of multiple models (e.g., MORBs), calculated solubilities for pure CO2 and pure H2O can deviate from one another by factors of >2 leading to 2x deviations in calculated saturation pressures (e.g., 5 to 10 kbar). The solubility of CO2 predicted by different rhyolitic models also differs substantially, overwhelming other sources of uncertainty such as analytical errors on measurements of volatile contents or uncertainties in crustal density profiles. This highlights the importance of model choice when drawing geological conclusions based on volatiles in magmas.
Kayla Iacovino↗