Search NASA⌕ Search

SEARCH · Search NASA

Results for “VOLCANIC THEORY”

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.

15 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Trapped bubbles keep pumice afloat and gas diffusion makes pumice sink

Pumice can float on water for months to years – long enough for pumice to travel across oceans and facilitate the spread of species. Long-lived pumice floatation is unexpected, however, because pumice pores are highly connected and water wets volcanic glass. As a result, observations of long floating times have not been reconciled with predictions of rapid sinking. We propose a mechanism to resolve this paradox – the trapping of gas bubbles by water within the pumice. Gas trapping refers to the isolation of gas by water within pore throats such that the gas becomes disconnected from the atmosphere and unable to escape. We use X-ray microtomography to image partially saturated pumice and demonstrate that non-condensable gas trapping occurs in both ambient temperature and hot (500 °C) pumice. Furthermore, we show that the size distribution of trapped gas clusters matches predictions of percolation theory. Finally, we propose that diffusion of trapped gas determines pumice floatation time. Experimental measurements of pumice floatation support a diffusion control on pumice buoyancy and we find that floatation time τ scales as τ∝L2Dθ2 where L is the characteristic length of pumice, D is the gas–water diffusion coefficient, and θ is pumice water saturation. A mechanistic understanding of pumice floatation is a step towards understanding how pumice is partitioned into floating and sinking components and provides an estimate for the lifetime of pumice rafts in the ocean.

58 GEOSCIENCES↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part I: Theory

We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.

East Pacific Ocean Islands↗
DOE OSTI2023

Modeling Wavelength Dependent Mid-infrared (5.5–25 μm) Optical Constants of Silicate Glasses: A Genetic Algorithm Approach

Wavelength-dependent mid-infrared (400–1,800 cm -1 ; 5.56–25 μm) optical constants (real and imaginary indices of refraction; n and k) are determined using reflectance spectra at a spectral sampling of 4 cm -1 for several silicate glasses of varying SiO 2 wt% which include (a) basaltic volcanic glass from Kilauea, Hawaii, (b) synthetic andesite, (c) two synthetic dacites, (d) obsidian volcanic glass from Mount Lassen, California, (e) synthetic rhyodacite, and (f) rhyolitic volcanic glass from Mexico. Because glasses are optically isotropic, no specific orientation was required for spectral measurements, and polished glass samples were measured at random orientations using micro-FTIR spectrometer. Lorentz-Lorenz dispersion theory and Fresnel reflectance model for high symmetry materials were used to model the optical constants by optimizing oscillator parameters in modeled spectra to match laboratory spectra. A genetic algorithm (GA) approach automatically finds the natural oscillators and their parameters for each spectrum, and then uses these parameters in a non-linear least squares optimization routine. The study compared spectral parameters such as Christiansen feature, reststrahlen bands (RBs), and the peaks centered around 860–1,100 (peak 1) and 400–480 cm -1 (peak 2) of n and k to their respective SiO 2 wt% of the glasses. CF, RBs, peak 1 of n and k, and peak 2 of n shift to higher wavenumbers with increased SiO 2 wt%, whereas peak 2 of k shifts to lower wavenumbers with increased SiO 2 wt%. Derived optical constants of these glasses will improve quantitative abundance mapping of volcanic materials on the surfaces of silicate targets in the Solar System.

58 GEOSCIENCES↗
DOE OSTI2015

High-temperature apparatus for chaotic mixing of natural silicate melts

A unique high-temperature apparatus was developed to trigger chaotic mixing at high-temperature (up to 1800 °C). This new apparatus, which we term Chaotic Magma Mixing Apparatus (COMMA), is designed to carry out experiments with high-temperature and high-viscosity (up to 10{sup 6} Pa s) natural silicate melts. This instrument allows us to follow in time and space the evolution of the mixing process and the associated modulation of chemical composition. This is essential to understand the dynamics of magma mixing and related chemical exchanges. The COMMA device is tested by mixing natural melts from Aeolian Islands (Italy). The experiment was performed at 1180 °C using shoshonite and rhyolite melts, resulting in a viscosity ratio of more than three orders of magnitude. This viscosity ratio is close to the maximum possible ratio of viscosity between high-temperature natural silicate melts. Results indicate that the generated mixing structures are topologically identical to those observed in natural volcanic rocks highlighting the enormous potential of the COMMA to replicate, as a first approximation, the same mixing patterns observed in the natural environment. COMMA can be used to investigate in detail the space and time development of magma mixing providing information about this fundamental petrological and volcanological process that would be impossible to investigate by direct observations. Among the potentials of this new experimental device is the construction of empirical relationships relating the mixing time, obtained through experimental time series, and chemical exchanges between the melts to constrain the mixing-to-eruption time of volcanic systems, a fundamental topic in volcanic hazard assessment.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The joy of transient chaos

