Reconstruction of topside density profile by using the topside sounder model profiler and digisonde data

To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (H_s), O⁺–H⁺ transition height (H_T), and their ratio Rt = H_s/H_T, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.     

DEMETER and DMSP satellite observations of the disturbed H/O ratio caused by Earth’s seismic activity in the Sumatra area during December 2004

In the present paper, plasma probe data taken from DEMETER and DMSP-F15 satellites were used to study the ion density and temperature disturbances in the morning topside ionosphere, caused by seismic activity at low latitudes. French DEMETER (Detection of Electro-Magnetic Emissions Transmitted from Earthquake Regions) micro-satellite mission had been especially designed to provide global scale observations in the topside ionosphere over seismically active regions. Onboard the DEMETER satellite, the thermal plasma instrument called “Instrument Analyser de Plasma” (IAP) provides ion mass and densities, ion temperature, three component ion drift and ion density irregularities measurements. As a part of “Defense Meteorological Satellite Program”, DMSP-F15 satellite is on orbit operation since 1999. It provides ionospheric plasma diagnostics by means of the “Special Sensor-Ion, Electron and Scintillations” (SSIES-2) instrument. We examined few examples of possible seismic effects in the equatorial ionosphere, probably associated with seismic activity during December month in the area of Sumatra Island, including main shock of giant Sumatra event. It is found that the localized topside ionospheric disturbances appear close to the epicenters of certain earthquakes in the Sumatra region. In two cases, ion H⁺/O⁺ ratio rises more than one hour before the main shock, due to the O⁺ density decrease at the winter side of the geomagnetic equator, with longitudinally closest location to the epicenter of the earthquakes. These anomalous depletions in O⁺ density do exist in all cases of SSIES-2 data. Particularly for Sumatra main event, more than one hour after the main shock, we observe large-scale depletion in O⁺ density northward of the geomagnetic equator at winter side hemisphere. Associated with O⁺ depletion, ion temperature latitudinal profile around the geomagnetic equator shows enhanced asymmetry with minimum at the summer side and maximum in positive Ti deviation from mean value at the winter side. This disturbance lasted for more than three hours, later in time observed at the same place by IAP/DEMETER.     

A topside investigation over a mid-latitude digisonde station in Cyprus

We examine the systematic differences between topside electron density measurements recorded by different techniques over the low-middle latitude operating European station in Nicosia, Cyprus (geographical coordinates: 35.14^oN, 33.2^oE), (magnetic coordinates 31.86^oN, 111.83 ^oE). These techniques include space-based in-situ data by Langmuir probes on board.European Space Agency (ESA) Swarm satellites, radio occultation measurements on board low Earth orbit (LEO) satellites from the COSMIC/FORMOSAT-3 mission and ground-based extrapolated topside electron density profiles from manually scaled ionograms. The measurements are also compared with International Reference Ionosphere Model (IRI-2016) topside estimations and IRI-corrected NeQuick topside formulation (method proposed by Pezzopane and Pignalberi (2019)). The comparison of Swarm and COSMIC observations with digisonde and IRI estimations verifies that in the majority of cases digisonde underestimates while IRI overestimates Swarm observations but in general, IRI provides a better topside representation than the digisonde. For COSMIC and digisonde profiles matched at the F layer peak the digisonde systematically underestimates topside COSMIC electron density values and the relative difference between COSMIC and digisonde increases with altitude (above hmF2), while IRI overestimates the topside COSMIC electron density but after a certain altitude (~150 km above hmF2) this overestimation starts to decrease with altitude. The IRI-corrected NeQuick underestimates the majority of topside COSMIC electron density profiles and relative difference is lower up to approximately 100 km (above the hmF2) and then it increases. The overall performance of IRI-corrected NeQuick improves with respect to IRI and digisonde.     

Positive ion composition and electron density in a combined auroral and NLC event

