In 2017, 2018 and 2019, comets 46P/Wirtanen, 45P/Honda-Mrkos-Pajdusakova, and 41P/Tuttle-Giacobini-Kresak all had perihelion passages. Their hydrogen comae were observed by the Solar Wind ANisotropies (SWAN) all-sky hydrogen Lyman-alpha Camera on the SOlar and Heliospheric Observer (SOHO) satellite: comet 46P for the fourth time and comets 45P and 41P for the third time each since 1997. Comet 46P/Wirtanen is one of a small class of so-called hyperactive comets whose gas production rates belie their small size. This comet was the original target comet of the Rosetta mission. The Solar Wind ANisotropies (SWAN) all-sky hydrogen Lyman-alpha camera on the SOlar and Heliospheric Observer (SOHO) satellite observed the hydrogen coma of comet 46P/Wirtanen during the apparitions of 1997, 2002, 2008 and 2018. Over the 22 years, the activity decreased and its variation with heliocentric distance has changed markedly in a way very similar to that of another hyperactive comet, 103P/Hartley 2. Comet 45P/Honda-Mrkos-Pajdusakova was observed by SWAN during its perihelion apparitions of 2001, 2011 and 2017. Over this time period the activity level has remained remarkably similar, with no long-term fading or abrupt decreases. Comet 41P/Tuttle-Giacobini-Kresak was observed by SWAN in its perihelion apparitions of 2001, 2006 and 2017 and has decreased in activity markedly over the same time period. In 1973 it was known for large outbursts, which continued during the 2001 (2 outbursts) and 2006 (1 outburst) apparitions. However, over the 2001 to 2017 time period covered by the SOHO/SWAN observations the water production rates have greatly decreased by factors of 10-30 over corresponding times during its orbit.

The Alice far-UV imaging spectrograph (700-2050 A) acquired over 70,000 spectral images during Rosetta's 2-year escort mission, including over 20,000 in the months surrounding perihelion when the comet activity level was highest. We have developed automated software to fit and remove ubiquitous H, O, C, S, and CO emissions from Alice spectra, along with reflected solar continuum and absorption from gaseous H2O in the comet's coma, which we apply to a "grand sum" of integrations taken near perihelion. We present upper limits on the presence of one ion and 17 neutral atomic species for this time period. These limits are compared to results obtained by other Rosetta instruments where possible, as well as to CI carbonaceous chondrites and solar photospheric abundances.

The Solar Wind ANisotropies (SWAN) all-sky hydrogen Lyman-alpha camera on the SOlar and Heliospheric Observer (SOHO) satellite observed the hydrogen coma of comet C/2017 S3 (PanSTARRS) for the last month of its activity from 2018 July 4 to August 4 and what appears to have been its final disintegration just 11 days before its perihelion on August 15. The hydrogen coma indicated water production had a small outburst on July 8 at a heliocentric distance of 1.1AU and then a much larger one on July 20 at 0.8 AU. Over the following two weeks the water production dropped by more than a factor of ten after which it was no longer detectable. The behavior is reminiscent of comet C/1999 S4 (LINEAR) in 2000, which had a few small outbursts on its inbound orbit and a major outburst at a heliocentric distance of about 0.8 AU, which was close to its perihelion, followed by its complete disintegration that was documented by several sets of observations including SWAN. C/2017 S3 (PanSTARRS) however had a much larger water production rate than C/ 1999 S4 (LINEAR). Here we estimate the size of the nucleus of C/2017 S3 just before its final outburst and apparent disintegration was estimated using the total amount of water produced during its last weeks for a range of values of the refractory/ice ratio in the nucleus. We also determine the size distribution of the disintegrating particles as the comet faded.

On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the one calculated using a thermo-physical model. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O) and 4200 nm (CO2). We measured the average signal in a circular annulus, to study spatial variation around the comet, and in a sector of the annulus, to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet centre. The spatial correlation between dust and water, both coming from the sun-lit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0deg sub-solar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production, with respect to the value calculated using a simple thermo-physical model, is observed when the head lobe and regions of the Southern hemisphere with strong seasonal variations are illuminated. A drastic decrease in dust production, when the water production (both measured and from the model) displays a maximum, happens when typical Northern consolidated regions are illuminated and the Southern hemisphere regions with strong seasonal variations are instead in shadow.

The southern hemisphere of comet 67P/Churyumov-Gerasimenko (67P) became observable by the Rosetta mission in March 2015, a few months before cometary southern vernal equinox. The Anhur region in the southern part of the comet's larger lobe was found to be highly eroded, enriched in volatiles, and highly active. We analyze high-resolution images of the Anhur region pre- and post-perihelion acquired by the OSIRIS imaging system on board the Rosetta mission. The Narrow Angle Camera is particularly useful for studying the evolution in Anhur in terms of morphological changes and color variations.}{Radiance factor images processed by the OSIRIS pipeline were coregistered, reprojected onto the 3D shape model of the comet, and corrected for the illumination conditions. We find a number of morphological changes in the Anhur region that are related to formation of new scarps; removal of dust coatings; localized resurfacing in some areas, including boulders displacements; and vanishing structures, which implies localized mass loss that we estimate to be higher than 50 million kg. The strongest changes took place in and nearby the Anhur canyon-like structure, where significant dust cover was removed, an entire structure vanished, and many boulders were rearranged. All such changes are potentially associated with one of the most intense outbursts registered by Rosetta during its observations, which occurred one day before perihelion passage. Moreover, in the niche at the foot of a new observed scarp, we also see evidence of water ice exposure that persisted for at least six months. The abundance of water ice, evaluated from a linear mixing model, is relatively high (> 20%). Our results confirm that the Anhur region is volatile-rich and probably is the area on 67P with the most pristine exposures near perihelion.

Following our previous detection of ubiquitous H2O and O2 absorption against the far-UV continuum of stars located near the nucleus of Comet 67P/Churyumov-Gerasimenko, we present a serendipitously observed stellar occultation that occurred on 2015 September 13, approximately one month after the comet's perihelion passage. The occultation appears in two consecutive 10-minute spectral images obtained by Alice, Rosetta's ultraviolet (700-2100 A) spectrograph, both of which show H2O absorption with column density $>10^{17.5} \mathrm{cm}^{-2}$ and significant O2 absorption ($\mathrm{O2/H2O} \approx 5$-10%). Because the projected distance from the star to the nucleus changes between exposures, our ability to study the H2O column density profile near the nucleus (impact parameters $<1$ km) is unmatched by our previous observations. We find that the H2O and O2 column densities decrease with increasing impact parameter, in accordance with expectations, but the O2 column decreases $\sim3$ times more quickly than H2O. When combined with previously published results from stellar appulses, we conclude that the O2 and H2O column densities are highly correlated, and O2/H2O decreases with increasing H2O column.

