Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “height estimation”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Marine Boundary Layer Cloud Boundaries and Phase Estimation Using Airborne Radar and In Situ Measurements During the SOCRATES Campaign over Southern Ocean

The Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) was an aircraft-based campaign (15 January–26 February 2018) that deployed in situ probes and remote sensors to investigate low-level clouds over the Southern Ocean (SO). A novel methodology was developed to identify cloud boundaries and classify cloud phases in single-layer, low-level marine boundary layer (MBL) clouds below 3 km using the HIAPER Cloud Radar (HCR) and in situ measurements. The cloud base and top heights derived from HCR reflectivity, Doppler velocity, and spectrum width measurements agreed well with corresponding lidar-based and in situ estimates of cloud boundaries, with mean differences below 100 m. A liquid water content–reflectivity (LWC-Z) relationship, LWC = 0.70Z0.29, was derived to retrieve the LWC and liquid water path (LWP) from HCR profiles. The cloud phase was classified using HCR measurements, temperature, and LWP, yielding 40.6% liquid, 18.3% mixed-phase, and 5.1% ice samples, along with drizzle (29.1%), rain (3.2%), and snow (3.7%) for drizzling cloud cases. The classification algorithm demonstrates good consistency with established methods. This study provides a framework for the boundary and phase detection of MBL clouds, offering insights into SO cloud microphysics and supporting future efforts in satellite retrievals and climate model evaluation.

MBL clouds over Southern Ocean↗
DOE OSTI2023

From Hydrometeor Size Distribution Measurements to Projections of Wind Turbine Blade Leading-Edge Erosion

Wind turbine blade leading-edge erosion (LEE) is a cause of increased operation and maintenance costs and decreased annual energy production. Thus, detailed, site-specific quantification of likely erosion conditions are critically needed to inform wind plant owner/operator decisions regarding mitigation strategies. Estimating the damage potential at a wind plant site requires accurate measurement of precipitation intensity, phase, droplet size distributions, wind speeds and their joint statistics. The current work quantifies the effect of disdrometer type on the characterization of LEE potential at a site in the US Southern Great Plains. using observations from three co-located disdrometers (an optical, an impact and a video disdrometer), along with hub-height wind-speed observations from a Doppler lidar and two LEE models: a kinetic energy model and the Springer model. Estimates of total kinetic energy of hydrometeor impacts over the four-year study period vary by as much as 38%, and coating lifetime derived from accumulated distance-to-failure estimates from the Springer model differ by an even greater amount, depending on disdrometer type. Damage potential at this site is concentrated in time, with 50% of impact kinetic energy occurring in 6–12 h per year, depending on which set of disdrometer observations is used. Rotor-speed curtailment during the most erosive 0.1–0.2% of 10 min periods is found to increase blade lifetimes and lead to the lowest levelized cost of energy.

17 WIND ENERGY↗
NASA NTRS1993

Comparison of dynamic height variations in the tropical Atlantic during 1987-1989 as viewed in sections hydrography and Geosat altimetry

During the 1980s the Soviet Union conducted a series of 25 hydrographic cruises in the tropical Atlantic under the Sections (acronym for Energetically Active Zones of the Ocean and Climate Variability) program. As one of its objectives, Sections has produced the most extensive hydrographic data set ever collected in the tropical Atlantic. Here we examine dynamic topography from seven cruises during the years 1987 through 1989 and compare them to corresponding maps of sea level from the Geosat satellite mission. The two independent estimates of sea level variation are well correlated, with the sea level fields explaining approximately half of the variance of dynamic height fields. There is a seasonal variation to the agreement, so that the correlation increases somewhat in summer and fall, as the variability of the height field increases. Wavenumber spectra of both dynamic height and sea level are quite red for wavelengths greater than 400-600 km. Smaller wavelengths cannot be resolved by the hydrographic sampling program. A comparison of the dynamic height and sea level spectra indicates that sea level has somewhat less variability than dynamic height at all wavenumbers.

Carton, James A.↗
DOE OSTI2025

