Assessment of International GNSS Service Global Ionosphere Map products over China region based on measurements from the Crustal Movement Observation Network of China

The global ionosphere maps(GIM)provided by the International GNSS Service(IGS)are extensively utilized for ionospheric morphology monitoring,scientific research,and practical application.Assessing the credibility of GIM products in data-sparse regions is of paramount importance.In this study,measurements from the Crustal Movement Observation Network of China(CMONOC)are leveraged to evaluate the suitability of IGS-GIM products over China region in 2013-2014.The indices of mean error(ME),root mean square error(RMSE),and normalized RMSE(NRMSE)are then utilized to quantify the accuracy of IGS-GIM products.Results revealed distinct local time and latitudinal dependencies in IGS-GIM errors,with substantially high errors at nighttime(NRMSE:39%)and above 40°latitude(NRMSE:49%).Seasonal differences also emerged,with larger equinoctial deviations(NRMSE:33.5%)compared with summer(20%).A preliminary analysis implied that the irregular assimilation of sparse IGS observations,compounded by China’s distinct geomagnetic topology,may manifest as error variations.These results suggest that modeling based solely on IGS-GIM observations engenders inadequate representations across China and that a thorough examination would proffer the necessary foundation for advancing regional total electron content(TEC)constructions.

The anisotropy of suprathermal electrons in the Martian ionosphere

Suprathermal electrons are an important population of the Martian ionosphere, either produced by photoionization of atmospheric neutrals or supplied from the Solar Wind (SW). This study is dedicated to an in-depth investigation of the pitch angle distribution of suprathermal electrons at two representative energies, 19−55 eV and 124−356 eV, using the extensive measurements made by the Solar Wind Electron Analyzer on board the Mars Atmosphere and Volatile Evolution. Throughout the study, we focus on the overall degree of anisotropy, defined as the standard deviation of suprathermal electron intensity among different directions which is normalized by the mean omni-directional intensity. The available data reveal the following characteristics: (1) In general, low energy electrons are more isotropic than high energy electrons, and dayside electrons are more isotropic than nightside electrons;(2) On the dayside, the anisotropy increases with increasing altitude at low energies but remains roughly constant at high energies, whereas on the nightside, the anisotropy decreases with increasing altitude at all energies;(3) Electrons tend to be more isotropic in strongly magnetized regions than in weakly magnetized regions, especially on the nightside. These observations indicate that the anisotropy is a useful diagnostic of suprathermal electron transport, for which the conversion between the parallel and perpendicular momenta as required by the conservation of the first adiabatic invariant, along with the atmospheric absorption at low altitudes, are two crucial factors modulating the observed variation of the anisotropy. Our analysis also highlights the different roles on the observed anisotropy exerted by suprathermal electrons of different origins.

The nature of electron density enhancement over a wide altitude range during ionosphere heating experiments at EISCAT

During the course of ionospheric heating experiments, researchers at the European Incoherent Scatter Scientific Association (EISCAT) observed an apparent electron density enhancement. The enhancement extended over a wide range of altitudes, above the reflection altitude of the high-frequency pump wave. However, whether this enhancement actually corresponds to a true enhancement in electron density remains an open question. When the dispersion relation of ion acoustic waves is followed, the frequency ratio of the enhanced ion line to the background ion line suggests that the profile of the effective ion mass may have remained unchanged. Furthermore, the solar radio flux and ion drift velocity indicate no significant changes in the ion species and their densities. In conclusion, the electron density enhancement observed at EISCAT should not, in fact, be considered a true enhancement.

Automatic recognition of tweek atmospherics and plasma diagnostics in the lower ionosphere with the machine learning method