In April 2016, the Rosetta spacecraft performed a low-altitude low-phase-angle flyby over the Imhotep-Khepry transition of 67P/Churyumov-Gerasimenko's nucleus. The OSIRIS/Narrow-Angle-Camera (NAC) acquired 112 images with mainly 3 broadband filters in the visible at a resolution of up to 0.53 m/px and for phase angles between 0.095{\deg} and 62{\deg}. Using those images, we have investigated the morphological and spectrophotometrical properties of this area. We assembled the images into coregistered color cubes. Using a 3D shape model, we produced the illumination conditions and georeference for each image. We projected the observations on a map to investigate its geomorphology. Observations were photometrically corrected using the Lommel-Seeliger disk law. Spectrophotometric analyses were performed on the coregistered color cubes. These data were used to estimate the local phase reddening. This region of the nucleus hosts numerous and varied types of terrains and features. We observe an association between a feature's nature, its reflectance, and its spectral slope. Fine material deposits exhibit an average reflectance and spectral slope, while terrains with diamictons, consolidated material, degraded outcrops, or features such as somber boulders, present a lower-than-average reflectance and higher-than-average spectral slope. Bright surfaces present here a spectral behavior consistent with terrains enriched in water-ice. We find a phase-reddening slope of 0.064{\pm}0.001{\%}/100nm/{\deg} at 2.7 au outbound, similarly to the one obtained at 2.3 au inbound during the February 2015 flyby. Identified as the source region of multiple jets and a host of water-ice material, the Imhotep-Khepry transition appeared in April 2016, close to the frost line, to further harbor several potential locations with exposed water-ice material among its numerous different morphological terrain units.

The phase function of the dust coma of comet 67P has been determined from Rosetta/OSIRIS images \citep{Bertini17}. This function show a deep minimum at phase angles near 100$^\circ$, and a strong backscattering enhancement. These two properties cannot be reproduced by regular models of cometary dust, most of them based on wavelength-sized and randomly-oriented aggregate particles. We show, however, that an ensamble of oriented elongated particles of a wide variety of aspect ratios, with radii $r \gtrsim$10 $\mu$m, and whose long axes are perpendicular to the direction of the solar radiation, are capable of reproducing the observed phase function. These particles must be absorbing, with an imaginary part of the refractive index of about 0.1 to match the expected geometric albedo, and with porosity in the 60-70\% range.

The Rosetta space probe accompanied comet 67P/Churyumov-Gerasimenko for more than two years, obtaining an unprecedented amount of unique data of the comet nucleus and inner coma. This work focuses identifying the source regions of faint jets and outbursts and on studying the spectrophotometric properties of some outbursts. We use observations acquired with the OSIRIS/NAC camera during July-October 2015, that is, close to perihelion. More than 200 jets of different intensities were identified directly on the nucleus. Some of the more intense outbursts appear spectrally bluer than the comet dark terrain in the vivible-to-near-infrared region. We attribute this spectral behavior to icy grains mixed with the ejected dust. Some of the jets have an extremely short lifetime. They appear on the cometary surface during the color sequence observations, and vanish in less than some few minutes after reaching their peak. We also report a resolved dust plume observed in May 2016 at a resolution of 55 cm/pixel, which allowed us to estimate an optical depth of $\sim$0.65 and an ejected mass of $\sim$ 2200 kg. We present the results on the location, duration, and colors of active sources on the nucleus of 67P from the medium-resolution (i.e., 6-10 m/pixel) images acquired close to perihelion passage. The observed jets are mainly located close to boundaries between different morphological regions. Jets depart not only from cliffs, but also from smooth and dust-covered areas, from fractures, pits, or cavities that cast shadows and favor the recondensation of volatiles. This study shows that faint jets or outbursts continuously contribute to the cometary activity close to perihelion passage, and that these events are triggered by illumination conditions. Faint jets or outbursts are not associated with a particular terrain type or morphology.

The Solar Wind Anisotropies (SWAN) instrument on the SOlar and Heliospheric Observatory (SOHO) satellite has observed 44 long period and new Oort cloud comets and 36 apparitions of 17 short period comets since its launch in December 1995. Water production rates have been determined from the over 3700 images producing a consistent set of activity variations over large parts of each comet's orbit. This has enabled the calculation of exponential power-law variations with heliocentric distance of these comets both before and after perihelion, as well as the absolute values of the water production rates. These various measures of overall water activity including pre- and post-perihelion exponents, absolute water production rates at 1 AU, active surface areas and their variations have been compared with a number of dynamical quantities for each comet including dynamical class, original semi-major axis, nucleus radius (when available), and compositional taxonomic class. Evidence for evolution of cometary nuclei is seen in both long-period and short-period comets.

The Alice ultraviolet spectrograph on the European Space Agency Rosetta spacecraft observed comet 67P/Churyumov-Gerasimenko in its orbit around the Sun for just over two years. Alice observations taken in 2015 October, two months after perihelion, show large increases in the comet's Ly-$\beta$, O I 1304, O I 1356, and C I 1657 $\AA$ atomic emission that initially appeared to indicate gaseous outbursts. However, the Rosetta Plasma Consortium instruments showed a coronal mass ejection (CME) impact at the comet coincident with the emission increases, suggesting that the CME impact may have been the cause of the increased emission. The presence of the semi-forbidden O I 1356 $\AA$ emission multiplet is indicative of a substantial increase in dissociative electron impact emission from the coma, suggesting a change in the electron population during the CME impact. The increase in dissociative electron impact could be a result of the interaction between the CME and the coma of 67P or an outburst coincident with the arrival of the CME. The observed dissociative electron impact emission during this period is used to characterize the O2 content of the coma at two peaks during the CME arrival. The mechanism that could cause the relationship between the CME and UV emission brightness is not well constrained, but we present several hypotheses to explain the correlation.

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.

The Alice far-ultraviolet imaging spectrograph onboard Rosetta observed emissions from atomic and molecular species from within the coma of comet 67P/Churyumov-Gerasimenko during the entire escort phase of the mission from 2014 August to 2016 September. The initial observations showed that emissions of atomic hydrogen and oxygen close to the surface were produced by energetic electron impact dissociation of H2O. Following delivery of the lander, Philae, on 2014 November 12, the trajectory of Rosetta shifted to near-terminator orbits that allowed for these emissions to be observed against the shadowed nucleus that, together with the compositional heterogeneity, enabled us to identify unique spectral signatures of dissociative electron impact excitation of H2O, CO2, and O2. CO emissions were found to be due to both electron and photoexcitation processes. Thus we are able, from far-ultraviolet spectroscopy, to qualitatively study the evolution of the primary molecular constituents of the gaseous coma from start to finish of the escort phase. Our results show asymmetric outgassing of H2O and CO2 about perihelion, H2O dominant before and CO2 dominant after, consistent with the results from both the in situ and other remote sensing instruments on Rosetta.

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.

Nine recently discovered long-period comets were observed by the Solar Wind Anisotropies (SWAN) Lyman-alpha all-sky camera on board the Solar and Heliosphere Observatory (SOHO) satellite during the period of 2013 to 2016. These were C/2012 K1 (PanSTARRS), C/2013 US10 (Catalina), C/2013 V5 (Oukaimeden), C/2013 R1 (Lovejoy), C/2014 E2 (Jacques), C/2014 Q2 (Lovejoy), C/2015 G2 (MASTER), C/2014 Q1 (PanSTARRS) and C/2013 X1 (PanSTARRS). Of these 9 comets 6 were long-period comets and 3 were possibly dynamically new. Water production rates were calculated from each of the 885 images using our standard time-resolved model that accounts for the whole water photodissociation chain, exothermic velocities and collisional escape of H atoms. For most of these comets there were enough observations over a broad enough range of heliocentric distances to calculate power-law fits to the variation of production rate with heliocentric distances for pre- and post-perihelion portions of the orbits. Comet C/2014 Q1 (PanSTARRS), with a perihelion distance of only ~0.3 AU, showed the most unusual variation of water production rate with heliocentric distance and the resulting active area variation, indicating that when the comet was within 0.7 AU its activity was dominated by the continuous release of icy grains and chunks, greatly increasing the active sublimation area by more than a factor of 10 beyond what it had at larger heliocentric distances. A possible interpretation suggests that a large fraction of the comet's mass was lost during the apparition.

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.

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.