We intend to show that transient chaos is a very appealing, but still not widely appreciated, subfield of nonlinear dynamics. Besides flashing its basic properties and giving a brief overview of the many applications, a few recent transient-chaos-related subjects are introduced in some detail. These include the dynamics of decision making, dispersion, and sedimentation of volcanic ash, doubly transient chaos of undriven autonomous mechanical systems, and a dynamical systems approach to energy absorption or explosion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Snake River Plain Play Fairway Analysis – Phase 1 Report

The Snake River volcanic province (SRP) overlies a thermal anomaly that extends deep into the mantle; it represents one of the highest heat flow provinces in North America. Our goals for this Phase 1 study are to: (1) adapt the methodology of Play Fairway Analysis for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal play fairway model for the SRP and identify the most promising plays, using software tools that are standard in the petroleum industry. The success of play fairway analysis in geothermal exploration depends critically on defining a systematic methodology that is grounded in theory (as developed within the petroleum industry over the last two decades) and within the geologic and hydrologic framework of real geothermal systems. Our preliminary assessment of the data suggests that important undiscovered geothermal resources may be located in several areas of the SRP, including the western SRP (associated with buried lineaments defined by gravity or magnetic anomalies, and capped by extensive deposits of lacustrine sediment), at lineament intersections in the central SRP (along the Banbury-Hagerman trend NW of Twin Falls, and along the northern margin of the Mt Bennett Hills-Camas Prairie area), and along the margins of the eastern SRP. Additional high temperature resources are likely associated with rhyolite domes and crypto-domes in the eastern SRP, but are masked by shallow groundwater flow leading to low upper crustal heat flow values. These blind resources may be exploitable with existing deep drilling technology. Groundwater modeling planned for later phases of the PFA project will address whether temperatures at viable producing depths are sufficient to support electricity production.

15 GEOTHERMAL ENERGY↗
DOE OSTI2015

Snake River Plain Play Fairway Analysis - Phase 1 Report

The Snake River volcanic province (SRP) overlies a thermal anomaly that extends deep into the mantle; it represents one of the highest heat flow provinces in North America. Our goals for this Phase 1 study are to: (1) adapt the methodology of Play Fairway Analysis for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal play fairway model for the SRP and identify the most promising plays, using software tools that are standard in the petroleum industry. The success of play fairway analysis in geothermal exploration depends critically on defining a systematic methodology that is grounded in theory (as developed within the petroleum industry over the last two decades) and within the geologic and hydrologic framework of real geothermal systems. Our preliminary assessment of the data suggests that important undiscovered geothermal resources may be located in several areas of the SRP, including the western SRP (associated with buried lineaments defined by gravity or magnetic anomalies, and capped by extensive deposits of lacustrine sediment), at lineament intersections in the central SRP (along the Banbury-Hagerman trend NW of Twin Falls, and along the northern margin of the Mt Bennett Hills-Camas Prairie area), and along the margins of the eastern SRP. Additional high temperature resources are likely associated with rhyolite domes and crypto-domes in the eastern SRP, but are masked by shallow groundwater flow leading to low upper crustal heat flow values. These blind resources may be exploitable with existing deep drilling technology. Groundwater modeling planned for later phases of the PFA project will address whether temperatures at viable producing depths are sufficient to support electricity production.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

The geometry of signal space: a case study of direct mapping between seismic signals and event distribution

SUMMARY Under favourable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal-space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection and scaling, with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning and possible application to studies of repeating events in induced, volcanic and glacial seismicity. Our results place limits on the widely held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗
DOE OSTI2021

LL18-Signal Analysis-NDD2Ad (Final Report)

Under favorable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially-distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection, and scaling with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning, and possible application to studies of repeating events in induced, volcanic, and glacial seismicity. Our results place limits on the widely-held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗
DOE OSTI2017

On the evolution of the plume function and entrainment in the near-source region of lazy plumes

Plumes occur in many natural and industrial settings, such as chimney smoke, volcanic eruptions and deep water oil spills. A plume function, $\unicode[STIX]{x1D6E4}$ , is used to characterize plumes and jets. The far-field behaviour of these flows has been studied in great detail while the near-field behaviour has not quite received the same attention. We examine near-field phenomena such as radial constriction, termed necking, and vortex structure formations with new high resolution direct numerical simulations. Four lazy plumes with increasing values of the source plume parameter, $\unicode[STIX]{x1D6E4}_{0}$ , are simulated. We study the evolution of entrainment and the plume function. The original assumptions, that Reynolds stresses dominate viscous shear stresses, do not hold for lazy plumes in the near field. Due to this, a deviation from self-similarity occurs initially and is corrected by a large entrainment coefficient caused by vortex stretching and compression. After correcting for the virtual origin, comparison between theory and simulations shows a monotonic decay of $\unicode[STIX]{x1D6E4}$ towards pure plume behaviour. The entrainment coefficient asymptotes to a widely accepted constant value for plumes.

Mechanics↗
DOE OSTI2015

Seismic energy data analysis of Merapi volcano to test the eruption time prediction using materials failure forecast method (FFM)

