Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “eruptions”

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.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Entropy-based feature selection for capturing impacts in Earth system models with abrupt forcing

This paper presents the development of a new entropy-based feature selection method for identifying and quantifying impacts. Here, impacts are defined as statistically significant differences in spatio-temporal fields when comparing datasets with and without an external forcing in an Earth system model. Temporal feature selection is performed by first computing the cross-fuzzy entropy to quantify similarity of patterns between two datasets and then applying changepoint detection to identify regions of statistically constant entropy. The method is used to capture temperate north surface cooling from a 9-member simulation ensemble of the Mt. Pinatubo volcanic eruption, which injected 10 Tg of SO 2 into the stratosphere. The results estimate a mean difference decrease in near surface air temperature of -0.560 K with a 99% confidence interval between -0.864 K and -0.257 K between April and November of 1992, one year following the eruption. A sensitivity analysis with decreasing SO 2 injection revealed that the impact is statistically significant at 5 Tg but not at 3 Tg. Using identified features, a dependency graph model based on a 9-day lag had significantly fewer nodes than a graph based on monthly means. Furthermore, this demonstrates our method’s ability to perform dimension reduction while still uncovering source-to-impact pathways.

Changepoint detection↗
DOE OSTI2020

Caltech Lab Experiments and the Insights They Provide Into Solar Corona Phenomena

A comprehensive overview of two decades of Caltech experiments relevant to solar corona physics is presented. The extent to which the experiments scale to the solar corona, the basic configurations and operation, and the importance of the magnetic force J × B common to all the experiments is discussed. Summaries are given of the various configurations used, the main observations, and interpretations of these observations, including new models developed to provide these interpretations. Topics include observations and explanations for flux rope self-collimation, axial flows along flux ropes, eruption of arched flux ropes, strapping magnetic fields that inhibit eruption, the torus instability, and effects such as X-ray emission of a kink-driven secondary Rayleigh-Taylor instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Narrow-Band Least-Squares Infrasound Array Processing

Infrasound data from arrays can be used to detect, locate, and quantify a variety of natural and anthropogenic sources from local to remote distances. However, many array processing methods use a single broad frequency range to process the data, which can lead to signals of interest being missed due to the choice of frequency limits or simultaneous clutter sources. In this work, we introduce a new open-source Python code that processes infrasound array data in multiple sequential narrow frequency bands using the least-squares approach. We test our algorithm on a few examples of natural sources (volcanic eruptions, mass movements, and bolides) for a variety of array configurations. Our method reduces the need to choose frequency limits for processing, which may result in missed signals, and it is parallelized to decrease the computational burden. Improvements of our narrow-band least-squares algorithm over broad-band least-squares processing include the ability to distinguish between multiple simultaneous sources if distinct in their frequency content (e.g., microbarom or surf vs. volcanic eruption), the ability to track changes in frequency content of a signal through time, and a decreased need to fine-tune frequency limits for processing. We incorporate a measure of planarity of the wavefield across the array (sigma tau, στ) as well as the ability to utilize the robust least trimmed squares algorithm to improve signal processing and insight into array performance. Our implementation allows for more detailed characterization of infrasound signals recorded at arrays that can improve monitoring and enhance research capabilities.

58 GEOSCIENCES↗
DOE OSTI2020

Chandra Resolves the Double FU Orionis System RNO 1B/1C in X-Rays

We present new Chandra X-ray observations of the close pair of young stars RNO 1B and 1C (6″ separation) located in the L1287 cloud. RNO 1B erupted in 1978–1990 and is classified as an FU Orionis star (FUor). RNO 1C also shows most of the properties of an FUor but no eruption has been seen yet. Only a few dozen FUors are known and the presence of two such objects with a small angular separation is rare, suggesting a common origin. Both stars were faintly detected by Chandra and we summarize their X-ray properties within the framework of other previously detected FUors. We also report other X-ray detections in L1287 including the deeply embedded young star RNO 1G, the jet-like radio source VLA 3, and an enigmatic hard flaring source with no Two Micron All Sky Survey counterpart that was only detected in the second of the two Chandra exposures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Observed and Projected Changes of Global Monsoons: Current Status and Future Perspectives

