In this work, we present results that simultaneously constrain the dust size distribution, dust-to-gas ratio, fraction of dust re-deposition, and total mass production rates for comet 67P/Churyumov-Gerasimenko. We use a 3D Direct Simulation Monte Carlo (DSMC) gas dynamics code to simulate the inner gas coma of the comet for the duration of the Rosetta mission. The gas model is constrained by ROSINA/COPS data. Further, we simulate for different epochs the inner dust coma using a 3D dust dynamics code including gas drag and the nucleus' gravity. Using advanced dust scattering properties these results are used to produce synthetic images that can be compared to the OSIRIS data set. These simulations allow us to constrain the properties of the dust coma and the total gas and dust production rates. We determined a total volatile mass loss of $(6.1 \pm 1.5) \cdot 10^9$~kg during the 2015 apparition. Further, we found that power-laws with $q=3.7^{+0.57}_{-0.078}$ are consistent with the data. This results in a total of $5.1^{+6.0}_{-4.9}\cdot10^9$~kg of dust being ejected from the nucleus surface, of which $4.4^{+4.9}_{-4.2}\cdot10^9$~kg escape to space and $6.8^{+11}_{-6.8}\cdot10^8$~kg (or an equivalent of $14^{+22}_{-14}$~cm over the smooth regions) is re-deposited on the surface. This leads to a dust-to-gas ratio of $0.73^{+1.3}_{-0.70}$ for the escaping material and $0.84^{+1.6}_{-0.81}$ for the ejected material. We have further found that the smallest dust size must be strictly smaller than $\sim30\mu$m and nominally even smaller than $\sim12\mu$m.

The European Space Agency's (ESA) Rosetta mission has returned a vast data set of measurements of the inner gas coma of comet 67P/Churyumov-Gerasimenko. These measurements have been used by different groups to determine the distribution of the gas sources at the nucleus surface. The solutions that have been found differ from each other substantially and illustrate the degeneracy of this issue. It is the aim of this work to explore the limitations that current gas models have in linking the coma measurements to the surface. In particular, we discuss the sensitivity of Rosetta's ROSINA/COPS, VIRTIS, and MIRO instruments to differentiate between vastly different spatial distributions of the gas emission from the surface. We have applied a state of the art 3D DSMC gas dynamics code to simulate the inner gas coma of different models that vary in the fraction of the surface that contains ice and in different sizes of active patches. These different distributions result in jet interactions that differ in their dynamical behaviour. We have found that ROSINA/COPS measurements by themselves cannot detect the differences in our models. While ROSINA/COPS measurements are important to constrain the regional inhomogeneities of the gas emission, they can by themselves not determine the surface emission distribution of the gas sources to a spatial accuracy of better than a few hundred metres (400 m). Any solutions fitting the ROSINA/COPS measurements is hence fundamentally degenerate, be it through a forward or inverse model. Only other instruments with complementary measurements can potentially lift this degeneracy as we show here for VIRTIS and MIRO. Finally, as a by-product, we have explored the effect of our activity distributions on lateral flow at the surface that may be responsible for some of the observed aeolian features.

The measured nitrogen-to-carbon ratio in comets is lower than for the Sun, a discrepancy which could be alleviated if there is an unknown reservoir of nitrogen in comets. The nucleus of comet 67P/Churyumov-Gerasimenko exhibits an unidentified broad spectral reflectance feature around 3.2 micrometers, which is ubiquitous across its surface. On the basis of laboratory experiments, we attribute this absorption band to ammonium salts mixed with dust on the surface. The depth of the band indicates that semivolatile ammonium salts are a substantial reservoir of nitrogen in the comet, potentially dominating over refractory organic matter and more volatile species. Similar absorption features appear in the spectra of some asteroids, implying a compositional link between asteroids, comets, and the parent interstellar cloud.

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.

We present new 3.6 and 4.5 micron secondary eclipse measurements for five cool (less than approximately 1000 K) transiting gas giant planets: HAT-P-15b, HAT-P-17b, HAT-P-18b, HAT-P-26b, and WASP-69b. We detect eclipses in at least one bandpass for all planets except HAT-P-15b. We confirm and refine the orbital eccentricity of HAT-P-17b, which is also the only planet in our sample with a known outer companion. We compare our measured eclipse depths in these two bands, which are sensitive to the relative abundances of methane versus carbon monoxide and carbon dioxide, respectively, to predictions from 1D atmosphere models for each planet. For planets with hydrogen-dominated atmospheres and equilibrium temperatures cooler than approximately 1000 K, this ratio should vary as a function of both atmospheric metallicity and the carbon-to-oxygen ratio. For HAT-P-26b, our observations are in good agreement with the low atmospheric metallicity inferred from transmission spectroscopy. We find that all four of the planets with detected eclipses are best matched by models with relatively efficient circulation of energy to the nightside. We see no evidence for a solar-system-like correlation between planet mass and atmospheric metallicity, but instead identify a potential (1.9 sigma) correlation between the inferred methane/(carbon monoxide + carbon dioxide) ratio and stellar metallicity. Our ability to characterize this potential trend is limited by the relatively large uncertainties in the stellar metallicity values. Our observations provide a first look at the brightness of these planets at wavelengths accessible to the James Webb Space Telescope, which will be able to resolve individual methane, carbon monoxide, and carbon dioxide bands and provide much stronger constraints on their atmospheric compositions.

This white paper proposes that AMBITION, a Comet Nucleus Sample Return mission, be a cornerstone of ESA's Voyage 2050 programme. We summarise some of the most important questions still open in cometary science after the successes of the Rosetta mission, many of which require sample analysis using techniques that are only possible in laboratories on Earth. We then summarise measurements, instrumentation and mission scenarios that can address these questions, with a recommendation that ESA select an ambitious cryogenic sample return mission. Rendezvous missions to Main Belt comets and Centaurs are compelling cases for M-class missions, expanding our knowledge by exploring new classes of comets. AMBITION would engage a wide community, drawing expertise from a vast range of disciplines within planetary science and astrophysics. With AMBITION, Europe will continue its leadership in the exploration of the most primitive Solar System bodies.

We examine the co-evolution of galaxies and supermassive black holes in the Simba cosmological hydrodynamic simulation. Simba grows black holes via gravitational torque-limited accretion from cold gas and Bondi accretion from hot gas, while feedback from black holes is modeled in radiative and jet modes depending on the Eddington ratio ($f_{Edd}$). Simba shows generally good agreement with local studies of black hole properties, such as the black hole mass--stellar velocity dispersion ($M_{BH}-\sigma$) relation, 2 the black hole accretion rate vs. star formation rate (BHAR--SFR), and the black hole mass function. $M_{BH}-\sigma$ evolves such that galaxies at a given $M_{BH}$ have higher $\sigma$ at higher redshift, consistent with no evolution in $M_{BH}-M_*$. For $M_{BH}<\sim 10^8 M_\odot$, $f_{Edd}$ is anti-correlated with $M_{BH}$ since the BHAR is approximately independent of $M_{BH}$, while at higher masses $f_{Edd}-M_{BH}$ flattens and has a larger scatter. BHAR vs. SFR is invariant with redshift, but $f_{Edd}$ drops steadily with time at a given $M_{BH}$, such that all but the most massive black holes are accreting in a radiatively efficient mode at $z>\sim 2$. The black hole mass function amplitude decreases with redshift and is locally dominated by quiescent galaxies for $M_{BH}>10^{8}M_{\odot}$, but for $z>\sim 1$ star forming galaxies dominate at all $M_{BH}$. The $z=0$ $f_{Edd}$ distribution is roughly lognormal with a peak at $f_{Edd}<\sim 0.01$ as observed, shifting to higher $f_{Edd}$ at higher redshifts. Finally, we study the dependence of black hole properties with \HI\ content and find that the correlation between gas content and star formation rate is modulated by black hole properties, such that higher SFR galaxies at a given gas content have smaller black holes with higher $f_{Edd}$

We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$\mu$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud.

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.

The polarimetry of the light scattered by planetary surfaces is a powerful tool to provide constraints on their microstructure. To improve the interpretation of polarimetric data from icy surfaces, we have developed the POLarimeter for ICE Samples (POLICES) complementing the measurement facilities of the Ice Laboratory at the University of Bern. The new setup uses a high precision Stokes polarimeter to measure the degree of polarization in the visible light scattered by surfaces at moderate phase angles (from 1.5 to 30{\deg}). We present the photometric and polarimetric phase curves measured on various surfaces made of pure water ice particles having well-controlled size and shape (spherical, crushed, frost). The results show how the amplitude and the shape of the negative polarization branch change with the particles sizes and the degree of metamorphism of the ice. We found that fresh frost formed by water condensation on cold surfaces has a phase curve characterized by resonances (Mie oscillations) indicating that frost embryos are transparent micrometer-sized particles with a narrow size distribution and spherical shape. Comparisons of these measurements with polarimetric observations of the icy satellites of the Solar System suggest that Europa is possibly covered by relatively coarser (~40-400 {\mu}m) and more sintered grains than Enceladus and Rhea, more likely covered by frost-like particles of few micrometers in average. The great sensitivity of polarization to grain size and degree of sintering makes it an ideal tool to detect hints of ongoing processes on icy planetary surfaces, such as cryovolcanism.

Linear polarization maps of the Carina Nebula were obtained at 250, 350, and 500 $\mu$m during the 2012 flight of the BLASTPol balloon-borne telescope. These measurements are combined with Planck 850 $\mu$m data in order to produce a submillimeter spectrum of the polarization fraction of the dust emission, averaged over the cloud. This spectrum is flat to within $\pm$15% (relative to the 350 $\mu$m polarization fraction). In particular, there is no evidence for a pronounced minimum of the spectrum near 350 $\mu$m, as suggested by previous ground-based measurements of other molecular clouds. This result of a flat polarization spectrum in Carina is consistent with recently-published BLASTPol measurements of the Vela C molecular cloud, and also agrees with a published model for an externally-illuminated, dense molecular cloud by Bethell and collaborators. The shape of the spectrum in Carina does not show any dependence on the radiative environment of the dust, as quantified by the Planck-derived dust temperature or dust optical depth at 353 GHz.

