Title: Interplay of Gravity Waves and Disturbance Electric Fields to the Abnormal Ionospheric Variations During the 11 May 2024 Superstorm
Abstract The strongest geomagnetic storm in the preceding two decades occurred in May 2024. Over these years, ground‐based observational capabilities have been significantly enhanced to monitor the ionospheric weather. Notably, the newly established Sanya incoherent scatter radar (SYISR) (Yue, Wan, Ning, & Jin, 2022,https://doi.org/10.1038/s41550‐022‐01684‐1), one of the critical infrastructures of the Chinese “Meridian Project,” provides multiple parameter measurements in the upper atmosphere at low latitudes over Asian longitudies. Unique ionospheric changes on superstorm day 11 May were first recorded by the SYISR experiments and the geostationary satellite (GEO) total electron content (TEC) network over the Asian sector. The electron density or TEC displayed wavelike structures rather than a regular diurnal pattern. Surprisingly, two humps, a common feature in the daytime equatorial ionization anomaly structure, disappeared. The SYISR observations revealed that multiple wind surges accompanied the downward phase propagation caused by atmospheric gravity waves (AGWs) originating from auroral zones. Meanwhile, strong upward and large downward drifts were respectively observed in the daytime and around sunset. The Thermosphere‐Ionosphere Electrodynamics Global Circulation Model (TIEGCM) simulations demonstrated that abnormal ionospheric changes were attributed to meridional wind disturbances associated with AGWs and recurrent penetration electric fields corresponding to largerBzsouthward excursions and disturbance dynamo. The complicated interplay between AGWs and disturbance electric fields contributed to this unique ionospheric variation.  more » « less
Award ID(s):
2033787 2149698 1952737
PAR ID:
10612069
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
AGU Advances
Volume:
6
Issue:
1
ISSN:
2576-604X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We introduce a new framework called Machine Learning (ML) based Auroral Ionospheric electrodynamics Model (ML‐AIM). ML‐AIM solves a current continuity equation by utilizing the ML model of Field Aligned Currents of Kunduri et al. (2020,https://doi.org/10.1029/2020JA027908), the FAC‐derived auroral conductance model of Robinson et al. (2020,https://doi.org/10.1029/2020JA028008), and the solar irradiance conductance model of Moen and Brekke (1993,https://doi.org/10.1029/92gl02109). The ML‐AIM inputs are 60‐min time histories of solar wind plasma, interplanetary magnetic fields (IMF), and geomagnetic indices, and its outputs are ionospheric electric potential, electric fields, Pedersen/Hall currents, and Joule Heating. We conduct two ML‐AIM simulations for a weak geomagnetic activity interval on 14 May 2013 and a geomagnetic storm on 7–8 September 2017. ML‐AIM produces physically accurate ionospheric potential patterns such as the two‐cell convection pattern and the enhancement of electric potentials during active times. The cross polar cap potentials (ΦPC) from ML‐AIM, the Weimer (2005,https://doi.org/10.1029/2004ja010884) model, and the Super Dual Auroral Radar Network (SuperDARN) data‐assimilated potentials, are compared to the ones from 3204 polar crossings of the Defense Meteorological Satellite Program F17 satellite, showing better performance of ML‐AIM than others. ML‐AIM is unique and innovative because it predicts ionospheric responses to the time‐varying solar wind and geomagnetic conditions, while the other traditional empirical models like Weimer (2005,https://doi.org/10.1029/2004ja010884) designed to provide a quasi‐static ionospheric condition under quasi‐steady solar wind/IMF conditions. Plans are underway to improve ML‐AIM performance by including a fully ML network of models of aurora precipitation and ionospheric conductance, targeting its characterization of geomagnetically active times. 
    more » « less
  2. Abstract We investigate Ionosphere‐Thermosphere (IT) responses to the March 2023 geomagnetic storm using GOLD and PFISR observations, along with TIEGCM simulations driven by data‐assimilated aurora and electric fields. A Lattice Kriging approach is implemented to assimilate auroral electron flux and characteristic energy from ground‐based (THEMIS/ASIs) and satellite (DMSP/SSUSIs) measurements, and the empirical model (Zhang & Paxton, 2008,https://doi.org/10.1016/j.jastp.2008.03.008). Electric field assimilation uses SuperDARN ion drifts, PFISR, and the SuperDARN SHF potential model. Compared to the default run, the assimilated simulation better captures GOLD observations, including an O/N2depletion by ∼0.5, and a neutral temperature enhancement of ∼170 K relative to a quiet day. These changes are primarily driven by enhanced Joule heating from assimilated aurora and electric fields, and the resultant upwelling. Results from different model runs implies that aurora modestly increases temperature, while electric fields further enhance it and modulate its spatial structure. The assimilated run also captures the local PFISR observations, including (a) F‐region elevated ion and electron temperatures (∼4000 K); (b) a negative ionospheric storm phase with ∼60% F‐region electron density reduction; and (c) E‐region electron density enhancements by 50%, relative to a quiet day. Electric fields more strongly elevate ion temperature and reduce F‐ region electron density, whereas aurora drives E‐region electron density enhancements. The combined effect of both drivers produces the largest height‐integrated Joule heating at the Poker Flat, reaching ∼120 mW/m2. Overall, we demonstrate that assimilating high‐latitude drivers into TIEGCM significantly improves its capability to simulate storm‐time IT dynamics with richer spatial structures and increased variability. 
    more » « less
  3. Abstract We compared the performance of DREAM3D simulations in reproducing the long‐term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low‐energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics insideL* = 6, with a better model performance at lowerμthan higherμ, whereμis the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data‐driven outer boundary condition is critical to the model performance, while adding in the data‐driven seed population doesn't further improve the performance. (b) The model shows comparable performance withDLLfrom Brautigam and Albert (2000,https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014,https://doi.org/10.1002/2013ja019204), and Liu et al. (2016,https://doi.org/10.1002/2015gl067398), while withDLLfrom Ali et al. (2016,https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data‐based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higherμusing the Wang et al. (2019,https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under‐reproduced in the current model. 
    more » « less
  4. Abstract In their article entitled “Trapdoor Fault Activation:A Step Toward Caldera Collapse at Sierra Negra,Galapagos,Ecuador” Shreve and Delgado (2023,https://doi.org/10.1029/2023jb026437) examine co‐eruptive deformation during the 2018 eruption of Sierra Negra Volcano. One of their major conclusions is that the 2018 eruption, and specifically co‐eruptive faulting, represents the initial stages of caldera collapse. They reach this conclusion because they focus their analysis solely on co‐eruptive deformation, and do not investigate the total (net) deformation for the 2005 to 2018 eruption cycle. Bell, La Famina, et al. (2021,https://doi.org/10.1038/s41467‐021‐21596‐4) investigated both the pre‐ and co‐eruptive phases of the 2018 eruption and showed that net deformation was one of caldera resurgence, not subsidence. In this comment, we demonstrate that the conclusion of collapse, or even initiation of collapse, is attributable to not accounting for pre‐eruptive deformation on the intra‐caldera Trapdoor Fault system and incorrectly assuming that the volcano‐tectonic dynamics of Sierra Negra mimic those of other basaltic calderas. 
    more » « less
  5. Abstract Most ionospheric models cannot sufficiently reproduce the observed electron density profiles in the E‐region ionosphere, since they usually underestimate electron densities and do not match the profile shape. Mitigation of these issues is often addressed by increasing the solar soft X‐ray flux which is ineffective for resolving data‐model discrepancies. We show that low‐resolution cross sections and solar spectral irradiances fail to preserve structure within the data, which considerably impacts radiative processes in the E‐region, and are largely responsible for the discrepancies between observations and simulations. To resolve data‐model inconsistencies, we utilize new high‐resolution (0.001 nm) atomic oxygen (O) and molecular nitrogen (N2) cross sections and solar spectral irradiances, which contain autoionization and narrow rotational lines that allow solar photons to reach lower altitudes and increase the photoelectron flux. This work improves upon Meier et al. (2007,https://doi.org/10.1029/2006gl028484) by additionally incorporating high‐resolution N2photoionization and photoabsorption cross sections in model calculations. Model results with the new inputs show increased O+production rates of over 500%, larger than those of Meier et al. (2007,https://doi.org/10.1029/2006gl028484) and total ion production rates of over 125%, while production rates decrease by ∼15% in the E‐region in comparison to the results obtained using the cross section compilation from Conway (1988,https://apps.dtic.mil/sti/pdfs/ADA193866.pdf). Low‐resolution molecular oxygen (O2) cross sections from the Conway compilation are utilized for all input cases and indicate that is a dominant contributor to the total ion production rate in the E‐region. Specifically, the photoionization contributed from longer wavelengths is a main contributor at ∼120 km. 
    more » « less