Ionospheric Current System Research Articles

Two distinct current systems in the ionosphere of Mars

When the solar wind interacts with the ionosphere of an unmagnetized planet, it induces currents that form an induced magnetosphere. These currents and their associated magnetic fields play a pivotal role in controlling the movement of charged particles, which is essential for understanding the escape of planetary ions. Unlike the well-documented magnetospheric current systems, the ionospheric current systems driven by solar wind and atmospheric neutral winds have not been quantitatively observed, which constrains the quantification of energy transfer from stars to these planets. Here, utilizing eight years of data from the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, we investigate the global distribution of ionospheric currents on Mars. We identify two distinct current systems in the ionosphere: one aligns with the solar wind electric field but exhibits hemispheric asymmetry perpendicular to the solar wind electric field direction; the other corresponds to the flow pattern of annually averaged neutral winds. We propose that these two current systems are driven by the solar wind and atmospheric neutral winds, respectively. Our findings reveal that Martian ionospheric dynamics are influenced by the neutral winds from below and the solar wind from above, highlighting the complex and intriguing nature of current systems on unmagnetized planets.

The Ionospheric Leg of the Substorm Current Wedge: Combining Iridium and Ground Magnetometers

AbstractUtilizing magnetic field measurements made by the Iridium satellites and by ground magnetometers in North America we calculate the full ionospheric current system and investigate the substorm current wedge. The current estimates are independent of ionospheric conductance, and are based on estimates of the divergence‐free (DF) ionospheric current from ground magnetometers and curl‐free (CF) ionospheric currents from Iridium. The DF and CF currents are represented using spherical elementary current systems (SECS), derived using a new inversion scheme that ensures the current systems' spatial scales are consistent. We present 18 substorm events and find a typical substorm current wedge (SCW) in 12 events. Our investigation of these substorms shows that during substorm expansion, equivalent field‐aligned currents (EFACs) derived with ground magnetometers are a poor proxy of the actual FAC. We also find that the intensification of the westward electrojet can occur without an intensification of the FACs. We present theoretical investigations that show that the observed deviation between FACs estimated with satellite measurements and ground‐based EFACs are consistent with the presence of a strong local enhancement of the ionospheric conductance, similar to the substorm bulge. Such enhancements of the auroral conductance can also change the ionospheric closure of pre‐existing FACs such that the ground magnetic field, and in particular the westward electrojet, changes significantly. These results demonstrate that attributing intensification of the westward electrojet to SCW current closure can yield false understanding of the ionospheric and magnetospheric state.

Response of total electron content to the October 25, 2022 partial solar eclipse from high to low latitudes in the Euro-Asian region

Solar eclipses (SEs) pertain to high-energy sources. They are capable of causing significant disturbances in all geoshells and geophysical fields. Despite the almost century-long history of studying the influence of eclipses on geoshells, the task of observing and analyzing disturbances in them remains relevant. This is explained by the fact that the response of the media significantly depends on the state of atmospheric and space weather, the location and time of observation, the SE magnitude, etc. The aim of this paper is to present the results of a statistical analysis of the features of the total electron content (TEC) variations during the October 25, 2022 SE from high to low latitudes in the Euro-Asian region. Data from Global Navigation Satellite System were used for the analysis. The error in estimating TEC in the vertical column did not exceed 0.1 TECU. For the first time, the response of the ionosphere to an SE has been studied on a global scale (from the Reykjavik station to the Novosibirsk station) using methods of statistical analysis, including the time moments after sunset. The decrease in TEC reached 1.0–19.0 TECU, and the relative decrease was approximately –0.32. Changes in all TEC disturbance parameters depended in a complex manner on the relative obscuration area of the solar disk. This can be explained by non-monotonic changes in electron density during the daytime, by the intensification of the ionosphere–plasmasphere coupling, the disturbance of the ionospheric current system, vertical plasma drift in crossed electric and magnetic fields, as well as the generation of wave disturbances during the SE period. The delay time of the maximum ionospheric response (maximum |δV| value) with respect to the maximum magnitude of the eclipse was 15–25 min. The duration of the ionospheric response to the SE slightly exceeded the duration of the perturbing source influence.

