Ions in the F region ionosphere at 150-400 km altitude consist mainly of molecular NO+ and O2+, and atomic O+. Incoherent scatter (IS) radars are sensitive to the molecular-to-atomic ion density ratio, but its effect to the observed incoherent scatter spectra is almost identical with that of the ion temperature. It is thus very difficult to fit both the ion temperature and the fraction of O+ ions to the observed spectra. In this paper, we introduce a novel combination of Bayesian filtering, smoothness priors, and chemistry modeling to solve for F1 region O+ ion fraction from EISCAT Svalbard IS radar (75.43° corrected geomagnetic latitude) data during the international polar year (IPY) 2007-2008. We find that the fraction of O+ ions in the F1 region ionosphere is controlled by ion temperature and electron production. The median value of the molecular-to-atomic ion transition altitude during IPY varies from 187 km at 16-17 MLT to 208 km at 04-05 MLT. The ion temperature has maxima at 05-06 MLT and 15-16 MLT, but the transition altitude does not follow the ion temperature, because photoionization lowers the transition altitude. A daytime transition altitude maximum is observed in winter, when lack of photoionization leads to very low daytime electron densities. Both ion temperature and the molecular-to-atomic ion transition altitude correlate with the Polar Cap North geomagnetic index. The annual medians of the fitted transition altitudes are 14-32 km lower than those predicted by the International Reference Ionosphere.

The most detrimental geomagnetically induced currents (GICs) documented to date have all taken place during geomagnetic storms. Yet, the probability of GICs throughout geomagnetic storms driven by different solar wind transients, such as high-speed streams/stream interaction regions (HSS/SIR) or interplanetary coronal mass ejection (ICME) sheaths and magnetic clouds (MC), is poorly understood. We present an algorithm to detect geomagnetic storms and storm phases, resulting in a catalog of 755 geomagnetic storms from January 1996 to June 2023 with the solar wind drivers. Using these storms and the IMAGE magnetometer network, we study the temporal and spatial evolution of spikes in the external dH\textsubscript{ext}/dt greater than 0.5 nT/s during geomagnetic storms driven by HSS/SIR, sheaths and MCs. Spikes occur more often toward the end of the storm main phase for HSS/SIR and MC-driven storms, while sheaths have spikes throughout the entire main phase. During the main phase most spikes occur in the morning sector around 05 magnetic local time (MLT) and the extent in MLT is narrowest for MCs and widest for sheaths. However, spikes in the pre-midnight sector during the main and recovery phases are most prominent for HSS/SIR-driven storms. During the storm sudden commencement (SSC), three MLT hotspots exist, the post-midnight at 04 MLT, pre-noon at 09 MLT and afternoon at 15 MLT. The pre-noon hotspot has the highest probability of spikes and the widest extent in magnetic latitude.

During geomagnetically active times, Joule heating in the Lower Thermosphere - Ionosphere is a significant energy source, greatly affecting density, temperature, composition and circulation. At the same time, Joule heating and the associated Pedersen conductivity are amongst the least known parameters in the upper atmosphere in terms of their quantification and spatial distribution, and their parameterization by geomagnetic parameters shows large discrepancies between estimation methodologies, primarily due to a lack of comprehensive measurements in the region where they maximize. In this work we perform a long-term statistical comparison of Joule heating as calculated by the NCAR Thermosphere - Ionosphere - Electrodynamics General Circulation Model (TIE-GCM) and as obtained through radar measurements by the European Incoherent Scatter Scientific Association (EISCAT). Statistical estimates of Joule heating and Pedersen conductivity are obtained from a simulation run over the 11 year period spanning from 2009 until 2019 and from radar measurements over the same period, during times of radar measurements. The results are statistically compared in different Magnetic Local Time sectors and Kp level ranges in terms of median values and percentiles of altitude profiles. It is found that Joule heating and Pedersen conductivity are higher on average in TIE-GCM than in EISCAT for low Kp and are lower than EISCAT for high Kp. It is also found that neutral winds cannot account for the discrepancies between TIE-GCM and EISCAT. Comparisons point towards the need for a Kp-dependent parameterization of Joule heating in TIE-GCM to account for the contribution of small scale effects.

