A different method to update monthly median hmF2 values

The height, h_mF2, and the electron density, N_mF2, of the F2 peak are key model parameters to characterize the actual state of the ionosphere. These parameters, or alternatively the propagation factor, M3000F2, and the critical frequency, f_oF2, of the F2 peak, which are related to h_mF2 and N_mF2, are used to anchor the electron density vertical profile computed with different models such as the International Reference Ionosphere ( Bilitza, 2002), as well as for radio propagation forecast purposes. Long time series of these parameters only exist in an inhomogeneous distribution of points over the surface of Earth, where dedicated instruments (typically ionosondes) have been working for many years. A commonly used procedure for representing median values of the aforementioned parameters all over the globe is the one recommended by the ITU-R ( ITU-R, 1997). This procedure, known as the Jones and Gallet mapping technique, was based on ionosondes measurements gathered from 1954 to 1958 by a global network of around 150 ionospheric stations (  and ). Even though several decades have passed since the development of that innovative work, only few efforts have been dedicated to establish a new mapping technique for computing h_mF2 and N_mF2 median values at global scale or to improve the old method using the increased observational database. Therefore, in this work three different procedures to describe the daily and global behavior of the height of the F2 peak are presented. All of them represent a different and simplified method to estimate h_mF2 and are based on different mathematical expressions. The advantages and disadvantages of these three techniques are analyzed, leading to the conclusion that the recommended procedure to represent h_mF2 is best characterized by a Spherical Harmonics expansion of degree and order equal to 15, since the differences between the h_mF2 values obtained with the Jones and Gallet technique and those obtained using the abovementioned procedure are of only 1%.     

The accuracy of data from ionosondes for the estimation of hmF2 and the identification of global change in the ionosphere

A long temporal series of simulated ionograms was generated with a superimposed secular variation of −14 km/century on the h_mF2 parameter. These ionograms were interpreted by the automatic scaling program Autoscala. By applying four different empirical formulas, four artificial series of h_mF2 were generated and then processed with the same methods used by other authors for real data sets. Data analysis of the simulated ionograms revealed the artificially imposed long-term trend. These results lead to the conclusion, that regardless of the empirical formula used, the accuracy of h_mF2 from ionosonde measurements would be adequate to observe a long-term trend of −14 km/century.     

Comparison of observed hmF2 and IRI 2007 model with M(3000)F2 estimation of hmF2 at low solar activity for an equatorial station

This is to investigate ways of improving the Equatorial F2-layer peak heights estimated from M(3000)F2 ionosonde data measured using the Ionospheric Prediction Service (IPS-42) sounder at Ouagadougou, Burkina Faso (Latitude +12.4°N, Longitude +1.5°W, Dip latitude +5.9°N) during a low solar activity year (1995). For this purpose, we have compared the observed h_mF2 (h_mF2_obs) deduced using an algorithm from scaled virtual heights of quiet day ionograms and the predicted h_mF2 values which is given by the IRI 2007 model (h_mF2_{IRI 2007}) with the ionosonde measured M(3000)F2 estimation of the h_mF2 values (h_mF2_est) respectively. The correlation coefficients R² for all the seasons were found to range from 0.259 to 0.692 for h_mF2_obs values, while it ranges from 0.551 to 0.875 for the h_mF2_{IRI 2007} values. During the nighttime, estimated h_mF2 (h_mF2_est) was found to be positively correlated with the h_mF2_obs values by the post-sunset peak representation which is also represented by the h_mF2_{IRI 2007} values. We also investigated the validity of the h_mF2_est values by finding the percentage deviations when compared with the h_mF2_obs and h_mF2_{IRI 2007}.     

Plasma chemical reactions in C2H2/N2, C2H4/N2, and C2H6/N2 gas mixtures of a laboratory dielectric barrier discharge