Tidal Composition Analysis of Global Sq Current System

AbstractIn this study, we examine the spatiotemporal variability of the global equivalent ionospheric current system estimated from geomagnetic Solar‐quiet (Sq) variations at 124 mid‐ to low‐latitude ground‐based magnetometer stations during 1–31 May 2020. In this period, geomagnetic activity was particularly low, that is, the hourly geomagnetic activity index Hp60 did not exceed 4o. The spherical harmonic analysis is performed on the Sq variations to estimate the equivalent current function at each universal time hour. Hourly maps of the global Sq current system are used to evaluate, for the first time, the tidal composition of mid‐latitude Sq currents and its temporal evolution. Although tides are known to be an important driver of Sq currents, the tidal composition analysis was difficult in previous studies that focused on Sq currents at particular longitude sectors. Our results show that the migrating tidal components (DW1, SW2, and TW3) are predominant in both hemispheres. This is as expected from the strong day‐night contrast in ionospheric conductivities. The day‐to‐day variations of the migrating tidal components are, however, not strongly correlated between the two hemispheres, suggesting that these variations arise not only from the global effect of solar radiation on ionospheric conductivities but also from the local effect of neutral winds. Variations associated with non‐migrating tides are also found. Especially, eastward‐propagating diurnal and semidiurnal tides with zonal wavenumber 1 (DE1 and SE1) are the largest non‐migrating components. Their production mechanisms remain to be understood.

Investigation of the influence of high-latitude ionospheric current systems on the results of impedance measurement in the Arctic on the example of the Kola Peninsula

We investigated the influence of high-latitude ionospheric current systems on the results of impedance measurements on the Kola Peninsula. Model calculations and joint analysis of satellite data and ground magnetometer data have shown the possibility of distortions in the interpretation of geophysical data due to the spatial vicinity of current systems in cases where rocks of the Earth’s crust have high electrical resistance.

Investigating the Effect of High-Latitude Ionospheric Current Systems on Results from Measuring Impedance in the Arctic, Using the Example of the Kola Peninsula

Investigating the Effect of High-Latitude Ionospheric Current Systems on Results from Measuring Impedance in the Arctic, Using the Example of the Kola Peninsula

Response of Global Ionospheric Currents to Solar Flares with Extreme Ultraviolet Late Phases

It is known that solar flares can affect the current system of the middle- and low-latitude ionosphere. Most earlier studies have focused on such effects during their impulsive phases. Recent studies have reported flares with a significant extreme ultraviolet (EUV) late phase, the effects of which on ionospheric currents have not yet been investigated. Here, we examine the solar quiet (Sq) currents and equatorial electrojets during two X-class flares with EUV late phases using data from more than 200 ground magnetometers. Our results indicate that the ionospheric currents could be significantly enhanced during the impulsive phase, while the effects of the EUV late phase may increase the global ionospheric currents, but are often weak and thus could be obscured by a change in solar wind conditions. In the X1.8 flare event on 2012 October 23, besides the solar flare effects, the currents were modulated by solar wind pressure. In the X1.3 flare event on 2014 April 25, the solar wind pressure was weak and stable, and the Sq currents were enhanced compared to nonflare conditions. We also found that even weak changes in the solar wind dynamic pressure, with magnitudes as low as ∼2 nPa, which are often ignored, may have an appreciable impact on the global ionospheric current system.

Machine Learning‐Based Emulator for the Physics‐Based Simulation of Auroral Current System

