Les comètes présentent un grand intérêt à la fois pour la planétologie et pour l'exobiologie. En effet, ces corps primitifs du fait de leur petite taille et de leurs réservoirs éloignés du soleil, n'ont pas ou que très peu évolué depuis leur formation. L'étude des comètes peut donc permettre de mieux comprendre les processus physico-chimiques ayant eu lieu lors de la formation du Système Solaire. D'autre part, les analyses menées en 1986 dans l'environnement de la comète 1P/Halley ont montré l'existence, dans les grains cométaires d'une phase solide riche en composés organiques. Ainsi, les comètes ont vraisemblablement pu apporter sur la Terre primitive des composés organiques, et favoriser ainsi l'apparition de la Vie. Néanmoins la nature de cette matière organique reste encore très largement méconnue. Ces composés organiques ont vraisemblablement été formés à partir des glaces observées dans le milieu interstellaire et qui sont soumises à différentes sources d'énergie. Les objectifs du travail expérimental mené au cours de cette thèse ont donc été de caractériser les différentes étapes conduisant à la synthèse des composés organiques complexes contenus dans les comètes à partir des glaces interstellaires. En particulier, j'ai étudié i.) la quantification des glaces présentes autour des étoiles naissantes, ii.) les processus de photolyse auxquels elles sont soumises et iii.) la nature des composés organiques qui peuvent être produits durant ces processus. Les observations infrarouges ont permis de détecter de nombreuses molécules en phase condensée autour des étoiles naissantes. Afin de préciser l'abondance de ces molécules, j'ai mesuré les sections efficaces intégrées, aussi appelées forces de bandes, pour huit d'entre elles (H2O, CO, CO2, CH3OH, NH3, CH4, HCOOH and H2CO). En effet, ce paramètre spectroscopique est nécessaire à la quantification des molécules et certaines des valeurs présentes dans la littérature affichaient une grande dispersion. Les nouvelles mesures que j'ai effectuées, basées sur une revue bibliographique exhaustive des masses volumiques et des indices optiques dans le visible, confirment pour certaines molécules (CO2, CO, CH4, NH3) les valeurs utilisées pour déterminer leur abondance. Néanmoins, j'ai pu montrer que les abondances d'autres molécules dans les milieux astrophysiques restent encore très incertaines, en particulier pour CH3OH, H2CO et HCOOH. En phase condensée, la dépendance en longueur d'onde des processus de photolyse est encore très largement négligée. Grâce à l'utilisation de deux lampes VUV, dont l'irradiance spectrale a été préalablement caractérisée, j'ai pu mesurer, dans deux gammes de longueurs d'onde différentes, les rendements quantiques de production de C2H6 et de CO lors de la photolyse de CH4 et CO2. J'ai ainsi pu confirmer que les rendements quantiques dépendent bien de la longueur d'onde de photolyse. L'extrapolation des résultats expérimentaux acquis en laboratoire aux différents milieux astrophysiques nécessite donc une bonne connaissance des spectres VUV mis en jeu. L'objectif final des simulations expérimentales est de prédire la nature de la matière organique cométaire en reproduisant au mieux la chimie pouvant se dérouler dans les glaces interstellaires. Or, le méthane a été détecté en phase condensée dans le milieu interstellaire, mais son influence sur la chimie se déroulant dans des mélanges de glaces contenant les principales molécules interstellaires a été très peu étudiée. J'ai donc soumis un mélange H2O : CH3OH : NH3 : CH4 (10 : 1 : 1 : 2) à une photolyse de 26 heures puis à un chauffage. L'influence du méthane se manifeste par la présence de C2H6 après la photolyse à basse température. Lors du chauffage, le méthane et ses photoproduits semblent se sublimer. Avec ou sans CH4, la chimie à des températures supérieures à 200 K apparaît très similaire. J'en conclu donc que la présence de méthane ne modifie pas notablement la chimie des glaces lors des simulations / Comets are very interesting for planetology as well as for exobiology. On one hand, held in the furthest and coldest regions of our solar system and due to their small size, they might not have been altered since their formation. The study of comets should allow a better understanding of the physic-chemical processes occurring during the Solar system formation. On the other hand, the analysis performed in 1986 on the environment of 1P/Halley showed the presence, in the cometary dust, of organic matter. Thus, comets might have brought organics on primitive Earth which might have contributed to the apparition of life. Nevertheless, the nature of these organics is still not well-known. Cometary organics might have been synthesized from the ices detected in interstellar medium which are submitted to different energetic processes. The aims of the experimental work performed during this thesis are to characterize the different steps of the synthesis of complex organic matter contained in comets from the interstellar ices. I studied : i) The quantification of interstellar ices detected around young stellar objects ii) The characterization of the photolysis process to which ices are submitted and iii.) The nature of the organic compounds produced during these processes. Once mixtures and energetic processes are under control, we can make cometary organic analogs. Infrared observations have revealed the presence of several molecules in the solid phase around young stellar objects. To precise their molecular abundances, I have measured the integrated cross sections, also called band strengths, of 8 molecules (H2O, CO, CO2, CH3OH, NH3, CH4, HCOOH and H2CO). Indeed, this spectroscopic parameter is required for the quantification of these molecules and some values presented in literature are scattered. The new measurements performed during this thesis, which are based on a bibliographic review of densities and optical indices in the visible range, confirm the values already used for the quantification of CO2, CO, CH4, NH3. But this work also underlines that abundances of CH3OH, H2CO and HCOOH in interstellar medium are still uncertain. In the solid phase, wavelength dependence of photolysis is often neglected. Thanks to two VUV lamps, for which the spectral irradiances have been characterized, I measured the production quantum yield, in two wavelength ranges, of C2H6 and CO, during photolysis of CH4 and CO2 respectively. Thanks to this study, I point out that quantum yield depends on the photolysis wavelength. Thus, the extrapolation of the experimental results to different astrophysical medium implies a good knowledge of VUV spectra. The final objective of experimental simulations is to foresee the nature of cometary organic matter by reproducing, as realistic as possible, the chemistry occurring in interstellar ices. Methane has been detected in the solid phase in the interstellar medium, but few studies implying methane have been undertaken. Thus, I have photolyzed a mixture composed of H2O: CH3OH: NH3: CH4 (10:1:1:2) during 26 hours at low temperature and then I applied a heating process. The influence of initial methane in the ice chemistry is demonstrated by the presence of its main photoproducts, C2H6, after photolysis. But while increasing temperature, methane and its photoproducts seem to sublimate. Therefore, with or without methane, chemistry occurring at temperature higher than 200K seems to be very similar. I conclude than methane does not have significant influence on ice chemistry