NANT Site - Microwave Radiometer Thermodynamic Retrievals TROPoe / Derived Data

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous brightness temperature measurements at 35 channels observed with a microwave radiometer MP3000A operated by NOAA Physical Sciences Laboratory on Nantucket Island for WFIP3. Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a sensor attached to the MWR housing. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗
DOE OSTI2025

NOAA PSL Microwave Radiometer Block Island / Thermodynamic retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous brightness temperature measurements at 35 channels observed with a microwave radiometer MP3000A operated by NOAA Physical Sciences Laboratory on the Block Island for WFIP3. Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a sensor attached to the MWR housing. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗
NASA NTRS1990

Laboratory simulation of the rocket motor thrust as a follower force

Ground tests of solid propellant rocket motors have shown that metal-containing propellants produce various amounts of slag (primarily aluminum oxide), which is trapped in the motor case causing a loss of specific impulse. Although not yet definitely established, the presence of a liquid pool of slag also may contribute to nutational instabilities that have been observed with certain spin-stabilized, upper-stage vehicles. Because of the rocket's axial acceleration - absent in the ground tests - estimates of in-flight slag mass have been very uncertain. Yet such estimates are needed to determine the magnitude of the control authority of the systems required for eliminating the instability. A test rig with an eccentrically mounted hemispherical bowl was designed and built that incorporates a follower force that properly aligns the thrust vector along the axis of spin. A program that computes the motion of a point mass in the spinning and precessing bowl was written. Using various rpm, friction factors, and initial starting conditions, plots were generated showing the trace of the point mass around the inside of the fuel tank. The apparatus will be used extensively during the 1990 to 1991 academic year and incorporate future design features such as a variable nutation angle and a film height measuring instrument. Data obtained on the nutational instability characteristics will be used to determine order-of-magnitude estimates of control authority needed to minimize the sloshing effect.

Source record↗
NASA NTRS2017

CME Velocity and Acceleration Error Estimates Using the Bootstrap Method

The bootstrap method is used to determine errors of basic attributes of coronal mass ejections (CMEs) visually identified in images obtained by the Solar and Heliospheric Observatory (SOHO) mission's Large Angle and Spectrometric Coronagraph (LASCO) instruments. The basic parameters of CMEs are stored, among others, in a database known as the SOHO/LASCO CME catalog and are widely employed for many research studies. The basic attributes of CMEs (e.g. velocity and acceleration) are obtained from manually generated height-time plots. The subjective nature of manual measurements introduces random errors that are difficult to quantify. In many studies the impact of such measurement errors is overlooked. In this study we present a new possibility to estimate measurements errors in the basic attributes of CMEs. This approach is a computer-intensive method because it requires repeating the original data analysis procedure several times using replicate datasets. This is also commonly called the bootstrap method in the literature. We show that the bootstrap approach can be used to estimate the errors of the basic attributes of CMEs having moderately large numbers of height-time measurements. The velocity errors are in the vast majority small and depend mostly on the number of height-time points measured for a particular event. In the case of acceleration, the errors are significant, and for more than half of all CMEs, they are larger than the acceleration itself.

Michalek, Grzegorz↗
NASA NTRS2006

Merging IceSAT GLAS and Terra MODIS Data in Order to Derive Forest Type Specific Tree Heights in the Central Siberian Boreal Forest

