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 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 four 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-1000 km horizontal wavelength for front shocks, going around the globe for three times, passing six times over the continental US in 100 hours since the eruption. TIDs following the shock fronts developed for ~8 hours 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 of 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.

This paper investigates the local and global ionospheric responses to the 2022 Tonga volcano eruption, using ground-based Global Navigation Satellite System (GNSS) total electron content (TEC), Swarm in-situ plasma density measurements, the Ionospheric Connection Explorer (ICON) Ion Velocity Meter (IVM) data, and ionosonde measurements. The main results are as follows: (1) A significant local ionospheric hole of more than 10 TECU depletion was observed near the epicenter ~45~min after the eruption, comprising of several cascading TEC decreases and quasi-periodic oscillations. Such a deep local plasma hole was also observed by space-borne in-situ measurements, with an estimated horizontal radius of 10-15 deg and persisted for more than 10 hours in ICON-IVM ion density profiles until local sunrise. (2) Pronounced post-volcanic evening equatorial plasma bubbles (EPBs) were continuously observed across the wide Asia-Oceania area after the arrival of volcano-induced waves; these caused a Ne decrease of 2-3 orders of magnitude at Swarm/ICON altitude between 450-575~km, covered wide longitudinal ranges of more than 140 deg and lasted around 12 hours. (3) Various acoustic-gravity wave modes due to volcano eruption were observed by accurate Beidou geostationary orbit (GEO) TEC, and the huge ionospheric hole was mainly caused by intense shock-acoustic impulses. TEC rate of change index revealed globally propagating ionospheric disturbances at a prevailing Lamb-wave mode of ~315 m/s; the large-scale EPBs could be seeded by acoustic-gravity resonance and coupling to less-damped Lamb waves, under a favorable condition of volcano-induced enhancement of dusktime plasma upward ExB drift and postsunset rise of the equatorial ionospheric F-layer.

We demonstrate a novel method for observing Large Scale Traveling Ionospheric Disturbances (LSTIDs) using high frequency (HF) amateur radio reporting networks, including the Reverse Beacon Network (RBN), Weak Signal Propagation Reporter Network (WSPRNet), and PSKReporter. LSTIDs are quasi-periodic variations in ionospheric densities with horizontal wavelengths > 1000 km and periods between 30 to 180 min. On 3 Nov 2017, LSTID signatures were observed simultaneously over the continental United States in amateur radio, SuperDARN HF radar, and GNSS Total Electron Content with a period of ~2.5 hr, propagation azimuth of ~163°, horizontal wavelength of ~1680 km, and phase speed of ~1200 km/hr. SuperMAG SME index enhancements and Poker Flat Incoherent Scatter Radar measurements suggest the LSTIDs were driven by auroral electrojet intensifications and Joule heating. This novel measurement technique has applications in future scientific studies and for assessing the impact of LSTIDs on HF communications.

We report on an extreme ionospheric plasma density enhancement and Global Positioning System (GPS) scintillation at dawn, observed within the expanding equatorial ionization anomaly (EIA). The total electron content (TEC) in central America reached 50~TECu at sunrise, the value almost twice as high as the normal afternoon peak. The enhanced EIA expanded poleward and westward from just below 20$^\circ$ magnetic latitude (MLAT) to beyond 30$^\circ$ MLAT at sunrise. The chief ramification of the enhanced EIA was strong GPS scintillation which was observed poleward of 30$^\circ$ northern MLAT and lasted until 8:00 local time. In total, the scintillation lasted for $\sim$5~hours at latitudes north of 20$^\circ$MLAT in central America.

Traveling Ionospheric Disturbances (TIDs) are propagating variations in ionospheric electron densities that affect radio communications and can help with understanding energy transport throughout the coupled magnetosphere-ionosphere-neutral atmosphere system. Large scale TIDs (LSTIDs) have periods T ≈30-180 min, horizontal phase velocities vH≈‍100-‍250 m/s, and horizontal wavelengths H>1000 km and are believed to be generated either by geomagnetic activity or lower atmospheric sources. TIDs create concavities in the ionospheric electron density profile that move horizontally with the TID and cause skip-distance focusing effects for high frequency (HF, 3-30 MHz) radio signals propagating through the ionosphere. The signature of this phenomena is manifest as quasi-periodic variations in contact ranges in HF amateur radio communication reports recorded by automated monitoring systems such as the Weak Signal Propagation Reporting Network (WSPRNet) and the Reverse Beacon Network (RBN). In this study, members of the Ham Radio Science Citizen Investigation (HamSCI) present a climatology of LSTID activity using RBN and WSPRNet observations on the 1.8, 3.5, 7, 14, 21, and 28 MHz amateur radio bands from 2017. Results will be organized as a function observation frequency, longitudinal sector (North America and Europe), season, and geomagnetic activity level. Connections to geospace are explored via SYM-H and Auroral Electrojet indexes, while neutral atmospheric sources are explored using NASA’s Modern-Era Retrospective Analysis for Research and Applications Version 2 (MERRA-2).

