Aerosol in the upper layer of earth’s atmosphere

The upper atmospheric layer of Venus, Mars, Jupiter, Saturn, and earth contains an aerosol layer. The meteorites, rings, and removal of small planetary particles may be responsible for its appearance. The observations from 1979–1992 have shown that the optical aerosol thickness over the earth’s polar regions varies from τ ≈ 0.0002 to 0.1 to λ = 1 μm. The highest τ value was in 1984 and 1992 and was preceded by intense activity of the El Chichon (1982) and Pinatubo (1991) volcanoes. We have shown that increase in τ of the stratospheric aerosol may lead to decrease in ozone layer registered in the 1970s. The nature of the stratospheric aerosol (a real part of the refraction index), effective size particles r, and latitudinal variation τ remain unknown. The analysis of phase dependence of the degree of polarization is effective among the distal methods of determination of n _r and r. The observation value of intensity and degree of polarization in the visible light are caused by the optical surface properties and optical atmospheric thickness, whose values varied with latitude, longitude, and in time. Thus, it is impossible to correctly distinguish the contribution of the stratospheric aerosol. In UV-rays (λ < 300 nm), the ozone layer stops the influence of the surface and earth’s atmosphere up to height of 20–25 km. In this spectrum area, the negative factors are emission of various depolarizating gases, horizontal heterogeneity of the effective optical height of the ozone layer, and oriented particles indicated by variation of the polarization plane.     

Estimation of secular density variations in the upper atmosphere from 1964–2007 satellite drag data

Orbital parameters of several artificial satellites of the Earth were analyzed for 1964–2007 and secular variations of the atmospheric density were estimated for the last 30–40 years. The analysis was based on the information about orbital parameters of 17 satellites and high-precision numerical integrations of the equations of motion with allowance for basic perturbing factors and spatiotemporal density variations, calculated from measured solar activity indices using the NRLMSISE-00 atmosphere model. The results demonstrate the presence of long-term variations in the atmospheric density not presented in modern atmosphere models. During solar-activity cycle 21, the atmospheric density became 0.4 to 19% higher (depending on height) than in cycle 20. It decreased by 1.0 to 11% (depending on height) in cycle 22 as compared to cycle 21. Both decreases and increases were observed in the atmospheric density during cycle 23, but with much smaller gradients. The results cannot be explained only by the growing concentration of greenhouse gases. Possible causes of the density variations and possible ways to take them into account in modern empirical and semiempirical atmospheric models are discussed.     

The density of the upper martian atmosphere measured by Lyman-α absorption with Mars Express SPICAM

Absorption of interplanetary Lyman-α emission by Mars’ nightside lower thermosphere was observed by Mars Express Spectrometer for Investigation of Characteristics of the Atmosphere of Mars (SPICAM), and is analyzed to derive the CO₂ density at 110 km during a martian year. The observed density seasonal variability is consistent with recent observations obtained by stellar occultations, proving that this method, though not as accurate as stellar occultations could be used complementary to them to characterize large variations of thermospheric density on Mars and provide a better spatial coverage by Lyman-α imagery.     

Influence of the hot oxygen corona on the satellite drag in the Earth’s upper atmosphere

Calculation results on the possible influence of the hot oxygen fraction on the satellite drag in the Earth’s upper atmosphere on the basis of the previously developed theoretical model of the hot oxygen geocorona are presented. Calculations have shown that for satellites with orbits above 500 km, the contribution from the corona is extremely important. Even for the energy flux Q ₀ = 1 erg cm⁻² s⁻¹, the contribution of the hot oxygen can reach tens of percent; and considering that real energy fluxes are usually higher, one can suggest that for extreme solar events, the contribution of hot oxygen to the atmospheric drag of the satellite will be dominant. For lower altitudes, the contribution of hot oxygen is, to a considerable degree, defined by the solar activity level. The calculations imply that for the daytime polar atmosphere, the change of the solar activity level from F _10.7 ∼ 200 to F _10.7 ∼ 70 leads to an increase in the ratio of the hot oxygen partial pressure to the thermal oxygen partial pressure by a factor of almost 30, from 0.85 to 25%. The transition from daytime conditions to nighttime conditions almost does not change the contribution from suprathermal particles. The decrease of the characteristic energy of precipitating particles, i.e., for the case of charged particles with a softer energy spectrum, leads to a noticeable increase of the contribution of the suprathermal fraction, by a factor of 1.5–2. It has been ascertained that electrons make the main contribution to the formation of the suprathermal fraction; and with the increase of the energy of precipitating electrons, the contribution of hot oxygen to the satellite drag also increases proportionally. Thus, for a typical burst, the contribution of the suprathermal fraction is 30% even at relatively high solar activity F _10.7 = 135.     

