The E ring in the vicinity of Enceladus: I. Spatial distribution and properties of the ring particles

Saturn's diffuse E ring is the largest ring of the Solar System and extends from about (Saturn radius RS=60,330 km) to at least encompassing the icy moons Mimas, Enceladus, Tethys, Dione, and Rhea. After Cassini's insertion into her saturnian orbit in July 2004, the spacecraft performed a number of equatorial as well as steep traversals through the E ring inside the orbit of the icy moon Dione. Here, we report about dust impact data we obtained during 2 shallow and 6 steep crossings of the orbit of the dominant ring source—the ice moon Enceladus. Based on impact data of grains exceeding 0.9 μm we conclude that Enceladus feeds a torus populated by grains of at least this size along its orbit. The vertical ring structure at agrees well with a Gaussian with a full-width-half-maximum (FWHM) of ∼4200 km. We show that the FWHM at is due to three-body interactions of dust grains ejected by Enceladus' recently discovered ice volcanoes with the moon during their first orbit. We find that particles with initial speeds between 225 and 235 m s⁻¹ relative to the moon's surface dominate the vertical distribution of dust. Particles with initial velocities exceeding the moon's escape speed of 207 m s⁻¹ but slower than 225 m s⁻¹ re-collide with Enceladus and do not contribute to the ring particle population. We find the peak number density to range between 16×¹⁰⁻² m⁻³ and 21×¹⁰⁻² m⁻³ for grains larger 0.9 μm, and 2.1×¹⁰⁻² m⁻³ and 7.6×¹⁰⁻² m⁻³ for grains larger than 1.6 μm. Our data imply that the densest point is displaced outwards by at least with respect of the Enceladus orbit. This finding provides direct evidence for plume particles dragged outwards by the ambient plasma. The differential size distribution for grains >0.9 μm is described best by a power law with slopes between 4 and 5. We also obtained dust data during ring plane crossings in the vicinity of the orbits of Mimas and Tethys. The vertical distribution of grains >0.8 μm at Mimas orbit is also well described by Gaussian with a FWHM of ∼5400 km and displaced southwards by ∼1200 km with respect to the geometrical equator. The vertical distribution of ring particles in the vicinity of Tethys, however, does not match a Gaussian. We use the FWHM values obtained from the vertical crossings to establish a 2-dimensional model for the ring particle distribution which matches our observations during vertical and equatorial traversals through the E ring.     

The planetary fourier spectrometer (PFS) onboard the European Venus Express mission

The planetary fourier spectrometer (PFS) for the Venus Express mission is an infrared spectrometer optimized for atmospheric studies. This instrument has a short wavelength (SW) channel that covers the spectral range from 1700 to 11400 cm⁻¹ (0.9–5.5 μm) and a long wavelength (LW) channel that covers 250–1700 cm⁻¹ (5.5–45 μm). Both channels have a uniform spectral resolution of 1.3 cm⁻¹. The instrument field of view FOV is about 1.6 ° (FWHM) for the short wavelength channel and 2.8 ° for the LW channel which corresponds to a spatial resolution of 7 and 12 km when Venus is observed from an altitude of 250 km. PFS can provide unique data necessary to improve our knowledge not only of the atmospheric properties but also surface properties (temperature) and the surface-atmosphere interaction (volcanic activity).PFS works primarily around the pericentre of the orbit, only occasionally observing Venus from larger distances. Each measurements takes 4.5 s, with a repetition time of 11.5 s. By working roughly 1.5 h around pericentre, a total of 460 measurements per orbit will be acquired plus 60 for calibrations. PFS is able to take measurements at all local times, enabling the retrieval of atmospheric vertical temperature profiles on both the day and the night side.The PFS measures a host of atmospheric and surface phenomena on Venus. These include the:(1) thermal surface flux at several wavelengths near 1 μm, with concurrent constraints on surface temperature and emissivity (indicative of composition); (2) the abundances of several highly-diagnostic trace molecular species; (3) atmospheric temperatures from 55 to 100 km altitude; (4) cloud opacities and cloud-tracked winds in the lower-level cloud layers near 50-km altitudes; (5) cloud top pressures of the uppermost haze/cloud region near 70–80 km altitude; and (6) oxygen airglow near the 100 km level. All of these will be observed repeatedly during the 500-day nominal mission of Venus Express to yield an increased understanding of meteorological, dynamical, photochemical, and thermo-chemical processes in the Venus atmosphere. Additionally, PFS will search for and characterize current volcanic activity through spatial and temporal anomalies in both the surface thermal flux and the abundances of volcanic trace species in the lower atmosphere.Measurement of the 15 μm CO₂ band is very important. Its profile gives, by means of a complex temperature profile retrieval technique, the vertical pressure-temperature relation, basis of the global atmospheric study.PFS is made of four modules called O, E, P and S being, respectively, the interferometer and proximity electronics, the digital control unit, the power supply and the pointing device.     