We have detected H$_2$O and O$_2$ absorption against the far-UV continuum of stars located on lines of sight near the nucleus of Comet 67P/Churyumov-Gerasimenko using the Alice imaging spectrograph on Rosetta. These stellar appulses occurred at impact parameters of $\rho=4$-20 km, and heliocentric distances ranging from $R_h=-1.8$ to 2.3 AU (negative values indicate pre-perihelion observations). The measured H$_2$O column densities agree well with nearly contemporaneous values measured by VIRTIS-H. The clear detection of O$_2$ independently confirms the initial detection by the ROSINA mass spectrometer; however, the relative abundance of O$_2$/H$_2$O derived from the stellar spectra (11%-68%, with a median value of 25%) is considerably larger than published values found by ROSINA. The cause of this difference is unclear, but potentially related to ROSINA measuring number density at the spacecraft position while Alice measures column density along a line of sight that passes near the nucleus.

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.

Aims. A recent debate on the decade-long stability of the interstellar He flow vector, and in particular the flow longitude, has prompted us to check for any variability in the interstellar H flow vector as observed by the SWAN instrument on board SOHO. Methods. We used a simple model-independent method to determine the interstellar H flow longitude, based on the parallax effects induced on the Lyman-{\alpha} intensity measured by SWAN following the satellite motion around the Sun. Results. Our results show that the interstellar H flow vector longitude does not vary significantly from an average value of 252.9$^{\circ}$ $\pm$ 1.4$^{\circ}$ throughout the 20-year span of the SWAN dataset, further strengthening the arguments for the stability of the interstellar gas flow.

We present the results of the photometric and spectrophotometric properties of the 67P/Churyumov-Gerasimenko nucleus derived with the OSIRIS instrument during the closest fly-by over the comet, which took place on 14 th February 2015 at a distance of {\~} 6 km from the surface. Several images covering the 0{\deg}-33{\deg} phase angle range were acquired, and the spatial resolution achieved was 11 cm/pxl. The flown-by region is located on the big lobe of the comet, near the borders of the Ash, Apis and Imhotep regions. Our analysis shows that this region features local heterogeneities at the decimetre scale. We observed difference of reflectance up to 40{\%} between bright spots and sombre regions, and spectral slope variations up to 50{\%}. The spectral reddening effect observed globally on the comet surface by Fornasier et al. (2015) is also observed locally on this region, but with a less steep behaviour. We note that numerous metre-sized boulders, which exhibit a smaller opposition effect, also appear spectrally redder than their surroundings. In this region, we found no evidence linking observed bright spots to exposed water-ice-rich material. We fitted our dataset using the Hapke 2008 photometric model. The region overflown is globally as dark as the whole nucleus (geometric albedo of 6.8{\%}) and it has a high porosity value in the uppermost-layers (86{\%}). These results of the photometric analysis at a decimetre scale indicate that the photometric properties of the flown-by region are similar to those previously found for the whole nucleus.

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.

Since the orbital insertion of the Rosetta spacecraft, comet 67P/Churyumov-Gerasimenko (67P/C-G) has been mapped by OSIRIS camera and VIRTIS spectro-imager, producing a huge quantity of images and spectra of the comet's nucleus. The aim of this work is to search for the presence of H$_2$O on the nucleus which, in general, appears very dark and rich in dehydrated organic material. After selecting images of the bright spots which could be good candidates to search for H$_2$O ice, taken at high resolution by OSIRIS, we check for spectral cubes of the selected coordinates to identify these spots observed by VIRTIS. The selected OSIRIS images were processed with the OSIRIS standard pipeline and corrected for the illumination conditions for each pixel using the Lommel-Seeliger disk law. The spots with higher I/F were selected and then analysed spectrophotometrically and compared with the surrounding area. We selected 13 spots as good targets to be analysed by VIRTIS to search for the 2 micron absorption band of water ice in the VIRTIS spectral cubes. Out of the 13 selected bright spots, eight of them present positive H$_2$O ice detection on the VIRTIS data. A spectral analysis was performed and the approximate temperature of each spot was computed. The H$_2$O ice content was confirmed by modeling the spectra with mixing (areal and intimate) of H$_2$O ice and dark terrain, using Hapke's radiative transfer modeling. We also present a detailed analysis of the detected spots.

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.

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.

Alice is a far-ultraviolet imaging spectrograph onboard Rosetta that, amongst multiple objectives, is designed to observe emissions from various atomic and molecular species from within the coma of comet 67P/Churyumov-Gerasimenko. The initial observations, made following orbit insertion in August 2014, showed emissions of atomic hydrogen and oxygen spatially localized close to the nucleus and attributed to photoelectron impact dissociation of H2O vapor. Weaker emissions from atomic carbon were subsequently detected and also attributed to electron impact dissociation, of CO2, the relative H I and C I line intensities reflecting the variation of CO2 to H2O column abundance along the line-of-sight through the coma. Beginning in mid-April 2015, Alice sporadically observed a number of outbursts above the sunward limb characterized by sudden increases in the atomic emissions, particularly the semi-forbidden O I 1356 multiplet, over a period of 10-30 minutes, without a corresponding enhancement in long wavelength solar reflected light characteristic of dust production. A large increase in the brightness ratio O I 1356/O I 1304 suggests O2 as the principal source of the additional gas. These outbursts do not correlate with any of the visible images of outbursts taken with either OSIRIS or the navigation camera. Beginning in June 2015 the nature of the Alice spectrum changed considerably with CO Fourth Positive band emission observed continuously, varying with pointing but otherwise fairly constant in time. However, CO does not appear to be a major driver of any of the observed outbursts.

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 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.

Aims. The Alice far-ultraviolet spectrograph onboard Rosetta is designed to observe emissions from various atomic and molecular species from within the coma of comet 67P/ Churyumov-Gerasimenko and to determine their spatial distribution and evolution with time and heliocentric distance. Methods. Following orbit insertion in August 2014, Alice made observations of the inner coma above the limbs of the nucleus of the comet from cometocentric distances varying between 10 and 80 km. Depending on the position and orientation of the slit relative to the nucleus, emissions of atomic hydrogen and oxygen were initially detected. These emissions are spatially localized close to the nucleus and spatially variable with a strong enhancement above the comet's neck at northern latitudes. Weaker emission from atomic carbon and CO were subsequently detected. Results. Analysis of the relative line intensities suggests photoelectron impact dissociation of H2O vapor as the source of the observed H I and O I emissions. The electrons are produced by photoionization of H2O. The observed C I emissions are also attributed to electron impact dissociation, of CO2, and their relative brightness to H I reflects the variation of CO2 to H2O column abundance in the coma.

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.