The wavelength calibration of spectrographs is an essential but challenging task in many disciplines. Calibration is traditionally accomplished by imaging the spectrum of a light source containing features that are known to appear at certain wavelengths and mapping them to their location on the sensor. This is typically required in conjunction with each scientific observation to account for mechanical and optical variations of the instrument over time, which may span years for certain projects. The method presented here investigates the usage of color itself instead of spectral features to calibrate a spectrograph. The primary advantage of such a calibration is that any broad-spectrum light source such as the sky or an incandescent bulb is suitable. This method allows for calibration using the full optical pathway of the instrument instead of incorporating separate calibration equipment that may introduce errors. This paper focuses on the potential for color calibration in the field of radial velocity astronomy, in which instruments must be finely calibrated for long periods of time to detect tiny Doppler wavelength shifts. This method is not restricted to radial velocity, however, and may find application in any field requiring calibrated spectrometers such as sea water analysis, cellular biology, chemistry, atmospheric studies, and so on. This paper demonstrates that color sensors have the potential to provide calibration with greatly reduced complexity.

Combining stellar atmospheric parameters, such as effective temperature, surface gravity, and metallicity, with barycentric radial velocity data provides insight into the chemo-dynamics of the MilkyWay and our local Galactic environment. We analyze 3075 stars with spectroscopic data from the Sloan Digital Sky Survey III (SDSS-III) MARVELS radial velocity survey and present atmospheric parameters for 2343 dwarf stars using the spectral indices method, a modified version of the equivalent width method. We present barycentric radial velocities for 2610 stars with uncertainties of $\sim$0.3 km s$^{-1}$. We determine stellar ages using two independent methods and calculate ages for 2343 stars with a maximum-likelihood isochronal age-dating method and for 1694 stars with a Bayesian age-dating method. Using previously published parallax data we compute Galactic orbits and space velocities for 1995 stars to explore stellar populations based on kinematic and age parameters. This study combines good ages and exquisite velocities to explore local chemo-kinematics of the Milky Way, which complements many of the recent studies of giant stars with the APOGEE survey, and we find our results to be in agreement with current chemo-dynamical models of the Milky Way. Particularly, we find from our metallicity distributions and velocity-age relations of a kinematically-defined thin disk that the metal rich end has stars of all ages, even after we clean the sample of highly eccentric stars, suggesting that radial migration plays a key role in the metallicity scatter of the thin disk. All stellar parameters and kinematic data derived in this work are catalogued and published online in machine-readable form.

We propose to search for biosignatures in the spectra of reflected light from about 100 Earth-sized planets that are already known to be orbiting in their habitable zones (HZ). For a sample of G and K type hosts, most of these planets will be between 25 and 50 milli-arcsec (mas) from their host star and 1 billion to 10 billion times fainter. To separate the planet's image from that of its host star at the wavelength (763nm) of the oxygen biosignature we need a telescope with an aperture of 16 metres. Furthermore, the intensity of the light from the host star at the position in the image of the exoplanet must be suppressed otherwise the exoplanet will be lost in the glare. This presents huge technical challenges. The Earth's atmosphere is turbulent which makes it impossible to achieve the required contrast from the ground at 763nm. The telescope therefore needs to be in space and to fit the telescope in the rocket fairing it must be a factor of 4 or more times smaller when folded than when operational. To obtain spectroscopy of the planet's biosignature at 763nm we need to use an integral field spectrometer (IFS) with a field of view (FOV) of 1000 x 1000 milli-arcsec (mas) and a spectral resolution of 100. This is a device that simultaneously takes many pictures of the exoplanet each at a slightly different wavelength which are then recorded as a data cube with two spatial dimensions and one wavelength dimension. In every data cube wavelength slice, the background light from the host star at the location of the planet image must be minimised. This is achieved via a coronagraph which blocks the light from the host star and active/adaptive optics techniques which continuously maintain very high accuracy optical alignment to make the images as sharp as possible. These are the technical challenges to be addressed in a design study.

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.

We present a study of an impacting descent probe that increases the science return of spacecraft orbiting or passing an atmosphere-less planetary body of the solar system, such as the Galilean moons of Jupiter. The descent probe is a carry-on small spacecraft (< 100 kg), to be deployed by the mother spacecraft, that brings itself onto a collisional trajectory with the targeted planetary body in a simple manner. A possible science payload includes instruments for surface imaging, characterisation of the neutral exosphere, and magnetic field and plasma measurement near the target body down to very low-altitudes (~1 km), during the probe's fast (~km/s) descent to the surface until impact. The science goals and the concept of operation are discussed with particular reference to Europa, including options for flying through water plumes and after-impact retrieval of very-low altitude science data. All in all, it is demonstrated how the descent probe has the potential to provide a high science return to a mission at a low extra level of complexity, engineering effort, and risk. This study builds upon earlier studies for a Callisto Descent Probe (CDP) for the former Europa-Jupiter System Mission (EJSM) of ESA and NASA, and extends them with a detailed assessment of a descent probe designed to be an additional science payload for the NASA Europa Mission.

CASTAway is a mission concept to explore our Solar System's main asteroid belt. Asteroids and comets provide a window into the formation and evolution of our Solar System and the composition of these objects can be inferred from space-based remote sensing using spectroscopic techniques. Variations in composition across the asteroid populations provide a tracer for the dynamical evolution of the Solar System. The mission combines a long-range (point source) telescopic survey of over 10,000 objects, targeted close encounters with 10 to 20 asteroids and serendipitous searches to constrain the distribution of smaller (e.g. 10 m) size objects into a single concept. With a carefully targeted trajectory that loops through the asteroid belt, CASTAway would provide a comprehensive survey of the main belt at multiple scales. The scientific payload comprises a 50 cm diameter telescope that includes an integrated low-resolution (R = 30 to 100) spectrometer and visible context imager, a thermal (e.g. 6 to 16 microns) imager for use during the flybys, and modified star tracker cameras to detect small (approx. 10 m) asteroids. The CASTAway spacecraft and payload have high levels of technology readiness and are designed to fit within the programmatic and cost caps for a European Space Agency medium class mission, whilst delivering a significant increase in knowledge of our Solar System.

On 3 July 2016, several instruments on board ESA's Rosetta spacecraft detected signs of an outburst event on comet 67P, at a heliocentric distance of 3.32 AU from the sun, outbound from perihelion. We here report on the inferred properties of the ejected dust and the surface change at the site of the outburst. The activity coincided with the local sunrise and continued over a time interval of 14 - 68 minutes. It left a 10m-sized icy patch on the surface. The ejected material comprised refractory grains of several hundred microns in size, and sub-micron-sized water ice grains. The high dust mass production rate is incompatible with the free sublimation of crystalline water ice under solar illumination as the only acceleration process. Additional energy stored near the surface must have increased the gas density. We suggest a pressurized sub-surface gas reservoir, or the crystallization of amorphous water ice as possible causes.

The M dwarf GJ 436 hosts a transiting warm Neptune known to experience atmospheric escape. Previous observations revealed the presence of a giant hydrogen exosphere transiting the star for more than 5 h, and absorbing up to 56% of the flux in the blue wing of the stellar Lyman-{\alpha} line of neutral hydrogen (H i Ly{\alpha}). The unexpected size of this comet-like exosphere prevented observing the full transit of its tail. In this Letter, we present new Ly{\alpha} observations of GJ 436 obtained with the Space Telescope Imaging Spectrograph (STIS) instrument onboard the Hubble Space Telescope. The stability of the Ly{\alpha} line over six years allowed us to combine these new observations with archival data sets, substantially expanding the coverage of the exospheric transit. Hydrogen atoms in the tail of the exospheric cloud keep occulting the star for 10-25 h after the transit of the planet, remarkably confirming a previous prediction based on 3D numerical simulations with the EVaporating Exoplanet code (EVE). This result strengthens the interpretation that the exosphere of GJ 436b is shaped by both radiative braking and charge exchanges with the stellar wind. We further report flux decreases of 15 +/- 2% and 47 +/- 10% in the red wing of the Ly{\alpha} line and in the line of ionised silicon (Si iii). Despite some temporal variability possibly linked with stellar activity, these two signals occur during the exospheric transit and could be of planetary origin. Follow-up observations will be required to assess the possibility that the redshifted Ly{\alpha} and Si iii absorption signatures arise from interactions between the exospheric flow and the magnetic field of the star.

We present a statistical analysis of the distribution of large scale topographic features on comet 67P/Churyumov-Gerasimenko. We observe that the cumulative cliff height distribution across the surface follows a power law with a slope equal to -1.69 +- 0.02. When this distribution is studied independently for each region, we find a good correlation between the slope of the power law and the orbital erosion rate of the surface. For instance, the northern hemisphere topography is dominated by structures on the 100~m scale while the southern hemisphere topography, illuminated at perihelion, is dominated by 10~m scale terrain features. Our study suggest that the current size of a cliff is controlled not only by material cohesion but by the dominant erosional process in each region. This observation can be generalized to other comets, where we argue that primitive nuclei are characterized by the presence of large cliffs with a cumulative height power index equal to or above -1.5, while older, eroded cometary surfaces have a power index equal to or below -2.3. In effect, our model shows that a measure of the topography provides a quantitative assessment of a comet's erosional history, i.e. its evolutionary age.