Plasma chemical reactions in C₂H₂/N₂, C₂H₄/N₂, and C₂H₆/N₂ gas mixtures have been studied by means of mass spectrometry at a medium pressure of 300 mbar in a laboratory dielectric barrier discharge. A major reaction scheme is production of larger hydrocarbons like C_nH_m with n up to 12 including formation of functional CN groups.     

The O2+ vibrational distribution in the Venusian ionosphere

Calculations are presented of the vibrational distribution of O₂⁺ in the Venusian ionosphere for a model atmosphere based on Pioneer Venus data. At 100 km, quenching precludes the survival of vibrationally excited O₂⁺. At the exobase, near 200 km, more than half are vibrationally excited. The effects of vibrationally excited O₂⁺ on the hot oxygen corona and the airglow are discussed.     

Effect of vertical plasma transport on ionospheric F2-region parameters at equatorial latitude

The variability of the F2-layer even during magnetically quiet times are fairly complex owing to the effects of plasma transport. The vertical E × B drift velocities (estimated from simplified electron density continuity equation) were used to investigate the seasonal effects of the vertical ion drifts on the bottomside daytime ionospheric parameters over an equatorial latitude in West Africa, Ibadan, Nigeria (Geographic: 7.4°N, 3.9°E, dip angle: 6°S) using 1 year of ionsonde data during International Geophysical Year (IGY) of 1958, that correspond to a period of high solar activity for quiet conditions. The variation patterns between the changes of the vertical ion drifts and the ionospheric F2-layer parameters, especially; f_oF₂ and h_mF₂ are seen remarkable. On the other hand, we observed strong anti-correlation between vertical drift velocities and h′F in all the seasons. We found no clear trend between N_mF₂ and h_mF₂ variations. The yearly average value of upward daytime drift at 300 km altitude was a little less than the generally reported magnitude of 20 ms⁻¹ for equatorial F-region in published literature, and the largest upward velocity was roughly 32 ms⁻¹. Our results indicate that vertical plasma drifts; ionospheric F2-layer peak height, and the critical frequency of F2-layer appear to be somewhat interconnected.     

Study of maximum electron density NmF2 at Karachi and Islamabad during solar minimum (1996) and solar maximum (2000) and its comparison with IRI

The monthly hourly medians of maximum electron density, N_mF₂, at two Pakistani ionospheric stations, Karachi and Islamabad, have been determined for solar minimum (1996) and solar maximum (2000) and compared with IRI predictions using the URSI coefficients. At night and pre-noon period the N_mF₂ values at both stations are almost equal during the 2 years. However, at post-noon the values at Karachi are considerably larger than those at Islamabad due to the equatorial or geomagnetic anomaly. Karachi (geomag. coord. 16.44°N, 139.08°E) lies near the region of the equatorial anomaly (+20 and −20 geomagnetic latitude), so most of the N_mF₂ values at Karachi are larger than those at Islamabad (geomag. coord. 24.46°N, 145.67°E). The maximum monthly values of N_mF₂ show a semi-annual variation at Karachi and Islamabad both during 1996 and 2000 as predicted by IRI.     

Implementation of a new characteristic parameter into the IRI sub-peak electron density profile

The shape of electron density profile in the International Reference Ionosphere could be improved significantly if the height h_g and electron density N_g of the F region sub-peak inflexion point were entered in the set of the profile standard parameters. To study variations of these important parameters, the N(h) analysis of the statistically-summarized ionograms at the latitudes of 40–80°N of the Eastern hemisphere has been carried out for the two-hours intervals of local time, three seasons (winter, equinox and summer) and two levels of solar activities characterized by Covington indicesF_10.7 = 100 and 200. It is shown that the parameters of the inflexion point can be expressed in most cases via the peak parameters of the F2 layer ash_g= 0.8 hmF2 and N_g= 0.5 NmF2.     

Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.     

Improvements in bottomside electron density definition in the Autoscala program