The thesis contains a survey of numerical tools for studying thermomechanical interactions of a two-phase system contained in a domain with an upper bound- ary that forms a free surface. The enthalpy diffused-interface formulation is used for an approximation of the phase change interface and the computing algorithm is benchmarked against an analytical solution of the Stefan problem. Arbitrary Lagrangian-Eulerian kinematical description of the continuum is applied to over- come the difficulty in the form of the free surface. The validity of the approach is examined on a thermal convection benchmark problem. 1

Chapter 1-2:The Cassini RADAR mapper has imaged elevated blocks and mountains on Titan we term ‘ridges’. Two unresolved problems regarding Titan's surface are still debated: what is the origin of its ridges and was there tectonic activity on Titan? To understand the processes that produced the ridges, in this study, (1) we analyze the distribution and orientation of ridges through systematic geomorphologic mapping and (2) we compare the location of the ridges to a new global topographic map to explore the correlation between elevation and ridges and the implications for Titan's surface evolution. Globally, the orientation of ridges is nearly E-W and the ridges are more common near the equator than at the poles, which suggests a tectonic origin for most of the ridges on Titan. In addition, the ridges are found to preferentially lie at higher-than-average elevations near the equator. We conclude the most reasonable formation scenario for Titan's ridges is that contractional tectonism built the ridges and thickened the icy lithosphere, causing regional uplift. The combination of global and regional tectonic events, likely contractional in nature, plus enhanced fluvial erosion and sedimentation near the poles, would have contributed to shaping Titan's tectonic landforms and surface morphology to what we see today. However, contractional structures (i.e. thrusts and folds) require large stresses (8~10 MPa), the sources of which probably do not exist on Titan. Liquid hydrocarbons in Titan's near subsurface must play a role similar to that of water on Earth and lead to fluid overpressures, which enable contractional deformation at smaller stresses (< 1MPa) by significantly reducing the shear strength of materials. We show that crustal conditions with enhanced pore fluid pressures on Titan favor the formation of thrust faults and related folds, in a contractional stress field. The production of folds, as on Earth, is facilitated by the presence of crustal liquids to weaken the crust. These hydrocarbon fluids have played a key role in Titan's tectonic evolutionary history, leaving it the only icy body on which strong evidence for contractional tectonism exists. Chapter 3: Arthur Wichmann's ‘Earthquakes of the Indian Archipelago’ documents several large earthquakes and tsunami throughout the Banda Arc region that can be interpreted as mega-thrust events. However, the source regions of these events are not known. One of the largest and well-documented events in the catalog is the great earthquake and tsunami affecting the Banda islands on 1 August 1629. It caused severe damage from a 15-meter tsunami that arrived at the Banda Islands about a half hour after violent shaking stopped. The earthquake was also recorded 230 km away in Ambon, but no tsunami is mentioned. This event was followed by at least 9 years of uncommonly frequent seismic activity in the region that tapered off with time, which can be interpreted as aftershocks. The combination of these observations indicates that the earthquake was most likely a mega-thrust event. We use an inverse modeling approach to numerically reconstruct the tsunami, which constrains the likely location and magnitude of the 1629 earthquake. Only linear numerical models are applied due to the low-resolution of bathymetry in the Banda Islands and Ambon. Therefore, we apply various wave amplification factors (1.5 to 4) derived from simulations of recent, well-constrained tsunami to bracket the upper and lower limits of earthquake moment magnitudes for the event. The closest major earthquake sources to the Banda Islands are the Tanimbar and Seram Troughs of the Banda subduction/collision zone. Other source regions are too far away for such a short arrival time of the tsunami after shaking. Moment magnitudes predicted by the models in order to produce a 15 m tsunami are Mw of 9.8 to 9.2 on the Tanimbar Trough and Mw 8.8 to 8.2 on the Seram Trough. The arrival times of these waves are 58 minutes for Tanimbar Trough and 30 minutes for Seram Trough. The model also predicts 5 meters run-up for Ambon from a Tanimbar Trough source, which is inconsistent with the historical records. Ambon is mostly shielded from a wave generated by a Seram Trough Source.We conclude that the most likely source of the 1629 mega-thrust earthquake is the Seram Trough. Only one earthquake > Mw 8.0 is recorded instrumentally from the eastern Indonesia region although high rates of strain (50-80 mm/a) are measured across the Seram section of the Banda subduction zone. Enough strain has already accumulated since the last major historical event to produce an earthquake of similar size to the 1629 event. Due to the rapid population growth in coastal areas in this region, it is imperative that the most vulnerable coastal areas prepare accordingly.

