Proposed reference models for CO2 and halogenated hydrocarbons

The vertical distribution of carbon dioxide, halocarbons and their sink products, HCl and HF, have become available, mainly by means of balloon measurements. Most measurements were made at northern mid-latitudes, but some constituents were measured at tropical latitudes and in the southern hemisphere as well. This report attempts to combine the available data for presentation of reference models for CO₂, CCl₄, CCl₃F, CCl₂F₂, CClF₃, CF₄, CCl₂F-CClF₂, CClF₂-CClF₂, CClF₂-CF₃, CF₃-CF₃, CH₃Cl, CHClF₂, CH₃-CCl₃, CBrClF₂, CBrF₃, HCl and HF.     

Vertical profiles of source gases at tropical latitudes

A first ISRO-DFVLR collaborative balloon flight of the MPAE cryogenic sampler was conducted at Hyderabad, India (17.5°N), on March 27 1985, and 15 air samples were collected between 10 and 35 km altitude. Vertical profiles of source gases analysed from these samples, such as CCl₃F, CCl₂F₂, CBrClF₂, and CH₄ are presented. Due to tropical upwelling, mixing ratios of source gases decrease less rapidly with altitude than at mid-latitudes.     

Atmospheric sampling

Four important sampling techniques are briefly reviewed: Selective sampling on impregnated filters for measuring acidic gases, the matrix isolation technique for measuring radicals, whole air grabsampling and whole air cryogenic sampling for measuring stable source gases.Vertical profiles of H₂, CH₄, CO, N₂O, CFCl₃ and CF₂Cl₂ resulting from gas chromatographic analysis of whole air samples collected with a cryogenic sampler are presented. Year-to-year variations are observed for H₂, CH₄ and N₂O above 25 km, while CFCl₃ and CF₂Cl₂ mixing ratios show a noticeable increase between 1977 and 1979 at almost every height level.The CO₂ mixing ratio is not constant with height but rather decreases from 332 ppmV at 10 km to 325 ppmV at 30 km.The vertical distribution of methyl chloride is characterized by a rapid decrease from 600 pptV in the troposphere to less than 10 pptV at 32 km in agreement with model results.     

Stratospheric trace gas distributions observed in different seasons

A newly developed balloon-borne cryogenic sampler has been used to collect large air samples in the midlatitude stratosphere (44° N). The samples were analysed for their content of several longlived trace gases. The vertical profiles of N₂O, CH₄, and CF₂Cl₂ derived from two balloon flights in winter 1982 and, as part of the MAP Globus Project, in fall 1983 indicate variations that are caused by large scale dynamic processes. Significant features of the profiles are discussed and compared with averaged profiles that were derived from a series of balloon sampling flights performed previously over the same location during the summer season. The temporal variation of the individual profiles is more detailed than predicted by model calculations.     

Origin and stability of lunar polar volatiles

Temperature regime at the LCROSS impact site is studied. All detected species in the Cabeus crater as well as CH₄ and CO clathrate hydrates except H₂, CO, and CH₄ are stable against evaporation at the LCROSS impact site. CO and CH₄ can be chemisorbed at the surface of the regolith particles and exist in the form of clathrate hydrates in the lunar cold traps. Flux rates of delivery of volatile species by asteroids, micrometeoroids, O-rich, C-rich, and low-speed comets into the permanently shadowed regions are estimated. Significant amounts of H₂O, CO, H₂, H₂S, SO₂, and CO₂ can be impact-produced during collisions between asteroids and O-rich comets with the Moon while CH₃OH, NH₃ and complex organic species survive during low-speed comet impacts as products of disequilibrium processes. C-rich comets are main sources of CH₄, and C₂H₄.     

Statistical global and regional atmospheric models of temperature and gas components over the northern hemisphere

Statistical mid-latitude models of altitude distribution of temperature, water vapor, ozone, carbon dioxide and trace gases (CO, CH₄, N₂O, NO, NO₂) are considered. The mean characteristics of altitude profiles of these parameters, as well as their time and space variability, have been taken into account. The statistical regional models were constructed using a temperature-humidy complex. The considered statistical mid-latitude models have been constructed as applied to solutions of the problems on remote sounding of the atmosphere and underlying surface from outer space.     

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%.     

An icy mineralogy package (IMP) for in-situ studies of Titan’s surface