Experimental reconstruction of Lomonosov’s discovery of Venus’s atmosphere with antique refractors during the 2012 transit of Venus

In 1761, the Russian polymath Mikhail Vasilievich Lomonosov (1711–1765) discovered the atmosphere of Venus during its transit over the Sun’s disc. In this paper we report on experimental reenactments of Lomonosov’s discovery with antique refractors during the transit of Venus June 5–6, 2012. We conclude that Lomonosov’s telescope was fully adequate to the task of detecting the arc of light around Venus off the Sun’s disc during ingress or egress provided proper experimental techniques as described by Lomonosov in his 1761 report are employed.     

First ever in situ observations of Venus’ polar upper atmosphere density using the tracking data of the Venus Express Atmospheric Drag Experiment (VExADE)

On its highly elliptical 24 h orbit around Venus, the Venus Express (VEX) spacecraft briefly reaches a periapsis altitude of nominally 250 km. Recently, however, dedicated and intense radio tracking campaigns have taken place in August 2008, October 2009, February and April 2010, for which the periapsis altitude was lowered to the 186–176 km altitude range in order to be able to probe the upper atmosphere of Venus above the North Pole for the first time ever in situ. As the spacecraft experiences atmospheric drag, its trajectory is measurably perturbed during the periapsis pass, allowing us to infer total atmospheric mass density at the periapsis altitude. A Precise Orbit Determination (POD) of the VEX motion is performed through an iterative least-squares fitting process to the Doppler tracking data, acquired by the VEX radioscience experiment (VeRa). The drag acceleration is modelled using an initial atmospheric density model (VTS3 model, Hedin, A.E., Niemann, H.B., Kasprzak, W.T., Seiff, A. [1983]. J. Geophys. Res. 88, 73–83). A scale factor of the drag acceleration is estimated for each periapsis pass, which scales Hedin’s density model in order to best fit the radio tracking data. Reliable density scale factors have been obtained for 10 passes mainly from the second (October 2009) and third (April 2010) VExADE campaigns, which indicate a lower density by a factor of about 1.8 than Hedin’s model predicts. These first ever in situ polar density measurements at solar minimum have allowed us to construct a diffusive equilibrium density model for Venus’ thermosphere, constrained in the lower thermosphere primarily by SPICAV-SOIR measurements and above 175 km by the VExADE drag measurements (Müller-Wodarg et al., in preparation). The preliminary results of the VExADE campaigns show that it is possible to obtain with the POD technique reliable estimates of Venus’ upper atmosphere densities at an altitude of around 175 km. Future VExADE campaigns will benefit from the planned further lowering of VEX pericenter altitude to below 170 km.     

Explosion of a comet nucleus fragment in the earth’s atmosphere

This paper analyzes data on thermal explosions of large meteoroids in the earth’s atmosphere. The cumulative function of flux of space bodies is corrected with regard to the explosion height, which is determined, according to our approach, by maximum braking. As a result, the integral function of flux in the work [Brown, P., Spalding, R.E., ReVelle, D.O., et al., The Flux of Small Near-Earth Objects Colliding with the Earth, Nature, 2002, vol. 420, pp. 314–316] is consistent with the one we derived earlier. It is found that at least one phenomenon of those discussed in the paper by Brown et al. is a result of explosion of a comet nucleus fragment. It is shown that the Tunguska phenomenon cannot be explained within a monolithic body model.     

