Comet Shoemaker-Levy 9 impacted Jupiter in July 1994, leaving its stratosphere with several new species, among them water vapor (H2O). With the aid of a photochemical model H2O can be used as a dynamical tracer in the jovian stratosphere. In this paper, we aim at constraining vertical eddy diffusion (Kzz) at the levels where H2O resides. We monitored the H2O disk-averaged emission at 556.936 GHz with the Odin space telescope between 2002 and 2019, covering nearly two decades. We analyzed the data with a combination of 1D photochemical and radiative transfer models to constrain vertical eddy diffusion in the stratosphere of Jupiter. The Odin observations show us that the emission of H2O has an almost linear decrease of about 40% between 2002 and 2019.We can only reproduce our time series if we increase the magnitude of Kzz in the pressure range where H2O diffuses downward from 2002 to 2019, i.e. from ~0.2 mbar to ~5 mbar. However, this modified Kzz is incompatible with hydrocarbon observations. We find that, even if allowance is made for the initially large abundances of H2O and CO at the impact latitudes, the photochemical conversion of H2O to CO2 is not sufficient to explain the progressive decline of the H2O line emission, suggestive of additional loss mechanisms. The Kzz we derived from the Odin observations of H2O can only be viewed as an upper limit in the ~0.2 mbar to ~5 mbar pressure range. The incompatibility between the interpretations made from H2O and hydrocarbon observations probably results from 1D modeling limitations. Meridional variability of H2O, most probably at auroral latitudes, would need to be assessed and compared with that of hydrocarbons to quantify the role of auroral chemistry in the temporal evolution of the H2O abundance since the SL9 impacts. Modeling the temporal evolution of SL9 species with a 2D model would be the next natural step.

Io's atmosphere is predominately SO2 sustained by a combination of volcanic outgassing and sublimation. The loss from the atmosphere is the main mass source for Jupiter's large magnetosphere. Previous studies attributed various transient phenomena in Io's environment and Jupiter's magnetosphere to a sudden change in the mass loss from the atmosphere supposedly triggered by a change in volcanic activity. Since the gas in volcanic plumes does not escape directly, such causal correlation would require a transient volcano-induced change in atmospheric abundance, which has never been observed so far. Here we report four observations of atmospheric SO2 and NaCl obtained with the IRAM NOEMA interferometer. These observations are compared to measurements of volcanic hot spots and Io's neutral and plasma environment. We find a stable NaCl column density in Io's atmosphere. The SO2 column density derived for December 2016 is about 30% lower compared to the period of March to April 2017. This increase in SO2 from December 2016 to March 2017 might be related to increasing volcanic activity observed at several sites in spring 2017, but the stability of the volcanic trace gas NaCl and resulting decrease in NaCl/SO2 ratio do not support this interpretation. Observed dimmings in both the sulfur ion torus and Na neutral cloud suggest rather a decrease in mass loading in the period of increasing SO2 abundance. The dimming Na brightness and stable atmospheric NaCl furthermore dispute an earlier suggested positive correlation of the sodium cloud and the hot spot activity at Loki Patara, which considerably increased in this period. The environment of Io overall appears to be in a quiescent state, preventing further conclusions. Only Jupiter's aurora morphology underwent several short-term changes, which are apparently unrelated to Io's quiescent environment or the relatively stable atmosphere.

The goal of this work is to determine the physical characteristics of resonant, detached and scattered disk objects in the transneptunian region, observed mainly in the framework of the "TNOs are Cool!" Herschel Open Time Key Program. Based on thermal emission measurements with the Herschel/PACS and Spitzer/MIPS instruments we determine size, albedo, and surface thermal properties for 23 objects using radiometric modelling techniques. This is the first analysis in which the physical properties of objects in the outer resonances are determined for a notable sample. In addition to the results for individual objects, we have compared these characteristics with the bulk properties of other populations of the transneptunian region. The newly analyzed objects show a large variety of beaming factors, indicating a diversity of thermal properties, and in general, they follow the albedo-colour clustering identified earlier for Kuiper belt objects and Centaurs, further strengthening the evidence for a compositional discontinuity in the young Solar System.

Saturn's moon Iapetus, which is in synchronous rotation, is covered by an optically dark material mainly on its leading side, while its trailing side is significantly brighter. Because longer wavelengths probe deeper into the subsurface, observing both sides at a variety of wavelengths brings to light possible changes in thermal, compositional, and physical properties with depth. We have observed Iapetus's leading and trailing hemispheres at 1.2 and 2.0 mm, using the NIKA2 camera mounted on the IRAM 30-m telescope, and compared our observations to others performed at mm to cm wavelengths. We calibrate our observations on Titan, which is simultaneously observed within the field of view. Due to the proximity of Saturn, it is sometimes difficult to separate Iapetus's and Titan's flux from that of Saturn, detected in the telescope's side lobes. Preliminary results show that the trailing hemisphere brightness temperatures at the two wavelengths are equal within error bars, unlike the prediction made by Ries (2012). On the leading side, we report a steep spectral slope of increasing brightness temperature (by 10 K) from 1.2 to 2.0 mm, which may indicate rapidly varying emissivities within the top few centimeters of the surface. Comparison to a diffuse scattering model and a thermal model will be necessary to further constrain the thermophysical properties of the subsurface of Iapetus's two faces.

From CCD observations carried out with different telescopes, we present short-term photometric measurements of the large trans-Neptunian object Varuna in 10 epochs, spanning around 19 years. We observe that the amplitude of the rotational light-curve has changed considerably during this period of time from 0.41 to 0.55 mag. In order to explain this variation, we constructed a model in which Varuna has a simple triaxial shape, assuming that the main effect comes from the change of the aspect angle as seen from Earth, due to Varuna's orbital motion in the 19-year time span. The best fits to the data correspond to a family of solutions with axial ratios b/a between 0.56 and 0.60. This constrains the pole orientation in two different ranges of solutions presented here as maps. Apart from the remarkable variation of the amplitude, we have detected changes in the overall shape of the rotational light-curve over shorter time scales. After the analysis of the periodogram of the residuals to a 6.343572 h double-peaked rotational light-curve fit, we find a clear additional periodicity. We propose that these changes in the rotational light-curve shape are due to a large and close-in satellite whose rotation induces the additional periodicity. The peak-to-valley amplitude of this oscillation is in the order of 0.04 mag. We estimate that the satellite orbits Varuna with a period of 11.9819 h (or 23.9638 h), assuming that the satellite is tidally locked, at a distance of ~ 1300 km (or ~ 2000 km) from Varuna, outside the Roche limit.

Extrasolar satellites are generally too small to be detected by nominal searches. By analogy to the most active body in the Solar System, Io, we describe how sodium (Na I) and potassium (K I) $\textit{gas}$ could be a signature of the geological activity venting from an otherwise hidden exo-Io. Analyzing $\sim$ a dozen close-in gas giants hosting robust alkaline detections, we show that an Io-sized satellite can be stable against orbital decay below a planetary tidal $\mathcal{Q}_p \lesssim 10^{11}$. This tidal energy is focused into the satellite driving a $\sim 10^{5 \pm 2}$ higher mass loss rate than Io's supply to Jupiter's Na exosphere, based on simple atmospheric loss estimates. The remarkable consequence is that several exo-Io column densities are on average $\textit{more than sufficient}$ to provide the $\sim$ 10$^{10 \pm 1}$ Na cm$^{-2}$ required by the equivalent width of exoplanet transmission spectra. Furthermore, the benchmark observations of both Jupiter's extended ($\sim 1000$ R$_J$) Na exosphere and Jupiter's atmosphere in transmission spectroscopy yield similar Na column densities that are purely exogenic in nature. As a proof of concept, we fit the "high-altitude" Na at WASP 49-b with an ionization-limited cloud similar to the observed Na profile about Io. Moving forward, we strongly encourage time-dependent ingress and egress monitoring along with spectroscopic searches for other volcanic volatiles.

Context. The origin of water in the stratospheres of Giant Planets has been an outstanding question ever since its first detection by ISO some 20 years ago. Water can originate from interplanetary dust particles, icy rings and satellites and large comet impacts. Analysis of Herschel Space Observatory observations have proven that the bulk of Jupiter's stratospheric water was delivered by the Shoemaker-Levy 9 impacts in 1994. In 2006, the Cassini mission detected water plumes at the South Pole of Enceladus, placing the moon as a serious candidate for Saturn's stratospheric water. Further evidence was found in 2011, when Herschel demonstrated the presence of a water torus at the orbital distance of Enceladus, fed by the moon's plumes. Finally, water falling from the rings onto Saturn's uppermost atmospheric layers at low latitudes was detected during the final orbits of Cassini's end-of-mission plunge into the atmosphere. Aims. In this paper, we use Herschel mapping observations of water in Saturn's stratosphere to identify its source. Methods. Several empirical models are tested against the Herschel-HIFI and -PACS observations, which were collected on December 30, 2010, and January 2nd, 2011 (respectively). Results. We demonstrate that Saturn's stratospheric water is not uniformly mixed as a function of latitude, but peaking at the equator and decreasing poleward with a Gaussian distribution. We obtain our best fit with an equatorial mole fraction 1.1 ppb and a half-width at half-maximum of 25{\deg}, when accounting for a temperature increase in the two warm stratospheric vortices produced by Saturn's Great Storm of 2010-2011. Conclusions. This work demonstrates that Enceladus is the main source of Saturn's stratospheric water.

From 2004 to 2017, the Cassini spacecraft orbited Saturn, completing 127 close flybys of its largest moon, Titan. Cassini's Composite Infrared Spectrometer (CIRS), one of 12 instruments carried on board, profiled Titan in the thermal infrared (7-1000 microns) throughout the entire 13-year mission. CIRS observed on both targeted encounters (flybys) and more distant opportunities, collecting 8.4 million spectra from 837 individual Titan observations over 3633 hours. Observations of multiple types were made throughout the mission, building up a vast mosaic picture of Titan's atmospheric state across spatial and temporal domains. This paper provides a guide to these observations, describing each type and chronicling its occurrences and global-seasonal coverage. The purpose is to provide a resource for future users of the CIRS data set, as well as those seeking to put existing CIRS publications into the overall context of the mission, and to facilitate future inter-comparison of CIRS results with those of other Cassini instruments, and ground-based observations.

The thermal emission of transneptunian objects (TNO) and Centaurs has been observed at mid- and far-infrared wavelengths - with the biggest contributions coming from the Spitzer and Herschel space observatories-, and the brightest ones also at sub-millimeter and millimeter wavelengths. These measurements allowed to determine the sizes and albedos for almost 180 objects, and densities for about 25 multiple systems. The derived very low thermal inertias show evidence for a decrease at large heliocentric distances and for high-albedo objects, which indicates porous and low-conductivity surfaces. The radio emissivity was found to be low ($\epsilon_r$=0.70$\pm$0.13) with possible spectral variations in a few cases. The general increase of density with object size points to different formation locations or times. The mean albedos increase from about 5-6% (Centaurs, Scattered-Disk Objects) to 15% for the Detached objects, with distinct cumulative albedo distributions for hot and cold classicals. The color-albedo separation in our sample is evidence for a compositional discontinuity in the young Solar System. The median albedo of the sample (excluding dwarf planets and the Haumea family) is 0.08, the albedo of Haumea family members is close to 0.5, best explained by the presence of water ice. The existing thermal measurements remain a treasure trove at times where the far-infrared regime is observationally not accessible.

