Region Of Ionosphere Research Articles

ACCURACY ASSESSMENT OF PRECISE POINT POSITIONING WITH MULTI-CONSTELLATION GNSS DATA UNDER STRONG SOLAR BURST EFECTS

Solar variations modify a layer of the Earth’s upper atmosphere known as the ionosphere. This is of particular concern for the aviation sector because of the way its communications and navigation systems can be affected. At the same time, one of the most complex atmospheric effects is the response of ionospheric regions to geomagnetic storms. The ionospheric response during the same storm can vary in time in different locations, which can introduce significant errors/displacement (meters) in single-frequency relative GNSS positioning (DGNSS technology). The residual effect can be somewhat mitigated by using dual- or multi-frequency GNSS, but dual frequency is not a guarantee against degradation of relative observations results, especially during significant geomagnetic storms. In this regard, PPP absolute positioning technology can be effective. However, another atmospheric effect – ionospheric scintillation can have a significant impact on the accuracy of both GNSS positioning approaches. The main goal of this study was to analyze the effect of second-order ionospheric delay during geomagnetic storms and ionospheric scintillations on GNSS positioning using the PPP method. GNSS data corrected and uncorrected for higher-order ionospheric delay, respectively, were processed by the static PPP-AR method using the PRIDE-PPPAR ver.2.2.6 software for the selected periods of geomagnetic storms. From the analysis of the influence of second-order ionospheric errors, it follows that their values can reach almost 4 cm for first-frequency signals under different states of ionospheric disturbances for the GPS constellation and almost an order of magnitude less for the GNSS quadroconstellation. The appearance of stronger geomagnetic storms increases the second-order ionospheric errors by several millimeters.

A ThingSpeak IOT on Real-Time Ionospheric Monitoring System

The primary factor for the ionospheric monitoring system that we take into consideration in this paper is the Total electron content (TEC). This TEC value varies over time, space, and as the ionospheric area absorbs solar heat. The ThingSpeak ionospheric monitoring system can be put into practice by taking into saved at several TEC receiver locations, such as Port Blair and Bengaluru. Along with the information on solar index values, these stations also have data on geometric index values. Additionally, all of the TEC values are taken into account and examined to develop the necessary monitoring system for the ionospheric region. By doing this, we will be able to check the report by using IOT to connect it to ThingSpeak, which can then be utilized or checked from anywhere using the API key.

Optimal Estimation Inversion of Ionospheric Electron Density from GNSS-POD Limb Measurements: Part I-Algorithm and Morphology

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in characterizing the global 3D distribution and variability of ionospheric electron density (Ne). In this study, we developed an optimal estimation (OE) method to retrieve Ne profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles (ε < 0°). Although both OE and onion-peeling (OP) methods use the Abel weighting function in the Ne inversion, they are significantly different in terms of performance in the lower ionosphere. The new OE results can overcome the large Ne oscillations, sometimes negative values, seen in the OP retrievals in the E-region ionosphere. In the companion paper in this Special Issue, the HmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006–2019), COSMIC-2 (2019–present), and Spire (2019–present) constellations, showing a consistent ionospheric Ne morphology. The unprecedented spatiotemporal sampling of the ionosphere from these constellations now allows a detailed analysis of the frequency–wavenumber spectra for the Ne variability at different heights. In the lower ionosphere (~150 km), we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2, and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4, and SW4. The co-existence of eastward- and westward-propagating wave4 components implies the presence of a stationary wave4 (SPW4), as suggested by other earlier studies. Further improvements to the OE method are proposed, including a tomographic inversion technique that leverages the asymmetric sampling about the tangent point associated with GNSS-LEO links.

Method of Development of a New Regional Ionosphere Model (RIM) to Improve Static Single-Frequency Precise Point Positioning (SF-PPP) for Egypt Using Bernese GNSS Software

