Simulation of Europa's water plume and structures related to energetic activities on solar system bodies from satellite images

This PhD thesis focuses on the analysis of different high energetic processes that affect the surface of planets, satellites and minor bodies as well as modify their surrounding environment. Specifically, this work concerns three main topics: (i) the simulation and analysis of one of the most geological energetic process, i.e. impact cratering; (ii) the investigation of the fragmentation processes that could have generated boulders on comet Churyumov-Gerasimenko 67P; (iii) the analysis of a possible transient plume originating from cryovolcanic events on Europa, the Jovian icy satellite, combined with an accurate characterisation of its exospheric background. The first topic addresses the investigation of the impact formation process through numerical modelling. Shocks code represent the most feasible method for studying impact craters, as they can simulate a large span of conditions beyond the reach of experiments (e.g., velocity, size). The iSALE hydrocode was used to simulate two different impact structures located on Mercury and Mars. On Mercury, the simulation allows to determine the genesis of a particular landform, i.e. a steep-sided cone with associated pyroclastic deposits, which was revealed by images acquired by MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft. On Mars, the simulation of the Firsoff crater in Arabia Terra permits a better understanding the subsequent geologic processes that led to crater post modification, defining which rheological structure is more likely in that region. In both cases, we conclude that the numerical modelling of impact process is a powerful tool to improve the comprehension of the Solar System. The second topic of the thesis has been developed after the Rosetta mission got inserted around the comet Churyumov-Gerasimenko 67P. We investigated the surface of comet Churyumov-Gerasimenko 67P focusing on the possible energetic events that lead to the formation of boulders; i.e. blocks that are ubiquitous on the surface of the comet. Different energetic formation processes were invoked to explain the presence of boulders, such as sublimation, fragmentation, outbursts and gravitational falls. Using images acquired by OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) and CIVA (Comet Infrared and Visible Analyser) cameras on board the spacecraft and the lander respectively, a quantitative analysis of different-size boulders has been performed in order to understand if the same energetic formation processes occur equally on different scales on the comet (m, cm and mm). Specifically, by means of different resolution images, we obtained several size-frequency distribution for: (i) boulders larger than 7 m, (ii) boulders larger than 1 m from higher resolution images used to analyse the Abydos site, the location where Philae is supposed to be, and (iii) pebbles (mm-scale structures) visible on CIVA images. The third topic is the icy satellite Europa in view of the future ESA/JUICE mission and because of our involvement in JANUS (Jovis, Amorum ac Natorum Undique Scrutator) visible camera. The presence of a subsurface ocean is a primary topic on Europa, in addition the recent discovery of a transient plume at the south pole by HST observations has raised many questions regarding the interaction between the subsurface/surface and the outer environment of Europa in terms of active processes affecting the icy surface. In this context, a possible plume deposit originating from cryovolcanic events was simulated to understand its detectability by JANUS camera during the Europa flyby phase of the JUICE mission. In addition, since the study of transient plumes has as a mandatory prerequisite an accurate characterisation of the exospheric background, a detailed study of the loss rates of Europa's tenuous atmosphere was performed. In particular, loss rates for electron impact dissociation and ionization processes, for charge-exchange (considering plasma torus, pick up and ionosphere ions) and for photo processes (for both cases of quite and active Sun) were calculated

Origin(s) of the local structures at the Philae landing site and possible implications on the formation and evolution of the 67P nucleus

