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.

Uranus and Neptune, and their diverse satellite and ring systems, represent the least explored environments of our Solar System, and yet may provide the archetype for the most common outcome of planetary formation throughout our galaxy. Ice Giants will be the last remaining class of Solar System planet to have a dedicated orbital explorer, and international efforts are under way to realise such an ambitious mission in the coming decades. In 2019, the European Space Agency released a call for scientific themes for its strategic science planning process for the 2030s and 2040s, known as Voyage 2050. We used this opportunity to review our present-day knowledge of the Uranus and Neptune systems, producing a revised and updated set of scientific questions and motivations for their exploration. This review article describes how such a mission could explore their origins, ice-rich interiors, dynamic atmospheres, unique magnetospheres, and myriad icy satellites, to address questions at the heart of modern planetary science. These two worlds are superb examples of how planets with shared origins can exhibit remarkably different evolutionary paths: Neptune as the archetype for Ice Giants, whereas Uranus may be atypical. Exploring Uranus' natural satellites and Neptune's captured moon Triton could reveal how Ocean Worlds form and remain active, redefining the extent of the habitable zone in our Solar System. For these reasons and more, we advocate that an Ice Giant System explorer should become a strategic cornerstone mission within ESA's Voyage 2050 programme, in partnership with international collaborators, and targeting launch opportunities in the early 2030s.

We seek to compile a uniform set of basic capabilities and needs to maximize the yield of Solar System science with future Astrophysics assets while allowing those assets to achieve their Astrophysics priorities. Within considerations of cost and complexity, inclusion of capabilities that make a particular platform useable to planetary science provide a critical advantage over platforms lacking such capabilities.

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.

Remote sensing observations suffer significant limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. This impacts our knowledge of the formation of these planets and the physics of their atmospheres. A remarkable example of the superiority of in situ probe measurements was 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 were only made available through in situ measurements by the Galileo probe. Here we describe the main scientific goals to be addressed by the future in situ exploration of Saturn, Uranus, and Neptune, placing the Galileo probe exploration of Jupiter in a broader context. An atmospheric entry probe targeting the 10-bar level would yield insight into two broad themes: i) the formation history of the giant planets and 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. An atmospheric probe could represent a significant ESA contribution to a future NASA New Frontiers or flagship mission to be launched toward Saturn, Uranus, and/or Neptune.

We present results of simulations with the Max Planck Institute general circulation model (MPI-MGCM) implementing a hydrological cycle scheme. The simulations reveal a seasonal water "pump" mechanism responsible for the upward transport of water vapor. This mechanism occurs in high latitudes above 60$^\circ$ of the southern hemisphere at perihelion, when the upward branch of the meridional circulation is particularly strong. A combination of the mean vertical flux with variations induced by solar tides facilitates penetration of water across the "bottleneck" at approximately 60 km. The meridional circulation then transports water across the globe to the northern hemisphere. Since the intensity of the meridional cell is tightly controlled by airborne dust, the water abundance in the thermosphere strongly increases during dust storms.

The MIRO instrument's remote sensing capability is integral in constraining water density, temperature and velocity fields in the coma of 67P/Churyumov-Gersimenko. Our aim is to quantify the contribution to the water density from all facets inside and outside the field-of-view (FOV) of MIRO, in both, nadir and limb geometries. This information is crucial for understanding the MIRO derived coma production rates and their relation to the nucleus characteristics, and inherent spatial resolution of the data. This study relies on a detailed 3D nucleus shape model, illumination conditions and the pointing information of the viewing geometry. With these parameters, we can evaluate the relative contribution of water density originating from facets directly inside the MIRO beam as well outside of the beam as a function of distance along the MIRO line-of-sight. We also calculate the ratio of in-beam versus out-of-beam number density. We demonstrate that despite the rather small MIRO field-of-view there is only a small fraction of molecules that originate from facets within the MIRO beam. This is true for nadir, but also translated into the limb observing geometry. The MIRO instrument cannot discriminate active from non-active regions directly from observations. This study also suggests that the beam averaged solar incidence angle, local time and mean normal vectors are not necessary related to molecules within the MIRO beam. These results also illustrate why the 1D spherical Haser model can be applied with relative success for analyzing the MIRO data (and generally any Rosetta measurements). The future possibilities of constraining gas activity distribution on the surface should use 3D codes extracting information from the MIRO spectral line shapes which contain additional information. The presented results are in fact applicable to all relevant instruments onboard Rosetta.

Our aim is to investigate early activity (July 2014) of 67P/CG with 3D coma and radiative transfer modeling of MIRO measurements, accounting for nucleus shape, illumination, and orientation of the comet. We investigate MIRO line shape information for spatial distribution of water activity on the nucleus during the onset of activity. During this period we show that MIRO line shape have enough information to clearly isolate contribution from Hapi and Inhotep independently, and compare it to the nominal case of activity from the entire illuminated surface. We also demonstrate that spectral line shapes differ from the 1D model for different viewing geometries and coma conditions relevant to this study. Specifically, line shapes are somewhat sensitive to the location of the terminator in the coma. At last, fitting the MIRO observations we show that the Imhotep region (possible extended source of H$_{2}$O due to CO$_{2}$ activity) contributes only a small fraction of the total number of water molecules into MIRO beam in the early activity. On the other hand, a strong enhancement of water activity from the Hapi region seems required to fit the MIRO line shapes. This is consistent with earlier Rosetta results. Nevertheless, within the assumption of our coma and surface boundary conditions, we cannot get a reasonable fit to all MIRO mapping observations in July 2014, which may illustrate that a more sophisticated coma model or more accurate temperature/velocity distribution is needed.