Due to the lack of coverage of IGS in Africa, especially over North Africa, and the construction revolution of infrastructure in Egypt, a geodetic CORS stations network was established in 2012. These CORS stations are operated by the Egyptian Surveying Authority (Egy. SA) and cover the whole of Egypt. The paper presents a fully developed regional ionosphere model (RIM) depending on the Egyptian CORS stations. The new model and the PPP solution were obtained using Bernese GNSS V. 5.2 software. An observation data series of eight days (DOY 201–208)/2019 was used in this study. Eighteen stations were used to develop the RIM model for each day; fifteen stations were used to validate the new RIM model. A static SF-PPP solution was obtained using the CODE-GIM and RIM models. Comparing the outcomes to the reference network solution, based on the recently developed RIM model, the solution showed a mean error of 0.06 m in the East direction, 0.13 m in the North direction, and 0.21 m in the height direction. In the East, North, and height directions, this solution improves the SF-PPP result achieved by the Global Ionosphere Maps (CODE-GIM) model by 60%, 68%, and 77%, respectively.

Variability of Vertical Total Electron Content Gradients in the Brazilian Sector

Ground Based Augmentation Systems are designed to meet demanding requirements during airport approach and landing sections of aircraft routes. It computes corrected and smoothed differential corrections and other information. These messages are transmitted to aircraft, which may use the corrections to improve their position estimation. However, residual correction errors remain, due to ionospheric spatial and temporal gradients between Ground Based Augmentation System reference receivers and approaching aircrafts. Large vertical ionospheric delay gradients detected during combinations of extreme geophysical environments indicate that there are problematic conditions for Ground Based Augmentation System operation over the Brazilian region. Motivated by the above conditions, data from Rede Brasileira de Monitoramento Contínuo, the Brazilian public network of dual-frequency Global Navigation Satellite System receivers operated by Instituto Brasileiro de Geografia e Estatística were analyzed to estimate ionospheric vertical Total Electron Content gradients. The effects from different locations (represented by the dip latitude) and geophysical conditions (represented by solar activity, season, and local time) on variations of ionospheric vertical Total Electron Content gradients have been statistically quantified. The results from the present work provide information to more detailed performance assessments of Ground Based Augmentation System procedures under equatorial and low-latitude ionospheric regions, aimed at flexibly delimiting the conditions for feasible operation in such regions.

On the nonlinear dynamics of large scale dust-acoustic solitary waves in a superthermal bi-ion dusty ionospheric magnetoplasma

The purpose of this study is to examine the properties of the dust-acoustic solitary waves in a complex magnetoplasma made up of negatively charged moving dust grains in the lower ionospheric region and inertialess electrons and ions obeying Maxwell and kappa distributions, respectively. In this context, the reductive perturbation technique is carried out to obtain the Zakharov–Kuznetsov (ZK) equation within the given framework. The obtained evolution equation, i.e., ZK equation is transformed to a planar dynamical system for studying the qualitative behavior of the solitary waves. The impact of important physical parameters, such as the dust number density, ion concentration, ion temperature, superthermality, and the background magnetic field, on the profile of the nonlinear structures is numerically investigated. The findings may be utilized to comprehend the low-frequency irregularities that are detected in the lower ionosphere.

The role of the storm-time prompt penetrating electric field on the net distribution of plasma density over the low latitude ionospheric regions