As a follow-up of a recent study, we challenge the claim that the flow of interstellar helium through the solar system has changed substantially over the last decades. We argue that only the IBEX-Lo 2009-2010 measurements are discrepant with older consensus values. Then we show that the probability of the claimed variations of longitude and velocity are highly unlikely (about 1 per cent), in view of the absence of change in latitude and absence of change in the (flow velocity, flow longitude) relation, while random values would be expected. Finally, we report other independent studies showing the stability of Helium flow and the Hydrogen flow over the years 1996-2012, consistent with the seventies earlier determinations of the interstellar flow.

Context: In the near future, new instruments such as ESPRESSO will arrive, allowing us to reach a precision in radial-velocity measurements on the order of 10 cm/s. At this level of precision, several noise sources that until now have been outweighed by photon noise will start to contribute significantly to the error budget. The telluric lines that are not neglected by the masks for the radial velocity computation, here called micro-telluric lines, are one such noise source. Aims: In this work we investigate the impact of micro-telluric lines in the radial velocities calculations. We also investigate how to correct the effect of these atmospheric lines on radial velocities. Methods: The work presented here follows two parallel lines. First, we calculated the impact of the micro-telluric lines by multiplying a synthetic solar-like stellar spectrum by synthetic atmospheric spectra and evaluated the effect created by the presence of the telluric lines. Then, we divided HARPS spectra by synthetic atmospheric spectra to correct for its presence on real data and calculated the radial velocity on the corrected spectra. When doing so, one considers two atmospheric models for the synthetic atmospheric spectra: the LBLRTM and TAPAS. Results: We find that the micro-telluric lines can induce an impact on the radial velocities calculation that can already be close to the current precision achieved with HARPS, and so its effect should not be neglected, especially for future instruments such as ESPRESSO. Moreover, we find that the micro-telluric lines' impact depends on factors, such as the radial velocity of the star, airmass, relative humidity, and the barycentric Earth radial velocity projected along the line of sight at the time of the observation.

We update two kinds of results obtained with the SWAN instrument on board SOHO. First, we use H cell data recorded in 2001 and derive the H flow direction in the same way we performed the study at solar minimum. We compare with the Helium flow and doing so we correct for the coordinate system change between the Ulysses and SOHO mission. The deflection plane we obtain is compatible with the previous result within error bars, confirming the predominant role of the interstellar magnetic field. Secondly, we extend the derivation of solar wind ionization temporal evolution as a function of heliolatitude. The pattern for the present solar minimum is strikingly different from the previous minimum, with a much wider slow solar wind equatorial belt which persists until at least 2008. Comparing with synoptic LASCO/C2 electron densities we infer from a preliminary study that the acceleration of the high speed solar wind occurs at a higher altitude during this minimum, a change expansion models should be able to explain.

We have revisited the series of observations recently used to infer a temporal variation of the interstellar helium flow over the last forty years. Concerning the recent IBEX-Lo direct detection of Helium neutrals, there are two types of precise and unambiguous measurements which do not rely on the exact response of the instrument: the count rate maxima as a function of the spin angle, which determines the ecliptic latitude of the flow, and the count rate maxima as a function of IBEX longitude, which determines a tight relationship between the ecliptic longitude of the flow and its velocity far from the Sun. These measurements provide parameters (and couples of parameters in the second case) remarkably similar to the canonical, old values. In contrast, the preferential choice of a lower velocity and higher longitude reported before from IBEX data is based only on the count rate variation (at each spin phase maximum) as a function of the satellite longitude, when drifting across the region of high fluxes. We have examined the consequences of dead time counting effects, and conclude that their inclusion at a realistic level is sufficient to reconcile the data with the old parameters, calling for further investigations. We discuss the analyses of the STEREO pickup ion (PUI) data and argue that the statistical method that has been preferred to infer the neutral flow longitude (instead of the more direct method based on the PUI maximum flux directions), is not appropriate. Moreover, transport effects may have been significant at the very weak solar activity level of 2007-2009, in which case the longitudes of the PUI maxima are only upper limits on the flow longitude. Finally, we found that the use of some flow longitude determinations based on UV glow data are not adequate. At variance with recent conclusions we find no evidence for a temporal variability of the interstellar helium flow.

The detectors of the UltraViolet Spectrographs (UVS) on Voyager 1/2 are recording a background that was earlier assigned to disintegrations in the RTG. We show that it arises instead from Galactic Cosmic Rays (GCRs). We show the 1992-2013 GCR flux measured by UVS on V1 and, by comparing with data from the GCR dedicated detectors, we estimate the energy range responsible for this UVS signal, around 300 MeV, and the response of UVS to the GCR anisotropy. After the abrupt jumps of May and August 2012 the count rate has been fluctuating only slightly around a constant value, but comparing with data from the LECP and the CRS instruments shows that those small variations are only responses to a varying anisotropy and not to a flux change. Taking advantage of the similarity in energy range to one of the products of the CRS instrument suite, we use the ratio between the two independent signals as a proxy for the temporal evolution of the GCR spectral slope around the 300 MeV range. We show that this slope has remained unchanged since Aug 2012 and find strong evidence that it will no longer vary, implying the end of the modulation at those energies and that V1 at this date is near or past the heliopause. The origin of this narrow and stagnating inner heliosheath is still unclear, and we discuss the potential effects of low solar wind speed episodes and subsequent self-amplified charge-exchange with interstellar neutrals, as a source of deceleration and collapse. We suggest that the quasi-static region encountered by V1 may be related to such effects, triggered by the strong post solar-maximum variability. This did not happen for V2 due to its trajectory at an angle further from the heliosphere axis and a later termination shock crossing. The existence on the upwind side of a mixing layer formed by charge transfer instead of a pure plasma contact discontinuity could explain various observations.

Spectra of astronomical targets acquired from ground-based instruments are affected by the atmospheric transmission. The authors and their institutes are developing a web-based service, TAPAS (Transmissions of the AtmosPhere for AStromomical data). This service, freely available, is developed and maintained within the atmospheric ETHER data center. TAPAS computes the atmospheric transmission in the line-of-sight to the target indicated by the user. The user files a request indicating the time, ground location, and either the equatorial coordinates of the target or the Zenith Angle of the line-of sight (LOS). The actual atmospheric profile (temperature, pressure, humidity, ozone content) at that time and place is retrieved from the ETHER atmospheric data base (from a combination of ECMWF meteorological field and other informations), and the atmospheric transmission is computed from LBLRTM software and HITRAN data base for a number of gases: O2, H2O, O3, CO2, and Rayleigh extinction. The first purpose of TAPAS output is to allow identification of observed spectral features as being from atmospheric or astrophysical origin. The returned transmission may also serve for characterizing the spectrometer in wavelength scale and Instrument Line Spectral Function (ILSF) by comparing one observed spectrum of an atmospheric feature to the transmission. Finally, the TOA (Top Of Atmosphere) spectrum may be obtained either by division of the observed spectrum by the computed transmission or other techniques developed on purpose. The obtention of transmissions for individual species allows more potentialities and better adjustments to the data. In this paper, we describe briefly the mechanism of computation of the atmospheric transmissions, and we show some results for O2 and H2O atmospheric absorption. The wavelength range is presently 500-2500 nm. The address is http://ether.ipsl.jussieu.fr/tapas/