This study presents an improved method to estimate differential energy flux, auroral power and field-aligned current of electron precipitation from incoherent scatter radar data. The method is based on a newly developed data analysis technique that uses Bayesian filtering to fit altitude profiles of electron density, electron temperature, and ion temperature to observed incoherent scatter spectra with high time and range resolutions. The electron energy spectra are inverted from the electron density profiles. Previous high-time resolution fits have relied on the raw electron density, which is calculated from the backscattered power assuming that the ion and electron temperatures are equal. The improved technique is applied to one auroral event measured by the EISCAT UHF radar and it is demonstrated that the effect of electron heating on electron energy spectra, auroral power and upward field-aligned current can be significant at times. Using the fitted electron densities instead of the raw ones may lead to wider electron energy spectra and auroral power up to 75% larger. The largest differences take place for precipitation that produces enhanced electron heating in the upper E region, and in this study correspond to fluxes of electrons with peak energies from 3 to 5 keV. Finally, the auroral power estimates are verified by comparison to the 427.8 nm auroral emission intensity, which show good correlation. The improved method makes it possible to calculate unbiased estimates of electron energy spectra with high time resolution and thereby to study rapidly varying aurora.

Incoherent scatter (IS) radars are invaluable instruments for ionospheric physics, since they observe altitude profiles of electron density (Ne), electron temperature (Te), ion temperature (Ti) and line-of-sight plasma velocity (Vi) from ground. However, the temperatures can be fitted to the observed IS spectra only when the ion composition is known, and resolutions of the fitted plasma parameters are often insufficient for auroral electron precipitation, which requires high resolutions in both range and time. The problem of unknown ion composition has been addressed by means of the full-profile analysis, which assumes that the plasma parameter profiles are smooth in altitude, or follow some predefined shape. In a similar manner, one could assume smooth time variations, but this option has not been used in IS analysis. We propose a plasma parameter fit technique based on Bayesian filtering, which we have implemented as an additional Bayesian Filtering Module (BAFIM) in the GUISDAP analysis package. BAFIM allows us to control gradients in both time and range directions for each plasma parameter separately. With BAFIM we can fit F1 region ion composition together with Ne, Te, Ti, and Vi, and we have reached 4 s/900 m time/range steps in four-parameter fits of Ne, Te, Ti and Vi in E region observations of auroral electron precipitation.

This study considers 28 geomagnetic storms with Dst $\leq-50$ nT driven by high-speed streams (HSSs) and associated stream interaction regions (SIRs) during 2010-2017. Their impact on ionospheric horizontal and field-aligned currents (FACs) have been investigated using superposed epoch analysis of SuperMAG and AMPERE data, respectively. The zero epoch ($t_0$) was set to the onset of the storm main phase. Storms begin in the SIR with enhanced solar wind density and compressed southward oriented magnetic field. The integrated FAC and equivalent currents maximise 40 and 58 min after $t_0$, respectively, followed by a small peak in the middle of the main phase ($t_0$+4h), and a slightly larger peak just before the Dst minimum ($t_0$+5.3h). The currents are strongly driven by the solar wind, and the correlation between the Akasofu $\varepsilon$ and integrated FAC is $0.90$. The number of substorm onsets maximises near $t_0$. The storms were also separated into two groups based on the solar wind dynamic pressure p_dyn in the vicinity of the SIR. High p_dyn storms reach solar wind velocity maxima earlier and have shorter lead times from the HSS arrival to storm onset compared with low p_dyn events. The high p_dyn events also have sudden storm commencements, stronger solar wind driving and ionospheric response at $t_0$, and are primarily responsible for the first peak in the currents after $t_0$. After $t_0+2$ days, the currents and number of substorm onsets become higher for low compared with high p_dyn events, which may be related to higher solar wind speed.