The impact of a weak south pole on thermal convection in Enceladus’ ice shell

Enceladus exhibits a strong hemispheric dichotomy of tectonism and heat flux, with geologically young, heavily tectonized terrains and a high heat flux in the South Polar Terrain (SPT) and relatively ancient terrains with presumably lower heat fluxes over the rest of the satellite. To understand the convective pattern and its relationship with surface tectonics, we present three-dimensional numerical models of convection in Enceladus’ ice shell including basal heating and tidal heating. Our thermal boundary conditions exhibit no north–south asymmetries, but because the tectonism at the SPT may weaken the ice there, we impose a mechanically weak lithosphere within the SPT. The weakening is parameterized by adopting a reduced viscosity contrast within the SPT. Without such a weak zone, convection (if any) resides in stagnant-lid mode and exhibits no hemispheric dichotomy. In the presence of such an SPT weak zone, however, we find vigorous convection in the ice underneath the SPT, with convective plumes rising close to the surface. In contrast, only stagnant lid convection, or no convection at all, occurs elsewhere over the satellite. Away from the SPT, the heat flux in our models is small (5–10 mW m^?2) and the surface strains are small enough to imply surface ages >10⁹ years. Within the SPT, however, our models yield peak heat fluxes of ～70–200 mW m^?2, implying heat flows integrated across the SPT of up to 5 GW, similar to that inferred from Cassini thermal observations. The surface strains in our models are high enough near the south pole to cause intense tectonism and imply surface ages of ～10⁶–10⁷ years, consistent with age estimates of the SPT.     

New wind measurements in Venus’ lower mesosphere from visible spectroscopy

The dynamics of Venus’ mesosphere (70–110 km) is characterized by the superposition of two different wind regimes: (1) Venus’ retrograde superrotation; (2) a sub-solar to anti-solar (SS–AS) flow pattern, driven by solar EUV heating on the sunlit hemisphere. Here, we report on new ground-based velocity measurements in the lower part of the mesosphere. We took advantage of two essentially symmetric Venus elongations in 2001 and 2002 to perform high-resolution Doppler spectroscopy (R=120,000) in ¹²C¹⁶O₂ visible lines of the 5ν₃ band and in a few solar Fraunhofer lines near 8700 Å. These measurements, mapped over several points on Venus’ illuminated hemisphere, probe the region of cloud tops. More precisely, the solar Fraunhofer lines sample levels a few kilometers below the UV features (i.e. near ∼67 km), while the CO₂ lines probe an altitude higher by about 7 km. The wind field over Venus’ disk is retrieved with an rms uncertainty of 15–25 m s⁻¹ on individual measurements. Kinematical fit to a one- or two-component circulation model indicates the dominance of the zonal retrograde flow with a mean equatorial velocity of ∼75 m s⁻¹, exhibiting very strong day-to-day variations (±65 m s⁻¹). Results are very consistent for the two kinds of lines, suggesting a negligible vertical wind shear over 67–74 km. The SS–AS flow is not detected in single-day observations, but combining the results from all data suggests that this component may invade the lower mesosphere with a ∼40 m s⁻¹ velocity.     

Penetration of nearby supernova dust in the inner solar system

We investigate the method by which nearby supernovae – within a few tens of pc of the solar system – can penetrate the solar system and deposit live radioactivities on earth. The radioactive isotopic signatures that could potentially leave an observable geological imprint are in the form of refractory metals; consequently, it is likely they would arrive in the form of supernova-produced dust grains. Such grains can penetrate into the solar system more easily than the bulk supernova plasma, which gets stalled and deflected near the solar system due to the solar wind plasma pressure. We therefore examine the motion of charged grains as they decouple from the supernova plasma and are influenced by the solar magnetic, radiation, and gravitational fields. We characterize the dust trajectories with analytical approximations which display the roles of grain size, initial velocity, and surface voltage. These results are verified with full numerical simulations for wide ranges of dust properties. We find that supernova dust grains traverse the inner solar system nearly undeflected, if the incoming grain velocity – which we take to be that of the incident supernova remnant – is comparable to the solar wind speeds and much larger than the escape velocity at 1 AU. Consequently, the dust penetration to 1 AU has essentially 100% transmission probability and the dust capture onto the earth should have a geometric cross section. Our results cast in a new light the terrestrial deposition of radioisotopes from nearby supernovae in the geological past. For explosions beyond ～10 pc from earth, dust grains can still deliver supernova ejecta to earth, and thus the amount of supernova material deposited is set by the efficiency of dust condensation and survival in supernovae. Turning the problem around, we use observations of live ⁶⁰Fe in both deep-ocean and lunar samples to infer a conservative lower bound iron condensation efficiency of M_dust,Fe/M_tot,Fe ? 4  × 10^?4 for the supernova which apparently produced these species 2–3 Myr ago.     