The global monsoon system, encompassing the Asian-Australian, African, and American monsoons, sustains two-thirds of the world’s population by regulating water resources and agriculture. Monsoon anomalies pose severe risks, including floods and droughts. Recent research associated with the implementation of the Global Monsoons Model Intercomparison Project under the umbrella of CMIP6 has advanced our understanding of its historical variability and driving mechanisms. Observational data reveal a 20th-century shift: increased rainfall pre-1950s, followed by aridification and partial recovery post-1980s, driven by both internal variability (e.g., Atlantic Multidecadal Oscillation) and external forcings (greenhouse gases, aerosols), while ENSO drives interannual variability through ocean-atmosphere interactions. Future projections under greenhouse forcing suggest long-term monsoon intensification, though regional disparities and model uncertainties persist. Models indicate robust trends but struggle to quantify extremes, where thermodynamic effects (warming-induced moisture rise) uniformly boost heavy rainfall, while dynamical shifts (circulation changes) create spatial heterogeneity. Volcanic eruptions and proposed solar radiation modification (SRM) further complicate predictions: tropical eruptions suppress monsoons, whereas high-latitude events alter cross-equatorial flows, highlighting unresolved feedbacks. The emergent constraint approach is booming in terms of correcting future projections and reducing uncertainty with respect to the global monsoons. Critical challenges remain. Model biases and sparse 20th-century observational data hinder accurate attribution. The interplay between natural variability and anthropogenic forcings, along with nonlinear extreme precipitation risks under warming, demands deeper mechanistic insights. Additionally, SRM’s regional impacts and hemispheric monsoon interactions require systematic evaluation. Addressing these gaps necessitates enhanced observational networks, refined climate models, and interdisciplinary efforts to disentangle multiscale drivers, ultimately improving resilience strategies for monsoon-dependent regions.

climate extreme events↗
DOE OSTI2021

Symmetry-breaking-induced rare fluctuations in a time-delay dynamic system

Inspired by the experimental and numerical findings, we study the dynamic instabilities of two coupled nonlinear delay differential equations that are used to describe the coherent oscillations between the top and bottom boundary layers in turbulent Rayleigh–Bénard convection. By introducing two sensitivity parameters for the instabilities of the top and bottom boundary layers, we find three different types of solutions, namely in-phase single-period oscillations, multi-period oscillations and chaos. The chaos solution contains rare but large amplitude fluctuations. The statistical properties of these fluctuations are consistent with those observed in the experiment for the massive eruption of thermal plumes, which causes random reversals of the large-scale circulation in turbulent Rayleigh–Bénard convection. Finally, our study thus provides new insights into the origin of rare massive eruptions and sudden changes of large-scale flow pattern that are often observed in convection systems of geophysical and astrophysical scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Seismic tomography of compressional wave velocity and attenuation structure for Makushin Volcano, Alaska

Located in the central portion of the Alaska-Aleutian subduction zone, Makushin Volcano is among the most active volcanoes in the United States and has been classified as a very high threat by the National Volcano Early Warning System based on its eruptive history and proximity to the City of Unalaska and international air routes. In 2015, five standalone seismic stations and three mini seismic arrays of 15 stations each were deployed on Unalaska Island to supplement the Alaska Volcano Observatory permanent seismic network. This temporary array was operational for one year. Taking advantage of the increased azimuthal coverage and the array's increased earthquake detection capability, we developed body-wave Vp and Vp/Vs seismic images of the velocity structure beneath the volcano. Our results show a complex structure with the upper 4 km of the crust. The shallow high-Vp features possibly delineate remnant magma pathways or conduits. Low-Vp regions are found east of the caldera at approximately 4–8 km depth. This is in agreement with previous seismic tomographic work and InSAR-based geodetic models, which had identified this region as a possible long-term source of magma. A three-dimensional compressional wave attenuation model complemented the velocity models. In general, we observe regions of corresponding low P-wave velocity and high attenuation east of the summit caldera between 5 and 7 km depth, supporting the presence of a melt-rich zone. The absolute Vp values (5.5–6.0 km/s) in this region are consistent with the presence of a magma body with a mafic signature. Here, we suggest that a triggering mechanism that involves mafic injections into a more evolved shallow magma chamber is representative of the current magmatic system below the volcano, similarly to what petrological studies have inferred for past eruptions.