We describe the scientific case for and preliminary design of an instrument whose primary goal is to determine the chemistry (element abundance) and mineralogy (compound identity and abundance) of Titan’s surface using a combination of energy dispersive X-ray fluorescence spectroscopy (EDXRF) and X-ray diffraction (XRD). XRD is capable of identifying any crystalline substance present on Titan’s surface at relative abundances greater than ∼1 wt%, allowing unambiguous identification of, for example, structure I and II clathrates (even in the presence of ice), and various organic solids, which may include C₂H₂, C₂H₄, C₄H₂, HCN, CH₃CN, HC₃N, and C₄N₂). The XRF component of the instrument will obtain elemental abundances for 16 < Z < 60 with minimum detection limits better than 10 ppm (including detection of atmospheric noble gas isotopes), and may achieve detection limits of 0.01–1% for lighter elements down to Z = 6 (carbon). The instrument is well suited to integration with other analytical tools as part of a light-weight surface chemistry and mineralogy package. Although considerably less sensitive to elemental abundance than GC–MS (10⁻² vs. 10⁻⁸) it is likely to be significantly lighter (<0.5 kg vs. 10 kg).     

A new inversion algorithm for HF sky-wave backscatter ionograms

HF sky-wave backscatter sounding system is capable of measuring the large-scale, two-dimensional (2-D) distributions of ionospheric electron density. The leading edge (LE) of a backscatter ionogram (BSI) is widely used for ionospheric inversion since it is hardly affected by any factors other than ionospheric electron density. Traditional BSI inversion methods have failed to distinguish LEs associated with different ionospheric layers, and simply utilize the minimum group path of each operating frequency, which generally corresponds to the LE associated with the F₂ layer. Consequently, while the inversion results can provide accurate profiles of the F region below the F₂ peak, the diagnostics may not be so effective for other ionospheric layers. In order to resolve this issue, we present a new BSI inversion method using LEs associated with different layers, which can further improve the accuracy of electron density distribution, especially the profile of the ionospheric layers below the F₂ region. The efficiency of the algorithm is evaluated by computing the mean and the standard deviation of the differences between inverted parameter values and true values obtained from both vertical and oblique incidence sounding. Test results clearly manifest that the method we have developed outputs more accurate electron density profiles due to improvements to acquire the profiles of the layers below the F₂ region. Our study can further improve the current BSI inversion methods on the reconstruction of 2-D electron density distribution in a vertical plane aligned with the direction of sounding.     

Rotational constants of linear and/or bent Cn+1H and CnN(n = 1–6): A DFT study

The geometries, dipole moments, and rotational constants for the linear and/or bent cations, C_n+1H⁺ and C_nN⁺(n = 1–6), were studied by the B3LYP method with the modest basis sets. For C_nH⁺(n = odd; 3, 5, 7) and C_nN⁺(n = even; 2, 4, 6), the theoretical rotational constants (B_es) of closed-shell singlet C₃H⁺, C₅H⁺, C₇H⁺, CCN⁺, C₄N⁺, and C₆N⁺ were calculated to be about 11,244, 2420, 885.2, 11,970, 2439, and 880.8 MHz, respectively. By contrast, the triplets are stable than the corresponding singlets for C_nH⁺(n = odd; 2, 4, 6) and C_nN⁺(n = even; 3, 5) except CN⁺.     

Ionospheric electron temperature at solar maximum

Langmuir probe measurements made at solar maximum from the Dynamics Explorer-2 satellite in 1981 and 1982 are employed to examine the latitudinal variation of electron temperature, T_e, at altitudes between 300 and 400 km and its response to 27 day variations of solar EUV. Comparison of these data with T_e models based on the solar minimum measurements from Atmosphere Explorer-C suggest that the daytime T_e does not change very much during the solar cycle, except at low latitudes where an especially large 27 day variation occurs. The 27 day component decreases from about 7°/F10.7 unit at the equator to 3°/F10.7 unit at 851V 3 middle and higher latitudes. From these DE-2 measurements, and those from AE-C, we conclude that the daytime T_e near the F₂ peak is more responsive to short-term (daily) variations in F10.7 than to any longer term changes that may occur between solar minimum and solar maximum. To investigate this sensitivity of the dayside ionosphere to solar activity we employ the inverse relationship of T_e and N_e, that was found at solar minimum, to see if it can be used to order the T_e behaviour at solar maximum. We introduce a simple quadratic correction for the F10.7 influence on T_e based on the entire daytime AE-C and DE-2 data base between 300 and 400 km. Although this equation may be found useful, the systematic deviations of the DE-2 data suggest that the solar minimum model does not accurately describe the T_e-N_e relationships at solar maximum, at least above 300 km where the DE-2 measurements were made. Future work with this data base should attempt to see if such a relationship exists.     

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.     

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.     

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₂.     

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.     

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.     

Bottom side profiles for two close stations at the southern crest of the EIA: Differences and comparison with IRI-2012 and NeQuick2 for low and high solar activity