Tweek atmospherics are extremely low frequency and very low frequency pulse signals with frequency dispersion characteristics that originate from lightning discharges and that propagate in the Earth–ionosphere waveguide over long distances.In this study,we developed an automatic method to recognize tweek atmospherics and diagnose the lower ionosphere based on the machine learning method.The differences(automatic−manual)in each ionosphere parameter between the automatic method and the manual method were−0.07±2.73 km,0.03±0.92 cm^(−3),and 91±1,068 km for the ionospheric reflection height(h),equivalent electron densities at reflection heights(Ne),and propagation distance(d),respectively.Moreover,the automatic method is capable of recognizing higher harmonic tweek sferics.The evaluation results of the model suggest that the automatic method is a powerful tool for investigating the long-term variations in the lower ionosphere.

Progress of Geodesy Related Ionosphere from Chinese Scientists in the Period of 2019—2023

The ionosphere is the ionized part of the upper atmosphere of the Earth,which plays an important role in atmospheric electricity and forms the inner edge of the magnetosphere.It influences radio propagation significantly,such as the Global Navigation Satellite System(GNSS).Meanwhile,the GNSS is also an essential technique for sensing the variation of ionosphere.During the years of 2019—2023,a large number of Chinese geodesy scientists devoted much efforts to the geodesy related ionosphere.Due to the very limited length,the achievements are carried out from the following six aspects,including:①The ionospheric correction models for BDS and BDSBAS;②Real-time global ionospheric monitoring and modeling;③The ionospheric 2D and 3D modeling based on GNSS and LEO satellites;④The ionospheric prediction based on artificial intelligence;⑤The monitoring and mitigation of ionospheric disturbances for GNSS users;⑥The ionospheric related data products and classical applications.

Statistical Characteristics of the Temporal Spectrum of Scattered Radiation in the Equatorial Ionosphere

On the basis of the solution of the space-time characteristic system by the method of geometric optics using symbolic calculations,analytical and numerical simulation of the propagation of the ordinary and extraordinary radio waves in the conducting equatorial ionospheric plasma was made considering the anisotropy of plasma irregularities and non-stationary nature of propagation medium.Broadening of the spectrum and the displacement of its maximum contain velocity of a turbulent plasma flow and parameters characterizing anisotropic plasmonic structures.Statistical moments of both radio waves do not depend on the absorption sign and are valid for both active and absorptive random media.Temporal pulsations and conductivity of a turbulent ionospheric plasma have an influence on the evaluation of the spectrum-varying propagation distances travelling by these waves.The new double-humped effect in the temporal spectrum has been revealed for the ordinary wave varying anisotropy coefficient and dip angle of stretched plasmonic structures.From a theoretical point of view,the algorithms developed in this work allow effective modelling of the propagation of both radio signals in the equatorial conductive ionospheric plasma,considering the external magnetic field,inhomogeneities of electron density in-homogeneities,as well as non-stationary.

Manifestations of the earthquake preparations in the ionosphere total electron content variations

Recent ionospheric observations report anomalous total electron content (TEC) deviations prior strong earthquakes. We discuss common fetures of the pre-earthquake TEC disturbances on the basis of statistics covering 50 strong seismic events during 2005-2006. The F2-layer ionospheric plasma drift under action of the electric fields of seismic origin is proposed as the main reason of producing TEC anomalies. The origin of such electric fields is discussed in terms of the lithosphere-atmosphere-ionosphere coupling system. This theory is supported by numerical simulations using global Upper Atmosphere Model (UAM). UAM calculations show that the vertical electric current with the density of about 20 - 40 nA/m2 flowing between the Earth and ionosphere over an area of about 200 by 2000 km is required to produce the TEC disturbances with the amplitude of about 30% - 50% relatively to the non-disturbed conditions. Ionosphere responses on the variations of the latitudinal position, direction and configuration of the vertical electric currents have been investigated. We show that not only the vertical component of the ionospheric plasma drift but also horizontal components play an important role in producing pre-earthquake TEC disturbances.

Heat-induced temperature enhancement in the polar summer ionosphere

Observation data recorded by the European Incoherent Scatter Scientific Association in TromsФ, Norway in August 2009 were analyzed to determine the heating effects in polar summer ionospheric modification experiments. There are two types of increases in electron temperature： large relative increases in a narrow range near 150 km and greater absolute increases in a wider range at 150-400 km. The percentage increase in temperature linearly increases with heating power, but the rate of increase decreases with increasing pump frequency. A clear two-dimensional distribution was found for the measurement made on August 15, and the beating effects are greater closer to the direction of the geomagnetic field. The heating effects obviously depend on the angle between the heating beam and geomagnetic field; as the angle increases, the heating effect decreases.

