Models of the high latitude stratosphere and mesosphere

Methods of calculation of interhemispheric and longitudinal differences are discussed. It is noted that in models one should take into account the solar activity influence. Network rocket sounding accuracy, devices and technique compatibility are analysed.     

Assessment of the ozone and temperature variability during 1979–1993 with the chemistry-climate model SOCOL

The chemistry-climate model SOCOL has been applied for the study of ozone and temperature anomalies during 1979–1993. Temperature and ozone anomalies have been obtained for a set of model runs forced by all major stratospheric forcing mechanisms. Forcings have been prescribed separately and together to assess their individual influence on stratospheric ozone and temperature. The results of these simulations have been compared to available satellite data. The model captures well ozone depletion and cooling in the upper stratosphere due to increases in the abundance of greenhouse gases and ozone depleting substances in the atmosphere. In the lower stratosphere, the model reproduces the warming over tropical and middle latitudes caused by the El-Chichon and Pinatubo eruptions. However, the simulated ozone response is overestimated in comparison with SAGE data. The best agreement with observations has been obtained for the run with all forcings included. This emphasizes the importance of the volcanic and solar forcings for the correct reproduction of observed trends. Comparison of near-global total ozone anomalies confirms an overestimation of ozone depletion just after volcanic eruptions, while the overall agreement with the model is fairly good.     

Wind structure and small-scale wind variability in the stratosphere and mesosphere during the November 1980 energy budget campaign

Between November 6 and December 1, 1980, a series of rocket observations obtained from two sites in northern Scandinavia as part of the Energy Budget Campaign indicated that significant vertical and temporal changes in the wind structure were present and were noted to coincide with different geomagnetic conditions, i.e., quiet and enhanced. This series of observations represents for the first time the largest amount of data ever gathered at high latitudes over such a short interval of time. It is observed that prior to November 16, the meridional wind component above 60 kilometers was found to be positive (southerly) while the magnitude of the zonal wind component increased with altitude. After November 16 the meridional component became negative (northerly) and the magnitude of the zonal wind component was noted to decrease with altitude. Time-sections of the perturbations of the zonal wind show the presence of vertically propagating waves which suggest gravity wave activity. These waves increase in wavelength from 3–4 kilometers near 40 kilometers to over 12 kilometers near 80 kilometers. The observational techniques employed at Andoya, Norway, and ESRANGE in Sweden, consisted of chaff foil, chemical trails, inflatable spheres, and parachutes.     

Height-latitude structure of stationary planetary waves in the stratosphere and lower mesosphere

Daily UK Met Office stratospheric assimilated data for the Northern and Southern Hemispheres, accumulated for the period from 2004 to 2012 and pressure range of 1000–0.1 hPa, are used in this paper. The paper presents and thoroughly discusses spatial–temporal distributions of stationary planetary wave (SPW) amplitudes and phases, calculated on the basis of geopotential height, temperature, zonal and meridional wind data for zonal wave numbers 1 and 2 (SPW1 and SPW2). The climatological planetary wave amplitudes and phases are calculated by extracting waves from three types of data: daily, monthly mean and climatological monthly mean. It has been established that magnitude of amplitudes and height-latitude distribution of amplitudes of SPW1 and SPW2 depend on data processing method for all parameters. It has also been established that height-latitude distribution amplitudes and phases significantly differ for geopotential height, temperature, zonal and meridional wind and depend on wave number and hemisphere. However, height-latitude distributions of phases are little different from each other for the used methods of data processing.     

A neutral gas mass spectrometer to measure the chemical composition of the stratosphere

The Polar Balloon Atmospheric Composition Experiment (P-BACE) is a new generation of neutral gas mass spectrometer based on the time-of-flight principle. P-BACE is the scientific experiment on the Mars Environment Analog Platform (MEAP) flown successfully on a balloon mission in summer 2008. The MEAP mission was flown with a 334,000 m³ helium balloon in the stratosphere on a semicircular trajectory from northern Sweden around the North Pole to Canada using the summer northern hemispheric wind current. The atmospheric conditions at an atmospheric altitude of 35–40 km are remarkably similar to those on the surface of Mars and thus the balloon mission was an ideal testbed for our mass spectrometer P-BACE. Originally this instrument was designed for in situ measurements of the chemical composition of the Martian atmosphere.P-BACE has a unique mass range from 0 to 1000 amu/q with a mass resolution m/Δm (FWHM) > 1000, and the dynamic range is at least six orders of magnitude. During this experiment, the acquisition of one mass spectrum is a sum of 65,535 single spectra, recorded in a time frame of 66 s.The balloon mission lasted 5 days and had successfully demonstrated the functionality of the P-BACE instrument during flight conditions. We had recorded more than 4500 mass spectra. With little modifications, P-BACE can be used on a planetary mission for Mars, but for example also for Venus or Mercury, if placed on a satellite.     