Spatial variability,composition and thickness of the seasonal north polar cap of Mars in mid-spring

The planetary fourier spectrometer (PFS) for the Mars express mission (MEX) is an infrared spectrometer operating in the wavelength range from 1.2 to 45 μm by means of two spectral channels, called SWC (short wavelength channel) and LWC (long wavelength channel), covering, respectively, 1.2–5.5 and 5.5–45 μm.The middle-spring Martian north polar cap (Ls∼40°) has been observed by PFS/MEX in illuminated conditions during orbit 452. The SWC spectra are here used to study the cap composition in terms of CO₂ ice, H₂O ice and dust content. Significant spectral variation is noted in the cap interior, and regions of varying CO₂ ice grain sizes, water frost abundance, CO₂ ice cover and dust contamination can be distinguished. In addition, we correlate the infrared spectra with an image acquired during the same orbit by the OMEGA imaging spectrometer and with the altimetry from MOLA data. Many of the spectra variations correlate with heterogeneities noted in the image, although significant spectral variations are not discernible in the visible. The data have been divided into five regions with different latitude ranges and strong similarities in the spectra, and then averaged. Bi-directional reflectance models have been run with the appropriate lighting geometry and used to fit the observed data, allowing for CO₂ ice and H₂O ice grain sizes, dust and H₂O ice contaminations in the form of intimate granular mixtures and spatial mixtures.A wide annulus of dusty water ice surrounds the recessing CO₂ seasonal cap. The inner cap exhibits a layered structure with a thin CO₂ layer with varying concentrations of dark dust, on top of an H₂O ice underneath ground. In the best-fits, the ices beneath the top layer have been considered as spatial mixtures. The results are still very good everywhere in the spectral range, except where the CO₂ ice absorption coefficients are such that even a thin layer is enough to totally absorb the incoming radiation (i.e. the band is saturated). This only happens around 3800 cm⁻¹, inside the strong 2.7-μm CO₂ ice absorption band. The effect of finite snow depth has been investigated through a layered albedo model. The thickness of the CO₂ ice deposits increases with latitude, ranging from 0.5–1 g cm⁻² within region II to 60–80 g cm⁻² within the highest-latitude (up to 84°N) region V.Region I is at the cap edge and extends from 65°N to 72°N latitude. No CO₂ ice is present in this region, which consists of relatively large grains of water ice (20 μm), highly contaminated by dust (0.15 wt%). The adjacent region II is a narrow region [76–79°N] right at the edge of the north residual polar cap. This region is very distinct in the OMEGA image, where it appears to surround the whole residual cap. The CO₂ ice features are barely visible in these spectra, except for the strong saturated 2.7 μm band. It basically consists of a thin layer of 5-mm CO₂ ice on top of an H₂O ice layer with the same composition as region I. A third interesting region III is found all along the shoulder of the residual cap [79–81°N]. It extends over 1.5 km in altitude and over only 2° of latitude and consists of CO₂ ice with a large dust content. It is an admixture of CO₂ ice (3–4 mm), with several tens of ppm by mass of water ice and more than 2 ppt by mass of dust. The surface temperatures have been retrieved from the LWC spectra for each observation. We found an increase in the surface temperature in this region, indicating a spatial mixture of cold CO₂ ice and warmer dust/H₂O ice. Region IV is close to the top of the residual cap [81–84°N]; it is much brighter than region III, with a dust content 10 times lower than the latter. The CO₂ grain size is 3 mm and strong CO₂ ice features are present in the data, indicating a thicker CO₂ ice layer than in region II (1–2 g cm⁻²). The final region V is right at the top of the residual cap (⩾84°N). It is “pure” CO₂ ice (no dust) of 5 mm grain sizes, with 30 ppm by weight of water ice. The CO₂ ice features are very pronounced and the 2.7 μm band is saturated. The optical thickness is close to the semi-infinite limit (30–40 g cm⁻²). Assuming a snowpack density of 0.5 g cm⁻³, we get a minimum thickness of 1–2 cm for the top-layer of regions II and III, 4–10 cm for region IV, and ⩾60–80 cm thickness for region V. These values are in close agreement with several recent results for the south seasonal polar cap.These results should provide new, useful constraints in models of the Martian climate system and volatile cycles.     