We collect observational evidence that supports the scheme of mass transfer on the nucleus of comet 67P/Churyumov-Gerasimenko. The obliquity of the rotation axis of 67P causes strong seasonal variations. During perihelion the southern hemisphere is four times more active than the north. Northern territories are widely covered by granular material that indicates back fall originating from the active south. Decimetre sized chunks contain water ice and their trajectories are influenced by an anti-solar force instigated by sublimation. OSIRIS observations suggest that up to 20 % of the particles directly return to the nucleus surface taking several hours of travel time. The back fall covered northern areas are active if illuminated but produce mainly water vapour. The decimetre chunks from the nucleus surface are too small to contain more volatile compounds such as CO 2 or CO. This causes a north-south dichotomy of the composition measurements in the coma. Active particles are trapped in the gravitational minimum of Hapi during northern winter. They are "shock frozen" and only reactivated when the comet approaches the sun after its aphelion passage. The insolation of the big cavity is enhanced by self-heating, i. e. reflection and IR radiation from the walls. This, together with the pristinity of the active back fall, explains the early observed activity of the Hapi region. Sobek may be a role model for the consolidated bottom of Hapi. Mass transfer in the case of 67P strongly influences the evolution of the nucleus and the interpretation of coma measurements.

The Southern hemisphere of the 67P/Churyumov-Gerasimenko comet has become visible from Rosetta only since March 2015. It was illuminated during the perihelion passage and therefore it contains the regions that experienced the strongest heating and erosion rate, thus exposing the subsurface most pristine material. In this work we investigate, thanks to the OSIRIS images, the geomorphology, the spectrophotometry and some transient events of two Southern hemisphere regions: Anhur and part of Bes. Bes is dominated by outcropping consolidated terrain covered with fine particle deposits, while Anhur appears strongly eroded with elongated canyon-like structures, scarp retreats, different kinds of deposits, and degraded sequences of strata indicating a pervasive layering. We discovered a new 140 m long and 10 m high scarp formed in the Anhur/Bes boundary during/after the perihelion passage, close to the area where exposed CO$_2$ and H$_2$O ices were previously detected. Several jets have been observed originating from these regions, including the strong perihelion outburst, an active pit, and a faint optically thick dust plume. We identify several areas with a relatively bluer slope (i.e. a lower spectral slope value) than their surroundings, indicating a surface composition enriched with some water ice. These spectrally bluer areas are observed especially in talus and gravitational accumulation deposits where freshly exposed material had fallen from nearby scarps and cliffs. The investigated regions become spectrally redder beyond 2 au outbound when the dust mantle became thicker, masking the underlying ice-rich layers.

Polarized emission from aligned dust is a crucial tool for studies of magnetism in the ISM and a troublesome contaminant for studies of CMB polarization. In each case, an understanding of the significance of the polarization signal requires well-calibrated physical models of dust grains. Despite decades of progress in theory and observation, polarized dust models remain largely underconstrained. During its 2012 flight, the balloon-borne telescope BLASTPol obtained simultaneous broad-band polarimetric maps of a translucent molecular cloud at 250, 350, and 500 microns. Combining these data with polarimetry from the Planck 850 micron band, we have produced a submillimeter polarization spectrum for a cloud of this type for the first time. We find the polarization degree to be largely constant across the four bands. This result introduces a new observable with the potential to place strong empirical constraints on ISM dust polarization models in a previously inaccessible density regime. Comparing with models by Draine and Fraisse (2009), our result disfavors two of their models for which all polarization arises due only to aligned silicate grains. By creating simple models for polarized emission in a translucent cloud, we verify that extinction within the cloud should have only a small effect on the polarization spectrum shape compared to the diffuse ISM. Thus we expect the measured polarization spectrum to be a valid check on diffuse ISM dust models. The general flatness of the observed polarization spectrum suggests a challenge to models where temperature and alignment degree are strongly correlated across major dust components.

There are many geometric calibration methods for "standard" cameras. These methods, however, cannot be used for the calibration of telescopes with large focal lengths and complex off-axis optics. Moreover, specialized calibration methods for the telescopes are scarce in literature. We describe the calibration method that we developed for the Colour and Stereo Surface Imaging System (CaSSIS) telescope, on board of the ExoMars Trace Gas Orbiter (TGO). Although our method is described in the context of CaSSIS, with camera-specific experiments, it is general and can be applied to other telescopes. We further encourage re-use of the proposed method by making our calibration code and data available on-line.

Images of the nucleus and the coma (gas and dust) of comet 67P/Churyumov- Gerasimenko have been acquired by the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras since March 2014 using both the Wide Angle Camera (WAC) and the Narrow Angle Camera (NAC). We use images from the NAC camera to study a bright outburst observed in the southern hemisphere on July 29, 2015. The high spatial resolution of the NAC is needed to localize the source point of the outburst on the surface of the nucleus. The heliocentric distance is 1.25 au and the spacecraft-comet distance is 186 km. Aiming to better understand the physics that led to the outgassing, we used the Direct Simulation Monte Carlo (DSMC) method to study the gas flow close to the nucleus and the dust trajectories. The goal is to understand the mechanisms producing the outburst. We reproduce the opening angle of the outburst in the model and constrain the outgassing ratio between the outburst source and the local region. The outburst is in fact a combination of both gas and dust, in which the active surface is approximately 10 times more active than the average rate found in the surrounding areas. We need a number of dust particles 7.83 $\times$ 10$^{11}$ - 6.90 $\times$ 10$^{15}$ (radius 1.97 - 185 {\mu}m), which corresponds to a mass of dust 220 - 21 $\times$ 10$^{3}$kg.

The ESA Rosetta spacecraft has been tracking its target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, in close vicinity for over two years. It hosts the OSIRIS instruments: the Optical, Spectroscopic, and Infrared Remote Imaging System composed of two cameras, see e.g. Keller et al. (2007). In some imaging sequences dedicated to observe dust particles in the comet's coma, the two cameras took images at the same time. The aim of this work is to use these simultaneous double camera observations to calculate the dust particles' distance to the spacecraft. As the two cameras are mounted on the spacecraft with an offset of 70 cm, the distance of particles observed by both cameras can be determined by a shift of the particles' apparent trails on the images. This paper presents first results of the ongoing work, introducing the distance determination method for the OSIRIS instrument and the analysis of an example particle. We note that this method works for particles in the range of about 500 m - 6000 m from the spacecraft.

Kepler-444 provides a unique opportunity to probe the atmospheric composition and evolution of a compact system of exoplanets smaller than the Earth. Five planets transit this bright K star at close orbital distances, but they are too small for their putative lower atmosphere to be probed at optical/infrared wavelengths. We used the Space Telescope Imaging Spectrograph instrument onboard the Hubble Space Telescope to search for the signature of the planet's upper atmospheres at six independent epochs in the Ly-$\alpha$ line. We detect significant flux variations during the transits of both Kepler-444e and f (~20%), and also at a time when none of the known planets was transiting (~40%). Variability in the transition region and corona of the host star might be the source of these variations. Yet, their amplitude over short time scales (~2-3 hours) is surprisingly strong for this old (11.2+-1.0Gyr) and apparently quiet main-sequence star. Alternatively, we show that the in-transits variations could be explained by absorption from neutral hydrogen exospheres trailing the two outer planets (Kepler-444e and f). They would have to contain substantial amounts of water to replenish such hydrogen exospheres, which would reveal them as the first confirmed ocean-planets. The out-of-transit variations, however, would require the presence of a yet-undetected Kepler-444g at larger orbital distance, casting doubt on the planetary origin scenario. Using HARPS-N observations in the sodium doublet, we derived the properties of two Interstellar Medium clouds along the line-of-sight toward Kepler-444. This allowed us to reconstruct the stellar Ly-$\alpha$ line profile and to estimate the XUV irradiation from the star, which would still allow for a moderate mass loss from the outer planets after 11.2Gyr. Follow-up of the system at XUV wavelengths will be required to assess this tantalizing possibility.

We present a large-scale combination of near-infrared (near-IR) interstellar polarization data from background starlight with polarized emission data at submillimeter (sub-mm) wavelengths for the Vela C molecular cloud. The near-IR data consist of more than 6700 detections probing a range of visual extinctions between $2$ and $20\,$mag in and around the cloud. The sub-mm data was collected in Antartica by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). This is the first direct combination of near-IR and sub-mm polarization data for a molecular cloud aimed at measuring the "polarization efficiency ratio" ($R_{\mathrm{eff}}$), a quantity that is expected to depend only on grain intrinsic physical properties. It is defined as $p_{500}/(p_{I}/\tau_{V})$, where $p_{500}$ and $p_{I}$ are polarization fractions at $500\,\mu$m and $I$-band, respectively, and $\tau_{V}$ is the optical depth. To ensure that the same column density of material is producing both polarization from emission and from extinction, we conducted a careful selection of near-background stars using 2MASS, $Herschel$ and $Planck$ data. This selection excludes objects contaminated by the Galactic diffuse background material as well as objects located in the foreground. Accounting for statistical and systematic uncertainties, we estimate an average $R_{\mathrm{eff}}$ value of $2.4\pm0.8$, which can be used to test the predictions of dust grain models designed for molecular clouds when such predictions become available. $R_{\mathrm{eff}}$ appears to be relatively flat as a function of the cloud depth for the range of visual extinctions probed.

We statistically evaluate the relative orientation between gas column density structures, inferred from Herschel submillimetre observations, and the magnetic field projected on the plane of sky, inferred from polarized thermal emission of Galactic dust observed by BLASTPol at 250, 350, and 500 micron, towards the Vela C molecular complex. First, we find very good agreement between the polarization orientations in the three wavelength-bands, suggesting that, at the considered common angular resolution of 3.0 arcminutes that corresponds to a physical scale of approximately 0.61 pc, the inferred magnetic field orientation is not significantly affected by temperature or dust grain alignment effects. Second, we find that the relative orientation between gas column density structures and the magnetic field changes progressively with increasing gas column density, from mostly parallel or having no preferred orientation at low column densities to mostly perpendicular at the highest column densities. This observation is in agreement with previous studies by the Planck collaboration towards more nearby molecular clouds. Finally, we find a correspondence between the trends in relative orientation and the shape of the column density probability distribution functions. In the sub-regions of Vela C dominated by one clear filamentary structure, or "ridges", we find a sharp transition from preferentially parallel or having no preferred relative orientation at low column densities to preferentially perpendicular at highest column densities. In the sub-regions of Vela C dominated by several filamentary structures with multiple orientations, or "nests", such a transition is also present, but it is clearly less sharp than in the ridge-like sub-regions. Both of these results suggest that the magnetic field is dynamically important for the formation of density structures in this region.

