Hybrid modelling the Pioneer Venus Orbiter magnetic field observations

This study presents comparisons between the Pioneer Venus Orbiter (PVO) magnetometer (OMAG) observations and the HYB-Venus hybrid simulation code. The comparisons are made near periapsides of four PVO orbits using the full resolution PVO/OMAG data. Also, the statistics of the solar wind and interplanetary magnetic field (IMF) conditions at Venus are studied using the PVO interplanetary dataset. The statistics include the histograms and the probability density maps of the selected upstream parameters. The confidence intervals derived from the upstream statistics demonstrate the nominal simulation input parameter space. Moreover, the probability density maps give the dependencies between the upstream parameters. The comparisons between the simulation code and the data along the spacecraft trajectory show that the basic, large scale, trends seen in the magnetic field can be understood by the current simulation runs. The discrepancies between the simulation and the data were found to arise at low altitudes close to the planetary ionosphere in the region which cannot be resolved in detail by the grid size of the runs.     

Unusually strong magnetic fields in Titan’s ionosphere: T42 case study

Observations of unusually large magnetic fields in the ionosphere indicate periods of maximum stress on Titan’s ionosphere and potentially of the strongest loss rates of ionospheric plasma. During Titan flyby T42, the observed magnetic field attained a maximum value of 37 nT between an altitude of 1200 and 1600 km, about 20 nT stronger than on any other Titan pass and close to five times greater in magnetic pressure. The strong fields occurred near the corotation-flow terminator rather than at the sub-flow point, suggesting that the flow which magnetized the ionosphere was from a direction far from corotation and possibly towards Saturn. Extrapolation of solar wind plasma conditions from Earth to Saturn using the University of Michigan MHD code predicts an enhanced solar wind dynamic pressure at Saturn close to this time. Cassini’s earlier exits from Saturn’s magnetosphere support this prediction because the Cassini Plasma Spectrometer instrument saw a magnetopause crossing three hours before the strong field observation. Thus it appears that Titan’s ionosphere was magnetized when the enhanced solar wind dynamic pressure compressed the Saturnian magnetosphere, and perhaps the magnetosheath magnetic field, against Titan. The solar wind pressure then decreased, leaving a strong fossil field in the ionosphere. When observed, this strong magnetic flux tube had begun to twist, further enhancing its strength.     

The Venus ionosphere

Physical properties of the Venus ionosphere obtained by experiments on the US Pioneer Venus and the Soviet Venera missions are presented in the form of models suitable for inclusion in the Venus International Reference Atmosphere. The models comprise electron density (from 120 km), electron and ion temperatures, and relative ion abundance in the altitude range from 150 km to 1000 km for solar zenith angles from 0° to 180°. In addition, information on ion transport velocities, ionopause altitudes, and magnetic field characteristics of the Venus ionosphere, are presented in tabular or graphical form. Also discussed is the solar control of the physical properties of the Venus ionosphere.     

Model for induced ionization by galactic cosmic rays in the Earth atmosphere and ionosphere

The ionization profiles produced by galactic cosmic rays in the Earth atmosphere and ionosphere are obtained on the basis of Monte Carlo simulations. Cascade processes in the atmosphere are simulated using CORSIKA 6.52 code with FLUKA 2006 and QGSJET II hadronic interaction subroutines. Proton induced showers are considered using a realistic atmospheric model (US Standard Atmosphere). The energy deposit from different components is taken into account, namely electromagnetic, hadron and muon components. The curvature of the atmosphere is considered in the computer code. On the basis of the computed ionization yield function the ion pair production rate in the atmosphere is obtained for different conditions and locations. The model is applicable to the entire atmosphere, from ground level to upper atmosphere. Several applications of the obtained results are discussed. The Monte Carlo simulation model considers nuclear interactions below the altitude of 35 km. It is compared with analytical–numerical electron production rate model. The latter model which takes into account the electromagnetic interactions above altitudes of 35 km has two main regions of application: above 50 km (thin target model) and between 35 and 50 km (intermediate target model). A good agreement between the CORSIKA results and analytical–numerical model results is found above altitude of 35 km.     

