Permittivity measurements of porous matter in support of investigations of the surface and interior of 67P/Churyumov-Gerasimenko

Aims. Permittivity measurements on porous samples of volcanic origin have been performed in the 0.05-190 GHz range under laboratory conditions in support of the Rosetta mission to comet 67P/Churyumov-Gerasimenko, specifically with the MIRO radiometric experiment and CONSERT radar experiment. Methods. The samples were split into several subsamples with different size ranges covering a few mu m to 500 mu m. Bulk densities of the subsamples were estimated to be in the 800 to 1500 kg/m(3) range. The porosities were in the range of 48% to 65%. From 50 MHz to 6 GHz and at 190 GHz, permittivity has been determined with a coaxial cell and with a quasi-optical bench, respectively. Results. Without taking into account the volume-scattering effect at 190 GHz, the real part of the permittivity, normalized by the bulk density, is in the range of 2.1 to 2.6. The results suggest that the real part of the permittivity of an ice-free dust mantle covering the nucleus is in the 1.5-2.2 range at 190 GHz. From these values, a lower limit for the absorption length for the millimeter receiver of MIRO has been estimated to be between 0.6 and 2 cm, in agreement with results obtained from MIRO in September 2014. At frequencies of interest for CONSERT experiment, the real part of the permittivity of a suspected ice-free dust mantle should be below 2.2. It may be in the range of 1.2 to 1.7 for the nucleus, in agreement with first CONSERT results, taking into account a mean temperature of 110 K and different values for the dust-to-ice volumetric ratio. Estimations of contributions of the different parameters to the permittivity variation may indicate that the porosity is the main parameter.

Orbital elements of the material surrounding comet 67P/Churyumov-Gerasimenko

Context. We investigate the dust coma within the Hill sphere of comet 67P/Churyumov-Gerasimenko. Aims. We aim to determine osculating orbital elements for individual distinguishable but unresolved slow-moving grains in the vicinity of the nucleus. In addition, we perform photometry and constrain grain sizes. Methods. We performed astrometry and photometry using images acquired by the OSIRIS Wide Angle Camera on the European Space Agency spacecraft Rosetta. Based on these measurements, we employed standard orbit determination and orbit improvement techniques. Results. Orbital elements and effective diameters of four grains were constrained, but we were unable to uniquely determine them. Two of the grains have light curves that indicate grain rotation. Conclusions. The four grains have diameters nominally in the range 0.14-0.50 m. For three of the grains, we found elliptic orbits, which is consistent with a cloud of bound particles around the nucleus. However, hyperbolic escape trajectories cannot be excluded for any of the grains, and for one grain this is the only known option. One grain may have originated from the surface shortly before observation. These results have possible implications for the understanding of the dispersal of the cloud of bound debris around comet nuclei, as well as for understanding the ejection of large grains far from the Sun.

Can we constrain the interior structure of rocky exoplanets from mass and radius measurements?

Aims. We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. Methods. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated uncertainties. Results. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii, and to exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision required to characterize the interior of exoplanets. Conclusions. Our main conclusions are (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are principal constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement accuracies, but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abundances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. We provide a general methodology of analyzing interior structures of exoplanets that may help to understand how interior models are distributed among star systems. The methodology we propose is sufficiently general to allow its future extension to more complex internal structures including hydrogen- and water-rich exoplanets.

Regional surface morphology of comet 67P/Churyumov-Gerasimenko from Rosetta/OSIRIS images

Aims. The OSIRIS camera onboard the Rosetta spacecraft has been acquiring images of the comet 67P/Churyumov-Gerasimenko (67P)'s nucleus at spatial resolutions down to similar to 0.17 m/px ever since Aug. 2014. These images have yielded unprecedented insight into the morphological diversity of the comet's surface. This paper presents an overview of the regional morphology of comet 67P. Methods. We used the images that were acquired at orbits similar to 20-30 km from the center of the comet to distinguish different regions on the surface and introduce the basic regional nomenclature adopted by all papers in this Rosetta special feature that address the comet's morphology and surface processes. We used anaglyphs to detect subtle regional and topographical boundaries and images from close orbit (similar to 10 km from the comet's center) to investigate the fine texture of the surface. Results. Nineteen regions have currently been defined on the nucleus based on morphological and/or structural boundaries, and they can be grouped into distinctive region types. Consolidated, fractured regions are the most common region type. Some of these regions enclose smooth units that appear to settle in gravitational sinks or topographically low areas. Both comet lobes have a significant portion of their surface covered by a dusty coating that appears to be recently placed and shows signs of mobilization by aeolian-like processes. The dusty coatings cover most of the regions on the surface but are notably absent from a couple of irregular large depressions that show sharp contacts with their surroundings and talus-like deposits in their interiors, which suggests that short-term explosive activity may play a significant role in shaping the comet's surface in addition to long-term sublimation loss. Finally, the presence of layered brittle units showing signs of mechanical failure predominantly in one of the comet's lobes can indicate a compositional heterogeneity between the two lobes.

