Primary Versus Secondary Gravity Wave Responses at F-Region Heights Generated by a Convective Source

A 2D nonlinear, compressible model is used to simulate the acoustic-gravity wave (AGW, i.e., encompassing the spectrum of acoustic and gravity waves) response to a thunderstorm squall-line type source. We investigate the primary and secondary neutral AGW response in the thermosphere, consistent with waves that can couple to the F-region ionospheric plasma, and manifest as Traveling Ionospheric Disturbances (TIDs). We find that primary waves at z = 240 km altitude have wavelengths and phase speeds in the range 170–270 km, and 180–320 m/s, respectively. The secondary waves generated have wavelengths ranging from ∼100 to 600 km, and phase speeds from 300 to 630 m/s. While there is overlap in the wave spectra, we find that the secondary waves (i.e., those that have been nonlinearly transformed or generated secondarily/subsequently from the primary wave) generally have faster phases than the primary waves. We also assess the notion that waves with fast phase speeds (that exceed proposed theoretical upper limits on passing from the mesosphere to thermosphere) observed at F-region heights must be secondary waves, for example, those generated in situ by wave breaking in the lower thermosphere, rather than directly propagating primary waves from their sources. We find that primary waves with phase speeds greater than this proposed upper limit can tunnel through a deep portion of the lower/middle atmosphere and emerge as propagating waves in the thermosphere. Therefore, comparing a TID\u27s/GWs phase speed with this upper limit is not a robust method of identifying whether an observed TID originates from a primary versus secondary AGW

Simulation of Infrasonic Acoustic Wave Imprints on Airglow Layers During the 2016 M7.8 Kaikoura Earthquake

Simulations of hypothesized but unobserved mesopause airglow (MA) disturbances generated by infrasonic acoustic waves (IAWs) during the 2016 M7.8 Kaikoura earthquake are performed. Realistic surface displacements are calculated in a forward seismic wave propagation model and incorporated into a 3-D nonlinear compressible neutral atmosphere model as a source of IAWs at the surface-air interface. Inchin et al. (2021), https://doi.org/10.1029/2020av000260 previously showed that Global Positioning System-based total electron content (TEC) observations can be used to constrain the finite-fault kinematics of the Kaikoura earthquake. However, due to limitations of Global Navigation Satellite System network coverage and coalescence of nonlinear IAW fronts, they pointed to the relative insensitivity of the observed near-zenith TEC perturbations to the rupture evolution on the Papatea fault (PF). Here, we demonstrate that MA observations may have been able to supplement the investigation of the PF, providing information on both the timing of rupture initiation and its direction of propagation. The amplitudes of perturbations of vertically integrated volume emission rates for the simulated hydroxyl (OH)(3,1) and atomic oxygen O( 1S) 557.7 nm reach ∼18% peak-to-peak, and ∼3.2% (5.8 K) peak-to-peak perturbations in OH(3,1) temperature. Our results suggest that observations of nighttime MA imprints of coseismic IAWs are feasible with ground-based imagers, and may supplement the study of finite-fault kinematics of large crustal earthquakes

Table: Spatial Extents of the Numerical Domains

The spatial extents of the numerical domains

Inferring the Evolution of a Large Earthquake from Its Acoustic Impacts on the Ionosphere

We investigate the possibility to constrain the evolution of the 2016 M7.8 Kaikoura earthquake evolution based on Global Positioning System signal-derived ionospheric total electron content (TEC) perturbations, that represent plasma responses to infrasonic acoustic waves (IAWs) generated by surface motion. This earthquake exhibited unusual complexity and some first-order aspects of its evolution remain unclear; for example, how and when the Papatea fault (PF) and the corresponding large surface deformation occurred. For various earthquake models, a seismic wave propagation code is used to simulate time-dependent surface deformations, which then excite IAWs in a 3D compressible nonlinear atmospheric model, coupled with a 2D nonlinear multispecies ionospheric plasma dynamic model. Our preferred finite-fault model reproduces the amplitudes, shapes, and time epochs of appearance of detected TEC perturbations well. Additionally, the incorporation of the PF, ruptured during the earthquake, results in the closest agreement between simulated and observed near-zenith vertical TEC perturbations, whereas its absence shows significant discrepancy. This supports the hypothesis that the PF was ruptured during the Kaikoura earthquake. Furthermore, the IAWs and resulting ionospheric plasma disturbances contain additional information on the PF rupture progression, including the timing of initiation and propagation direction, indicating new opportunities to further constrain the PF rupture with low elevation angle “slant” TEC data. The results confirm the ability for TEC measurements to constrain evolutions of large crustal earthquakes to provide new insight beyond traditional seismic and geodetic data sets

The Dynamics of Nonlinear Atmospheric Acoustic-Gravity Waves Generated by Tsunamis Over Realistic Bathymetry

The investigation of atmospheric tsunamigenic acoustic and gravity wave (TAGW) dynamics, from the ocean surface to the thermosphere, is performed through the numerical computations of the 3D compressible nonlinear Navier-Stokes equations. Tsunami propagation is first simulated using a nonlinear shallow water model, which incorporates instantaneous or temporal evolutions of initial tsunami distributions (ITD). Ocean surface dynamics are then imposed as a boundary condition to excite TAGWs into the atmosphere from the ground level. We perform a case study of a large tsunami associated with the 2011 M9.1 Tohuku-Oki earthquake and parametric studies with simplified and demonstrative bathymetry and ITD. Our results demonstrate that TAGW propagation, controlled by the atmospheric state, can evolve nonlinearly and lead to wave self-acceleration effects and instabilities, followed by the excitation of secondary acoustic and gravity waves (SAGWs), spanning a broad frequency range. The variations of the ocean depth result in a change of tsunami characteristics and subsequent tilt of the TAGW packet, as the wave\u27s intrinsic frequency spectrum is varied. In addition, focusing of tsunamis and their interactions with seamounts and islands may result in localized enhancements of TAGWs, which further indicates the crucial role of bathymetry variations. Along with SAGWs, leading long-period phases of the TAGW packet propagate ahead of the tsunami wavefront and thus can be observed prior to the tsunami arrival. Our modeling results suggest that TAGWs from large tsunamis can drive detectable and quantifiable perturbations in the upper atmosphere under a wide range of scenarios and uncover new challenges and opportunities for their observations