A group of trans-Neptunian objects (TNO) are dynamically related to the dwarf planet 136108 Haumea. Ten of them show strong indications of water ice on their surfaces, are assumed to have resulted from a collision, and are accepted as the only known TNO collisional family. Nineteen other dynamically similar objects lack water ice absorptions and are hypothesized to be dynamical interlopers. We have made observations to determine sizes and geometric albedos of six of the accepted Haumea family members and one dynamical interloper. Ten other dynamical interlopers have been measured by previous works. We compare the individual and statistical properties of the family members and interlopers, examining the size and albedo distributions of both groups. We also examine implications for the total mass of the family and their ejection velocities. We use far-infrared space-based telescopes to observe the target TNOs near their thermal peak and combine these data with optical magnitudes to derive sizes and albedos using radiometric techniques. We determine the power-law slope of ejection velocity as a function of effective diameter. The detected Haumea family members have a diversity of geometric albedos $\sim$ 0.3-0.8, which are higher than geometric albedos of dynamically similar objects without water ice. The median geometric albedo for accepted family members is $p_V=0.48_{-0.18}^{+0.28}$, compared to 0.08$_{-0.05}^{+0.07}$ for the dynamical interlopers. In the size range $D=175-300$ km, the slope of the cumulative size distribution is $q$=3.2$_{-0.4}^{+0.7}$ for accepted family members, steeper than the $q$=2.0$\pm$0.6 slope for the dynamical interlopers with D$< $500 km. The total mass of Haumea's moons and family members is 2.4% of Haumea's mass. The ejection velocities required to emplace them on their current orbits show a dependence on diameter, with a power-law slope of 0.21-0.50.

Winds in Titan's lower and middle atmosphere have been determined by a variety of techniques, including direct measurements from the Huygens Probe over 0-150 km, Doppler shifts of molecular spectral lines in the optical, thermal infrared and mm ranges, probing altogether the ~100-450 km altitude range, and inferences from thermal field over 10 mbar - 10 -3 mbar (i.e. ~100-500 km) and from central flashes in stellar occultation curves. These measurements predominantly indicated strong prograde winds, reaching maximum speeds of ~150-200 m/s in the upper stratosphere, with important latitudinal and seasonal variations. However, these observations provided incomplete atmospheric sounding; in particular, the wind regime in Titan's upper mesosphere and thermosphere (500- 1200 km) has remained unconstrained so far. Here we report direct wind measurements based on Doppler shifts of six molecular species observed with ALMA. We show that unlike expectations, strong prograde winds extend up to the thermosphere, with the circulation progressively turning into an equatorial jet regime as altitude increases, reaching ~340 m/s at 1000 km. We suggest that these winds may represent the dynamical response of forcing by waves launched at upper stratospheric/mesospheric levels and/or magnetospheric-ionospheric interaction. We also demonstrate that the HNC distribution is restricted to Titan's thermosphere above ~870 km altitude.

Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere.

A recent occultation measurement of the dwarf planet Haumea (Ortiz et al. 2017) revealed an elongated shape with the longest axis comparable to Pluto's mean diameter. The chords also indicate a ring around Haumea's equatorial plane, where its largest moon, Hi'iaka, is also located. The Haumea occultation size estimate (equivalent diameter 1595 km) is larger than previous radiometric solutions (in the range between 1150 and 1350 km), which lowers the object's density to about 1.8 gcm-3, a value closer to the densities of other large TNOs. We present unpublished and also reprocessed Herschel and Spitzer MIR and FIR measurements. We compare 100 and 160um thermal lightcurve amplitudes - originating from Haumea itself - with models of the total measured system fluxes (ring, satellite, Haumea) from 24-350um. The combination with results derived from the occultation measurements show that Haumea's crystalline water ice surface must have a thermal inertia of ~5 SIu. We also have indications that the satellites (at least Hi'iaka) must have high geometric albedos >=0.5, otherwise the derived thermal amplitude would be inconsistent with the total measured system fluxes at 24, 70, 100, 160, 250, and 350um. The high albedos imply sizes of about 300 and 150 km for Hi'iaka and Namaka, respectively, indicating unexpectedly high densities >1.0 gcm-3 for TNOs this small, and the assumed collisional formation from Haumea's icy crust. We also estimated the thermal emission of the ring for the time period 1980-2030, showing that the contribution during the Spitzer and Herschel epochs was small, but not negligible. Due to the progressive opening of the ring plane, the ring emission will be increasing in the next decade when JWST is operational. In the MIRI 25.5um band it will also be possible to obtain a very high-quality thermal lightcurve to test the derived Haumea properties.

We report on results of an observing campaign to support the Juno mission. At the beginning of 2016, using TEXES (Texas Echelon cross-dispersed Echelle Spectrograph), mounted on the NASA Infrared Telescope Facility (IRTF), we obtained data cubes of Jupiter in the 1930--1943 cm$^{-1}$ spectral ranges (around 5 $\mu$m), which probe the atmosphere in the 1--4 bar region, with a spectral resolution of $\approx$ 0.15 cm$^{-1}$ and an angular resolution of $\approx$ 1.4". This dataset is analysed by a code that combines a line-by-line radiative transfer model with a non-linear optimal estimation inversion method. The inversion retrieves the vertical abundance profiles of NH$_3$ - which is the main contributor at these wavelengths - with a maximum sensitivity at $\approx$ 1--3 bar, as well as the cloud transmittance. This retrieval is performed on over more than one thousand pixels of our data cubes, producing maps of the disk, where all the major belts are visible. We present our retrieved NH$_3$ abundance maps which can be compared with the distribution observed by Juno's MWR (Bolton et al., 2017; Li et al., 2017) in the 2 bar region and discuss their significance for the understanding of Jupiter's atmospheric dynamics. We are able to show important latitudinal variations - such as in the North Equatorial Belt (NEB), where the NH$_3$ abundance is observed to drop down to 60 ppmv at 2 bar - as well as longitudinal variability. In the zones, we find the NH$_3$ abundance to increase with depth, from 100 $\pm$ 15 ppmv at 1 bar to 500 $\pm$ 30 ppmv at 3 bar. We also display the cloud transmittance--NH$_3$ abundance relationship, and find different behaviour for the NEB, the other belts and the zones. Using a simple cloud model, we are able to fit this relationship, at least in the NEB, including either NH$_3$-ice or NH$_4$SH particles with sizes between 10 and 100 $\mu$m.

The sensitivity of ALMA makes it possible to detect thermal mm/submm emission from small/distant Solar System bodies at the sub-mJy level. Measured fluxes are primarily sensitive to the objects' diameters, but deriving precise sizes is somewhat hampered by the uncertain effective emissivity at these wavelengths. Following Brown and Butler (2017) who presented ALMA data for four binary TNOs, we report ALMA 1.29 mm measurements of four Centaurs (2002 GZ$_{32}$, Bienor, Chiron, Chariklo) and two TNOs (Huya and Makemake), sampling a range of size, albedo and composition. These thermal fluxes are combined with mid/far-infrared fluxes to derive the relative emissivity at radio (mm/submm) wavelengths, using NEATM and thermophysical models. We reassess earlier thermal measurements of these and other objects -- including Pluto/Charon and Varuna -- exploring effects due to non-spherical shape and varying apparent pole orientation, and show that those can be key for reconciling previous diameter determinations and correctly estimating the spectral emissivities. We also evaluate the possible contribution to thermal fluxes of established (Chariklo) or claimed (Chiron) ring systems. As a general conclusion, all the objects, except Makemake, have radio emissivities significantly lower than unity. Although the emissivity values show diversity, we do not find any significant trend with physical parameters such as diameter, composition, beaming factor, albedo, or color, but we suggest that the emissivity could be correlated with grain size. The mean relative radio emissivity is found to be 0.70$\pm$0.13, a value that we recommend for the analysis of further mm/submm data.

We use data from five stellar occultations observed between 2013 and 2016 to constrain Chariklo's size and shape, and the ring reflectivity. We consider four possible models for Chariklo (sphere, Maclaurin spheroid, tri-axial ellipsoid and Jacobi ellipsoid) and we use a Bayesian approach to estimate the corresponding parameters. The spherical model has a radius $R=129\pm3$ km. The Maclaurin model has equatorial and polar radii $a=b=143^{+3}_{-6}$ km and $c=96^{+14}_{-4}$ km, respectively, with density $970^{+300}_{-180}$ kg m$^{-3}$. The ellipsoidal model has semiaxes $a=148^{+6}_{-4}$ km, $b=132^{+6}_{-5}$ km and $c=102^{+10}_{-8}$ km. Finally, the Jacobi model has semiaxes $a$=157$\pm$4 km, $b$=139$\pm$ 4 km and $c$=86$\pm$1 km, and density $796^{+2}_{-4}$ kg m$^{-3}$ . Depending on the model, we obtain topographic features of 6-11 km, typical of Saturn icy satellites with similar size and density. We constrain Chariklo's geometric albedo between 3.1\% (sphere) and 4.9\% (ellipsoid), while the ring $I/F$ reflectivity is less constrained between 0.6\% (Jacobi) and 8.9\% (sphere). The ellipsoid model explains both the optical light curve and the long-term photometry variation of the system, giving a plausible value for the geometric albedo of the ring particles of $10-15\%$. The derived Chariklo's mass of 6-8$\times10^{18}$ kg places the rings close to the 3:1 resonance between the ring mean motion and Chariklo's rotation period.

The ice giants Uranus and Neptune are the least understood class of planets in our solar system but the most frequently observed type of exoplanets. Presumed to have a small rocky core, a deep interior comprising ~70% heavy elements surrounded by a more dilute outer envelope of H2 and He, Uranus and Neptune are fundamentally different from the better-explored gas giants Jupiter and Saturn. Because of the lack of dedicated exploration missions, our knowledge of the composition and atmospheric processes of these distant worlds is primarily derived from remote sensing from Earth-based observatories and space telescopes. As a result, Uranus's and Neptune's physical and atmospheric properties remain poorly constrained and their roles in the evolution of the Solar System not well understood. Exploration of an ice giant system is therefore a high-priority science objective as these systems (including the magnetosphere, satellites, rings, atmosphere, and interior) challenge our understanding of planetary formation and evolution. Here we describe the main scientific goals to be addressed by a future in situ exploration of an ice giant. An atmospheric entry probe targeting the 10-bar level, about 5 scale heights beneath the tropopause, would yield insight into two broad themes: i) the formation history of the ice giants and, in a broader extent, that of the Solar System, and ii) the processes at play in planetary atmospheres. The probe would descend under parachute to measure composition, structure, and dynamics, with data returned to Earth using a Carrier Relay Spacecraft as a relay station. In addition, possible mission concepts and partnerships are presented, and a strawman ice-giant probe payload is described. An ice-giant atmospheric probe could represent a significant ESA contribution to a future NASA ice-giant flagship mission.