Spectrophotometric properties of the nucleus of comet 67P/Churyumov-Gerasimenko from the OSIRIS instrument onboard the ROSETTA spacecraft

Context. 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. Aims. 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 observations presented here were performed during July and the beginning of August 2014, during the approach phase, when OSIRIS was mapping the surface of the comet with several filters at different phase angles (1.3 degrees-54 degrees). The resolution reached up to 2.1 m/px. Methods. The OSIRIS images were processed with the OSIRIS standard pipeline, then converted into I/F. radiance factors and corrected for the illumination conditions at each pixel using the Lommel-Seeliger disk law. Color cubes of the surface were produced by stacking registered and illumination-corrected images. Furthermore, photometric analysis was performed both on disk-averaged photometry in several filters and on disk-resolved images acquired with the NAC orange filter, centered at 649 ran, using Hapke modeling. Results. The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13 +/- 0.01 in the HG system formalism and an absolute magnitude H-v(1, 1, 0) = 15.74 +/- 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 similar to 290 rim that is possibly due to SO2 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 degrees-54 degrees phase angle range. The geometric albedo of the comet is 6.5 +/- 0.2% at 649 nm, with local variations of up to similar to 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.

Gravitational slopes, geomorphology, and material strengths of the nucleus of comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Aims. 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 (tensile, shear, and compressive strengths). Methods. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface (Imhotep, Ash, Seth, Hathor, and Agilkia), 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 (overhangs, collapsed structures, boulders, cliffs, and Philae's footprint) and mechanical considerations. Results. 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 degrees) are covered by a fine material and contain a few large (>10 m) and isolated boulders; ii) intermediate-slope terrains (20-45 degrees) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from <1m to 10m for the majority of them; and iii) high-slope terrains (45-90 degrees) 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.

Temporal morphological changes in the Imhotep region of comet 67P/Churyumov-Gerasimenko

Aims. 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. Methods. We used images of the OSIRIS cameras onboard Rosetta to follow the temporal changes from 24 May 2015 to 11 July 2015. Results. The morphological changes observed on the surface are visible in the form of roundish features that are growing in size from a given location in a preferential direction at a rate of 5.6-8.1 x 10(-5) m s(-1) during the observational period. The location where the changes started and the contours of the expanding features are bluer than the surroundings, which suggests that ices (H2O and/or CO2) are 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.

Insolation, erosion, and morphology of comet 67P/Churyumov-Gerasimenko

Context. The complex shape of comet 67P and its oblique rotation axis cause pronounced seasonal effects. Irradiation and hence activity vary strongly. Aims. We investigate the insolation of the cometary surface in order to predict the sublimation of water ice. The strongly varying erosion levels are correlated with the topography and morphology of the present cometary surface and its evolution. Methods. The insolation as a function of heliocentric distance and diurnal (spin dependent) variation is calculated using >10(5) facets of a detailed digital terrain model. Shading, but also illumination and thermal radiation by facets in the field of view of a specific facet are iteratively taken into account. We use a two-layer model of a thin porous dust cover above an icy surface to calculate the water sublimation, presuming steady state and a uniform surface. Our second model, which includes the history of warming and cooling due to thermal inertia, is restricted to a much simpler shape model but allows us to test various distributions of active areas. Results. Sublimation from a dirty ice surface yields maximum erosion. A thin dust cover of 50 pm yields similar rates at perihelion. Only about 6% of the surface needs to be active to match the observed water production rates at perihelion. A dust layer of 1 mm thickness suppresses the activity by a factor of 4 to 5. Erosion on the south side can reach more than 10 m per orbit at active spots. The energy input to the concave neck area (Hapi) during northern summer is enhanced by about 50% owing to self-illumination. Here surface temperatures reach maximum values along the foot of the Hathor wall. Integrated over the whole orbit this area receives the least energy input. Based on the detailed shape model, the simulations identify "hot spots" in depressions and larger pits in good correlation with observed dust activity. Three-quarters of the total sublimation is produced while the sub-solar latitude is south, resulting in a distinct dichotomy in activity and morphology. Conclusions. The northern areas display a much rougher morphology than what is seen on Imhotep, an area at the equator that will be fully illuminated when 67P is closer to the Sun. Self-illumination in concave regions enhance the energy input and hence erosion. This explains the early activity observed at Hapi. Cliffs are more prone to erosion than horizontal, often dust covered, areas, which leads to surface planation. Local activity can only persist if the forming cliff walls are eroding. Comet 67P has two lobes and also two distinct sides. Transport of material from the south to the north is probable. The morphology of the Imhotep plain should be typical for the terrains of the yet unseen southern hemisphere.

Geomorphology and spectrophotometry of Philae's landing site on comet 67P/Churyumov-Gerasimenko

Context. On 12 November 2014 the European mission Rosetta succeeded in delivering a lander, named Philae, on the surface of one of the smallest, low-gravity and most primitive bodies of the solar system, the comet 67P/Churyumov-Gerasimenko (67P). Aims. The aim of this paper is to provide a comprehensive geomorphological and spectrophotometric analysis of Philae's landing site (Agilkia) to give an essential framework for the interpretation of its in situ measurements. Methods. OSIRIS images, coupled with gravitational slopes derived from the 3D shape model based on stereo-photogrammetry were used to interpret the geomorphology of the site. We adopted the Hapke model, using previously derived parameters, to photometrically correct the images in orange filter (649.2 nm). The best approximation to the Hapke model, given by the Akimov parameter-less function, was used to correct the reflectance for the effects of viewing and illumination conditions in the other filters. Spectral analyses on coregistered color cubes were used to retrieve spectrophotometric properties. Results. The landing site shows an average normal albedo of 6.7% in the orange filter with variations of similar to 15% and a global featureless spectrum with an average red spectral slope of 15.2%/100 nm between 480.7 nm (blue filter) and 882.1 nm (near-IR filter). The spatial analysis shows a well-established correlation between the geomorphological units and the photometric characteristics of the surface. In particular, smooth deposits have the highest reflectance a bluer spectrum than the outcropping material across the area. Conclusions. The featureless spectrum and the redness of the material are compatible with the results by other instruments that have suggested an organic composition. The observed small spectral variegation could be due to grain size effects. However, the combination of photometric and spectral variegation suggests that a compositional differentiation is more likely. This might be tentatively interpreted as the effect of the efficient dust-transport processes acting on 67P. High-activity regions might be the original sources for smooth fine-grained materials that then covered Agilkia as a consequence of airfall of residual material. More observations performed by OSIRIS as the comet approaches the Sun would help interpreting the processes that work at shaping the landing site and the overall nucleus.

Large-scale dust jets in the coma of 67P/Churyumov-Gerasimenko as seen by the OSIRIS instrument onboard Rosetta

Context. During the most recent perihelion passage in 2009 of comet 67P/Churyumov-Gerasimenko (67P), ground-based observations showed an anisotropic dust coma where jet-like features were detected at similar to 1.3 AU from the Sun. The current perihelion passage is exceptional as the Rosetta spacecraft is monitoring the nucleus activity since March 2014, when a clear dust coma was already surrounding the nucleus at 4.3 AU from the Sun. Subsequently, the OSIRIS camera also witnessed an outburst in activity between April 27 and 30, and since mid-July, the dust coma at rh similar to 3.7-3.6 AU preperihelion is clearly non-isotropic, pointing to the existence of dust jet-like features. Aims. We aim to ascertain on the nucleus surface the origin of the dust jet-like features detected as early as in mid-July 2014. This will help to establish how the localized comet nucleus activity compares with that seen in previous apparitions and will also help following its evolution as the comet approaches its perihelion, at which phase most of the jets were detected from ground-based observations. Determining these areas also allows locating them in regions on the nucleus with spectroscopic or geomorphological distinct characteristics. Methods. Three series of dust images of comet 67P obtained with the Wide Angle Camera (WAC) of the OSIRIS instrument onboard the Rosetta spacecraft were processed with different enhancement techniques. This was made to clearly show the existence of jet-like features in the dust coma, whose appearance toward the observer changed as a result of the rotation of the comet nucleus and of the changing observing geometry from the spacecraft. The position angles of these features in the coma together with information on the observing geometry, nucleus shape, and rotation, allowed us to determine the most likely locations on the nucleus surface where the jets originate from. Results. Geometrical tracing of jet sources indicates that the activity of the nucleus of 67P gave rise during July and August 2014 to large-scale jet-like features from the Hapi, Hathor, Anuket, and Aten regions, confirming that active regions may be present on the nucleus localized at 60. northern latitude as deduced from previous comet apparitions. There are also hints that large-scale jets observed from the ground are possibly composed, at their place of origin on the nucleus surface, of numerous small-scale features.

Morphology and dynamics of the jets of comet 67P/Churyumov-Gerasimenko: Early-phase development

Aims. The OSIRIS camera onboard the Rosetta spacecraft obtained close-up views of the dust coma of comet 67P. The jet structures can be used to trace their source regions and to examine the possible effect of gas-surface interaction. Methods. We analyzed the wide-angle images obtained in the special dust observation sequences between August and September 2014. The jet features detected in different images were compared to study their time variability. The locations of the potential source regions of some of the jets are identified by ray tracing. We used a ring-masking technique to calculate the brightness distribution of dust jets along the projected distance. Results. The jets detected between August and September 2014 mostly originated in the Hapi region. Morphological changes appeared over a timescale of several days in September. The brightness slope of the dust jets is much steeper than the background coma. This might be related to the sublimation or fragmentation of the emitted dust grains. Interaction of the expanding gas flow with the cliff walls on both sides of Hapi could lead to erosion and material down-fall to the nucleus surface.

The HARPS search for southern extra-solar planets XXXVII. Five new long-period giant planets and a system update

We describe radial-velocity time series obtained by HARPS on the 3.60 m telescope in La Silla (ESO, Chile) over ten years and report the discovery of five new giant exoplanets in distant orbits; these new planets orbit the stars HD 564, HD 30669, HD 108341, and BD -114672. Their periods range from 492 to 1684 days, semi-major axes range from 1.2 to 2.69 AU, and eccentricities range from 0 to 0.85. Their minimum mass ranges from 0.33 to 3.5 MJup. We also refine the parameters of two planets announced previously around HD 113538, based on a longer series of measurements. The planets have a period of 663 ± 8 and 1818 ± 25 days, orbital eccentricities of 0.14 ± 0.08 and 0.20 ± 0.04, and minimum masses of 0.36 ± 0.04 and 0.93 ± 0.06 MJup. Finally, we report the discovery of a new hot-Jupiter planet around an active star, HD 103720; the planet has a period of 4.5557 ± 0.0001 days and a minimum mass of 0.62 ± 0.025 MJup. We discuss the fundamental parameters of these systems and limitations due to stellar activity in quiet stars with typical 2 m s-1 radial velocity precision.

Size-frequency distribution of boulders ≥7 m on comet 67P/Churyumov-Gerasimenko

Aims. We derive for the first time the size-frequency distribution of boulders on a comet, 67P/Churyumov-Gerasimenko (67P), computed from the images taken by the Rosetta/OSIRIS imaging system. We highlight the possible physical processes that lead to these boulder size distributions. Methods. We used images acquired by the OSIRIS Narrow Angle Camera, NAC, on 5 and 6 August 2014. The scale of these images (2.44−2.03 m/px) is such that boulders ≥7 m can be identified and manually extracted from the datasets with the software ArcGIS. We derived both global and localized size-frequency distributions. The three-pixel sampling detection, coupled with the favorable shadowing of the surface (observation phase angle ranging from 48° to 53°), enables unequivocally detecting boulders scattered all over the illuminated side of 67P. Results. We identify 3546 boulders larger than 7 m on the imaged surface (36.4 km2), with a global number density of nearly 100/km2 and a cumulative size-frequency distribution represented by a power-law with index of −3.6 +0.2/−0.3. The two lobes of 67P appear to have slightly different distributions, with an index of −3.5 +0.2/−0.3 for the main lobe (body) and −4.0 +0.3/−0.2 for the small lobe (head). The steeper distribution of the small lobe might be due to a more pervasive fracturing. The difference of the distribution for the connecting region (neck) is much more significant, with an index value of −2.2 +0.2/−0.2. We propose that the boulder field located in the neck area is the result of blocks falling from the contiguous Hathor cliff. The lower slope of the size-frequency distribution we see today in the neck area might be due to the concurrent processes acting on the smallest boulders, such as i) disintegration or fragmentation and vanishing through sublimation; ii) uplifting by gas drag and consequent redistribution; and iii) burial beneath a debris blanket. We also derived the cumulative size-frequency distribution per km2 of localized areas on 67P. By comparing the cumulative size-frequency distributions of similar geomorphological settings, we derived similar power-law index values. This suggests that despite the selected locations on different and often opposite sides of the comet, similar sublimation or activity processes, pit formation or collapses, as well as thermal stresses or fracturing events occurred on multiple areas of the comet, shaping its surface into the appearance we see today.