AbstractUsing a machine learning technique called echo state network (ESN), we have developed an emulator to model the physics‐based global magnetohydrodynamic simulation results of REPPU (REProduce Plasma Universe) code. The inputs are the solar wind time series with date and time, and the outputs are the time series of the ionospheric auroral current system in the form of two‐dimensional (2D) patterns of field‐aligned current, potential, and conductivity. We mediated a principal component analysis for a dimensionality reduction of the 2D map time series. In this study, we report the latest upgraded Surrogate Model for REPPU Auroral Ionosphere version 2 (SMRAI2) with significantly improved resolutions in time and space (5 min in time, ∼1° in latitude, and 4.5° in longitude), where the dipole tilt angle is also newly added as one of the input parameters to reproduce the seasonal dependence. The fundamental dependencies of the steady‐state potential and field‐aligned current patterns on the interplanetary magnetic field directions are consistent with those obtained from empirical models. Further, we show that the ESN‐based emulator can output the AE index so that we can evaluate the performance of the dynamically changing results, comparing with the observed AE index. Since the ESN‐based emulator runs a million times faster than the REPPU simulation, it is promising that we can utilize the emulator for the real‐time space weather forecast of the auroral current system as well as to obtain large‐number ensembles to achieve future data assimilation‐based forecast.

Determination of the Effective Collision Frequency of Electrons in the E and D Regions of the High-Latitude Ionosphere from Analysis of Radio Occultation Measurements

Collisions between electrons and neutral molecules are of special interest for the physics of the Earth’s ionosphere, in particular, for determining the ionospheric conductivity and current systems in the lower ionosphere of the planet, as well as elucidating the role they play in attenuating radio waves propagating inside the D and E regions of the ionosphere. The effective collision frequency of electrons can be estimated from laboratory studies of electron mobility in atmospheric gases in combination with rocket measurements of temperature and particle density in the Earth’s upper atmosphere, or it can be determined independently from analysis of radio occultation data. We have developed a method for reconstructing the vertical profiles of the absorption coefficient of decimeter (wavelength ~19 cm) radio waves by solving the inverse problem of signal absorption in the D and E regions of the Earth’s ionosphere. Based on the analysis of radio occultation data from the FORMOSAT-3/COSMIC satellites, the altitude profiles of the absorption coefficient of decimeter (DM) radio waves in the planet’s ionosphere during the geomagnetic storm on June 22–23, 2015, were determined. It is known that the absorption coefficient at a given fixed frequency is directly proportional to both the electron density and the collision frequency of electrons with ions and neutrals. Using the data on the vertical profiles of the absorption coefficient of DM radio waves and the electron density reconstructed from the analysis of FORMOSAT-3/COSMIC radio occultations, we estimated the effective collision frequency of electrons in the D and E regions of the Earth’s high-latitude ionosphere. The practical significance of studying the frequency of electron collisions and the effects of radio wave absorption in the D and E regions of the planet’s ionosphere is associated with maintaining the uninterrupted operation of space radio communication and navigation systems.

Quantifying the ability of magnetohydrodynamic models to reproduce observed Birkeland current and auroral electrojet magnitudes

Although global magnetohydrodynamic (MHD) models have increased in sophistication and are now at the forefront of modeling Space Weather, there is still no clear understanding of how well these models replicate the observed ionospheric current systems. Without a full understanding and treatment of the ionospheric current systems, global models will have significant shortcomings that will limit their use. In this study we focus on reproducing observed seasonal interhemispheric asymmetry in ionospheric currents using the Space Weather Modeling Framework (SWMF). We find that SWMF does reproduce the linear relationship between the electrojets and the FACs, despite the underestimation of the currents’ magnitudes. Quantitatively, we find that at best SWMF is only capturing approximately 60% of the observed current. We also investigate how varying F10.7 effects the ionospheric potential and currents during the summer and winter. We find that simulations ran with higher F10.7 result in lower ionospheric potentials. Additionally, we find that the models do not always replicate the expected behavior of the currents with varying F10.7. This work points to a needed improvement in ionospheric conductance models.

Physical reasons for regular discrepancies in values of the PCN and PCS indices characterizing magnetic activity in the northern and southern polar caps