Time series observations of the dwarf planet Haumea and the Plutinos 2003VS2 and 2003AZ84 with Herschel/PACS are presented in this work. Thermal emission of these trans-Neptunian objects were acquired as part of the TNOs are Cool Herschel Space Observatory key programme. We search for the thermal light curves at 100 and 160um of Haumea and 2003AZ84, and at 70 and 160um for 2003VS2 by means of photometric analysis of the PACS data. The goal of this work is to use these thermal light curves to obtain physical and thermophysical properties of these icy Solar System bodies. Haumea's thermal light curve is clearly detected at 100 and 160um. The effect of the reported dark spot is apparent at 100um. Different thermophysical models were applied to these light curves, varying the thermophysical properties of the surface within and outside the spot. Although no model gives a perfect fit to the thermal observations, results imply an extremely low thermal inertia (< 0.5 MKS) and a high phase integral (> 0.73) for Haumea's surface. We note that the dark spot region appears to be only weakly different from the rest of the object, with modest changes in thermal inertia and/or phase integral. The thermal light curve of 2003VS2 is not firmly detected at 70 and at 160um but a thermal inertia of 2+/-0.5 MKS can be derived from these data. The thermal light curve of 2003AZ84 is not firmly detected at 100um. We apply a thermophysical model to the mean thermal fluxes and to all the Herschel/PACS and Spitzer/MIPS thermal data of 2003AZ84, obtaining a close to pole-on orientation as the most likely for this TNO. For the three TNOs, the thermal inertias derived from light curve analyses or from the thermophysical analysis of the mean thermal fluxes confirm the generally small or very small surface thermal inertias of the TNO population, which is consistent with a statistical mean value of 2.5+/-0.5 MKS.

Measurements of flux density are described for five planets, Mars, Jupiter, Saturn, Uranus, and Neptune, across the six Planck High Frequency Instrument frequency bands (100-857 GHz) and these are then compared with models and existing data. In our analysis, we have also included estimates of the brightness of Jupiter and Saturn at the three frequencies of the Planck Low Frequency Instrument (30, 44, and 70 GHz). The results provide constraints on the intrinsic brightness and the brightness time-variability of these planets. The majority of the planet flux density estimates are limited by systematic errors, but still yield better than 1% measurements in many cases. Applying data from Planck HFI, the Wilkinson Microwave Anisotropy Probe (WMAP), and the Atacama Cosmology Telescope (ACT) to a model that incorporates contributions from Saturn's rings to the planet's total flux density suggests a best fit value for the spectral index of Saturn's ring system of $\beta _\mathrm{ring} = 2.30\pm0.03$ over the 30-1000 GHz frequency range. The average ratio between the Planck-HFI measurements and the adopted model predictions for all five planets (excluding Jupiter observations for 353 GHz) is 0.997, 0.997, 1.018, and 1.032 for 100, 143, 217, and 353 GHz, respectively. Model predictions for planet thermodynamic temperatures are therefore consistent with the absolute calibration of Planck-HFI detectors at about the three-percent-level. We compare our measurements with published results from recent cosmic microwave background experiments. In particular, we observe that the flux densities measured by Planck HFI and WMAP agree to within 2%. These results allow experiments operating in the mm-wavelength range to cross-calibrate against Planck and improve models of radiative transport used in planetary science.

Observations of the Pluto-Charon system, acquired with the ALMA interferometer on June 12-13, 2015, have yielded a detection of the CO(3-2) and HCN(4-3) rotational transitions from Pluto, providing a strong confirmation of the presence of CO, and the first observation of HCN, in Pluto's atmosphere. The CO and HCN lines probe Pluto's atmosphere up to ~450 km and ~900 km altitude, respectively. The CO detection yields (i) a much improved determination of the CO mole fraction, as 515+/-40 ppm for a 12 ubar surface pressure (ii) clear evidence for a well-marked temperature decrease (i.e., mesosphere) above the 30-50 km stratopause and a best-determined temperature of 70+/-2 K at 300 km, in agreement with recent inferences from New Horizons / Alice solar occultation data. The HCN line shape implies a high abundance of this species in the upper atmosphere, with a mole fraction >1.5x10-5 above 450 km and a value of 4x10-5 near 800 km. The large HCN abundance and the cold upper atmosphere imply supersaturation of HCN to a degree (7-8 orders of magnitude) hitherto unseen in planetary atmospheres, probably due to the slow kinetics of condensation at the low pressure and temperature conditions of Pluto's upper atmosphere. HCN is also present in the bottom ~100 km of the atmosphere, with a 10-8 - 10-7 mole fraction; this implies either HCN saturation or undersaturation there, depending on the precise stratopause temperature. The HCN column is (1.6+/-0.4)x10^14 cm-2, suggesting a surface-referred net production rate of ~2x10^7 cm-2s-1. Although HCN rotational line cooling affects Pluto's atmosphere heat budget, the amounts determined in this study are insufficient to explain the well-marked mesosphere and upper atmosphere's ~70 K temperature. We finally report an upper limit on the HC3N column density (< 2x10^13 cm-2) and on the HC15N / HC14N ratio (< 1/125).

We compare the absolute gain photometric calibration of the Planck/HFI and Herschel/SPIRE instruments on diffuse emission. The absolute calibration of HFI and SPIRE each relies on planet flux measurements and comparison with theoretical far-infrared emission models of planetary atmospheres. We measure the photometric cross calibration between the instruments at two overlapping bands, 545 GHz / 500 $\mu$m and 857 GHz / 350 $\mu$m. The SPIRE maps used have been processed in the Herschel Interactive Processing Environment (Version 12) and the HFI data are from the 2015 Public Data Release 2. For our study we used 15 large fields observed with SPIRE, which cover a total of about 120 deg^2. We have selected these fields carefully to provide high signal-to-noise ratio, avoid residual systematics in the SPIRE maps, and span a wide range of surface brightness. The HFI maps are bandpass-corrected to match the emission observed by the SPIRE bandpasses. The SPIRE maps are convolved to match the HFI beam and put on a common pixel grid. We measure the cross-calibration relative gain between the instruments using two methods in each field, pixel-to-pixel correlation and angular power spectrum measurements. The SPIRE / HFI relative gains are 1.047 ($\pm$ 0.0069) and 1.003 ($\pm$ 0.0080) at 545 and 857 GHz, respectively, indicating very good agreement between the instruments. These relative gains deviate from unity by much less than the uncertainty of the absolute extended emission calibration, which is about 6.4% and 9.5% for HFI and SPIRE, respectively, but the deviations are comparable to the values 1.4% and 5.5% for HFI and SPIRE if the uncertainty from models of the common calibrator can be discounted. Of the 5.5% uncertainty for SPIRE, 4% arises from the uncertainty of the effective beam solid angle, which impacts the adopted SPIRE point source to extended source unit conversion factor (Abridged)

We present results from a multi-chord Pluto stellar occultation observed on 29 June 2015 from New Zealand and Australia. This occurred only two weeks before the NASA New Horizons flyby of the Pluto system and serves as a useful comparison between ground-based and space results. We find that Pluto's atmosphere is still expanding, with a significant pressure increase of 5$\pm$2\% since 2013 and a factor of almost three since 1988. This trend rules out, as of today, an atmospheric collapse associated with Pluto's recession from the Sun. A central flash, a rare occurrence, was observed from several sites in New Zealand. The flash shape and amplitude are compatible with a spherical and transparent atmospheric layer of roughly 3~km in thickness whose base lies at about 4~km above Pluto's surface, and where an average thermal gradient of about 5 K~km$^{-1}$ prevails. We discuss the possibility that small departures between the observed and modeled flash are caused by local topographic features (mountains) along Pluto's limb that block the stellar light. Finally, using two possible temperature profiles, and extrapolating our pressure profile from our deepest accessible level down to the surface, we obtain a possible range of 11.9-13.7~$\mu$bar for the surface pressure.

Thermal observations of the Pluto-Charon system acquired by the Herschel Space Observatory in February 2012 are presented. They consist of photometric measurements with the PACS and SPIRE instruments (nine visits to the Pluto system each), covering six wavelengths from 70 to 500 $\mu$m altogether. The thermal light curve of Pluto-Charon is observed in all filters, albeit more marginally at 160 and especially 500 $\mu$m. Putting these data into the context of older ISO, Spitzer and ground-based observations indicates that the brightness temperature (T$_B$) of the system (rescaled to a common heliocentric distance) drastically decreases with increasing wavelength, from $\sim$53 K at 20 $\mu$m to $\sim$35 K at 500 $\mu$m, and perhaps ever less at longer wavelengths. Considering a variety of diurnal and/or seasonal thermophysical models, we show that T$_B$ values of 35 K are lower than any expected temperature for the dayside surface or subsurface of Pluto and Charon, implying a low surface emissivity. Based on multiterrain modeling, we infer a spectral emissivity that decreases steadily from 1 at 20-25 $\mu$m to $\sim$0.7 at 500~$\mu$m. This kind of behavior is usually not observed in asteroids (when proper allowance is made for subsurface sounding), but is found in several icy surfaces of the solar system. We tentatively identify that a combination of a strong dielectric constant and a considerable surface material transparency (typical penetration depth $\sim$1 cm) is responsible for the effect. Our results have implications for the interpretation of the temperature measurements by REX/New Horizons at 4.2 cm wavelength.

In this paper we investigate the opportunities provided by the James Webb Space Telescope (JWST) for significant scientific advances in the study of solar system bodies and rings using stellar occultations. The strengths and weaknesses of the stellar occultation technique are evaluated in light of JWST's unique capabilities. We identify several possible JWST occultation events by minor bodies and rings, and evaluate their potential scientific value. These predictions depend critically on accurate a priori knowledge of the orbit of JWST near the Sun-Earth Lagrange-point 2 (L2). We also explore the possibility of serendipitous stellar occultations by very small minor bodies as a by-product of other JWST observing programs. Finally, to optimize the potential scientific return of stellar occultation observations, we identify several characteristics of JWST's orbit and instrumentation that should be taken into account during JWST's development.

The James Webb Space Telescope (JWST), scheduled for launch in 2018, is the successor to the Hubble Space Telescope (HST) but with a significantly larger aperture (6.5 m) and advanced instrumentation focusing on infrared science (0.6-28.0 $\mu$m ). In this paper we examine the potential for scientific investigation of Titan using JWST, primarily with three of the four instruments: NIRSpec, NIRCam and MIRI, noting that science with NIRISS will be complementary. Five core scientific themes are identified: (i) surface (ii) tropospheric clouds (iii) tropospheric gases (iv) stratospheric composition and (v) stratospheric hazes. We discuss each theme in depth, including the scientific purpose, capabilities and limitations of the instrument suite, and suggested observing schemes. We pay particular attention to saturation, which is a problem for all three instruments, but may be alleviated for NIRCam through use of selecting small sub-arrays of the detectors - sufficient to encompass Titan, but with significantly faster read-out times. We find that JWST has very significant potential for advancing Titan science, with a spectral resolution exceeding the Cassini instrument suite at near-infrared wavelengths, and a spatial resolution exceeding HST at the same wavelengths. In particular, JWST will be valuable for time-domain monitoring of Titan, given a five to ten year expected lifetime for the observatory, for example monitoring the seasonal appearance of clouds. JWST observations in the post-Cassini period will complement those of other large facilities such as HST, ALMA, SOFIA and next-generation ground-based telescopes (TMT, GMT, EELT).