Planet searches using the radial velocity technique show a paucity of companions to solar-type stars within ~5 AU in the mass range of ~10 - 80 M$_{\text{Jup}}$. This deficit, known as the brown dwarf desert, currently has no conclusive explanation. New substellar companions in this region help asses the reality of the desert and provide insight to the formation and evolution of these objects. Here we present 10 new brown dwarf and two low-mass stellar companion candidates around solar-type stars from the Multi-object APO Radial-Velocity Exoplanet Large-Area Survey (MARVELS) of the Sloan Digital Sky Survey III (SDSS-III). These companions were selected from processed MARVELS data using the latest University of Florida Two Dimensional (UF2D) pipeline, which shows significant improvement and reduction of systematic errors over previous pipelines. The 10 brown dwarf companions range in mass from ~13 to 76 M$_{\text{Jup}}$ and have orbital radii of less than 1 AU. The two stellar companions have minimum masses of ~98 and 100 M$_{\text{Jup}}$. The host stars of the MARVELS brown dwarf sample have a mean metallicity of [Fe/H] = 0.03 $\pm$ 0.08 dex. Given our stellar sample we estimate the brown dwarf occurrence rate around solar-type stars with periods less than ~300 days to be ~0.56%.

The visible reflectance spectrum of many solar system bodies changes with changing viewing geometry for reasons not fully understood. It is often observed to redden (increasing spectral slope) with increasing solar phase angle, an effect known as phase reddening. Only once, in an observation of the Martian surface by the Viking 1 lander, was reddening observed up to a certain phase angle with bluing beyond, making the reflectance ratio as a function of phase angle shaped like an arch. However, in laboratory experiments this arch-shape is frequently encountered. To investigate this, we measured the bidirectional reflectance of particulate samples of several common rock types in the 400-1000 nm wavelength range and performed ray-tracing simulations. We confirm the occurrence of the arch for surfaces that are forward scattering, i.e. are composed of semi-transparent particles and are smooth on the scale of the particles, and for which the reflectance increases from the lower to the higher wavelength in the reflectance ratio. The arch shape is reproduced by the simulations, which assume a smooth surface. However, surface roughness on the scale of the particles, such as the Hapke and van Horn (1963) fairy castles that can spontaneously form when sprinkling a fine powder, leads to monotonic reddening. A further consequence of this form of microscopic roughness (being indistinct without the use of a microscope) is a flattening of the disk function at visible wavelengths, i.e. Lommel-Seeliger-type scattering. The experiments further reveal monotonic reddening for reflectance ratios at near-IR wavelengths. The simulations fail to reproduce this particular reddening, and we suspect that it results from roughness on the surface of the particles. Given that the regolith of atmosphereless solar system bodies is composed of small particles, our results indicate that the prevalence of monotonic (...)

We use a laboratory facility to study the sputtering properties of centimeter-thick porous water ice subjected to the bombardment of ions and electrons to better understand the formation of exospheres of the icy moons of Jupiter. Our ice samples are as similar as possible to the expected moon surfaces but surface charging of the samples during ion irradiation may distort the experimental results. We therefore monitor the time scales for charging and dis- charging of the samples when subjected to a beam of ions. These experiments allow us to derive an electric conductivity of deep porous ice layers. The results imply that electron irradiation and sputtering play a non-negligible role for certain plasma conditions at the icy moons of Jupiter. The observed ion sputtering yields from our ice samples are similar to previous experiments where compact ice films were sputtered off a micro-balance.

DXL (Diffuse X-rays from the Local Galaxy) is a sounding rocket mission designed to quantify and characterize the contribution of Solar Wind Charge eXchange (SWCX) to the Diffuse X-ray Background and study the properties of the Local Hot Bubble (LHB). Based on the results from the DXL mission, we quantified and removed the contribution of SWCX to the diffuse X-ray background measured by the ROSAT All Sky Survey (RASS). The "cleaned" maps were used to investigate the physical properties of the LHB. Assuming thermal ionization equilibrium, we measured a highly uniform temperature distributed around kT=0.097 keV+/-0.013 keV (FWHM)+/-0.006 keV (systematic). We also generated a thermal emission measure map and used it to characterize the three-dimensional (3D) structure of the LHB which we found to be in good agreement with the structure of the local cavity measured from dust and gas.

During its two years mission around comet 67P/Churyumov-Gerasimenko, ESA's Rosetta spacecraft had the unique opportunity to follow closely a comet in the most active part of its orbit. Many studies have presented the typical features associated to the activity of the nucleus, such as localized dust and gas jets. Here we report on series of more energetic transient events observed during the three months surrounding the comet's perihelion passage in August 2015. We detected and characterized 34 outbursts with the Rosetta cameras, one every 2.4 nucleus rotation. We identified 3 main dust plume morphologies associated to these events: a narrow jet, a broad fan, and more complex plumes featuring both previous types together. These plumes are comparable in scale and temporal variation to what has been observed on other comets. We present a map of the outbursts source locations, and discuss the associated topography. We find that the spatial distribution sources on the nucleus correlates well with morphological region boundaries, especially in areas marked by steep scarps or cliffs. Outbursts occur either in the early morning or shortly after the local noon, indicating two potential processes: Morning outbursts may be triggered by thermal stresses linked to the rapid change of temperature, afternoon events are most likely related to the diurnal or seasonal heat wave reaching volatiles buried under the first surface layer. In addition, we propose that some events can be the result of a completely different mechanism, in which most of the dust is released upon the collapse of a cliff.

Beginning in March 2014, the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras began capturing images of the nucleus and coma (gas and dust) of comet 67P/Churyumov-Gerasimenko using both the wide angle camera (WAC) and the narrow angle camera (NAC). The many observations taken since July of 2014 have been used to study the morphology, location, and temporal variation of the comet's dust jets. We analyzed the dust monitoring observations shortly after the southern vernal equinox on May 30 and 31, 2015 with the WAC at the heliocentric distance Rh = 1.53 AU, where it is possible to observe that the jet rotates with the nucleus. We found that the decline of brightness as a function of the distance of the jet is much steeper than the background coma, which is a first indication of sublimation. We adapted a model of sublimation of icy aggregates and studied the effect as a function of the physical properties of the aggregates (composition and size). The major finding of this article was that through the sublimation of the aggregates of dirty grains (radius a between 5 microm and 50 microm) we were able to completely reproduce the radial brightness profile of a jet beyond 4 km from the nucleus. To reproduce the data we needed to inject a number of aggregates between 8.5 x $10^{13}$ and 8.5 x $10^{10}$ for a = 5 microm and 50 microm respectively, or an initial mass of $H_2O$ ice around 22kg.

We present OSIRIS/NAC observations of decimetre-sized, likely ice-containing aggregates ejected from a confined region on the surface of comet 67P/Churyumov-Gerasimenko. The images were obtained in January 2016 when the comet was at 2 AU from the Sun out-bound from perihelion. We measure the acceleration of individual aggregates through a two-hour image series. Approximately 50% of the aggregates are accelerated away from the nucleus, and 50% towards it, and likewise towards either horizontal direction. The accelerations are up to one order of magnitude stronger than local gravity, and are most simply explained by the combined effect of gas drag accelerating all aggregates upwards, and the recoil force from asymmetric outgassing, either from rotating aggregates with randomly oriented spin axes and sufficient thermal inertia to shift the temperature maximum away from an aggregate's subsolar region, or from aggregates with variable ice content. At least 10% of the aggregates will escape the gravity field of the nucleus and feed the comet's debris trail, while others may fall back to the surface and contribute to the deposits covering parts of the northern hemisphere. The rocket force plays a crucial role in pushing these aggregates back towards the surface. Our observations show the future back fall material in the process of ejection, and provide the first direct measurement of the acceleration of aggregates in the innermost coma (<2km) of a comet, where gas drag is still significant.

We report the detections of a giant planet (MARVELS-7b) and a brown dwarf candidate (MARVELS-7c) around the primary star in the close binary system, HD 87646. It is the first close binary system with more than one substellar circum-primary companion discovered to the best of our knowledge. The detection of this giant planet was accomplished using the first multi-object Doppler instrument (KeckET) at the Sloan Digital Sky Survey (SDSS) telescope. Subsequent radial velocity observations using ET at Kitt Peak National Observatory, HRS at HET, the "Classic" spectrograph at the Automatic Spectroscopic Telescope at Fairborn Observatory, and MARVELS from SDSS-III confirmed this giant planet discovery and revealed the existence of a long-period brown dwarf in this binary. HD 87646 is a close binary with a separation of $\sim22$ AU between the two stars, estimated using the Hipparcos catalogue and our newly acquired AO image from PALAO on the 200-inch Hale Telescope at Palomar. The primary star in the binary, HD 87646A, has Teff = 5770$\pm$80K, log(g)=4.1$\pm$0.1 and [Fe/H] = $-0.17\pm0.08$. The derived minimum masses of the two substellar companions of HD 87646A are 12.4$\pm$0.7M$_{\rm Jup}$ and 57.0$\pm3.7$M$_{\rm Jup}$. The periods are 13.481$\pm$0.001 days and 674$\pm$4 days and the measured eccentricities are 0.05$\pm$0.02 and 0.50$\pm$0.02 respectively. Our dynamical simulations show the system is stable if the binary orbit has a large semi-major axis and a low eccentricity, which can be verified with future astrometry observations.