PC index characterizes magnetic activity generated in polar caps by R1 FAC system related to the solar wind electric field EKL. The PC index is calculated separately by data on magnetic disturbance values δF at stations Thule (PCN) and Vostok (PCS) with use of the “unified PC derivation method” taking into account statistically justified relationships between values of EKL and δF at these stations. Nevertheless, discrepancy between PCN and PCS indices is regularly observed: the value of difference ΔPC = PCN-PCS demonstrates the obvious seasonal variation and strongly increases in periods of high magnetic activity, the ΔPC values being the highest (>1 mV/m) in epochs of solar activity maximum and being the lowest (<0.5 mV/m) in the quiet periods. Analysis of discrepancies ΔPC, carried out by data on the 1-min definitive PCN and PCS indices for 22 year (1997–2022), showed that difference between the daily PCN and PCS values is caused by influence of the azimuthal BY IMF component, which generates the specific BY FAC system with reverse polarity of field-aligned currents in the northern and southern hemispheres. Effect of these currents, which is opposite in northern and southern polar caps, critically distorts the regular ionospheric current systems related to EKL field influence, resulting in discrepancy of corresponding PCN and PCS indices.

Modeling the Climatology of Low‐ and Mid‐Latitude F‐Region Ionospheric Currents Using the Swarm Constellation

AbstractWe present a new empirical model of quiet‐time F‐region ionospheric currents and associated magnetic fields. This model is designed to accurately represent these currents and fields at low and mid latitudes. For each individual Swarm satellite, the preprocessed data is represented as a non‐potential toroidal magnetic field using the Mie representation in a thin‐shell and spherical harmonic expansions. This approach allows to fully separate spatial and climatological variations as well as to assess the robustness of the model with respect to both measurement errors and data sampling with local time. The obtained model describes the toroidal magnetic fields and the associated radial poloidal electric currents at two distinct altitudes in the ionosphere F region. Clear signatures of low‐ and mid‐latitude interhemispheric field‐aligned currents (IHFACs) are identified. The model reproduces well‐known characteristics of the climatology of IHFACs and provides new insights, for example, on the average daily variations of IHFACs during winter in the Northern Hemisphere. It also well recovers the variations of IHFACs with longitude. The potential driving mechanisms of these variations, such as longitudinal variations of the main field and modulation by upward propagating atmospheric tides, are discussed. The new model can be used to analyze the relationship between atmospheric tides and IHFACs. It can also be used to investigate the connection between the magnetic fields and electric currents from the ionospheric E and F regions in order to improve the separation of these fields as well as our understanding of the overall ionospheric electric current system.

PHYSICAL EFFECTS OF THE POWERFUL TONGA VOLCANO EXPLOSION IN THE EARTH – ATMOSPHERE – IONOSPHERE – MAGNETOSPHERE SYSTEM ON JANUARY 15, 2022

PHYSICAL EFFECTS OF THE POWERFUL TONGA VOLCANO EXPLOSION IN THE EARTH – ATMOSPHERE – IONOSPHERE – MAGNETOSPHERE SYSTEM ON JANUARY 15, 2022

Ionospheric Currents in the Equatorial and Low Latitudes of Africa

The magnetometer data obtained for 2008 from geomagnetic stations installed across Africa by magnetic data acquisition set (MAGDAS) have been used to study the ionospheric Sq current system in the equatorial and lowlatitudes of Africa. The aim of this work is to separate the quiet-day feld variations obtained in the equatorial and low latitude regions of Africa into their external and internal feld contributions and then to use the paired external and internal coeffcients of the SHAto determine the source current and induced currents. The method used involved a spherical harmonic analysis (SHA). This was applied in the separation of the internal and external field/current contribution to the Sq variations. The result shows that the variation in the currents is seen to be a dawn-to-dusk phenomenon with the variation in the external currents different from that of the internal currents both in amplitude and in phase. Furthermore, the seasonal variation in the external current maximizes during the March equinox and minimizes during the December solstice. The maximum current observed in AAB and ILR is due to the Equatorial Electrojet Current present in the AAB and ILR stations. Seasonal variation was observed in the geomagnetic component variations as well as in the currents. This is attributed to the position of the sun with respect to the earth at different months of the year. The equinoctial maximum is observed in external current intensity which occurred mostly during the March Equinox.