Here I study the “plasma interactions with icy bodies in the solar system”, that is, my quest to understand the fundamental processes that govern such interactions. By using numerical modelling combined with in situ observations, one can infer the internal structure of icy bodies and their plasma environments. After a broad overview of the laws governing space plasmas a more detailed part follows. This contains the method on how to model the interaction between space plasmas and icy bodies. Numerical modelling of space plasmas is applied to the icy bodies Callisto (a satellite of Jupiter), the dwarf planet Ceres (located in the asteroid main belt) and the comet 67P/Churyumov-Gerasimenko. The time-varying magnetic field of Jupiter induces currents inside the electrically conducting moon Callisto. These create magnetic field perturbations thought to be related to conducting subsurface oceans. The flow of plasma in the vicinity of Callisto is greatly affected by these magnetic field perturbations. By using a hybrid plasma solver, the interaction has been modelled when including magnetic induction and agrees well with magnetometer data from flybys (C3 and C9) made by the Galileo spacecraft. The magnetic field configuration allows an inflow of ions onto Callisto’s surface in the central wake. Plasma that hits the surface knocks away matter (sputtering) and creates Callisto’s tenuous atmosphere. A long term study of solar wind protons as seen by the Rosetta spacecraft was conducted as the comet 67P/Churyumov-Gerasimenko approached the Sun. Here, extreme ultraviolet radiation from the Sun ionizes the neutral water of the comet’s coma. Newly produced water ions get picked up by the solar wind flow, and forces the solar wind protons to deflect due to conservation of momentum. This effect of mass-loading increases steadily as the comet draws closer to the Sun. The solar wind is deflected, but does not lose much energy. Hybrid modelling of the solar wind interaction with the coma agrees with the observations; the force acting to deflect the bulk of the solar wind plasma is greater than the force acting to slow it down. Ceres can have high outgassing of water vapour, according to observations by the Herschel Space Observatory in 2012 and 2013. There, two regions were identified as sources of water vapour. As Ceres rotates, so will the source regions. The plasma interaction close to Ceres depends greatly on the source location of water vapour, whereas far from Ceres it does not. On a global scale, Ceres has a comet-like interaction with the solar wind, where the solar wind is perturbed far downstream of Ceres. / Här studerar jag “plasmaväxelverkan med isiga kroppar i solsystemet”, det vill säga, min strävan är att förstå de grundläggande processerna som styr sådana interaktioner. Genom att använda numerisk modellering i kombination med observationer på plats vid himlakropparna kan man förstå sig på deras interna strukturer och rymdmiljöer. Efter en bred översikt över de fysiska lagar som styr ett rymdplasma följer en mer detaljerad del. Denna innehåller metoder för hur man kan modellera växelverkan mellan rymdplasma och isiga kroppar. Numerisk modellering av rymdplasma appliceras på de isiga himlakropparna Callisto (en måne kring Jupiter), dvärgplaneten Ceres (lokaliserad i asteroidbältet mellan Mars och Jupiter) och kometen 67P/Churyumov-Gerasimenko. Det tidsvarierande magnetiska fältet kring Jupiter inducerar strömmar inuti den elektriskt ledande månen Callisto. Dessa strömmar skapar magnetfältsstörningar som tros vara relaterade till ett elektriskt ledande hav under Callistos yta. Plasmaflödet i närheten av Callisto påverkas i hög grad av dessa magnetfältsstörningar. Genom att använda en hybrid-plasma-lösare har växelverkan modellerats, där effekten av magnetisk induktion har inkluderats. Resultaten stämmer väl överens med magnetfältsdata från förbiflygningarna av Callisto (C3 och C9) som gjordes av den obemannade rymdfarkosten Galileo i dess bana kring Jupiter. Den magnetiska konfigurationen som uppstår möjliggör ett inflöde av laddade joner på Callistos baksida. Plasma som träffar ytan slår bort materia och skapar Callistos tunna atmosfär. En långtidsstudie av solvindsprotoner sett från rymdfarkosten Rosetta utfördes då kometen 67P/Churyumov-Gerasimenko närmade sig solen. Ultraviolett strålning från solen joniserar det neutrala vattnet i kometens koma (kometens atmosfär). Nyligt joniserade vattenmolekyler plockas upp av solvindsflödet och tvingar solvindsprotonernas banor att böjas av, så att rörelsemängden bevaras. Denna effekt ökar stadigt då kometen närmar sig solen. Solvinden böjs av kraftigt, men förlorar inte mycket energi. Hybridmodellering av solvindens växelverkan bekräftar att kraften som verkar på solvinden till störst del får den att böjas av, medan kraften som verkar till att sänka dess fart är mycket lägre. Ceres har enligt observationer av rymdteleskopet Herschel under 2012 och 2013 haft högt utflöde av vattenånga från dess yta. Där har två regioner identifierats som källor för vattenångan. Eftersom Ceres roterar kommer källornas regioner göra det också. Plasmaväxelverkan i närheten av Ceres beror i hög grad på vattenångskällans placeringen, medan det inte gör det långt ifrån Ceres. På global nivå har Ceres en kometliknande växelverkan med solvinden, där störningar i solvinden propagerar långt nedströms från Ceres.