Mapping of boreal forest's type, biomass, and other structural parameters are critical for understanding of the boreal forest's significance in the carbon cycle, its response to and impact on global climate change. We believe the nature of the forest structure information available from MISR and GLAS can be used to help identify forest type, age class, and estimate above ground biomass levels beyond that now possible with MODIS alone. The ground measurements will be used to develop relationships between remote sensing observables and forest characteristics and provide new information for understanding forest changes with respect to environmental change. Lidar is a laser altimeter that determines the distance from the instrument to the physical surface by measuring the time elapsed between the pulse emission and the reflected return. Other studies have shown that the returned signal may identify multiple returns originating from trees, building and other objects and permits the calculation of their height. Studies using field data have shown that lidar data can provide estimates of structural parameters such as biomass, stand volume and leaf area index and allows remarkable differentiation between primary and secondary forest. NASA's IceSAT Geoscience Laser Altimeter System (GLAS) was launched in January 2003 and collected data during February and September of that year. This study used data acquired over our study sites in central Siberia to examine the GLAS signal as a source of forest height and other structural characteristics. The purpose of our Siberia project is to improve forest cover maps and produce above-ground biomass maps of the boreal forest in Northern Eurasia from MODIS by incorporating structural information inherent in the Terra MISR and ICESAT Geoscience Laser Altimeter System (GLAS) instruments. A number of forest cover classifications exist for the boreal forest. We believe the limiting factor in these products is the lack of structural information, particularly in the vertical dimension. The emphasis of this project is to improve upon satellite maps of boreal forest structure parameters (i.e. height and biomass) using temporal, multi-angle, and vertical profile information of GLAS data. The existing and near future lidar data is useful for demonstrating these techniques and pursuing current estimates. Future lidar missions may be several years in the future, so we will work other new data sets that may aide in biomass estimates such as ALOS PALSAR We will continue this work to produce an accurate map of current above ground forest phytomass/carbon storage possible for the study area. We plan to develop, test, and integrate remote sensing methods for extracting forest canopy structure measures. We are compiling our field measurements and will compare them with the remote sensing methods where possible. We also be able to produce a realistic error bound on the remotely sensed carbon estimates.

Ranson, K. Jon↗
NASA NTRS2011

Video Meteor Fluxes

The flux of meteoroids, or number of meteoroids per unit area per unit time, is critical for calibrating models of meteoroid stream formation and for estimating the hazard to spacecraft from shower and sporadic meteors. Although observations of meteors in the millimetre to centimetre size range are common, flux measurements (particularly for sporadic meteors, which make up the majority of meteoroid flux) are less so. It is necessary to know the collecting area and collection time for a given set of observations, and to correct for observing biases and the sensitivity of the system. Previous measurements of sporadic fluxes are summarized in Figure 1; the values are given as a total number of meteoroids striking the earth in one year to a given limiting mass. The Gr n et al. (1985) flux model is included in the figure for reference. Fluxes for sporadic meteoroids impacting the Earth have been calculated for objects in the centimeter size range using Super-Schmidt observations (Hawkins & Upton, 1958); this study used about 300 meteors, and used only the physical area of overlap of the cameras at 90 km to calculate the flux, corrected for angular speed of meteors, since a large angular speed reduces the maximum brightness of the meteor on the film, and radiant elevation, which takes into account the geometric reduction in flux when the meteors are not perpendicular to the horizontal. They bring up corrections for both partial trails (which tends to increase the collecting area) and incomplete overlap at heights other than 90 km (which tends to decrease it) as effects that will affect the flux, but estimated that the two effects cancelled one another. Halliday et al. (1984) calculated the flux of meteorite-dropping fireballs with fragment masses greater than 50 g, over the physical area of sky accessible to the MORP fireball cameras, counting only observations in clear weather. In the micron size range, LDEF measurements of small craters on spacecraft have been used to estimate the flux (Love & Brownlee, 1993); here the physical area of the detector is well known, but the masses depend strongly on the unknown velocity distribution. In the same size range, Thomas & Netherway (1989) used the narrow-beam radar at Jindalee to calculate the flux of sporadics. In between these very large and very small sizes, a number of video and photographic observations were reduced by Ceplecha (2001). These fluxes were calculated (details are given in Ceplecha, 1988) taking the Halliday et al. (1984) MORP fireball fluxes, slightly corrected in mass, as a calibration, and adjusting the flux of small cameras to overlap with the number/mass relation from that work.

Campbell-Brown, M. D.↗
DOE OSTI2025