This work conducts a statistical study of the subauroral polarization stream (SAPS) feature in the North American sector using Millstone Hill incoherent scatter radar measurements from 1979 to 2019, which provides a comprehensive SAPS climatology using a significantly larger database of radar observations than was used in seminal earlier works. Key features of SAPS and associated Ne/Ti/Te are investigated using a superposed epoch analysis method. The characteristics of these parameters are investigated with respect to magnetic local time, season, geomagnetic activity, solar activity, and interplanetary magnetic field orientation, respectively. The main results are as follows: (1) Conditions for SAPS are more favorable for dusk than near midnight, for winter compared to summer, for active geomagnetic periods compared to quiet time, for solar minimum compared to solar maximum, and for IMF conditions with negative By and negative Bz. (2) SAPS is usually associated with a midlatitude trough of 15–20\% depletion in the background density. The SAPS-related trough is more pronounced in the postmidnight sector and near the equinoxes. (3) Subauroral ion and electron temperatures exhibit a 3–8\% (50–120 K) enhancement in SAPS regions, which tend to have higher percentage enhancement during geomagnetically active periods and at midnight. Ion temperature enhancements are more favored during low solar activity periods, while the electron temperature enhancement remains almost constant as a function of the solar cycle. (4) The electron thermal content, Te \times Ne, in the SAPS associated region is strongly dependent on 1/Ne, with Te exhibiting a negative correlation with respect to $Ne$.

We report new findings of total electron content (TEC) perturbations in the southern hemisphere at conjugate locations to the northern eclipse on 21 August 2017. We identified a persistent conjugate TEC depletion by 10-15\% during eclipse time, elongating along magnetic latitudes with $\sim$5$^\circ$ latitudinal width, moving equatorward, and becoming most pronounced at lower magnetic latitudes ($<$20$^\circ$S) when the ionosphere in the northern low latitudes were masked. This depletion was coincident with a weakening of the southern crest of the equatorial ionization anomaly (EIA), while the northern EIA crest stayed almost undisturbed or was slightly enhanced. We suggest these conjugate perturbations were associated with dramatic eclipse initiated plasma pressure reductions in the flux tubes, with a large portion of shorter tubes located at low latitudes underneath the Moon’s shadow. These small L-shell tubes transversed the F region ionosphere at low and equatorial latitudes. The plasma pressure gradient was markedly skewed northward in the flux tubes at low and equatorial latitudes, as was the neutral pressure. These effects caused a general northward motion tendency for plasma within the flux tubes, and inhibited normal southward diffusion of equatorial fountain plasma into the southern EIA region.

A newly-released, novel ionospheric dataset of global gridded vertical total electron content (VTEC) is introduced in this paper. This VTEC dataset, provided by Los Alamos National Laboratory (LANL), is derived from very-high frequency (VHF; defined as 30-300 MHz) broadband radio-frequency (RF) measurements of lightning made by U.S. Department of Defense sensing systems on board Global Positioning System (GPS) satellites. This paper presents the new dataset (LANL VTEC), discusses the errors inherent in VHF TEC estimation due to ionospheric dispersion, and compares the LANL VTEC to two community standard VTEC gridded products: Jet Propulsion Laboratory’s Global Ionospheric Model (JPL GIM) and the CEDAR community’s Open Madrigal VTEC gridded measurements of L-band GNSS (global navigation satellite systems) TEC. We find that the LANL VTEC data has an offset of 3 TECU from CEDAR Madrigal GNSS VTEC, and a full-width-half-maximum (FWHM) of 6 TECU. In comparison, the offset between LANL VTEC and the JPL GIM model is -3 TECU, but with a FWHM of 5 TECU. We also compare to Jason-3 VTEC measurements over the ocean, finding an offset of less than 0.5 TECU and a FWHM of < 5 TECU. Because this technique uses a completely different methodology to determine TEC, the sources of errors are distinct from the typical ground-based GNSS L-band (GHz) TEC measurements. Also, because it is derived from RF lightning signals, this dataset provides measurements in regions that are not well covered by ground-based GPS measurements, such as over oceans and over central Africa.

The unambiguous 2-dimensional (2D) maps of OI 135.6 nm radiance retrieved from the Global Observation of Limb and Disk (GOLD) after sunset are compared with the total electron content (TEC) maps measured by GPS receivers in the America sector. The OI 135.6 nm radiance observed by GOLD is an indicator of the peak electron density of the ionosphere, while the TEC depends on the total electron density in the column. The comparisons show that both of them are able to capture the large structures in the equatorial ionization anomaly (EIA) well, and sometimes they both also capture bubbles. Both show that the ionosphere after sunset is quite dynamic and has strong day-to-day variability. A statistical study has also been carried out to check the occurrence rate of bubbles and the apparent EIA structure between Oct 17, 2018 to May 31, 2019. GOLD is able to image the areas where it is difficult to situate GPS receivers such as the ocean, while TEC data covers the full-diurnal cycle. In all, the GOLD and TEC have valuable synergy to allow us to gain a better understanding of the equatorial ionosphere.