(Abridged) Searching for planets around stars with different masses probes the outcome of planetary formation for different initial conditions. This drives observations of a sample of 102 southern nearby M dwarfs, using a fraction of our guaranteed time on the ESO/HARPS spectrograph (Feb. 11th, 2003 to Apr. 1st 2009). This paper makes available the sample's time series, presents their precision and variability. We apply systematic searches and diagnostics to discriminate whether the observed Doppler shifts are caused by stellar surface inhomogeneities or by the radial pull of orbiting planets. We recover the planetary signals corresponding to 9 planets already announced by our group (Gl176b, Gl581b, c, d & e, Gl674b, Gl433b, Gl 667Cb and c). We present radial velocities that confirm GJ 849 hosts a Jupiter-mass planet, plus a long-term radial-velocity variation. We also present RVs that precise the planetary mass and period of Gl 832b. We detect long-term RV changes for Gl 367, Gl 680 and Gl 880 betraying yet unknown long-period companions. We identify candidate signals in the radial-velocity time series and demonstrate they are most probably caused by stellar surface inhomogeneities. Finally, we derive a first estimate of the occurrence of M-dwarf planets as a function of their minimum mass and orbital period. In particular, we find that giant planets (m sin i = 100-1,000 Mearth) have a low frequency (e.g. f<1% for P=1-10 d and f=0.02^{+0.03}_{-0.01} for P=10-100 d), whereas super-Earths (m sin i = 1-10 Mearth) are likely very abundant (f=0.36^{+0.25}_{-0.10} for P=1-10 d and f=0.35^{+0.45}_{-0.11} for P=10-100 d). We also obtained eta_earth=0.41^{+0.54}_{-0.13}, the frequency of habitable planets orbiting M dwarfs (1<m sin i<10 Mearth). For the first time, eta_earth is a direct measure and not a number extrapolated from the statistic of more massive and/or shorter-period planets.

We report on the results of an 8-year survey carried out at the La Silla Observatory with the HARPS spectrograph to detect and characterize planets in the super-Earth and Neptune mass regime. The size of our star sample and the precision achieved with HARPS have allowed the detection of a sufficiently large number of low-mass planets to study the statistical properties of their orbital elements, the correlation of the host-star metallicity with the planet masses, as well as the occurrence rate of planetary systems around solar-type stars. A robust estimate of the frequency of systems shows that more than 50% of solar-type stars harbor at least one planet of any mass and with period up to 100 days. Different properties are observed for the population of planets less massive than about 30M-Earth compared to the population of gaseous giant planets. The mass distribution of Super-Earths and Neptune-mass planets (SEN) is strongly increasing between 30 and 15M-Earth. The SEN occurence rate does not exhibit a preference for metal rich stars. Most of the SEN planets belong to multi-planetary systems. The orbital eccentricities of the SEN planets seems limited to 0.45. At the opposite, the occurence rate of gaseous giant planets is growing with the logarithm of the period, and is strongly increasing with the host-star metallicity. About 14% of solar-type stars have a planetary companion more massive than 50M-Earth? on an orbit with a period shorter than 10 years. Orbital eccentricities of giant planets are observed up to 0.9 and beyond. The precision of HARPS-type spectrographs opens the possibility to detect planets in the habitable zone of solar-type stars. Identification of a significant number of super-Earths orbiting solar-type of the Sun vicinity is achieved by Doppler spectroscopy. 37 newly discovered planets are announced in the Appendix of this paper, among which 15 Super-Earths.

The European Space Agency's Rosetta spacecraft, en route to a 2014 encounter with comet 67P/Churyumov-Gerasimenko, made a gravity assist swing-by of Mars on 25 February 2007, closest approach being at 01:54UT. The Alice instrument on board Rosetta, a lightweight far-ultraviolet imaging spectrograph optimized for in situ cometary spectroscopy in the 750-2000 A spectral band, was used to study the daytime Mars upper atmosphere including emissions from exospheric hydrogen and oxygen. Offset pointing, obtained five hours before closest approach, enabled us to detect and map the HI Lyman-alpha and Lyman-beta emissions from exospheric hydrogen out beyond 30,000 km from the planet's center. These data are fit with a Chamberlain exospheric model from which we derive the hydrogen density at the 200 km exobase and the H escape flux. The results are comparable to those found from the the Ultraviolet Spectrometer experiment on the Mariner 6 and 7 fly-bys of Mars in 1969. Atomic oxygen emission at 1304 A is detected at altitudes of 400 to 1000 km above the limb during limb scans shortly after closest approach. However, the derived oxygen scale height is not consistent with recent models of oxygen escape based on the production of suprathermal oxygen atoms by the dissociative recombination of O2+.

Global water production rates were determined from the Lyman-{\alpha} emission of hydrogen around comet 103P/Hartley 2, observed with the SWAN (Solar Wind ANisotropies) all-sky camera on the SOHO spacecraft from September 14 through December 12, 2010. This time period included the November 4 flyby by the EPOXI spacecraft. Water production was 3 times lower than during the 1997 apparition also measured by SWAN. In 2010 it increased by a factor of ~2.5 within one day on September 30 with a similar corresponding drop between November 24 and 30. The total surface area of sublimating water within {\pm}20 days of perihelion was ~0.5 km^2, about half of the mean cross section of the nucleus. Outside this period it was ~0.2 km^2. The peak value was 90%, implying a significant water production by released nucleus icy fragments.

Context. Low-mass extrasolar planets are presently being discovered at an increased pace by radial velocity and transit surveys, opening a new window on planetary systems. Aims. We are conducting a high-precision radial velocity survey with the HARPS spectrograph which aims at characterizing the population of ice giants and super-Earths around nearby solar-type stars. This will lead to a better understanding of their formation and evolution, and yield a global picture of planetary systems from gas giants down to telluric planets. Methods. Progress has been possible in this field thanks in particular to the sub-m/s radial velocity precision achieved by HARPS. We present here new high-quality measurements from this instrument. Results. We report the discovery of a planetary system comprising at least five Neptune-like planets with minimum masses ranging from 12 to 25 M_Earth, orbiting the solar-type star HD 10180 at separations between 0.06 and 1.4 AU. A sixth radial velocity signal is present at a longer period, probably due to a 65-M_Earth object. Moreover, another body with a minimum mass as low as 1.4 M_Earth may be present at 0.02 AU from the star. This is the most populated exoplanetary system known to date. The planets are in a dense but still well-separated configuration, with significant secular interactions. Some of the orbital period ratios are fairly close to integer or half-integer values, but the system does not exhibit any mean-motion resonances. General relativity effects and tidal dissipation play an important role to stabilize the innermost planet and the system as a whole. Numerical integrations show long-term dynamical stability provided true masses are within a factor ~3 from minimum masses. We further note that several low-mass planetary systems exhibit a rather "packed" orbital architecture with little or no space left for additional planets. (Abridged)

We report the detection of three new extrasolar planets orbiting the solar type stars HD 85390, HD 90156 and HD 103197 with the HARPS spectrograph mounted on the ESO 3.6-m telescope at La Silla observatory. HD 85390 has a planetary companion with a projected intermediate mass (42.0 Earth masses) on a 788-day orbit (a=1.52 AU) with an eccentricity of 0.41, for which there is no analogue in the solar system. A drift in the data indicates the presence of another companion on a long period orbit, which is however not covered by our measurements. HD 90156 is orbited by a warm Neptune analogue with a minimum mass of 17.98 Earth masses (1.05 Neptune masses), a period of 49.8 days (a=0.25 AU) and an eccentricity of 0.31. HD 103197 has an intermediate mass planet on a circular orbit (P=47.8 d, Msini=31.2 Earth masses). We discuss the formation of planets of intermediate mass (about 30-100 Earth masses) which should be rare inside a few AU according to core accretion formation models.

