This article investigates reconstructing the internal mass density of a numerical asteroid model using the gradient of a simulated gravity field as synthetic measurement data. Our goal is to advance the mathematical inversion methodology and find feasibility constraints for the resolution, noise and orbit selection for future space missions. We base our model on the shape of the asteroid Itokawa as well as on the recent observations and simulation studies which suggest that the internal density varies, increasing towards the center, and that the asteroid may have a detailed structure. We introduce randomized multiresolution scan algorithm which might provide a robust way to cancel out bias and artifact effects related to the measurement noise and numerical discretization. In this scheme, the inverse algorithm can reconstruct details of various sizes without fixing the exact resolution {\em a priori}, and the randomization minimizes the effect of discretization on the solution. We show that the adopted methodology provides an advantageous way to diminish the surface bias of the inverse solution. The results also suggest that a noise level below 80 Eotvos will be sufficient for the detection of internal voids and high density anomalies, if a sparse set of measurements can be obtained from a close-enough distance to the target.

In this paper, we investigate full-wave computed radar tomography (CRT) using a rubble-pile asteroid model in which a realistic shape (Itokawa) is coupled with a synthetic material composition and structure model. The aim is to show that sparse bistatic radar measurements can distinguish details inside a complex-structured rubble-pile asteroid. The results obtained suggest that distinct local permittivity distribution changes such as surface layers, voids, low-permittivity anomalies, high-permittivity boulders, and cracks can be detected with bistatic CRT, when the total noise level in the data is around -10 dB with respect to the signal amplitude. Moreover, the bistatic measurement set-up improves the robustness of the inversion compared to the monostatic case. Reconstructing the smooth Gaussian background distribution was found to be difficult with the present approach, suggesting that complementary techniques, such as gravimetry, might be needed to improve the reliability of the inference in practice.

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.

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.

We have performed an initial stage conceptual design study for the Deep Interior Scanning CubeSat (DISCUS), a tandem 6U CubeSat carrying a bistatic radar as main payload. DISCUS will be operated either as an independent mission or accompanying a larger one. It is designed to determine the internal macroporosity of a 260-600 m diameter Near Earth Asteroid (NEA) from a few kilometers distance. The main goal will be to achieve a global penetration with a low-frequency signal as well as to analyze the scattering strength for various different penetration depths and measurement positions. Moreover, the measurements will be inverted through a computed radar tomography (CRT) approach. The scientific data provided by DISCUS would bring more knowledge of the internal configuration of rubble pile asteroids and their collisional evolution in the Solar System. It would also advance the design of future asteroid deflection concepts. We aim at a single-unit (1U) radar design equipped with a half-wavelength dipole antenna. The radar will utilize a stepped-frequency modulation technique the baseline of which was developed for ESA's technology projects GINGER and PIRA. The radar measurements will be used for CRT and shape reconstruction. The CubeSat will also be equipped with an optical camera system and laser altimeter to sup- port navigation and shape reconstruction. We provide the details of the measurement methods to be applied along with the requirements derived of the known characteristics of rubble pile asteroids.

This study aims at advancing mathematical and computational techniques for reconstructing the interior structure of a near earth asteroid (NEA) via Computed Radar Tomography (CRT). We introduce a far-field model for full-wave CRT and validate it numerically for sparse and limited-angle data and an orbiting distance of about 5 km. As a potential future application of the proposed method, we consider the Deep Interior Scanning CUbeSat (DISCUS) concept in which two small spacecraft carrying a bistatic radar localize internal macropores inside a rubble pile asteroid.

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.

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.

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.

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.

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.

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.

P/2013 P5 PANSTARRS was discovered in Aug. 2013, displaying a cometary tail, but with orbital elements typical for a member of the inner asteroid Main Belt. We monitored the object from 2013 Aug. 30 until Oct. 05 using the CFHT, NTT, CA 1.23m, Perkins 1.8m (Lowell), and the 0.6m TRAPPIST telescopes. We measured its nuclear radius to be r < 0.25-0.29km, and its colours g-r = 0.58+/-0.05 and r-i = 0.23+/-0.06, typical for an S-class asteroid. We failed to detect any rotational light curve, with an amplitude < 0.05mag and a double-peaked rotation period < 20h. A detailed Finson-Probstein analysis of deep NTT and CFHT images indicated that the object was active since at least late January 2013 until the time of the latest observations in 2013 September, with at least two peaks of activity around 2013 June 14+/-10d and 2013 July 22+/-3d. The changes of activity level and the activity peaks were extremely sharp and short, shorter than the temporal resolution of our observations (about 1d). The dust distribution was similar during these two events, with dust grains covering at least the 1-1000{\mu}m range. The total mass ejected in grains <1mm was estimated to be 3.0 10$^6$kg and 2.6 10$^7$kg around the two activity peaks. Rotational disruption cannot be ruled out as the cause of the dust ejection. We also propose that the components of a contact binary might gently rub and produce the observed emission. Volatile sublimation might also explain what appears as cometary activity over a period of 8 months. However, while Main Belt comets best explained by ice sublimation are found in the outskirts of the Main Belt, where water ice is believed to be able to survive buried in moderately large objects for the age of the solar system deeply, the presence of volatiles in an object smaller than 300m in radius would be very surprising in the inner asteroid belt.

We report new mid-eclipse times of the short-period sdB/dM binary HW Vir, which differ substantially from the times predicted by a previous model. The proposed orbits of the two planets in that model are found to be unstable. We present a new secularly stable solution, which involves two companions orbiting HW VIr with periods of 12.7 yr and 55 +/-15 yr. For orbits coplanar with the binary, the inner companion is a giant planet with mass M_3 sin i_3 = 14 M_Jup and the outer one a brown dwarf or low-mass star with a mass of M_4 sin i_4 = 30-120 M_Jup. Using the mercury6 code, we find that such a system would be stable over more than 10^7 yr, in spite of the sizeable interaction. Our model fits the observed eclipse-time variations by the light-travel time effect alone, without invoking any additional process, thereby providing support for the planetary hypothesis of the eclipse-time variations in close binaries. The signature of non-Keplerian orbits may be visible in the data.