Aims. We use four observational data sets, mainly from the Rosetta mission, to constrain the activity pattern of the nucleus of comet 67P/Churyumov-Gerasimenko. Methods. We develop a numerical model that computes the production rate and non-gravitational acceleration of the nucleus of comet 67P as a function of time, taking into account its complex shape with a shape model reconstructed from OSIRIS imagery. We use this model to fit three observational data sets: the trajectory data from flight dynamics; the rotation state, as reconstructed from OSIRIS imagery; and the water production measurements from ROSINA, of 67P. The two key parameters of our model, adjusted to fit the three data sets all together, are the activity pattern and the momentum transfer efficiency (i.e., the so-called "$\eta$ parameter" of the non-gravitational forces). Results. We find an activity pattern able to successfully reproduce the three data sets simultaneously. The fitted activity pattern exhibits two main features: a higher effective active fraction in two southern super-regions ($\sim 10$~\%) outside perihelion compared to the northern ones ($< 4$~\%), and a drastic rise of the effective active fraction of the southern regions ($\sim 25-35$~\%) around perihelion. We interpret the time-varying southern effective active fraction by cyclic formation and removal of a dust mantle in these regions. Our analysis supports moderate values of the momentum transfer coefficient $\eta$ in the range $0.6-0.7$; values $\eta\leq0.5$ or $\eta\geq0.8$ degrade significantly the fit to the three data sets. Our conclusions reinforce the idea that seasonal effects linked to the orientation of the spin axis play a key role in the formation and evolution of dust mantles, and in turn largely control the temporal variations of the gas flux.

Carbon dioxide (CO$_2$) ice clouds have been routinely observed in the middle atmosphere of Mars. However, there are still uncertainties concerning physical mechanisms that control their altitude, geographical, and seasonal distributions. Using the Max Planck Institute Martian General Circulation Model (MPI-MGCM), incorporating a state-of-the-art whole atmosphere subgrid-scale gravity wave parameterization [Yi\u{g}it et al., 2008], we demonstrate that internal gravity waves generated by lower atmospheric weather processes have wide reaching impact on the Martian climate. Globally, GWs cool the upper atmosphere of Mars by $\sim$10 % and facilitate high-altitude CO$_2$ ice cloud formation. CO$_2$ ice cloud seasonal variations in the mesosphere and the mesopause region appreciably coincide with the spatio-temporal variations of GW effects, providing insight into the observed distribution of clouds. Our results suggest that GW propagation and dissipation constitute a necessary physical mechanism for CO$_2$ ice cloud formation in the Martian upper atmosphere during all seasons.

We directly measure twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimate the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around one Pa or less (one to five Pa, when scaled to a metre length). The presence of eroded material at the base of most overhangs, as well as the observed collapse of two features and implied previous collapse of another, suggests that they are prone to failure and that true material strengths are close to these lower limits (although we only consider static stresses and not dynamic stress from, for example, cometary activity). Thus, a tensile strength of a few pascals is a good approximation for the tensile strength of 67P's nucleus material, which is in agreement with previous work. We find no particular trends in overhang properties with size, over the $\sim10-100$ m range studied here, or location on the nucleus. There are no obvious differences, in terms of strength, height or evidence of collapse, between the populations of overhangs on the two cometary lobes, suggesting that 67P is relatively homogenous in terms of tensile strength. Low material strengths are supportive of cometary formation as a primordial rubble pile or by collisional fragmentation of a small (tens of km) body.

We simulate the stresses induced by temperature changes in a putative hard layer near the surface of comet 67P/Churyumov--Gerasimenko with a thermo-viscoelastic model. Such a layer could be formed by the recondensation or sintering of water ice (and dust grains), as suggested by laboratory experiments and computer simulations, and would explain the high compressive strength encountered by experiments on board the Philae lander. Changes in temperature from seasonal insolation variation penetrate into the comet's surface to depths controlled by the thermal inertia, causing the material to expand and contract. Modelling this with a Maxwellian viscoelastic response on a spherical nucleus, we show that a hard, icy layer with similar properties to Martian permafrost will experience high stresses: up to tens of MPa, which exceed its material strength (a few MPa), down to depths of centimetres to a metre. The stress distribution with latitude is confirmed qualitatively when taking into account the comet's complex shape but neglecting thermal inertia. Stress is found to be comparable to the material strength everywhere for sufficient thermal inertia ($\gtrsim50$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) and ice content ($\gtrsim 45\%$ at the equator). In this case, stresses penetrate to a typical depth of $\sim0.25$ m, consistent with the detection of metre-scale thermal contraction crack polygons all over the comet. Thermal fracturing may be an important erosion process on cometary surfaces which breaks down material and weakens cliffs.

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.