The asteroid (21) Lutetia is the target of a planned close encounter by the Rosetta spacecraft in July 2010. To prepare for that flyby, Lutetia has been extensively observed by a variety of astronomical facilities. We used the Hubble Space Telescope (HST) to determine the albedo of Lutetia over a wide wavelength range, extending from ~150 nm to ~700 nm. Using data from a variety of HST filters and a ground-based visible light spectrum, we employed synthetic photometry techniques to derive absolute fluxes for Lutetia. New results from ground-based measurements of Lutetia's size and shape were used to convert the absolute fluxes into albedos. We present our best model for the spectral energy distribution of Lutetia over the wavelength range 120-800 nm. There appears to be a steep drop in the albedo (by a factor of ~2) for wavelengths shorter than ~300 nm. Nevertheless, the far ultraviolet albedo of Lutetia (~10%) is considerably larger than that of typical C-chondrite material (~4%). The geometric albedo at 550 nm is 16.5 +/- 1%. Lutetia's reflectivity is not consistent with a metal-dominated surface at infrared or radar wavelengths, and its albedo at all wavelengths (UV-visibile-IR-radar) is larger than observed for typical primitive, chondritic material. We derive a relatively high FUV albedo of ~10%, a result that will be tested by observations with the Alice spectrograph during the Rosetta flyby of Lutetia in July 2010.

We update the SOHO/SWAN H Lyman-alpha brightness analysis to cover the 1996-2008 time interval. A forward model applied to the intensity maps provides the latitude and time dependence of the interstellar Hydrogen ionisation rate over more than a full solar cycle. The hydrogen ionisation, being almost entirely due to charge-exchange with solar wind ions, reflects closely the solar wind flux. Our results show that the solar wind latitudinal structure during the present solar minimum is strikingly different from the previous minimum, with a much wider slow solar wind equatorial belt which persists until at least the end of 2008. We compute absolute values of the in-ecliptic H ionisation rates using OMNI solar wind data and use them to calibrate our ionisation rates at all heliographic latitudes. We then compare the resulting fluxes with the synoptic LASCO/C2 electron densities at 6 solar radii. The two time-latitude patterns are strikingly similar over all the cycle. This comparison shows that densities at 6 solar radii can be used to infer the solar wind type close to its source, with high (resp. low) densities tracing the slow (resp. fast) solar wind, simply because the density reflects at which altitude occurs the acceleration. The comparison between the two minima suggests that the fast polar wind acceleration occurs at larger distance during the current minimum compared to the previous one. This difference, potentially linked to the magnetic field decrease or(and) the coronal temperature decrease should be reproduced by solar wind expansion models.

We have analyzed data from a multi-site campaign to observe oscillations in the F5 star Procyon. The data consist of high-precision velocities that we obtained over more than three weeks with eleven telescopes. A new method for adjusting the data weights allows us to suppress the sidelobes in the power spectrum. Stacking the power spectrum in a so-called echelle diagram reveals two clear ridges that we identify with even and odd values of the angular degree (l=0 and 2, and l=1 and 3, respectively). We interpret a strong, narrow peak at 446 muHz that lies close to the l=1 ridge as a mode with mixed character. We show that the frequencies of the ridge centroids and their separations are useful diagnostics for asteroseismology. In particular, variations in the large separation appear to indicate a glitch in the sound-speed profile at an acoustic depth of about 1000 s. We list frequencies for 55 modes extracted from the data spanning 20 radial orders, a range comparable to the best solar data, which will provide valuable constraints for theoretical models. A preliminary comparison with published models shows that the offset between observed and calculated frequencies for the radial modes is very different for Procyon than for the Sun and other cool stars. We find the mean lifetime of the modes in Procyon to be 1.29 +0.55/-0.49 days, which is significantly shorter than the 2-4 days seen in the Sun.

Precise radial-velocity measurements for data acquired with the HARPS spectrograph infer that three planets orbit the M4 dwarf star GJ876. In particular, we confirm the existence of planet "d", which orbits every 1.93785 days. We find that its orbit may have significant eccentricity (e=0.14), and deduce a more accurate estimate of its minimum mass of 6.3 Earth masses. Dynamical modeling of the HARPS measurements combined with literature velocities from the Keck Observatory strongly constrain the orbital inclinations of the "b" and "c" planets. We find that i_b = 48.9 degrees and i_c = 48.1 degrees, which infers the true planet masses of M_b = 2.64 Jupiter masses and M_c = 0.83 Jupiter masses, respectively. Radial velocities alone, in this favorable case, can therefore fully determine the orbital architecture of a multi-planet system, without the input from astrometry or transits. The orbits of the two giant planets are nearly coplanar, and their 2:1 mean motion resonance ensures stability over at least 5 Gyr. The libration amplitude is smaller than 2 degrees, suggesting that it was damped by some dissipative process during planet formation. The system has space for a stable fourth planet in a 4:1 mean motion resonance with planet "b", with a period around 15 days. The radial velocity measurements constrain the mass of this possible additional planet to be at most that of the Earth.

The GJ 581 planetary system was already known to harbour three planets, including two super-Earths planets which straddle its habitable zone. We report here the detection of an additional planet -- GJ 581e -- with a minimum mass of 1.9 M_earth. With a period of 3.15 days, it is the innermost planet of the system and has a ~5% transit probability. We also correct our previous confusion of the orbital period of GJ 581d (the outermost planet) with a one-year alias, thanks to an extended time span and many more measurements. The revised period is 66.8 days, and locates the semi-major axis inside the habitable zone of the low mass star. The dynamical stability of the 4-planet system imposes an upper bound on the orbital plane inclination. The planets cannot be more massive than approximately 1.6 times their minimum mass.

Precise radial-velocity measurements with the HARPS spectrograph reveal the presence of two planets orbiting the solar-type star HD45364. The companion masses are 0.187 Mjup and 0.658 Mjup, with semi-major axes of 0.681 AU and 0.897 AU, and eccentricities of 0.168 and 0.097, respectively. A dynamical analysis of the system further shows a 3:2 mean motion resonance between the two planets, which prevents close encounters and ensures the stability of the system over 5 Gyr. This is the first time that such a resonant configuration has been observed for extra-solar planets, although there is an analogue in our Solar System formed by Neptune and Pluto. This singular planetary system may provide important constraints on planetary formation and migration scenarios.

This paper reports on the detection of two new multiple planet systems around solar-like stars HD47186 and HD181433. The first system includes a hot Neptune of 22.78 M_Earth at 4.08-days period and a Saturn of 0.35 M_Jup at 3.7-years period. The second system includes a Super-Earth of 7.5 M_Earth at 9.4-days period, a 0.64 M$_Jup at 2.6-years period as well as a third companion of 0.54 M_Jup with a period of about 6 years. These detections increase to 20 the number of close-in low-mass exoplanets (below 0.1 M_Jup) and strengthen the fact that 80% of these planets are in a multiple planetary systems.