Assessing the long-term variability of Venus winds at cloud level from VIRTIS–Venus Express

The Venus Express (VEX) mission has been in orbit to Venus for more than 4 years now. The Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) instrument onboard VEX observes Venus in two channels (visible and infrared) obtaining spectra and multi-wavelength images of the planet that can be used to sample the atmosphere at different altitudes. Day-side images in the ultraviolet range (380 nm) are used to study the dynamics of the upper cloud at 66–72 km while night-side images in the near infrared (1.74 μm) map the opacity of the lower cloud deck at 44–48 km. Here we present a long-term analysis of the global atmospheric dynamics at these levels using a large selection of orbits from the VIRTIS-M dataset covering 860 Earth days that extends our previous work (Sánchez-Lavega, A. et al. [2008]. Geophys. Res. Lett. 35, L13204) and allows studying the variability of the global circulation at the two altitude levels. The atmospheric superrotation is evident with equatorial to mid-latitudes westward velocities of 100 and 60 m s^?1 in the upper and lower cloud layers. These zonal velocities are almost constant in latitude from the equator to 50°S. From 50°S to 90°S the zonal winds at both cloud layers decrease steadily to zero at the pole. Individual cloud tracked winds have errors of 3–10 m s^?1 with a mean of 5 m s^?1 and the standard deviations for a given latitude of our zonal and meridional winds are 9 m s^?1. The zonal winds in the upper cloud change with the local time in a way that can be interpreted in terms of a solar tide. The zonal winds in the lower cloud are stable at mid-latitudes to the tropics and present variability at subpolar latitudes apparently linked to the activity of the South polar vortex. While the upper cloud presents a net meridional motion consistent with the upper branch of a Hadley cell with peak velocity v = 10 m s^?1 at 50°S, the lower cloud meridional motions are less organized with some cloud features moving with intense northwards and southwards motions up to v = ±15 m s^?1 but, on average, with almost null global meridional motions at all latitudes. We also examine the long-term behavior of the winds at these two vertical layers by comparing our extended wind tracked data with results from previous missions.     

Spectral behavior of the symbiotic nova AG Pegasi observed with IUE and HST

Ultraviolet spectra from the International Ultraviolet Explorer (IUE) and from the Hubble Space Telescope (HST) of the symbiotic novae AG Peg during the period 1978–1996 are analyzed. Some spectra showing the modulations of spectral lines at different times are presented. We determined the reddening from the 2200 Å feature, finding that E(B−V) = 0.10 ± 0.02. We studied N IV] at 1486 Å, C IV 1550 Å, and O III] at 1660 Å, produced in the fast wind from the hot white dwarf. The mean wind velocity of the three emission lines is 1300 km s⁻¹ (FWHM). The mean wind mass loss rate is ∼6 × 10⁻⁷ M_⊙ yr⁻¹. The mean temperature is ∼6.5 × 10⁵ K. The mean ultraviolet luminosity is ∼5 × 10³³ erg s⁻¹. The modulations of line fluxes in the emitting region at different times are attributed to the variations of density and temperature of the ejected matter as a result of variations in the rate of mass loss.     

Study of large solar energetic particle events with halo coronal mass ejections and their associated solar flares

The purpose of the present study is to investigate the association of solar energetic particle (SEP) events with halo coronal mass ejections (CME) and with their associated solar flares during the period 1997–2014 (solar cycle 23 and 24). We have found that halo CMEs are more effective in producing SEP events. The occurrence probability and peak fluxes of SEPs strongly depend on the halo CMEs speed (V) as follows. The highest associations, 56% for occurrence probability and 90% for average peak fluxes, are found for the halo CMEs with V> 1400 km s⁻¹ but the lowest associations, 20% for occurrence probability and 5% for average peak fluxes, are found for halo CMEs with speed range 600 ≤ V ≤ 1000 km s⁻¹. We have also examined the relationship between SEP events and halo CME associated solar flares and found that 73% of events are associated with western solar flares while only 27% are with eastern solar flares. For longitudinal study, 0–20° belt is found to be more dominant for the SEP events. The association of SEP events with latitudinal solar flares is also examined in the study. 51% of events are associated with those halo CMEs associated solar flares which occur in the southern hemisphere of the Sun while 49% are with those solar flares that occur in the northern hemisphere of the Sun. Also, 10–20° latitudinal belt is found to be likely associated with the SEP events. Further, 45% of SEP events are associated with M-class solar flares while 44% and 11% are with X and C-class respectively. Maximum number of SEP events are found for the fast halo CME associated X- class solar flares (68%) than M and C- class solar flares.     