Thermochemical models have been used in the past to constrain the deep oxygen abundance in the gas and ice giant planets from tropospheric CO spectroscopic measurements. Knowing the oxygen abundance of these planets is a key to better understand their formation. These models have widely used dry and/or moist adiabats to extrapolate temperatures from the measured values in the upper troposphere down to the level where the thermochemical equilibrium between H$_2$O and CO is established. The mean molecular mass gradient produced by the condensation of H$_2$O stabilizes the atmosphere against convection and results in a vertical thermal profile and H$_2$O distribution that departs significantly from previous estimates. We revisit O/H estimates using an atmospheric structure that accounts for the inhibition of the convection by condensation. We use a thermochemical network and the latest observations of CO in Uranus and Neptune to calculate the internal oxygen enrichment required to satisfy both these new estimates of the thermal profile and the observations. We also present the current limitations of such modeling.

The European Far-Infrared (FIR) Space Roadmap focuses on fundamental, yet still unresolved, astrophysical questions that can only be answered through a far-infrared space mission and gives an overview of the technology required to answer them. The document discusses topics ranging from Solar System and Planet Formation, Our Galaxy and nearby Galaxies and Distant Galaxies and Galaxy Evolution. The FIR Roadmap was open to comments from the wider astronomical community following a presentation during EWASS 2016.

We present a model of cometary dust capable of simulating the dynamics within the first few tens of km of the comet surface. Recent measurements by the GIADA and COSIMA instruments on Rosetta show that the nucleus emits fluffy dust particles with porosities above 50% and sizes up to at least mm (Schulz et al. 2015, Rotundi et al. 2015, Fulle et al. 2015). Retrieval of the physical properties of these particles requires a model of the effective forces governing their dynamics. Here, we present a model capable of simulating realistic, large and porous particles using hierarchical aggregates, which shows previous extrapolations to be inadequate. The main strengths of our approach are that we can simulate very large (mm-scale) non-spherical agglomerates and can accurately determine their 1) effective cross-section and ratio of cross-section to mass, 2) gas drag coefficient, and 3) light scattering properties. In practical terms, we find that a more detailed treatment of the dust structure results in 3-5 times higher velocities for large dust particles in the inner coma than previously estimated using spherical particles of the same mass. We apply our model to the dynamics of dust in the vicinity of the nucleus of comet 67P and successfully reproduce the dust speeds reported early on when the comet was roughly 3.5 AU from the Sun. At this stage, we employ a simple spherical comet nucleus, we model activity as constant velocity gas expansion from a uniformly active surface, and use Mie scattering. We discuss pathways to improve on these simplifications in the future.

Global characteristics of the small-scale gravity wave (GW) field in the Martian atmosphere obtained from a high-resolution general circulation model (GCM) are presented for the first time. The simulated GW-induced temperature variances are in a good agreement with available radio occultation data in the lower atmosphere between 10 and 30 km. The model reveals a latitudinal asymmetry with stronger wave generation in the winter hemisphere, and two distinctive sources of GWs: mountainous regions and the meandering winter polar jet. Orographic GWs are filtered while propagating upward, and the mesosphere is primarily dominated by harmonics with faster horizontal phase velocities. Wave fluxes are directed mainly against the local wind. GW dissipation in the upper mesosphere generates body forces of tens of m~s$^{-1}$~sol$^{-1}$, which tend to close the simulated jets. The results represent a realistic surrogate for missing observations, which can be used for constraining GW parameterizations and validating GCM simulations.

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 the first general circulation model simulations that quantify and reproduce patches of extremely cold air required for CO$_2$ condensation and cloud formation in the Martian mesosphere. They are created by subgrid-scale gravity waves (GWs) accounted for in the model with the interactively implemented spectral parameterization. Distributions of GW-induced temperature fluctuations and occurrences of supersaturation conditions are in a good agreement with observations of high-altitude CO$_2$ ice clouds. Our study confirms the key role of GWs in facilitating CO$_2$ cloud formation, discusses their tidal modulation, and predicts clouds at altitudes higher than have been observed to date.