The properties of the low altitude magnetic belt in the Venus ionosphere

When the solar wind dynamic pressure is high, the Venus ionosphere usually contains a belt of steady magnetic field at the very lowest altitudes to which Pioneer Venus probes. The current layer that flows on the high altitude side of this low altitude belt is centered at an altitude which ranges from 170 to 190 km with a most probable altitude of 182 km. This altitude is independent of solar zenith angle and hence the current system is flowing horizontally rather than vertically as proposed by Cloutier and co-workers. The lower edge of the magnetic belt was probed only on the lowest altitude passes of Pioneer Venus. This boundary is even more stable in location. The belt has decayed to 90% of its maximum strength usually by 162 km and to 50% of its maximum strength by 155 km. We interpret these data to indicate that the observed magnetic structure of the Venus ionosphere is a product of temporal evolution rather than of spacecraft motion through a spatially varying static structure.     

Neutral atmosphere and crustal magnetization as factors determining differences in plasma convection and magnetic fields distribution in the ionospheres of Mars and Venus

The pattern of the magnetic field/plasma convection can be, to some extent, recovered from the magnetic field measurements by employing either theoretical or numerical models. We use the MAG/ER day-time measurements of the magnetic field at the altitudes from 90 to 180 km during the elliptical orbits of MGS. Analysis of the altitude variation of the characteristics of the large-scale magnetic fields, which were measured some distance away from strong crustal magnetic anomalies, is summarized. The low density of the Martian atmosphere together with the crustal magnetization result in critical differences in plasma convection which are followed by remarkable differences of the magnetic field features within the ionosphere of Venus and Mars (even in its northern hemisphere where the crustal magnetization is, on the average, low) and distribution of currents.     

Effects of solar activity variations on the low latitude topside nighttime ionosphere

We have studied the topside nighttime ionosphere of the low latitude region using data obtained from DMSP F15, ROCSAT-1, KOMPSAT-1, and GUVI on the TIMED satellite for the period of 2000–2004, during which solar activity decreased from its maximum. As these satellites operated at different altitudes, we were able to discriminate altitude dependence of several key ionospheric parameters on the level of solar activity. For example, with intensifying solar activity, electron density was seen to increase more rapidly at higher altitudes than at lower altitudes, implying that the corresponding scale height also increased. The density increased without saturation at all observed altitudes when plotted against solar EUV flux instead of F10.7. The results of the present study, as compared with those of previous studies for lower altitudes, indicate that topside vertical scale height increases with altitude and that, when solar activity increases, topside vertical scale height increases more rapidly at higher altitudes than at lower altitudes. Temperature also increased more rapidly at higher altitudes than at lower altitudes as solar activity increased. In addition, the height of the F2 peak was seen to increase with increasing solar activity, along with the oxygen ion fraction measured above the F2 peak. These results confirm that the topside ionosphere rises and expands with increasing solar activity.     

Changes in total electron content (TEC) during the annular solar eclipse of 15 January 2010

The solar eclipse of 15 January 2010 was an annular eclipse of the Sun with a maximum magnitude of 0.96 at 1.62°N, 69.29°E. To study the effect of this solar eclipse on the ionosphere the GPS data recorded at three different Indian stations Varanasi (Geographic latitude 25°, 16′N, longitude 82°, 59′E), Hyderabad (Geographic latitude 17°, 20′N, longitude 78°, 30′E) and Bengaluru (Geographic latitude 12°, 58′N, longitude 77°, 33′E) have been used to retrieve ionospheric total electron content (TEC). The ionospheric response to this rare event has been studied in terms of GPS-derived TEC observed at all the three Indian stations. A significant reduction in TEC reflected by all PRNs at all the three stations has been observed. The magnitude of the reduction in VTEC compared to quiet mean VTEC depends on latitude as well as longitude. The amount of reduction observed from different satellites (PRN) is different and depends on the location of the satellite from the solar eclipse path.     

GPS电离层反演方法研究及其在地震方面的应用