Global Characteristics of Improved Interhemispheric Field‐Aligned Currents and of F‐Region Meridional Currents Observed by the Swarm Dual‐Spacecraft

AbstractThe Swarm satellite constellation provides an excellent opportunity to explore ionospheric current systems. In this study, we have reanalyzed the ionospheric currents derived by the Swarm dual‐satellite approach, to investigate the characteristic of inter‐hemispheric field‐aligned currents (IHFAC). One major improvement is that the influence of the ambient magnetic field on the IHFAC intensity has been taken into account, and the derived IHFAC densities are normalized to their ionospheric E‐region footprints. In addition, we have extended the analysis to the middle latitudes within ±60° MLat, which show IHFAC features different from those at low latitudes. For the first time the tidal features at mid‐latitudes are studied. They are dominated by longitudinal wave‐1 and wave‐2 patterns. A superposition of these tidal components reflects a confinement of the IHFAC tidal modulation to daytime in the western hemisphere. The tidal signatures at middle latitudes are better organized in universal time than in local time. The strongest IHFACs are appearing around 18–19 UT, near noon in the American sector. This is related to the overlap with the South Atlantic Anomaly. For the first time we define and analyze the F‐region meridional currents (FMC), interpreting the remaining part of the ionospheric radial current (IRC), which is not covered by IHFAC. The mean FMCs are found to be highly symmetric with respect to the magnetic equator and only amounts to one‐third or half of the IHFAC density, however, they show neither typical longitudinal, diurnal variations, or tidal signatures.

Signatures of wedgelets over Fennoscandia during the St Patrick’s Day Storm 2015

During the long main phase of the St Patrick’s Day storm on March 17, 2015, we found three separate enhancements of the westward electrojet. These enhancements are observed in the ionospheric equivalent currents computed using geomagnetic data over Fennoscandia. Using data from the IMAGE magnetometer network, we identified localised field-aligned current (FAC) systems superimposed on the pre-existing ionospheric current system. We suggest that these localised current systems are wedgelets and that they can potentially contribute to a larger-scale structure of a substorm current wedge (SCW). Each wedgelet is associated with a negative BXspike. Each spike is recorded at a higher latitude than the former one and all three are very localised over Fennoscandia. The first spike occurred at 17:34 UT and was observed at Lycksele, Rørvik and Nurmijärvi, the second spike was recorded at 17:41 UT and located at Lycksele and Rørvik, whereas the last spike occurred at 17:47 UT and was observed at Kevo and Abisko. Simultaneous optical auroral data and electron injections at the geosynchronous orbit indicate that one or more substorms took place in the polar ionosphere at the time of the wedgelets. This study demonstrates the occurrence of small and short-lived structures such as wedgelets at different locations over a short time scale, 15 min in this case.

The time derivative of the geomagnetic field has a short memory

Abstract. Solar eruptions and other types of space weather effects can pose a hazard to the high voltage power grids via geomagnetically induced currents (GICs). In worst cases, they can even cause large-scale power outages. GICs are a complex phenomenon, closely related to the time derivative of the geomagnetic field. However, the behavior of the time derivative is chaotic and has proven to be tricky to predict. In our study, we look at the dynamics of the geomagnetic field during active space weather. We try to characterize the magnetic field behavior, to better understand the drivers behind strong GIC events. We use geomagnetic data from the IMAGE (International Monitor for Auroral Geomagnetic Effect) magnetometer network between 1996 and 2018. The measured geomagnetic field is primarily produced by currents in the ionosphere and magnetosphere, and secondarily by currents in the conducting ground. We use the separated magnetic field in our analysis. The separation of the field means that the measured magnetic field is computationally divided into external and internal parts corresponding to the ionospheric and telluric origin, respectively. We study the yearly directional distributions of the baseline subtracted, separated horizontal geomagnetic field, ΔH, and its time derivative, dΔH/dt. The yearly distributions do not have a clear solar cycle dependency. The internal field distributions are more scattered than the external field. There are also clear, station-specific differences in the distributions related to sharp conductivity contrasts between continental and ocean regions or to inland conductivity anomalies. One of our main findings is that the direction of dΔH/dt has a very short “reset time“, around 2 min, but ΔH does not have this kind of behavior. These results hold true even with less active space weather conditions. We conclude that this result gives insight into the time scale of ionospheric current systems, which are the primary driver behind the time derivative's behavior. It also emphasizes a very short persistence of dΔH/dt compared to ΔH, and highlights the challenges in forecasting dΔH/dt (and GIC).