Observations show that the lower thermosphere of Mars ($\sim$100--140 km) is up to 40 K colder than the current general circulation models (GCMs) can reproduce. Possible candidates for physical processes missing in the models are larger abundances of atomic oxygen facilitating stronger CO$_2$ radiative cooling, and thermal effects of gravity waves. Using two state-of-the-art Martian GCMs, the Laboratoire de M\'et\'eorologie Dynamique and Max Planck Institute models that self-consistently cover the atmosphere from the surface to the thermosphere, these physical mechanisms are investigated. Simulations demonstrate that the CO$_2$ radiative cooling with a sufficiently large atomic oxygen abundance, and the gravity wave-induced cooling can alone result in up to 40 K colder temperature in the lower thermosphere. Accounting for both mechanisms produce stronger cooling at high latitudes. However, radiative cooling effects peak above the mesopause, while gravity wave cooling rates continuously increase with height. Although both mechanisms act simultaneously, these peculiarities could help to further quantify their relative contributions from future observations.

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.

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard, however, is showing its age. Developed in the late 1970s, the FITS authors made a number of implementation choices that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not anticipate the challenges which we are facing today in astronomical computing. Difficulties we now face include, but are not limited to, addressing the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets, and capturing significantly more complex metadata and data relationships. There are members of the community today who find some or all of these limitations unworkable, and have decided to move ahead with storing data in other formats. If this fragmentation continues, we risk abandoning the advantages of broad interoperability, and ready archivability, that the FITS format provides for astronomy. In this paper we detail some selected important problems which exist within the FITS standard today. These problems may provide insight into deeper underlying issues which reside in the format and we provide a discussion of some lessons learned. It is not our intention here to prescribe specific remedies to these issues; rather, it is to call attention of the FITS and greater astronomical computing communities to these problems in the hope that it will spur action to address them.

We aimed to measure the H2O and dust production rates in C/2006 W3 (Christensen) with the Herschel Space Observatory at a heliocentric distance of ~ 5 AU. We have searched for emission in the H2O and NH3 ground-state rotational transitions at 557 GHz and 572 GHz, simultaneously, with HIFI onboard Herschel on UT 1.5 September 2010. Photometric observations of the dust coma in the 70 and 160 {\mu}m channels were acquired with the PACS instrument on UT 26.5 August 2010. A tentative 4-{\sigma} H2O line emission feature was found in the spectra obtained with the HIFI wide-band and high-resolution spectrometers, from which we derive a water production rate of $2.0(5) \times 10^{27}$ molec. s$^{-1}$. A 3-{\sigma} upper limit for the ammonia production rate of <$1.5 \times 10^{27}$ molec. s$^{-1}$ is obtained taking into account the contribution from all hyperfine components. The blueshift of the water line detected by HIFI suggests preferential emission from the subsolar point. However, it is also possible that water sublimation occurs in small ice-bearing grains that are emitted from an active region on the nucleus surface at a speed of ~ 0.2 km s$^{-1}$. The dust thermal emission was detected in the 70 and 160 {\mu}m filters, with a more extended emission in the blue channel. The dust production rates, obtained for a dust size distribution index that explains the fluxes at the photocenters of the PACS images, lie in the range 70-110 kg s$^{-1}$. Scaling the CO production rate measured post-perihelion at 3.20 and 3.32 AU, these values correspond to a dust-to-gas production rate ratio in the range 0.3-0.4. The dust production rates derived in August 2010 are roughly one order of magnitude lower than in September 2009, suggesting that the dust-to-gas production rate ratio remained approximately constant during the period when the activity became increasingly dominated by CO outgassing.

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.

HCl and HF are expected to be the main reservoirs of fluorine and chlorine wherever hydrogen is predominantly molecular. They are found to be strongly depleted in dense molecular clouds, suggesting freeze-out onto grains in such cold environments. We can then expect that HCl and HF were also the major carriers of Cl and F in the gas and icy phases of the outer solar nebula, and were incorporated into comets. We aimed to measure the HCl and HF abundances in cometary ices as they can provide insights on the halogen chemistry in the early solar nebula. We searched for the J(1-0) lines of HCl and HF at 626 and 1232 GHz, respectively, using the HIFI instrument on board the Herschel Space Observatory. HCl was searched for in comets 103P/Hartley 2 and C/2009 P1 (Garradd), whereas observations of HF were conducted in comet C/2009 P1. In addition, observations of H$_2$O and H$_2^{18}$O lines were performed in C/2009 P1 to measure the H$_2$O production rate. Three lines of CH$_3$OH were serendipitously observed in the HCl receiver setting. HCl is not detected, whereas a marginal (3.6-$\sigma$) detection of HF is obtained. The upper limits for the HCl abundance relative to water are 0.011% and 0.022%, for 103P and C/2009 P1, respectively, showing that HCl is depleted with respect to the solar Cl/O abundance by a factor more than 6$^{+6}_{-3}$ in 103P, where the error is related to the uncertainty in the chlorine solar abundance. The marginal HF detection obtained in C/2009 P1 corresponds to an HF abundance relative to water of (1.8$\pm$0.5) $\times$ 10$^{-4}$, which is approximately consistent with a solar photospheric F/O abundance. The observed depletion of HCl suggests that HCl was not the main reservoir of chlorine in the regions of the solar nebula where these comets formed. HF was possibly the main fluorine compound in the gas phase of the outer solar nebula.

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 line emission in the coma was measured in the comets C/2001 Q4 (NEAT) and C/2002 T7 (LINEAR), that were observed on five consecutive nights, 7-11 May 2004, at heliocentric distances of 1.0 and 0.7 AU, respectively, by means of high-resolution spectroscopy using the 10-m Submillimeter Telescope (SMT). We present a search for six parent- and product-volatile species (HCN, H2CO, CO, CS, CH3OH, and HNC) in both comets. Multiline observations of the CH3OH J = 5-4 series allow us to estimate the rotational temperature using the rotation diagram technique. We derive rotational temperatures of 54(9) K for C/2001 Q4 (NEAT) and 119(34) K for C/2002 T7 (LINEAR) that are roughly consistent with observations of other comets at similar distances from the Sun. The gas production rates of material are computed using a spherically symmetric molecular excitation code that includes collisions between neutrals and electrons. We find an HCN production rate of 2.96(5)e26 molec.s-1 for comet C/2001 Q4 (NEAT), corresponding to a mixing ratio with respect to H2O of 1.12(2)e-3. The mean HCN production rate during the observing period is 4.54(10)e26 molec.s-1 for comet C/2002 T7 (LINEAR), which gives a Q_HCN/Q_H2O mixing ratio of 1.51(3)e-3. With systematically lower mixing ratios in comet C/2001 Q4 (NEAT), production rate ratios of the observed species with respect to H2O lie within the typical ranges of dynamically new comets in both objects. We find a relative low abundance of CO in C/2001 Q4 (NEAT) compared to the observed range in other comets based on millimeter/submillimeter observations, and a significant upper limit on the CO production in C/2002 T7 (LINEAR) is derived. Depletion of CO suggests partial evaporation from the surface layers during previous visits to the outer Solar System and agrees with previous measurements of dynamically new comets.

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.