NWTC Site 4.0 - NREL ASSIST (SN10) / Thermodynamic retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe v0.12 (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NREL ASSIST-II (SN 10) infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a Vaisala CL51 ceilometer. The full pipeline for running the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. Met data were not ingested. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Denver, CO.

17 WIND ENERGY↗
DOE OSTI2025

NWTC Site 3.2 - NREL ASSIST (SN12) / Thermodynamic retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe v0.12 (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NREL ASSIST-II (SN 12) infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a Vaisala CL51 ceilometer. The full pipeline for running the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. Met data were not ingested. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Denver, CO.

17 WIND ENERGY↗
DOE OSTI2025

Title NWTC Site 3.2 - NREL ASSIST (SN11) / Thermodynamic retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe v0.12 (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NREL ASSIST-II (SN 11) infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a Vaisala CL51 ceilometer. The full pipeline for running the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. Met data were not ingested. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Denver, CO.

17 WIND ENERGY↗
NASA NTRS2000

Apse-Alignment of the Uranian Rings

An explanation of the dynamical mechanism for apse-alignment of the eccentric Uranian rings is necessary before observations can be used to determine properties such as ring masses, particle sizes, and elasticities. The leading model relies on the ring self-gravity to accomplish this task, yet it yields equilibrium masses which are not in accord with Voyager radio measurements. We explore possible solutions such that the self-gravity and the collisional terms are both involved in the process of apse-alignment. We consider limits that correspond to a hot and a cold ring, and show that pressure terms may play a significant role in the equilibrium conditions for the narrow Uranian rings. In the cold ring case, where the scale height of the ring near periapse is comparable to the ring particle size, we introduce a new pressure correction pertaining to a region of the ring where the particles are locked in their relative positions and jammed against their neighbors, and the velocity dispersion is so low that the collisions are nearly elastic. In this case, we find a solution such that the ring self-gravity maintains apse-alignment against both differential precession (m = 1 mode) and the fluid pressure. We apply this model to the Uranian alpha ring, and show that, compared to the previous self-gravity model, the mass estimate for this ring increases by an order of magnitude. In the case of a hot ring, where the scale height can reach a value as much as fifty times larger than a particle size, we find velocity dispersion profiles that result in pressure forces which act in such a way as to alter the ring equilibrium conditions, again leading to a ring mass increase of an order of magnitude; however, such a velocity dispersion profile would require a different mechanism than is currently envisioned for establishing heating/cooling balance in a finite-sized, inelastic particle ring. Finally, we introduce an important correction to the model of Chiang and Goldreich.

Mosqueira, I.↗
DOE OSTI2026

FC Site 4.0 - NLR Thermodynamic profiler (ASSIST II-11) Thermodynamic Retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe v0.19 (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NLR ASSIST II infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height (CBH) from co-located scanning lidar. The full pipeline for running the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. Met data was not ingested. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches in Denver, CO.

17 WIND ENERGY↗
DOE OSTI2022

Tallo: A global tree allometry and crown architecture database

Abstract Data capturing multiple axes of tree size and shape, such as a tree's stem diameter, height and crown size, underpin a wide range of ecological research—from developing and testing theory on forest structure and dynamics, to estimating forest carbon stocks and their uncertainties, and integrating remote sensing imagery into forest monitoring programmes. However, these data can be surprisingly hard to come by, particularly for certain regions of the world and for specific taxonomic groups, posing a real barrier to progress in these fields. To overcome this challenge, we developed the Tallo database, a collection of 498,838 georeferenced and taxonomically standardized records of individual trees for which stem diameter, height and/or crown radius have been measured. These data were collected at 61,856 globally distributed sites, spanning all major forested and non‐forested biomes. The majority of trees in the database are identified to species (88%), and collectively Tallo includes data for 5163 species distributed across 1453 genera and 187 plant families. The database is publicly archived under a CC‐BY 4.0 licence and can be access from: https://doi.org/10.5281/zenodo.6637599 . To demonstrate its value, here we present three case studies that highlight how the Tallo database can be used to address a range of theoretical and applied questions in ecology—from testing the predictions of metabolic scaling theory, to exploring the limits of tree allometric plasticity along environmental gradients and modelling global variation in maximum attainable tree height. In doing so, we provide a key resource for field ecologists, remote sensing researchers and the modelling community working together to better understand the role that trees play in regulating the terrestrial carbon cycle.

Jucker, Tommaso↗
NASA NTRS2013

Modeling of 2008 Kasatochi Volcanic Sulfate Direct Radiative Forcing: Assimilation of OMI SO2 Plume Height Data and Comparison with MODIS and CALIOP Observations

Volcanic SO2 column amount and injection height retrieved from the Ozone Monitoring Instrument (OMI) with the Extended Iterative Spectral Fitting (EISF) technique are used to initialize a global chemistry transport model (GEOS-Chem) to simulate the atmospheric transport and lifecycle of volcanic SO2 and sulfate aerosol from the 2008 Kasatochi eruption, and to subsequently estimate the direct shortwave, top-of-the-atmosphere radiative forcing of the volcanic sulfate aerosol. Analysis shows that the integrated use of OMI SO2 plume height in GEOS-Chem yields: (a) good agreement of the temporal evolution of 3-D volcanic sulfate distributions between model simulations and satellite observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Cloud-Aerosol Lidar with Orthogonal Polarisation (CALIOP), and (b) an e-folding time for volcanic SO2 that is consistent with OMI measurements, reflecting SO2 oxidation in the upper troposphere and stratosphere is reliably represented in the model. However, a consistent (approx. 25 %) low bias is found in the GEOS-Chem simulated SO2 burden, and is likely due to a high (approx.20 %) bias of cloud liquid water amount (as compared to the MODIS cloud product) and the resultant stronger SO2 oxidation in the GEOS meteorological data during the first week after eruption when part of SO2 underwent aqueous-phase oxidation in clouds. Radiative transfer calculations show that the forcing by Kasatochi volcanic sulfate aerosol becomes negligible 6 months after the eruption, but its global average over the first month is -1.3W/sq m, with the majority of the forcing-influenced region located north of 20degN, and with daily peak values up to -2W/sq m on days 16-17. Sensitivity experiments show that every 2 km decrease of SO2 injection height in the GEOS-Chem simulations will result in a approx.25% decrease in volcanic sulfate forcing; similar sensitivity but opposite sign also holds for a 0.03 m increase of geometric radius of the volcanic aerosol particles. Both sensitivities highlight the need to characterize the SO2 plume height and aerosol particle size from space. While more research efforts are warranted, this study is among the first to assimilate both satellite-based SO2 plume height and amount into a chemical transport model for an improved simulation of volcanic SO2 and sulfate transport.

Wang, J.↗
DOE OSTI2022

Avalanche gain and its effect on energy resolution in GEM-based detectors

Here, we present avalanche gain and associated resolution measurements recorded with a 4 He:CO 2 (70:30) gas mixture and pure SF 6 , a Negative Ion (NI) gas. SF 6 is of particular interest to the directional dark matter detection community, as its low thermal diffusion helps to retain recoil ionization track features over long drift lengths. With the aid of a general form of the reduced first Townsend coefficient (RFTC), multiple GEM-based detector data sets are used to study the high-gain behavior of the 4 He:CO 2 gas mixture. The high-gain data is well described purely in terms of the reduced electric field strength and the number of GEMs, and the robust relationship between the RFTC and the average, reduced, electric field strength across the GEMs is emphasized. The associated (pulse-height) resolution measurements are used to discuss the variance of the avalanche distribution and to describe and estimate the lower limits of energy resolution one should expect to measure using a simple relationship with the RFTC. In the end, a description of avalanche gain, its effect on energy resolution, and the contributing experimental parameters in GEM-based detectors is developed over a broad parameter space for further use.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Episodic Magma Hammers for the 15 January 2022 Cataclysmic Eruption of Hunga Tonga‐Hunga Ha'apai

Understanding the forces and magma system dynamics on timescales of seconds to minutes remains challenging. In the January 2022 phreatoplinian Hunga Tonga-Hunga Ha'apai eruption, four remarkably similar seismic subevents within a 5-min interval occurred during the intensifying early eruptive phase. The subevents are similar in waveforms and durations (~25 s each). Each subevent begins with an unusual negative P-wave polarity which is inferred, using full-wave seismic modeling, to be caused by an upward single-force mechanism at the volcano created by a magma hammer likely in response to magma flow blockage/constriction during the early part of the eruption as discharge rapidly increased over orders of magnitude with concomitant conduit geometry evolution and instability. Our proposed episodic magma hammer model is consistent with thermodynamic and phase properties of the magmatic mixture, and yields an estimate of conduit mass flow in agreement with vent discharge rates derived from satellite imagery of plume heights.

58 GEOSCIENCES↗