Q-Bursts Generated by the Contact of the Interplanetary Shock Wave with the Magnetosphere

The study presents the results of the global registration of a short-term train of damped oscillations in the frequency range 0,5–30,0 Hz (Q-burst), observed at the initial phase of the preliminary impulse of the sudden onset SC on 09.09.2011. The spatial features of the spectrum of train oscillations are studied. We discovered A spectral-resonance structure of bursts, which is observed globally. It is assumed that bursts are generated as a result of powerful lightning discharges that create red sprites in the ionosphere. The electromagnetic radiation of sprites is captured in the earth-ionosphere waveguide and propagates globally. We assumed that one of the reasons for the increase in the power of the Q-signal is the coincidence of the beginning of its generation with the amplification of the ionospheric current system of the preliminary impulse of the sudden onset of SC.

Pressure-Gradient Current at High Latitude from Swarm Measurements

The pressure-gradient current is among the weaker ionospheric current systems arising from plasma pressure variations. It is also called diamagnetic current because it produces a magnetic field which is oriented oppositely to the ambient magnetic field, causing its reduction. The magnetic reduction can be revealed in measurements made by low-Earth orbiting satellites flying close to ionospheric plasma regions where rapid changes in density occur. Using geomagnetic field, plasma density and electron temperature measurements recorded on board ESA Swarm A satellite from April 2014 to March 2018, we reconstruct the flow patterns of the pressure-gradient current at high-latitude ionosphere in both hemispheres, and investigate their dependence on magnetic local time, geomagnetic activity, season and solar forcing drivers. Although being small in amplitude these currents appear to be a ubiquitous phenomenon at ionospheric high latitudes characterized by well defined flow patterns, which can cause artifacts in the main field models. Our findings can be used to correct magnetic field measurements for diamagnetic current effect, to improve modern magnetic field models, as well as to understand the impact of ionospheric irregularities on ionospheric dynamics at small-scale sizes of a few tens of kilometers.

Natural Orthogonal Component Analysis of Daily Magnetic Variations at the Martian Surface: InSight Observations

AbstractDistinguishing different sources of magnetic field variations at InSight is important to understand dynamic processes in the Martian ionosphere as well as the coupling between the solar wind and the Martian induced magnetosphere. Recent studies based on magnetic field measurements from InSight have suggested that the daily and seasonal variations in the magnetic field at the Martian surface are at least partially the result of neutral wind‐driven ionospheric dynamo currents and their seasonal variations. However, the sources of the daily variations with different time scales need be further investigated. In this paper, magnetic field variations in a sol as well as during nearly a whole Martian year from InSight observations were decomposed into their natural orthogonal components. We found that the first eigenmode shows the previously identified early to midmorning peak, and varies with season. This corresponds to the solar quiet variations. The second and higher eigenmodes manifest the quasi‐Carrington and sub‐Carrington rotation periodicity represent disturbed components that may be stimulated by variations in the draped interplanetary magnetic field and/or the Martian ionospheric electron density. Different from their counterparts at the Earth, the amplitude of the first eigenmode is comparable with the sum of second to fifth ones, showing that the quiet and disturbed diurnal variations contribute similarly to the total diurnal variation. Decomposition of Martian surface magnetic field variations could provide monitoring of the Martian ionospheric current system as well as the solar wind conditions in the near‐Mars space, which will be greatly enhanced when combined with Zhurong Martian surface field measurements in the future.