skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Ionospheric Disturbances Generated by the 2015 Calbuco Eruption: Comparison of GITM‐R Simulations and GNSS Observations
Abstract Volcanic eruptions provide broad spectral forcing to the atmosphere and understanding the primary mechanisms that are relevant to explain the variety in waveform characteristics in the Ionosphere‐Thermosphere (IT) is still an important open question for the community. In this study, Global Navigation Satellite System (GNSS) Total Electron Content (TEC) data are analyzed and compared to simulations performed by the Global Ionosphere‐Thermosphere Model with Local Mesh Refinement (GITM‐R) for the first phase of the 2015 Calbuco eruption that occurred on 22 April. A simplified source representation and spectral acoustic‐gravity wave (AGW) propagation model are used to specify the perturbation at the lower boundary of GITM‐R at 100 km altitude. Two assumptions on the propagation structure, Direct Spherical (DS) and Ground Coupled (GC), are compared to the GNSS data and these modeling specifications show good agreement with different aspects of the observations for some waveform characteristics. Most notably, GITM‐R is able to reproduce the relative wave amplitude of AGWs as a function of radial distance from the vent, showing acoustic dominant forcing in the near field (<500 km) and gravity dominant forcing in the far‐field (>500 km). The estimated apparent phase speeds from GITM‐R simulations are consistent with observations with ∼10% difference from observation for both acoustic wave packets and a trailing gravity mode. The relevance of the simplifications made in the lower atmosphere to the simulated IT response is then discussed.  more » « less
Award ID(s):
2149698 2033787 1952737
PAR ID:
10511945
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Space Weather
Volume:
22
Issue:
2
ISSN:
1542-7390
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The response of the thermosphere ionosphere system to an X1.3 class solar flare is studied using observations of the total electron content (TEC) and the Global Ionosphere Thermosphere Model (GITM) simulations. The solar flare erupted from the active region AR12975 on 30 March 2022. Owing to the absence of accompanying severe geomagnetic activity, it was possible to isolate the effects of the flare on the upper atmosphere. TEC data are processed for Continental USA (CONUS), employing filtering and binning techniques to create 2D variation maps. The spectral content of the TEC variations is analyzed using a wavelet coherence method. The immediate response of the solar flare exhibited broad similarities, while notable differences were observed during the recovery period between the East and West sides of the CONUS. GITM is used to explore the East–West asymmetry of the key T‐I parameters. Simulation results reveal that the coinciding interplanetary magnetic field southward turning had a greater influence on these parameters compared to the solar flare, while their nonlinear interaction introduced complex variations. Additional investigation reveals gravity wave damping also contributes to the asymmetric solar flare response. 
    more » « less
  2. The Tonga volcano eruption at 04:14:45 UT on 2022-01-15 released enormous amounts of energy into the atmosphere, triggering very significant geophysical variations not only in the immediate proximity of the epicenter but also globally across the whole atmosphere. This study provides a global picture of ionospheric disturbances over an extended period for at least 4 days. We find traveling ionospheric disturbances (TIDs) radially outbound and inbound along entire Great-Circle loci at primary speeds of ∼300–350 m/s (depending on the propagation direction) and 500–1,000 km horizontal wavelength for front shocks, going around the globe for three times, passing six times over the continental US in 100 h since the eruption. TIDs following the shock fronts developed for ∼8 h with 10–30 min predominant periods in near- and far- fields. TID global propagation is consistent with the effect of Lamb waves which travel at the speed of sound. Although these oscillations are often confined to the troposphere, Lamb wave energy is known to leak into the thermosphere through channels such as atmospheric resonance at acoustic and gravity wave frequencies, carrying substantial wave amplitudes at high altitudes. Prevailing Lamb waves have been reported in the literature as atmospheric responses to the gigantic Krakatoa eruption in 1883 and other geohazards. This study provides substantial first evidence of their long-duration imprints up in the global ionosphere. This study was enabled by ionospheric measurements from 5,000+ world-wide Global Navigation Satellite System (GNSS) ground receivers, demonstrating the broad implication of the ionosphere measurement as a sensitive detector for atmospheric waves and geophysical disturbances. 
    more » « less
  3. Abstract During the second recovery phase of the 13–14 March 2022 storm, intense high‐latitude neutral mass density spikes are detected by satellites at ∼500 km in both hemispheres. These density spikes, accurately modeled by the Global Ionospheric Thermosphere Model (GITM), are identified as high‐latitude neutral mass density anomalies (HDAs). The GITM simulation indicates that these HDAs, which extends over the polar region, are important features in high‐latitude neutral density. Furthermore, GITM reveals that these HDAs are manifestations of transpolar traveling atmospheric disturbances triggered on the dawn side. Moreover, GITM also reveals significant interhemispheric asymmetries (IHAs) in the magnitude, propagation speed, and propagation direction of HDAs, which are linked to the IHAs in the distribution and magnitude of Joule heating deposited as well as the thermospheric background conditions. This study presents a dynamic perspective on the IHA of storm‐time high‐latitude neutral density variations that is particularly helpful to the proper interpretation of satellite observations. 
    more » « less
  4. Abstract In this paper, we simulate an observed mountain wave event over central Europe and investigate the subsequent generation, propagation, phase speeds and spatial scales, and momentum deposition of secondary waves under three different tidal wind conditions. We find the mountain wave breaks just below the lowest critical level in the mesosphere. As the mountain wave breaks, it extends outwards along the phases and fluid associated with the breaking flows downstream of its original location by 500–1,000 km. The breaking generates a broad range of secondary waves with horizontal scales ranging from the mountain wave instability scales (20–300 km), to multiples of the mountain wave packet scale (420 km+) and phase speeds from 40 to 150 m/s in the lower thermosphere. The secondary wave morphology consists of semi‐concentric patterns with wave propagation generally opposing the local tidal winds in the mesosphere. Shears in the tidal winds cause breaking of the secondary waves and local wave forcing which generates even more secondary waves. The tidal winds also influence the dominant wavelengths and phase speeds of secondary waves that reach the thermosphere. The secondary waves that reach the thermosphere deposit their energy and momentum over a broad area of the thermosphere, mostly eastward of the source and concentrated between 110 and 130 km altitude. The secondary wave forcing is significant and will likely be very important for the dynamics of the thermosphere. A large portion of this forcing comes from nonlinearly generated secondary waves at relatively small‐scales which arise from the wave breaking processes. 
    more » « less
  5. This study explores the meteorological source and vertical propagation of gravity waves (GWs) that drive daytime traveling ionospheric disturbances (TIDs), using the specified dynamics version of the SD-WACCM-X (Whole Atmosphere Community Climate Model with thermosphere-ionosphere eXtension) and the SAMI3 (Sami3 is Also a Model of the Ionosphere) simulations driven by SD-WACCM-X neutral wind and composition. A cold weather front moved over the northern-central USA (90–100°W, 35–45°N) during the daytime of 20 October 2020, with strong upward airflow. GWs with ~500–700 km horizontal wavelengths propagated southward and northward in the thermosphere over the north-central USA. Also, the perturbations were coherent from the surface to the thermosphere; therefore, the GWs were likely generated by vertical acceleration associated with the cold front over Minnesota and South Dakota. The convectively generated GWs had almost infinite vertical wavelength below ~100 km due to being evanescent. This implies that the GWs tunneled through their evanescent region in the middle atmosphere (where a squared vertical wavenumber is equal to or smaller than 0) and became freely propagating in the thermosphere and ionosphere. Medium-scale TIDs (MSTIDs) also propagated southward with the GWs, suggesting that the convectively generated GWs created MSTIDs. 
    more » « less