Excitations of nonmigrating diurnal tides in the mesosphere and lower thermosphere simulated by the Kyushu-GCM

Excitation mechanisms of nonmigrating diurnal tides in the MLT region simulated by the Kyushu-GCM are examined. It is shown that the westward propagating diurnal tide with zonal wavenumber s = 2 is mainly excited by nonlinear interactions between the migrating diurnal tide and the stationary planetary wave with zonal wavenumber s = 1, while the nonlinear excitation of the standing diurnal tide with zonal wavenumber s = 0 is less important than the excitation by tropospheric heating. Nonlinear interactions between the migrating diurnal tide and the stationary planetary wave with zonal wavenumber s = 2 are not dominant to excite the westward propagating diurnal tide with zonal wavenumber s = 3, and it is shown that the excitation by tropospheric heating is comparable to the nonlinear excitation. It is also shown that other nonmigrating diurnal tides are excited by tropospheric heating.     

The influence of interplanetary shocks on solar protons measured in the stratosphere.

Since the beginning of the 22nd solar cycle twenty solar proton events were observed by the regular balloon measurements of cosmic rays. Temporal changes of intensities and energy spectra of solar protons with energy 100-500 MeV were obtained. The strong influence of interplanetary shock waves on the proton flux characteristics near the Earth was observed. Possible effects of solar proton transport in the vicinity of shock fronts are discussed to explain the observational data.     

The structure of the solar cycle maximum phase in the galactic cosmic ray 11-year variation

Two phenomena connected with the maximum phase of the 11-year solar cycle in the galactic cosmic ray intensity – the change in the energy dependence of the intensity variations and the double-peak structure in the intensity modulation time profile – are considered for the last five solar cycles (Nos. 19–23). The distinct 22-year cycle in the magnitude of the so called energy hysteresis is observed.The periods of the solar cycle maximum phase in the galactic cosmic ray intensity, characterized by the specific energy dependence of the intensity, are estimated. It is found that the double-peak structures belonging to the solar cycle maximum phase and those around it are very similar both in the amplitude and in its energy dependence.     

Features of the 11-year variation of galactic cosmic rays in different periods of solar magnetic cycles

Variations of galactic cosmic ray intensity have been studied based on the neutron monitors and interplanetary magnetic field experimental data for different ascending and descending epochs of solar activity. The dependence of the diffusion coefficient on the cosmic ray particles rigidity R is stronger in the maxima epoch than in the minima epoch of solar activity. For the period of 1977–1981 (qA > 0) the diffusion coefficient for the minimum epoch is, χ_min ∝ R^{0.7 ± 0.04} and for the maximum χ_max ∝ R^{1.3 ± 0.05}; for the period of 1987–1990 (qA < 0), χ_min ∝ R^{0.8 ± 0.05} and χ_max ∝ R^{1.1 ± 0.04}. The exponents ν_y and ν_z of the power spectral density of the B_y and B_z components of the IMF in the region of the frequencies (10⁻⁶– 4 × 10⁻⁶) Hz are larger for the minimum epoch of 1987 (ν_y ≈ 2.0 and ν_z ≈ 1.93) than for the maximum epoch of 1990 (ν_y ≈ 1.43 and ν_z ≈ 1.27).     

Possible influence of solar radiation variability on the stratospheric temperature

Using daily temperature data available from radio-sonde measurements over Barajas (Madrid), La Coruña and Palma de Mallorca stations, for the time span 1971–1982 and an altitude range 100-30 mb, a study is made comparing temperatures at differents levels with the 10.7 cm flux in order to check whether radiation variability must be included in lower stratospheric models. At the latitude studied, stratospheric temperatures are uninfluenced by sudden warming phenomena avoiding difficulties of masking found in previous studies.     

Application of a modified generalized sail model on a solar sail with wrinkles

In this paper, we present an analysis of effect of wrinkles on the solar sail performance. We describe different analytical, semi-analytical and numerical approaches to the calculation of general large-scale curvature of a solar sail as well as parameters of so-called wrinkled domains, and introduce the impact of such wrinkles on the thrust and torque of the solar sail. Finally, we present a model of an optically-orthotropic surface for such non-ideal sail, providing a connection with the Generalized Sail Model, and other solar sail thrust models.     