The Wide Angle Camera of the OSIRIS instrument on board the Rosetta spacecraft is equipped with several narrowband filters that are centered on the emission lines and bands of various fragment species. These are used to determine the evolution of the production and spatial distribution of the gas in the inner coma of comet 67P with time and heliocentric distance, here between 2.6 - 1.3 AU pre-perihelion. Our observations indicate that the emission observed in the OH, OI, CN, NH, and NH2 filters is mostly produced by dissociative electron impact excitation of different parent species. We conclude that CO2 rather than H2O is a significant source of the [OI] 630 nm emission. A strong plume-like feature observed in the in CN and [OI] filters is present throughout our observations. This plume is not present in OH emission and indicates a local enhancement of the CO2/H2O ratio by as much as a factor of 3. We observed a sudden decrease in intensity levels after March 2015, which we attribute to decreased electron temperatures in the first kilometers above the nucleus surface.

DXL (Diffuse X-ray emission from the Local Galaxy) is a sounding rocket mission designed to estimate the contribution of Solar Wind Charge eXchange (SWCX) to the Diffuse X-ray Background (DXB) and to help determine the properties of the Local Hot Bubble (LHB). The detectors are large-area thin-window proportional counters with a spectral response similar to that of the PSPC used in the ROSAT All Sky Survey (RASS). A direct comparison of DXL and RASS data for the same part of the sky viewed from quite different vantage points in the Solar system and the assumption of approximate isotropy for the Solar wind allowed us to quantify the SWCX contribution to all 6 RASS bands (R1-R7, excepting R3). We find that the SWCX contribution at l=140 deg, b=0 deg, where the DXL path crosses the Galactic plane is 33%+-6% (statistical)+-12%(systematic) for R1, 44%+-\%+-5% for R2, 18%+-12%+-11% for R4, 14%+-11%+-9% for R5, and negligible for R6 and R7 bands. Reliable models for the distribution of neutral H and He in the Solar system permit estimation of the contribution of interplanetary SWCX emission over the the whole sky and correction of the RASS maps. We find that the average SWCX contribution in the whole sky is 26%+-6%+-13% for R1, 30%+-4%+-4% for R2, 8%+-5%+-5% for R4, 6%+-4%+-4% for R5, and negligible for R6 and R7.

Polarization maps of the Vela C molecular cloud were obtained at 250, 350, and 500um during the 2012 flight of the balloon-borne telescope BLASTPol. These measurements are used in conjunction with 850um data from Planck to study the submillimeter spectrum of the polarization fraction for this cloud. The spectrum is relatively flat and does not exhibit a pronounced minimum at \lambda ~350um as suggested by previous measurements of other molecular clouds. The shape of the spectrum does not depend strongly on the radiative environment of the dust, as quantified by the column density or the dust temperature obtained from Herschel data. The polarization ratios observed in Vela C are consistent with a model of a porous clumpy molecular cloud being uniformly heated by the interstellar radiation field.

We present results for Vela C obtained during the 2012 flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). We mapped polarized intensity across almost the entire extent of this giant molecular cloud, in bands centered at 250, 350, and 500 {\mu}m. In this initial paper, we show our 500 {\mu}m data smoothed to a resolution of 2.5 arcminutes (approximately 0.5 pc). We show that the mean level of the fractional polarization p and most of its spatial variations can be accounted for using an empirical three-parameter power-law fit, p = p_0 N^(-0.4) S^(-0.6), where N is the hydrogen column density and S is the polarization-angle dispersion on 0.5 pc scales. The decrease of p with increasing S is expected because changes in the magnetic field direction within the cloud volume sampled by each measurement will lead to cancellation of polarization signals. The decrease of p with increasing N might be caused by the same effect, if magnetic field disorder increases for high column density sightlines. Alternatively, the intrinsic polarization efficiency of the dust grain population might be lower for material along higher density sightlines. We find no significant correlation between N and S. Comparison of observed submillimeter polarization maps with synthetic polarization maps derived from numerical simulations provides a promising method for testing star formation theories. Realistic simulations should allow for the possibility of variable intrinsic polarization efficiency. The measured levels of correlation among p, N, and S provide points of comparison between observations and simulations.

The ESA Rosetta spacecraft, currently orbiting around comet 67P, has already provided in situ measurements of the dust grain properties from several instruments, particularly OSIRIS and GIADA. We propose adding value to those measurements by combining them with ground-based observations of the dust tail to monitor the overall, time-dependent dust-production rate and size distribution. To constrain the dust grain properties, we take Rosetta OSIRIS and GIADA results into account, and combine OSIRIS data during the approach phase (from late April to early June 2014) with a large data set of ground-based images that were acquired with the ESO Very Large Telescope (VLT) from February to November 2014. A Monte Carlo dust tail code has been applied to retrieve the dust parameters. Key properties of the grains (density, velocity, and size distribution) were obtained from Rosetta observations: these parameters were used as input of the code to considerably reduce the number of free parameters. In this way, the overall dust mass-loss rate and its dependence on the heliocentric distance could be obtained accurately. The dust parameters derived from the inner coma measurements by OSIRIS and GIADA and from distant imaging using VLT data are consistent, except for the power index of the size-distribution function, which is $\alpha$=--3, instead of $\alpha$=--2, for grains smaller than 1 mm. This is possibly linked to the presence of fluffy aggregates in the coma. The onset of cometary activity occurs at approximately 4.3 au, with a dust production rate of 0.5 kg/s, increasing up to 15 kg/s at 2.9 au. This implies a dust-to-gas mass ratio varying between 3.8 and 6.5 for the best-fit model when combined with water-production rates from the MIRO experiment.

We analyze the physical properties and dynamical origin of a curved jet of comet 67P/Churyumov-Gerasimenko that was observed repeatedly in several nucleus rotations starting on May 30 and persisting until early August, 2015. We simulated the motion of dust grains ejected from the nucleus surface under the influence of the gravity and viscous drag effect of the expanding gas flow from the rotating nucleus. The formation of the curved jet is a combination of the size of the dust particles (~0.1-1 mm) and the location of the source region near the nucleus equator. This enhances the spiral feature of the collimated dust stream after the dust is accelerated to a terminal speed on the order of m/s.

Dust jets, i.e. fuzzy collimated streams of cometary material arising from the nucleus, have been observed in-situ on all comets since the Giotto mission flew by comet 1P/Halley in 1986. Yet their formation mechanism remains unknown. Several solutions have been proposed, from localized physical mechanisms on the surface/sub-surface (see review in Belton (2010)) to purely dynamical processes involving the focusing of gas flows by the local topography (Crifo et al. 2002). While the latter seems to be responsible for the larger features, high resolution imagery has shown that broad streams are composed of many smaller features (a few meters wide) that connect directly to the nucleus surface. We monitored these jets at high resolution and over several months to understand what are the physical processes driving their formation, and how this affects the surface. Using many images of the same areas with different viewing angles, we performed a 3-dimensional reconstruction of collimated jets, and linked them precisely to their sources on the nucleus. Results.We show here observational evidence that the Northern hemisphere jets of comet 67P arise from areas with sharp topographic changes and describe the physical processes involved. We propose a model in which active cliffs are the main source of jet-like features, and therefore the regions eroding the fastest on comets. We suggest that this is a common mechanism taking place on all comets.

We use an existing laboratory facility for space hardware calibration in vacuum to study the impact of energetic ions on water ice. The experiment is intended to simulate the conditions on the surface of Jupiter's icy moons. We present first results of ion sputtering in a sample of porous ice, including the first experimental results for sulphur ion sputtering of ice. The results confirm theoretical predictions and extrapolations from previous sputtering experiments obtained at different impact angles for non-porous water ice.

We investigate the thermal evolution of comet 67P/Churyumov-Gerasimenko's subsurface in the Seth_01 region, where active pits have been observed by the ESA/Rosetta mission. Our simulations show that clathrate destabilization and amorphous ice crystallization can occur at depths corresponding to those of the observed pits in a timescale shorter than 67P/Churyumov-Gerasimenko's lifetime in the comet's activity zone in the inner solar system. Sublimation of crystalline ice down to such depths is possible only in the absence of a dust mantle, which requires the presence of dust grains in the matrix small enough to be dragged out by gas from the pores. Our results are consistent with both pits formation via sinkholes or subsequent to outbursts, the dominant process depending on the status of the subsurface porosity. A sealed dust mantle would favor episodic and disruptive outgassing as a result of an increasing gas pressure in the pores, while a high porosity should allow the formation of large voids in the subsurface due to the continuous escape of volatiles. We finally conclude that the subsurface of 67P/Churyumov-Gerasimenko is not uniform at a spatial scale of 100-200~m.

We report on the first major temporal morphological changes observed on the surface of the nucleus of comet 67P/Churyumov-Gerasimenko, in the smooth terrains of the Imhotep region. We use images of the OSIRIS cameras onboard Rosetta to follow the temporal changes from 24 May 2015 to 11 July 2015. The morphological changes observed on the surface are visible in the form of roundish features, which are growing in size from a given location in a preferential direction, at a rate of 5.6 - 8.1$\times$10$^{-5}$ m s$^{-1}$ during the observational period. The location where changes started and the contours of the expanding features are bluer than the surroundings, suggesting the presence of ices (H$_2$O and/or CO$_2$) exposed on the surface. However, sublimation of ices alone is not sufficient to explain the observed expanding features. No significant variations in the dust activity pattern are observed during the period of changes.

We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface, using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies and mechanical considerations. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 deg) are covered by a fine material and contain a few large ($>$10 m) and isolated boulders, ii) intermediate-slope terrains (20-45 deg) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from $<$1 m to 10 m for the majority of them, and iii) high-slope terrains (45-90 deg) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.

The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13$\pm$0.01 in the HG system formalism and an absolute magnitude $H_v(1,1,0)$ = 15.74$\pm$0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at $\sim$ 290 nm that is possibly due to SO$_2$ ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11\%/(100 nm) to 16\%/(100 nm) in the 1.3$^{\circ}$--54$^{\circ}$ phase angle range. The geometric albedo of the comet is 6.5$\pm$0.2\% at 649 nm, with local variations of up to $\sim$ 16\% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