利用地基GPS数据计算了电离层单球壳模型穿刺点上的垂向总电子含量(VTEC), 根据VTEC和卫星及测站接收机的差分码间偏差DCB的不同时变特性采用了复弧法, 将VTEC作为局部变量(每30 min 一组, 可调), DCB作为一天的全局量进行解算. 在解算的过程中, 充分考虑VTEC的空间分布特性, 利用变异函数通过Kriging插值法建立电离层VTEC的二维格网模型, 并给出了卫星和接收机的差分码间偏差DCB. 通过与IGS结果的比较, 发现其结果可靠, 且时空分辨率和稳定性都有较大提高. 同时, 基于简化的三元样条插值基函数对电离层电子密度进行三维展开, 利用乘型代数重建技术MART算法构建了同批数据的四维层析成像结果, 获得了电离层电子密 度的四维分布. 其结果与CHAMP无线电掩星结果非常一致. 利用上面两种算法又分别对2008年5月长江三角洲地区地基GPS数据进行处理, 简要分析了该时段该地区上空电离层总电子含量和电子密度的变化情况及其对汶川地震的响应.      

Planetary wave coupling of the atmosphere–ionosphere system during the Northern winter of 2008/2009

This paper presents the global spatial (latitude and altitude) structure and temporal variability of the ∼23-day ionospheric zonally symmetric (s = 0) planetary wave (PW) seen in the Northern winter of 2008/2009 (October 2008–March 2009). It is shown that these ∼23-day ionospheric oscillations are forced from PWs propagating from below. The COSMIC ionospheric parameters f_oF2 and h_mF2 and electron density at fixed altitudes and the SABER temperatures were utilized in order to define the waves which are present simultaneously in the atmosphere and ionosphere. The long-period PWs from the two data sets have been extracted through the same data analysis method. The similarity between the lower thermospheric ∼23-day (s = 0) temperature PW and its ionospheric electron density response provides valuable and strong experimental evidence for confirming the paradigm of atmosphere–ionosphere coupling.     

World grid of cosmic ray vertical cut-off rigidity for the last decade

Cosmic ray cut-off rigidity tables and maps over the world concerning the epochs 2010, 2015 and the current one 2020 have been constructed. These maps display the effective cut-off rigidity in every five degrees in latitude and in longitude at the altitude of 20 km above the surface of the international reference ellipsoid. The values of the geomagnetic cut-off rigidity were calculated in every 5° in latitude and in every 15° in longitude applying the well-known method of particle trajectory calculations resulted from the theory of the particle motion in the Earth's magnetic field. The applied software employed the 12th Generation of the International Geomagnetic Reference Field (IGRF 12) and trajectories were calculated at 0.01 GV intervals in order to determine the vertical cut-off rigidity for each location. Beyond the use of the calculated cut-off rigidity values as a basic reference of charged particle access to different geographical locations during quiet and/or more intense geomagnetic periods, these results can be used for a long- term forecasting of the geomagnetic conditions variations.     

不同上边界条件下的极区电离层数值模拟

利用一维自洽的极区电离层模型,研究了沿磁力线方向不同电离层-磁层耦合条件下极区电离层的响应.此模型在110-610km的电离层空间区域内,综合求解描述极区电离层的连续性方程、动量方程和能量方程,以得到电离层数值解.研究发现,上边界条件在200 km以上的高度能显著地影响电离层参量的形态.较高的O+上行速度对应较低的F层峰值和较高的电子温度.不同边界O+上行速度对应的温度高度剖面完全不同.200km以上电子温度高度剖面不但由来自磁层的热流通量所控制,同时还受到场向O+速度的影响.对利用电离层模型研究电离层内部物理过程提出了建议.      

Dependence of Venus ionopause altitude and ionospheric magnetic field on solar wind dynamic pressure

The shape of the dayside Venus ionopause, and its dependence on solar wind parameters, is examined using Pioneer Venus Orbiter field and particle data. The ionopause is defined here as the altitude of pressure equality between magnetosheath magnetic pressure and ionospheric thermal pressure; its typical altitudes range from ～300 km near the subsolar point to ～900 km near the terminator. A strong correlation between ionopause altitude and magnetosheath magnetic pressure is demonstrated; correlation between magnetic pressure and the normally incident component of solar wind dynamic pressure is also evident. The data support the hypothesis of control of the ionopause altitude by solar wind dynamic pressure, manifested in the sheath as magnetic pressure. The presence of large scale magnetic fields in the ionosphere is observed primarily when dynamic pressure is high and the ionopause is low.     