58 GEOSCIENCES↗
DOE OSTI2021

Recent seismicity at Ceboruco Volcano (Mexico)

Magma and related hydrothermal fluid movement, and their interaction with solid rock, in active volcanic regions, generate a wide variety of seismic waves whose characterization can mitigate the risk of a potential eruption. Located in the western region of the Trans-Mexican Volcanic Belt, Ceboruco Volcano, whose last eruptive period was 1870–1875, is considered to be one of the most hazardous volcanoes in Mexico. We have conducted a detailed study of the seismicity in the surroundings of Ceboruco's volcanic edifice to assess the current state of this volcano. A dense temporary seismic network with 25 seismic stations in an area of 16 km × 16 km was deployed between November 2016 and July 2017, as part of the P-24 project of the CeMIEGeo consortium; this effort has allowed the detection of 81 earthquakes concentrated beneath the crater with depths between 4 and 8 km. In this study, we observe that the recorded seismicity occurs in swarms, and we specifically identify four sequences that we characterize in detail via the first focal mechanisms available for this volcano. Our results suggest a change in the local seismicity distribution compared to earlier observations, which reported seismic activity near the volcano edifice associated with fluid migration along zones of weakness related to the extensional stresses of the Tepic-Zacoalco rift. The changes in seismic patterns and obtained focal mechanisms are consistent with observed fluid effects at many geothermal sites worldwide, but also could suggest resumption of activity at this currently dormant volcano.

58 GEOSCIENCES↗
DOE OSTI2022

Contrasting Stratospheric Smoke Mass and Lifetime From 2017 Canadian and 2019/2020 Australian Megafires: Global Simulations and Satellite Observations

Stratospheric injections of carbonaceous aerosols and combustion gases by extreme wildfires have become increasingly common. Recent “megafires,” particularly large and intense fires, delivered particulate burdens to the lower stratosphere comparable to those of moderate volcanic eruptions. The 2017 Canadian megafire generated four large Pyrocumulonimbi (pyroCbs), injecting up to ≈0.3 Tg of smoke in the lower stratosphere. Even more extreme, the 2019/2020 Australian event produced a pyroCb activity resulting in stratospheric smoke intrusions of ≈1 Tg. To understand their contrasting behavior, we present global climate simulations of the atmospheric response to these events, applying smoke burdens informed by remote observations. Model outcomes, compared to satellite data of smoke transport, reproduce reasonably well the initial plume rise, at 0.2–0.3 km/day, attaining heights of ≈20 km in Canada and above 30 km in Australia. Global dispersal of the plume occurs within about 3 weeks in both cases, consistent with observations. Smoke removal timescales, ≈5 months for the Canadian megafire, agree with remote measurements. During the Australian megafire, observations indicate stratospheric injections three times as large, and models predict comparatively longer smoke lifetimes, ≈16 months. After the latter event, atmospheric optical depths and radiative cooling achieved values close to those measured following the Pinatubo eruption. Sensitivity tests of model assumptions indicate, in accord with prior studies, that smoke burden, injection heights, and black carbon content can determine plume evolution and possible climate impacts. An empirical relation between peak heights of stratospheric plumes and lifetimes is derived that can help assess megafire impacts on the stratosphere, climate and the Earth system

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Solar, Atmospheric, and Volcanic Impacts on 10 Be Depositions in Greenland and Antarctica During the Last 100 Years