In this document, we summarize the main capabilities of the James Webb Space Telescope (JWST) for performing observations of Mars. The distinctive vantage point of JWST at the Sun-Earth Lagrange point (L2) will allow sampling the full observable disk, permitting the study of short-term phenomena, diurnal processes (across the East-West axis) and latitudinal processes between the hemispheres (including seasonal effects) with excellent spatial resolutions (0.07 arcsec at 2 {\mu}m). Spectroscopic observations will be achievable in the 0.7-5 {\mu}m spectral region with NIRSpec at a maximum resolving power of 2700, and with 8000 in the 1-1.25 {\mu}m range. Imaging will be attainable with NIRCam at 4.3 {\mu}m and with two narrow filters near 2 {\mu}m, while the nightside will be accessible with several filters in the 0.5 to 2 {\mu}m. Such a powerful suite of instruments will be a major asset for the exploration and characterization of Mars. Some science cases include the mapping of the water D/H ratio, investigations of the Martian mesosphere via the characterization of the non-LTE CO$_2$ emission at 4.3 {\mu}m, studies of chemical transport via observations of the O$_2$ nightglow at 1.27 {\mu}m, high cadence mapping of the variability dust and water ice clouds, and sensitive searches for trace species and hydrated features on the Martian surface. In-flight characterization of the instruments may allow for additional science opportunities.

We present THERMAP, a mid-infrared (8-16 {\mu}m) spectro-imager for space missions to small bodies in the inner solar system, developed in the framework of the MarcoPolo-R asteroid sample return mission. THERMAP is very well suited to characterize the surface thermal environment of a NEO and to map its surface composition. The instrument has two channels, one for imaging and one for spectroscopy: it is both a thermal camera with full 2D imaging capabilities and a slit spectrometer. THERMAP takes advantage of the recent technological developments of uncooled microbolometers detectors, sensitive in the mid-infrared spectral range. THERMAP can acquire thermal images (8-18 {\mu}m) of the surface and perform absolute temperature measurements with a precision better than 3.5 K above 200 K. THERMAP can acquire mid-infrared spectra (8-16 {\mu}m) of the surface with a spectral resolution {\Delta}{\lambda} of 0.3 {\mu}m. For surface temperatures above 350 K, spectra have a signal-to-noise ratio >60 in the spectral range 9-13 {\mu}m where most emission features occur.

We analyze two multi-chord stellar occultations by Pluto observed on July 18th, 2012 and May 4th, 2013, and monitored respectively from five and six sites. They provide a total of fifteen light-curves, twelve of them being used for a simultaneous fit that uses a unique temperature profile, assuming a clear (no-haze) and pure N_2 atmosphere, but allowing for a possible pressure variation between the two dates. We find a solution that fits satisfactorily (i.e. within the noise level) all the twelve light-curves, providing atmospheric constraints between ~1,190 km (pressure ~ 11 \mubar) and ~ 1,450 km (pressure ~0.1 \mubar) from Pluto's center. Our main results are: (1) the best-fitting temperature profile shows a stratosphere with strong positive gradient between 1,190 km (at 36 K, 11 \mubar) and r = 1,215 km (6.0 \mubar), where a temperature maximum of 110 K is reached; above it is a mesosphere with negative thermal gradient of -0.2 K/km up to ~ 1,390 km (0.25 \mubar), where, the mesosphere connects itself to a more isothermal upper branch around 81 K; (2) the pressure shows a small (6 %) but significant increase (6-\sigma level) between the two dates; (3) without troposphere, Pluto's radius is found to be R_P = 1,190 +/- 5km. Allowing for a troposphere, R_P is constrained to lie between 1,168 and 1,195 km; (4) the currently measured CO abundance is too small to explain the mesospheric negative thermal gradient. Cooling by HCN is possible, but only if this species is largely saturated; Alternative explanations like zonal winds or vertical compositional variations of the atmosphere are unable to explain the observed mesospheric trend.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

The Planck High Frequency Instrument (HFI) has observed the full sky at six frequencies (100, 143, 217, 353, 545, and 857 GHz) in intensity and at four frequencies in linear polarization (100, 143, 217, and 353 GHz). In order to obtain sky maps, the time-ordered information (TOI) containing the detector and pointing samples must be processed and the angular response must be assessed. The full mission TOI is included in the Planck 2015 release. This paper describes the HFI TOI and beam processing for the 2015 release. HFI calibration and map-making are described in a companion paper. The main pipeline has been modified since the last release (2013 nominal mission in intensity only), by including a correction for the non-linearity of the warm readout and by improving the model of the bolometer time response. The beam processing is an essential tool that derives the angular response used in all the Planck science papers and we report an improvement in the effective beam window function uncertainty of more than a factor 10 relative to the 2013 release. Noise correlations introduced by pipeline filtering function are assessed using dedicated simulations. Angular cross-power spectra using datasets that are decorrelated in time are immune to the main systematic effects.

This paper describes the processing applied to the Planck High Frequency Instrument (HFI) cleaned, time-ordered information to produce photometrically calibrated maps in temperature and (for the first time) in polarization. The data from the entire 2.5 year HFI mission include almost five independent full-sky surveys. HFI observes the sky over a broad range of frequencies, from 100 to 857 GHz. To obtain the best accuracy on the calibration over such a large range, two different photometric calibration schemes have been used. The 545 and 857 GHz data are calibrated using models of planetary atmospheric emission. The lower frequencies (from 100 to 353 GHz) are calibrated using the time-variable cosmological microwave background dipole, which we call the "orbital dipole". This source of calibration only depends on the satellite velocity with respect to the solar system. Using a CMB temperature of 2.7255 +/- 0.0006 K, it permits an independent measurement of the amplitude of the CMB solar dipole (3364.3 +/- 1.5 \mu K) which is approximatively 1\sigma\ higher than the WMAP measurement with a direction that is consistent between both experiments. We describe the pipeline used to produce the maps of intensity and linear polarization from the HFI timelines, and the scheme used to set the zero level of the maps a posteriori. We also summarize the noise characteristics of the HFI maps in the 2015 Planck data release and present some null tests to assess their quality. Finally, we discuss the major systematic effects and in particular the leakage induced by flux mismatch between the detectors that leads to spurious polarization signal.

In this letter we explore the environment of Pluto and Charon in the far infrared with the main aim to identify the signs of any possible dust ring, should it exist in the system. Our study is based on observations performed at 70 um with the PACS instrument onboard the Herschel Space Observatory at 9 epochs between March 14 and 19, 2012. The far-infrared images of the Pluto-Charon system are compared to those of the point spread function (PSF) reference quasar 3C454.3. The deviation between the observed Pluto-Charon and reference PSFs are less then 1 sigma indicating that clear evidence for an extended dust ring around the system was not found. Our method is capable of detecting a hypothetical ring with a total flux of ~3.3 mJy at a distance of ~153 000 km (~8.2 Pluto-Charon distances) from the system barycentre. We place upper limits on the total disk mass and on the column density in a reasonable disk configuration and analyse the hazard during the flyby of NASAs New Horizons in July 2015. This realistic model configuration predicts a column density of 8.7x10^(-10) gcm^(-2) along the path of the probe and an impactor mass of 8.7x10^(-5) g.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

We report on high-resolution and spatially-resolved spectra of Io in the 4.0 {\mu}m region, recorded with the VLT/CRIRES instrument in 2008 and 2010, which provide the first detection of the {\nu}1 + {\nu}3 band of SO2 in Io's atmosphere. Data are analyzed to constrain the latitudinal, longitudinal, and diurnal distribution of Io's SO2 atmosphere as well as its characteristic temperature. equatorial SO2 column densities clearly show longitudinal asymmetry, but with a maximum of around 1.5e17 cm-2 at central meridian longitude L = 200-220 and a minimum of around 3e16 cm-2 at L = 285-300, the longitudinal pattern somewhat differs from earlier inferences from Ly {\alpha} and thermal IR measurements. Within the accuracy of the measurements, no evolution of the atmospheric density from mid-2008 to mid-2010 can be distinguished. The decrease of the SO2 column density towards high latitude is apparent, and the typical latitudinal extent of the atmosphere found to be (+-) 40{\deg} at half-maximum. The data show moderate diurnal variations of the equatorial atmosphere, which is evidence for a partially sublimation-supported atmospheric component. Compared to local noon, factor of 2 lower densities are observed around 40{\deg} before and 80{\deg} after noon. Best-fit gas temperatures range from 150 to 220 K, with a weighted mean value of 170 (+-) 20 K, which should represent the column-weighted mean kinetic temperature of Io's atmosphere. Finally, although the data include clear thermal emission due to Pillan (in outburst in July 2008) and Loki, no detectable enhancements in the SO2 atmosphere above these volcanic regions are found, with an upper limit of 4e16 cm-2 at Pillan and 1e17 cm-2 at Loki.

We present an analysis of the VIMS solar occultations dataset, which allows us to extract vertically resolved information on the characteristics of Titan's atmosphere between 100-700 km with a characteristic vertical resolution of 10 km. After a series of data treatment procedures, 4 occultations out of 10 are retained. This sample covers different seasons and latitudes of Titan. The transmittances show clearly the evolution of the haze and detect the detached layer at 310 km in Sept. 2011 at mid-northern latitudes. Through the inversion of the transmission spectra with a line-by-line radiative transfer code we retrieve the vertical distribution of CH$_4$ and CO mixing ratio. The two methane bands at 1.4 and 1.7 {\mu}m are always in good agreement and yield an average stratospheric abundance of $1.28\pm0.08$%. This is significantly less than the value of 1.48% obtained by the GCMS/Huygens instrument. The analysis of the residual spectra after the inversion shows that there are additional absorptions which affect a great part of the VIMS wavelength range. We attribute many of these additional bands to gaseous ethane, whose near-infrared spectrum is not well modeled yet. Ethane contributes significantly to the strong absorption between 3.2-3.5 {\mu}m that was previously attributed only to C-H stretching bands from aerosols. Ethane bands may affect the surface windows too, especially at 2.7 {\mu}m. Other residual bands are generated by stretching modes of C-H, C-C and C-N bonds. In addition to the C-H stretch from aliphatic hydrocarbons at 3.4 {\mu}m, we detect a strong and narrow absorption at 3.28 {\mu}m which we tentatively attribute to the presence of PAHs in the stratosphere. C-C and C-N stretching bands are possibly present between 4.3-4.5 {\mu}m. Finally, we obtain the CO mixing ratio between 70-170 km. The average result of $46\pm16$ ppm is in good agreement with previous studies.

Until now, rings have been detected in the Solar System exclusively around the four giant planets. Here we report the discovery of the first minor-body ring system around the Centaur object (10199) Chariklo, a body with equivalent radius 124$\pm$9 km. A multi-chord stellar occultation revealed the presence of two dense rings around Chariklo, with widths of about 7 km and 3 km, optical depths 0.4 and 0.06, and orbital radii 391 and 405 km, respectively. The present orientation of the ring is consistent with an edge-on geometry in 2008, thus providing a simple explanation for the dimming of Chariklo's system between 1997 and 2008, and for the gradual disappearance of ice and other absorption features in its spectrum over the same period. This implies that the rings are partially composed of water ice. These rings may be the remnants of a debris disk, which were possibly confined by embedded kilometre-sized satellites.