La glace de H2O est présente dans de nombreux environnements planétaires, et notamment sous forme de polymorphe de haute pression au sein des satellites de glaces ainsi que dans le manteau des planètes extrasolaires, dites planètes océan. La diversité des conditions thermodynamiques prédite au sein de ces corps planétaires a souligné le besoin de nouvelles données de laboratoire et de calculs sur les glaces de H2O afin de pouvoir modéliser leur évolution et leur structure interne.Si les propriétés structurales et spectroscopiques des pôles purs de ces glaces sont déjà relativement bien connues, une description pétrologique plus réaliste des solutions solides et des phases riches en impureté, manque encore à la communauté. Ce travail de thèse s’est concentré sur l’étude de la fusion des glaces VI et VII dans le binaire H2O-NaCl grâce aux techniques de cellules à enclumes en diamants et la spectroscopie vibrationelle Raman. Ces données ont été complétées par des mesures du fractionnement du sel analogue RbI entre les glace VI et VII et le fluide aqueux en utilisant la cartographie de fluorescence X et de diffraction des rayons X réalisées à l’European Synchrotron Research Facility (Grenoble). Ceci as permis de mettre en évidence une inversion de densité entre le fluide riche en sel et la glace VI et de révéler une forte différence de partage du sel entre la glace VI et la glace VII avec un coefficient de partage du RbI estimé à Kd(VI-VII)=4.5(±2.7)10-2.Au sein des plus gros corps riches en H2O appelés planète océan, le manteau de glace potentiellement épais de plus de 1000 km abrite un type de glace de ultra haute pression appelé glace X. Cette phase de la glace d’eau est unique de part sa structure cristallographique ionique, contrairement aux autres glaces de plus basse pression, toutes de structure moléculaire. Cette caractéristique structurale et l’absence de données concernant ses propriétés mécaniques ont motivé l’étude de ses propriétés élastiques et plastiques. Ainsi à partir de calcul ab initio et du modèle de Peierls Nabarro, j’ai pu déterminer une forte variation de l’anisotropie élastique avec la pression, les différentes structures de cœurs des dislocations vis et coin et les systèmes de glissement préférentiels au sein de la glace X dans son champ de stabilité de 100 à 350 GPa. Nos calculs suggèrent que la déformation de la glace X est toujours localisée sur le plan {110} et que le système <110>{110} contrôle la déformation plastique en dessous de 250 GPa et que le système <100>{110} est dominant à plus haute pression. Nos résultats montrent aussi que si l’anisotropie élastique augmente rapidement avec la pression, la plasticité de la glace X devient quasi-isotrope à 350 GPa. / H2O ice is found in a variety of planetary environments, notably in the form of high pressure polymorphs inside icy moons and extrasolar ocean planets. The great diversity of thermodynamic conditions predicted inside such planetary bodies, reveals the need for new experimental and computational data to allow modeling of their internal structure and dynamics.Structural and spectral properties of H2O pure ices have been intensively studied, but surprisingly there is a lack of petrological data on impurities rich ice solid solutions. This Ph.D. thesis work focused on the study of ice VI and ice VII fusion curves in the H2O-NaCl binary, using diamond anvil cell and Raman spectroscopy. We later determined the partitioning of the NaCl analog salt, RbI, between ice VI and VII and the aqueous fluid using X- ray fluorescence and X-ray diffraction techniques at the European Synchrotron Research Facility (Grenoble). Our results enable us to observe a density inversion between ice VI and the salty fluid, and to measure a strong difference in salt partitioning between ice VI and ice VII with a partition coefficient of Kd(VI-VII)=4.5(±2.7)10-2. Inside the largest H2O rich planetary bodies, called ocean planets, the icy mantle, putatively more than 1000 km thick, shelters an ultra high pressure ice form called ice X. This H2O ice phase is unique because of its ionic crystallographic structure, in contrast with lower pressure ices polymorphs, all being molecular solids. This characteristic coupled with the fact that no data are available yet on its mechanical properties, encouraged us to study its elastic and plastic properties. Using ab initio calculations and the Peierls Nabarro model, I showed the strong variation of elastic anisotropy with increasing pressure and determined the dominant slip system inside the structure of ice X over its entire pressure stability range from 100 to 350 GPa. Our calculations suggest that plasticity in ice X is dominated by displacement always occurring on the {110} glide plane. Also, it reveals that the <110>{110} glide system is dominant below 250 GPa and that the <100>{110} slip system controls the plasticity of ice X. Our results also show that, if elastic anisotropy of ice X is strongly increasing with increasing pressure, the plasticity becomes almost isotropic at 350 GPa.