Penetration characteristics of VLF wave from atmosphere into lower ionosphere

The factors affecting the reflection and transmission coefficient of the ionosphere have been analyzed.These factors include wave frequency,incident angle,geomagnetic inclination,electron density and collision frequency in the ionosphere.The ionosphere refractive index is also analyzed.The ionosphere above 70 km is considered to be homogeneous and anisotropic,and the reflection and transmission coefficient matrix is calculated using matrix method.Simultaneously the Booker quartic equation is solved to get the refractive index in the ionosphere.The results show that when the wave frequency is higher,it is easier to penetrate into the ionosphere from its bottom boundary and the propagation attenuation in the ionosphere is smaller.TE（traverse electric） wave and TM（traverse magnetic） wave can both penetrate into the ionosphere with a small incident angle,while TE wave can hardly transmit into the ionosphere when the incident angle is large.The transmission coefficient decreases as the geomagnetic inclination increases.TE and TM wave cannot penetrate into the ionosphere at magnetic equator.When the electron collision frequency is higher,it is easier for VLF wave to penetrate into the ionosphere and the attenuation of ordinary wave is weaker,which may be caused by the energy transportation between the waves and the particles.The ordinary（O） wave experiences severer attenuation than extraordinary（X） wave,and X wave is a penetration mode whereas O wave is a non-penetration mode in the ionosphere.All the results indicate that VLF wave with higher frequency is easier to penetrate into the ionosphere and to be recorded by the satellites at high latitude.It is hard for ULF and the lower frequency VLF wave to transmit into the ionosphere directly for the severe reflection and attenuation.It may transmit into the ionosphere with a small incident angle due to the nonlinear effect,for example,the interaction between the waves and the particles or cross modulation,and then propagate along the whistle duct with small attenuation.This work may be a preliminary theoretical exploration for the future calculation on the response of ground based VLF artificial transmitter in the ionosphere and further study on the seismic ionosphere coupling model.

The earthquake-related disturbances in ionosphere and project of the first China seismo-electromagnetic satellite

Based on the case studies and statistical analysis of earthquake-related ionospheric disturbances mainly from DEMETER satellite, ground-based GPS and ionosounding data, this paper summarizes the statisw tical characteristics of earthquake-related ionospheric disturbances, including electromagnetic emissions, plasma perturbations and variation of energetic particle flux. According to the main results done by Chinese scientists, fusing with the existed study from global researches, seismo-ionospheric disturbances usually occurred a few days or hours before earthquake occurrence. Paralleling to these case studies, lithosphere-atmosphere-ionosphere （LAI） coupling mechanisms are checked and optimized. A thermo-electric model was proposed to explain the seismo-electromagnetic effects before earthquakes. A propagation model was put forward to explain the electromagnetic waves into the ionosphere. According to the requirement of earthquake prediction research, China seismo-electromagnetic satellite, the first space-based platform of Chinese earthquake stereoscopic observation system, is proposed and planned to launch at about the end of 2014. It focuses on checking the LAI model and distinguishing earthquake-related ionospheric disturbance. The preliminary design for the satellite will adopt CAST-2000 platform with eight payloads onboard. It is believed that the satellite will work together with the ground monitoring network to improve the capability to capture seismo-electromagnetic information, which is beneficial for earthquake monitoring and prediction researches.

Novel adaptive Hatch filter to mitigate the effects of ionosphere and multipath on LAAS

