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.     

The Ashima/MIT Mars GCM and argon in the martian atmosphere

We investigate the ability of modern general circulation models (GCMs) to simulate transport in the martian atmosphere using measurements of argon as a proxy for the transport processes. Argon provides the simplest measure of transport as it is a noble gas with no sinks or sources on seasonal timescales. Variations in argon result solely from ‘freeze distillation’, as the atmosphere condenses at the winter poles, and from atmospheric transport. Comparison of all previously published models when rescaled to a common definition of the argon enhancement factor (EF) suggest that models generally do a poor job in predicting the peak enhancement in southern winter over the winter pole – the time when the capability of the model transport approaches are most severely tested. Despite observed peak EF values of ～6, previously published model predictions peaked at EF values of only 2–3. We introduce a new GCM that provides a better treatment of mass conservation within the dynamical core, includes more sophisticated tracer transport approaches, and utilizes a cube–sphere grid structure thus avoiding the grid-point convergence problem at the pole that exists for most current Mars GCMs. We describe this model – the Ashima Research/Massachusetts Institute of Technology Mars General Circulation Model (Ashima/MIT Mars GCM) and use it to demonstrate the significant sensitivity of peak EF to the choices of transport approach for both tracers and heat. We obtain a peak EF of 4.75 which, while over 50% higher than any prior model, remains well short of the observed value. We show that the polar EF value in winter is primarily determined by the competition between two processes: (1) mean meridional import of lower-latitude air not enriched in argon and (2) the leakage of enriched argon out of the polar column by eddies in the lowest atmospheric levels. We suggest possibilities for improving GCM representation of the CO₂ cycle and the general circulation that may further improve the simulation of the argon cycle. We conclude that current GCMs may be insufficient for detailed simulation of transport-sensitive problems like the water cycle and potentially also the dust cycle.     

Small-scale variability in Saturn’s lower ionosphere

We perform and present a wavelet analysis on all 31 Cassini electron density profiles published to date (Nagy, A.F. et al. [2006]. J. Geophys. Res. 111 (A6), CiteID A06310; Kliore, A.J. et al. [2009]. J. Geophys. Res. 114 (A4), CiteID A04315). We detect several discrete scales of variability present in the observations. Small-scale variability (S < 700 km) is observed in almost all data sets at different latitudes, both at dawn and dusk conditions. The most typical scale of variability is 300 km with scales between 200 km and 450 km being commonly present in the vast majority of the profiles. A low latitude dawn/dusk asymmetry is noted in the prevalent scales with the spectrum peaking sharply at the 300 km scale at dusk conditions and being broader at dawn conditions. Compared to dawn conditions the dusk ionosphere also shows more significant variability at the 100 km scale. The 300 km vertical scale is also present in the few available profiles from the northern hemisphere. Early observations from 2005 show a dominant scale at 350 km whereas later in 2007–2008 the spectrum shifts to the shorter scales with the most prominent scale being 300 km. The performed wavelet analysis and the obtained results are independent of assumptions about the nature of the layers and do not require a definition for a “background” electron density profile.In the second part of the paper we present a gravity wave propagation/dissipation model for Saturn’s upper atmosphere and compare the wave properties to the characteristics of the observed electron density variability at different scales. The general features observed in the data are consistent with gravity waves being present in the lower ionosphere and causing layering of the ions and the electrons. The wave-driving mechanism provides a simultaneous explanation for several of the properties of the observed variability: (i) lack of variability in the electron density above the predicted region of wave dissipation; (ii) in most cases the peak amplitude of variability occurs within the altitude range for dissipation of gravity waves or below; (iii) shorter scales have smaller amplitudes than the longer scales; (iv) shorter scales are present at lower altitudes whereas longer scales persist to higher altitudes; and (v) several layers often form a system of equally spaced maxima and minima that can be traced over a large altitude range.     

N-body integrators with individual time steps from Hierarchical splitting

We review the implementation of individual particle time-stepping for N-body dynamics. We present a class of integrators derived from second order Hamiltonian splitting. In contrast to the usual implementation of individual time-stepping, these integrators are momentum conserving and show excellent energy conservation in conjunction with a symmetrized time step criterion. We use an explicit but approximate formula for the time symmetrization that is compatible with the use of individual time steps. No iterative scheme is necessary. We implement these ideas in the HUAYNO¹ code and present tests of the integrators and show that the presented integration schemes shows good energy conservation, with little or no systematic drift, while conserving momentum and angular momentum to machine precision for long term integrations.     

Magnetic field measurements of CP stars from hydrogen line cores