A 10.24 days Neptune-mass planet was recently announced to orbit the nearby M2 dwarf Gl 176, based on 28 radial velocities measured with the HRS spectrograph on the Hobby-Heberly Telescope (HET). We obtained 57 radial velocities of Gl 176 with the ESO 3.6m telescope and the HARPS spectrograph, which is known for its sub-m/s stability. The median photon-noise standard error of our measurements is 1.1 m/s, significantly lower than the 4.7 m/s of the HET velocities, and the 4 years period over which they were obtained has much overlap with the epochs of the HET measurements. The HARPS measurements show no evidence for a signal at the period of the putative HET planet, suggesting that its detection was spurious. We do find, on the other hand, strong evidence for a lower mass 8.4 Mearth planet, in a quasi-circular orbit and at the different period of 8.78 days. The host star has moderate magnetic activity and rotates on a 39-days period, which we confirm through modulation of both contemporaneous photometry and chromospheric indices. We detect that period as well in the radial velocities, but it is well removed from the orbital period and no cause for confusion. This new detection of a super-Earth (2 Mearth < M sin(i) < 10 Mearth) around an M dwarf adds to the growing evidence that such planets are common around very low mass stars: a third of the 20 known planets with M sin(i) < 0.1 Mjup and 3 of the 7 known planets with M sin(i) < 10 Mearth orbit an M dwarf, in contrast to just 4 of the ~300 known Jupiter-mass planets.

We have carried out a multi-site campaign to measure oscillations in the F5 star Procyon A. We obtained high-precision velocity observations over more than three weeks with eleven telescopes, with almost continuous coverage for the central ten days. This represents the most extensive campaign so far organized on any solar-type oscillator. We describe in detail the methods we used for processing and combining the data. These involved calculating weights for the velocity time series from the measurement uncertainties and adjusting them in order to minimize the noise level of the combined data. The time series of velocities for Procyon shows the clear signature of oscillations, with a plateau of excess power that is centred at 0.9 mHz and is broader than has been seen for other stars. The mean amplitude of the radial modes is 38.1 +/- 1.3 cm/s (2.0 times solar), which is consistent with previous detections from the ground and by the WIRE spacecraft, and also with the upper limit set by the MOST spacecraft. The variation of the amplitude during the observing campaign allows us to estimate the mode lifetime to be 1.5 d (+1.9/-0.8 d). We also find a slow variation in the radial velocity of Procyon, with good agreement between different telescopes. These variations are remarkably similar to those seen in the Sun, and we interpret them as being due to rotational modulation from active regions on the stellar surface. The variations appear to have a period of about 10 days, which presumably equals the stellar rotation period or, perhaps, half of it. The amount of power in these slow variations indicates that the fractional area of Procyon covered by active regions is slightly higher than for the Sun.

This paper reports on the detection of a planetary system with three Super-Earths orbiting HD40307. HD40307 is a K2V metal-deficient star at a distance of only 13 parsec, part of the HARPS GTO high-precision planet-search programme. The three planets on circular orbits have very low minimum masses of respectively 4.2, 6.9 and 9.2 Earth masses and periods of 4.3, 9.6 and 20.5 days. The planet with the shortest period is the lightest planet detected to-date orbiting a main sequence star. The detection of the correspondingly low amplitudes of the induced radial-velocity variations is completely secured by the 135 very high-quality HARPS observations illustrated by the radial-velocity residuals around the 3-Keplerian solution of only 0.85 m/s. Activity and bisector indicators exclude any significant perturbations of stellar intrinsic origin, which supports the planetary interpretation. Contrary to most planet-host stars, HD40307 has a marked sub-solar metallicity ([Fe/H]=-0.31), further supporting the already raised possibility that the occurrence of very light planets might show a different dependence on host star's metallicity compared to the population of gas giant planets. In addition to the 3 planets close to the central star, a small drift of the radial-velocity residuals reveals the presence of another companion in the system the nature of which is still unknown.

This Letter reports on the detection of two super-Earth planets in the Gl581 system, already known to harbour a hot Neptune. One of the planets has a mass of 5 M_Earth and resides at the ``warm'' edge of the habitable zone of the star. It is thus the known exoplanet which most resembles our own Earth. The other planet has a 7.7 M_Earth mass and orbits at 0.25 AU from the star, close to the ``cold'' edge of the habitable zone. These two new light planets around an M3 dwarf further confirm the formerly tentative statistical trend for i) many more very low-mass planets being found around M dwarfs than around solar-type stars and ii) low-mass planets outnumbering Jovian planets around M dwarfs.

Context: How planet properties depend on stellar mass is a key diagnostic of planetary formation mechanisms. Aims: This motivates planet searches around stars which are significantly more massive or less massive than the Sun, and in particular our radial velocity search for planets around very-low mass stars. Methods: As part of that program, we obtained measurements of GJ 674, an M2.5 dwarf at d=4.5 pc, which have a dispersion much in excess of their internal errors. An intensive observing campaign demonstrates that the excess dispersion is due to two superimposed coherent signals, with periods of 4.69 and 35 days. Results: These data are well described by a 2-planet Keplerian model where each planet has a ~11 Mearth minimum mass. A careful analysis of the (low level) magnetic activity of GJ 674 however demonstrates that the 35-day period coincides with the stellar rotation period. This signal therefore originates in a spot inhomogeneity modulated by stellar rotation. The 4.69-day signal on the other hand is caused by a bona-fide planet, GJ 674b. Conclusion: Its detection adds to the growing number of Neptune-mass planets around M-dwarfs, and reinforces the emerging conclusion that this mass domain is much more populated than the jovian mass range. We discuss the metallicity distributions of M dwarf with and without planets and find a low 11% probability that they are drawn from the same parent distribution. Moreover, we find tentative evidence that the host star metallicity correlates with the total mass of their planetary system.

Over the past two years, the search for low-mass extrasolar planets has led to the detection of seven so-called 'hot Neptunes' or 'super-Earths' around Sun-like stars. These planets have masses 5-20 times larger than the Earth and are mainly found on close-in orbits with periods of 2-15 days. Here we report a system of three Neptune-mass planets with periods of 8.67, 31.6 and 197 days, orbiting the nearby star HD 69830. This star was already known to show an infrared excess possibly caused by an asteroid belt within 1 AU (the Sun-Earth distance). Simulations show that the system is in a dynamically stable configuration. Theoretical calculations favour a mainly rocky composition for both inner planets, while the outer planet probably has a significant gaseous envelope surrounding its rocky/icy core; the outer planet orbits within the habitable zone of this star.

Interplanetary Lyman alpha line profiles are derived from the SWAN H cell data measurements. The measurements cover a 6-year period from solar minimum (1996) to after the solar maximum of 2001. This allows us to study the variations of the line profiles with solar activity. These line profiles were used to derive line shifts and line widths in the interplanetary medium for various angles of the LOS with the interstellar flow direction. The SWAN data results were then compared to an interplanetary background upwind spectrum obtained by STIS/HST in March 2001. We find that the LOS upwind velocity associated with the mean line shift of the IP \lya line varies from 25.7 km/s to 21.4 km/s from solar minimum to solar maximum. Most of this change is linked with variations in the radiation pressure. LOS kinetic temperatures derived from IP line widths do not vary monotonically with the upwind angle of the LOS. This is not compatible with calculations of IP line profiles based on hot model distributions of interplanetary hydrogen. We also find that the line profiles get narrower during solar maximum. The results obtained on the line widths (LOS temperature) show that the IP line is composed of two components scattered by two hydrogen populations with different bulk velocities and temperature. This is a clear signature of the heliospheric interface on the line profiles seen at 1 AU from the sun.