The positive ion composition and electron density were measured in the lower ionosphere above Kiruna in salvo A of CAMP (Cold Arctic Mesopause Project). The CAMP/P (S37/P) payload carrying a magnetic ion spectrometer, positive ion and electron probes, and propagation experiments was launched on 3 August 1982 2332 UT during extended Noctilucent Clouds (NLC) and auroral activities over Kiruna. The measured electron density was 5×10³cm^?3 at 80 km and 2.5×10⁵cm^?3 at 90 km. The increase of ion and electron densities in the D- and E-region during twilight was caused by precipitating auroral particles. The height distribution of the positive ions measured by the mass spectrometer in the mass range 19–280 amu is different from a winter flight with similar auroral conditions. Below 85.5 km proton hydrates H⁺(H₂O)₃ ? H⁺(H₂O)₈ were the dominant ions. The heaviest proton hydrates H⁺(H₂O)₇ and H⁺(H₂O)₈ were most abundant at 82–85.5 km, the altitude of visible NLC. Above 85.5 km O₂⁺ and NO⁺ became dominant. A small metal ion layer was observed between 90.5–93 km with a maximum ion density of 10% of the total positive ion density at 91 km altitude. The metal ion density disappeared within about a km below 90.5 km.     

Comparison of topside ionosphere scale height modeled by the Topside Sounder Model and incoherent scatter radar ionospheric model

The topside ionosphere scale height extracted from two empirical models are compared in the paper. The Topside Sounder Model (TSM) provides directly the scale height (H_T), while the incoherent scatter radar ionospheric model (ISRIM) provides electron density profiles and its scale height (H_R) is determined by the lowest gradient in the topside part of the profile. H_T and H_R are presented for 7 ISR locations along with their dependences on season, local time, solar flux F10.7, and geomagnetic index ap. Comparison reveals that H_T values are systematically lower than respective H_R values as the average offset for all 7 stations is 55 km. For the midlatitude stations Arecibo, Shigaraki, and Millstone Hill this difference is reduced to 43 km. The range of variations of H_R is much larger than that of H_T, as the H_T range overlaps the lower part of the H_R range. Dependences on ap, DoY and LT are much stronger in the ISRIM than in TSM. This results in much larger values of H_R at higher ap. Diurnal amplitude of H_R is much larger than that of H_T, with large maximum of H_R at night. The present comparison yields the conclusion that the ISR measurements provide steeper topside Ne profiles than that provided by the topside sounders.     

Field-aligned plasmaspheric flows at moderate latitudes

In a model of the plasmasphere, coupled time-dependent continuity, momentum and energy equations are solved for thermal O⁺, H⁺ and electrons. The field-aligned mass flow coupling and thermal coupling of conjugate ionsopheres via the protonosphere are studied. For solstice conditions, thermal coupling between conjugate hemispheres gives rise to very strong upward flows of O⁺ in the topside ionosphere of the summer hemisphere at the time of sunrise in the conjugate (winter) ionosphere; a less marked effect (but with downward flow) occurs in the summer ionosphere at winter sunset. In addition, there are strong upward and downward flows of O⁺ at local sunrise and sunset, respectively, in both hemispheres. At both L = 1.5 and L = 3, the 24-hour time-integrated interhemispheric H⁺ flux is in the summer - winter direction. At L = 1.5 its magnitude is in good agreement with the magnitude of the time- integrated field-aligned plasma (O⁺ + H⁺) flux at 1000 km altitude; there is no such agreement at L = 3.     

Positive ion distributions in the morning auroral zone: Local acceleration and drift effects

We report on the typical structure of the large scale ion precipitation in the morning sector of the auroral zone and associated low frequency electromagnetic waves. Data obtained during near radial passes of the AUREOL-3 satellite point to a distinction between two main precipitation regions: 1) In the poleward part of the auroral zone the latitudinal variation of the average energy (or temperature) of the precipitated ions (mainly H⁺) indicate that they are adiabatically accelerated in the outer magnetosphere. This “high energy” (? 3 to > 20 keV) precipitation is usually associated with a low energy (E < 110 eV) upward flowing 0⁺ and H⁺ component, and 2) near the boundary between discrete and diffuse electron aurorae a drastic change in the ion characteristics is observed. The flux of energetic precipitated H⁺ ions is sharply reduced, which suggests the formation of an Alfvén layer. However, intense fluxes of precipitated H⁺, O⁺, and He⁺ ions with energies < 3 keV are observed equatorward of the Alfvén layer, in coincidence with the diffuse aurora and in association with quasi-monochromatic electromagnetic waves with frequencies around the proton gyrofrequency. As the characteristic convection and bounce times of the low energy upward flowing ion component are comparable (τ > 3 hours) we suggest that the precipitation of ionospheric ions inside the diffuse aurora results from convection and corotation of the ions accelerated to suprathermal energies at higher latitudes.     