Cosmogenic radionuclides (e.g., 10 Be) from ice cores are a powerful tool for solar reconstructions back in time. However, superimposed on the solar signal, other factors like weather/climate and volcanic influences on 10 Be can complicate the interpretation of 10 Be data. A comprehensive study of 10 Be records over the recent period, when atmospheric 10 Be production and meteorological conditions are relatively well-known, can improve our interpretation of 10 Be records. Here we conduct a systematic study of the production and climate/volcanic signals in Antarctica and Greenland 10 Be records, including a new 10 Be record from the East GReenland Ice-core Project site. Greenland and Antarctica records show significant decreasing trends (5%–6.5%/decade) for 1900–1950, which is comparable with the expected production rate inferred from sunspot observations. By comparing 10 Be records with reanalysis data and atmospheric circulation patterns, 10 Be records from Southern/Southeastern Greenland are significantly correlated with the Scandinavia pattern. Stacking 10 Be records from different locations can enhance the production signal. However, this approach is not always straightforward as uncertainties in some records can lead to a weaker solar signal. A strategy can be employed to select records for the bipolar stack by comparing Greenland records with Antarctica records, assuming the shared signal is a production signal. Finally, we observe significant increases (36%–64%) in 10 Be depositions in Greenland related to the Agung eruption. This large increase in Greenland 10 Be records after the Agung eruption, could be partly explained by the enhanced air mass transport from mid-latitudes coinciding with the decreased precipitation en-route.

10Be↗
DOE OSTI2023

Perturbation of Tropical Stratospheric Ozone Through Homogeneous and Heterogeneous Chemistry Due To Pinatubo

Abstract The Pinatubo eruption in 1991 injected 10–20 Tg SO 2 into the stratosphere, which formed sulfate aerosols through oxidation. Our modeling results show that volcanic heating significantly perturbs the heterogeneous and homogeneous chemistry including NO x and HO x catalytic cycles in the tropical stratosphere. The simulated tropical chemical ozone tendency is positive at 20 mb while negative at 10 mb in the tropics. The simulated ozone chemical tendency is of the same magnitude as the dynamical ozone tendency caused by the accelerated tropical upwelling, but with the opposite sign. Our study finds that the tropical ozone chemical tendency due to homogeneous chemistry becomes more important than heterogeneous chemistry 3 months after eruption. Sensitivity simulations further suggest that the tropical ozone tendency through heterogeneous chemistry is saturated when the injected amount exceeds 2 Tg.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Space‐Time Causal Discovery in Earth System Science: A Local Stencil Learning Approach

Causal discovery tools enable scientists to infer meaningful relationships from observational data, spurring advances in fields as diverse as biology, economics, and climate science. Despite these successes, the application of causal discovery to space-time systems remains immensely challenging due to the high-dimensional nature of the data. For example, in climate sciences, modern observational temperature records over the past few decades regularly measure thousands of locations around the globe. To address these challenges, we introduce Causal Space-Time Stencil Learning (CaStLe), a novel meta-algorithm for discovering causal structures in complex space-time systems. CaStLe leverages regularities in local space-time dependencies to learn governing global dynamics. This local perspective eliminates spurious confounding and drastically reduces sample complexity, making space-time causal discovery practical and effective. For causal discovery, CaStLe flexibly accepts any appropriately adapted time series causal discovery algorithm to recover local causal structures. These advances enable causal discovery of geophysical phenomena that were previously unapproachable, including non-periodic, transient phenomena such as volcanic eruption plumes. Regularities in local space-time dependencies are transformed into informative spatial replicates, which actually improve CaStLe's performance when applied to ever-larger spatial grids. We successfully apply CaStLe to discover the atmospheric dynamics governing the climate response to the 1991 Mount Pinatubo volcanic eruption. We provide validation experiments to demonstrate the effectiveness of CaStLe over existing causal-discovery frameworks on a range of geophysics-inspired benchmarks while identifying the method's limitations and domains where its assumptions may not hold.

