It is well-known that equatorial plasma bubbles (EPBs) are highly correlated to the post-sunset rise of the ionosphere on a climatological basis. However, when proceeding to the daily EPB development, what controls the day-to-day/longitudinal variability of EPBs remains a puzzle. In this study, we investigate the underlying physics responsible for the day-to-day/longitudinal variability of EPBs using the Sami3 is A Model of the Ionosphere (SAMI3) and the Whole Atmosphere Community Climate Model with thermosphere-ionosphere eXtension (WACCM-X). Simulation results on October 20, 22, and 24, 2020 were presented. SAMI3/WACCM-X self-consistently generated midnight EPBs on October 20 and 24, displaying irregular and regular spatial distributions, respectively. However, EPBs are absent on October 22. We investigate the role of gravity waves on upwelling growth and EPB development and discuss how gravity waves contribute to the distributions of EPBs. Of particular significance is that we found the westward wind associated with solar terminator waves and gravity waves causes midnight vertical drift enhancement and collisional shear instability, which provides conditions favorable for upwelling growth and EPB development. The converging and diverging winds associated with solar terminator waves and midnight temperature maximum also affect the longitudinal distribution of EPBs. The absence of EPBs on October 22 is related to the weak upward drift induced by weak westward wind associated with solar terminator waves.

Assessing space weather modeling capability is a key element in improving existing models and developing new ones. In order to track improvement of the models and investigate impacts of forcing, from the lower atmosphere below and from the magnetosphere above, on the performance of ionosphere-thermosphere models, we expand our previous assessment for 2013 March storm event [Shim et al., 2018]. In this study, we evaluate new simulations from upgraded models (Coupled Thermosphere Ionosphere Plasmasphere Electrodynamics (CTIPe) model version 4.1 and Global Ionosphere Thermosphere Model (GITM) version 21.11) and from NCAR Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension (WACCM-X) version 2.2 including 8 simulations in the previous study. A simulation of NCAR Thermosphere-Ionosphere-Electrodynamics General Circulation Model version 2 (TIE-GCM 2) is also included for comparison with WACCM-X. TEC and foF2 changes from quiet-time background are considered to evaluate the model performance on the storm impacts. For evaluation, we employ 4 skill scores: Correlation coefficient (CC), root-mean square error (RMSE), ratio of the modeled to observed maximum percentage changes (Yield), and timing error(TE). It is found that the models tend to underestimate the storm-time enhancements of foF2 (F2-layer critical frequency) and TEC (Total Electron Content) and to predict foF2 and/or TEC better in the North America but worse in the Southern Hemisphere. The ensemble simulation for TEC is comparable to results from a data assimilation model (Utah State University-Global Assimilation of Ionospheric Measurement (USU-GAIM)) with differences in skill score less than 3% and 6% for CC and RMSE, respectively.

The explosive eruption of the Hunga-Tonga volcano in the southwest Pacific at 0415UT on 15 January 2022 triggered gigantic atmospheric disturbances with surface air pressure wave propagating around the globe in Lamb mode. In space, concentric traveling ionosphere disturbances (CTIDs) are also observed as a manifestation of air pressure acoustic waves in New Zealand ~0500UT and Australia ~0630UT. As soon as the wave reached central Australia ~0800UT, CTIDs appeared simultaneously in the northern hemispheres through magnetic field line conjugate effect, which is much earlier than the arrival of the air pressure wave to Japan after 1100UT. Combining observations over Australia and Japan between 0800-1000UT, CTIDs with characteristics of phase velocities of 320-390 m/s are observed, matching with the dispersion relation of Lamb mode. The arrival of atmospheric Lamb wave to Japan later created in situ CTIDs showing the same Lamb mode characteristics as the earlier arriving CTIDs.

The Tonga volcano eruption of 15 January 2022 unleashed a variety of atmospheric perturbations, coinciding with the recovery phase of a geomagnetic storm. The ensuing thermospheric variations created rare display of extreme poleward-expanding conjugate plasma bubbles seen in the rate of total electron content index (ROTI) over 100-150°E. This is associated with fluctuations in FORMOSAT-7/COSMIC-2 (F7/C2) ion-density measurements and spread-F signatures in ionograms, reaching ~40°N geographic latitude. This was preceded by an unusually strong pre-reversal enhancement (PRE) in the global ionospheric specification (GIS) electron density profiles derived from F7/C2 observations. The GIS further revealed a decrease of equatorial ionization anomaly (EIA) crest density due to the storm impact. A sharp decrease of E-region conductivity by volcano-induced waves, combined with enhanced F-region wind over EIA with less ion-drag apparently intensified the PRE. The strong PRE and seed perturbations from the volcano-induced waves likely further triggered super plasma bubble activity.

We report on a new method to calibrate the FORMOSAT-7/COSMIC-2 (F7/C2) ion velocity meter (IVM) in-situ ion density data using the Tri GNSS (Global Navigation Satellite System) Radio-occultation System (TGRS) differential total electron content (TEC) data. This calibration is made using collocated measurements from the IVM and TGRS instruments on the same satellite. We found that the IVM ion density is about 8-15% lower than the TGRS derived density at the insertion orbit (~ 710 km) and 5% higher at the mission operation orbit (~ 540 km). Using a linear correction specific to orbit altitude and satellite, an adjustment is implemented for the IVM density data. These corrections remove the offsets between the IVM and TGRS density measurements. We believe the corrected densities eliminate any inter-spacecraft discrepancy in the IVM density data, making it suitable for use in multi-satellite scientific investigations of longitudinal and local time variations and space weather operational applications.