Bottom side electron density profiles for two stations at the southern crest of the Equatorial Ionization Anomaly (EIA), São José dos Campos (23.1°S, 314.5°E, dip latitude 19.8°S; Brazil) and Tucumán (26.9°S, 294.6°E, dip latitude 14.0°S; Argentina), located at similar latitude and separated by only 20° in longitude, have been compared during equinoctial, winter and summer months under low (year 2008, minimum of the solar cycle 23/24) and high solar activity (years 2013–2014, maximum of the solar cycle 24) conditions. An analysis of parameters describing the bottom side part of the electron density profile, namely the peak electron density NmF2, the height hmF2 at which it is reached, the thickness parameter B₀ and the shape parameter B₁, is carried out. Further, a comparison of bottom side profiles and F-layer parameters with the corresponding outputs of IRI-2012 and NeQuick2 models is also reported. The variations of NmF2 at both stations reveal the absence of semi-annual anomaly for low solar activity (LSA), evidencing the anomalous activity of the last solar minimum, while those related to hmF2 show an uplift of the ionosphere for high solar activity (HSA). As expected, the EIA is particularly visible at both stations during equinox for HSA, when its strength is at maximum in the South American sector. Despite the similar latitude of the two stations upon the southern crest of the EIA, the anomaly effect is more pronounced at Tucumán than at São José dos Campos. The differences encountered between these very close stations suggest that in this sector relevant longitudinal-dependent variations could occur, with the longitudinal gradient of the Equatorial Electrojet that plays a key role to explain such differences together with the 5.8° separation in dip latitude between the two ionosondes. Furthermore at Tucumán, the daily peak value of NmF2 around 21:00 LT during equinox for HSA is in temporal coincidence with an impulsive enhancement of hmF2, showing a kind of “elastic rebound” under the action of the EIA. IRI-2012 and NeQuick2 bottom side profiles show significant deviations from ionosonde observations. In particular, both models provide a clear underestimation of the EIA strength at both stations, with more pronounced differences for Tucumán. Large discrepancies are obtained for the parameter hmF2 for HSA during daytime at São José dos Campos, where clear underestimations made by both models are observed. The shape parameter B₀ is quite well described by the IRI-2012 model, with very good agreement in particular during equinox for both stations for both LSA and HSA. On the contrary, the two models show poor agreements with ionosonde data concerning the shape parameter B₁.     

A southern hemisphere error in the IRI

I would like to call attention to the fact that the IRI computes erroneously the F2-layer semithickness parameter B₀ at southern hemisphere locations. The values of B₀, based on northern latitude observations, have a seasonal variation which must be preserved at southern latitudes.The error was found in the course of a study to develop a new ionospheric model for radio-propagation predictions. We observed at southern latitudes major discrepancies between the IRI and the Bradley-Dudeney (1973) model in relation to the F2-layer semithickness. This is estimated in the latter model as the difference between the height of maximum electron concentration (hmF2) and the ionospheric characteristic h′F,F2, the minimum observed virtual height of reflection from the F2-layer, corrected taking into account underlying ionization. The profiles for both models were drawn using the same values of foF2 and hmF2. Then, our analysis served also to test the IRI model with h′F,F2 data obtained from CCIR maps but not used as primary inputs by the IRI.     

The infrared synthetic spectrum of comet Halley

In order to prepare infrared sounding of comet Halley from the flyby VEGA probes, we have computed the synthetic spectrum between 2.5 and 15 μ of a typical comet at a heliocentric distance of ～ 0.8 AU. The present paper is particularly devoted to the contribution from the cometary gases. For a selection of 20 possible parent molecules, the most efficient excitation process is resonant fluorescence by the solar radiation field. The H₂O, CO, CO₂, CH₄, NH₃ and H₂CO molecules are the best candidates for detection by the IKS infrared spectrometers aboard the VEGA probes. For the water molecule, collisions are too rare to ensure thermal equilibrium in the whole coma ; therefore a limited number of fluorescence lines are expected to be present in the H₂O vibrational bands.     

The influence of radiative cooling and turbulence on the heat budget of the thermosphere

Numerical models of the thermal budget of the Earth's upper atmosphere in the height range of 90–500km are developed. The main sources and sinks of energy including infra-red radiative cooling by vibrational-rotational bands of NO, CO₂, OH and O₃ as well as heating and cooling arising from dissipation of turbulent energy and eddy heat transport are taken into account. The calculated temperature and density height profiles are in good agreement with the respective profiles from CIRA 72 and Jacchia 1977 models. It is shown for the models considered that IR-radiative cooling by CO₂ and NO in the 15μ and 5.3μ bands, not eddy turbulence provides the major loss of heat from 90 to 180km.