skip to main content

Title: The Response of Middle Thermosphere (∼160 km) Composition to the November 20 and 21, 2003 Superstorm

TIMED/GUVI limb measurements and first‐principles simulations from the Thermosphere Ionosphere Electrodynamics Global Circulation Model (TIEGCM) are used to investigate thermospheric atomic oxygen (O) and molecular nitrogen (N2) responses in the middle thermosphere on a constant pressure surface (∼160 km) to the November 20 and 21, 2003 superstorm. The consistency between GUVI observations and TIEGCM simulated composition changes allows us to utilize TIEGCM outputs to investigate the storm‐time behaviors of O and N2systematically. Diagnostic analysis shows that horizontal and vertical advection are the two main processes that determine the storm‐induced perturbations in the middle thermosphere. Molecular diffusion has a relatively smaller magnitude than the two advection processes, acting to compensate for the changes caused by the transport partly. Contributions from chemistry and eddy diffusion are negligible. During the storm initial and main phases, composition variations at high latitudes are determined by both horizontal and vertical advection. At middle‐low latitudes, horizontal advection is the main driver for the composition changes where O mass mixing ratiodecreases (N2mass mixing ratioincreases); whereas horizontal and vertical advection combined to dominate the changes in the regions whereincreases (decreases). Over the entire storm period, horizontal advection plays a significant role in transporting high‐latitude composition perturbations globally. Our results also demonstrate that storm‐time temperature changes are not the direct cause of the composition perturbations on constant pressure surfaces.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Space Physics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    A thermospheric O and N2column density ratio (∑O/N2) depletion with long‐duration (>16 hr) was observed by the Global‐scale Observations of the Limb and Disk at the Atlantic longitudes (75W–20W) and middle latitudes (20N–50N) during the recovery phase of the 8 June 2019 geomagnetic storm. The National Center for Atmospheric Research Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM) simulations reproduced the ∑O/N2depletion patterns with a similar magnitude, and indicated that the composition recovery at middle latitudes began several hours after the beginning of the recovery phase of the geomagnetic storm. The TIEGCM simulations enable quantitative analysis of the physical mechanisms driving the middle‐latitude composition changes during the storm recovery phase. This analysis indicates that vertical advection and molecular diffusion dominated the initial recovery of composition perturbations at middle latitudes. Horizontal advection was also a main driver in the initial recovery of composition, but its contribution decreased rapidly. In the late recovery phase, the composition recovery was mainly determined by horizontal advection. In comparison, vertical advection and molecular diffusion played a much less important role.

    more » « less
  2. Abstract

    In this study, we analyze the thermospheric density data provided by the Gravity Field and Steady‐State Ocean Circulation Explorer during June–August 2010–2013 at ∼260 km altitude and the Challenging Minisatellite Payload during June–August 2004–2007 at ∼370 km altitude to study high latitude traveling atmospheric disturbances (TADs) in austral winter. We extract the TADs along the satellite tracks from the density for varyingKp, and linearly extrapolate the TAD distribution toKp = 0; we call these the geomagnetic “quiet time” results here. We find that the quiet time spatial distribution of TADs depends on the spatial scale (along‐track horizontal wavelength) and altitude. Atz∼ 260 km, TADs with ≤ 330 km are seen mainly around and slightly downstream of the Southern Andes‐Antarctic region, while TADs with > 800 km are distributed fairly evenly around the geographic South pole at latitudes ≥60°S. Atz∼ 370 km, TADs with ≤ 330 km are relatively weak and are distributed fairly evenly over Antarctica, while TADs with > 330 km make up a bipolar distribution. For the latter, the larger size lobe is centered at ∼60°S, and is located around, downstream and somewhat upstream of the Andes/Antarctic Peninsula, while the smaller lobe is located over the Antarctic continent at 90°–150°E. We also find that the TAD morphology forKp ≥ 2 and > 330 km depends strongly on geomagnetic activity, likely due to auroral activity, with greatly enhanced TAD amplitudes with increasingKp.

    more » « less
  3. 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
  4. Abstract

    Based on the observations from the balloon‐borne instrument High‐altitude Interferometer WIND experiment (HIWIND) and the simulations from the Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM), the Grid Agnostic MHD Environment for Research Applications (GAMERA)‐TIEGCM (GT), and the GAMERA‐TIEGCM‐RCM (GTR), we investigate the variations of summer high‐latitude thermospheric winds and their physical mechanisms from 25 to 30 June, 2018. HIWIND observations show that the meridional winds were the largest at midnight and exhibited strong day‐to‐day variations during the 6‐day period, which were generally reproduced by those three models. The day‐to‐day variations of winds were mainly associated with the interplanetary magnetic field (IMF)perturbations, while the magnetic latitude variations also contributed to the large day‐to‐day variations of the winds seen in the observations. Meanwhile, the zonal winds were mostly westward during the daytime, and the wind speed became large, especially in the afternoon, which is related to the westward ion drift velocity. The observed meridional winds tend to turn equatorward during the daytime on some days, while the simulated winds blow mostly poleward except for simulations by the GTR model on 26 June. The GTR model revealed that the equatorward meridional winds on 26 June were associated with strong and negative IMFconditions, which tilts the convection pattern to the prenoon sector. The simulations also revealed that the ring current could contribute to affecting the neutral wind variations, especially under geomagnetically active conditions.

    more » « less
  5. Abstract

    We present a statistical investigation of the effects of interplanetary magnetic field (IMF) on hemispheric asymmetry in auroral currents. Nearly 6 years of magnetic field measurements from Swarm A and C satellites are analyzed. Bootstrap resampling is used to remove the difference in the number of samples and IMF conditions between the local seasons and the hemispheres. Currents are stronger in Northern Hemisphere (NH) than Southern Hemisphere (SH) for IMF Bin NH (Bin SH) in most local seasons under both signs of IMF B. For Bin NH (Bin SH), the hemispheric difference in currents is small except in local winter when currents in NH are stronger than in SH. During Band Bin NH (Band Bin SH), the largest hemispheric asymmetry occurs in local winter and autumn, when the NH/SH ratio of field aligned current (FAC) is 1.180.09 in winter and 1.170.09 in autumn. During Band Bin NH (Band Bin SH), the largest asymmetry is observed in local autumn with NH/SH ratio of 1.160.07 for FAC. We also find an explicit Beffect on auroral currents in a given hemisphere: on average Bin NH and Bin SH causes larger currents than vice versa. The explicit Beffect on divergence‐free current during IMF Bis in very good agreement with the Beffect on the cross polar cap potential from the Super Dual Auroral Radar Network dynamic model except at SH equinox and NH summer.

    more » « less