Mesospheric H2O and H2O2 densities inferred from in situ positive ion composition measurement

The measurements of positive ion composition in the high latitude D-region have revealed an excess of 34⁺ under distrubed conditions which has been interpreted as H₂O₂⁺. At the same altitude range near the transition height oxonium ions were measured as well. This paper presents a new model for the production and loss of oxonium ions with their production from H₂O₂⁺ + H₂O → H₃O⁺ + HO₂ and their loss by attachment of N₂ and/or CO₂. A reaction constant of 8.5×10^?28 (300/T)⁴ cm⁶s^?1 has been obtained for the three body attachment H₃O⁺ + CO₂ + M → H₃O⁺.CO₂ + M from the measured density profile of 63⁺ in flight 18.1020. Mesospheric H₂O and H₂O₂ densities are inferred from measurements of four high latitude ion compositions based on the oxonium model. The mixing ratios of hydrogen peroxide are up to two orders of magnitude higher compared to previous model calculations. In order to explain the missing production of odd hydrogen, we consider larger O(¹D) densities, surface reactions of O(³P) on particles, and cathalytic photodissociation of water vapor on aerosol particles.     

Comparing topside and bottomside-measured characteristics of the F2 layer peak

The ionospheric characteristics of the F2 layer peak have been measured with ionosondes from the ground or with satellites from space. The most common characteristics are the F2-peak density NmF2 and peak height hmF2. In addition to these two parameters this paper studies the F2-peak scale height. Comparing the median values of hmF2 and NmF2 obtained from topside and bottomside sounding shows good agreement in general. The Chapman scale height values for the F2 layer peak derived from topside profiles, H_m,top, are generally several times larger than H_m,bot derived from bottomside profiles.     

Modelling bottom and topside electron density and TEC with profile data from topside ionograms

The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.     

Relative contributions of terrestrial and solar wind ions in the plasma sheet

A major uncertainty concerning the origins of plasma sheet ions is due to the fact that terrestrial H⁺ can have similar fluxes and energies as H⁺ from the solar wind. The situation is especially ambiguous during magnetically quiet conditions (AE < 60γ) when H⁺ typically contributes more than 90% of the plasma sheet ion population. In this study we examine that problem using a large data set obtained by the ISEE-1 Plasma Composition Experiment. The data suggest that one component of the H⁺ increases in energy with increasing activity, roughly in proportion to $\text{1}\text{4}$ the energy of the He⁺⁺, whereas the other H⁺ component has about the same energy at all activity levels, as do the O⁺ and the He⁺. If we can assume that the H⁺ of solar wind origin on the average has about the same energy-per-nucleon as the He⁺⁺, which is presumably almost entirely from the solar wind, then the data imply that as much as 20–30% of the H⁺ can be of terrestrial origin even during quiet conditions.     

Storm time behavior of topside scale height inferred from the ionosphere–plasmasphere model driven by the F2 layer peak and GPS-TEC observations

The International Reference Ionosphere model extended to the plasmasphere, IRI-Plas, presents global electron density profiles and total electron content, TECiri, up to the altitude of the GPS satellites (20,000 km). The model code is modified by input of GPS-derived total electron content, TECgps, so that the topside scale height, Hsc, is obtained minimizing in one step the difference between TECiri and TECgps observation. The topside basis scale height, Hsc, presents the distance in km above the peak height at which the peak plasma density, NmF2, decays by a factor of e (∼2.718). The ionosonde derived F2 layer peak density and height and GPS-derived TECgps data are used with IRI-Plas code during the main phase of more than 100 space weather storms for a period of 1999–2006. Data of seven stations are used for the analysis, and data from five other stations served as testing database. It is found that the topside basis scale height is growing (depressing) when the peak electron density (critical frequency foF2) and electron content are decreased (increased) compared to the median value, and vice versa. Relative variability of the scale height, rHsc, and the instantaneous Hsc are inferred analytically in a function of the instantaneous foF2, median fmF2 and median Hmsc avoiding a reference to geomagnetic indices. Results of validation suggest reliability of proposed algorithm for implementation in an operational mode.     