The relative role of the prompt penetration of electric fields and neutral winds in the distribution of ionospheric plasma during geomagnetic storms over the Indian ionosphere region is studied in this paper. Four geomagnetic storms with similar onset times but of varying strengths have been analyzed with the help of data from ground-based magnetometers, and the ACE, GPS, and TIMED satellites. The strength and magnitude of the prompt penetration electric field were inferred from the magnetic field variations of the horizontal component of the magnetic field (ΔH) at Tirunelveli, a dip-equatorial station, using the ground-based magnetometers. The electron density distribution over the Indian ionospheric region was quantified using GPS-derived vertical total electron content (VTEC). Signatures of prompt penetration electric field were present for all the storms considered whose main phase was during the daytime. It is noted that as the strength of the penetration electric field increases, the latitudinal extent to which the electron density gets distributed also increases. The gradient in electron density extends even up to Shimla (a mid-latitude station) during a severe geomagnetic storm. The gradient was less for moderate geomagnetic storms. Since during the geomagnetic storm, a factor that controls electron density modulation is the neutral wind, its contribution is analyzed with the help of the thermospheric O/N2 ratio, measured by the global ultraviolet imager (GUVI) instrument, onboard NASA’s TIMED satellite. The observed increase in the O/N2 ratio does not correspond with the increase in the VTEC. We surmise that during the first day of the storm, the penetrating electric field is an essential factor that controls the electron density distribution in the low-latitude region. Later the plasma density distribution is controlled mainly by the thermospheric wind.

Forecasting Equatorial Ionospheric Convective Instability With ICON Satellite Measurements

AbstractMeasurements from the Ionospheric Connections Explorer satellite (ICON) form the basis of direct numerical forecast simulations of plasma convective instability in the postsunset equatorial F region ionosphere. ICON data are selected and used to initialize and force the simulations and then to test the results one orbit later when the satellite revisits the same longitude. Data from the IVM plasma density and drifts instrument and the MIGHTI red‐line thermospheric winds instrument are used to force the simulation. Data from IVM are also used to test for irregularities (electrically polarized plasma depletions). Fourteen datasets from late March 2022, were examined. The simulations correctly predicted the occurrence or non‐occurrence of irregularities 12 times while producing one false positive and one false negative. This demonstrates that the important telltales of instability are present in the ICON state variables and that the important mechanisms for irregularity formation are captured by the simulation code. Possible refinements to the forecast strategy are discussed.

The effects of atmospheric dust and solar radiation on the dayside ionosphere of Mars derived from 17 years of Mars Express radio science observations

This work combines 17 years of Mars Express radio science (MaRS) observations with proxies for insolation and local/global atmospheric dust to investigate the combined and individual effects on the dayside ionosphere of Mars from the top down to the ionospheric base.The increase in insolation from orbital apocenter to pericenter in combination with Mars‘ dust cycle causes an average rise of the whole photochemically dominated region of the dayside ionosphere, ranging from 13 km at the ionospheric base up to 22 km above the main peak during conditions without a global dust storm. The declining phase of the 2018 global dust storm was observed by MaRS on the southern hemisphere and close to pericenter. The observed lifting effect on the whole photochemically dominated region of the ionosphere from the increased insolation and the high local and global atmospheric dust levels exceeds that seen by MaRS from similar seasons during years without a global dust storm.The average ionospheric peak altitude at the subsolar point rises for increasing levels of local atmospheric dust until a maximum elevation is reached. This maximum depends on the available insolation at the top of the planetary atmosphere. Further increases of the local atmospheric dust levels do not lead to a further rise of the average ionospheric peak altitude in the investigated data set. This indicates a limit for the warming/expansion of the lower neutral atmosphere and the consecutive lifting of the ionosphere based on the available insolation and explains why regional dust storms can cause a similar lifting of the ionospheric main peak region as global dust storms.

Statistical analysis of the regional and global ionosphere model on intense geomagnetic storm

Statistical analysis of the regional and global ionosphere model on intense geomagnetic storm

Troubles in Mars Ionosphere Modeling