We present Herschel observations of water isotopologues in the atmosphere of the Jupiter-family comet 45P/Honda-Mrkos-Pajdusakova. No HDO emission is detected, with a 3 sigma upper limit of 2.0 10-4 for the D/H ratio. This value is consistent with the earlier Herschel measurement in the Jupiter-family comet 103P/Hartley 2. The canonical value of 3 10-4 measured pre-Herschel in a sample of Oort-cloud comets can be excluded at a 4.5 sigma level. The observations presented here further confirm that a diversity of D/H ratios exists in the comet population and emphasize the need for additional measurements with future ground-based facilities, such as CCAT, in the post-Herschel era.

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.

We aim to determine the production rates of several parent and product volatiles and the 12C/13C isotopic carbon ratio in the long-period comet C/2004 Q2 (Machholz), which is likely to originate from the Oort Cloud. The line emission from several molecules in the coma was measured with high signal-to-noise ratio in January 2005 at heliocentric distance of 1.2 AU by means of high-resolution spectroscopic observations using the Submillimeter Telescope (SMT). We have obtained production rates of several volatiles (CH3OH, HCN, H13CN, HNC, H2CO, CO and CS) by comparing the observed and simulated line-integrated intensities. Furthermore, multiline observations of the CH3OH (7-6) series allow us to estimate the rotational temperature using the rotation diagram technique. We find that the CH3OH population distribution of the levels sampled by these lines can be described by a rotational temperature of 40 \pm 3 K. Derived mixing ratios relative to hydrogen cyanide are CO/CH3OH/H2CO/CS/HNC/H13CN/HCN = 30.9/24.6/4.8/0.57/0.031/0.013/1 assuming a pointing offset of 8" due to the uncertain ephemeris at the time of the observations and the telescope pointing error. The measured relative molecular abundances in C/2004 Q2 (Machholz) are between low- to typical values of those obtained in Oort Cloud comets, suggesting that it has visited the inner solar system previously and undergone thermal processing. The HNC/HCN abundance ratio of ~3.1% is comparable to that found in other comets, accounting for the dependence on the heliocentric distance, and could possibly be explained by ion-molecule chemical processes in the low-temperature atmosphere. From a tentative H13CN detection, the measured value of 97 \pm 30 for the H12CN/H13CN isotopologue pair is consistent with a telluric value.

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.

Methods: We used the dual-band receiver GREAT on board the SOFIA airborne telescope to perform observations of the [C II] 158 {\mu}m fine-structure line at the postitions of two giant molecular clouds (GMC) in the center of IC 342 (GMCs C and E) and compared the spectra with corresponding ground-based data for low- and mid-J CO and [C I]. We performed model calculations assuming a clumpy photo-dissociation region (PDR) environment using the KOSMA-tau PDR model code to derive physical parameters of the local medium. Results: The [C II] 158 {\mu}m emission resembles the spectral signature of ground-based atomic and molecular lines, which indicates a common origin. The emission from GMC E can be decomposed into a cool, molecular component with weak far-ultraviolet (FUV) fields and low, mean densities of 103 cm^-3 and a strongly excited starburst/PDR region with higher densities of 104 cm^-3 and FUV intensities of 250-300 Draine fields. The emission from GMC C is consistent with gas densities of 5000 cm^-3, FUV intensities of a few Draine fields and total gas masses of 20\times10^6 M$_\odot$. Conclusions: The high spectral resolution of the GREAT receiver allowed us to decompose the [C II] emission of the GMC E into a strongly excited gas component resembling a PDR/starburst environment and a quieter, less excited gas component and to analyze the different components within a single beam individually.

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.