F-Region ion composition modeling

Measurements of the principal ion species of the F₁- and F₂- regions have been used to develop an empirical model of the ion composition for altitudes between 150 and 500 km. The species measured by the S3-1 satellite include N⁺, O⁺, N₂⁺, NO⁺ and O₂⁺. The data were obtained near the minimum of the solar cycle, thus limited information on the ionospheric variation with solar flux is available. However, the range of latitude, altitude, local time and geomagnetic activity does provide a useful basis for modeling the F-region. The ion composition measurements have been used to provide a model for relative ion composition which is compatible with the total ion density from the International Reference Ionosphere model.     

Nightglow studies with the world's largest optical telescope

By making use of 0.02-nm-resolution sky spectra from the HIRES echelle spectrograph at the W.M. Keck(I) observatory on Mauna Kea, HI, obtained during normal astronomical observations, we have shown that 650–870 nm emission from the vibrationally-excited levels of the O₂(b¹Σ⁺_g) state is a significant component of the terrestrial nightglow, with a total average intensity from levels v = 1–15 of approximately 150 R. The b¹Σ⁺_g state vibrational distribution is bimodal, with peaks at v = 3, 4 and v = 12, and a deep minimum at v = 8 containing only ∼5% of the v = 3 population. The b-X 0-0 band (the Franhofer A-band) is discernible via isotopic emission, with an intensity comparable to that of the strongest of the new b-X bands. There are indications that the vibrationally-excited b¹Σ⁺_g state emissions correlate temporally with the OH Meinel band emission, and satellite measurements appear to show that the relatively strong emission from the O₂(b-X) 4-3 band originates at an alitutde near 87 km, i.e. in the OH region.In an 8-minute observation on 20 Nov 1999, the Keck telescope was pointed at Venus, and the first high-resolution nightglow spectrum was recorded in the visible spectral region. The oxygen green line was revealed, with an intensity near 150 R, approximately the same as in the terrestrial atmosphere. The only previous observation, from the Venera 9/10 orbiters, did not detect this emission.     

Temporal and spatial variations observed in the ionospheric composition of Venus: Implications for empirical modelling

In situ measurements of the thermal ion composition of the ionosphere of Venus have been obtained for a period of two Venus years from the Bennett rf ion mass spectrometer on the Pioneer Venus Orbiter. Ion measurements within an altitude interval of 160 to 300 kilometers, corresponding to an overall latitude interval of about ?4° to 34°N, are assembled from the interval December 1978 to March 1980. This time interval corresponds to two revolutions of Venus about the Sun, designated as two “diurnal cycles”. The distributions of several ion species in this data base have been sorted to identify temporal and spatial variations, and to determine the feasibility of an analytical representation of the experimental results. The first results from the sorting of several prominent ions including O⁺, O₂⁺, and H⁺ and several minor ions including CO₂⁺, C⁺, and H₂⁺ reveal significant diurnal variations, with superimposed modulation associated with solar activity and solar wind variations. The diurnal variation consists of strong day to night contrast in the ion concentrations, with differences of one to two orders of magnitude, depending upon ion mass and altitude. The concentrations of O₂⁺, O⁺, CO₂⁺ and C⁺ peak throughout the dayside decreasing sharply at the terminators to nightside levels, lower by one to two orders of magnitude relative to the dayside. The diurnal variations of the light ions H⁺ and H₂⁺ peak during the night, exhibiting asymmetric nightside bulges favoring the pre-dawn sector, near 0400 solar hour angle. Superimposed upon the diurnal distributions are modulation signatures which correlate well with modulation in the F_10.7 index, indicating a strong influence of solar variability on the ion production and distribution. The influence of solar wind perturbations upon the ion distributions are also indicated, by a significant increase in the scatter of the observations with increasing altitude as higher altitudes, approaching 300 kilometers, are sampled. Together, these temporal and spatial variations make the task of modelling the ionosphere of Venus both very interesting and challenging.     

Towards a new reference model of hmF2 for IRI