Cluster observations of Pc 1–2 waves and associated ion distributions during the October and November 2003 magnetic storms

Unusual wave activity in the Pc 1–2 frequency band (0.1–5 Hz) was observed by the Cluster spacecraft in association with the two large geomagnetic storms of late 2003. During the onset of the Hallowe’en storm on October 29, 2003, intense broadband activity between ∼0.1 and 0.6 Hz appeared at all 4 spacecraft on both sides of the magnetic equator at perigee (near 1400 UT and 08:45 MLT). Power was especially strong and more structured in frequency in the compressional component: a minimum in wave power was observed at 0.38 Hz, corresponding to the oxygen ion cyclotron frequency. Poynting vector calculations indicated that wave power was primarily directed radially inward rather than along the magnetic field. Narrowband purely compressional waves near 0.15 Hz appeared at higher dayside latitudes in the southern hemisphere. CIS ion spectrometer data during this pass revealed that O⁺ was the dominant energetic ion. During the recovery phase of the November storm, on November 22, 2003, predominantly transverse 1.8 Hz waves with peak-to-peak amplitude of 10 nT were observed by all four spacecraft near perigee at L=4.4. During this more typical Pc 1 event, wave power was directed along B, toward the northern ionosphere. An unusually polarized 2.3 Hz emission (with power in the radial and compressional, but not azimuthal directions) was observed at L=5.4–5.9, 10–15° south of the magnetic equator. We infer that this wave event may have been generated on lower L shells and propagated obliquely to Cluster's location. Consistent with other recent observations, anisotropic plasma sheet/ring current proton distributions appeared to be a necessary condition for occurrence of waves during both passes, but was not always a sufficient condition. The transverse waves of November 22 occurred in regions which also contained greatly increased fluxes of cool ions (E<1 keV). On both days, Cluster observed features not previously reported, and we note that the purely compressional nature of the October 29 events was not anticipated in previous theoretical studies. The fact that these unusually polarized waves occurred in association with very intense geomagnetic storms suggests that they are likely to be extremely rare.     

The heat flow during the formation of ribbon terrains on Venus

Ribbons are regularly spaced, between 2 and 6 km, troughs that exist on venusian tesserae, which are mainly located in, and characterize to, venusian crustal plateaus. Independent of the geological or temporal relations with other features, regularly and similarly spaced ribbons on several tesserae strongly suggest a thermal control on the thickness of the deformed layer. This can be used to constraint the heat flow at the time of ribbon formation, which holds important implications for the viability of the hypotheses that address the origin and evolution of crustal plateaus. For a brittle–ductile transition ∼1–3 km deep (as proposed from ribbon spacing), realistic strain rates, and a present-day surface temperature of 740 K, the implied heat flow is very high, 130–780 mW m⁻². If Venus has experienced higher surface temperatures due to climate forcing by massive volcanism, then the heat flow could be greatly reduced. For surface temperatures of 850–900 K the heat flow is 190–560, 60–230 and 20–130 mW m⁻² for brittle–ductile transition depths of 1, 2 and 3 km, respectively. Heat flow values around 80–100 mW m⁻² are reasonable for venusian hotspots, based on terrestrial analogs, but hardly consistent with coldspot settings. High surface temperatures are also required to maintain the crustal solidus deeper than a few kilometers during the formation of ribbon terrains. For the obtained heat flows, a solidus deeper than ∼30 km (the likely mean value for the crustal thickness) is difficult to achieve. This suggests that a substantial proportion of the crust beneath crustal plateaus was emplaced subsequently to the time when ribbon terrains were formed. Alternatively, at that time a magma reservoir inside the crust could have existed.     

The symbiotic system AG Draconis: Soft X-ray bremsstrahlung from the nebulae

The modeling of UV and optical spectra emitted from the symbiotic system AG Dra, adopting collision of the winds, predicts soft X-ray bremsstrahlung from nebulae downstream of the reverse shock with velocities > 150 km s^?1 and intensities comparable to those of the WD black body flux. At outbursts, the envelop of debris, which corresponds to the nebula downstream of the high velocity shocks (700–1000 km s^?1) accompanying the blast wave, absorbs the black body soft X-ray flux from the WD, explains the broad component of the H and He lines, and leads to low optical-UV–X-ray continuum fluxes. The high optical-UV flux observed at the outbursts is explained by bremsstrahlung downstream of the reverse shock between the stars. The depletion of C, N, O, and Mg relative to H indicates that they are trapped into dust grains and/or into diatomic molecules, suggesting that the collision of the wind from the WD with the dusty shells, ejected from the RG with a ～1 year periodicity, leads to the U-band fluctuations during the major bursts.     

