A novel Python module called PyIRTAM was developed in addition to the previously released PyIRI model. The PyIRI computational approach was extended to work with Real-Time Assimilative Model (IRTAM) coefficients. This enables global approach to the density specification, utilizing computation of the global and diurnal functions and their multiplication with coefficients in the matrix form. PyIRTAM tool is made publicly available through PyPI and GitHub.

A novel model called PyIRI was recently developed. It offers a fully Python alternative to the widely used FORTRAN International Reference Ionosphere (IRI) model to construct the ionospheric electron density for the entire day and on the entire global grid in one computation, which has a significantly lower computational overhead. PyIRI introduced a novel approach to the computation of the global and diurnal functions and their matrix multiplication with Consultative Committee on International Radio (CCIR) coefficients, that enabled this global approach for the density specification. Since the IRI-based Real-Time Assimilative Model (IRTAM) produces coefficients in a similar format as CCIR coefficients, the PyIRI software was extended to be able to work with IRTAM coefficients using the same computationally efficient approach. This technical note describes the PyIRTAM software and shows examples of how to use it. PyIRTAM is 21 times more efficient than the classical IRTAM. PyIRTAM tool is made publicly available.

ANCHOR is a novel assimilative model developed at the U.S. Naval Research Laboratory. It extracts ionospheric parameters from RO and ionosonde data and assimilates them as point measurements into the maps of the background parameters using Kalman Filter approach. This paper introduces the ANCHOR algorithm, discusses its coordinate system and background, explains the background covariance formation, discusses the extraction of the ionospheric parameters from the data and the assimilation process, and, finally, shows the results of the observing system simulation experiment.

The EXospheric TEMeratures on a PoLyhedrAl gRid (EXTEMPLAR) method predicts the neutral densities in the thermosphere. The performance of this model has been evaluated through a comparison with the Air Force High Accuracy Satellite Drag Model (HASDM). The Space Environment Technologies (SET) HASDM database that was used for this test spans the 20 years 2000 through 2019, containing densities at 3 hour time intervals at 25 km altitude steps, and a spatial resolution of 10 degrees latitude by 15 degrees longitude. The upgraded EXTEMPLAR that was tested uses the newer Naval Research Laboratory MSIS 2.0 model to convert global exospheric temperature values to neutral density as a function of altitude. The revision also incorporated time delays that varied as a function of location, between the total Poynting flux in the polar regions and the exospheric temperature response. The density values from both models were integrated on spherical shells at altitudes ranging from 200 to 800 km. These sums were compared as a function of time. The results show an excellent agreement at temporal scales ranging from hours to years. The EXTEMPLAR model performs best at altitudes of 400 km and above, where geomagnetic storms produce the largest relative changes in neutral density. In addition to providing an effective method to compare models that have very different spatial resolutions, the use of density totals at various altitudes presents a useful illustration of how the thermosphere behaves at different altitudes, on time scales ranging from hours to complete solar cycles.

The formation of layers at mid-latitudes has been related to neutral winds activity at altitudes below 130km in the Mesosphere and Lower Thermosphere (MLT). Recent SAMI3 simulations by Krall et al. (2020) of ionospheric metallic layers at Arecibo suggest that forces induced by the meridional winds cause low altitude layers near 100 km. However, the classic mechanism, originally proposed by Whitehead (1961), correctly states that zonal wind shear has a bigger effect than meridional wind shear in the lower E region. Haldoupis and Shalimov (2021), referring to observations of ionosonde-based sporadic E statistics and radio occultation sporadic E measurements using low Earth orbiting satellites, support the idea that zonal winds dominate layer formation at these altitudes, apparently disputing the findings of Krall et al. (2020). Perhaps the latest technique to continuously measure mid-latitude MLT daytime neutral winds was developed by Hysell et al. (2014). That technique used a unique configuration of the Arecibo radar dual-beam. Unfortunately, since Arecibo lost the capability of the dual-beam in 2017 (when one antenna was destroyed by Hurricane Maria), there are only few valuable data sets that can help elucidate the origin of the lower altitude layers at Arecibo. We present Arecibo neutral wind data correlated with lower altitude layers. While not disputing current theory, we find that, near 100 km, meridional neutral wind shear can be much stronger than zonal wind shear when a layer is present, with the meridional shear correctly positioned to support the layer. We also present a complete analysis of the vertical ion drift, including declination, where the meridional winds become more important and with a reversed mechanism for altitudes below 115km for Arecibo conditions. References: Haldoupis et al. (2021), https://doi.org/10.1016/j.jastp.2021.105537 Hysell et al. (2014), http://doi.org/10.1002/2013JA019621 Krall et al. (2020), https://doi.org/10.1029/2019JA027297

The Naval Research Laboratory (NRL) Sami3 is Also a Model of the Ionosphere (SAMI3) ionosphere/plasmasphere code is used to examine the physics of metallic layers at altitudes from 80 to 160 km. Results are presented near the simulated location of the Arecibo observatory (18N, 66W). We find that simulations, using winds from the empirical horizontal wind model (HWM14), produce layers consistent with those observed at Arecibo. Specifically, we find upper semidiurnal and lower diurnal traces similar to those identified in previous observational surveys. While metallic layers are shaped by meridional winds, zonal winds, and electric fields, much of the observed structure is found if only meridional wind forces are included in the model. Stratification below 110 km, where the ions are very weakly magnetized, is supported mainly by meridional wind shear.