Theoretical models of planetary system formation II. Post-formation evolution

Aims. We extend the results of planetary formation synthesis by computing the long-term evolution of synthetic systems from the clearing of the gas disk into the dynamical evolution phase. Methods. We use the symplectic integrator SyMBA to numerically integrate the orbits of planets for 100 Myr, using populations from previous studies as initial conditions. Results. We show that within the populations studied, mass and semimajor axis distributions experience only minor changes from post-formation evolution. We also show that, depending upon their initial distribution, planetary eccentricities can statistically increase or decrease as a result of gravitational interactions. We find that planetary masses and orbital spacings provided by planet formation models do not result in eccentricity distributions comparable to observed exoplanet eccentricities, requiring other phenomena, such as stellar fly-bys, to account for observed eccentricities.

Comet 67P/Churyumov-Gerasimenko: Constraints on its origin from OSIRIS observations

Context. One of the main aims of the ESA Rosetta mission is to study the origin of the solar system by exploring comet 67P/Churyumov-Gerasimenko at close range. Aims. In this paper we discuss the origin and evolution of comet 67P/Churyumov-Gerasimenko in relation to that of comets in general and in the framework of current solar system formation models. Methods. We use data from the OSIRIS scientific cameras as basic constraints. In particular, we discuss the overall bi-lobate shape and the presence of key geological features, such as layers and fractures. We also treat the problem of collisional evolution of comet nuclei by a particle-in-a-box calculation for an estimate of the probability of survival for 67P/Churyumov-Gerasimenko during the early epochs of the solar system. Results. We argue that the two lobes of the 67P/Churyumov-Gerasimenko nucleus are derived from two distinct objects that have formed a contact binary via a gentle merger. The lobes are separate bodies, though sufficiently similar to have formed in the same environment. An estimate of the collisional rate in the primordial, trans-planetary disk shows that most comets of similar size to 67P/Churyumov-Gerasimenko are likely collisional fragments, although survival of primordial planetesimals cannot be excluded. Conclusions. A collisional origin of the contact binary is suggested, and the low bulk density of the aggregate and abundance of volatile species show that a very gentle merger must have occurred. We thus consider two main scenarios: the primordial accretion of planetesimals, and the re-accretion of fragments after an energetic impact onto a larger parent body. We point to the primordial signatures exhibited by 67P/Churyumov-Gerasimenko and other comet nuclei as critical tests of the collisional evolution.