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 We study the dynamical and thermal roles of internal gravity waves generated in the troposphere and above using the Coupled Middle Atmosphere Thermosphere‐2 General Circulation Model. This model incorporates the whole atmosphere nonlinear spectral gravity wave parameterization and its extension to include non‐tropospheric sources. We conducted model experiments for northern summer solstice conditions, first including only tropospheric sources, then including sources localized at 50 and 90 km, and uniformly distributed over all heights. The simulated differences in mean temperature and horizontal winds demonstrate that gravity waves produce the greatest dynamical and thermal changes in the latter case compared to the localized sources. While the gravity wave drag is longitudinally uniform in the lower thermosphere, it is more localized in the upper thermosphere in all the simulations. Waves from uniformly distributed sources increase the longitudinal variability of zonal winds in the thermosphere up to 150 km. Gravity wave effects exhibit different local time variations in the lower thermosphere (100–140 km) than in the upper thermosphere. In the upper thermosphere, gravity wave effects are stronger during the day than at night. In contrast, nighttime gravity wave effects are stronger than the daytime ones in the lower thermosphere. Finally, a comparison with ICON‐MIGHTI observations shows that the model reproduces the basic structure of thermospheric winds, performing better with zonal winds than with meridional winds. Adding non‐tropospheric wave sources modifies wind structures in wave‐breaking regions, but does not improve the global statistical comparison. 
    more » « less
  4. 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
  5. 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