We analyze albedo data obtained using the Herschel Space Observatory that reveal the existence of two distinct types of surface among midsized transneptunian objects. A color-albedo diagram shows two large clusters of objects, one redder and higher albedo and another darker and more neutrally colored. Crucially, all objects in our sample located in dynamically stable orbits within the classical Kuiper belt region and beyond are confined to the bright-red group, implying a compositional link. Those objects are believed to have formed further from the Sun than the dark-neutral bodies. This color-albedo separation is evidence for a compositional discontinuity in the young solar system.

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations.

Aims: We investigate the composition of Titan's stratosphere from new medium-resolution far-infrared observations performed with the full range of Herschel's Photodetector Array Camera and Spectrometer (PACS) (51-220 $\mu$m at a resolution $\lambda$/$\Delta \lambda$ ranging from 950 to 5500 depending on wavelength and grating order). Methods: Using PACS, we obtained the spectral emission of several features of the Titan's stratosphere. We used a line-by-line radiative transfer code and the least-squares fitting technique to infer the abundances of the trace constituents. Results: Numerous spectral features attributable to CH$_4$, CO, HCN, and H$_2$O are present. From the flux density spectrum measured and by a detailed comparison with synthetic spectra, we constrain the stratospheric abundance of CH4, which is assumed to be constant with altitude, to be 1.29 $\pm$ 0.03%. Similarly, we constrain the abundance of CO to be 50 $\pm$ 2 ppm, and the HCN vertical distribution consistent with an increase from 40 ppb at $\sim$100 km to 4 ppm at $\sim$200 km, which is an altitude region where the HCN signatures are sensitive. Measurements of three H$_2$O rotational lines confirm the H$_2$O distribution profile recently obtained with Herschel. Furthermore, we determine the isotopic ratios $^{12}$C/$^{13}$C in CO and HCN to be 124 $\pm$ 58, and 66 $\pm$ 35, respectively. Comparisons between our results and the values derived with other instruments show that our results are consistent with the vertical distributions and isotopic ratios in previous studies, except for the HCN distribution obtained with Cassini/CIRS, which does not agree with the PACS lines at the 1-sigma confidence interval.

The classical Kuiper belt contains objects both from a low-inclination, presumably primordial, distribution and from a high-inclination dynamically excited population. Based on a sample of classical TNOs with observations at thermal wavelengths we determine radiometric sizes, geometric albedos and thermal beaming factors as well as study sample properties of dynamically hot and cold classicals. Observations near the thermal peak of TNOs using infra-red space telescopes are combined with optical magnitudes using the radiometric technique with near-Earth asteroid thermal model (NEATM). We have determined three-band flux densities from Herschel/PACS observations at 70.0, 100.0 and 160.0 $\mu$m and Spitzer/MIPS at 23.68 and 71.42 $\mu$m when available. We have analysed 18 classical TNOs with previously unpublished data and re-analysed previously published targets with updated data reduction to determine their sizes and geometric albedos as well as beaming factors when data quality allows. We have combined these samples with classical TNOs with radiometric results in the literature for the analysis of sample properties of a total of 44 objects. We find a median geometric albedo for cold classical TNOs of 0.14 and for dynamically hot classical TNOs, excluding the Haumea family and dwarf planets, 0.085. We have determined the bulk densities of Borasisi-Pabu (2.1 g/cm^3), Varda-Ilmare (1.25 g/cm^3) and 2001 QC298 (1.14 g/cm^3) as well as updated previous density estimates of four targets. We have determined the slope parameter of the debiased cumulative size distribution of dynamically hot classical TNOs as q=2.3 +- 0.1 in the diameter range 100<D<500 km. For dynamically cold classical TNOs we determine q=5.1 +- 1.1 in the diameter range 160<D<280 km as the cold classical TNOs have a smaller maximum size.

High-resolution spectra of Pluto in the 1.66 um region, recorded with the VLT/CRIRES instrument in 2008 (2 spectra) and 2012 (5 spectra), are analyzed to constrain the spatial and vertical distribution of methane in Pluto's atmosphere and to search for mid-term (4 year) variability. A sensitivity study to model assumptions (temperature structure, surface pressure, Pluto's radius) is performed. Results indicate that (i) no variation of the CH4 atmospheric content (column-density or mixing ratio) with Pluto rotational phase is present in excess of 20 % (ii) CH4 column densities show at most marginal variations between 2008 and 2012, with a best guess estimate of a ~20 % decrease over this time frame. As stellar occultations indicate that Pluto's surface pressure has continued to increase over this period, this implies a concomitant decrease of the methane mixing ratio (iii) the data do not show evidence for an altitude-varying methane distribution; in particular, they imply a roughly uniform mixing ratio in at least the first 22-27 km of the atmosphere, and high concentrations of low-temperature methane near the surface can be ruled out. Our results are also best consistent with a relatively large (> 1180 km) Pluto radius. Comparison with predictions from a recently developed global climate model GCM indicates that these features are best explained if the source of methane occurs in regional-scale CH4 ice deposits, including both low latitudes and high Northern latitudes, evidence for which is present from the rotational and secular evolution of the near-IR features due to CH4 ice. Our "best guess" predictions for the New Horizons encounter in 2015 are: a 1184 km radius, a 17 ubar surface pressure, and a 0.44 % CH4 mixing ratio with negligible longitudinal variations.

Centaurs are the transitional population between trans-Neptunian objects (TNOs) and Jupiter-family comets. For this reason it is possible to access the smaller ones, which is more difficult to do with the TNO population. The goal of this work is to characterize a set of 16 Centaurs in terms of their size, albedo, and thermal properties. We study the correlations, for a more extended sample obtained from the literature, of diameter, albedo, orbital parameters, and spectral slopes. We performed three-band photometric observations using Herschel-PACS and used a consistent method for the data reduction and aperture photometry of this sample to obtain monochromatic flux densities at 70, 100, and 160 $\mu$m. Additionally, we used Spitzer-MIPS flux densities at 24 and 70 $\mu$m when available. We also included in our Centaur sample scattered disk objects (SDOs), a dynamical family of TNOS, using results previously published by our team, and some Centaurs observed only with the Spitzer/MIPS instrument. We have determined new radiometric sizes and albedos of 16 Centaurs. The first conclusion is that the albedos of Centaur objects are not correlated with their orbital parameters. Similarly, there is no correlation between diameter and orbital parameters. Most of the objects in our sample are dark (pv < 7%) and most of them are small (D < 120km). However, there is no correlation between albedo and diameter, in particular for the group of the small objects as albedo values are homogeneously distributed between 4 to 16%. The correlation with the color of the objects showed that red objects are all small (mean diameter 65 km), while the gray ones span a wide range of sizes (mean diameter 120 km). Moreover, the gray objects tend to be darker, with a mean albedo of 5.6%, compared with a mean of 8.5% (ranging from 5 to 15%) for the red objects.

This paper describes the processing applied to the HFI cleaned time-ordered data to produce photometrically calibrated maps. HFI observes the sky over a broad range of frequencies, from 100 to 857 GHz. To get the best accuracy on the calibration on such a large range, two different photometric calibration schemes have to be used. The 545 and 857 \GHz\ data are calibrated using Uranus and Neptune flux density measurements, compared with models of their atmospheric emissions to calibrate the data. The lower frequencies (below 353 GHz) are calibrated using the cosmological microwave background dipole.One of the components of this anisotropy results from the orbital motion of the satellite in the Solar System, and is therefore time-variable. Photometric calibration is thus tightly linked to mapmaking, which also addresses low frequency noise removal. The 2013 released HFI data show some evidence for apparent gain variations of the HFI bolometers' detection chain. These variations were identified by comparing observations taken more than one year apart in the same configuration. We developed an effective correction to limit its effect on calibration, and assess its accuracy. We present several methods used to estimate the precision of the photometric calibration. We distinguish relative (from one detector to another, or from one frequency to another) and absolute uncertainties. In both cases, we found that these uncertainties range from a few $10^{-3}$ to several per cents from 100 to 857 GHz. We describe the pipeline producing the maps from the HFI timelines, based on the photometric calibration parameters and we detail the scheme used to a posteriori set the zero level of the maps. We also briefly discuss the cross-calibration between HFI and the SPIRE instrument on board Herschel. We finally summarize the basic characteristics of the set of the HFI maps from the 2013 Planck data release.

Context. Carbon monoxide (CO) has been detected in all Giant Planets and its origin is both internal and external in Jupiter and Neptune. Despite its first detection in Uranus a decade ago, the magnitude of its internal and external sources remains unconstrained. Aims. We targeted CO lines in Uranus in the submillimeter range to constrain its origin. Methods. We recorded disk-averaged spectra of Uranus with a very high spectral resolution at the frequencies of CO rotational lines in the submillimeter range in 2011-2012. We used empirical and diffusion models of the atmosphere of Uranus to constrain the origin of CO. We also used a thermochemical model of its troposphere to derive an upper limit on the O/H ratio in the deep atmosphere of Uranus. Results. We have detected the CO(8-7) rotational line for the first time, with Herschel-HIFI. Both empirical and diffusion model results show that CO has an external origin. An empirical profile in which CO is constant above the 100mbar level with a mole fraction of 7.1-9.0E-9 (depending on the adopted stratospheric thermal structure) reproduces the data. Sporadic and steady source models cannot be differentiated with our data. Taking Teanby & Irwin's internal source model upper limit of a mole fraction of 2.1E-9 [Teanby & Irwin 2013. ApJ, 775, L49], the deep O/H ratio of Uranus is lower than 500 times solar according to our thermochemical computations. Conclusions. Our work shows that the average mole fraction of CO decreases from the stratosphere to the troposphere and thus strongly advocates for an external source of CO in Uranus. Photochemical modeling of oxygen species in the atmosphere of Uranus and more sensitive observations are needed to reveal the nature of the external source.

The "TNOs are Cool!: A survey of the trans-Neptunian region" is a Herschel Open Time Key Program that aims to characterize planetary bodies at the outskirts of the Solar System using PACS and SPIRE data, mostly taken as scan-maps. In this paper we summarize our PACS data reduction scheme that uses a modified version of the standard pipeline for basic data reduction, optimized for faint, moving targets. Due to the low flux density of our targets the observations are confusion noise limited or at least often affected by bright nearby background sources at 100 and 160\,$\mu$m. To overcome these problems we developed techniques to characterize and eliminate the background at the positions of our targets and a background matching technique to compensate for pointing errors. We derive a variety of maps as science data products that are used depending on the source flux and background levels and the scientific purpose. Our techniques are also applicable to a wealth of other Herschel solar system photometric observations, e.g. comets and near-Earth asteroids. The principles of our observing strategies and reduction techniques for moving targets will also be applicable for similar surveys of future infrared space projects.