It has been proven that carrier smoothing and differential global positioning system （DGPS） are effective to improve the accuracy of pseudorange by reducing the noise in it and eliminating almost all the common mode errors between the ground station and user. However, another issue coming with local area augmentation system （LAAS） is how to find an adaptive smoothing window width to minimize the error on account of ionosphere delay and multipath. Based on the errors analysis in carrier smoothing process, a novel algorithm is formulated to design adaptive Hatch filter whose smoothing window width flexibly varies with the characteristic of ionosphere delay and multipath in the differential carrier smoothing process. By conducting the simulation in LAAS and after compared with traditional Hatch filers, it reveals that not only the accuracy of differential correction, but also the accuracy and the robustness of positioning results are significantly improved by using the designed adaptive Hatch filter.

Online computation of International Reference Ionosphere Extended to Plasmasphere(IRI-Plas) model for space weather

Ionosphere is the most challenging part of Space Weather with its spatio-temporal variability and dispersive nature. Ionospheric models are very important in reducing positioning error in GNSS system.International Reference Ionosphere(IRI) is an empirical, deterministic and climatic model of ionosphere up to 2000 km in height. Recently, IRI Extended to Plasmasphere(IRI-Plas) model has been developed to extend the interest region of IRI to the GPS orbital height of 20,000 km. Both IRI and IRI-Plas provide ionospheric parameters such as electron density, electron and ion temperatures according to their height profiles. In order to update the model to current ionospheric conditions, IRI-Plas can input F2 layer critical frequency(foF2), maximum ionization height(hmF2), and also Total Electron Content(TEC).Online IRI-Plas is developed for the ionospheric community to run multiple tasks at various locations,dates and times with optional foF2, hmF2 and TEC inputs in a user-friendly manner. In this paper, we are going to present the capabilities of the Online IRI-Plas service and provide some comparisons between IRI-Plas outputs and ionosonde measurements. The comparison between online IRI-Plas foF2 outputs and ionosonde foF2 measurements indicates that the model with TEC input can significantly improve the representation of the current ionospheric state, which is very successful especially in the geomagnetically disturbed days.

Temporal distribution characteristics of GNSS ionospheric occultation data and its effects in earthquake-ionosphere anomaly detection

The temporal distribution characteristics of COSMIC occultation data are analyzed in detail, and the limitations in earthquake-ionosphere anomaly detection caused by the temporal distribution characteristics of COSMIC occultation data are discussed using the example of the Wenchuan earthquake. The results demonstrate that there is no fixed temporal resolution for COSMIC occultation data when compared with other ionospheric observation techniques. Therefore, occultation data cannot currently be independently utilized in research studies but can only be used as a complement to other ionospheric observation techniques for applications with high temporal resolution demands, such as earthquake-ionosphere anomaly detection.

An Introduction to Magnetosphere-Ionosphere-Thermosphere Coupling Small Satellite Constellation

A future Chinese mission is introduced to study the coupling between magnetosphere,ionosphere and thermosphere,i.e.the Magnetosphere-Ionosphere-Thermosphere Coupling Small Satellite Constellation(MIT).The scientific objective of the mission is to focus on the outflow ions from the ionosphere to the magnetosphere.The constellation is planning to be composed of four small satellites;each small satellite has its own orbit and crosses the polar region at nearly the same time but at different altitude.The payloads onboard include particle detectors,electromagnetic payloads,auroral imagers and neutral atom imagers.With these payloads,the mission will be able to investigate acceleration mechanism of the upflow ions at different altitudes.Currently the orbits have been determined and prototypes of some have also been completed.Competition for next phase selection is scheduled in late 2015.

Propagation characteristics of power line harmonic radiation in the ionosphere

Theoretical model and solutions on power line harmonic radiation （PLHR） propagating in the ground, air, and anisotropic homogeneous ionosphere are presented, The theoretical model is verified by the PLHR events observed by the DEMETER satellite. Some propagation characteristics of PLHR based on the model are obtained. This paper is bene- ficial to quantitatively interpret the formation mechanism of PLHR phenomenon.

Manifestation of earthquake preparation zone in the ionosphere before 2021 Sonitpur,Assam earthquake revealed by GPS-TEC data