Hydrogen density in the dayside venusian exosphere derived from Lyman-α observations by SPICAV on Venus Express

A series of observations of the venusian hydrogen corona made by SPICAV on Venus Express are analyzed to estimate the amount of hydrogen in the exosphere of Venus. These observations were made between November 2006 and July 2007 at altitudes from 1000 km to 8000 km on the dayside. The Lyman-α brightness profiles derived are reproduced by the sum of a cold hydrogen population dominant below ～2000 km and a hot hydrogen population dominant above ～4000 km. The temperature (～300 K) and hydrogen density at 250 km (～10⁵ cm^?3) derived for the cold populations, near noon, are in good agreement with previous observations. Strong dawn–dusk exospheric asymmetry is observed from this set of observations, with a larger exobase density on the dawn side than on the dusk side, consistent with asymmetry previously observed in the venusian thermosphere, but with a lower dawn/dusk contrast. The hot hydrogen density derived is very sensitive to the sky background estimate, but is well constrained near 5000 km. The density of the hot population is reproduced by the exospheric model from Hodges (Hodges, R.R. [1999]. J. Geophys. Res. 104, 8463–8471) in which the hot population is produced by neutral–ions interactions in the thermosphere of Venus.     

Dynamical properties of the Venus mesosphere from the radio-occultation experiment VeRa onboard Venus Express

The dynamics of Venus’ mesosphere (60–100 km altitude) was investigated using data acquired by the radio-occultation experiment VeRa on board Venus Express. VeRa provides vertical profiles of density, temperature and pressure between 40 and 90 km of altitude with a vertical resolution of few hundred meters of both the Northern and Southern hemisphere. Pressure and temperature vertical profiles were used to derive zonal winds by applying an approximation of the Navier–Stokes equation, the cyclostrophic balance, which applies well on slowly rotating planets with fast zonal winds, like Venus and Titan. The main features of the retrieved winds are a midlatitude jet with a maximum speed up to 140 ± 15 m s^?1 which extends between 20°S and 50°S latitude at 70 km altitude and a decrease of wind speed with increasing height above the jet. Cyclostrophic winds show satisfactory agreement with the cloud-tracked winds derived from the Venus Monitoring Camera (VMC/VEx) UV images, although a disagreement is observed at the equator and near the pole due to the breakdown of the cyclostrophic approximation. Knowledge of both temperature and wind fields allowed us to study the stability of the atmosphere with respect to convection and turbulence. The Richardson number Ri was evaluated from zonal field of measured temperatures and thermal winds. The atmosphere is characterised by a low value of Richardson number from ～45 km up to ～60 km altitude at all latitudes that corresponds to the lower and middle cloud layer indicating an almost adiabatic atmosphere. A high value of Richardson number was found in the region of the midlatitude jet indicating a highly stable atmosphere. The necessary condition for barotropic instability was verified: it is satisfied on the poleward side of the midlatitude jet, indicating the possible presence of wave instability.     

An Accuracy Analysis of the SGP4/SDP4 Model

Based on the latest release of the SGP4/SDP4 (Simplified General Perturbation Version 4/ Simplified Deep-space Perturbation Version 4) model, in this paper we have designed an orbit determination program. Through calculations for the 1120 objects with various types and orbital elements selected from the space objects database, we have obtained the accuracies of the orbit determination prediction dealt with various types of space objects by the SGP4/SDP4 model. The results show that the accuracies of the near-earth objects are in the order of magnitude of 100 meters; the averages of the orbit determination accuracies of the semi-synchronous and geosynchronous orbits are, respectively, 0.7 and 1.9 km. The orbit determination accuracies of the elliptical orbit objects are related to their eccentricities. Except for few elliptical orbit objects with e > 0.8, the orbit determination errors of the vast majority of the elliptical orbit objects are all less than 10 km. By using the SGP4/SDP4 model to make 3 days predictions for near-earth objects, 30 days for semi-synchronous orbit objects, 15 days for geosynchronous orbit objects and 1 day for elliptical orbit objects, the errors of prediction generally don’t exceed 40 km.     

Characterizing Io’s Pele,Tvashtar and Pillan plumes: Lessons learned from Hubble