The composition of Io's tenuous atmosphere is poorly constrained. Only the major species SO2 and a handful of minor species have been positively identified, but a variety of other molecular species should be present, based on thermochemical equilibrium models of volcanic gas chemistry and the composition of Io's environment. This paper focuses on the spectral search for expected yet undetected molecular species (KCl, SiO, S2O) and isotopes (34SO2). We analyze a disk-averaged spectrum of a potentially line-rich spectral window around 345 GHz, obtained in 2010 at the APEX-12m antenna (Atacama Pathfinder EXperiment). Using different models assuming either extended atmospheric distributions or a purely volcanically-sustained atmosphere, we tentatively measure the KCl relative abundance with respect to SO2 and derive a range of 4x10^{-4}-8x10^{-3}. We do not detect SiO or S2O and present new upper limits on their abundances. We also present the first measurement of the 34S/32S isotopic ratio in gas phase on Io, which appears to be twice as high as the Earth and ISM reference values. Strong lines of SO2 and SO are also analyzed to check for longitudinal variations of column density and relative abundance. Our models show that, based on their predicted relative abundance with respect to SO2 in volcanic plumes, both the tentative KCl detection and SiO upper limit are compatible with a purely volcanic origin for these species.

Transneptunian objects (TNOs) are bodies populating the Kuiper Belt and they are believed to retain the most pristine and least altered material of the solar system. The Herschel Open Time Key Program entitled "TNOs are Cool: A survey of the trans-Neptunian region" has been awarded 373 h to investigate the albedo, size distribution and thermal properties of TNOs and Centaurs. Here we focus on the brightest targets observed by both the PACS and SPIRE multiband photometers: the dwarf planet Haumea, six TNOs (Huya, Orcus, Quaoar, Salacia, 2002 UX25, and 2002 TC302), and two Centaurs (Chiron and Chariklo). Flux densities are derived from PACS and SPIRE instruments using optimised data reduction methods. The spectral energy distribution obtained with the Herschel PACS and SPIRE instruments over 6 bands (centred at 70, 100, 160, 250, 350, and 500 $\mu$m), and with Spitzer-MIPS at 23.7 and 71.4 $\mu$m has been modelled with the NEATM thermal model in order to derive the albedo, diameter, and beaming factor. For the Centaurs Chiron and Chariklo and for the 1000 km sized Orcus and Quaoar, a thermophysical model was also run to better constrain their thermal properties. We derive the size, albedo, and thermal properties, including thermal inertia and surface emissivity, for the 9 TNOs and Centaurs. Several targets show a significant decrease in their spectral emissivity longwards of $\sim$300 $\mu$m and especially at 500 $\mu$m. Using our size estimations and the mass values available in the literature, we also derive the bulk densities for the binaries Quaoar/Weywot (2.18$^{+0.43}_{-0.36}$ g/cm$^3$), Orcus/Vanth (1.53$^{+0.15}_{-0.13}$ g/cm$^3$), and Salacia/Actea (1.29$^{+0.29}_{-0.23}$ g/cm$^3$). Quaoar's density is similar to that of the other dwarf planets Pluto and Haumea, and its value implies high contents of refractory materials mixed with ices.

Herschel-PACS measurements of the rotational R(0) and R(1) HD lines in the atmospheres of Uranus and Neptune are analyzed in order to derive a D/H ratio with improved precision for both planets. The derivation of the D/H ratio includes also previous measurements of the R(2) line by the Short Wavelength Spectrometer on board the Infrared Space Observatory (ISO). The available spectroscopic line information of the three rotational transitions is discussed and applied in the radiative transfer calculations. The best simultaneous fit of all three lines requires only a minor departure from the Spitzer temperature profile of Uranus and a departure limited to 2K from the Voyager temperature profile of Neptune (both around the tropopause). The resulting and remarkably similar D/H ratios for Uranus and Neptune are found to be (4.4$\pm$0.4)$\times10^{-5}$ and (4.1$\pm$0.4)$\times10^{-5}$ respectively. Although the deuterium enrichment in both atmospheres compared to the protosolar value is confirmed, it is found to be lower compared to previous analysis. Using the interior models of Podolak et al. (1995), Helled et al. (2011) and Nettelmann et al. (2013), and assuming that complete mixing of the atmosphere and interior occured during the planets history, we derive a D/H in protoplanetary ices between (5.75--7.0)$\times10^{-5}$ for Uranus and between (5.1--7.7)$\times10^{-5}$ for Neptune. Conversely, adopting a cometary D/H for the protoplanetary ices between (15-30)$\times10^{-5}$, we constrain the interior models of both planets to have an ice mass fraction of 14-32%, i.e. that the two planets are rock-dominated.

176P/LINEAR is a member of the new cometary class known as main-belt comets (MBCs). It displayed cometary activity shortly during its 2005 perihelion passage that may be driven by the sublimation of sub-surface ices. We have therefore searched for emission of the H2O 110-101 ground state rotational line at 557 GHz toward 176P/LINEAR with the Heterodyne Instrument for the Far Infrared (HIFI) on board the Herschel Space Observatory on UT 8.78 August 2011, about 40 days after its most recent perihelion passage, when the object was at a heliocentric distance of 2.58 AU. No H2O line emission was detected in our observations, from which we derive sensitive 3-sigma upper limits for the water production rate and column density of < 4e25 molec/s and of < 3e10 cm^{-2}, respectively. From the peak brightness measured during the object's active period in 2005, this upper limit is lower than predicted by the relation between production rates and visual magnitudes observed for a sample of comets by Jorda et al. (2008) at this heliocentric distance. Thus, 176P/LINEAR was likely less active at the time of our observation than during its previous perihelion passage. The retrieved upper limit is lower than most values derived for the H2O production rate from the spectroscopic search for CN emission in MBCs.

The D/H ratio in cometary water is believed to be an important indicator of the conditions under which icy planetesimals formed and can provide clues to the contribution of comets to the delivery of water and other volatiles to Earth. Available measurements suggest that there is isotopic diversity in the comet population. The Herschel Space Observatory revealed an ocean-like ratio in the Jupiter-family comet 103P/Hartley 2, whereas most values measured in Oort-cloud comets are twice as high as the ocean D/H ratio. We present here a new measurement of the D/H ratio in the water of an Oort-cloud comet. HDO, H_2O, and H_2^18O lines were observed with high signal-to-noise ratio in comet C/2009 P1 (Garradd) using the Herschel HIFI instrument. Spectral maps of two water lines were obtained to constrain the water excitation. The D/H ratio derived from the measured H_2^16O and HDO production rates is 2.06+/-0.22 X 10**-4. This result shows that the D/H in the water of Oort-cloud comets is not as high as previously thought, at least for a fraction of the population, hence the paradigm of a single, archetypal D/H ratio for all Oort-cloud comets is no longer tenable. Nevertheless, the value measured in C/2009 P1 (Garradd) is significantly higher than the Earth's ocean value of 1.558 X 10**-4. The measured H_2^16O/H_2^18O ratio of 523+/-32 is, however, consistent with the terrestrial value.

Comet 17P/Holmes underwent a dramatic outburst in October 2007, caused by the sudden fragmentation of its nucleus and the production of a large quantity of grains scattering sunlight. We report on 90 GHz continuum observations carried out with the IRAM Plateau de Bure interferometer on 27.1 and 28.2 October 2007 UT, i.e., 4-5 days after the outburst. These observations probed the thermal radiation of large dust particles, and therefore provide the best constraints on the mass in the ejecta debris. The thermal emission of the debris was modelled and coupled to a time-dependent description of their expansion after the outburst. The analysis was performed in the Fourier plane. Visibilities were computed for the two observing dates and compared to the data to measure their velocity and mass. Optical data and 250-GHz continuum measurements published in the literature were used to further constrain the dust kinematics and size distribution. Two distinct dust components in terms of kinematic properties are identified in the data. The large-velocity component, with typical velocities V0 of 50-100 m/s for 1 mm particles, displays a steep size distribution with a size index estimated to q = -3.7 (\pm0.1), assuming a minimum grain size of 0.1 \mum. It corresponds to the fast expanding shell observed in optical images. The slowly-moving "core" component (V0 = 7-9 m/s) detected near the nucleus has a size index |q| < 3.4 and contains a higher proportion of large particles than the shell. The dust mass in the core is in the range 0.1-1 that of the shell. Using optical constants pertaining to porous grains (50% porosity) made of astronomical silicates mixed with water ice (48% in mass), the total dust mass Mdust injected by the outburst is estimated to 4-14 x 10**11 kg, corresponding to 3-9% the nucleus mass.

We present estimates of the basic physical properties (size and albedo) of (90377) Sedna, a prominent member of the detached trans-Neptunian object population and the recently discovered scattered disk object 2010 EK139, based on the recent observations acquired with the Herschel Space Observatory, within the "TNOs are Cool!" key programme. Our modeling of the thermal measurements shows that both objects have larger albedos and smaller sizes than the previous expectations, thus their surfaces might be covered by ices in a significantly larger fraction. The derived diameter of Sedna and 2010 EK139 are 995 +/- 80 km and 470 +35/-10 km, while the respective geometric albedos are pV 0.32 +/- 0.06 and 0.25 +0.02/-0.05. These estimates are based on thermophysical model techniques.

Trans-Neptunian objects (TNO) represent the leftovers of the formation of the Solar System. Their physical properties provide constraints to the models of formation and evolution of the various dynamical classes of objects in the outer Solar System. Based on a sample of 19 classical TNOs we determine radiometric sizes, geometric albedos and beaming parameters. Our sample is composed of both dynamically hot and cold classicals. We study the correlations of diameter and albedo of these two subsamples with each other and with orbital parameters, spectral slopes and colors. We have done three-band photometric observations with Herschel/PACS and we use a consistent method for data reduction and aperture photometry of this sample to obtain monochromatic flux densities at 70.0, 100.0 and 160.0 \mu m. Additionally, we use Spitzer/MIPS flux densities at 23.68 and 71.42 \mu m when available, and we present new Spitzer flux densities of eight targets. We derive diameters and albedos with the near-Earth asteroid thermal model (NEATM). As auxiliary data we use reexamined absolute visual magnitudes from the literature and data bases, part of which have been obtained by ground based programs in support of our Herschel key program. We have determined for the first time radiometric sizes and albedos of eight classical TNOs, and refined previous size and albedo estimates or limits of 11 other classicals. The new size estimates of 2002 MS4 and 120347 Salacia indicate that they are among the 10 largest TNOs known. Our new results confirm the recent findings that there are very diverse albedos among the classical TNOs and that cold classicals possess a high average albedo (0.17 +/- 0.04). Diameters of classical TNOs strongly correlate with orbital inclination in our sample. We also determine the bulk densities of six binary TNOs.