Context. The HII region/PDR/molecular cloud complex S106 is excited by a single O-star. The full extent of the warm and dense gas close to the star has not been mapped in spectrally resolved high-J CO or [CII] lines, so the kinematics of the warm, partially ionized gas, are unknown. Whether the prominent dark lane bisecting the hourglass-shaped nebula is due solely to the shadow cast by a small disk around the exciting star or also to extinction in high column foreground gas was an open question until now. Aims. To disentangle the morphology and kinematics of warm neutral and ionized gas close to the star, study their relation to the bulk of the molecular gas, and to investigate the nature of the dark lane. Methods. We use the heterodyne receiver GREAT on board SOFIA to observe velocity resolved spectral lines of [C II] and CO 11-10 in comparison with so far unpublished submm continuum data at 350 micron (SHARC-II) and complementary molecular line data. Results. The high angular and spectral resolution observations show a very complex morphology and kinematics of the inner S106 region, with many different components at different excitation conditions contributing to the observed emission. The [C II] lines are found to be bright and very broad, tracing high velocity gas close to the interface of molecular cloud and HII region. CO 11-10 emission is more confined, both spatially and in velocity, to the immediate surroundings of S106 IR showing the presence of warm, high density (clumpy) gas. Our high angular resolution submm continuum observations rule out the scenario where the dark lane separating the two lobes is due solely to the shadow cast by a small disk close to the star. The lane is clearly seen also as warm, high column density gas at the boundary of the molecular cloud and HII region.

The hydroxyl radical (OH) is found in various environments within the interstellar medium (ISM) of the Milky Way and external galaxies, mostly either in diffuse interstellar clouds or in the warm, dense environments of newly formed low-mass and high-mass stars, i.e, in the dense shells of compact and ultracompact HII regions (UCHIIRs). Until today, most studies of interstellar OH involved the molecule's radio wavelength hyperfine structure (hfs) transitions. These lines are generally not in LTE and either masing or over-cooling complicates their interpretation. In the past, observations of transitions between different rotational levels of OH, which are at far-infrared wavelengths, have suffered from limited spectral and angular resolution. Since these lines have critical densities many orders of magnitude higher than the radio wavelength ground state hfs lines and are emitted from levels with more than 100 K above the ground state, when observed in emission, they probe very dense and warm material. We probe the warm and dense molecular material surrounding the UCHIIR/OH maser sources W3(OH), G10.62-0.39 and NGC 7538 IRS1 by studying the $^2\Pi_{{1/2}}, J = {3/2} - {1/2}$ rotational transition of OH in emission and, toward the last source also the molecule's $^2\Pi_{3/2}, J = 5/2 - 3/2$ ground-state transition in absorption. We used the Stratospheric Observatory for Infrared Astronomy (SOFIA) to observe these OH lines, which are near 1.84 THz ($163 \mu$m) and 2.51 THz ($119.3 \mu$m). We clearly detect the OH lines, some of which are blended with each other. Employing non-LTE radiative transfer calculations we predict line intensities using models of a low OH abundance envelope versus a compact, high-abundance source corresponding to the origin of the radio OH lines.

We describe the design and construction of GREAT, the German REceiver for Astronomy at Terahertz frequencies operated on the Stratospheric Observatory for Infrared Astronomy (SOFIA). GREAT is a modular dual-color heterodyne instrument for highresolution far-infrared (FIR) spectroscopy. Selected for SOFIA's Early Science demonstration, the instrument has successfully performed three Short and more than a dozen Basic Science flights since first light was recorded on its April 1, 2011 commissioning flight. We report on the in-flight performance and operation of the receiver that - in various flight configurations, with three different detector channels - observed in several science-defined frequency windows between 1.25 and 2.5 THz. The receiver optics was verified to be diffraction-limited as designed, with nominal efficiencies; receiver sensitivities are state-of-the-art, with excellent system stability. The modular design allows for the continuous integration of latest technologies; we briefly discuss additional channels under development and ongoing improvements for Cycle 1 observations. GREAT is a principal investigator instrument, developed by a consortium of four German research institutes, available to the SOFIA users on a collaborative basis.

Context. Observing supernova remnants (SNRs) and modelling the shocks they are associated with is the best way to quantify the energy SNRs re-distribute back into the Interstellar Medium (ISM). Aims. We present comparisons of shock models with CO observations in the F knot of the W28 supernova remnant. These comparisons constitute a valuable tool to constrain both the shock characteristics and pre-shock conditions. Methods. New CO observations from the shocked regions with the APEX and SOFIA telescopes are presented and combined. The integrated intensities are compared to the outputs of a grid of models, which were combined from an MHD shock code that calculates the dynamical and chemical structure of these regions, and a radiative transfer module based on the 'large velocity gradient' (LVG) approximation. Results. We base our modelling method on the higher J CO transitions, which unambiguously trace the passage of a shock wave. We provide fits for the blue- and red-lobe components of the observed shocks. We find that only stationary, C-type shock models can reproduce the observed levels of CO emission. Our best models are found for a pre-shock density of 104 cm-3, with the magnetic field strength varying between 45 and 100 {\mu}G, and a higher shock velocity for the so-called blue shock (\sim25 km s-1) than for the red one (\sim20 km s-1). Our models also satisfactorily account for the pure rotational H2 emission that is observed with Spitzer.

