Volcanic eruption source parameters may be estimated from acoustic pressure recordings dominant at infrasonic frequencies (< 20 Hz), yet uncertainties may be high due in part to poorly understood propagation dynamics. Linear acoustic propagation of volcano infrasound is commonly assumed, but nonlinear processes such as wave steepening may distort waveforms and obscure the sourcing process in recorded waveforms. Here we use a previously developed frequency-domain nonlinearity indicator to quantify spectral changes due to nonlinear propagation primarily in 80 signals from explosions at Yasur Volcano, Vanuatu. We find evidence for
Sound waves generated by erupting volcanoes can be used to infer important source dynamics, yet acoustic source‐time functions may be distorted during propagation, even at local recording distances (
- NSF-PAR ID:
- 10449120
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 125
- Issue:
- 3
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract $$\le$$ -
more » « less
Abstract The biggest volcanic eruption since 1991 happened on 15 January 2022 on the island of Hunga Tonga‐Hunga Haʻapai (20.6°S; 175.4°W) in the South Pacific between 4:00 and 4:16 UT. The updrafts from the eruption reached 58 km height. In order to observe its ionospheric effects, approximately 750 GNSS receivers in New Zealand and Australia were used to calculate the detrended total electron content (dTEC). Traveling ionospheric disturbances (TIDs) were observed over New Zealand 1.0–1.5 hr after the volcano eruption, with a horizontal wavelength (
) of 1,525 km, horizontal phase velocity ( ) of 635 m/s, period ( τ ) of 40 min, and azimuth (α ) of 214°. On the other hand, TIDs were observed 2–3 hr after the eruption in Australia with, , τ , andα of 922 km, 375 m/s, 41 min, and 266°, respectively. Using reverse ray tracing, we found that these GWs originated atz > 100 km at a location ∼500 km south of Tonga, in agreement with model results for the location of a large amplitude body force created from the breaking of primary GWs from the eruption. Thus, we found that these fast GWs were secondary, not primary GWs from the Tonga eruption. -
more » « less
Abstract Infrasound observations are increasingly used to constrain properties of volcanic eruptions. In order to better interpret infrasound observations, however, there is a need to better understand the relationship between eruption properties and sound generation. Here we perform two‐dimensional computational aeroacoustic simulations where we solve the compressible Navier‐Stokes equations for pure‐air with a large‐eddy simulation approximation. We simulate idealized impulsive volcanic eruptions where the exit velocity is specified and the eruption is pressure‐balanced with the atmosphere. Our nonlinear simulation results are compared with the commonly used analytical linear acoustics model of a compact monopole source radiating acoustic waves isotropically in a half space. The monopole source model matches the simulations for low exit velocities (
m/s or where is the Mach number); however, the two solutions diverge as the exit velocity increases with the simulations developing lower peak amplitude, more rapid onset, and anisotropic radiation with stronger infrasound signals recorded above the vent than on Earth's surface. Our simulations show that interpreting ground‐based infrasound observations with the monopole source model can result in an underestimation of the erupted volume for eruptions with sonic or supersonic exit velocities. We examine nonlinear effects and show that nonlinear effects during propagation are relatively minor for the parameters considered. Instead, the dominant nonlinear effect is advection by the complex flow structure that develops above the vent. This work demonstrates the need to consider anisotropic radiation patterns and jet dynamics when interpreting infrasound observations, particularly for eruptions with sonic or supersonic exit velocities. -
more » « less
Abstract The Whitham equation was proposed as a model for surface water waves that combines the quadratic flux nonlinearity
of the Korteweg–de Vries equation and the full linear dispersion relation of unidirectional gravity water waves in suitably scaled variables. This paper proposes and analyzes a generalization of Whitham's model to unidirectional nonlinear wave equations consisting of a general nonlinear flux function and a general linear dispersion relation . Assuming the existence of periodic traveling wave solutions to this generalized Whitham equation, their slow modulations are studied in the context of Whitham modulation theory. A multiple scales calculation yields the modulation equations, a system of three conservation laws that describe the slow evolution of the periodic traveling wave's wavenumber, amplitude, and mean. In the weakly nonlinear limit, explicit, simple criteria in terms of general and establishing the strict hyperbolicity and genuine nonlinearity of the modulation equations are determined. This result is interpreted as a generalized Lighthill–Whitham criterion for modulational instability. -
more » « less
Abstract A number of interdependent conditions and processes contribute to ionospheric‐origin energetic (
10 eV to several keV) ion outflows. Due to these interdependences and the associated observational challenges, energetic ion outflows remain a poorly understood facet of atmosphere‐ionosphere‐magnetosphere coupling. Here we demonstrate the relationship between east‐west magnetic field fluctuations ( ) and energetic outflows in the magnetosphere‐ionosphere transition region. We use dayside cusp region FAST satellite observations made near apogee ( 4,180‐km altitude) near fall equinox and solstices in both hemispheres to derive statistical relationships between ion upflow and spectral power as a function of spacecraft frame frequency bands between 0 and 4 Hz. Identification of ionospheric‐origin energetic ion upflows is automated, and the spectral power in each frequency band is obtained via integration of power spectral density. Derived relationships are of the form for upward ion flux at 130‐km altitude, with the mapped upward ion flux for a nominal spectral power nT . The highest correlation coefficients are obtained for spacecraft frame frequencies 0.1–0.5 Hz. Summer solstice and fall equinox observations yield power law indices 0.9–1.3 and correlation coefficients , while winter solstice observations yield 0.4–0.8 with . Mass spectrometer observations reveal that the oxygen/hydrogen ion composition ratio near summer solstice is much greater than the corresponding ratio near winter. These results reinforce the importance of ion composition in outflow models. If observed perturbations result from Doppler‐shifted wave structures with near‐zero frequencies, we show that spacecraft frame frequencies 0.1–0.5 Hz correspond to perpendicular spatial scales of several to tens of kilometers.