The spatio-spectral localization approach to modeling VTEC over the western part of the USA using GPS observations

The precise ionosphere modeling is crucial and remains a challenge for GPS positioning and navigation, as well as many other Earth observation systems. In this research, two approaches have been proposed to model the vertical total electron content (VTEC) of the ionosphere using spherical slepian function. For the two-dimensional case, VTEC has been modeled in a sun-fixed reference frame and the three-dimensional approach based on the system of the three-dimensional base functions has been defined as the tensor product of the spherical slepian function for the longitude and latitude in an Earth-fixed reference frame, and the polynomial B-spline function for time. Rather than the spherical harmonics, the spherical slepian functions can be employed to produce the locally and globally orthogonal bases to optimally represent the data in any arbitrary region up to a given degree. The spherical slepian functions have been applied to the real data obtained from the ground-based GPS observation across the western part of the USA.     

Ion precipitation into the ionosphere during geomagnetic storms

This review presents numerous recent examples of interesting variations in the composition and intensity of the hot ion flux (10 eV - 15 keV/e) provided by the AUREOL-3 satellite as a function of latitude and local time during periods of magnetic activity. In particular, these results reveal that although H⁺ is the most abundant ion during magnetically quiet periods, the ion composition of hot plasma at ionospheric altitudes is quite variable, and depends strongly on magnetic activity; results obtained during main and recovery phases of several magnetic storms demonstrate clearly (below 15 keV/Q) the great importance of the low altitude ionospheric source (H⁺, O⁺, and to a lesser degree He⁺) particularly at low latitudes (L ～ 3 - 4) where the flux of O⁺ ions becomes very large and even dominates. The results of the AUREOL-3 ion spectrometers establish the fact that upflowing suprathermal ionospheric ions (E_i < 100 eV/e) appear over large regions of the auroral ionosphere, the polar caps, and the polar cusp, as well as in or at the boundary of the plasmasphere during magnetospheric substorms or magnetic storms, and may consequently contribute significantly to the plasma sheet and to the inner storm time ring current. Most of the properties of the storm time ring current found by the GEOS, SCATHA, and ISEE satellites apply to lower altitudes, although the role of the ionospheric and/or plasmaspheric source appears accentuated.     

Dayside ionospheric uplift during strong geomagnetic storms as detected by the CHAMP,SAC-C,TOPEX and Jason-1 satellites

Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.     

Negative VTEC with spherical harmonic expansion model under different solar activity conditions

The Earth’s ionosphere can be described by a spherical harmonic (SH) expansion up to a specific degree. However, there exist negative vertical total electron content (VTEC) values in the global ionosphere map (GIM) with the SH expansion model. In this contribution, we specifically investigated the negative VTEC values that are induced by the SH expansion model and validated the performance of the inequality-constrained least squares (ICLS) method in eliminating the negative VTEC values. The GPS data from 2004 to 2017 was selected to cover one solar cycle and the experiments under different solar activity conditions were analyzed. The results in our work show that the occurrence of the negative VTEC values is attributed to the deficiency of the SH expansion model when the VTEC itself is small instead of the unevenly distribution of the GNSS stations. The negative VTEC values appear periodically in the temporal domain, showing apparently one year and half year periods. During one year, two peaks in June and December can be observed in the time series of the negative VTEC values. The number of negative VTEC values in June is obvious larger than that in December. During one solar cycle, the number of negative VTEC values under quiet solar activity condition is obvious larger than that under strong solar activity condition. In the spatial domain, the appearance of the negative VTEC values is strongly related with the movement of the subsolar point. In the latitude of the subsolar point has the largest magnitude, the negative values will appear on the opposite hemisphere and the further from the subsolar point the more negative values. The maximum number of the negative VTEC values in the southern hemisphere appears in June, while the peak value in the northern hemisphere appears in December. The maximum number of negative VTEC values in the southern hemisphere is generally larger than that in the northern hemisphere. In addition, the negative VTEC values are distributed both at middle latitude and high latitude in the southern hemisphere, while they are mainly distributed at high latitude in the northern hemisphere. When the ICLS method is used, the negative VTEC values can be eliminated efficiently and it has nearly no influence on the positive VTEC values. The ICLS method can also improve the receiver’s differential code bias (DCB) and significantly decrease the unreasonable negative slant TEC (STEC) values along the lines of sight. Using the final GIM product of the Jet Propulsion Laboratory (JPLG) as a reference, the root mean square (RMS) of the ICLS solution shows maximum 25%, 20% and 45% improvement relative to the least squares (LS) solution at northern high latitude, southern middle latitude and southern high latitude, respectively.     