Mesopause Airglow Disturbances Driven by Nonlinear Infrasonic Acoustic Waves Generated by Large Earthquakes

Near-epicentral mesopause airglow perturbations, driven by infrasonic acoustic waves (AWs) during a nighttime analog of the 2011 M9.1 Tohoku-Oki earthquake, are simulated through the direct numerical computation of the 3D nonlinear Navier-Stokes equations. Surface dynamics from a forward seismic wave propagation simulation, initialized with a kinematic slip model and performed with the SPECFEM3D_GLOBE model, are used to excite AWs into the atmosphere from ground level. Simulated mesopause airglow perturbations include steep oscillations and persistent nonlinear depletions up to 50% and 70% from the background state, respectively, for the hydroxyl OH(3,1) and oxygen O(1S) 557.7-nm emissions. Results suggest that AWs excited near a large earthquake\u27s epicenter may be strong enough to drive fluctuations in mesopause airglow, some which may persist after the AWs have passed, that could be readily detectable with ground- and/or satellite-based imagers. Synthetic data demonstrate that future airglow observations may be used for the characterization of earthquake mechanisms and surface seismic waves propagation, potentially complementing tsunami early-warning systems based on total electron content (TEC) observations

Modeling of Ionospheric Responses to Atmospheric Acoustic and Gravity Waves Driven by the 2015 Nepal M w 7.8 Gorkha Earthquake

Near- and far-field ionospheric responses to atmospheric acoustic and gravity waves (AGWs) generated by surface displacements during the 2015 Nepal 7.8 Gorkha earthquake are simulated. Realistic surface displacements driven by the earthquake are calculated in three-dimensional forward seismic waves propagation simulation, based on kinematic slip model. They are used to excite AGWs at ground level in the direct numerical simulation of three-dimensional nonlinear compressible Navier-Stokes equations with neutral atmosphere model, which is coupled with a two-dimensional nonlinear multifluid electrodynamic ionospheric model. The importance of incorporating earthquake rupture kinematics for the simulation of realistic coseismic ionospheric disturbances (CIDs) is demonstrated and the possibility of describing faulting mechanisms and surface deformations based on ionospheric observations is discussed in details. Simulation results at the near-epicentral region are comparable with total electron content (TEC) observations in periods ( 3.3 and 6-10 min for acoustic and gravity waves, respectively), propagation velocities ( 0.92 km/s for acoustic waves) and amplitudes (up to 2 TECu). Simulated far-field CIDs correspond to long-period ( 4 mHz) Rayleigh waves (RWs), propagating with the same phase velocity of 4 km/s. The characteristics of modeled RW-related ionospheric disturbances differ from previously-reported observations based on TEC data; possible reasons for these differences are discussed

Comparison of vTEC perturbations from simulations without coseismic rupturing process on Papatea fault

Based on separate simulations, we investigate the sensitivity of vTEC perturbations to the background ionospheric plasma state. The Figure supports the discussion of simulation uncertainties in Section 5.1 of the manuscript

Multi-Layer Evolution of Acoustic-Gravity Waves and Ionospheric Disturbances Over the United States After the 2022 Hunga Tonga Volcano Eruption

e Hunga-Tonga Hunga-Ha\u27apai volcano underwent a series of large-magnitude eruptions that generated in the atmosphere. We investigate the spatial and temporal evolutions of fluctuations driven by atmospheric acoustic-gravity waves (AGWs) and, in particular, the Lamb wave modes in high spatial resolution data sets measured over the Continental United States (CONUS), complemented with data over the Americas and the Pacific. Along with \u3e800 barometer sites, tropospheric observations, and Total Electron Content data from \u3e3,000 receivers, we report detections of volcano-induced AGWs in mesopause and ionosphere-thermosphere airglow imagery and Fabry-Perot interferometry. We also report unique AGW signatures in the ionospheric D-region, measured using Long-Range Navigation pulsed low-frequency transmitter signals. Although we observed fluctuations over a wide range of periods and speeds, we identify Lamb wave modes exhibiting 295–345 m s−1 phase front velocities with correlated spatial variability of their amplitudes from the Earth\u27s surface to the ionosphere. Results suggest that the Lamb wave modes, tracked by our ray-tracing modeling results, were accompanied by deep fluctuation fields coupled throughout the atmosphere, and were all largely consistent in arrival times with the sequence of eruptions over 8 hr. The ray results also highlight the importance of winds in reducing wave amplitudes at CONUS midlatitudes. The ability to identify and interpret Lamb wave modes and accompanying fluctuations on the basis of arrival times and speeds, despite complexity in their spectra and modulations by the inhomogeneous here, suggests opportunities for analysis and modeling to understand their signals to constrain features of azardous events

The details of forward seismic wave simulation method

The description of forward seismic wave simulation methodology and a comparison of observed and simulated vertical displacements