The in situ images of the 67P/Churyumov-Gerasimenko nucleus acquired by the CIVA cameras on-board PHILAE revealed a rough, irregular and inhomogeneous terrains dominated by fractures and agglomerates of consolidated materials. While the composition of these materials is unknown, they provide unique structures to constrain the conditions prevailing at the surface of a comet and also possibly to the nucleus formation. A quantitative analysis of some microscopic structures (namely fractures and grains that look like pebbles) will be presented using a manual extraction from the CIVA data set with the software ArcGIS. Fractures/cracks are rather ubiquitous at various spatial scales with network and size (from sub-cm to 10 cm) well correlated to the texture of the landscape. The pebble size distribution are reasonably well fitted by power-laws having different cumulative indexes. The nature of the landscape of the landing site will be then discussed in relation to both endogenic and exogenic processes that could have sculpted it. The block seen in CIVA#1 is interpreted to be a close-up of fractured boulders/cliff belonging to the boulder field identified from the orbit near Abydos, this boulder field being itself the result of gravitational regressive erosion due to sublimation (Lucchetti et al. 2016). The observed fractures are best explained by thermal insolation leading to thermal fatigue and/or to loss of volatile materials (e.g., desiccation). This surficial fragmentation (up to >10 cm length) could generate macroscopic erosion that is also visible at larger scale from the orbit. While the pebbles are difficult to be formed by any current physical processes, there is at least an intriguing possibility that they are remnants of primordial accretion processes, as there are several lines of evidence that the nucleus could be primordial (Davidsson et al., 2016), and not a collisional rubble piles of a large body (Morbidelli and Rickman A&A, 2015). We thus speculate that the Abydos landscape could be in favour of pebble accretion model instead of runaway coagulation model with a formation location in the outer region of the Solar System

New simulation of Phobos Stickney crater

In this work we model the Phobos Stickney impact crater using the iSALE hydrocode and considering different scenarios that could form the well-studied crater

The heterogeneous ice shell thickness of Enceladus

Saturn's moon Enceladus is the smallest Solar System body that presents an intense geologic activity on its surface. Plumes erupting from Enceladus' South Polar terrain (SPT) provide direct evidence of a reservoir of liquid below the surface. Previous analysis of gravity data determined that the ice shell above the liquid ocean must be 30-40 km thick from the South Pole up to 50° S latitude (Iess et al., 2014), however, understand the global or regional nature of the ocean beneath the ice crust is still challenging. To infer the thickness of the outer ice shell and prove the global extent of the ocean, we used the self-similar clustering method (Bonnet et al., 2001; Bour et al., 2002) to analyze the widespread fractures of the Enceladus's surface. The spatial distribution of fractures has been analyzed in terms of their self-similar clustering and a two-point correlation method was used to measure the fractal dimension of the fractures population (Mazzarini, 2004, 2010). A self-similar clustering of fractures is characterized by a correlation coefficient with a size range defined by a lower and upper cut-off, that represent a mechanical discontinuity and the thickness of the fractured icy crust, thus connected to the liquid reservoir. Hence, this method allowed us to estimate the icy shell thickness values in different regions of Enceladus from SPT up to northern regions.We mapped fractures in ESRI ArcGis environment in different regions of the satellite improving the recently published geological map (Crow-Willard and Pappalardo, 2015). On these regions we have taken into account the fractures, such as wide troughs and narrow troughs, located in well-defined geological units. Firstly, we analyzed the distribution of South Polar Region fracture patterns finding an ice shell thickness of ~ 31 km, in agreement with gravity measurements (Iess et al., 2014). Then, we applied the same approach to other four regions of the satellite inferring an increasing of the ice shell thickness from 31 to 70 km from the South Pole to northern regions. By these findings, we prove the global extent of the ocean underneath the ice crust of the satellite

Fractal analysis on Enceladus: a global ocean underneath the icy crust

Plumes of water have been observed erupting from Enceladus' south polar terrain providing direct evidence of a reservoir of liquid below the surface, that could be considered global or just a small body of water concentrated at its south pole. Gravity data collected during the spacecraft's several close flyby over the south polar region determined that the icy shell above the liquid ocean must be 30-40 km thick extending from the south pole up to 50°S (Iess et al. 2014). The hypothesis of a global ocean beneath the icy crust has been raised even in a recent paper of Thomas et al. (2015) thanks to the measurements of the very slight wobble that Enceladus displays as it orbits Saturn. In this work we support the hypothesis of the presence of an ocean layer using the fractal percolation theory. This method allowed us to estimate the icy shell thickness values in different regions of Enceladus from the south polar terrain up to the north pole. The spatial distribution of fractures on Enceladus has been analyzed in terms of their self-similar clustering and a two-point correlation method was used to measure the fractal dimension of the fractures population (Mazzarini, 2004, 2010). A self-similar clustering of fractures is characterized by a correlation coefficient with a size range defined by a lower and upper cut-off, that represent a mechanical discontinuity and the thickness of the fractured icy crust, thus connected to the liquid reservoir. We mapped the fractures on Enceladus surface based on April 2010 global mosaic from Cassini mission and applied the fractal method firstly to the south polar terrain finding indeed a fractal correlation of fractures and providing an ice shell thickness of ~40 km. Then, we analyzed fractures of four different regions around the equator and around the north pole inferring an overall ice shell thickness ranging from 35 to 45 km. Our results are in agreement with the gravity observations (Iess et al., 2014) and the mechanical models (Yin and Pappalardo, 2015) of the south polar region and, in addition, they are consistent with the hypothesis of a global ocean beneath the icy shell