The dynamics of the upper mesosphere of Venus (~85-115 km) have been characterized as a combination of a retrograde superrotating zonal wind (RSZ) with a subsolar-to-antisolar flow (SSAS). Numerous mm-wave single-dish observations have been obtained and could directly measure mesospheric line-of-sight winds by mapping Doppler-shifts on CO rotational lines, but their limited spatial resolution makes their interpretation difficult. By using interferometric facilities, one can obtain better resolution on Doppler-shifts maps, allowing in particular to put firmer constraints on the respective contributions of the SSAS and RSZ circulations to the global mesospheric wind field. We report on interferometric observations of the CO(1-0) line obtained with the IRAM-Plateau de Bure interferometer in November 2007 and June 2009, that could map the upper mesosphere dynamics on the morning hemisphere with a very good spatial resolution (3.5-5.5"). All the obtained measurements show, with a remarkably good temporal stability, that the wind globally flows in the (sky) East-West direction, corresponding in the observed geometry either to an unexpected prograde zonal wind or a SSAS flow. A very localized inversion of the wind direction, that could correspond to a RSZ wind, is also repeatedly detected in the night hemisphere. The presence of significant meridional winds is not evidenced. Using models with different combinations of zonal and SSAS winds, we find that the data is best reproduced by a dominant SSAS flow with a maximal velocity at the terminator of ~200 m/s, displaying large diurnal and latitudinal asymmetries, combined with an equatorial RSZ wind of 70-100 m/s, overall indicating a wind-field structure consistent with but much more complex than the usual representation of the mesospheric dynamics.

We present Herschel PACS photometry of 18 Plutinos and determine sizes and albedos for these objects using thermal modeling. We analyze our results for correlations, draw conclusions on the Plutino size distribution, and compare to earlier results. Flux densities are derived from PACS mini scan-maps using specialized data reduction and photometry methods. In order to improve the quality of our results, we combine our PACS data with existing Spitzer MIPS data where possible, and refine existing absolute magnitudes for the targets. The physical characterization of our sample is done using a thermal model. Uncertainties of the physical parameters are derived using customized Monte Carlo methods. The correlation analysis is performed using a bootstrap Spearman rank analysis. We find the sizes of our Plutinos to range from 150 to 730 km and geometric albedos to vary between 0.04 and 0.28. The average albedo of the sample is 0.08 \pm 0.03, which is comparable to the mean albedo of Centaurs, Jupiter Family comets and other Trans-Neptunian Objects. We were able to calibrate the Plutino size scale for the first time and find the cumulative Plutino size distribution to be best fit using a cumulative power law with q = 2 at sizes ranging from 120-400 km and q = 3 at larger sizes. We revise the bulk density of 1999 TC36 and find a density of 0.64 (+0.15/-0.11) g cm-3. On the basis of a modified Spearman rank analysis technique our Plutino sample appears to be biased with respect to object size but unbiased with respect to albedo. Furthermore, we find biases based on geometrical aspects and color in our sample. There is qualitative evidence that icy Plutinos have higher albedos than the average of the sample.

Physical characterization of Trans-Neptunian objects, a primitive population of the outer solar system, may provide constraints on their formation and evolution. The goal of this work is to characterize a set of 15 scattered disk (SDOs) and detached objects, in terms of their size, albedo, and thermal properties. Thermal flux measurements obtained with the Herschel-PACS instrument at 70, 100 and 160 \mu m, and whenever applicable, with Spitzer-MIPS at 24 and 70 \mu m, are modeled with radiometric techniques, in order to derive the objects' individual size, albedo and when possible beaming factor. Error bars are obtained from a Monte-Carlo approach. We look for correlations between these and other physical and orbital parameters. Diameters obtained for our sample range from 100 to 2400 km, and the geometric albedos (in V band) vary from 3.8 % to 84.5 %. The unweighted mean V geometric albedo for the whole sample is 11.2 % (excluding Eris); 6.9 % for the SDOs, and 17.0 % for the detached objects (excluding Eris). We obtain new bulk densities for three binary systems: Ceto/Phorcys, Typhon/Echidna and Eris/Dysnomia. Apart from correlations clearly due to observational bias, we find significant correlations between albedo and diameter (more reflective objects being bigger), and between albedo, diameter and perihelion distance (brighter and bigger objects having larger perihelia). We discuss possible explanations for these correlations.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

The goal is to determine the composition of Pluto's atmosphere and to constrain the nature of surface-atmosphere interactions. We perform high--resolution spectroscopic observations in the 2.33--2.36 $\mu$m range, using CRIRES at the VLT. We obtain (i) the first detection of gaseous methane in this spectral range, through lines of the $\nu_3$ + $\nu_4$ and $\nu_1$ + $\nu_4$ bands (ii) strong evidence (6-$\sigma$ confidence) for gaseous CO in Pluto. For an isothermal atmosphere at 90 K, the CH$_4$ and CO column densities are 0.75 and 0.07 cm-am, within factors of 2 and 3, respectively. Using a physically--based thermal structure model of Pluto's atmosphere also satisfying constraints from stellar occultations, we infer CH$_4$ and CO mixing ratios q$_{CH_4}$= 0.6$^{+0.6}_{-0.3}$% (consistent with results from the 1.66 $\mu$m range) and q$_{CO}$ = 0.5$^{+1}_{-0.25}$$\times10^{-3}$. The CO atmospheric abundance is consistent with its surface abundance. As for Triton, it is probably controlled by a thin, CO-rich, detailed balancing layer resulting from seasonal transport and/or atmospheric escape.

Once completed, the Atacama Large Millimeter Array (ALMA) will be the most powerful (sub)millimeter interferometer in terms of sensitivity, spatial resolution and imaging. This paper presents the capabilities of ALMA applied to the observation of Centaurs and Trans-Neptunian Objects, and their possible output in terms of physical properties. Realistic simulations were performed to explore the performances of the different frequency bands and array configurations, and several projects are detailed along with their feasibility, their limitations and their possible targets. Determination of diameters and albedos via the radiometric method appears to be possible on ~500 objects, while sampling of the thermal lightcurve to derive the bodies' ellipticity could be performed at least 30 bodies that display a significant optical lightcurve. On a limited number of objects, the spatial resolution allows for direct measurement of the size or even surface mapping with a resolution down to 13 milliarcseconds. Finally, ALMA could separate members of multiple systems with a separation power comparable to that of the HST. The overall performance of ALMA will make it an invaluable instrument to explore the outer solar system, complementary to space-based telescopes and spacecrafts.

Little is known about the physical properties of cometary nuclei. Apart from space mission targets, measuring the thermal emission of a nucleus is one of the few means to derive its size, independently of its albedo, and to constrain some of its thermal properties. This emission is difficult to detect from Earth but space telescopes (Infrared Space Observatory, Spitzer Space Telescope, Herschel Space Observatory) allow reliable measurements in the infrared and the sub-millimetre domains. We aim at better characterizing the thermal properties of the nucleus of comet 8P/Tuttle using multi-wavelentgh space- and ground-based observations, in the visible, infrared, and millimetre range. We used the Plateau de Bure Interferometer to measure the millimetre thermal emission of comet 8P/Tuttle at 240 GHz (1.25 mm) and analysed the observations with the shape model derived from Hubble Space Telescope observations and the nucleus size derived from Spitzer Space Telescope observations. We report on the first detection of the millimetre thermal emission of a cometary nucleus since comet C/1995 O1 Hale-Bopp in 1997. Using the two contact spheres shape model derived from Hubble Space Telescope observations, we constrained the thermal properties of the nucleus. Our millimetre observations are best match with: i) a thermal inertia lower than ~10 J K-1 m-2 s-1/2, ii) an emissivity lower than 0.8, indicating a non-negligible contribution of the colder sub-surface layers to the outcoming millimetre flux.

We report on the initial analysis of Herschel/HIFI carbon monoxide (CO) observations of the Martian atmosphere performed between 11 and 16 April 2010. We selected the (7-6) rotational transitions of the isotopes ^{13}CO at 771 GHz and C^{18}O at 768 GHz in order to retrieve the mean vertical profile of temperature and the mean volume mixing ratio of carbon monoxide. The derived temperature profile agrees within less than 5 K with general circulation model (GCM) predictions up to an altitude of 45 km, however, show about 12-15 K lower values at 60 km. The CO mixing ratio was determined as 980 \pm 150 ppm, in agreement with the 900 ppm derived from Herschel/SPIRE observations in November 2009.

We report on an initial analysis of Herschel/HIFI observations of hydrogen chloride (HCl), hydrogen peroxide (H_2O_2), and molecular oxygen (O_2) in the Martian atmosphere performed on 13 and 16 April 2010 (L_s ~ 77{\deg}). We derived a constant volume mixing ratio of 1400 +/- 120 ppm for O_2 and determined upper limits of 200 ppt for HCl and 2 ppb for H_2O_2. Radiative transfer model calculations indicate that the vertical profile of O_2 may not be constant. Photochemical models determine the lowest values of H_2O_2 to be around L_s ~ 75{\deg} but overestimate the volume mixing ratio compared to our measurements.

The high spectral resolution and sensitivity of Herschel/HIFI allows for the detection of multiple rotational water lines and accurate determinations of water production rates in comets. In this letter we present HIFI observations of the fundamental 110-101 (557 GHz) ortho and 111-000 (1113 GHz) para rotational transitions of water in comet 81P/Wild 2 acquired in February 2010. We mapped the extent of the water line emission with five point scans. Line profiles are computed using excitation models which include excitation by collisions with electrons and neutrals and solar infrared radiation. We derive a mean water production rate of $1.0 \times 10^{28}$ molecules s$^{-1}$ at a heliocentric distance of 1.61 AU about 20 days before perihelion, in agreement with production rates measured from the ground using observations of the 18-cm OH lines. Furthermore, we constrain the electron density profile and gas kinetic temperature, and estimate the coma expansion velocity by fitting the water line shapes.

We report on the initial analysis of a Herschel/PACS full range spectrum of Neptune, covering the 51-220 micrometer range with a mean resolving power of ~ 3000, and complemented by a dedicated observation of CH4 at 120 micrometers. Numerous spectral features due to HD (R(0) and R(1)), H2O, CH4, and CO are present, but so far no new species have been found. Our results indicate that (i) Neptune's mean thermal profile is warmer by ~ 3 K than inferred from the Voyager radio-occultation; (ii) the D/H mixing ratio is (4.5+/-1) X 10**-5, confirming the enrichment of Neptune in deuterium over the protosolar value (~ 2.1 X 10**-5); (iii) the CH4 mixing ratio in the mid stratosphere is (1.5+/-0.2) X 10**-3, and CH4 appears to decrease in the lower stratosphere at a rate consistent with local saturation, in agreement with the scenario of CH4 stratospheric injection from Neptune's warm south polar region; (iv) the H2O stratospheric column is (2.1+/-0.5) X 10**14 cm-2 but its vertical distribution is still to be determined, so the H2O external flux remains uncertain by over an order of magnitude; and (v) the CO stratospheric abundance is about twice the tropospheric value, confirming the dual origin of CO suspected from ground-based millimeter/submillimeter observations.

Thermal emission from Kuiper Belt object (136108) Haumea was measured with Herschel-PACS at 100 and 160 micrometers for almost a full rotation period. Observations clearly indicate a 100-micrometer thermal lightcurve with an amplitude of a factor of ~ 2, which is positively correlated with the optical lightcurve. This confirms that both are primarily due to shape effects. A 160-micrometer lightcurve is marginally detected. Radiometric fits of the mean Herschel- and Spitzer- fluxes indicate an equivalent diameter D ~ 1300 km and a geometric albedo p_v ~ 0.70-0.75. These values agree with inferences from the optical lightcurve, supporting the hydrostatic equilibrium hypothesis. The large amplitude of the 100-micrometer lightcurve suggests that the object has a high projected a/b axis ratio (~ 1.3) and a low thermal inertia as well as possible variable infrared beaming. This may point to fine regolith on the surface, with a lunar-type photometric behavior. The quality of the thermal data is not sufficient to clearly detect the effects of a surface dark spot.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