Gamma-ray bursts and the production of cosmogenic radionuclides in the Earth’s atmosphere

An explanation is offered for the impulsive increase in the concentration of cosmogenic radiocarbon in annual tree rings (Δ¹⁴C ～ 12‰) from AD ?775. A possible cause of such an increase could be the high-energy emission from a Galactic gamma-ray burst. It is shown that such an event should not lead to an increase in the total production of ¹⁰Be in the atmosphere, as distinct from the effect of cosmic-ray fluxes on the atmosphere. At the same time, the production of an appreciable amount of ³⁶Cl, which can be detected in Greenland and Antarctica ice samples of the corresponding age, should be expected. This allows the effects caused by a gamma-ray burst and anomalously powerful proton events to be distinguished.     

Neutral atmosphere near the icy surface of Jupiter’s moon Ganymede

The paper discusses the formation and dynamics of the rarefied gas envelope near the icy surface of Jupiter’s moon Ganymede. Being the most massive icy moon, Ganymede can form a rarefied exosphere with a relatively dense near-surface layer. The main parent component of the gas shell is water vapor, which enters the atmosphere due to thermal degassing, nonthermal radiolysis, and other active processes and phenomena on the moon’s icy surface. A numerical kinetic simulation is performed to investigate, at the molecular level, the formation, chemical evolution, and dynamics of the mainly H₂O- and O₂-dominant rarefied gas envelopes. The ionization processes in these rarefied gas envelopes are due to exposure to ultraviolet radiation from the Sun and the magnetospheric plasma. The chemical diversity of the icy moon’s gas envelope is attributed to the primary action of ultraviolet solar photons and plasma electrons on the rarefied gas in the H₂O- or O₂-dominant atmosphere. The model is used to calculate the formation and development of the chemical diversity in the relatively dense near-surface envelope of Ganymede, where an important contribution comes from collisions between parent molecules and the products of their photolysis and radiolysis.     

Dynamics of the Earth’s core from VLBI observations

The Earth’s rotation is accompanied by free circadian oscillations of its liquid core in the inner cavity of the lower mantle, which perturb the angular momentum of the entire Earth and produce an additional free nutation of the celestial pole called free core nutation (FCN). Since this nutation causes resonances in the diurnal tides and in the expansions of luni—solar nutation, its study, especially an improvement of the FCN period, is of fundamental importance for the theory of the Earth’s rotation. We have determined the FCN parameters from a joint analysis of equidistant series of coordinates of the celestial pole obtained from the combined processing of VLBI observations on global networks of stations for the interval 1984.0–2008.4 by IERS (International Earth Rotation and Reference System Service, Paris, France) and NEOS (National Earth Orientation Service, Washington, USA). Applying a moving least-squares filter (MLSF) to these data has shown that the FCN period averaged over this time interval differs significantly from the theoretical one and its phase varies over a wide range. Using the mean quadratic collocation (MQC) method, we have obtained a new, more accurate stochastic FCN model. Its analysis by the envelope method has revealed long-term linear phase trends, calling into question not only the adopted FCN period but also its stability and, hence, the stability of the resonant effects in the Earth’s luni—solar nutation.     

Effects of the solar eclipse of August 1, 2008, on the earth’s lower atmosphere

The paper presents results of optical observations and analysis of dynamics of effects on the earth’s lower atmosphere of the partial solar eclipse (of magnitude 42%) of August 1, 2008, near Kharkov. This is compared with the effects induced by the partial solar eclipses on August 11, 1999, and October 3, 2005. All three eclipses occurred around midday. The standard deviation of the solar-limb displacement σ _S during the eclipses on October 3, 2005, August 1, 2008, and August 11, 1999, was established to decrease by 0.13, 0.30, and 0.68″ at the maximum of the solar obscuration function 0.13, 0.31, and 0.73, respectively, so that the temperature drop in the earth’s lower atmosphere t _a was 1.3, 2.0, and 7.3 K. The time lags of decreases of σ _S and t _a was found to be 15 and 5 minutes.     