We report on spectroscopic observations of periodic comet 9P/Tempel 1 by the Alice ultraviolet spectrograph on the Rosetta spacecraft in conjunction with NASA's Deep Impact mission. Our objectives were to measure an increase in atomic and molecular emissions produced by the excavation of volatile sub-surface material. We unambiguously detected atomic oxygen emission from the quiescent coma but no enhancement at the 10% (1-sigma) level following the impact. We derive a quiescent water production rate of 9 x 10^27 molecules per second with an estimated uncertainty of 30%. Our upper limits to the volatiles produced by the impact are consistent with other estimates.

The mu Ara planetary system is rather complex: It contains two already known planets, mu Ara b with P=640 days and mu Ara c with P=9.64 days, and a third companion on a wide but still poorly defined orbit. Even with three planets in the system, the data points keep anomalously high dispersion around the fitted solution. The high residuals are only partially due to the strong p-mode oscillations of the host star. We have therefore studied in this paper the possible presence of a fourth planet in the system. During the past years we have carried out additional and extremely precise radial-velocity measurements with the HARPS spectrograph. We provide in this paper a full orbital solution of the planetary system around mu Ara. It turns out to be the second system known to harbor 4 planetary companions. Thanks to the new data points acquired with HARPS we can confirm the presence of mu Ara c at P=9.64 days, which produces a coherent RV signal over more than two years. The new orbital fit sets the mass of mu Ara c to 10.5 M_Earth. Furthermore, we present the discovery of mu Ara d, a new planet on an almost circular 310 days-period and with a mass of 0.52 M_Jup. Finally, we give completely new orbital parameters for the longest-period planet, mu Ara e. It is the first time that this companion is constrained by radial-velocity data into a dynamical stable orbit, which leaves no doubt about its planetary nature. (Abridged).

We describe the design, performance and scientific objectives of the NASA-funded ALICE instrument aboard the ESA Rosetta asteroid flyby/comet rendezvous mission. ALICE is a lightweight, low-power, and low-cost imaging spectrograph optimized for cometary far-ultraviolet (FUV) spectroscopy. It will be the first UV spectrograph to study a comet at close range. It is designed to obtain spatially-resolved spectra of Rosetta mission targets in the 700-2050 A spectral band with a spectral resolution between 8 A and 12 A for extended sources that fill its ~0.05 deg x 6.0 deg field-of-view. ALICE employs an off-axis telescope feeding a 0.15-m normal incidence Rowland circle spectrograph with a concave holographic reflection grating. The imaging microchannel plate detector utilizes dual solar-blind opaque photocathodes (KBr and CsI) and employs a 2 D delay-line readout array. The instrument is controlled by an internal microprocessor. During the prime Rosetta mission, ALICE will characterize comet 67P/Churyumov-Gerasimenko's coma, its nucleus, and the nucleus/coma coupling; during cruise to the comet, ALICE will make observations of the mission's two asteroid flyby targets and of Mars, its moons, and of Earth's moon. ALICE has already successfully completed the in-flight commissioning phase and is operating normally in flight. It has been characterized in flight with stellar flux calibrations, observations of the Moon during the first Earth fly-by, and observations of comet Linear T7 in 2004 and comet 9P/Tempel 1 during the 2005 Deep Impact comet-collision observing campaign

We present here the discovery and characterisation of a very light planet around HD4308. The planet orbits its star in 15.56 days. The circular radial-velocity variation presents a tiny semi-amplitude of 4.1 m/s that corresponds to a planetary minimum mass m2sin(i)=14.1 Earth masses. The planet was unveiled by high-precision radial-velocity measurements obtained with the HARPS spectrograph on the ESO 3.6-m telescope. The radial-velocity residuals around the Keplerian solution are 1.3 m/s, demonstrating the very high quality of the HARPS measurements. Activity and bisector indicators exclude any significant perturbations of stellar intrinsic origin, which supports the planetary interpretation. Contrary to most planet-host stars, HD4308 has a marked sub-solar metallicity ([Fe/H]=-0.31), raising the possibility that very light planet occurrence might show a different coupling with the parent star's metallicity than do giant gaseous extra-solar planets. Together with Neptune-mass planets close to their parent stars, the new planet occupies a position in the mass-separation parameter space that is constraining for planet-formation and evolution theories. The question of whether they can be considered as residuals of evaporated gaseous giant planets, ice giants, or super-earth planets is discussed in the context of the latest core-accretion models.

We report the discovery of a Neptune-mass planet around Gl 581 (M3V, M = 0.31 Msol), based on precise Doppler measurements with the HARPS spectrograph at La Silla Observatory. The radial velocities reveal a circular orbit of period P = 5.366 days and semi-amplitude K1 = 13.2 m/s. The resulting minimum mass of the planet (m2 sin i) is only 0.052 Mjup = 0.97 Mnep = 16.6 Mearth making Gl 581b one of the lightest extra-solar planet known to date. The Gl 581 planetary system is only the third centered on an M dwarf, joining the Gl 876 three-planet system and the lone planet around Gl 436. Its discovery reinforces the emerging tendency of such planets to be of low mass, and found at short orbital periods. The statistical properties of the planets orbiting M dwarfs do not seem to match a simple mass scaling of their counterparts around solar-type stars.

We report on the detection of three Saturn-mass planets discovered with the HARPS instrument. HD 93083 shows radial-velocity (RV) variations best explained by the presence of a companion of 0.37 M_Jup orbiting in 143.6 days. HD 101930 b has an orbital period of 70.5 days and a minimum mass of 0.30 M_Jup. For HD 102117, we present the independent detection of a companion with m2sini = 0.14 M_Jup and orbital period P = 20.7 days. This planet was recently detected by Tinney et al. (2004). Activity and bisector indicators exclude any significant RV perturbations of stellar origin, reinforcing the planetary interpretation of the RV variations. The radial-velocity residuals around the Keplerian fits are 2.0, 1.8 and 0.9 m/s respectively, showing the unprecedented RV accuracy achieved with HARPS. A sample of stable stars observed with HARPS is also presented to illustrate the long-term precision of the instrument. All three stars are metal-rich, confirming the now well-established relation between planet occurrence and metallicity. The new planets are all in the Saturn-mass range, orbiting at moderate distance from their parent star, thereby occupying an area of the parameter space which seems difficult to populate according to planet formation theories. A systematic exploration of these regions will provide new constraints on formation scenarios in the near future.

In this letter we present the discovery of a very light planetary companion to the star mu Ara (HD160691). The planet orbits its host once every 9.5days, and induces a sinusoidal radial velocity signal with a semi-amplitude of 4.1 m/s, the smallest Doppler amplitude detected so far. These values imply a mass of m2 sini = 14 earth-masses. This detection represents the discovery of a planet with a mass slightly smaller than that of Uranus, the smallest ``ice giant" in our Solar System. Whether this planet can be considered an ice giant or a super-earth planet is discussed in the context of the core-accretion and migration models.