We present a statistical investigation of the effects of interplanetary magnetic field (IMF) on hemispheric asymmetry in auroral currents. Nearly six years of magnetic field measurements from Swarm A and C satellites are analyzed. Bootstrap resampling is used to remove the difference in the number of samples and IMF conditions between the local seasons and the hemispheres. Currents are stronger in Northern Hemisphere (NH) than Southern Hemisphere (SH) for IMF B$y^+$ in NH (B$y^-$ in SH) in most local seasons under both signs of IMF B$z$. For B$y^-$ in NH (B$y^+$ in SH), the hemispheric difference in currents is small except in local winter when currents in NH are stronger than in SH. During B$y^+$ and B$z^+$ in NH (B$y^-$ and B$z^+$ in SH), the largest hemispheric asymmetry occurs in local winter and autumn when the NH/SH ratio of field-aligned current (FAC) is 1.18$\pm$0.09 in winter and 1.17$\pm$0.09 in autumn. During B$y^+$ and B$z^-$ in NH (B$y^-$ and B$z^-$ in SH), the largest asymmetry is observed in local autumn with NH/SH ratio of 1.16$\pm$0.07 for FAC. We also find an explicit B$y$ effect on auroral currents in a given hemisphere: on average B$y^+$ in NH and B$y^-$ in SH causes larger currents than vice versa. The explicit B$y$ effect on divergence-free (DF) current during IMF B$z^+$ is in very good agreement with the B$y$ effect on the cross polar cap potential (CPCP) from the Super Dual Auroral Radar Network (SuperDARN) dynamic model except at SH equinox and NH summer.

We present a new analysis technique for estimating 2D neutral wind pattern using data from a single Scanning Doppler Imager (SDI) or a combination of SDIs, incoherent scatter radars (ISR) and Fabry-Perot interferometers (FPI) within overlapping field-of-views. Neutral wind plays an important role in ionospheric electrodynamics and Ionosphere-Thermosphere coupling, by for example affecting the Joule heating rates and plasma transport. However, reliable and extensive measurements of the neutral wind are rather difficult to obtain. Pointwise measurements can be obtained with ISRs or FPIs, but these measurements can not provide 2D latitude-longitude maps of the neutral wind pattern needed in mesospheric studies. A Scanning Doppler Imager can measure the line-of-sight (LOS) component of the neutral wind in dozens of directions simultaneously. However, further modeling is needed to convert the LOS velocities into 2D velocity maps. Unfortunately these maps are far from unique, as perpendicular velocities (e.g. rotation around the measurement site) are not visible in the LOS data. This can be mitigated by combining data from several nearby SDIs, or a combination of SDIs, FPIs and ISRs. Our analysis technique is based on fitting the LOS data with special vector basis functions called Spherical Elementary Current Systems (SECS). In this approach the wind is naturally divided into curl-free and divergence-free components, and there is no need to provide any explicit boundary conditions on the wind pattern. We present several synthetic test scenarios as well as first results using data from SDIs located in Alaska. Using the synthetic test scenarios we further estimate optimal locations for 2 or 3 SDIs that could be located around the future EISCAT_3D radar system in northern Scandinavia.

We present a statistical investigation of the seasonal effect on hemispheric asymmetry in the auroral currents during low (Kp $<$ 2) and high (Kp $\geq$ 2) geomagnetic activity. Five years of magnetic data from the Swarm satellites has been analysed by applying the spherical elementary current system (SECS) method. Bootstrap resampling has been used to remove the seasonal differences between the hemispheres in the dataset. In general, the currents are larger in the Northern Hemisphere (NH) than in the Southern Hemisphere (SH). Asymmetry is larger during low than high Kp, and during winter and autumn than summer and spring. The NH/SH ratio for FACs in winter, autumn, spring and summer are 1.17 $\pm$ 0.05, 1.14 $\pm$ 0.05, 1.07 $\pm$ 0.04 and 1.02 $\pm$ 0.04, respectively. The largest asymmetry is observed during low Kp winter, when the excess in the NH currents is 21$\pm$5\% in FAC, 14 $\pm$ 3\% in curl-free (CF), and 10$\pm$3\% in divergence-free (DF) current. We also find that evening sector (13-24 MLT) contributes more to the high NH/SH ratio than the morning (01-12 MLT) sector. The physical mechanisms producing the hemispheric asymmetry are not presently understood. We calculated the background ionospheric conductances during low Kp conditions from the IRI, NRLMSISE and CHAOS models. The results indicate that only a small part of the hemispheric asymmetry can be explained by variations in the solar induced conductances.