On the probable change of the radius and nature of aerosol particles in the deep layers of Jupiter’s atmosphere

When analyzing the pressure dependences of the aerosol volume scattering coefficient calculated from the measurement data on the geometric albedo of Jupiter obtained in 1993 in the methane absorption bands at 619, 727, and 842 nm, the signs of probable changes in the parameters of aerosol particles in the deep atmospheric layers were detected and the first estimates of the magnitude of these changes were obtained. It has been found that, in the pressure interval from 4 to 14 bar, the effective radius of particles may increase twofold and more (larger than 0.73 μm) and the real part of the refractive index may grow by 10% (from 1.44 and higher) relative to the values of these parameters in the upper atmosphere. If we take into account these changes, we find no signs of aerosol deep in the atmosphere of Jupiter.     

A Nonequilibrium Figure of Saturn’s Satellite Iapetus and the Origin of the Equatorial Ridge on Its Surface

The structure, dynamical equilibrium, and evolution of Saturn’s moon Iapetus are studied. It has been shown that, in the current epoch, the oblateness of the satellite ε₂ ≈ 0.046 does not correspond to its angular velocity of rotation, which causes the secular spherization behavior of the ice shell of Iapetus. To study this evolution, we apply a spheroidal model, containing a rock core and an ice shell with an external surface ε₂, to Iapetus. The model is based on the equilibrium finite-difference equation of the Clairaut theory, while the model parameters are taken from observations. The mean radius of the rock core and the oblateness of its level surface, ε₁ ≈ 0.028, were determined. It was found that Iapetus is covered with a thick ice shell, which is 56.6% of the mean radius of the figure. We analyze a role of the core in the evolution of the shape of a gravitating figure. It was determined that the rock core plays a key part in the settling of the ice masses of the equatorial bulge, which finally results in the formation of a large circular equatorial ridge on the surface of the satellite. From the known mean altitude of this ice ridge, it was found that, in the epoch of its formation, the rotation period of Iapetus was 166 times shorter than that at present, as little as T ≈ 11^h27^m. This is consistent with the fact that a driving force of the evolution of the satellite in our model was its substantial despinning. The model also predicts that the ice ridge should be formed more intensively in the leading (dark and, consequently, warmer) hemisphere of the satellite, where the ice is softer. This inference agrees with the observations: in the leading hemisphere of Iapetus, the ridge is actually high and continuous everywhere, while it degenerates into individual ice peaks in the opposite colder hemisphere.     

Determination of coronal magnetic fields from 10 to 26R ⊙ using the density compression ratios of CME-driven shocks

Recently, the estimation of coronal magnetic field using new methods, such as standoff distance method or density compression ratio method has been reported. In the present work, we utilized the density compression ratio of CME-driven shocks for 10 events at 29 different locations in the upper solar corona (10–26R _⊙) and determined the coronal magnetic field for two different adiabatic indices (γ=4/3 and 5/3). In addition, radial dependence of shock parameters in the corona is studied. It is found that the magnetic field estimated in the above range agree with the general trend. In addition, we obtained a radial profile of magnetic field [B(R)=623R ^?1.4] in the entire upper corona (3–30R _⊙) by combining the magnetic field estimated by Kim et al. (Astrophys. J. 746:118, 2012) in the range 3–15R _⊙ and that estimated in the present study in the range (10–26R _⊙). The power-law indices are nearly in agreement with recent results of CME-driven shocks reported in the literature. The results are discussed with the comparison of newly reported coronal magnetic field values obtained by different techniques and found that the power-law relation closely follow the literature values.     

Corrected spectral dependence of the imaginary part of the refractive index of aerosol in Jupiter’s atmosphere in the short-wavelength spectral range