Latitudinal features of Total Electron Content over the African and European longitude sector following the St. Patrick’s day storm of 2015

GNSS TEC values have been obtained from 18 stations distributed from the magnetic equator to nearly 80°N magnetic dip in the African and west-European longitude sector corresponding to the March 17–18, 2015 geomagnetic storm. Significantly depleted ionosphere have been observed at stations north of 50°N geographic on March 18, 2015 following the above storm over a longitude swath 11.9°–21°E covering the Eastern Africa and Western European longitude sector. High ROTI values were noted on March 17^th at locations around 80°N magnetic dip. Two prominent peaks in PCN were noted around 09:00 UT and 14:00 UT on March 17, 2015 and around 15:00 UT on March 18, 2015. Daytime thermospheric (O/N₂) ratio was markedly less on March 18th at latitudes above 60°N geographic which is suggested to be the major driver behind depleted high latitude ionosphere during the recovery phase of the storm on March 18, 2015.     

Ionospheric reflection of the magnetic activity described by the index η

Differences in the external part of the vertical geomagnetic component point to the existence of local inhomogeneities in the magnetosphere or the ionosphere. Usually used magnetic indices are not sufficient to express the state of ionosphere, the common used global Kp index derived in the three-hour interval does not indicate much more rapidly changes appearing in ionosphere. Magnetic index η reflects ionospheric disturbances when other indices show very quiet conditions. Data of ionospheric characteristics (foE, foEs, h’E, h’F2) during 28-day long quiet day conditions (Kp = 0–2) in 2004 were analyzed. The correlations between strong local disturbances in ionosphere during very quiet days and high values of magnetic index η were found. The most sensitive to magnetic influence – ionospheric E layer data (foE characteristic) – reaches median deviations up to (+0.8 MHz and −0.8 MHz) during very low magnetic activity (Kp = 0–1). The high peaks (2–2.7) of the magnetic index η correlate in time with large local median deviations of foE. Such local deviations can suggest local inhomogeneities (vertical drifts) in the ionosphere. The correlation in space is not trivial. The strong peak of η is situated between the positive and negative deviations of foE. Additional observation is connected with correlation in time of the high η value with the negative median deviations of h’F2 (in some cases up to −90 km). The analysis was based on one-minute data recorded at each of 20 European Magnetic Observatories working in the INTERMAGNET network and from 19 ionosondes for 2004. Ionospheric data are sparse in time and in space in opposite to the magnetic data. The map of the magnetic indices can suggest the behavior of ionospheric characteristics in the areas where we have no data.     

Studies of electron and ion temperatures at 500 km altitude during sunrise using Indian SROSS C2 satellite

Solar dependence of electron and ion temperatures (T_e and T_i) in the ionosphere is studied using RPA data onboard SROSS C2 at an altitude of ∼500 km and 77°E longitude during early morning hours (04:00–07:00 LT) for three solar activities: solar minimum, moderate and maximum during winter, summer and equinox months in 10°S–20°N geomagnetic latitude. In winter the morning overshoot phenomenon is observed around 06:00 LT (T_e enhances to ∼4000 K) during low-solar activity and to T_e ∼ 3800 K, during higher solar activity. In summer, it is observed around 05:30 LT, but the rate of T_e enhancement is higher during moderate solar activity (∼2700 K/hr) than the low-solar activity (∼1700 K/hr). During equinox, this phenomenon is delayed and is observed around 06:00 LT (∼4200 K) during all three activities.