Properties of the mesosphere and thermosphere and comparison with CIRA 72

Exospheric temperatures of several reference atmosphere are reviewed and a recommendation is made for the exospheric temperature of a proposed mean CIRA. One of the deficiencies of CIRA 72 and other present thermospheric models is the representation of density changes with geomagnetic activity. This deficiency is illustrated with samples of data. The data show the effects of geomagnetic activity, particle precipitation, a solar proton event, and gravity waves. An empirical model developed from the unique AFGL satellite density data bank using multiple linear regression is reviewed. The present model is for low to moderate solar flux and quiet geomagnetic conditions, but it is planned to extend the model to active conditions. Good progress has been made since CIRA 72 was specified in our knowledge and understanding of the properties of the lower thermosphere, although there are still some unresolved problems. The biggest progress has been made in the theory of tidal effects and of particulate energy deposition and of electrojet heating. On the other hand, it is still not possible to define adequately the systematic variations of the lower boundary conditions of thermospheric models. This is due to lack of knowledge of the systematic variations of the structure properties in the 100 to 120 km altitude region and inadequate information on the mesospheric turbulence profile and variations in the turbopause altitude.     

Influence of optical parameters on a solar sail motion

Detailed dynamic modeling of a solar sail requires recording of solar radiation pressure influence. A photon-solar sail is determined by the thrust value and the direction. We define the solar sail’s reflectivity depending on the film materials, the sail design and temperature, the thickness of multiple layers, and degradation factor, with a reasonable degree of accuracy. Thus, this work is devoted to the identification of optical characteristics of thin multilayer films in space flight conditions, i.e. to finding its reflectance, absorbance, and transmittance. In particular, the paper asks whether the solar sail simulates by a mathematical model of the optical characteristics of a multilayer epitaxial thin film. The temperature change effect and optical properties of solar sail degradation are considered as well. Solar sail flight from Earth to Mercury is designed as a simulation of the flight change in optical parameters.     

A global model of circulation and temperature for the upper mesosphere - lower thermosphere region

The zonal mean model of zonal and meridional wind for the Northern and Southern Hemisphere based the analysis of meteor radar wind and partial reflection drift data for the 70–110-km height interval is constructed. The height-latitudinal cross-sections of vertical wind are calculated from data on the latitudinal structure of a meridional wind field using the continuity equation. The temperature field cross-sections from the zonal wind model using the thermal wind equation are derived.     

Influence of curved thin-film device on deformation of a solar sail

A spinning solar sail IKAROS’s membrane is estimated to unexpectedly deform into an inverted pyramid shape due to thin-film devices with curvature, such as thin-film solar cells and steering devices on the membrane. It is important to investigate the deformation caused by the curved thin-film devices and predict the sail shape because the out-of-plane deformation greatly affects solar radiation pressure (SRP) and SRP torque. The purpose of this paper is to clarify the relationship between the global shape and orientation and position of curved thin-film devices and to evaluate SRP torque on the global shape using finite element analysis. The global shape is evaluated based on the out-of-plane displacement and the SRP torque. When the curved thin-film devices make the membrane shrink in the circumferential, diagonal, and radial direction, the sail deforms into a pyramid shape, an inverted pyramid one, and a saddle one, respectively. The saddle shape is more desirable for solar sails than the inverted pyramid shape and the pyramid one from the viewpoint of shape stability to SRP and control of SRP torque in the normal direction of the sail (windmill torque). The position of the thin-film device tends to increase the absolute value of windmill torque when it is biased circumferentially from the petal central axis. The suggested design principles for the arrangement of thin-film devices is that the curved thin-film devices should be directed so that the sail shrinks in the radial direction in order to deform the sail into a saddle shape with high shape stability, and the position of the thin-film devices should be biased in the circumferential direction paying attention to the absolute value of windmill torque to determine the direction of windmill torque.     

The physical mechanism of the solar variability influence on electrical and climatic characteristics of the troposphere