We investigate the gas dynamics and the physical properties of photodissociation regions (PDRs) in IC1396A, which is an illuminated bright-rimmed globule with internal structures created by young stellar objects. Our mapping observations of the [CII] emission in IC1396A with GREAT onboard SOFIA revealed the detailed velocity structure of this region. We combined them with observations of the [CI] 3P_1 - 3P_0 and CO(4-3) emissions to study the dynamics of the different tracers and physical properties of the PDRs. The [CII] emission generally matches the IRAC 8 micron, which traces the polycyclic aromatic hydrocarbon (PAH) emissions. The CO(4-3) emission peaks inside the globule, and the [CI] emission is strong in outer regions, following the 8 micron emission to some degree, but its peak is different from that of [CII]. The [CII] emitting gas shows a clear velocity gradient within the globule, which is not significant in the [CI] and CO(4-3) emission. Some clumps that are prominent in [CII] emission appear to be blown away from the rim of the globule. The observed ratios of [CII]/[CI] and [CII]/CO(4-3) are compared to the KOSMA-tau PDR model, which indicates a density of 10^4-10^5 cm-3.

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.

The Jupiter-family comet 10P/Tempel 2 was observed during its 2010 return with the Herschel Space Observatory. We present here the observation of the (J, K) = (1, 0)-(0, 0) transition of ammonia at 572 GHz in this comet with the Heterodyne Instrument for the Far Infrared (HIFI) of Herschel. We also report on radio observations of other molecules (HCN, CH3OH, H2S and CS) obtained during the 1999 return of the comet with the CSO telescope and the JCMT, and during its 2010 return with the IRAM 30-m telescope. Molecular abundances relative to water are 0.09%, 1.8%, 0.4%, and 0.08% for HCN, CH3OH, H2S, and CS, respectively. An abundance of 0.5% for NH3 is obtained, which is similar to the values measured in other comets. The hyperfine structure of the ammonia line is resolved for the first time in an astronomical source. Strong anisotropy in the outgassing is present in all observations from 1999 to 2010 and is modelled to derive the production rates.

We present Herschel observations of the water molecule in the massive star-forming region W3 IRS5. The o-H17O 110-101, p-H18O 111-000, p-H2O 22 202-111, p-H2O 111-000, o-H2O 221-212, and o-H2O 212-101 lines, covering a frequency range from 552 up to 1669 GHz, have been detected at high spectral resolution with HIFI. The water lines in W3 IRS5 show well-defined high-velocity wings that indicate a clear contribution by outflows. Moreover, the systematically blue-shifted absorption in the H2O lines suggests expansion, presumably driven by the outflow. No infall signatures are detected. The p-H2O 111-000 and o-H2O 212-101 lines show absorption from the cold material (T ~ 10 K) in which the high-mass protostellar envelope is embedded. One-dimensional radiative transfer models are used to estimate water abundances and to further study the kinematics of the region. We show that the emission in the rare isotopologues comes directly from the inner parts of the envelope (T > 100 K) where water ices in the dust mantles evaporate and the gas-phase abundance increases. The resulting jump in the water abundance (with a constant inner abundance of 10^{-4}) is needed to reproduce the o-H17O 110-101 and p-H18O 111-000 spectra in our models. We estimate water abundances of 10^{-8} to 10^{-9} in the outer parts of the envelope (T < 100 K). The possibility of two protostellar objects contributing to the emission is discussed.

We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate.

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

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.

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the Herschel/HIFI instrument. No strong lines are detected down to 3sigma levels in 0.5 km/s channels of 4.2 mK for the 1_{10}--1_{01} line and 12.6 mK for the 1_{11}--0_{00} line. We report a very tentative detection, however, of the 1_{10}--1_{01} line in the Wide Band Spectrometer, with a strength of T_{mb}=2.7 mK, a width of 5.6 km/s and an integrated intensity of 16.0 mK km/s. The latter constitutes a 6sigma detection. Regardless of the reality of this tentative detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane.

We report a detection of the fundamental rotational transition of hydrogen fluoride in absorption towards Orion KL using Herschel/HIFI. After the removal of contaminating features associated with common molecules ("weeds"), the HF spectrum shows a P-Cygni profile, with weak redshifted emission and strong blue-shifted absorption, associated with the low-velocity molecular outflow. We derive an estimate of 2.9 x 10^13 cm^-2 for the HF column density responsible for the broad absorption component. Using our best estimate of the H2 column density within the low-velocity molecular outflow, we obtain a lower limit of ~1.6 x 10^-10 for the HF abundance relative to hydrogen nuclei, corresponding to 0.6% of the solar abundance of fluorine. This value is close to that inferred from previous ISO observations of HF J=2--1 absorption towards Sgr B2, but is in sharp contrast to the lower limit of 6 x 10^-9 derived by Neufeld et al. (2010) for cold, foreground clouds on the line of sight towards G10.6-0.4.

