In this study, we use measurements from over 4,735 globally distributed Global Navigation Satellite System receivers to track the progression of traveling ionospheric disturbances (TIDs) associated with the 15 January 2022 Hunga Tonga‐Hunga Ha'apai submarine volcanic eruption. We identify two distinct Large Scale traveling ionospheric disturbances (LSTIDs) and several subsequent Medium Scale traveling ionospheric disturbances (MSTIDs) that propagate radially outward from the eruption site. Within 3,000 km of epicenter, LSTIDs of >1,600 km wavelengths are initially observed propagating at speeds of ∼950 and ∼555 ms−1, before substantial slowing to ∼600 and ∼390 ms−1, respectively. MSTIDs with speeds of 200–400 ms−1 are observed for 6 hrs following eruption, the first of which comprises the dominant global ionospheric response and coincides with the atmospheric surface pressure disturbance associated with the eruption. These are the first results demonstrating the global impact of the Tonga eruption on the ionospheric state.

In this study, we present a topside model representation to be used by the Empirical Canadian High Arctic Ionospheric Model (E‐CHAIM). In the process of this, we also present a comprehensive evaluation of the NeQuick's, and by extension the International Reference Ionosphere's, topside electron density model for middle and high latitudes in the Northern Hemisphere. Using data gathered from all available incoherent scatter radars, topside sounders, and Global Navigation Satellite System Radio Occultation satellites, we show that the current NeQuick parameterization suboptimally represents the shape of the topside electron density profile at these latitudes and performs poorly in the representation of seasonal and solar cycle variations of the topside scale thickness. Despite this, the simple, one variable, NeQuick model is a powerful tool for modeling the topside ionosphere. By refitting the parameters that define the maximum topside scale thickness and the rate of increase of the scale height within the NeQuick topside model function, r and g, respectively, and refitting the model's parameterization of the scale height at the F region peak, H0, we find considerable improvement in the NeQuick's ability to represent the topside shape and behavior. Building on these results, we present a new topside model extension of the E‐CHAIM based on the revised NeQuick function. Overall, root‐mean‐square errors in topside electron density are improved over the traditional International Reference Ionosphere/NeQuick topside by 31% for a new NeQuick parameterization and by 36% for a newly proposed topside for E‐CHAIM.

In this study, we present a bottomside model representation to be used by the Empirical Canadian High Arctic Ionospheric Model (E-CHAIM). This model features a new approach to modeling the bottomside electron density; namely, instead of modelling electron density directly, E-CHAIM models the altitude profile of the scale thickness of a single bottomside layer. In this approach, the curvature in the bottomside associated with the E region and F1 layer is represented in the scale thickness domain as a peak function centered at the layer peak altitude. The use of this approach ensures the production of explicitly doubly differentiable bottomside electron density profiles and directly avoids issues known to exist within current standards, such as the International Reference Ionosphere (IRI), which has discontinuities in space, time, and in the vertical electron density gradient. In terms of performance, after removing the impacts of hmF2 and NmF2, the new E-CHAIM profile function generally performs comparably to the IRI, with bottomside TEC from both models within 2.0 TECU (1 TECU = 10 16 e/m 3 ) of observations. More specifically, the E-CHAIM bottomside is demonstrated to outperform the IRI bottomside function in the F region during low solar activity periods with respect to incoherent scatter radar observations. At high latitudes, E-CHAIM tends to outperform the IRI during winter months by between 10% and 40% of NmF2 while being outperformed by the IRI by between 10% and 25% of NmF2 during summer periods, mainly during the daytime at high solar activity.

The high latitude ionosphere has historically been a challenging system to model (Buchert, 2020;Lockwood et al., 1990;Rasmussen et al., 1986). A rich collection of external drivers and interactions drive ionospheric behavior, including strong electric fields, magnetospheric coupling via particle precipitation and current systems, and rapid changes in the thermospheric state. These dynamic conditions, paired with a lack of high latitude observations when compared to mid and low latitudes, present a substantial problem for operational ionospheric modeling. With increased interest in polar ionospheric monitoring (Thayaparan et al., 2018) and High Frequency (HF) communications, it is now imperative that a near-real-time, operational model of high latitude electron density be developed and deployed for use in this region.

String cavitation has been studied in an optical automotive size fuel injector with true-scale flow geometry at injection pressures of up to 2050 bar. The multihole nozzle geometry studied allowed observation of the hole-to-hole vortex interaction and, in particular, that of a bridging vortex in the sac region between the holes. A dependency on Reynolds number was observed in the formation of the visible, vapor filled vortex cores. Above a threshold Reynolds number, their formation and appearance during a 2 ms injection event was repeatable and independent of upstream pressure and cavitation number.