Some improvements introduced in the Autoscala program are presented. They include improvements in E valley modeling of the electron density profile N_e(h), and in the link between the E valley and bottom-side F regions. An abrupt variation in N_e(h) generated by the previous version of Autoscala under night conditions has been eliminated.A series of ionograms recorded by the Millstone Hill digisonde (42.6°, 288.5°) were automatically interpreted by the previous version of Autoscala and by the new one. Data from Incoherent Scatter Radar (ISR) were used to comparatively assess the performance of the two versions. For this purpose, the root mean square errors (RMSEs) of the N_e(h) provided by Autoscala were calculated relative to the corresponding values provided by ISR.A more accurate overall modeling of N_e(h) was achieved by the new Autoscala version (RMSE = 0.51 MHz for the new version against RMSE = 0.67 MHz for the previous one).     

Equatorial M(3000)F2 estimation of F2-layer models peak heights during low solar activity

M(3000)F2 estimation of h_mF2 based on four different formulated models viz: (1) Shimazaki (1955) (2) Bradley and Dudeney (1973), (3) Dudeney (1974) and (4) Bilitza et al. (1979) at an equatorial station in West Africa during low solar activity period (1995) are used to validate its conformity with observed and International Reference Ionosphere (IRI) model. Local time analyses of data from fifteen (15) selected days during the January and July solstices and April and October equinoxes are used. The results obtained show that the M(3000)F2 estimation of h_mF2 from the ionosonde-measured values using the Ionospheric Prediction Service (IPS-42) sounder compared to the observed values which were deduced using an algorithm from scaled virtual heights of quiet day ionograms are highly correlated with Bilitza model. International Reference Ionosphere (IRI 2007) model for the equatorial region also agrees with the formulation developed by Bilitza et al. (1979) for the four different seasons of the year. h_mF2 is highest (425 km) in summer (June solstice) season and lowest (386 km) in autumn (September equinox) season with daytimes peaks occurring at 1100–1200 LT during the solstices and at 1000 LT during the equinoxes respectively. Also, the post-sunset peaks are highest (362 km) at the spring (March equinox) and lowest (308 km) at the summer (June solstice) both occurring between 1800 and 2000 LT.     

Estimation of the elastic Earth parameters (h2, l2) using SLR data

We present results for the global elastic parameters h₂ and l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two satellites LAGEOS 1 and LAGEOS 2 observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analysed using different approaches. The analysis was done separately for the two satellites and approaches to estimate the two elastic parameters independently and together were performed. We do a sequential analysis and study the stability of the estimates as a function of length of the data set used. The adjusted final values for h₂ equal to 0.6151 ± 0.0008 and 0.6152 ± 0.0008, and those for l₂ equal to 0.0886 ± 0.0003 and 0.0881 ± 0.0003 for LAGEOS 1 and LAGEOS 2 tracking data are compared to other independently derived estimates. These parameters and their errors achieve stability at about the 24 and 27 month time interval for h₂ and l₂, respectively.     

Comparison of the CHAMP radio occultation data with the Canadian advanced digital ionosonde in the Polar Regions

In this paper we present the results of the comparison of the retrieved electron density profiles of the Ionospheric Radio Occultation (IRO) experiment on board CHAMP (CHAllenging Minisatellite Payload), with the ground ionosonde profiles for the Polar Regions. IRO retrieved electron density profiles from CHAMP are compared with Canadian Advanced Digital Ionosonde (CADI) measurements at two vertical sounding stations well within the Polar Cap, Eureka (geog. 80°13′ N; 86°11′ W) and Resolute Bay (geog. 74°41′ N; 94°54′ W). We compared the ionospheric parameters such as the peak electron density of the F-layer (NmF₂) and the peak height of the F-layer (h_mF₂) for a 3-year period, 2004–2006. CHAMP derived NmF₂ shows reasonable agreement with the ionosonde retrieved NmF₂ for both the stations (0.76 and 0.71 correlation coefficient, for Eureka and Resolute Bay, respectively) whereas the h_mF₂ agreement is not that acceptable (0.25 and 0.37 correlation coefficient, respectively). The h_mF₂ from vertical sounding showed less spread than the CHAMP h_mF₂.     