We present initial results from the Herschel GT key program: Herschel observations of EXtra-Ordinary Sources (HEXOS) and outline the promise and potential of spectral surveys with Herschel/HIFI. The HIFI instrument offers unprecedented sensitivity, as well as continuous spectral coverage across the gaps imposed by the atmosphere, opening up a largely unexplored wavelength regime to high-resolution spectroscopy. We show the spectrum of Orion KL between 480 and 560 GHz and from 1.06 to 1.115 THz. From these data, we confirm that HIFI separately measures the dust continuum and spectrally resolves emission lines in Orion KL. Based on this capability we demonstrate that the line contribution to the broad-band continuum in this molecule-rich source is ~20-40% below 1 THz and declines to a few percent at higher frequencies. We also tentatively identify multiple transitions of HD18O in the spectra. The first detection of this rare isotopologue in the interstellar medium suggests that HDO emission is optically thick in the Orion hot core with HDO/H2O ~ 0.02. We discuss the implications of this detection for the water D/H ratio in hot cores.

We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds show `spiral arm' cloud characteristics, while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr~B2 envelope. The observations show that each `spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have internal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H$_2$ column densities found at lower $A_V$ by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which probably means that CH is depleted in the denser cores of these clouds.

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.

Observations of the water inventory as well as other chemically important species on Jupiter will be performed in the frame of the guaranteed time key project of the Herschel Space Observatory entitled "Water and related chemistry in the Solar system". Among other onboard instruments, PACS (Photodetector Array Camera and Spectrometer) will provide new data of the spectral atlas in a wide region covering the far-infrared and submillimetre domains, with an improved spectral resolution and a higher sensitivity compared to previous observations carried out by Cassini/CIRS (Composite InfraRed Spectrometer) and by ISO (Infrared Space Observatory). In order to optimise the observational plan and to prepare for the data analysis, we have simulated the expected spectra of PACS Jupiter observations. Our simulation shows that PACS will promisingly detect several H2O emission lines. As PACS is capable of spatially resolving the Jovian disk, we will be able to discern the external oxygen sources in the giant planets by exploring the horizontal distribution of water. In addition to H2O lines, some absorption lines due to tropospheric CH4, HD, PH3 and NH3 lines will be observed with PACS. Furthermore, owing to the high sensitivity of the instrument, the current upper limit on the abundance of hydrogen halides such as HCl will be also improved.

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.

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.

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

The origin of the atmosphere of the largest moon of Saturn, Titan, is poorly understood and its chemistry is rather complicated. Ground-based millimeter/sub-millimeter heterodyne spectroscopy resolves line shapes sufficiently to determine information in Titan's atmospheric composition (on vertical profiles and isotopic ratios). We test the capabilities of the Swedish Heterodyne Facility Instrument (SHFI), Receiver APEX-1, together with the Atacama Pathfinder EXperiment APEX 12-m telescope for Titan's atmospheric observations. In particular we present sub-millimeter observations of the CO(2-1) and HCN(3-2) lines of the Titan stratosphere with APEX, and with SHFI taken during the Science Verification (SV) instrument phase on March and June 2008. With the help of appropriate radiative transfer calculations we investigate the possibility to constrain the chemical concentrations and optimize the performance of the APEX-1 instrument for inferring vertical profiles of molecular components of the atmosphere of Titan.

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.

We present submillimeter observations of 12CO J=3-2 and J=2-1, and 13CO J = 2-1 lines of the Venusian mesosphere and lower thermosphere with the Heinrich Hertz Submillimeter Telescope (HHSMT) taken around the second MESSENGER flyby of Venus on 5 June 2007. The observations cover a range of Venus solar elongations with different fractional disk illuminations. Preliminary results like temperature and CO abundance profiles are presented.These data are part of a coordinated observational campaign in support of the ESA Venus Express mission. Furthermore, this study attempts to contribute to cross-calibrate space- and ground-based observations, to constrain radiative transfer and retrieval algorithms for planetary atmospheres, and to a more thorough understanding of the global patters of circulation of the Venusian atmosphere.

We report vertical thermal structure and wind velocities in the Venusian mesosphere retrieved from carbon monoxide (12CO J=2-1 and 13CO J=2-1) spectral line observations obtained with the Heinrich Hertz Submillimeter Telescope (HHSMT). We observed the mesosphere of Venus from two days after the second Messenger flyby of Venus (on June 5 2007 at 23:10 UTC) during five days. Day-to-day and day-to-night temperature variations and short-term fluctuations of the mesospheric zonal flow were evident in our data. The extensive layer of warm air detected recently by SPICAV at 90 - to 100 km altitude is also detected in the temperature profiles reported here. These data were part of a coordinated ground-based Venus observational campaign in support of the ESA Venus Express mission. Furthermore, this study attempts to cross-calibrate space- and ground-based observations, to constrain radiative transfer and retrieval algorithms for planetary atmospheres, and to contribute to a more thorough understanding of the global patterns of circulation of the Venusian atmosphere.