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.     

Modelling studies of ionospheric convection in northern and southern polar regions

The realistic model of Quegan et al. has been used to investigate the convection paths of ionospheric plasma at 300 km altitude, for different polar cap radii and in both hemispheres. Taking the Northern magnetic dip pole to be at a co-latitude of 11° and the Southern magnetic dip pole at a co-latitude of 23°, these paths are presented in a Sun-Earth frame, with the position of the Earth's axis fixed as it is on 21 March, as polar plots centred on the magnetic pole. There are marked hemispheric differences between 13 and 23 L.T., particularly near the stagnation region at 18 to 21 L.T., but only minor differences between 00 and 12 L.T., when the radius of the polar cap exceeds 12°. For a smaller polar cap, the differences between the hemispheres are small at all local times. The time taken to perform a complete circuit is most dependent on the polar cap radius, and most variable - between 15 and 36 h - for convection paths starting near 60° latitude. The time that plasma convecting from noon to near midnight across the Northern polar cap spends within the 10° co-latitude circle increases from 6 h, for a polar cap radius of 10°, to 11.5 h at 17°. These results are compared and contrasted with other model calculation results and with some ground-based and satellite observations of plasma densities at high latitudes.     

The D-region ion composition and electron density response to strong solar proton events (model simulations)

Fluxes of energetic solar protons penetrate deep into the Earth’s polar cap middle atmosphere. Interacting with molecules of the air they cause additional dissociation and ionization, and the formed NO_x, OH_y and ions enter chemical and ion-molecular reactions. Induced changes of the ionospheric D-layer are modeled by a 1D model of lower ionosphere with chemistry, using neutral species concentrations calculated by a 1D photochemical time-dependent model. Changes of the electron and ion densities, and the most important ionospheric parameters are calculated after SPE with the onset on July 14, 2000 and the results are compared with our results obtained previously for the October 19, 1989 SPE. It is shown that not only electron density increases after SPE, but also the amount of clusters. It is found that the magnitude of the ionospheric response depends on season.     

Relationship between electron density and electron temperature as a function of solar activity

This paper attempts to examine the control of electron density and solar activity on the F-region electron temperature. This is achieved by obtaining coefficients relating electron temperature with electron density and solar activity by using incoherent scatter radar measurements at Arecibo for the period August 1966 to May 1977. These coefficients are then used to construct an empirical model of F-region electron temperature. The model values are compared with measurements at other locations and show reasonable agreement.     

Theoretical model of electron temperature and electron density: Comparisons with IRI and MSIS

An earlier theoretical model (UW-87) accurately predicted the electron temperature in the daytime F-region but suggested N₂ concentrations significantly greater than the predictions of MSIS-86. This discrepancy is resolved when the model is developed to include the effects of vibrationally excited nitrogen molecules and electronically excited oxygen ions on the F-region recombination rate. The revised model (UW-92) continues to predict electron temperatures close to the layer peak with great accuracy but it is now more closely consistent with MSIS. However, the electron temperatures predicted by this model, which are in close agreement with EISCAT observations, are significantly higher than the values predicted by the international Reference Ionosphere.     

Electric field and electron density measurements in the equatorial E-region

A centaure rocket, with payloads of Langmuir probe and Electric field probe, was launched from Thumba (8° 31'N, O° 47'S dip), India on February 12, 1981 at 1057 Hrs IST. The aim of the experiment was to study the role of localised electric fields in the generation of plasma density irregularities through cross field instability and the two-stream instability mechanism. The rocket was launched at a time when Type I irregularities were observed with VHF radar at Thumba.     

On long-term changes of electron density in the upper mesosphere

Recent review study done jointly by 19 experts of 17 institutes shows zero trend of temperature in the upper mesosphere. In the light of this latest development, we have examined the long-term changes in electron density, [e], in this region. The study has been concentrated at 80 km. At this altitude, electrons are mainly produced by the interaction of nitric oxide, NO, by solar Ly-α. Any long-term change in this flux will affect trend of [e]. Considering this flux proportional to 10.7 cm solar flux, analysis of available 10.7 cm solar flux data from 1948 to 2003 has been made. A decreasing trend up to about 1970 and then an increasing trend are found. The over-all increasing trend of Ly-α flux during the past five decades is ∼0.17% per year. This increase also gives a ∼0.17% increasing trend per year in [e]. This non-anthropogenic increase is much less compared to the observed increase in [e] which is reported to be >0.7% per year. The observed increase in [e] of this magnitude will then, predominantly, be due to the anthropogenic effect. In zero trend in temperature, significant change in electron loss coefficient, α_eff, and [NO] are unlikely to take place to cause a significant change in [e]. The increase in [e] > 0.7% per year then can be explained by considering a decreasing trend in [O₂].     