Satellite measurements of atmospheric composition: Three years' observations of CH4 and N2O

The stratospheric and mesospheric sounder (SAMS) was launched in October 1978 on the NIMBUS 7 satellite. Between then and its eventual failure in June 1983 the instrument was used to collect over four years of radiance data from which atmospheric temperature and the abundances of a number of minor constituents have been derived. The paper will present fields of CH₄ and N₂O between 50S and 70N derived from SAMS data for the period 1979–1981. Global distributions of CH₄ and N₂O will be presented in various forms and the observed seasonal changes and interannual variability will be described. The paper will compare the SAMS CH₄ and N₂O data with model predictions and will comment upon some other areas of interest.     

Thermosphere–ionosphere response to a severe magnetic storm: A case study

This paper reports the response of the ionosphere–thermosphere system to an intense geomagnetic storm. For that, data taken by instruments on board Dynamic Explorer 2 at heights of the F2-layer (molecular nitrogen N₂ and atomic oxygen O compositions, neutral temperature Tg and electron density Ne) were used. The ionospheric response is characterized by a negative storm effect expanding from mid–high to low latitude. It is observed during this severe geomagnetic storm that negative effects were caused mainly by an increase in molecular nitrogen composition N₂ and almost no changes in atomic oxygen composition O.     

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.     

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.     

Occurrence features of simultaneous H+- and He+-band EMIC emissions in the outer radiation belt

As an important loss mechanism of radiation belt electrons, electromagnetic ion cyclotron (EMIC) waves show up as three distinct frequency bands below the hydrogen (H⁺), helium (He⁺), and oxygen (O⁺) ion gyrofrequencies. Compared to O⁺-band EMIC waves, H⁺- and He⁺-band emissions generally occur more frequently and result in more efficient scattering removal of <～5?MeV relativistic electrons. Therefore, knowledge about the occurrence of these two bands is important for understanding the evolution of the relativistic electron population. To evaluate the occurrence pattern and wave properties of H⁺- and He⁺-band EMIC waves when they occur concurrently, we investigate 64 events of multi-band EMIC emissions identified from high quality Van Allen Probes wave data. Our quantitative results demonstrate a strong occurrence dependence of the multi-band EMIC emissions on magnetic local time (MLT) and L-shell to mainly concentrate on the dayside region of L?=?～4–6. We also find that the average magnetic field amplitude of H⁺-band waves is larger than that of He⁺-band waves only when L?<?4.5 and AE¹?<?300?nT, and He⁺-band emissions are more intense under all other conditions. In contrast to 5 events that have average H⁺-band amplitude over 2 nT, 19 events exhibit >2 nT He⁺-band amplitude, indicating that the He⁺-band waves can be more easily amplified than the H⁺-band waves under the same circumstances. For simultaneous occurrences of the two EMIC wave bands, their frequencies vary with L-shell and geomagnetic activity: the peak wave frequency of H⁺-band emissions varies between 0.25 and 0.8 f_cp with the average between 0.25 and 0.6 f_cp, while that of He⁺-band emissions varies between 0.03 and 0.23 f_cp with the average between 0.05 and 0.15 f_cp. These newly observed occurrence features of simultaneous H⁺- and He⁺-band EMIC emissions provide improved information to quantify the overall contribution of multi-band EMIC waves to the loss processes of radiation belt electrons.     

Estimation of the Love and Shida numbers: h2, l2 using SLR data for the low satellites

In this paper we present results for the global elastic parameters: Love number h₂ and Shida number l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two low satellites STELLA (H = 800 km) and STARLETTE (H = 810 km) observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analyzed. The analysis was done separately for the two satellites. We do a sequential analysis and study the stability and convergence of the estimates as a function of length of the data set used.