We have obtained the first continuous disk averaged spectrum of Mars from 450 to 1550 Ghz using the Herschel-SPIRE Fourier Transform Spectrometer. The spectrum was obtained at a constant resolution of 1.4 GHz across the whole band. The flux from the planet is such that the instrument was operated in "bright source" mode to prevent saturation of the detectors. This was the first successful use of this mode and in this work we describe the method used for observing Mars together with a detailed discussion of the data reduction techniques required to calibrate the spectrum. We discuss the calibration accuracy obtained and describe the first comparison with surface and atmospheric models. In addition to a direct photometric measurement of the planet the spectrum contains the characteristic transitions of 12CO from J 5-4 to J 13-12 as well as numerous H2O transitions. Together these allow the comparison to global atmospheric models allowing the mean mixing ratios of water and 12CO to be investigated. We find that it is possible to match the observed depth of the absorption features in the spectrum with a fixed water mixing ratio of 1 x 10-4 and a 12CO mixing ratio of 9 x 10-4

High-resolution far-infrared and sub-millimetre spectroscopy of water lines is an important tool to understand the physical and chemical properties of cometary atmospheres. We present observations of several rotational ortho- and para-water transitions in comet C/2008 Q3 (Garradd) performed with HIFI on Herschel. These observations have provided the first detection of the 2_{12}-1_{01} (1669 GHz) ortho and 1_{11}-0_{00} (1113 GHz) para transitions of water in a cometary spectrum. In addition, the ground-state transition 1_{10}-1_{01} at 557 GHz is detected and mapped. By detecting several water lines quasi-simultaneously and mapping their emission we can constrain the excitation parameters in the coma. Synthetic line profiles are computed using excitation models which include excitation by collisions, solar infrared radiation, and radiation trapping. We obtain the gas kinetic temperature, constrain the electron density profile, and estimate the coma expansion velocity by analyzing the map and line shapes. We derive water production rates of 1.7-2.8 x 10^{28} s^{-1} over the range r_h = 1.83-1.85 AU.

The goal of the Herschel Open Time Key programme "TNOs are Cool!" is to derive the physical and thermal properties for a large sample of Centaurs and trans-Neptunian objects (TNOs), including resonant, classical, detached and scattered disk objects. We present results for seven targets either observed in PACS point-source, or in mini scan-map mode. Spitzer-MIPS observations were included for three objects. The sizes of these targets range from 100 km to almost 1000 km, five have low geometric albedos below 10%, (145480) 2005 TB190 has a higher albedo above 15%. Classical thermal models driven by an intermediate beaming factor of $\eta$=1.2 or $\eta$-values adjusted to the observed colour temperature fit the multi-band observations well in most cases. More sophisticated thermophysical models give very similar diameter and albedo values for thermal inertias in the range 0-25 Jm-2s-0.5K-1, consistent with very low heat conductivities at temperatures far away from the Sun. The early experience with observing and model strategies will allow us to derive physical and thermal properties for our complete Herschel TNO sample of 140 targets as a benchmark for understanding the solar system debris disk, and extra-solar ones as well.

Comet C/2006 W3 (Christensen) was observed in November 2009 at 3.3 AU from the Sun with Herschel. The PACS instrument acquired images of the dust coma in 70- and 160-micrometers filters, and spectra covering several H2O rotational lines. Spectra in the range 450-1550 GHz were acquired with SPIRE. The comet emission continuum from 70 to 672 micrometers was measured, but no lines were detected. The spectral energy distribution indicates thermal emission from large particles and provides a measure of the size distribution index and dust production rate. The upper limit to the water production rate is compared to the production rates of other species (CO, CH3OH, HCN, H2S, OH) measured with the IRAM 30-m and Nancay telescopes. The coma is found to be strongly enriched in species more volatile than water, in comparison to comets observed closer to the Sun. The CO to H2O production rate ratio exceeds 220%. The dust to gas production rate ratio is on the order of 1.

Triton possesses a thin atmosphere, primarily composed of nitrogen, sustained by the sublimation of surface ices. The goal is to determine the composition of Triton's atmosphere and to constrain the nature of surface-atmosphere interactions. We perform high-resolution spectroscopic observations in the 2.32-2.37 $\mu$m range, using CRIRES at the VLT. From this first spectroscopic detection of Triton's atmosphere in the infrared, we report (i) the first observation of gaseous methane since its discovery in the ultraviolet by Voyager in 1989 and (ii) the first ever detection of gaseous CO in the satellite. The CO atmospheric abundance is remarkably similar to its surface abundance, and appears to be controlled by a thin, CO-enriched, surface veneer resulting from seasonal transport and/or atmospheric escape. The CH$_4$ partial pressure is several times larger than inferred from Voyager. This confirms that Triton's atmosphere is seasonally variable and is best interpreted by the warming of CH$_4$-rich icy grains as Triton passed southern summer solstice in 2000. The presence of CO in Triton's atmosphere also affects its temperature, photochemistry and ionospheric composition. An improved upper limit on CO in Pluto's atmosphere is also reported.

Context: Pluto possesses a thin atmosphere, primarily composed of nitrogen, in which the detection of methane has been reported. Aims: The goal is to constrain essential but so far unknown parameters of Pluto's atmosphere such as the surface pressure, lower atmosphere thermal stucture, and methane mixing ratio. Methods: We use high-resolution spectroscopic observations of gaseous methane, and a novel analysis of occultation light-curves. Results: We show that (i) Pluto's surface pressure is currently in the 6.5-24 microbar range (ii) the methane mixing ratio is 0.5+/-0.1 %, adequate to explain Pluto's inverted thermal structure and ~100 K upper atmosphere temperature (iii) a troposphere is not required by our data, but if present, it has a depth of at most 17 km, i.e. less than one pressure scale height; in this case methane is supersaturated in most of it. The atmospheric and bulk surface abundance of methane are strikingly similar, a possible consequence of the presence of a CH4-rich top surface layer.

Characterizing the Martian atmosphere is an essential objective to understand its meteorology and its climate. The lower atmosphere (< 40 km) and middle atmosphere (40-80 km) of Mars appear dynamically coupled at much higher levels than in the case of the Earth. The vertical extension of the weather phenomena is considerable with for example Hadley's cells reaching the top of the neutral atmosphere (120 km). The circulation in the middle atmosphere modifies the meteorology of the lower atmosphere, affecting the transport and climatic processes Observations of the CO rotational lines at millimeter (mm) wavelengths (Clancy et al 1990) have strongly contributed in the study of the vertical distribution of this compound and the thermal profile in the atmosphere of Mars over 0-70 km. Singledish observations of the CO Doppler lineshift have allowed direct wind measurements in the martian middle atmosphere near 50 km altitude (Lellouch et al 1991), but at a low spatial resolution (12'') enabling only an essentially hemispheric resolution of the martian disk. The use of mm interferometry has allowed us to better spatially resolve the Martian disk, in order to obtain wind maps of the middle atmosphere (Moreno et al 2001).

We report on coordinated thermal and optical measurements of trans-Neptunian object (20000) Varuna obtained in January-February 2002, respectively from the IRAM 30-m and IAA 1.5 m telescopes. The optical data show a lightcurve with a period of 3.176+/-0.010 hr, a mean V magnitude of 20.37+/-0.08 and a 0.42+/-0.01 magnitude amplitude. They also tentatively indicate that the lightcurve is asymmetric and double-peaked. The thermal observations indicate a 1.12+/-0.41 mJy flux, averaged over the object's rotation. Combining the two datasets, we infer that Varuna has a mean 1060(+180/-220) km diameter and a mean 0.038(+0.022/-0.010) V geometric albedo, in general agreement with an earlier determination using the same technique.

We have analysed the ISO-SWS spectrum of Jupiter in the 12-16 micron range, where several hydrocarbons exhibit rovibrational bands. Using temperature information from the methane and hydrogen emissions, we derive the mixing ratios (q) of acetylene and ethane at two independent pressure levels. For acetylene, we find $q=(8.9^{+1.1}_{-0.6})\times10^{-7}$ at 0.3 mbar and $q=(1.1^{+0.2}_{-0.1})\times10^{-7}$ at 4 mbar, giving a slope $-d\ln q / d\ln P=0.8\pm0.1$, while for ethane $q=(1.0\pm0.2)\times10^{-5}$ at 1 mbar and $q=(2.6^{+0.5}_{-0.6})\times10^{-6}$ at 10 mbar, giving $-d\ln q / d\ln P=0.6\pm0.2$. The ethane slope is consistent with the predictions of Gladstone et al. (1996), but that predicted for acetylene is larger than we observe. This disagreement is best explained by an overestimation of the acetylene production rate compared to that of ethane in the Gladstone et al. (1996) model. At 15.8 micron, methylacetylene is detected for the first time at low jovian latitudes, and a stratospheric column density of $(1.5\pm0.4)\times10^{15}$ molecule.cm-2 is inferred. We also derive an upper limit for the diacetylene column density of $7\times10^{13}$ molecule.cm-2.

The impact of the vapor pressure difference between deuterated and nondeuterated condensing molecules in planetary atmospheres is quantitatively assessed. This difference results in a loss of deuterium in the vapor phase above the condensation level. In Titan, Uranus and Neptune, the effect on CH3D is too subtle to alter current D/H ratio determinations. In Mars, the effect can induce a large depletion of HDO, starting about one scale height above the condensation level. Although the current infrared measurements of the D/H ratio appear to be almost unaffected, the intensity of disk-averaged millimetric HDO lines can be modified by about 10%. The effect is much stronger in limb sounding, and can be easily detected from orbiter observations.

We present the results of the Infrared Space Observatory Short Wavelength Spectrometer (ISO-SWS) observations of Jupiter related to ammonia. We focus on two spectral regions; the first one (the 10-micron region), probes atmospheric levels between 1 and 0.2 bar, while the second one (the 5-micron window), sounds the atmosphere between 8 and 2 bars. The two spectral windows cannot be fitted with the same ammonia vertical distribution. From the 10-micron region we infer an ammonia distribution of about half the saturation profile above the 1-bar level, where the N/H ratio is roughly solar. A totally different picture is derived from the 5-micron window, where we determine an upper limit of 3.7E-5 at 1 bar and find an increasing NH3 abundance at least down to 4 bar. This profile is similar to the one measured by the Galileo probe (Folkner et al. 1998). The discrepancy between the two spectral regions most likely arises from the spatial heterogeneity of Jupiter, the 5-micron window sounding dry areas unveiled by a locally thin cloud cover (the 5-micron hot spots), and the 10-micron region probing the mean jovian atmosphere above 1 bar. The 15NH3 mixing ratio is measured around 400 mbar from nu_2 band absorptions in the 10-micron region. We find the atmosphere of Jupiter highly depleted in 15N at this pressure level (15N/14N)_{Jupiter}=(1.9^{+0.9}_{-1.))E-3, while (15N/14N)_{earth}=3.68E-3. It is not clear whether this depletion reveals the global jovian (15N/14N) ratio. Instead an isotopic fractionation process, taking place during the ammonia cloud condensation, is indicated as a possible mechanism. A fractionation coefficient alpha higher than 1.08 would explain the observed isotopic ratio, but the lack of laboratory data does not allow us to decide unambiguously on the origin of the observed low (15N/14N) ratio.