Nichol, J. Jake [Univ. of New Mexico, Albuquerque,↗
DOE OSTI2022

Experimental measurements of the viscosity and melt structure of alkali basalts at high pressure and temperature

Abstract Volcanic eruptions are shallow phenomena that represent the final stage of density- and viscosity- driven processes of melt migration from source rocks at upper mantle depths. In this experimental study, we investigated the effect of pressure (0.7–7.0 GPa) and temperature (1335–2000 °C) on the viscosity and the atomic melt structure of a synthetic anhydrous primitive alkaline basalt, an analogue of the pre-eruptive magma that likely feeds the Campi Flegrei Volcanic District at present day. Obtained viscosities (0.5–3.0 Pa s), mobility (0.1–0.4 g cm 3 Pa −1 s −1 ) and ascent velocity (1.5–6.0 m yr −1 ) are presented to support geochemical and geophysical observations of Campi Flegrei as a critical volcanic district currently undergoing gradual magma recharge at depth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Non-ideal MHD growth of current interchange tearing modes at plasma edge and response to externally-imposed flow

Current interchange tearing modes (CITMs) have been proposed to explain an intermittent eruption of electric current toward a divertor at an edge region of a tokamak device. Extended magnetohydrodynamic (MHD) simulations combined with the Scrape-Off-Layer diffusion model, which has been developed by the authors [Miura et al., Phys. Plasmas 24, 092111 (2017)], are carried out in order to study the growth of a CITM under influences of two-fluid and gyro-viscous effects, or non-ideal MHD effects. Numerical simulations both with and without an externally imposed azimuthal flow show that the growth of a CITM is observed for a relatively small flow velocity, and that a CITM can be a candidate mechanism for the intermittent eruption of the current in a tokamak, whether the growth is under influences of non-ideal MHD effects or not. Here, the growth can be suppressed when an externally imposed azimuthal flow is sufficiently large to cause a finite radial displacement because such a displacement prevents a transition from an interchange to a tearing mode that is the key process of a CITM. Furthermore, it is also found that a stripe pattern similar to so-called streamers is formed in the course of the CITM growth. The basic nature of this streamer-like structure is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Creating and studying a scaled interplanetary coronal mass ejection

The Sun, being an active star, undergoes eruptions of magnetized plasma that reach the Earth and cause the aurorae near the poles. These eruptions, called coronal mass ejections (CMEs), send plasma and magnetic fields out into space. CMEs that reach planetary orbits are called interplanetary coronal mass ejections (ICMEs) and are a source of geomagnetic storms, which can cause major damage to our modern electrical systems with limited warning. To study ICME propagation, we devised a scaled experiment using the Big Red Ball (BRB) plasma containment device at the Wisconsin Plasma Physics Laboratory. These experiments inject a compact torus of plasma as an ICME through an ambient plasma inside the BRB, which acts as the interplanetary medium. Magnetic and temperature probes provide three-dimensional magnetic field information in time and space, as well as temperature and density as a function of time. Using this information, we can identify features in the compact torus that are consistent with those in real ICMEs. We also identify the shock, sheath, and ejecta similar to the structure of an ICME event. This experiment acts as a first step to providing information that can inform predictive models, which can give us time to shield our satellites and large electrical systems in the event that a powerful ICME were to strike.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The possible disappearance of a massive star in the low-metallicity galaxy PHL 293B

ABSTRACT We investigate a suspected very massive star in one of the most metal-poor dwarf galaxies, PHL 293B. Excitingly, we find the sudden disappearance of the stellar signatures from our 2019 spectra, in particular the broad H lines with P Cygni profiles that have been associated with a massive luminous blue variable (LBV) star. Such features are absent from our spectra obtained in 2019 with the Echelle Spectrograph for Rocky Exoplanet- and Stable Spectroscopic Observation and X-shooter instruments of the European Southern Observatory’s Very Large Telescope. We compute radiative transfer models using cmfgen, which fit the observed spectrum of the LBV and are consistent with ground-based and archival Hubble Space Telescope photometry. Our models show that during 2001–2011, the LBV had a luminosity L* = 2.5–3.5 × 106 L⊙, a mass-loss rate $\dot{M} = 0.005{-}0.020 ~{\rm M}_{\odot }$ yr−1, a wind velocity of 1000 km s−1, and effective and stellar temperatures of Teff = 6000–6800 and T* = 9500–15 000 K. These stellar properties indicate an eruptive state. We consider two main hypotheses for the absence of the broad emission components from the spectra obtained since 2011. One possibility is that we are seeing the end of an LBV eruption of a surviving star, with a mild drop in luminosity, a shift to hotter effective temperatures, and some dust obscuration. Alternatively, the LBV could have collapsed to a massive black hole without the production of a bright supernova.

Allan, Andrew P.↗
DOE OSTI2022

Experimental calibration of an Fe3+/Fe2+-in-amphibole oxybarometer and its application to shallow magmatic processes at Shiveluch Volcano, Kamchatka

Abstract Oxygen fugacity is an important but difficult parameter to constrain for primitive arc magmas. In this study, the partitioning behavior of Fe3+/Fe2+ between amphibole and glass synthesized in piston-cylinder and cold-seal apparatus experiments is developed as an oxybarometer, applicable to magmas ranging from basaltic to dacitic composition. The partitioning of Fe2+ is strongly dependent on melt polymerization; the relative compatibility of Fe2+ in amphibole decreases with increasing polymerization. The Fe2+/Mg distribution coefficient between amphibole and melt is a relatively constant value across all compositions and is, on average, 0.27. The amphibole oxybarometer is applied to amphibole in mafic enclaves, cumulates, and basaltic tephra erupted from Shiveluch volcano in Kamchatka with measured Fe3+/FeTotal. An average Fe3+/Fe2+ amphibole-glass distribution coefficient for basalt is used to convert the Fe3+/FeTotal of amphibole in samples from Shiveluch to magmatic oxygen fugacity relative to NNO. The fO2 of primitive melts at the volcano is approximately NNO+2 and is faithfully recorded in amphibole from an amphibole-rich cumulate and the basaltic tephra. Apparently, higher fO2 recorded by amphibole in mafic enclaves likely results from partial dehydrogenation of amphibole during residence in a shallow andesite storage region. We identify three pulses of mafic magma recharge within two weeks of, a month before, and two to three months before the eruption and find that, at each of these times, the host andesite was recharged by at least two magmas at varying stages of differentiation. Application of the amphibole oxybarometer not only gives insight into magmatic fO2 but also potentially details of shallow magmatic processes.

Geochemistry & Geophysics↗
DOE OSTI2024

CLimate Impact: Determining Etiology thRough pAthways (CLDERA)

Climate impacts have broad economic, health, political, and national security ramifications. Societally relevant impacts are typically farther downstream, are the product of multiple interacting processes, and can arise over small regions and timeframes because their sources are short-term and localized. Short-term forcings (as can be seen in volcanic eruptions, climatic tipping points (e.g., the collapse of rainforests or the disappearance of sea ice), or in increasingly plausible climate interventions) fundamentally possess low signal-to-noise and could benefit from accounting for the multiple conditional processes through which a downstream impact arises. Under the Grand Challenge LDRD CLDERA (CLimate impacts: Discovering Etiology thRough pAthways), we have developed tools to enable downstream impact attribution from geographically and temporally localized source forcings in the climate. CLDERA developed methods that can distinguish how a localized source drives the climate system to respond with particular impacts. The how is embodied in pathways – the spatio-temporally evolving chain of physical processes that connects a source to a series of increasingly distant impacts. Novel analytic methods in pursuit of downstream impact attribution were developed and demonstrated on simulations and observations of the 1991 eruption of Mt. Pinatubo in the Philippines. As described within this report we have • developed stratospheric expertise and aerosol modeling capabilities in E3SM, • created original methods to detect and model pathways from source-to-impact, and • advanced climate attribution through novel methods, cases, and approaches. Further, CLDERA developed a tiered verification process consisting of controlled datasets to prototype, verify, and refine the original method development. CLDERA increased Sandia’s footprint in the climate analytics community and developed new climate collaborations whilst also creating a cadre of climate analysts at Sandia. The products from CLDERA have been extensive with a total of 9 journal articles published, 12 articles submitted and under review, and an additional 8 articles in preparation. We have produced 1750 simulated years and developed 9 code-bases. This report details these accomplishments and serves as a summary of the work completed during the CLDERA Grand Challenge.

54 ENVIRONMENTAL SCIENCES↗