Hubble Space Telescope/Wide Field and Planetary Camera 2 (HST/WFPC2) images of Io obtained between 1995 and 2007 between 0.24 and 0.42 μm led to the detection of the Pele plume in reflected sunlight in 1995 and 1999; imaging of the Pele plume via absorption of jovian light in 1996 and 1999; detection of the Prometheus-type Pillan plume in reflected sunlight in 1997; and detection of the 2007 Pele-type Tvashtar plume eruption in reflected sunlight and via absorption of jovian light. Based on a detailed analysis of these observations we characterize and compare the gas and dust properties of each of the detected plumes. In each case, the brightness of the plumes in reflected sunlight is less at 0.26 μm than at 0.33 μm. Mie scattering analysis of the wavelength dependence of each plume’s reflectance signature suggests that range of particle sizes within the plumes is quite narrow. Assuming a normal distribution of particle sizes, the range of mean particle sizes is ～0.035–0.12 μm for the 1997 Pillan eruption, ～0.05–0.08 μm for the 1999 Pele and 2007 Tvasthar plumes, and ～0.05–0.11 μm for the 1995 Pele plume, and in each case the standard deviation in the particle size distribution is <15%. The Mie analysis also suggests that the 2007 Tvashtar eruption released ～10⁹ g of sulfur dust, the 1999 Pele eruption released ～10⁹ g of SO₂ dust, the 1997 Pillan eruption released ～10¹⁰ g of SO₂ dust, and the 1995 Pele plume may have released ～10¹⁰ g of SO₂ dust. Analysis of the plume absorption signatures recorded in the F255W filter bandpass (0.24–0.28 μm) indicates that the opacity of the 2007 Tvashtar plume was 2× that of the 1996 and 1999 Pele plume eruptions. While the sulfur dust density estimated for the Tvashtar from the reflected sunlight data could have produced 61% of the observed plume opacity, <10% of the 1999 Pele F255W plume opacity could have resulted from the SO₂ dust detected in the eruption. Accounting for the remaining F255W opacity level of the Pele and Tvasthar plumes based on SO₂ and S₂ gas absorption, the SO₂ and S₂ gas density inferred for each plume is almost equivalent corresponding to ～2–6 × 10¹⁶ cm^?2 and 3–5 × 10¹⁵ cm^?2, respectively, producing SO₂ and S₂ gas resurfacing rates ～0.04–0.2 cm yr^?1 and 0.007–0.01 cm yr^?1; and SO₂ and S₂ gas masses ～1–4 × 10¹⁰ g and ～2–3 × 10⁹ g; for a total dust to gas ratio in the plumes ～10^?1–10^?2. The 2007 Tvashtar plume was detected by HST at ～380 ± 40 km in both reflected sunlight and absorbed jovian light; in 1999, the detected Pele plume altitude was 500 km in absorbed jovian light, but in reflected sunlight the detected height was ～2× lower. Thus, for the 1999 Pele plume, similar to the 1979 Voyager Pele plume observations, the most efficient dust reflections occurred in the region closest to the plume vent. The 0.33–0.42 μm brightness of the 1997 Pillan plume was 10–20× greater than the Pele or Tvashtar plumes, exceeding by a factor of 3 the average brightness levels observed within 200 km of 1979 Loki eruption vent. But, the 0.26 μm brightness of the 1997 Pillan plume in reflected sunlight was significantly lower than would be predicted by the dust scattering model. Presuming that the 0.26 μm brightness of the 1997 Pillan plume was attenuated by the eruption plume’s gas component, then an SO₂ gas density ～3–6 × 10¹⁸ cm^?2 is inferred from the data (for S₂/SO₂ ratios ?4%), comparable to the 0.3–2 × 10¹⁸ cm^?2 SO₂ density detected at Loki in 1979 (Pearl, J.C. et al. [1979]. Nature 280, 755; Lellouch et al., 1992), and producing an SO₂ gas mass ～3–8 × 10¹¹ g and an SO₂ resurfacing rate ～8–23 cm yr^?1. These results confirm the connection between high (?10¹⁷ cm^?2) SO₂ gas content and plumes that scatter strongly at nearly blue wavelengths, and it validates the occurrence of high density SO₂ gas eruptions on Io. Noting that the SO₂ gas content inferred from a spectrum of the 2003 Pillan plume was significantly lower ～2 × 10¹⁶ cm^?2 (Jessup, K.L., Spencer, J., Yelle, R. [2007]. Icarus 192, 24–40); and that the Pillan caldera was flooded with fresh SO₂ frost/slush just prior to the 1997 Pillan plume eruption (Geissler, P., McEwen, A., Phillips, C., Keszthelyi, L., Spencer, J. [2004a]. Icarus 169, 29–64; Phillips, C.B. [2000]. Voyager and Galileo SSI Views of Volcanic Resurfacing on Io and the Search for Geologic Activity at Europa. Ph.D. Thesis, Univ. of Ariz., Tucson); we propose that the density of SO₂ gas released by this volcano is directly linked to the local SO₂ frost abundance at the time of eruption.     