The production factor, or broad band averaged cross-section, for solar wind charge-exchange with hydrogen producing emission in the ROSAT 1/4 keV (R12) band is $3.8\pm0.2\times10^{-20}$ count degree$^{-2}$ cm$^4$. This value is derived from a comparison of the Long-Term (background) Enhancements in the ROSAT All-Sky Survey with magnetohysdrodynamic simulations of the magnetosheath. This value is 1.8 to 4.5 times higher than values derived from limited atomic data, suggesting that those values may be missing a large number of faint lines. This production factor is important for deriving the exact amount of 1/4 keV band flux that is due to the Local Hot Bubble, for planning future observations in the 1/4 keV band, and for evaluating proposals for remote sensing of the magnetosheath. The same method cannot be applied to the 3/4 keV band as that band, being composed primarily of the oxygen lines, is far more sensitive to the detailed abundances and ionization balance in the solar wind. We also show, incidentally, that recent efforts to correlate XMM-Newton observing geometry with magnetosheath solar wind charge-exchange emission in the oxygen lines have been, quite literally, misguided. Simulations of the inner heliosphere show that broader efforts to correlate heliospheric solar wind charge-exchange with local solar wind parameters are unlikely to produce useful results.

The Solar neighborhood is the closest and most easily studied sample of the Galactic interstellar medium, an understanding of which is essential for models of star formation and galaxy evolution. Observations of an unexpectedly intense diffuse flux of easily-absorbed 1/4 keV X rays, coupled with the discovery that interstellar space within ~100 pc of the Sun is almost completely devoid of cool absorbing gas led to a picture of a "local cavity" filled with X-ray emitting hot gas dubbed the local hot bubble. This model was recently upset by suggestions that the emission could instead be produced readily within the solar system by heavy solar wind ions charge exchanging with neutral H and He in interplanetary space, potentially removing the major piece of evidence for the existence of million-degree gas within the Galactic disk. Here we report results showing that the total solar wind charge exchange contribution is 40% +/- 5% (stat) +/- 5% (sys) of the 1/4 keV flux in the Galactic plane. The fact that the measured flux is not dominated by charge exchange supports the notion of a million-degree hot bubble of order 100 pc extent surrounding the Sun.

An attitude determination system for balloon-borne experiments is presented. The system provides pointing information in azimuth and elevation for instruments flying on stratospheric balloons over Antarctica. In-flight attitude is given by the real-time combination of readings from star cameras, a magnetometer, sun sensors, GPS, gyroscopes, tilt sensors and an elevation encoder. Post-flight attitude reconstruction is determined from star camera solutions, interpolated by the gyroscopes using an extended Kalman Filter. The multi-sensor system was employed by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol), an experiment that measures polarized thermal emission from interstellar dust clouds. A similar system was designed for the upcoming flight of SPIDER, a Cosmic Microwave Background polarization experiment. The pointing requirements for these experiments are discussed, as well as the challenges in designing attitude reconstruction systems for high altitude balloon flights. In the 2010 and 2012 BLASTPol flights from McMurdo Station, Antarctica, the system demonstrated an accuracy of <5' rms in-flight, and <5" rms post-flight.

We present the thermal model of the Balloon-borne Large-Aperture Submillimeter Telescope for Polarimetry (BLASTPol). This instrument was successfully flown in two circumpolar flights from McMurdo, Antarctica in 2010 and 2012. During these two flights, BLASTPol obtained unprecedented information about the magnetic field in molecular clouds through the measurement of the polarized thermal emission of interstellar dust grains. The thermal design of the experiment addresses the stability and control of the payload necessary for this kind of measurement. We describe the thermal modeling of the payload including the sun-shielding strategy. We present the in-flight thermal performance of the instrument and compare the predictions of the model with the temperatures registered during the flight. We describe the difficulties of modeling the thermal behavior of the balloon-borne platform and establish landmarks that can be used in the design of future balloon-borne instruments.

We present the second generation BLASTbus electronics. The primary purposes of this system are detector readout, attitude control, and cryogenic housekeeping, for balloon-borne telescopes. Readout of neutron transmutation doped germanium (NTD-Ge) bolometers requires low noise and parallel acquisition of hundreds of analog signals. Controlling a telescope's attitude requires the capability to interface to a wide variety of sensors and motors, and to use them together in a fast, closed loop. To achieve these different goals, the BLASTbus system employs a flexible motherboard-daughterboard architecture. The programmable motherboard features a digital signal processor (DSP) and field-programmable gate array (FPGA), as well as slots for three daughterboards. The daughterboards provide the interface to the outside world, with versions for analog to digital conversion, and optoisolated digital input/output. With the versatility afforded by this design, the BLASTbus also finds uses in cryogenic, thermometry, and power systems. For accurate timing control to tie everything together, the system operates in a fully synchronous manner. BLASTbus electronics have been successfully deployed to the South Pole, and flown on stratospheric balloons.

Turbulence and magnetic fields are expected to be important for regulating molecular cloud formation and evolution. However, their effects on subparsec to 100 parsec scales, leading to the formation of starless cores, is not well understood. We investigate the prestellar core structure morphologies obtained from analysis of the Herschel-SPIRE 350 $\mu$m maps of the Lupus I cloud. This distribution is first compared on a statistical basis to the large scale shape of the main filament. We find the distribution of the elongation position angle of the cores to be consistent with a random distribution, which means no specific orientation of the morphology of the cores is observed with respect to a large-scale filament shape model for Lupus I, or relative to a large-scale bent filament model. This distribution is also compared to the mean orientation of the large-scale magnetic fields probed at 350 $\mu$m with the Balloon-borne Large Aperture Telescope for Polarimetry (BLASTPol) during its 2010 campaign. Here again we do not find any correlation between the core morphology distribution and the average orientation of the magnetic fields on parsec scales. Our main conclusion is that the local filament dynamics - including secondary filaments that often run orthogonally to the primary filament - and possibly small-scale variations in the local magnetic field direction, could be the dominant factors for explaining the final orientation of each core.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was created by adding polarimetric capability to the BLAST experiment that was flown in 2003, 2005, and 2006. BLASTPol inherited BLAST's 1.8 m primary and its Herschel/SPIRE heritage focal plane that allows simultaneous observation at 250, 350, and 500 {\mu}m. We flew BLASTPol in 2010 and again in 2012. Both were long duration Antarctic flights. Here we present polarimetry of the nearby filamentary dark cloud Lupus I obtained during the 2010 flight. Despite limitations imposed by the effects of a damaged optical component, we were able to clearly detect submillimeter polarization on degree scales. We compare the resulting BLASTPol magnetic field map with a similar map made via optical polarimetry (The optical data were published in 1998 by J. Rizzo and collaborators.). The two maps partially overlap and are reasonably consistent with one another. We compare these magnetic field maps to the orientations of filaments in Lupus I, and we find that the dominant filament in the cloud is approximately perpendicular to the large-scale field, while secondary filaments appear to run parallel to the magnetic fields in their vicinities. This is similar to what is observed in Serpens South via near-IR polarimetry, and consistent with what is seen in MHD simulations by F. Nakamura and Z. Li.

Ground based radial velocity (RV) searches continue to discover exoplanets below Neptune mass down to Earth mass. Furthermore, ground based transit searches now reach milli-mag photometric precision and can discover Neptune size planets around bright stars. These searches will find exoplanets around bright stars anywhere on the sky, their discoveries representing prime science targets for further study due to the proximity and brightness of their host stars. A mission for transit follow-up measurements of these prime targets is currently lacking. The first ESA S-class mission CHEOPS (CHaracterizing ExoPlanet Satellite) will fill this gap. It will perform ultra-high precision photometric monitoring of selected bright target stars almost anywhere on the sky with sufficient precision to detect Earth sized transits. It will be able to detect transits of RV-planets by photometric monitoring if the geometric configuration results in a transit. For Hot Neptunes discovered from the ground, CHEOPS will be able to improve the transit light curve so that the radius can be determined precisely. Because of the host stars' brightness, high precision RV measurements will be possible for all targets. All planets observed in transit by CHEOPS will be validated and their masses will be known. This will provide valuable data for constraining the mass-radius relation of exoplanets, especially in the Neptune-mass regime. During the planned 3.5 year mission, about 500 targets will be observed. There will be 20% of open time available for the community to develop new science programmes.

The European Space Agency's Rosetta spacecraft passed by the main belt asteroid (21) Lutetia the 10th July 2010. With its ~100km size, Lutetia is one of the largest asteroids ever imaged by a spacecraft. During the flyby, the on-board OSIRIS imaging system acquired spectacular images of Lutetia's northern hemisphere revealing a complex surface scarred by numerous impact craters, reaching the maximum dimension of about 55km. In this paper, we assess the cratering history of the asteroid. For this purpose, we apply current models describing the formation and evolution of main belt asteroids, that provide the rate and velocity distributions of impactors. These models, coupled with appropriate crater scaling laws, allow us to interpret the observed crater size-frequency distribution (SFD) and constrain the cratering history. Thanks to this approach, we derive the crater retention age of several regions on Lutetia, namely the time lapsed since their formation or global surface reset. We also investigate the influence of various factors -like Lutetia's bulk structure and crater obliteration- on the observed crater SFDs and the estimated surface ages. From our analysis, it emerges that Lutetia underwent a complex collisional evolution, involving major local resurfacing events till recent times. The difference in crater density between the youngest and oldest recognized units implies a difference in age of more than a factor of 10. The youngest unit (Beatica) has an estimated age of tens to hundreds of Myr, while the oldest one (Achaia) formed during a period when the bombardment of asteroids was more intense than the current one, presumably around 3.6Gyr ago or older.