A payload for the study of electric fields and electron density in the equatorial region

A rocket borne payload for simultaneous measurement of the electric field along and perpendicular to the rocket spin axis and the electron density in the medium was developed and flown from Thumba (8° 31′N, 0° 47′S dip) onboard two Centaure rockets for the study of plasma dynamcis in the equatorial E-region. The arrangement of sensors in this payload allows near continuous measurements of some of these parameters to be made.     

Seismo-ionospheric anomalies in total electron content of the GIM and electron density of DEMETER before the 27 February 2010 M8.8 Chile earthquake

In this paper we examine the pre-earthquake ionospheric anomalies by the total electron content (TEC) extracted from GIM (global ionospheric map) and the electron density (Ne) observed by the DEMETER (Detection of Electro-Magnetic Emissions Transmitted from Earthquake Regions) satellite during the 2010 M8.8 Chile earthquake. Temporal variations show the nighttime TEC and Ne simultaneously increase 9–19 days before the earthquake. A cross-comparison of data recorded during the period of 1 February to 3 March in 2006–2010 confirms the above temporal anomalies specifically appear in 2010. The spatial analyses show that the anomalies tend to appear over the epicenter.     

Local ionospheric electron density profile reconstruction in real time from simultaneous ground-based GNSS and ionosonde measurements

The purpose of the LIEDR (local ionospheric electron density profile reconstruction) system is to acquire and process data from simultaneous ground-based total electron content (TEC) and digital ionosonde measurements, and subsequently to deduce the vertical electron density distribution above the ionosonde’s location. LIEDR is primarily designed to operate in real time for service applications and, for research applications and further development of the system, in a post-processing mode. The system is suitable for use at sites where collocated TEC and digital ionosonde measurements are available. Developments, implementations, and some preliminary results are presented and discussed in view of possible applications.     

Regional 4-D modeling of the ionospheric electron density from satellite data and IRI

Accurate knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying ionosphere physics. During the last decade Global Navigation Satellite Systems (GNSS) have become a promising tool for monitoring ionospheric parameters such as the total electron content (TEC). In this contribution we present a four-dimensional (4-D) model of the electron density consisting of a given reference part, i.e., the International Reference Ionosphere (IRI), and an unknown correction term expanded in terms of multi-dimensional base functions. The corresponding series coefficients are calculable from the satellite measurements by applying parameter estimation procedures. Since satellite data are usually sampled between GPS satellites and ground stations, finer structures of the electron density are modelable just in regions with a sufficient number of ground stations. The proposed method is applied to simulated geometry-free GPS phase measurements. The procedure can be used, for example, to study the equatorial anomaly.     

Midday bottomside electron density profiles during moderate solar activity and comparison with IRI-2001

Bottomside electron density (Ne-h) profiles during midday (10–14 h) are analyzed using modern digital ionosonde observations at a low-middle latitude station, New Delhi (28.6N, 77.2E, dip 42.4N), for the period from January 2003 to December 2003, pertaining to moderate solar activity (MSA). Each individual profile is normalized with respect to the peak height and density (hmF2, NmF2) of the F2-region. These profiles are compared with those obtained from the International Reference Ionosphere (IRI-2001) model. Bilitza [Bilitza, D. International Reference Ionosphere 2000. Radio Sci. 36 (2), 261–275, 2001] using both the options namely: Gulyaeva’s model [Gulyaeva, T.L. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7 (6) 39–48, 1987] and B0 Tab. option [Bilitza, D., Radicella, S.M., Reinisch, B.W., Adeniyi, J.O., Mosert Gonzalez, M.E., Zhang, S.R., Obrou, O. New B0 and B1 models for IRI. Adv. Space Res. 25 (1), 89–95, 2000]. The study reveals that during summer and equinox, the IRI model with B0 Tab. option in general, produces better agreement with the observed median profiles, while the IRI predictions using Gulyaeva’s option, overestimate the electron density distribution at all the heights below the F2-peak. However, during winter, in general, the IRI model, using both the options, reveals shows fairly good agreement with the observations.     