A numerical model of the peak height of the F2 layer, hmF2_top, is derived from the topside sounding database of 90,000 electron density profiles for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The model of regular hmF2 variations is produced in terms of local time, season, geomagnetic latitude, geodetic longitude and solar radio flux. No geomagnetic activity trends were discernible in the topside sounding data. The nighttime peak of hmF2_top evident for mid-latitudes disappears near the geomagnetic equator where a maximum of hmF2_top occurs at sunset hours when it can exceed 500 km at solar maximum. The hmF2 given by the IRI exceeds hmF2_top at the low solar activities. The hmF2_top, obtained by extrapolation of the first derivative of the topside profile to zero shows saturation similar to foF2 the greater the solar activity. The proposed model differs from hmF2 given by IRI based on M(3000)F2 to hmF2 conversion by empirical relationships in terms of foF2, foE and R12 with these quantities mapped globally by the ITU-R (former CCIR) from ground-based ionosonde data. The differences can be attributed to the different techniques of the peak height derivation, different epochs and different global distribution of the source data as well as the different mathematical functions involved in the maps and the model presentation.     

Ion composition behaviour in low and mid-latitudes during high solar activity

In the present analysis, the mass spectrometer data from the ISS-b satellite, available in the form of contour plots at an average height of 1100 km for every alternate hour, is used. This analysis showed some interesting results in terms of the diurnal variation of the transition heights; at times dominance of He⁺ over the other ions, and the seasonal variations of different constituents. In the development of ion composition models, it is suggested that this type of result should be taken into consideration.     

Modeling of sub-ionospheric VLF signal perturbations associated with total solar eclipse, 2009 in Indian subcontinent

During the total solar eclipse of 2009, a week-long campaign was conducted in the Indian sub-continent to study the low-latitude D-region ionosphere using the very low frequency (VLF) signal from the Indian Navy transmitter (call sign: VTX3) operating at 18.2 kHz. It was observed that in several places, the signal amplitude is enhanced while in other places the amplitude is reduced. We simulated the observational results using the well known Long Wavelength Propagation Capability (LWPC) code. As a first order approximation, the ionospheric parameters were assumed to vary according to the degree of solar obscuration on the way to the receivers. This automatically brought in non-uniformity of the ionospheric parameters along the propagation paths. We find that an assumption of 4 km increase of lower ionospheric height for places going through totality in the propagation path simulate the observations very well at Kathmandu and Raiganj. We find an increase of the height parameter by h^′=+3.0$h^{\prime }=+3.0$ km for the VTX-Malda path and h^′=+1.8$h^{\prime }=+1.8$ km for the VTX-Kolkata path. We also present, as an example, the altitude variation of electron number density throughout the eclipse time at Raiganj.     

Analysis of electron content variations over Japan during solar minimum: Observations and modeling

Comparative analysis of GPS TEC data and FORMOSAT-3/COSMIC radio occultation measurements was carried out for Japan region during period of the extremely prolonged solar minimum of cycle 23/24. COSMIC data for different seasons corresponded to equinox and solstices of the years 2007–2009 were analyzed. All selected electron density profiles were integrated up to the height of 700 km (altitude of COSMIC satellites), the monthly median estimates of Ionospheric Electron Content (IEC) were retrieved with use of spherical harmonics expansion. Monthly medians of TEC values were calculated from diurnal variations of GPS TEC estimates during considered month. Joint analysis of GPS TEC and COSMIC data allows us to extract and estimate electron content corresponded to the ionosphere (its bottom and topside parts) and the plasmasphere (h > 700 km) for different seasons of 2007–2009. Percentage contribution of ECpl to GPS TEC indicates the clear dependence from the time and varies from a minimum of about 25–50% during day-time to the value of 50–75% at night-time. Contribution of both bottom-side and topside IEC has minimal values during winter season in compare with summer season (for both day- and night-time). On average bottom-side IEC contributes about 5–10% of GPS TEC during night and about 20–27% during day-time. Topside IEC contributes about 15–20% of GPS TEC during night and about 35–40% during day-time. The obtained results were compared with TEC, IEC and ECpl estimates retrieved by Standard Plasmasphere–Ionosphere Model that has the plasmasphere extension up to 20,000 km (GPS orbit).     

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.