We identify and investigate the nature of the 20 brightest 250-micron sources detected by the BLAST within the central 150 sq. arcmin of the GOODS-South field. Aided by the available deep VLA radio imaging, reaching S_1.4 = 30 micro-Jy, we have identified radio counterparts for 17/20 of the 250-micron sources. The resulting enhanced positional accuracy of ~1 arcsec has then allowed us to exploit the deep multi-frequency imaging of GOODS-South to establish secure galaxy counterparts for the 17 radio-identified sources, and plausible galaxy candidates for the 3 radio-unidentified sources. Confusion is a serious issue for this deep BLAST 250-micron survey, due to the large size of the beam. Nevertheless, we argue that our chosen counterparts are significant, and often dominant contributors to the measured BLAST flux densities. For all of these 20 galaxies we have been able to determine spectroscopic (8) or photometric (12) redshifts. The result is the first `complete' redshift distribution for a deep 250-micron selected galaxy sample. This reveals that 250-micron surveys reaching detection limits of ~30 mJy contain not only low-redshift spirals/LIRGs, but also the extreme z~2 dust-enshrouded starburst galaxies previously discovered at sub-millimetre wavelengths. Based on their IRAC colours, we find that virtually all of the BLAST galaxy identifications appear better described as analogues of the M82 starburst galaxy, or Sc star-forming discs rather than highly obscured ULIRGs. Inspection of the LABOCA 870-micron imaging of the GOODS-South field yields detections of 7/11 of the z>1 BLAST sources, and reveals 250/870 flux-density ratios consistent with a standard 40K modified black-body fit with a dust emissivity index beta=1.5.

We present results from a survey carried out by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) on a 9 deg^2 field near the South Ecliptic Pole at 250, 350 and 500 {\mu}m. The median 1{\sigma} depths of the maps are 36.0, 26.4 and 18.4 mJy, respectively. We apply a statistical method to estimate submillimeter galaxy number counts and find that they are in agreement with other measurements made with the same instrument and with the more recent results from Herschel/SPIRE. Thanks to the large field observed, the new measurements give additional constraints on the bright end of the counts. We identify 132, 89 and 61 sources with S/N>4 at 250, 350, 500 {\mu}m, respectively and provide a multi-wavelength combined catalog of 232 sources with a significance >4{\sigma} in at least one BLAST band. The new BLAST maps and catalogs are available publicly at http://blastexperiment.info.

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 Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLAST-Pol) is a suborbital mapping experiment designed to study the role played by magnetic fields in the star formation process. BLAST-Pol is the reconstructed BLAST telescope, with the addition of linear polarization capability. Using a 1.8 m Cassegrain telescope, BLAST-Pol images the sky onto a focal plane that consists of 280 bolometric detectors in three arrays, observing simultaneously at 250, 350, and 500 um. The diffraction-limited optical system provides a resolution of 30'' at 250 um. The polarimeter consists of photolithographic polarizing grids mounted in front of each bolometer/detector array. A rotating 4 K achromatic half-wave plate provides additional polarization modulation. With its unprecedented mapping speed and resolution, BLAST-Pol will produce three-color polarization maps for a large number of molecular clouds. The instrument provides a much needed bridge in spatial coverage between larger-scale, coarse resolution surveys and narrow field of view, and high resolution observations of substructure within molecular cloud cores. The first science flight will be from McMurdo Station, Antarctica in December 2010.

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.

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

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 investigated whether the predictions and results of Stanishev et al (2002) concerning a possible relationship between eclipse depths in PX And and its retrograde disc precession phase, could be confirmed in long term observations made by SuperWASP. In addition, two further CVs (DQ Her and V795 Her) in the same SuperWASP data set were investigated to see whether evidence of superhump periods and disc precession periods were present and what other, if any, long term periods could be detected. The results of our period and eclipse analysis on PX And confirm that the negative superhump period is 0.1417 \pm 0.0001d. We find no evidence of positive superhumps in our data suggesting that PX And may have been in a low state during our observations. We improve on existing estimates of the disc precession period and find it to be 4.43 \pm 0.05d. Our results confirm the predictions of Stanishev et al (2002). We find that DQ Her does not appear to show a similar variation for we find no evidence of negative superhumps or of a retrograde disc precession. We also find no evidence of positive superhumps or of a prograde disc precession and we attribute the lack of positive superhumps in DQ Her to be due to the high mass ratio of this CV. We do however find evidence for a modulation of the eclipse depth over a period of 100 days which may be linked with solar-type magnetic cycles which give rise to long term photometric variations. The periodogram analysis for V795 Her detected the likely positive superhump period 0.1165d, however, neither the 0.10826d orbital period or the prograde 1.53d disc precession period were seen. Here though we have found a variation in the periodogram power function at the positive superhump period, over a period of at least 120 days.

We present observations at 250, 350, and 500 um of the nearby galaxy cluster Abell 3112 (z=0.075) carried out with BLAST, the Balloon-borne Large Aperture Submillimeter Telescope. Five cluster members are individually detected as bright submillimetre sources. Their far-infrared SEDs and optical colours identify them as normal star-forming galaxies of high mass, with globally evolved stellar populations. They all have B-R colours of 1.38+/-0.08, transitional between the blue, active population and the red, evolved galaxies that dominate the cluster core. We stack to determine the mean submillimetre emission from all cluster members, which is determined to be 16.6+/-2.5, 6.1+/-1.9, and 1.5+/-1.3 mJy at 250, 350, and 500 um, respectively. Stacking analyses of the submillimetre emission of cluster members reveal trends in the mean far-infrared luminosity with respect to cluster-centric radius and Ks-band magnitude. We find that a large fraction of submillimetre emission comes from the boundary of the inner, virialized region of the cluster, at cluster-centric distances around R_500. Stacking also shows that the bulk of the submillimetre emission arises in intermediate-mass galaxies (L<L*), with Ks magnitude ~1 mag fainter than the giant ellipticals. The results and constraints obtained in this work will provide a useful reference for the forthcoming surveys to be conducted on galaxy clusters by Herschel.

The Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) carried out a 250, 350 and 500 micron survey of the galactic plane encompassing the Vela Molecular Ridge, with the primary goal of identifying the coldest dense cores possibly associated with the earliest stages of star formation. Here we present the results from observations of the Vela-D region, covering about 4 square degrees, in which we find 141 BLAST cores. We exploit existing data taken with the Spitzer MIPS, IRAC and SEST-SIMBA instruments to constrain their (single-temperature) spectral energy distributions, assuming a dust emissivity index beta = 2.0. This combination of data allows us to determine the temperature, luminosity and mass of each BLAST core, and also enables us to separate starless from proto-stellar sources. We also analyze the effects that the uncertainties on the derived physical parameters of the individual sources have on the overall physical properties of starless and proto-stellar cores, and we find that there appear to be a smooth transition from the pre- to the proto-stellar phase. In particular, for proto-stellar cores we find a correlation between the MIPS24 flux, associated with the central protostar, and the temperature of the dust envelope. We also find that the core mass function of the Vela-D cores has a slope consistent with other similar (sub)millimeter surveys.

We describe the application of a statistical method to estimate submillimeter galaxy number counts from confusion limited observations by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). Our method is based on a maximum likelihood fit to the pixel histogram, sometimes called 'P(D)', an approach which has been used before to probe faint counts, the difference being that here we advocate its use even for sources with relatively high signal-to-noise ratios. This method has an advantage over standard techniques of source extraction in providing an unbiased estimate of the counts from the bright end down to flux densities well below the confusion limit. We specifically analyse BLAST observations of a roughly 10 sq. deg. map centered on the Great Observatories Origins Deep Survey South (GOODS-S) field. We provide estimates of number counts at the three BLAST wavelengths, 250, 350, and 500 microns; instead of counting sources in flux bins we estimate the counts at several flux density nodes connected with power-laws. We observe a generally very steep slope for the counts of about -3.7 at 250 microns and -4.5 at 350 and 500 microns, over the range ~0.02-0.5 Jy, breaking to a shallower slope below about 0.015 Jy at all three wavelengths. We also describe how to estimate the uncertainties and correlations in this method so that the results can be used for model-fitting. This method should be well-suited for analysis of data from the Herschel satellite.

We detect correlations in the cosmic far-infrared background due to the clustering of star-forming galaxies in observations made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST, at 250, 350, and 500 microns. We perform jackknife and other tests to confirm the reality of the signal. The measured correlations are well fit by a power law over scales of 5-25 arcminutes, with Delta I/I = 15.1 +/- 1.7%. We adopt a specific model for submillimeter sources in which the contribution to clustering comes from sources in the redshift ranges 1.3 <= z <= 2.2, 1.5 <= z <= 2.7, and 1.7 <= z <= 3.2, at 250, 350, and 500 microns, respectively. With these distributions, our measurement of the power spectrum, P(k_theta), corresponds to linear bias parameters, b = 3.8 +/- 0.6, 3.9 +/- 0.6 and 4.4 +/- 0.7, respectively. We further interpret the results in terms of the halo model, and find that at the smaller scales, the simplest halo model fails to fit our results. One way to improve the fit is to increase the radius at which dark matter halos are artificially truncated in the model, which is equivalent to having some star-forming galaxies at z >= 1 located in the outskirts of groups and clusters. In the context of this model we find a minimum halo mass required to host a galaxy is log (M_min / M_sun) = 11.5 (+0.4/-0.1), and we derive effective biases $b_eff = 2.2 +/- 0.2, 2.4 +/- 0.2, and 2.6 +/- 0.2, and effective masses log (M_eff / M_sun) = 12.9 +/- 0.3, 12.8 +/- 0.2, and 12.7 +/- 0.2, at 250, 350, and 500 microns, corresponding to spatial correlation lengths of r_0 = 4.9, 5.0, and 5.2 +/- 0.7 h^-1 Mpc, respectively. Finally, we discuss implications for clustering measurement strategies with Herschel and Planck.

Over the course of two flights, the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) made resolved maps of seven nearby (<25 Mpc) galaxies at 250, 350, and 500 microns. During its June 2005 flight from Sweden, BLAST observed a single nearby galaxy, NGC 4565. During the December 2006 flight from Antarctica, BLAST observed the nearby galaxies NGC 1097, NGC 1291, NGC 1365, NGC 1512, NGC 1566, and NGC 1808. We fit physical dust models to a combination of BLAST observations and other available data for the galaxies observed by Spitzer. We fit a modified blackbody to the remaining galaxies to obtain total dust mass and mean dust temperature. For the four galaxies with Spitzer data, we also produce maps and radial profiles of dust column density and temperature. We measure the fraction of BLAST detected flux originating from the central cores of these galaxies and use this to calculate a "core fraction," an upper limit on the "AGN fraction" of these galaxies. We also find our resolved observations of these galaxies give a dust mass estimate 5-19 times larger than an unresolved observations would predict. Finally, we are able to use these data to derive a value for the dust mass absorption co-efficient of kappa = 0.29 +/-0.03 m^2 kg^-1 at 250 microns. This study is an introduction to future higher-resolution and higher-sensitivity studies to be conducted by Herschel and SCUBA-II.

We directly measure redshift evolution in the mean physical properties (far-infrared luminosity, temperature, and mass) of the galaxies that produce the cosmic infrared background (CIB), using measurements from the Balloon-borne Large Aperture Sub-millimeter Telescope (BLAST), and Spitzer which constrain the CIB emission peak. This sample is known to produce a surface brightness in the BLAST bands consistent with the full CIB, and photometric redshifts are identified for all of the objects. We find that most of the 70 micron background is generated at z <~ 1 and the 500 micron background generated at z >~ 1. A significant growth is observed in the mean luminosity from ~ 10^9 - 10^12 L_sun, and in the mean temperature by 10 K, from redshifts 0< z < 3. However, there is only weak positive evolution in the comoving dust mass in these galaxies across the same redshift range. We also measure the evolution of the far-infrared luminosity density, and the star-formation rate history for these objects, finding good agreement with other infrared studies up to z ~1, exceeding the contribution attributed to optically-selected galaxies.

We present first results from an unbiased 50 deg^2 submillimeter Galactic survey at 250, 350, and 500 micron from the 2006 flight of the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). The map has resolution ranging from 36 arcsec to 60 arcsec in the three submillimeter bands spanning the thermal emission peak of cold starless cores. We determine the temperature, luminosity, and mass of more than one thousand compact sources in a range of evolutionary stages and an unbiased statistical characterization of the population. From comparison with C^(18)O data, we find the dust opacity per gas mass, kappa r = 0.16 cm^2 g^(-1) at 250 micron, for cold clumps. We find that 2% of the mass of the molecular gas over this diverse region is in cores colder than 14 K, and that the mass function for these cold cores is consistent with a power law with index alpha = -3.22 +/- 0.14 over the mass range 14 M_sun < M < 80 M_sun. Additionally, we infer a mass-dependent cold core lifetime of t_c(M) = 4E6 (M/20 M_sun)^(-0.9) years - longer than what has been found in previous surveys of either low or high mass cores, and significantly longer than free fall or likely turbulent decay times. This implies some form of non-thermal support for cold cores during this early stage of star formation.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has made one square degree, deep, confusion limited maps at three different bands, centered on the Great Observatories Origins Deep Survey South field. By calculating the covariance of these maps with catalogs of 24 micron sources from the Far-Infrared Deep Extragalactic Legacy Survey (FIDEL), we have determined that the total submillimeter intensities are 8.60 +/- 0.59, 4.93 +/- 0.34, and 2.27 +/- 0.20 nW m^2 sr^(-1) at 250, 350, and 500 micron, respectively. These numbers are more precise than previous estimates of the cosmic infrared background (CIB) and are consistent with 24 micron-selected galaxies generating the full intensity of the CIB. We find that the fraction of the CIB that originates from sources at z >= 1.2 increases with wavelength, with 60% from high redshift sources at 500 micron. At all BLAST wavelengths, the relative intensity of high-z sources is higher for 24 micron-faint sources than it is for 24 micron-bright sources. Galaxies identified as active galactic nuclei (AGN) by their Infrared Array Camera (IRAC) colors are 1.6-2.6 times brighter than the average population at 250-500 micron, consistent with what is found for X-ray-selected AGN. BzK-selected galaxies are found to be moderately brighter than typical 24 micron-selected galaxies in the BLAST bands. These data provide high precision constraints for models of the evolution of the number density and intensity of star forming galaxies at high redshift.

We present the 250, 350, and 500 micron detection of bright submillimeter emission in the direction of the Bullet Cluster measured by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). The 500 micron centroid is coincident with an AzTEC 1.1 mm point-source detection at a position close to the peak lensing magnification produced by the cluster. However, the 250 micron and 350 micron centroids are elongated and shifted toward the south with a differential shift between bands that cannot be explained by pointing uncertainties. We therefore conclude that the BLAST detection is likely contaminated by emission from foreground galaxies associated with the Bullet Cluster. The submillimeter redshift estimate based on 250-1100 micron photometry at the position of the AzTEC source is z_phot = 2.9 (+0.6 -0.3), consistent with the infrared color redshift estimation of the most likely IRAC counterpart. These flux densities indicate an apparent far-infrared luminosity of L_FIR = 2E13 Lsun. When the amplification due to the gravitational lensing of the cluster is removed, the intrinsic far-infrared luminosity of the source is found to be L_FIR <= 10^12 Lsun, consistent with typical luminous infrared galaxies.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully during a 250-hour flight over Antarctica in December 2006 (BLAST06). As part of the calibration and pointing procedures, the red hypergiant star VY CMa was observed and used as the primary calibrator. Details of the overall BLAST06 calibration procedure are discussed. The 1-sigma absolute calibration is accurate to 10, 12, and 13% at the 250, 350, and 500 micron bands, respectively. The errors are highly correlated between bands resulting in much lower error for the derived shape of the 250-500 micron continuum. The overall pointing error is <5" rms for the 36, 42, and 60" beams. The performance of the optics and pointing systems is discussed.

We have identified radio and/or mid-infrared counterparts to 198 out of 350 sources detected at >=5 sigma over ~ 9 square degrees centered on the Chandra Deep Field South (CDFS) by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) at 250, 350 and 500 um. We have matched 114 of these counterparts to optical sources with previously derived photometric redshifts and fitted SEDs to the BLAST fluxes and fluxes at 70 and 160 um acquired with the Spitzer Space Telescope. In this way, we have constrained dust temperatures, total far-infrared/sub-millimeter luminosities and star formation rates for each source. Our findings show that on average, the BLAST sources lie at significantly lower redshifts and have significantly lower rest-frame dust temperatures compared to submm sources detected in surveys conducted at 850 um. We demonstrate that an apparent increase in dust temperature with redshift in our sample arises as a result of selection effects. Finally, we provide the full multi-wavelength catalog of >= 5 sigma BLAST sources contained within the complete ~ 9 square degree survey area.

Submillimetre surveys during the past decade have discovered a population of luminous, high-redshift, dusty starburst galaxies. In the redshift range 1 <= z <= 4, these massive submillimetre galaxies go through a phase characterized by optically obscured star formation at rates several hundred times that in the local Universe. Half of the starlight from this highly energetic process is absorbed and thermally re-radiated by clouds of dust at temperatures near 30 K with spectral energy distributions peaking at 100 microns in the rest frame. At 1 <= z <= 4, the peak is redshifted to wavelengths between 200 and 500 microns. The cumulative effect of these galaxies is to yield extragalactic optical and far-infrared backgrounds with approximately equal energy densities. Since the initial detection of the far-infrared background (FIRB), higher-resolution experiments have sought to decompose this integrated radiation into the contributions from individual galaxies. Here we report the results of an extragalactic survey at 250, 350 and 500 microns. Combining our results at 500 microns with those at 24 microns, we determine that all of the FIRB comes from individual galaxies, with galaxies at z >= 1.2 accounting for 70 per cent of it. As expected, at the longest wavelengths the signal is dominated by ultraluminous galaxies at z > 1.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a sub-orbital experiment designed to study the process of star formation in local galaxies (including the Milky Way) and in galaxies at cosmological distances. Using a 2-m Cassegrain telescope, BLAST images the sky onto a focal plane, which consists of 270 bolometric detectors split between three arrays, observing simultaneously in 30% wide bands, centered at 250, 350, and 500 microns. The diffraction-limited optical system provides a resolution of 30" at 250 microns. The pointing system enables raster-like scans with a positional accuracy of ~30", reconstructed to better than 5" rms in post-flight analysis. BLAST had two successful flights, from the Arctic in 2005, and from Antarctica in 2006, which provided the first high-resolution and large-area (~0.8-200 deg^2) submillimeter surveys at these wavelengths. As a pathfinder for the SPIRE instrument on Herschel, BLAST shares with the ESA satellite similar focal plane technology and scientific motivation. A third flight in 2009 will see the instrument modified to be polarization-sensitive (BLAST-Pol). With its unprecedented mapping speed and resolution, BLAST-Pol will provide insights into Galactic star-forming nurseries, and give the necessary link between the larger, coarse resolution surveys and the narrow, resolved observations of star-forming structures from space and ground based instruments being commissioned in the next 5 years.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a sub-orbital surveying experiment designed to study the evolutionary history and processes of star formation in local galaxies (including the Milky Way) and galaxies at cosmological distances. The BLAST continuum camera, which consists of 270 detectors distributed between 3 arrays, observes simultaneously in broad-band (30%) spectral-windows at 250, 350, and 500 microns. The optical design is based on a 2m diameter telescope, providing a diffraction-limited resolution of 30" at 250 microns. The gondola pointing system enables raster mapping of arbitrary geometry, with a repeatable positional accuracy of ~30"; post-flight pointing reconstruction to ~5" rms is achieved. The on-board telescope control software permits autonomous execution of a pre-selected set of maps, with the option of manual override. In this paper we describe the primary characteristics and measured in-flight performance of BLAST. BLAST performed a test-flight in 2003 and has since made two scientifically productive long-duration balloon flights: a 100-hour flight from ESRANGE (Kiruna), Sweden to Victoria Island, northern Canada in June 2005; and a 250-hour, circumpolar-flight from McMurdo Station, Antarctica, in December 2006.