In July 2031, a new chapter in the exploration of the Jovian system will begin with the arrival of the Jupiter Icy Moons Explorer (Juice) at Jupiter. Launched on April 14 2024 as part of ESA’s Cosmic Vision programme, the mission aims to study Jupiter and its icy Galilean moons Callisto, Europa, and Ganymede. Juice carries a whole suite of instruments for in-situ and remote ground observations, one of which is the Jovian plasma Dynamics and Composition analyser (JDC). As a part of the Particle Environment Package (PEP), the particle detector will measure the energy, mass, charge and arrival direction of ions and electrons in the Jovian magnetosphere. Spacecraft charging caused by interactions between the spacecraft and its surrounding plasma environment poses a significant problem for JDC because the electrostatic potential of the spacecraft accelerates/decelerates charged particles, resulting in distorted measurements, particularly for the lowest energy particles.  In this report, we show the results of spacecraft charging and instrument simulations performed in the Spacecraft Plasma Interaction System (SPIS) for the solar wind and Earth’s magnetosheath—two environments that Juice will encounter at the start of the cruise phase. We found that the conductive surfaces that cover the majority of the spacecraft become positively charged as a result of a large photoelectron current in both the solar wind and magnetosheath environments. We show that these surfaces are expected to reach potentials of 9 V in the solar wind and 4 V in the magnetosheath. The four radiators on Juice that are covered with dielectric paint and shaded by the sun shield become negatively charged in both simulated environments. The radiator potentials can be as low as -40 V in the solar wind and -100 V in the magnetosheath. We also conclude that due to blocking by the spacecraft main body, the ion population cannot be sampled in the solar wind unless a spacecraft roll is performed. Furthermore, due to the high ion f low energy, spacecraft charging will not influence JDC measurements in this environment.  In the magnetosheath, the ion population can be sampled by JDC, and we identified three distortion mechanisms: (1) repulsion by the main body, (2) attraction by two of the radiators, and (3) repulsion by the MAG boom. Of all the distortion modes, the one originating from a negatively charged (-67.8 V) radiator close to JDC is the strongest, affecting ions with energies above 80 eV. The least powerful but most prevalent mode is the repulsion of ions by the main body. Our results can be compared with future in-situ measurements to identify distortion mechanisms well ahead of the science phase in which the scientifically important measurements will be carried out.