We present the results of measurements of magnetic fields of chemically peculiar (CP) stars, performed from the shifts between the circularly polarized components of metal and hydrogen lines (the Babcock method). The observations are carried out with an analyzer of circular polarization at the 6‐m telescope of the SAO RAS. We found that for the absolute majority of the objects studied (in 22 CP stars out of 23), the magnetic fields, determined from the Zeeman shifts in the hydrogen line cores, are significantly lower than those obtained from metal lines in the same spectra. This disparity varies between the stars. We show that instrumental effects can not produce the above features, and discuss the possible causes of the observed effect. The discovered condition reveals a more complicated structure of magnetic fields of CP stars than a simple dipole, in particular, a reduction of the field strength in the upper atmosphere with the vertical gradient, significantly higher than the dipole (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A technique for processing of planetary images with heterogeneous characteristics for estimating geodetic parameters of celestial bodies with the example of Ganymede

The new technique for generation of coordinate control point networks based on photogrammetric processing of heterogeneous planetary images (obtained at different time, scale, with different illumination or oblique view) is developed. The technique is verified with the example for processing the heterogeneous information obtained by remote sensing of Ganymede by the spacecraft Voyager-1, -2 and Galileo. Using this technique the first 3D control point network for Ganymede is formed: the error of the altitude coordinates obtained as a result of adjustment is less than 5 km. The new control point network makes it possible to obtain basic geodesic parameters of the body (axes size) and to estimate forced librations. On the basis of the control point network, digital terrain models (DTMs) with different resolutions are generated and used for mapping the surface of Ganymede with different levels of detail (Zubarev et al., 2015b).     

Are galaxies with active galactic nuclei a transition population?

We present the results of an analysis of a well-selected sample of galaxies with active and inactive galactic nuclei from the Sloan Digital Sky Survey, in the range  0.01 < z < 0.16  . The SDSS galaxy catalogue was split into two classes of active galaxies, Type 2 active galactic nuclei (AGN) and composites, and one set of inactive, star-forming/passive galaxies. For each active galaxy, two inactive control galaxies were selected by matching redshift, absolute magnitude, inclination, and radius. The sample of inactive galaxies naturally divides into a red and a blue sequence, while the vast majority of AGN hosts occur along the red sequence. In terms of Hα equivalent width (EW), the population of composite galaxies peaks in the valley between the two modes, suggesting a transition population. However, this effect is not observed in other properties such as the colour–magnitude space or colour–concentration plane. Active galaxies are seen to be generally bulge-dominated systems, but with enhanced Hα emission compared to inactive red-sequence galaxies. AGN and composites also occur in less dense environments than inactive red-sequence galaxies, implying that the fuelling of AGN is more restricted in high-density environments. These results are therefore inconsistent with theories in which AGN host galaxies are a 'transition' population. We also introduce a systematic 3D spectroscopic imaging survey, to quantify and compare the gaseous and stellar kinematics of a well-selected, distance-limited sample of up to 20 nearby Seyfert galaxies, and 20 inactive control galaxies with well-matched optical properties. The survey aims to search for dynamical triggers of nuclear activity and address outstanding controversies in optical/infrared imaging surveys.     

Reconstruction of fire spread within wildland fire events in Northern Eurasia from the MODIS active fire product

Russian boreal forests have been reshaped by wildland fire for millennia. While fire is a natural component of boreal ecosystems, it impacts various aspects of the environment and affects human well-being. Often fires occur over large remote areas with limited access, which makes their ground-based observation difficult. A significant progress has been made in mapping burned area from satellite imagery, which provides consistent and fairly unbiased estimates of fire impact on areas of interest at multiple scales. Although the information provided by burned area products is highly important, the spatio-temporal dynamics of individual fire events and their impact are less known. In high northern latitudes of Northern Eurasia, MODIS (Moderate Resolution Imaging Spectroradiometer) makes up to four daily observations from each of the Terra and Aqua satellites providing consistent data on fire development with high temporal frequency. Here we introduce an approach to reconstruct the development of fire events based on active fire detections from MODIS. Fire Spread Reconstruction (FSR) provides a means for characterization of fire occurrence over large territories from remotely sensed data. Individual fire detections are clustered within a GIS environment based on a set of rules determining proximity between fire observations in space and time. FSR determines the number of fire events, their approximate size, duration, and fire spread rate and allows for the analysis of fire occurrence and spread as a function of vegetation, fire season, fire weather and other parameters. FSR clusters were compared to burned scars mapped from Landsat7/ETM+ imagery over Yakutia (Russia). While some smaller burn scars were found to be formed through a continuous burning of a single fire event, large burned areas in Siberia were created by a constellation of fire events incorporating over 100 individual fire clusters. Geographic regions were found to have a stronger influence on the rates of fire activity in the area compared to vegetation zones. In addition, fire spread rates do not directly correlate with the intensity of a given fire season. FSR is also used to identify the points of ignition for individual fire events in spatio-temporal domain for fire danger and fire threat modeling. This approach presents another step towards the more complete characterization of fire events from remotely sensed data.