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
Ionospheric Disturbances Generated by the 2015 Calbuco Eruption: Comparison of GITM‐R Simulations and GNSS Observations
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- PAR ID:
- 10511945
- 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
-
-
more » « less
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
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
Abstract The global 3‐dimensional structure of the concentric traveling ionospheric disturbances (CTIDs) triggered by 2022 Tonga volcano was reconstructed by using the 3‐dimensional computerized ionospheric tomography (3DCIT) technique and extensive global navigation satellite system (GNSS) observations. This study provides the first estimation of the CTIDs vertical wavelengths, ∼736 km, which was much larger than the gravity wave (GW) vertical wavelength, 240–400 km, estimated using ICON neutral wind observations. Notable trend with the variation of azimuth was also found in horizontal speeds at 200 and 500 km altitudes and differences between them. These results imply that (a) the global propagation of Lamb waves determined the arrival time of local ionospheric disturbances, and (b) the arriving Lamb waves caused vertical atmospheric perturbations that are not typical of GWs, resulting in local thermospheric horizontal wave propagation which is faster than the Lamb wave propagation at lower altitudes.
-
more » « less
We present the first simultaneous observations of a traveling ionosphere wave (TID) event, measuring electron concentration (Ne), vertical plasma drift (Vz), and ion and electron temperatures (Ti, Te) using the Arecibo incoherent scatter radar. A TID with a period of 135 min was evident in all four state variables in the thermosphere. The amplitudes of Vz and relative Ti fluctuations show only small height variations from 200 to 500 km and their vertical wavelengths increase with altitude. The Te fluctuation shows different characteristics from EISCAT in both phase and amplitude. When the geomagnetic dip angle is 45°, half of the driving gravity wave’s (GW’s) equatorward velocity is mapped to Vz. This meridional-to-vertical velocity coupling amplifies GW’s effect in Ne through vertical transport. The amplifying and anisotropic effects of the geomagnetic field explain the ubiquitous presence of TIDs and their preferred equatorward propagation direction in the geomagnetic mid-latitudes, as well as the midnight collapse phenomenon observed at Arecibo.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1029/2023SW003502