Global Positioning System(GPS)Continuously Operating Reference Station(CORS)data analysis shows that the ionosphere’s electron density variability is linked to the deformation and stress accumulation in the Earth’s crust.Anomalies in ionosphere total electron content(TEC)variability before 2021 M6.4 Sonitpur,Assam earthquake were detected using L1 and L2 GPS frequencies that showed three distinct abnormalities on April 3,9,10,2021.Pearson’s correlation coefficient(r)of TEC decreases in the CORS that lies away from the earthquake epicenter,indicating the possibilities of a positive relationship between TEC variability and earthquake epicenter.TEC concentration also decreases towards the epicenter within the earthquake preparation zone(EPZ).It is also observed that the Pearson’s correlation coefficient(r)of TEC decreases linearly near the EPZ.The study demonstrates the possibilities of determining the TEC anomalous zone in the ionosphere that coincides with the EPZ in the crustal rocks.The research indicated the possibilities of magnitude estimation of an impending earthquake based on the TEC anomalous zone in the ionosphere using closely spaced dense CORS network data.

Magnetic Storm Effects in the Auroral Ionosphere Observed with EISCAT Radar-Two Case Studies

Storm-time changes of main plasma parameters in the auroral ionosphere are analyzed for two intense storms occurring on May 15, 1997 and Sept. 25, 1998, with emphasis on their relationship to the solar wind dynamic pressure and the IMFB z component. Strong hard particle precipitation occurred in the initial phase for both storms, associated with high solar wind dynamical pressure. During the recovery phase of the storms, some strong particle precipitation was neither concerned with high solar wind pressure nor southward IMFB z. Severe negative storm effects depicted by electron density depletion appeared in theF-region during the main and recovery phase of both storms, caused by intensive electric field-related strong Joule/frictional heating when IMF was largely southward. The ion temperature behaved similarly inE-andF-region, but the electron temperature did quite different, with a strong increase in the lowerE-region relating to plasma instability excited by strong electric field and a slight decrease in theF-region probably concerning with a cooling process. The field-aligned ion velocity was high and apparently anticorrelated with the northward component of the ion convection velocity.

EISCAT Observations of Main Ionization Troughs in the High-Latitude Ionosphere

F-region electron density depletions associated with main ionization troughs in the high-latitude ionosphere are studied using EISCAT CP3 data of meridian scanning experiments. The troughs in our observations are found to appear mainly in dusk sector, extending from late afternoon to pre-midnight, with higher occurrence rate during equinox and winter. Simultaneous ion drift velocity inF-region shows that the main trough minimum is mostly located at the equator ward edge of the plasma convection flow, rather than in the region where the largest ion flow are observed.

Comparison of Diurnal,Seasonal and Solar Cycle Variations of High-latitude,Mid-latitude and Low-latitude Ionosphere

Comparison of regular(diurnal,seasonal and solar cycle)variations of high-latitude,mid-latitude and low-latitude ionospheric characteristics has been provided on basis of local empirical models of the peak electron density and the peak height.The local empirical models were derived from the hand-scaled ionogram data recorded by DPS-4 digisondes located at Norilsk(69°N,88°E),Irkutsk(52°N,104°E)and Hainan(19°N,109°E)for a 6-year period from December,2002 to December,2008.The technique used to build the local empirical model is described.The primary focus is diurnal-seasonal behavior under low solar activity and its change with increasing solar activity.Both common and specific features of the high-latitude(Norilsk),mid-latitude(Irkutsk)and low-latitude(Hainan)regular variations were revealed using their local empirical models.

On ionosphere-delay processing methods for single-frequency precise-point positioning

In single-frequency precise-point positioning of a satellite,ionosphere delay is one of the most important factors impacting the accuracy. Because of the instability of the ionosphere and uncertainty of its physical properties, the positioning accuracy is seriously limited when using a precision-limited model for correction. In order to reduce the error, we propose to introduce some ionosphere parameter for real-time ionosphere-delay estimation by applying various mapping functions. Through calculation with data from the IGS（ International GPS Service） tracking station and comparison among results of using several different models and mapping functions, the feasibility and effectiveness of the new method are verified.