AbstractThe Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft has acquired comprehensive measurements of the properties that determine the density and composition of the ionosphere of Mars, namely solar irradiance, thermospheric density and composition, and electron temperature. Ionospheric models that use these measured properties as inputs ought, in principle, to accurately reproduce corresponding MAVEN measurements of ionospheric density and composition. Here we test whether or not several Mars ionosphere models do so. We focus on the optically thin portion of the photochemical equilibrium region of the dayside ionosphere, 160–200 km altitude, where ionospheric modeling is most straight‐forward. The predictions of Fox et al. (2021, https://doi.org/10.1016/j.icarus.2020.114186) clearly delivered better model‐data agreement than the other models considered herein, although the model of Mayyasi et al. (2019, https://doi.org/10.1029/2019ja026481) also performed well. For the other models assessed herein, model‐data discrepancies in ion densities of at least a factor of two were present in at least one major ion species.

Simulating the Upper Hybrid Instability as a Cause for 150 km Echoes

AbstractOne Hundred and Fifty kilometer echoes are a type of strong radar echo observed in the valley region of the equatorial ionosphere, whose origin was long standing mystery. Recently, a new upper hybrid (UH) instability theory, driven by high energy photoelectrons, has been proposed to explain most features of 150 km echoes. However, this instability excites high frequency electron modes, whereas radars observe low frequency ion modes. To explain 150 km echoes, the UH instability must ultimately excite ion acoustic modes. This paper describes a set of particle‐in‐cell simulations used to study how photoelectrons interacting with a cold background plasma generate UH waves, and how these drive ion acoustic waves measured by radars. We implement a new electron‐N2 collision algorithm to better model the bump on tail features of the photoelectron distribution in the valley region. These simulations show that photoelectrons drive unstable UH waves that strongly enhance ion acoustic waves. We also show a strong frequency dependence on power in the ion line, which explains observational differences between radars, most notably the lack of echoes at ALTAIR. The photoelectron driven UH instability successfully reproduces most features of 150 km echoes associated with naturally enhanced incoherent scattering.

Comparative Analysis of Global and Regional Ionospheric Responses during Two Geomagnetic Storms on 3 and 4 February 2022

The present work examines the spatial and temporal distribution of positive and negative TEC anomalies on the global and regional scale. To study the local response of the ionosphere, foF2 data from ground ionosonde stations and TEC data from Madrigal and CODE databases have been used. The relative deviation, which also determines the type of TEC response during geomagnetic storms on 3 and 4 February 2022, is considered. In the present study, the regions of positive and negative TEC anomalies and their evolution during storms are examined in detail. As a result of the study, estimates of the following were obtained: (i) the location of the sectors of the polar regions, in where the particle precipitation from the solar wind is observed, (ii) the mid-latitude regions, in which the mechanism of influence of the O/N2 ratio dominates, and (iii) the region around the equator, in which the influence of the electric field dominates. An attempt was made to determine which mechanism of influence of geomagnetic storms on the ionospheric electron density is dominant in different regions. The following main mechanisms are considered: (a) the additional ionization from the particles’ precipitation, (b) the change of the ratio of atomic oxygen (O) to molecular nitrogen (N2) due to the heating of the neutral air, and (c) the influence on the equatorial ionospheric anomaly.

Three-dimensional ionospheric tomography based on compressed sensing

The insufficient number of low-elevation observations is a limitation of the three-dimensional ionospheric computer tomography (CT) based on the global navigation satellite system (GNSS). To solve this problem, accurate prior information on the regional ionosphere must be obtained. However, it is difficult to explicitly and accurately express prior ionospheric information. This study uses compressed sensing (CS) for ionospheric tomography for the first time. Specifically, the electron density obtained from the international reference ionosphere is used to build a dictionary to fully integrate the prior information into the dictionary. Then, the electron density is reconstructed by using the compressive sampling matching pursuit method. Subsequently, the GNSS data of China (Region I) and Europe (Region II) were utilized to validate this proposed method, and the results are compared with ionosonde observations. The mean and standard deviation (SD) of the difference with respect to the ionosonde result are 41 and 22 km, respectively. The mean and SD of relative deviation were 16% and 9%, respectively. In Region II, the mean and SD of the deviation between the reversed peak electron density and the result of the ionosonde were 1.9 × 1010 m−3 and 8.1 × 1010 m−3, respectively. The mean and SD of the relative deviation were 3% and 13%, respectively. The mean and SD of the peak height deviation were 33 and 19 km, and the mean and SD of the relative deviation were 11% and 7%. The electron density distribution and variation in these two regions showed a local time dependence, and the horizontal gradient of the electron density in the latitude was greater than that in the longitude. Moreover, CT by CS is efficient, taking about 6 s per inversion based on an desktop computer with 16 GB RAM and Intel (R) Core (TM) i7-8700 CPU.

MODELING OF SPATIAL-TEMPORAL VARIATIONS OF DYNAMIC AND THERMAL PROCESS PARAMETERS IN GEOSPACE OVER UKRAINE DURING THE MINIMUM OF 24-TH CYCLE OF SOLAR ACTIVITY (2009, 2019)

We have performed the modelling of spatiotemporal variations of parameters of dynamic and thermal processes in ionospheric plasma at the phases of the minimum of the 24-th cycle of solar activity according to the Kharkiv radar of incoherent scattering. The diurnal dependences of parameters of the processes in the ionospheric plasma at altitudes from 210 to 450 km are constructed for typical geophysical periods (vernal and autumn equinoxes, summer and winter solstices). In the paper, the analysis of spatial and temporal variations of parameters of dynamic and thermal processes in the ionosphere is presented. We determined the value of the plasma transfer velocity due to ambipolar diffusion, the density of the full plasma flux, and the flux of charged particles due to ambipolar diffusion, the value of the energy supplied to the electron gas, the density of the heat flux transferred by electrons from the plasmasphere to the ionosphere, as well as the velocity of the equivalent neutral wind, and the meridional component of the neutral wind velocity. We found that weak variations in space weather do not lead to significant changes in spatiotemporal variations of the parameters of dynamic and thermal processes in the ionosphere for most of the studied periods. Quantitative and qualitative characteristics of most of these parameters and their diurnal variations were typical for the considered seasons. On the contrary, the velocity of the equivalent neutral wind changed significantly (up to 2—2.5 times) even with a weak increase in geomagnetic activity. The reasons for such changes may be the strengthening of horizontal thermospheric winds and the penetration of zonal magnetospheric electric fields into midlatitudes during the equinoxes. The obtained results of calculations can be used in basic studies of solar-terrestrial relations and geospace, for the solution of applied problems related to the ability to predict the state of space weather, as well as for further development of the regional ionosphere model CERIM IION. The object of research: physical processes in the ionospheric plasma. The subject of research: spatiotemporal dependences of the main parameters of ionospheric plasma, which were obtained using incoherent scattering radar. Research methods — terrestrial radiophysical method of incoherent scatter of radio waves, statistical analysis of observation results, semi-empirical modelling of parameters of dynamic and thermal processes.

Near real-time global ionospheric total electron content modeling and nowcasting based on GNSS observations

For the purposes of routinely providing reliable and low-latency Global Ionosphere Maps (GIMs), a method of estimating hourly updated near real-time GIM with a time latency of about 1–2 h based on a 24-h data sliding window of Global Navigation Satellite System (GNSS) near real-time observations and real-time data streams was presented. On the basis of the implementation of near real-time GIM estimation, an hourly updated GIM nowcasting method was further proposed to improve the accurate of short-term total electron content (TEC) prediction. We estimated the Shanghai Astronomical Observatory near real-time GIM (SHUG) and nowcasting GIM (SHPG) in the solar relatively active year (2014) and quiet year (2021), and employed GIMs provided by the International GNSS Service, the Global Positioning System (GPS) differential slant TECs (dSTECs) extracted from global independent GNSS stations, and the vertical TECs (VTECs) inverted from satellite altimetry as the references to validate the estimated results. The GPS dSTECs evaluation results show that SHUG behaves fairly consistent with the rapid GIMs, with a discrepancy of less than 1 TEC unit (TECu) overall. The standard deviations (STDs) of SHUG with respect to Jason-2/-3 VTECs are no more than 10% over the majority of rapid GIMs due to the instability of observations. The performance of 1-h nowcasting SHPG is significantlybetter than the Center for Orbit Determination in Europe (CODE) 1-day predicted GIM (C1PG). GPS dSTEC validation results indicate that 1-h nowcasting SHPG is 1 to 2 TECu more reliable than C1PG in eventful ionospheric electron activity regions, and it outperforms the C1PG by 10% overall versus Jason-2/-3 VTECs. The hourly updated SHUG and SHPG have relatively high reliability and low time latency, and thus can provide excellent service for (near) real-time users and offer more accurate TEC background information than daily predicted GIM for real-time GIM estimation.

A Regional Ionospheric Storm Forecasting Method Using a Deep Learning Algorithm: LSTM

AbstractAn ionospheric storm forecasting method was proposed using a deep learning algorithm, LSTM (long short‐term memory). We used the perturbation index to denote the level of an ionospheric storm, deduced from foF2 data, and helped to remove most of the local time and seasonal variations in the ionosphere. In constructing the model, a number of correlated factors were used as inputs, including the properties of coronal mass ejections, solar flare bursts, interplanetary conditions, and geomagnetic and ionospheric states, and the output was whether an ionospheric storm occurred locally in the next 24 hr. Data sets from 2007 to 2014 were used to train the model, and those from 2015 to 2016 were used for validation. The results showed that the model behaved well in most events. The mean precision rate, recall rate, accuracy, and F1 score of the model were 71.7%, 59.7%, 92.7%, and 65.0% in northern China and 78.9%, 56.3%, 96.3% and 65.0% in southern China, respectively. The LSTM forecasting model performed better than other models such as persistence, multiple‐layer perceptron and support vector machine models. Case studies also showed good performance during geomagnetic storms of different strengths. We believe that this model can be beneficial for functional ionospheric storm operation.

Modeling of the cold electron plasma density for radiation belt physics

This review focusses strictly on existing plasma density models, including ionospheric source models, empirical density models, physics-based and machine-learning density models. This review is framed in the context of radiation belt physics and space weather codes. The review is limited to the most commonly used models or to models recently developed and promising. A great variety of conditions is considered such as the magnetic local time variation, geomagnetic conditions, ionospheric source regions, radial and latitudinal dependence, and collisional vs. collisionless conditions. These models can serve to complement satellite observations of the electron plasma density when data are lacking, are for most of them commonly used in radiation belt physics simulations, and can improve our understanding of the plasmasphere dynamics.

Effects of Energetic Particle Precipitation on Electron Temperature in the E-Region of Ionosphere: an Experimental Study

Effects of Energetic Particle Precipitation on Electron Temperature in the E-Region of Ionosphere: an Experimental Study

Imprint of Mesospheric Quasi 2‐Day Wave in the Ground Geomagnetic Field Variations at Low Latitudes

AbstractThe strong quasi 2‐day wave (Q2DW) event that occurred during January 2015 was studied by employing medium frequency and meteor radar observations of winds from three different locations at low latitudes with considerable longitudinal separation. The large‐amplitude wave activity that occurred in mid‐January is the focus of this study. Using the advantage offered by the longitudinal separation among the stations, the zonal structure of the Q2DW has been examined. The presence of a dominant westward propagating zonal wavenumber 3 (W3) mode was revealed in the analysis. The availability of geomagnetic field data from a number of stations at low latitudes enabled us to observe the wave event at ionospheric dynamo region heights. The phase‐longitude cross‐section of the quasi 2‐day oscillation in the Y component of the geomagnetic field reiterates the presence of the dominant W3 mode observed at lower altitudes. Similar zonal characteristics and periods of the wave signatures found at mesospheric altitudes and in the manifestation of quiet‐time ionospheric currents in geomagnetic field are in conformity with the notion of the wave coupling of the atmosphere‐ionosphere system.