Possible mechanisms of solar–climatic connections, which may be of importance over short and long time intervals, are discussed. The variations of energetic balance of Earth’s climatic system for the last 50 years are estimated. It is ascertained that the imbalance between the flux of solar energy that comes to the Earth and radiates to space is of 0.1% for the last ten years. The significance is analyzed for the possible influence of variations of solar constant upon the energetic balance of the atmosphere.The physical mechanism of the influence of solar activity on climatic characteristics and the atmospheric circulation is suggested and theoretically substantiated. The mechanism is based on the redistribution in lower-troposphere of condensation nuclei by the vertical electric field. This electric field is determined by the ionosphere–Earth electric potential, which in the Polar Regions is controlled not only by tropical thunderstorms and by the galactic cosmic-ray intensity but also by solar cosmic-ray fluxes. The height redistribution in the atmosphere of condensation nuclei with a change of the electric field of the atmosphere is accompanied by a change in total latent heat (phase transition of water vapor), by changes in radiation balance, and by subsequent changes of the thermobaric field of troposphere. The results of analysis of thermobaric field variations for the periods of invasion of abnormally powerful solar cosmic ray fluxes and magnetic storms confirm the reality of manifestation of heliogeophysical disturbances.     

On the diurnal variations in the temperature and composition: a three-dimensional model with superrotation

Based on a simplified theoretical interpretation of the composition measurements with the ONMS and OIMS experiments on Pioneer Venus, the conclusion was drawn that the rotation rate of the thermosphere should be close (within a factor of two) to that of the lower atmosphere. A more realistic three-dimensional model of the thermosphere dynamics is now being developed, considering non-linear processes, higher order modes and collisional momentum exchange between the major species CO₂, CO and O, which describes the diurnal variations in temperature and composition (Niemann et al., $\text{JGR}$, 1980). The computed horizontal winds are about 300 m/sec near the terminators and poles. Results are also presented from a two-dimensional (quasi-axisymmetric) spectral model which describes the four day superrotation in the lower atmosphere of Venus.     

Patterns of carbon monoxide in the middle atmosphere and effects of solar variability

We determine the spatial-time patterns of zonally averaged carbon monoxide (CO) in the middle atmosphere by applying Principle Component Analysis to the CO data obtained from the Microwave Limb Sounder (MLS) measurements on the Aura satellite in 2004–2012. The first two principal components characterize more than 90% of the CO variability. Both principal components are localized in the low thermosphere near the mesopause. The first principal component is asymmetric relative to the poles. It has opposite signs in the Northern and Southern Hemisphere at mid to high latitudes and strongly oscillates with an annual periodicity. The second principal component has the same sign in both hemispheres and oscillates mainly with a semi-annual frequency. Both principal components are modulated by the 11-year solar cycle and display short-term variations. To test possible correlations of these variations with the short term solar ultraviolet (UV) variability we use the simultaneous measurements of the UV solar radiance from the Solar-Stellar Irradiance Comparison Experiment (SOLSTICE) on the Solar Radiation and Climate Experiment (SORCE) satellite to investigate the correlation between CO in the middle atmosphere and solar UV in 2004–2012. Using a wavelet coherence technique a weak, intermittent 27-day signal is detected in high-frequency parts of the CO principal components.     

Time step size and model stiffness in the simulated slew of a tow of square sails

In order to assess space tow solar sail stability and control feasibility, slew simulations are performed for a simplified but dynamically representative km-class tow-like sail of sixteen 25 m square units (10,000 m² total area and 110 kg gross mass) with a 250 kg payload. It is seen that, for the dimensions considered, the space tow concept is structurally sound and its control is feasible. While observed instabilities are identified as numeric in nature and are eliminated accordingly, their very occurrence highlights the need for a refinement of the model for future studies. The analyses are carried out with custom software implementing non-standard implicit-iterative time integration with innovative elements. A new damping model, specifically tailored for the analysis of truly gossamer systems such as solar sails, is also proposed.     

Effect of current sheets on the power spectrum of the solar wind magnetic field using a cell model

A puzzling observation of solar wind MHD turbulence is the often seen Kolmogorov scaling of k^-5/3$k^{-5/3}$, even though the solar wind MHD turbulence is dominated by Alfvénic fluctuations. Recently Li et al. (2011) proposed that the presence of current sheets may be the cause of the Kolmogorov scaling. Here, using a cell model of the solar wind we examine the effect of current sheets on the power spectrum of the solar wind magnetic field. We model the solar wind as multiple cells separated by current sheets. We prescribe the spectra of turbulent magnetic field in individual cells as IK-like and examine the spectra along trajectories that cross multiple boundaries. We find that these spectra become softer and are consistent with the Kolmogorov-scaling. Our finding supports our recent proposal of Li et al. (2011).