An empirical model of electron density for low and middle latitudes below 200 km

Our empirical model of electron density (n_e) for quiet and weakly disturbed geomagnetic conditions (Kp not greater 4) takes account of comparative analysis of existing models and of experimental data obtained by rockets and incoherent scatter radar. The model describes the n_e distribution in the 80 to 200 km height range at low and middle latitudes, and to some extent, in the subauroral region. It is presented in analytical form thus allowing one to calculate electron density profiles for any time. The electron density distribution at 140 km depends on the season (day of the year) and on the solar zenith angle. Profile variations during the day are for one season shown. Different from other models, ours specifies the variations during sunrise and sunset and reflects the particular profile shape at night admitting the occurrence of an intermediate layer.     

Ionospheric electron density over Resolute Bay according to E-CHAIM model and RISR radar measurements

In this study, predictions of the E-CHAIM ionospheric model are compared with measurements by the incoherent scatter radars RISR at Resolute Bay, Canada, in the northern polar cap. Reasonable coverage was available for all seasons except winter for which no conclusions were drawn. It is shown that ratios of the model-to measured electron densities are close to unity in the central part of the F layer, around its peak. This is particularly evident for summer daytime. Distributions of the ratios are wider for other seasons indicating larger number of cases when the model underestimates or overestimates. E-CHAIM underestimates the electron density at ionospheric topside and bottomside by ～ 10–20 %. At the bottomside, the underestimations are strongest in summer and equinoctial nighttime. At the topside, the underestimations are strongest in autumn nighttime. Model overestimations are noticeable in the middle part of the F layer during dawn hours in autumn. Overall, the model tends to not predict highest-observed peak electron densities and the largest-observed heights of the peak.     

Comparison between d- and lower e-region electron density profiles and IRI-79

The D and E-region electron density profiles obtained by different techniques are compared with the IRI-79 model to see how they fit. The rocket data showed good agreement. However discrepancies between the observed and model values were found especially for solar zenith angle greater than 50 degrees.     

Reanalysis of radiation belt electron phase space density using various boundary conditions and loss models

Data assimilation is becoming an increasingly important tool for understanding the near Earth hazardous radiation environments. Reanalysis of the radiation belts can be used to identify the electron acceleration mechanism and distinguish local acceleration from radial diffusion. However, for any practical applications we need to determine how reliable is reanalysis, and how significant is the dependence of the results on the assumptions of the code and choice of boundary conditions. We present the sensitivity of reanalysis of the radiation belt electron phase space density (PSD) to the assumed location of the outer boundary, using the VERB code and a Kalman filter. We analyze the sensitivity of reanalysis to changes in the electron-loss throughout the domain, and the sensitivity to the assumed boundary condition and its effect on the innovation vector. All the simulations presented in this study for all assumed loss models and boundary conditions, show that peaks in the phase space density of relativistic electrons build up between 4.5 and 6 R_E during relativistic electron flux enhancements in the outer radiation belt. This clearly shows that peaks build up in the heart of the electron radiation belt independent of the assumptions in the model, and that local acceleration is operating there. The work here is also an important step toward performing reanalysis using observations from current and future missions.     

Linkage of the ionospheric peak electron density and height deduced from the topside sounding data

Topside sounding electron density profiles are analyzed to explore interrelations of the F2 layer critical frequency and the peak height for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The foF2 and hmF2 are delivered with exponential extrapolation of electron density profile to zero of its 1st derivative. It is shown that the linear regression exists between foF2 and hmF2 under different conditions. The linkage between the two parameters amended to the empirical model of the peak height [Gulyaeva, T.L., Bradley, P.A., Stanislawska, I., Juchnikowski, G. Towards a new reference model of hmF2 for IRI. Adv. Space Res. 42, 666–672, doi:10.1016/j.asr.2008.02.021, 2008] results in an empirical model of the both foF2 and hmF2 expressed by superposition of functions in terms of local-time, season, geodetic longitude, modified dip latitude and solar activity. For the solar activity we use a proxy Fsp index averaged from the mean solar radio flux F10.7s for the past 81 days (3 solar rotations) and F10.7 value for 1 day prior the day of observation. Impact of geomagnetic activity is not discernible with the topside sounding data due to mixed positive and negative storm-time effects. Appreciable differences have been revealed between IRI-CCIR predictions and outcome of the new model which might be attributed to the different techniques of the peak electron density and 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.