The test of eruption time prediction is an effort to prepare volcanic disaster mitigation, especially in the volcano’s inhabited slope area, such as Merapi Volcano. The test can be conducted by observing the increase of volcanic activity, such as seismicity degree, deformation and SO2 gas emission. One of methods that can be used to predict the time of eruption is Materials Failure Forecast Method (FFM). Materials Failure Forecast Method (FFM) is a predictive method to determine the time of volcanic eruption which was introduced by Voight (1988). This method requires an increase in the rate of change, or acceleration of the observed volcanic activity parameters. The parameter used in this study is the seismic energy value of Merapi Volcano from 1990 – 2012. The data was plotted in form of graphs of seismic energy rate inverse versus time with FFM graphical technique approach uses simple linear regression. The data quality control used to increase the time precision employs the data correlation coefficient value of the seismic energy rate inverse versus time. From the results of graph analysis, the precision of prediction time toward the real time of eruption vary between −2.86 up to 5.49 days.

42 ENGINEERING↗
DOE OSTI2021

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

Geochemistry & Geophysics↗
DOE OSTI2022

Estimating Explosion Yields Using Moment Tensor Solutions and Seismic Moment

We report seismic moment, a measurable and well-understood quantity of seismic sources, is used to estimate the yield of explosions. Application of such a method in the past, as in the manner of m b -derived yields, has been complicated by the effect of variations in the explosion working point, depth, and secondary source effects (such as spalling and tectonic release) on the observed moment. We start using the full (six-element) moment tensor solution, which can capture the relevant source physics and, at least in theory, better isolate the primary explosion source. The moment-to-yield ratio is then estimated using an explosion source model which, provided with emplacement conditions, can relate the two parameters. We discuss the major sources of uncertainty associated with the method, and calibrate it with chemical and nuclear explosions at the Nevada National Security Site. We then apply the method to published moment tensor solutions for the six declared North Korean nuclear explosions that occurred between 2006 and 2017. The results are mostly consistent with other yield estimates made using a variety of high-frequency methods. This technique is a new approach to estimating explosive yield and simple to implement, as much of the complexity is captured by the source models.

58 GEOSCIENCES↗
DOE OSTI2017

Source partitioning of anthropogenic groundwater nitrogen in a mixed-use landscape, Tutuila, American Samoa

This study presents a modeling framework for quantifying human impacts and for partitioning the sources of contamination related to water quality in the mixed-use landscape of a small tropical volcanic island. On Tutuila, the main island of American Samoa, production wells in the most populated region (the Tafuna-Leone Plain) produce most of the island’s drinking water. However, much of this water has been deemed unsafe to drink since 2009. Tutuila has three predominant anthropogenic non-point-groundwater-pollution sources of concern: on-site disposal systems (OSDS), agricultural chemicals, and pig manure. These sources are broadly distributed throughout the landscape and are located near many drinking-water wells. Water quality analyses show a link between elevated levels of total dissolved groundwater nitrogen (TN) and areas with high non-point-source pollution density, suggesting that TN can be used as a tracer of groundwater contamination from these sources. The modeling framework used in this study integrates land-use information, hydrological data, and water quality analyses with nitrogen loading and transport models. The approach utilizes a numerical groundwater flow model, a nitrogen-loading model, and a multi-species contaminant transport model. Nitrogen from each source is modeled as an independent component in order to trace the impact from individual land-use activities. Model results are calibrated and validated with dissolved groundwater TN concentrations and inorganic δ{sup 15}N values, respectively. Results indicate that OSDS contribute significantly more TN to Tutuila’s aquifers than other sources, and thus should be prioritized in future water-quality management efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Constraining the Magmatic System at Mount St. Helens (2004-2008) Using Bayesian Inversion With Physics-Based Models Including Gas Escape and Crystallization

Physics-based models of volcanic eruptions track conduit processes as functions of depth and time. When used in inversions, these models permit integration of diverse geological and geophysical data sets to constrain important parameters of magmatic systems. We develop a 1-D steady state conduit model for effusive eruptions including equilibrium crystallization and gas transport through the conduit and compare with the quasi-steady dome growth phase of Mount St. Helens in 2005. Viscosity increase resulting from pressure-dependent crystallization leads to a natural transition from viscous flow to frictional sliding on the conduit margin. Erupted mass flux depends strongly on wall rock and magma permeabilities due to their impact on magma density. Including both lateral and vertical gas transport reveals competing effects that produce nonmonotonic behavior in the mass flux when increasing magma permeability. Using this physics-based model in a Bayesian inversion, we link data sets from Mount St. Helens such as extrusion flux and earthquake depths with petrological data to estimate unknown model parameters, including magma chamber pressure and water content, magma permeability constants, conduit radius, and friction along the conduit walls. Even with this relatively simple model and limited data, we obtain improved constraints on important model parameters. We find that the magma chamber had low (<5 wt %) total volatiles and that the magma permeability scale is well constrained at ~10 –11.4 m 2 to reproduce observed dome rock porosities. Here, compared with previous results, higher magma overpressure and lower wall friction are required to compensate for increased viscous resistance while keeping extrusion rate at the observed value.

58 GEOSCIENCES↗