Search for relativistic magnetic monopoles with the Baikal neutrino telescope

We present results of a search for relativistic magnetic monopoles with the Baikal neutrino telescope NT200, using data taken between April 1998 and February 2003. No monopole candidates have been found. We set an upper limit 4.6 × 10⁻¹⁷ cm⁻² s⁻¹ sr⁻¹ for the flux of monopoles with β_m = 1. This is a factor of 20 below the Chudakov–Parker bound which is inferred from the very existence of large-scale galactic magnetic fields.     

Direct wind measurements from November 2007 in Venus’ upper atmosphere using ground-based heterodyne spectroscopy of CO2 at 10 μm wavelength

Between November 23 and 28, 2007, the Cologne Tuneable Heterodyne Infrared Spectrometer THIS was installed at the McMath-Pierce Solar Telescope (Kitt Peak, Arizona, USA) to determine zonal wind velocities and to estimate the subsolar-to-antisolar flow. We investigate dynamics in the upper atmosphere of Venus by measuring the Doppler shift of fully-resolved non-LTE CO₂ emission lines at 959.3917 cm^?1 (10.423 μm), which probe a narrow altitude region in Venus’ atmosphere around 110 ± 10 km (～1 μbar). The results show no significant zonal wind velocity at the equator. An increase with latitude up to 43 ± 13 m/s at a latitude of 33°N was observed. This confirms the deduction of a minor influence of Venus superrotation at an altitude of 110 km from previous measurements in May 2007 (Sornig et al., 2008). The specific observing geometry enables estimating the maximum cross terminator velocity of the subsolar-to-antisolar flow at 72 ± 47 m/s.     

Planetesimal growth in turbulent discs before the onset of gravitational instability

It is difficult to imagine a planet formation model that does not at some stage include a gravitationally unstable disc. Initially unstable gas–dust discs may form planets directly, but the high surface density required has motivated the alternative that gravitational instability occurs in a dust sub-layer only after grains have grown large enough by electrostatic sticking. Although such growth up to the instability stage is efficient for laminar discs, previous research concluded that realistic disc turbulence catastrophically increases the settling time, thereby requiring additional processes to facilitate planet formation on the needed time scales. We develop a different model for the influence of turbulence on the collisional velocity of grains and on the scale height of the dust layer and find that the earlier conclusions must be revisited. The model produces a disc-radius dependent time scale to reach a gravitationally unstable phase of planet formation. For a range of dust sticking and disc parameters, we find that for viscosity parameters α ⩽ 10⁻³, this time scale is short enough over a significant range in radii R that turbulence does not catastrophically slow the early phases of planet formation, even in the absence of agglomeration enhancement agents like vortices.     

Spatial and temporal variations in Titan's surface temperatures from Cassini CIRS observations

We report a wide-ranging study of Titan's surface temperatures by analysis of the Moon's outgoing radiance through a spectral window in the thermal infrared at 19 μm (530 cm^?1) characterized by lower atmospheric opacity. We begin by modeling Cassini Composite Infrared Spectrometer (CIRS) far infrared spectra collected in the period 2004–2010, using a radiative transfer forward model combined with a non-linear optimal estimation inversion method. At low-latitudes, we agree with the HASI near-surface temperature of about 94 K at 10°S (Fulchignoni et al., 2005). We find a systematic decrease from the equator toward the poles, hemispherically asymmetric, of ～1 K at 60° south and ～3 K at 60° north, in general agreement with a previous analysis of CIRS data (Jennings et al., 2009), and with Voyager results from the previous northern winter. Subdividing the available database, corresponding to about one Titan season, into 3 consecutive periods, small seasonal changes of up to 2 K at 60°N became noticeable in the results. In addition, clear evidence of diurnal variations of the surface temperatures near the equator are observed for the first time: we find a trend of slowly increasing temperature from the morning to the early afternoon and a faster decrease during the night. The diurnal change is ～1.5 K, in agreement with model predictions for a surface with a thermal inertia between 300 and 600 J m^?2 s^?1/2 K^?1. These results provide important constraints on coupled surface–atmosphere models of Titan's meteorology and atmospheric dynamic.