To correctly determine the relative contribution of aerosol to the scattering properties of a gas–aerosol medium in the continuum, we propose a method that allows more reliable values of the imaginary part of the refractive index n _i to be obtained for Jupiter’s atmosphere in the short-wavelength spectral range. We considered the measurement data on the spectral values of the geometric albedo of Jupiter acquired in 1993 and used the model of homogeneous spherical aerosol particles. The obtained values of n _i are 0.00378, 0.00309, 0.00254, 0.00175, 0.00123, 0.00084, 0.00064, 0.00045, 0.00031, 0.00033, 0.00013, and 0.00008 at wavelengths λ = 320, 350, 375, 400, 420, 450, 470, 500, 520, 550, 606, and 631 nm, respectively.     

Io’s atmosphere: Constraints on sublimation support from density variations on seasonal timescales using NASA IRTF/TEXES observations from 2001 to 2010

We present an analysis of 19 μm spectra of Io’s SO₂ atmosphere from the TEXES mid-infrared high spectral resolution spectrograph on NASA’s Infrared Telescope Facility, incorporating new data taken between January 2005 and June 2010 and a re-analysis of earlier data taken from November 2001 to January 2004. This is the longest set of contiguous observations of Io’s atmosphere using the same instrument and technique thus far. We have fitted all 16 detected blended absorption lines of the ν₂ SO₂ vibrational band to retrieve the subsolar values of SO₂ column abundance and the gas kinetic temperature. By incorporating an existing model of Io’s surface temperatures and atmosphere, we retrieve sub-solar column densities from the disk-integrated data. Spectra from all years are best fit by atmospheric temperatures <150 K. Best-fit gas kinetic temperatures on the anti-Jupiter hemisphere, where SO₂ gas abundance is highest, are low and stable, with a mean of 108 (±18) K. The sub-solar SO₂ column density between longitudes of 90–220° varies from a low of 0.61 (±0.145) × 10^?17 cm^?2, near aphelion in 2004, to a high of 1.51 (±0.215) × 10¹⁷ cm^?2 in 2010 when Jupiter was approaching its early 2011 perihelion. No correlation in the gas temperature was seen with the increasing SO₂ column densities outside the errors.Assuming that any volcanic component of the atmosphere is constant with time, the correlation of increasing SO₂ abundance with decreasing heliocentric distance provides good evidence that the atmosphere is at least partially supported by frost sublimation. The SO₂ frost thermal inertias and albedos that fit the variation in atmospheric density best are between 150–1250 W m^?2 s^?1/2 K^?1 and 0.613–0.425 respectively. Photometric evidence favors albedos near the upper end of this range, corresponding to thermal inertias near the lower end. This relatively low frost thermal inertia produces larger amplitude seasonal variations than are observed, which in turn implies a substantial additional volcanic atmospheric component to moderate the amplitude of the seasonal variations of the total atmosphere on the anti-Jupiter hemisphere. The seasonal thermal inertia we measure is unique both because it refers exclusively to the SO₂ frost surface component, and also because it refers to relatively deep subsurface layers (few meters) due to the timescales of many years, while previous studies have determined thermal inertias at shallower levels (few centimeters), relevant for timescales of ～2 h (eclipse) or ～2 days (diurnal curves).     

Distribution of HCN in Titan’s upper atmosphere from Cassini/VIMS observations at 3 μm