Loss rates of Europa's tenuous atmosphere

Loss processes in Europa's tenuous atmosphere are dominated by plasma-neutral interactions. Based on the updated data of the plasma conditions in the vicinity of Europa (Bagenal et al. 2015), we provide estimations of the atmosphere loss rates for the H2O, O2 and H2 populations. Due to the high variability of the plasma proprieties, we perform our investigation for three sample plasma environment cases identified by Bagenal et al. as hot/low density, cold/high density, and an intermediate case. The role of charge-exchange interactions between atmospheric neutrals and three different plasma populations, i.e. magnetospheric, pickup, and ionospheric ions, is examined in detail. Our assumptions related to the pickup and to the ionospheric populations are based on the model by Sittler et al. (2013). We find that O2-O2+ charge-exchange is the fastest loss process for the most abundant atmospheric species O2, though this loss process has been neglected in previous atmospheric models. Using both the revised O2 column density obtained from the EGEON model (Plainaki et al., 2013) and the current loss rate estimates, we find that the upper limit for the volume integrated loss rate due to O2-O2+ charge exchange is in the range (13-51)×1026 s-1, depending on the moon's orbital phase and illumination conditions. The results of the current study are relevant to the investigation of Europa's interaction with Jupiter's magnetospheric plasma

Planetary space weather: scientific aspects and future perspectives

International audienceIn this paper, we review the scientific aspects of planetary space weather at different regions of our Solar System, performing a comparative planetology analysis that includes a direct reference to the circum-terrestrial case. Through an interdisciplinary analysis of existing results based both on observational data and theoretical models, we review the nature of the interactions between the environment of a Solar System body other than the Earth and the impinging plasma/radiation, and we offer some considerations related to the planning of future space observations. We highlight the importance of such comparative studies for data interpretations in the context of future space missions (e.g. ESA JUICE; ESA/JAXA BEPI COLOMBO). Moreover, we discuss how the study of planetary space weather can provide feedback for better understanding the traditional circum-terrestrial space weather. Finally, a strategy for future global investigations related to this thematic is proposed

Simulation of Europa's water plume

Plumes on Europa would be extremely interesting science and mission targets, particularly due to the unique opportunity to obtain direct information on the subsurface composition, thereby addressing Europa's potential habitability. The existence of water plume on the Jupiter's moon Europa has been long speculated until the recent discover. HST imaged surpluses of hydrogen Lyman alpha and oxygen emissions above the southern hemisphere in December 2012 that are consistent with two 200 km high plumes of water vapor (Roth et al. 2013). In previous works ballistic cryovolcanism has been considered and modeled as a possible mechanism for the formation of low-albedo features on Europa's surface (Fagents et al. 2000). Our simulation agrees with the model of Fagents et al. (2000) and consists of icy particles that follow ballistic trajectories. The goal of such an analysis is to define the height, the distribution and the extension of the icy particles falling on the moon's surface as well as the thickness of the deposited layer. We expect to observe high albedo regions in contrast with the background albedo of Europa surface since we consider that material falling after a cryovolcanic plume consists of snow. In order to understand if this phenomenon is detectable we convert the particles deposit in a pixel image of albedo data. We consider also the limb view of the plume because, even if this detection requires optimal viewing geometry, it is easier detectable in principle against sky. Furthermore, we are studying the loss rates due to impact electron dissociation and ionization to understand how these reactions decrease the intensity of the phenomenon. We expect to obtain constraints on imaging requirements necessary to detect potential plumes that could be useful for ESA's JUICE mission, and in particular for the JANUS camera (Palumbo et al. 2014)