Cassini/VIMS limb observations have been used to retrieve vertical profiles of hydrogen cyanide (HCN) from its 3 μm emission in the region from 600 to 1100 km altitude at daytime. While the daytime emission is large up to about 1100 km, it vanishes at nighttime at very low altitudes, suggesting that the daytime emission originates under non-LTE conditions. The spectrally integrated radiances around 3.0 μm shows a monotonically decrease with tangent altitude, and a slight increase with solar zenith angle in the 40-80° interval around 800 km.A sophisticated non-LTE model of HCN energy levels has been developed in order to retrieve the HCN abundance. The population of the HCN 0 0⁰ 1 energy level, that contributes mostly to the 3.0 μm limb radiance, has been shown to change significantly with the solar zenith angle (SZA) and HCN abundance. Also its population varies with the collisional rate coefficients, whose uncertainties induced errors in the retrieved HCN of about 10% at 600-800 km and about 5% above. HCN concentrations have been retrieved from a set of spectra profiles, covering a wide range of latitudes and solar zenith angles, by applying a line-by-line inversion code. The results show a significant atmospheric variability above ∼800 km with larger values for weaker solar illumination. The HCN shows a very good correlation with solar zenith angles, irrespective of latitude and local time, suggesting that HCN at these high altitudes is in or close to photochemical equilibrium. A comparison with UVS and UVIS measurements show that these are close to the lower limit (smaller SZAs) of the VIMS observations above 750 km. However, they are in reasonable agreement when combining the rather large UV measurement errors and the atmospheric variability observed in VIMS. A comparison of the mean profile derived here with the widely used profile reported by Yelle and Griffith (Yelle R.V., Griffith, C.A. [2003]. Icarus 166, 107-115) shows a good agreement for altitudes ranging from 850 to 1050 km, while below these altitudes our result exhibits higher concentrations.     

Theoretical calculation of the interannual variability of the Earth’s insolation with daily resolution

Based on the astronomical ephemerides DE-406, theoretical calculations have been performed of the interannual variability of the Earth’s insolation related to celestial-mechanical processes for 365 points of a tropical year in the time period from 1900 to 2050. It has been determined that the average amplitude of variations of the interannual insolation is 0.310 W/m² (0.023% of the solar constant). The calculated variations are characterized by strict periodicity that corresponds with the length of a synodic month. Connection between the extreme values of the calculated insolation variability and syzygies has been defined. The average amplitude of the calculated variability exceeds by 1.7 times (0.01% of the solar constant) the amplitude of the interannual variability in the 11-year variation of the total Earth’s insolation.     

Influence of magnetic field on propagation of five-minute oscillations in the sun’s atmosphere: Phase shifts

From results of spectral (in Ba II λ 455.4-nm line) and spectropolarimetric (in Fe I λλ 1564.3–1565.8-nm lines) observations of the active region (an isolated faculae at the solar disk center) with the German vacuum tower telescope (VTT) at the Institute of Astrophysics on the Canary Islands, the peculiarities of propagation of five-minute oscillations from the photosphere base (h = 0 km) to the lower chromosphere (h = 650 km) were investigated. At the height of the continuum formation (h = 0 km), the nature of wave propagation in the active region does not differ much from that in the quiet region: 80–90% of the investigated areas are occupied by waves moving up and down. In the lower chromosphere (h = 650 km), differences in the behavior of the waves are fundamental. In a quiet area, the waves become standing for 90% of the cases. In contrast to this, in the presence of moderate and strong magnetic fields (B = 30–180 mT), in 47% of the cases, the waves are running upward, which gives the principal possibility to heat the active region. The investigations revealed the presence of the waves in the active region, for which the phase shift Φ _T,V of the temperature and velocity oscillations is between ?90° and 0°. These waves cannot propagate in a quiet atmosphere.     

Probing Milky Way’s hot gas halo density distribution using the dispersion measure of pulsars

A number of recent studies indicates a significant amount of ionized gas in a form of the hot gas halo around the Milky Way. The halo extends over the region of 100 kpc and may be acountable for the missing baryon mass. In this paper we calculate the contribution of the proposed halo to the dispersion measure (DM) of the pulsars. The Navarro, Frenk, and White (NFW), Maller and Bullock (MB), and Feldmann, Hooper, and Gnedin (FHG) density distibutions are considered for the gas halo. The data set includes pulsars with the distance known independently from the DM, e.g., pulsars in globular clusters, LMC, SMC and pulsars with known parallax. The results exclude the NFW distribution for the hot gas, while the more realisticMB and FHG models are compatible with the observed dispersion measure.