OPC applications

TP1 | Atmospheres and Exospheres of Terrestrial Bodies

EPSC2024-231 | ECP | Posters | TP1

Validation of Martian dust storm trajectories in the Mars PCM using observational datasets 

Demetrius Ramette, Ehouarn Millour, Tanguy Bertrand, Antoine Bierjon, and Kerstin Schepanski
Mon, 09 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P1

Dust is a key component affecting the radiative budget of the Martian atmosphere. In order to achieve accurate numerical modelling of the Martian weather and climate, it is therefore crucial to understand the dynamics of Martian dust storms throughout their whole life cycle, from the injection to the deposition of dust [1]. However, our understanding of dust storm dynamics is still incomplete and simulations of the Martian dust cycle come with large uncertainties, especially concerning the transport of dust by model winds. The latter are difficult to validate, as in-situ or remote-sensing measurements of winds in the Martian atmosphere are rare or inexistent. How well do dust storm trajectories predicted by a Mars Global Climate Model match observed trajectories?

In this study, we use recently published observations [2, 3] of dust storm trajectories and dust storm simulations by the Mars PCM6 [4], in order to compare the direction and velocity of dust transport of individual dust storms. This will allow to identify strengths and weaknesses of the Mars PCM6 with respect to dust transport.

The Mars PCM in its newest version 6 [5, 6] is able to initiate a dust storm with a given dust scenario – based on maps of the observed column dust optical depth from [7] – and then let it freely evolve, following the large scale winds and other relevant physical processes. The Mars Dust Activity Database (MDAD, [2]) and the Mars Dust Storm Sequence Dataset (MDSSD, [3]) contain compilations of observed dust storms, based on high-resolution (0.1° x 0.1°) daily mosaics of wide-angle images from orbital cameras. Fig. 1 shows an exemplary trajectory of a dust storm member from the MDAD.

This new combination of model results and observations may further open opportunities to better understand the dynamics of Martian dust storms, especially concerning emission and deposition processes and, based on that, to improve dust cycle representation in the Mars PCM.

Fig. 1: Plotting of the trajectory of an exemplary dust storm member during Mars Year 24 on a map of Mars, around Ls=185° (here: dust storm member m03_126 from the MDAD [2]). The size of the colored circles scales with the surface area of the dust storm and the color indicates the time of observation (solar longitude Ls, in degrees). The vertical black lines indicate the maximum latitudinal extension of the dust storm. The black broken line connects the individual positions of the dust storm centroid and thus indicates the evolution of its position over time. The contour lines (dashed and solid lines) represent the Martian topography with a line interval of 1500m. Here, we see that the dust storm originates in the Northern Hemisphere west of Acidalia Planitia and ends in Utopia Planitia. In this study, we want to compare these kind of tracks with model simulations.

 

[1] Kahre, M. A., Murphy, J. R., Newman, C. E., Wilson, R. J., Cantor, B. A., Lemmon, M. T., & Wolff, M. J. (2017). The Mars dust cycle. The atmosphere and climate of Mars, 18, 295.

[2] Battalio, M., & Wang, H. (2021). The Mars Dust Activity Database (MDAD): A comprehensive statistical study of dust storm sequences. Icarus, 354, 114059.

[3] Wang, H., Saidel, M., Richardson, M. I., Toigo, A. D., & Battalio, J. M. (2023). Martian dust storm distribution and annual cycle from Mars daily global map observations. Icarus, 394, 115416.

[4] Forget, F., Hourdin, F., Fournier, R., Hourdin, C., Talagrand, O., Collins, M., Lewis, S.R., Read P.L. & Huot, J. P. (1999). Improved general circulation models of the Martian atmosphere from the surface to above 80 km. Journal of Geophysical Research: Planets, 104(E10), 24155-24175.

[5] Forget, F., Millour, E., Bierjon, A., Delavois, A., Fan, S., Lange, L., Liu, J., Mathe, C., Naar, J., Pierron, T., Vandemeulebrouck, R., Spiga, A., Montabone, L., Chaufray, J.-Y., Lefèvre, F. Määttänen, A., Montmessin, F., Rossi, L., Vals, M., Gonzalez-Galindo, F., Lopez-Valverde, Wolff, M.J., Young, R., Lewis, S.R. & Read, P. L. (2022). Challenges in Mars Climate Modelling with the LMD Mars Global Climate Model, Now Called the Mars “Planetary Climate Model”(PCM). In Seventh international workshop on the Mars atmosphere: Modelling and observations.

[6] Millour, E., Bierjon, A., Forget, F., Spiga, A., Wang, C. and the Mars PCM Team, "Improving the vertical distribution of airborne dust in the Mars PCM", EPSC 2024

[7] Montabone, L., Forget, F., Millour, E., Wilson, R. J., Lewis, S. R., Cantor, B., Kass, D., Kleinböhl, A., Lemmon, L.T., Smith, M.D. & Wolff, M. J. (2015). Eight-year climatology of dust optical depth on Mars. Icarus, 251, 65-95.

How to cite: Ramette, D., Millour, E., Bertrand, T., Bierjon, A., and Schepanski, K.: Validation of Martian dust storm trajectories in the Mars PCM using observational datasets, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-231, https://doi.org/10.5194/epsc2024-231, 2024.

TP2 | Mars Surface and Interior

EPSC2024-748 | ECP | Posters | TP2 | OPC: evaluations required

Water and Sediment Transport Processes in Jezero Crater 

Anastasiia Ovchinnikova, Ralf Jaumann, Sebastian H. G. Walter, Christoph Gross, Wilhelm Zuschneid, and Frank Postberg
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P7

Introduction: The current NASA's Mars 2020 mission is exploring the Jezero crater that was once filled with water. Since it is widely acknowledged that access to liquid water is essential for life, studying the fluvial activity in Jezero can aid in the crater’s habitability assessment. We examined ten water-related processes using water and sediment transport models by [3]: 1) the western inlet valley carving, 2) the northern inlet valley carving, 3) crater flooding by only northern inlet, 4) by both northern and western inlets, 5) erosion of the western rim by the western inlet, 6) erosion of the eastern rim due to the outlet, 7) water outflow from the crater, 8) outlet valley carving, 9) western delta deposition, 10) northern delta deposition. The goal of our study was to calculate the minimum timescales for each event and estimate the minimum volume of water provided/released during each event. Relative comparison of timescales and of water amount which we calculated based on new geomorphological observations introduces a deeper understanding of the water history in Jezero.

Data: Measurements of channel sizes, valleys, deltas, eroded rims, and outflowed water volumes were based on Mars 2020 Science Investigation CTX DEM Mosaic [5] and HRSC Mars Chart (HMC) DTM and corresponding orthomosaics [2].

Geomorphological observations: The northern inlet was most likely involved in the crater flooding because it has terraces at the same height as the breaching terraces in the eastern rim (breaching happened in 3 phases, as shown in [7]). The western inlet, in contrast, has no terraces at the breaching heights. This implies that either it was not involved in crater flooding, or its terraces were eroded.

Mapping: Both deltas were mapped using three potential extents: minimum, medium, and maximum. Valleys were mapped based on two morphological features: 1) initial valley, borders of which are not visible on HMC ortho-mosaics and could only be recognized on slope and profile curvature rasters, calculated from HMC DTM; 2) last incision valley, which was mapped on HMC ortho-rectified image mosaic.

Measurements: Dimensions of the channels were derived from longitudinal and cross-sectional profiles on CTX DTM. Water volume to fill the crater before breaching and the amount of sediments, transported from valleys and deposited in deltas were estimated using ArcGIS Tools „Surface Volume“ and „CutFill“.

Methodology: Flow discharge and sediment transport models [3] are used to calculate water and sediment transport timescales under constant bank-full discharge when most of erosion occurs. The models do not include the climate and non-bank-full conditions to constrain minimum timescales. Main input parameters include: Median Grain Size, Channel depth, width and slope, Sediment porosity, Shields criterion for incipient motion, Sediment density, Martian gravity, and Water density. Data from the Perseverance rover [1] were taken for the Median Grain Size estimation. Sediment porosity, Sediment density, and Shields criterion for incipient motion were taken based on previous research [3], [4], [6]. Several scenarios were modeled with varying values of input parameters in the most expectable ranges. Eastern and western rim breechings were modelled both in a catastrophic scenario and in a long-term erosion under constant flow scenario.

Results: A comparison of the timescales of the last incised valleys carving and delta depositions showed that deltas were deposited during the last incision of the corresponding valleys. For the northern delta, the medium extent is the most probable; for the western delta, the maximum extent is the most probable.

Most modelled cases, in both long-term and catastrophic scenarios of the rim breaching, showed that the outlet valley carving lasted longer than the eastern rim erosion. Therefore, Jezero was an open-basin lake after breaching.

The eastern rim erosion and water outflow during breaching showed different results depending on the scenario. In the long-term scenario, the eastern rim erosion lasted longer than the outflow of water which was stored in the crater before breaching. That means, that this amount of water (236 km3) was not enough to carve the breach. In the catastrophic scenario there are overlapping timescales, therefore, the eastern rim breaching, and water outflow could happen simultaneously.

Multiplying the timescales by corresponding discharges allows us to calculate the minimum water volume provided/released during each event. Dividing the minimum water volume by the volume of the crater before breaching (446 km3) shows that the northern inlet as well as the western inlet could alone flood the crater (last column in Table 1).

A comparison of the minimum amount of water discharged after the breach (236 km3) with the amount of water needed to carve the whole outlet valley (1000 – 14800 km3, Table 1) confirms that Jezero must have been an open-basin lake after the breaching.

Table 1. Minimum water volume provided/released during the carving of the valleys. Timescales presented for the total valley carving (initial valley and last incised valley together).

 

Timescale, Earth years

Discharge, km3/day

Minimum volume of provided/released water, km3

 How many times Jezero could be filled before breaching (basin volume = 446 km3)

Northern Valley Total

70 – 632

~0.2

4600 – 42000

10 – 94

Western Valley Total

505 – 3792

~0.2

33500 – 253000

75 – 567

Outlet Valley Total

1.4 – 21.2

~1.9

1000 – 14800

 

References:

[1] Farley K., and Stack K. (February 15, 2023). Mars 2020 reports, Volume 2. https://mars.nasa.gov/internal_resources/1656/

[2] Gwinner, K. et al. (2016). The High Resolution Stereo Camera (HRSC) of Mars Express and its approach to science analysis and mapping for Mars and its satellites. Planetary and Space Science, 126, 93–138. https://doi.org/10.1016/j.pss.2016.02.014

[3] Kleinhans, M. G. (2005). Flow discharge and sediment transport models for estimating a minimum timescale of hydrological activity and channel and delta formation on Mars. Journal of Geophysical Research: Planets, 110(12), 1–23. https://doi.org/10.1029/2005JE002521

[4] Kleinhans, M. G., van de Kasteele, H. E., & Hauber, E. (2010). Palaeoflow reconstruction from fan delta morphology on Mars. Earth and Planetary Science Letters, 294(3–4), 378–392. https://doi.org/10.1016/j.epsl.2009.11.025

[5] Malin, M. C. et al. (2007). Context Camera Investigation on board the Mars Reconnaissance Orbiter. Journal of Geophysical Research: Planets, 112(5). https://doi.org/10.1029/2006JE002808

[6] Roda, M. et al. (2014). Catastrophic ice lake collapse in Aram Chaos, Mars. Icarus, 236, 104–121. https://doi.org/10.1016/j.icarus.2014.03.023

[7] Salese, F. et al. (2020). Estimated Minimum Life Span of the Jezero Fluvial Delta (Mars). Astrobiology, 20(8), 977–993. https://doi.org/10.1089/ast.2020.2228

How to cite: Ovchinnikova, A., Jaumann, R., Walter, S. H. G., Gross, C., Zuschneid, W., and Postberg, F.: Water and Sediment Transport Processes in Jezero Crater, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-748, https://doi.org/10.5194/epsc2024-748, 2024.

EPSC2024-777 | ECP | Posters | TP2

Polyphase Amazonian floods in the Olympica – Jovis Fossae channel system. 

Anita Zambrowska, Daniel Mège, and Sam Poppe
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P14
  • Centrum Badań Kosmicznych Polskiej Akademii Nauk (CBK PAN), Bartycka 18A, 00-716 Warszawa, Poland
IntroductionIntroductio

Introduction

Geomorphological evidence of deglaciation was discovered on Mars at mid-latitudes [1, 2], in response to planetary obliquity changes [3]. Many craters contain evidence of subsurface ice remaining from Late Amazonian deglaciation [4]. Martian climate models predict the presence of ice accumulation in the mid-latitudes, especially on the western flanks of the Tharsis montes [5]. During the Amazonian epochs, the early Tharsis reservoirs were partly withdrawn to the surface, through outflow channels [6]. The largest of them formed between 3.7 and 3 billion years ago. There is evidence, however, that the crustal water reservoir was not exhausted after this time [7][8]. Many channels, including the largest Amazonian system Olympica – Jovis Fossae, formed at the top of the Tharsis bulge in lava flows that may be only a few hundred million years old [9]. To investigate the formation of such channel systems, we present a new geological map of Olympica-Jovis Fossae (Figure 1).

Data and methods

To map morphological and geological features, we used visible imagery from the Mars Reconnaissance Orbiter (MRO) Context Camera (CTX) with a resolution of 6 meters per pixel [10][11] and MRO High-Resolution Imaging Science Experiment (HiRISE) camera with a resolution of 25 cm per pixel [12]. To constrain the intersecting relationships between channels we used the newest DEM topography data (2021 – 2023) from the Mars Express High Resolution Stereo Camera (HRSC) with a resolution of ~20-40 meters per pixel [11][13] and optical imagery data from the ExoMars Trace Gas Orbiter, Colour and Stereo Surface Imaging System (CaSSIS) with a resolution of 4.8 meters per pixel ([14]).

Results

Our new geological map includes fluvial units (Figure 2), volcanic units, volcanic edifices and vents, glacial lineations, regional fracture systems, and impact craters. We identified 28 fluvial units that could be attributed to six classes, based on discriminating geomorphological features. Observed channel morphologies include the main channel (set of troughs and channels extending from the triple junction to Jovis Fossae), braided channels with densely interconnected channels with dendritic and sinuous geometry, terraced channels with eroded floors, step-like landforms, meandering channels, tectonically controlled channels, and outflow channels with streamlined islands. The oldest in the sequence of fluvial channels are terraced channels, representing multiple flooding events, followed by braided channels, meandering channels, outflow channels that originate at

Jovis Fossae system fissures, and tectonically controlled channels (attributed to extensional stress, subsequently developed and filled by water and lava). Investigation of the mapped units and correlation of channel formation reveal at least three bigger flooding events. We determined the relative ages and relationships between fluvial channels from cross-cutting relationships (Figure 3).

Figure 1. Geologic map of the Olympica Fossae and Jovis Fossae hydrologic system.

Figure 2. Zoom on two regions of the map. Legend: clt – collapsed lava tube, Atc – Amazonian terraced channel, Amc – Amazonian meandering channel, Abc – Amazonian braided channel, Atcc – Amazonian tectonically controlled channel.

Figure 3. Chart of cross-section relationships between 28 fluvial channels. Arrowheads indicate younger flows.

Discussion and conclusion

We present the first detailed geological map and relative dating of the Olympica Fossae and Jovis Fossae, providing a basis for understanding the recent hydrologic and volcanic activity at the top of the Tharsis dome. Our map includes 28 individual fluvial units, glacial lineations, impact craters, volcanic units, edifices, and vents, as well as geological structures. We found a series of complex cross-cutting relationships between channels. General investigation of the channel’s morphology and their correlation with volcanic edifices reveal that the Olympica Fossae and Jovis Fossae formed by repeated overflows from various fissures related to an intensive period of distributed volcanism. Olympica – Jovis Fossae shares similar fluvial morphology with other outflow channels on Mars. Since our geologic map reveals six types of channels, future detailed morphologic analysis will be crucial to better understand the general conditions for catastrophic groundwater flooding in a volcanic context. Further efforts are needed to identify the chronological succession of deposition and flows and interfaces between floods, lava flows, and airfall deposits, based on crater retention age determination.

References:

[1] Baker, D. M. H., & Head, J. W. (2015). Icarus, 260, 269–288.

[2] Head, J. W., et al. (2006). Geophysical Research Letters, 33(8).

[3] Head, J. W., et al. (2006). Earth and Planetary Science Letters, 241(3-4), 663–671.

[4] Dickson, J.L., et al. (2010). Earth and Planetary Science Letters, 294(3-4), 332–342.

[5] Shean, D. E. (2005). Journal of Geophysical Research, 110(E5).

[6] Cassanelli, J.P., Head, J.W. (2019). Planetary and Space Science 169, 45-69.

[7] Rodriguez, J. A. P., et al. (2015). Icarus, 257, 387–395.

[8] Hiatt, E., et al. (2024). Icarus, 408, 115774

[9] Keske, A.L., et al. (2015). Icarus, 245, 333-347.

[10] Fergason, R. L., et al. (2021). 51st  Lunar Planet. Sci. Conf. Abstract #2020.

[11] Jaumann, R., et al. (2007). Planetary and Space Science, 55(7), 928–952.

[12] McEwen, A. S., et al. (2007). Journal of Geophysical Research, 112(E5), E05S02.

[13] Neukum, G., Jaumann, R. (2004). Eur. Space Agency Spec. Publ., ESA-SP 1240, 17–35.

[14] Thomas, N., et al. (2017). Space Sci Rev 212, 1897–1944.

How to cite: Zambrowska, A., Mège, D., and Poppe, S.: Polyphase Amazonian floods in the Olympica – Jovis Fossae channel system., Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-777, https://doi.org/10.5194/epsc2024-777, 2024.

EPSC2024-1077 | ECP | Posters | TP2

Evolution of fluvial, and possible lacustrine or marine activity in Nilosyrtis, Mars -  A geomorphological & chronostratigraphic analysis. 

Emelie Saar, Cynthia Sassenroth, and Andreas Johnsson
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P12

 Introduction

Today, Mars’ surface is defined as a cold hyper-arid desert, but geomorphological and geochemical observations have shown that in the past, vast amounts of liquid water were stable at the planet’s surface (Di Achille et al. 2006, Erkeling et al. 2012, Hauber et al. 2009). With the increased coverage of high-resolution remote sensing datasets, investigations of possible remnants of fluvial processes in the Martian landscape have significantly increased during the last couple of decades. The spatiotemporal evolution of hydrological events is however not well understood (Erkeling et al. 2012). In this master thesis study, an impact crater in the region of Nilosyrtis Mensae has been mapped to identify geomorphological evidence that can increase the understanding of the spatiotemporal evolution of landforms at the dichotomy boundary on Mars.

Method & Data

The mapping was made using ArcGIS Pro version 3.1.0. with inspiration taken from the FGDC Digital Cartographic Standard for Geologic Map Symbolization (USGS, 2006) as well as the Italian cartographic standard Geomorphological Map of Italy at 1:50,000 scale by ISPRA (2021). The Context Camera [CTX] mosaic & raw images (5-6 m/px) and High Resolution Imaging Science Experiment [HiRISE] images (0.3-0.25m/px) as well as the THEMIS IR mosaic dataset of 100m/px were used to map surface features. 3D models were used in this study to support interpretations, such as the HRSC DEM product mosaic of 50m/px. By using the MarsSI application (marssi.univ-lyon1.fr) overlapping CTX raw images have been processed for HiRISE & CTX imagery to create DEM products with 12 m/px and 0.5 m/px respectively.

Geological setting

The study crater is 80+ km wide and is located at the dichotomy boundary. The crater lacks primary impact morphology and ejecta. The northeastern wall is completely degraded and opens up towards the northern plains. Moreover, the crater wall is incised by several channels of which two channels connects to adjacent craters of similar size, located east and west of the study area. The crater floor forms a palimpsest landscape of different landforms suggesting a very complex history of sedimentation, exhumation and erosion.   

Result

The main result is the geomorphological map produced by this study (Fig.1). Key findings include a set of fan-shaped landforms that partially cover the crater floor. Some fans are associated with valley incisions and one of the larger fans is heavily affected by impact cratering. The fan shapes have a variety in surface texture, morphology, size and stratigraphic positions. Multiple elongated fan shapes combined with streamlined plains, alluvial deposits have been identified at the termination of the southwestern valley inlets (Fig. 2).

Figure 1. A geomorphological map of an impact crater basin in Nilosyrtis Mensae, Mars, at approximate
29°N, 72°E draped over the CTX satellite mosaic (NASA/JPL/MSSS/The Murray lab
Spatial reference: CGS Mars 2000 Sphere. Datum: Mars 2000 sphere. Map units: Degree.

Figure 2. A) 3D view of the valley inlet entering the study crater from the southeast. B) An overview of the mapped units of possible fluvial genesis in the western region of the impact crater basin. Fan shape 1 and 2 are deposited in a close spatial proximity to the channel inlet, while fan shape unit 3 is found further out, divided by a depression of possible alluvial deposition. Fan shape unit 4 covers a large part of the crater floor and does exhibit unique geomorphological features at the top of the streamlined plains marked at 5. Image credit: NASA/JPL/MSSS/The Murray lab.

 Discussion & Conclusion

The diverse classification of landforms and interpreted geneses shows that the studied crater has experienced a long history of degradation, deposition by various processes and erosion since its formation. The superposition between geological units suggests a scenario where fluvial deposition have occurred episodically. Bamberg et. al., (2014) describes the geological history of the floor-fractured crater which is connected to the studied crater by a channel inlet from the northwest. The possible fluvial landforms in the studied crater suggests a close relationship to the hypothesized scenario by Bamberg et al., (2014) of a major fluvial event. Lesser fluvial activity followed, predominantly controlled by conditions in the study crater.  In conclusion, the geomorphology suggests a long and complex history of fluvial events. An impact of possible lacustrine or marine conditions could not be confirmed nor refuted. The framework of the thesis did not allow for more detailed studies using spectral data. Work that includes CRISM may provide further insights into the temporal evolution in the studied crater.   

References

Bamberg et.al., (2014). Floor-Fractured Craters on Mars – Observations and Origin. Planetary and Space Science, 98, 146–162. https://doi.org/10.1016/j.pss.2013.09.017

Di Achille et.al., (2006) Geological evolution of the Tyras Vallis paleolacustrine system, Mars. Journal of Geophysical Research, 111(E4), E04003-N/a.

Erkeling, et.al., (2012). Valleys, paleolakes and possible shorelines at the Libya Montes/Isidis boundary; implications for the hydrologic evolution of Mars. Icarus 219(1), 393-413.

Hargitai, (2019). Planetary Cartography and GIS. Springer.

Hauber, et.al., (2009). Sedimentary deposits in Xanthe Terra: Implications for the ancient climate on Mars. Planetary and Space Science, 57(8), 944-957

ISPRA. (2021). Geomorphological Map of Italy at 1:50,000 scale - Update and additions to the Guidelines of the Geomorphological Map of Italy at 1:50,000 scale and Geomorphological Database (Version 2.0). Volume 13, issue 1 Volume 13, Part I - The CARG Project: changes and additions to the Quaderno n. 4/19944/2022. https://www.isprambiente.gov.it/resolveuid/45c17b95e93f45cda3e1c2fa0f5e06fa

USGS. (2006). FGDC Digital Cartographic Standard for Geologic Map Symbolization (PostScript Implementation. U.S. Geological Survey Techniques and Methods 11-A2. Derived 2023-12-29 from:https://pubs.usgs.gov/tm/2006/11A02/

 

How to cite: Saar, E., Sassenroth, C., and Johnsson, A.: Evolution of fluvial, and possible lacustrine or marine activity in Nilosyrtis, Mars -  A geomorphological & chronostratigraphic analysis., Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1077, https://doi.org/10.5194/epsc2024-1077, 2024.

EPSC2024-1115 | Posters | TP2 | OPC: evaluations required

Exploring Polygonal Patterned Grounds in the Hyper-arid Atacama Desert: Insights into Formation Mechanisms and Implications for Martian Analogues 

Josephine Anderson, Christof Sager, Alessandro Airo, Andrea Miedtank, Franziska Schwonke, and Jenny Feige
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P18

Understanding the formation mechanisms of polygonal patterned grounds in hyper-arid environments like the Atacama Desert is crucial for unraveling their significance as an analogue for extraterrestrial landscapes, notably Mars. Unlike their extensively studied periglacial counterparts, the genesis of hyper-arid polygons requires more detailed exploration, since multiple processes, such as, desiccation [1] (including gypsum dehydration), thermal contraction [2], and haloturbation [3] are currently under debate. Therefore, we excavated a trench in polygonated ground in the Yungay area of the Atacama Desert and systematically collected over 80 soil samples covering an entire polygon, its adjacent sand wedges, and underlying sediments (190 x 190 cm). Sedimentological and geochemical analyses indicate a continuous interplay of atmospheric salt deposition, dissolution, precipitation, and inducing clast heave and shattering, under extreme dry conditions. Our findings emphasize that polygons in the Atacama Desert provide a promising analogue to Martian saline polygons, emphasizing the relevance of terrestrial studies for advancing our understanding of extraterrestrial landscapes.

References:

1 Ewing 2006

² Sager 2021

³ Zineabelidin 2024 (preprint)

How to cite: Anderson, J., Sager, C., Airo, A., Miedtank, A., Schwonke, F., and Feige, J.: Exploring Polygonal Patterned Grounds in the Hyper-arid Atacama Desert: Insights into Formation Mechanisms and Implications for Martian Analogues, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1115, https://doi.org/10.5194/epsc2024-1115, 2024.

EPSC2024-1355 | ECP | Posters | TP2 | OPC: evaluations required

Spatiotemporal development of two stepped fans in Xanthe Terra and Terra Sirenum, Mars 

Cynthia Sassenroth, Ernst Hauber, Maria Cristina Salvatore, and Carlo Baroni
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P9

INTRODUCTION

Deltas on Mars are crucial markers for reconstructing the planet's past climatic and hydrological conditions. The detailed study of Martian deltas, particularly through the interpretation of high-resolution imaging and geomorphological mapping, reveals insights into the environmental changes the planet has undergone. This study investigates stepped-fan deposits on Mars. The research utilized remote sensing techniques and landscape analysis to produce geomorphological maps of two key sites: Picardi crater in Terra Sirenum and Dukhan crater in Xanthe Terra (Fig. 1). These sites showcase extremely well-preserved stepped-fan deposits, possibly have formed during groundwater sapping events in the early Amazonian [1], [2].

 

Figure 1: A) Location of stepped-fan deposits in Picardi (Terra Sirenum) and Dukhan Crater (Xanthe Terra). B) Stepped-fan in Picardi crater. C) Stepped-fan deposits in Dukhan crater. 

 

METHODS

We used data from NASA's MRO mission (HiRISE and CTX imagery and DEMs) to produced detailed geomorphological maps of the study sites. The maps were generated using a combination of two geomorphological keys. The planetary symbolization key by the Federal Geographic Data Committee of the United States (U.S.G.S., 2006) lacks the representation of small-scale features like alluvial fans, deltas, and various surface textures observed at the sites. Therefore, a synergistic mapping approach was adopted by also integrating the symbology of the cartographic standard proposed by the Italian Geological Survey and the Italian Association of Physical Geography and Geomorphology (Gruppo Nazionale Geografia Fisica e Geomorfologia CNR, 1986).

RESULTS

 

Picardi crater showcases an exceptionally well-developed  stepped-fan deposit [3], [4]. A multitude of geomorphological features and deposits, related to impact cratering, paleochannel activity, aeolian erosion and deposition, mass wasting and tectonism, have been identified . Most intriguing are the stepped, fan-shaped deposits in the southeastern part of Picardi crater, which have been further subdivided into individual geomorphic units, each differentiated by their specific morphometrical and sedimentological characteristics, and potential origin. The fan is connected to a small, steep-walled amphitheater-headed channel. Seven distinct fan units, originating from a single amphitheater-headed valley, were deposited radially outward from the apex (Fig. 2). 

 

Figure 2: Detailed HiRISE images of the Picardi crater fan deposit. A) Fan Unit 1. Wrinkle ridge and smooth, sparsely cratered sediments of Fan Unit 1 distal to the main fan-geological units. B) Fan Unit 5. Finger-like extensions at the topographically low, distal part of the main fan. C) Secondary, isolated fan deposit. D) Upper part of the fan deposits with Fan Units 6 and 7 covering the highly eroded, semicircular Fan Unit 4. North is up in all panels.

Dukhan crater exhibits a diverse array of features, including plains, plateaus, tectonic structures, impact craters, and various depositional features such as fluvial, gravitational, and aeolian deposits (see figure 1c and 3). The stepped-fan deposits in Dukhan crater consist of distinct units with varied characteristics. Fan Unit 1, the most distal, has a wedge-like shape with sparse cratering and large boulders on its flat top. Fan Unit 2 is partially covered by Unit 4 and displays a slightly higher albedo with erosion on its distal side. Fan Unit 3, with a smooth surface and minor disruptions, forms the bulk of the fan. Unit 4, at the uppermost part, features asymmetrical geometry and erosion on its eastern side. The fan is connected to an amphitheater-headed inlet valley, aligned with a wrinkle ridge system on the highland plateau, suggesting underlying structural control. Along the crater walls, multiple channel-like features indicate the presence of alluvial deposits and paleochannels, gradually transitioning from rugged terrain to linear paths with changes in albedo.

Figure 3: Detailed HiRISE images of the Picardi crater deposit. A) Fan Unit 1. Wrinkle ridge and smooth, sparsely cratered sediments of Fan Unit 1 distal to the main fan-geological units. B) Fan Unit 5. Finger-like extensions at the topographically low, distal part of the main fan. C) Secondary, isolated fan deposit. D) Upper part of the fan deposits with Fan Units 6 and 7 covering the highly eroded, semicircular Fan Unit 4. North is up in all panels.

DISCUSSON

The stepped-fan deposits in Picardi crater exhibit a complex evolution, transitioning from deltaic depositional in the lower fan to alluvial dominance in the upper fan. In the middle section Unit 4 (fig. 2d), characteristics of glacial activity were recognised. The formation mechanisms of these units suggest a dynamic interplay of environmental factors, including shifts in climate and hydrological activity as well as glacial activity.

In Dukhan crater, the lowermost unit displays large boulders on its surface with diameters of >1 m and displays a massive sedimentary unit. It likely formed as a result of a landslide, and initiated the formation of the amphitheater-headed channel as a result of groundwater sapping. The upper units, resembling conical, layered deposits, have characteristics of both alluvial fans and deltaic deposits, making their differentiation challenging. While they share similarities with deltaic deposits found in other Martian regions, their lack of typical features suggests an alternative formation process, possibly through alluvial mechanisms driven by groundwater aquifers.

 

CONCLUSIONS

The study reveals a more complex formation history for Martian stepped-fans than previously recognized. Both Picardi and Dukhan craters show evidence of multiple processes, including deltaic, alluvial, and glacial and mass wasting processes. These findings enhance our understanding of Mars' climatic and hydrological evolution, indicating shifts from water-rich to glacial and arid conditions over time.

REFERENCES

[1]             E. Hauber et al., ‘Asynchronous formation of Hesperian and Amazonian-aged deltas on Mars and implications for climate’, Journal of Geophysical Research: Planets, vol. 118, no. 7, pp. 1529–1544, 2013, doi: 10.1002/jgre.20107.

[2]             E. Hauber et al., ‘Old or not so old: That is the question for deltas and fans in Xanthe Terra, mars’, presented at the Third Conference on Early Mars, 2012.

[3]             E. R. Kraal, M. van Dijk, G. Postma, and M. G. Kleinhans, ‘Martian stepped-delta formation by rapid water release’, Nature, vol. 451, no. 7181, pp. 973–976, Feb. 2008, doi: 10.1038/nature06615.

[4]             G. G. Ori, L. Marinangeli, and A. Baliva, ‘Terraces and Gilbert-type deltas in crater lakes in Ismenius Lacus and Memnonia (Mars)’, Journal of Geophysical Research: Planets, vol. 105, no. E7, pp. 17629–17641, 2000, doi: 10.1029/1999JE001219.

How to cite: Sassenroth, C., Hauber, E., Salvatore, M. C., and Baroni, C.: Spatiotemporal development of two stepped fans in Xanthe Terra and Terra Sirenum, Mars, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1355, https://doi.org/10.5194/epsc2024-1355, 2024.

TP3 | Planetary field analogues for Space Research

EPSC2024-124 | ECP | Posters | TP3 | OPC: evaluations required

Rivers' classification: integrating Deep Learning and statistical techniques for terrestrial and extraterrestrial drainage networks analysis 

Mariarca D'Aniello, Maria Raffaella Zampella, Andrea Dosi, Alvi Rownok, Michele Delli Veneri, Adriano Ettari, Stefano Cavuoti, Luca Sannino, Massimo Brescia, Carlo Donadio, and Giuseppe Longo
Tue, 10 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P10

Various approaches have been proposed to describe the geomorphology of drainage networks and the intricate relationships between abiotic/biotic factors and their surrounding environment. There is an intrinsic complexity of the explicit qualification of the morphological variations in response to various types of control factors and the difficulty of expressing the cause-effect links. Traditional methods of drainage network classification are based on the manual extraction of key characteristics, subsequently applied to pattern recognition frameworks. These attitudes, however, have low predictive and uniform ability. For this reason, we present a different approach, based on the data-driven supervised learning by images, extended also to extraterrestrial cases. Using deep learning models, the extraction and classification phases are integrated within a more objective, analytical, and automatic toolkit (Donadio et al., 2021).

 

Pre - processing of satellite and topographical images through image segmentation methods

Extraterrestrial and terrestrial drainage pattern analysis is a topic of central interest since it allows scientists to understand the hydrogeological and geomorphological past of planets and satellites. Earth, Mars, Venus, and Titan’s patterns have been taken into consideration in this work. The extraction process, necessary to obtain an outline of the river, can be done through image segmentation, considering different mathematical methods whose efficiency varies according to the conditions and properties of the images. The segmentation methods used can reliably identify objects’ contours if in the foreground, separating them from the background. In the context of image processing, this is addressed as Edge Detection.

Two types of images were addressed, respectively, topographic and satellite. In both, an extensive pre-processing phase has been carried out, to reduce background noise with computationally efficient and optimized algorithms. This makes the profiles of the drainage networks stand out from the rest of the image, minimizing the loss of important information and the need for human intervention.

This aims to make the preparatory phase of the images (pre-processing) as self-consistent as possible, to be effectively applied to large volumes of images, allowing the generation of a valid training set for the classification of drainage patterns using self-adaptive methods, based on machine and deep learning paradigms.

In the final work, a good trade-off has been achieved between efficiency and effectiveness of the edge-detection methods. As will be discussed later, the need for an expert’s intervention is extremely limited, in most cases not needed at all.

 

River Zoo survey and classification based on Deep Learning models

This work introduces an innovative approach to river hydrographic basins classification within the River Zoo Survey project. The main goal is to perform a statistical evaluation of the classification of terrestrial and extraterrestrial drainage networks by human experts to be subsequently used as base of knowledge to train supervised Artificial Intelligence (AI) methods.

The idea is to analyze the degree of reliability of class assignment to drainage samples, driven domain expert decisions, based on visual inspection of images and the identification of the right pattern type. Through the analysis of Earth, Mars and Titan’s rivers, experts were asked to classify rivers into one of ten distinct patterns, further categorized into two macro-classes: dendritic and non-dendritic.

Figure 1 – Different classes of drainage patterns: a) dendritic; b) sub-dendritic; c) pinnate; d) parallel; e) radial; f) rectangular; g) trellis; h) angular; i) annular; j) contorted. (a)–(c) patterns are related to dendritic forms (D), (d)–(j) to non-dendritic ones (ND).

The purpose of this study is to establish an objective classification system for rivers, improving the understanding of terrestrial and extra-terrestrial rivers drainage networks. Using statistical techniques, the study explores methods to reduce noise in human based classification, thus providing a robust classification system for the automatic processing of river data with a detail much better (10 classes) than the current two class classification systems (dendritic and non-dendritic).

This work focuses on the methodology and objectives of the research, highlighting its interdisciplinary nature and potential contributions to a better comprehension of river morphologies across different planetary bodies.

 

Classification of drainage patterns using properties of fractals

Machine Learning (ML) models often require a differentiated and big enough training sequence so that they can output a good prediction rule, a predictor, to then use in labeling unseen elements belonging to the testing set. While experts can give their opinions on a given river in relation to the ten classes here considered, it is more of a subjective truth than a ground truth. To minimize the model’s bias, the training sequence should contain data as accurately labeled as possible, thus having a great probability of minimizing the true error.

Introducing fractal geometry, involving self-similar objects with a fixed degree of complexity, often referred to as Hausdorff Dimension (HD). By computing the HD for a given set of rivers, they can be grouped into classes by defining step thresholds determining the belonging to any of the ten categories. The classification can be further refined by considering the Horton-Strahler number, which is a way of establishing a hierarchy between tributaries in a drainage network. This makes it possible to keep track of the branching of rivers, along with their intrinsic fractal complexity.

To compute the HDs, different methods will be used, leveraging the flexibility and the capabilities of the Python programming language. The best algorithm to compute the HD on rivers, Box Counting (BC), will be analyzed, and compared to the experts’ classification, to further comprehend the thought process of a human mind when presented with a classification task involving complex and branched structures.

Results show that fractal analysis is reliable in the context of geomorphology and river patterns, allowing for the creation of a ground truth for the RiverZoo images, and laying the basis for the development of advanced ML algorithms used for classification purposes.

How to cite: D'Aniello, M., Zampella, M. R., Dosi, A., Rownok, A., Delli Veneri, M., Ettari, A., Cavuoti, S., Sannino, L., Brescia, M., Donadio, C., and Longo, G.: Rivers' classification: integrating Deep Learning and statistical techniques for terrestrial and extraterrestrial drainage networks analysis, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-124, https://doi.org/10.5194/epsc2024-124, 2024.

EPSC2024-695 | ECP | Posters | TP3 | OPC: evaluations required

Mapping of Martian geomorphologies as a contribution to construct the first artificial analog in Colombia 

Laura Romero, Javier Suarez, Oscar Ojeda, Yael Mendez, and Luis Ochoa
Tue, 10 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P2

1. Introduction

The development and validation of specific procedures for scientific activities in planetary analog sites, whether natural or artificial, are key to ensure the success, effectiveness, and safety in future extraterrestrial human exploration missions (Foucher et al. 2021). However, this has been developed in facilities and sites generally inaccessible for the Latin American scientific community and constitutes an obstacle for contributing to the study in analogs.

This work presents the geomorphological assessment of three representative landforms on Mars: a lava flow field at the southwestern base of Olympus Mons (13.18° N, 134.95 W), a lacustrine deposit in a paleolake within an impact crater near Maja Valles (5.33° N, 58.58° W) and the impact crater Catota in Acidalia Planitia (51.65° N, 25.98° W); which were used to test the analogy of morphological simulations of these landforms created in the first artificial Mars analog facility in Latin America, called "Rock Yard" within the Simulated Analog Space Exploration Habitat in Colombia (HAdEES-C).

2. Data

In order to carry out the geomorphological mapping and interpretation, images from the High Resolution Imaging Science Experiment (HiRISE), Mars Reconnaissance Orbiter (MRO) Context Camera (CTX) and the High Resolution Stereo Camera (HRSC) as well as digital terrain models available were collected and processed for each of the Mars’s sites, in addition to basemaps of Mars.

3. Results

3.1. Mapping of landforms

For each of the real sites on Mars, a geomorphological map was made and interpreted.

The area in the southernmost part at the base of Olympus Mons forms a flow field in an anastomosing network (Morris & Tanaka, 1995). Flows are broadly divided into three categories according to their morphology and stratigraphy (Peters et al. 2021). Corrugated flows (Cof) are the younger flows characterized for their rough surface texture resembling ‘a’a¯ flows in Hawaii (Hargitai & Kereszturi, 2015). Underlying are the Channelized flows (Caf), numerous and narrower flows with a smooth-textured central channel and rough-textured levees. And Ancient corrugated flows (Acf) are the oldest flows. These define at least three different episodes of volcanic activity with different emplacement styles. Figure 1 shows this map.

The possible paleolake is in a ~37 km complex impact crater and forms a closed-basin lake (Cabrol & Grin, 1999) with a >30 km long incised Inlet valley (Iv) (Goudge et al. 2015) but no outlet valley. It has an inlet continuous delta observed extending over 4 km, where different depositional events were identified based on the overlapping Deltaic slopes (Ds). The basin floor is nearly level given that the lacustrine deposits filled the topographic irregularities, but only resurfacing units are visible: Volcanic mantles (Vm) at east, and dark and finer Eolian deposits (Ed) at west.

The impact crater has a diameter of 1.3 km, is bowl-shaped and its rim is raised. The inner flank shows a Spur-and-gully (Sg) morphology (Watters et al. 2015) in the upper part and talus deposits in the lower part. The outer flank is covered by the preserved ejecta that divides in Continuous ejecta (Ec) and Discontinuous ejecta (Ed).

Figure 1: Geomorphological map of lava flows at the base of Olympus Mons.

3.2. Analogue model constructed

Following a low-cost and open science concept, the Rock Yard was built with two 4 m2, one 8 m2 and 18 m2 interconnected modules, each one with a specific landform made with easily accessible local construction materials to mechanically simulate the Martian soil. The cost did not exceed 400€ and it was built in 8 days. Figure 2 contains a photograph of the complete Rock Yard.

Figure 2: Rock Yard at HAdEES-C.

4. Discussions and conclusions

The recreated landforms represent approximately 44% of the Martian surface compared to the latest defined geological units on Mars defined by Tanaka et al. (2014). Several of the defined geomorphological units can indeed be correlated with the structures present in the simulations of the Rock Yard, and the main differences correspond to inconsistencies due to the scale, and the impossibility to recreate secondary surface processes on Mars that have modified the original landforms. Therefore, the Rock Yard constitutes a characterized facility whose main parameters are controlled with a geomorphologic analogy that is coherent within the limits of the scale of representation and the composition of the materials.

Work is currently underway to expand the area to include simulations of other landforms and generate interactions between them that represent a realistic geological history, to make this space a more rigorous scientific representation of the structures and processes on Mars. In any case, while being a first approximation to the design and construction of these facilities in Colombia, the Rock Yard is already a compact, representative, and accessible geologic analog. One that allows research, outreach, and educational activities that will contribute to the local production of scientific knowledge, and to bring the study of analogues to a wider public.

5. Acknowledgements

This project has received funding from Cydonia Foundation and GMAS and initiated as an undergraduate geology dissertation in the National University of Colombia.

6. References

  • Cabrol, N., & Grin, E. (1999). Distribution, Classification, and Ages of Martian Impact Crater Lakes. Icarus, 142, 160–172.
  • Foucher, F., Hickman-Lewis, K., Hutzler, A., Joy, K. H., Folco, L., Bridges, J. C., Wozniakiewicz, P., Martínez-Frías, J., Debaille, V., Zolensky, …, Westall, F. (2021). Definition and use of functional analogues in planetary exploration. Planetary and Space Science, 197.
  • Goudge, T., Aureli, K., Head, J., Fassett, C., & Mustard, J. (2015). Classification and analysis of candidate impact crater-hosted closed-basin lakes on Mars. Icarus, 260, 346-367.
  • Hargitai, H., & Kereszturi, Á. (2015). Encyclopedia of Planetary Landforms. Springer. doi:10.1007/978-1-4614-3134-3.
  • Morris, E., & Tanaka, K. (1995). Geologic Maps of the Olympus Mons Region of Mars. USGS Astrogeology Science Center.
  • Peters, S. I., Christensen, P. R., & Clarke, A. (2021). Lava Flow Eruption Conditions in the Tharsis Volcanic. Journal of Geophysical Research: Planets, 126.
  • Watters, W., Geiger, L., Fendrock, M., & Gibson, R. (2015). Morphometry of small recent impact craters on Mars: Size and terrain dependence, short-term modification. J. Geophys. Res. Planets, 120, 226–254.

How to cite: Romero, L., Suarez, J., Ojeda, O., Mendez, Y., and Ochoa, L.: Mapping of Martian geomorphologies as a contribution to construct the first artificial analog in Colombia, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-695, https://doi.org/10.5194/epsc2024-695, 2024.

EPSC2024-914 | ECP | Posters | TP3 | OPC: evaluations required

Towards Low Earth Orbit Exposure Experiments on the ISS -Designing a Simulation Setup for Mars Like Conditions 

Ruben Nitsche, Severin Wipf, Lucas Bourmancé, Adrienne Kish, and Andreas Elsaesser
Tue, 10 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P12

Abstract
In the search for life, Mars is considered to be a major target due to its similarity and relative
proximity to Earth, which makes it accessible for scientific investigation. Considering the past
chemical, geological and physical environment on Mars, the planets surface might have been
habitable to life during the so called Noachian¹. The quest to identify complex organic molecules on
the surface of Mars is an ongoing effort using instruments like SHERLOC onboard the NASA
Perseverance Rover² or the SAM³ and CheMin⁴ instruments onboard the NASA Curiosity Rover.
Other investigations of possible biosignatures on Mars are focusing on the search for chemical
processes exclusive to life. Potential indicators are specific atmospheric gases like Methane or
mineralogical signatures in composition or morphology that indicate past or present presence life.
Recent measurements indicate not only the presence of frozen but also subsurface liquid water on
Mars. Such water could only be stable on the planet in the form of highly concentrated brines.
Halophilic organisms, that are known to survive environments with high salinity, have therefore
become a focus for Astrobiology in the context of Mars⁵.
The atmosphere of Mars mostly consists of CO2 (95%), but oxygen (0.174%) and water vapor (0.03%,
variable) are also present⁶ at an atmospheric pressure of around 6 mbar. The high abundance of CO2
blocks UV radiation below around 200 nm while any UV light at higher wavelengths reaches the
surface of Mars. This is in clear contrast to solar radiation on the surface of Earth where UV light
below around 300 nm is blocked from reaching ground level due to the higher concentration of
oxygen and ozone. Additionally, the absence of magnetic shielding around Mars means that
energetic particle radiation can reach the surface of Mars. The average surface temperature of Mars
is considered to be around −63 ℃, reaching up to 20 ℃ in the equatorial regions and go as low as
−153 ℃ at the poles, with daily variations often exceeding 80 ℃⁷. The surface of Mars is covered in a
fine, unconsolidated regolith mostly originating from eroded volcanic rocks exhibiting a distinct red
color caused by high abundances of iron oxides. Varying amounts of phyllosilicates have been found
indicating the past presence of water⁸. To investigate the photochemistry of possible biosignatures in
a laboratory or space born context it is necessary to reproduce these extreme conditions as
accurately as possible.
A number of radiation exposure experiments under Mars-like conditions in Low Earth Orbit (LEO)
involving organic molecules and other astrobiological samples have been performed or are currently
under development. Considering high costs and limited availability of space born experiments we
have developed a laboratory based Mars simulation setup. Our setup partly reuses concepts of LEO
experiments while adding simulation parameters that are not yet possible to recreate in LEO due to
their technical complexity. Specialized reaction cells have been developed for the NASA O/OREOS
cube satellite experiments⁹. They hold samples, applied as thin films, in a sealed gas volume while
being transparent to irradiation and spectroscopy measurements. These reaction cells are also
planned to be used in the upcoming LEO experiments ExoCube Chem and OREOCube¹⁰ outside the
International Space Station (ISS). The reaction cells consist of a central stainless steel ring, sealed
using indium rings with a sample window on either side The window materials are chosen to allow
both irradiation and transmission spectroscopy from the UV up to the IR range. Using FTIR
spectroscopy, we can show that the reaction cells lose less than 60% of CO2 gas content over a span
of 18 months. The Radiation Background on Mars is complex and can not fully be recreated
accurately. We therefore focus on simulating electromagnetic radiation as found on the surface of
Mars. To do so we use a Xenon Arc lamp that produces a wide spectrum of light similar to solar
radiation. It also produces significant amounts of UV radiation below 300 nm so that it can be used
as an adequate radiation source for Mars Simulation. The setup has space for up to 10 reaction cells
placed in a ring for the most uniform irradiation. The irradiance was checked at each sample spot
with a relative variation in irradiance of less than 5%. The custom made sample holder can be cooled
using liquid nitrogen (LN2). An off the shelf solenoid valve is used to control the flow of LN2 while
the temperature is controlled using a PT1000 temperature probe in the same form factor as the
reaction cells. In practice this system can be used to cool samples to temperatures between room
temperatures and about −150 ℃. The custom PID control is not limited to a fixed temperature but
also allows to perform temperature protocols (e.g. diurnal cycles). The variance in temperature from
the setpoint using this temperature control is typically below 1 ℃.
FTIR spectroscopy is performed using the ARCoptix OEM FT-IR module which is also planned to be
used in the ExoCube Chem LEO experiment. For UV-VIS measurements we use an Ocean Insight
Flame-S UV-VIS spectrometer, which is planned to be used in the OREOCube LEO experiment. Both
spectroscopy setups are placed on a xy-stage, measuring individual samples in transmission during
the irradiation. The spectroscopic measurements are fully automated, such that only the exchange of
liquid nitrogen has to be performed manually. The setup will be used in the context of the ExoCube
Halo project to investigate photochemical processes involving halophilic organisms exposed to
extreme Mars like conditions.
Acknowledgements
This work was supported by the Deutsche Forschungsgemeinschaft (DFG, grant number 490702919)
and the Volkswagen Foundation and its Freigeist Program.
References
1 https://doi.org/10.1089/ast.2013.1106
2 https://doi.org/10.1007/s11214-021-00812-z
3 https://doi.org/10.1016/j.pss.2016.06.007.
4 https://doi.org/10.1016/j.chemer.2020.125605.
5 https://doi.org/10.1038/s41550-020-1080-9
6 https://doi.org/10.1016/j.pss.2017.01.014.
7 https://doi.org/10.1029/1999JE001095
8 https://doi.org/10.1016/j.icarus.2018.08.019
9 https://doi.org/10.1016/j.actaastro.2012.09.009.
10 https://ui.adsabs.harvard.edu/abs/2022cosp…44.2758W

How to cite: Nitsche, R., Wipf, S., Bourmancé, L., Kish, A., and Elsaesser, A.: Towards Low Earth Orbit Exposure Experiments on the ISS -Designing a Simulation Setup for Mars Like Conditions, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-914, https://doi.org/10.5194/epsc2024-914, 2024.

EPSC2024-982 | ECP | Posters | TP3

Charge Distribution of Ejected Particles after Impact Splash on Mars: A Laboratory Approach 

Tim Becker, Florence Chioma Onyeagusi, Jens Teiser, Teresa Jardiel, Marco Peiteado, Olga Munoz, Julia Martikkainen, Juan Carlos Gomez Martin, Jonathan Merrison, and Gerhard Wurm
Tue, 10 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P4

We present laboratory measurements for the charge distribution of dust to sand sized particles that are ejected in small-scale sand grain impacts, simulating saltation bombardement events on the Martian surface. Using a high speed camera, tracking particles within an electric field, we determined charges of individual grains characterizing single grain impacts. While the charge of small particles seems to have random polarity, larger particles show some bias toward positive charges. Such preference could lesd to a charge separation in the free air stream and thus to the establishment of an electric field near the surface, aiding further lifting as e.g. proposed by Renno&Kok 2008* and Holstein-Rathlou et al. 2010**.

*Renno, N. O., & Kok, J. F. 2008, SSRv, 137, 419, doi: 10.1007/s11214-008-9377-5
**Holstein-Rathlou, C., Gunnlaugsson, H. P., Merrison, J. P., et al. 2010, Journal of Geophysical Research: Planets, 115, doi: 10.1029/2009JE003411

How to cite: Becker, T., Onyeagusi, F. C., Teiser, J., Jardiel, T., Peiteado, M., Munoz, O., Martikkainen, J., Gomez Martin, J. C., Merrison, J., and Wurm, G.: Charge Distribution of Ejected Particles after Impact Splash on Mars: A Laboratory Approach, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-982, https://doi.org/10.5194/epsc2024-982, 2024.

TP4 | Planetary Dynamics: Shape, Gravity, Orbit, Tides, and Rotation from Observations and Models

EPSC2024-952 | ECP | Posters | TP4 | OPC: evaluations required

Satellite gravity-rate observations to uncover Martian plume-lithosphere dynamics 

Riva Alkahal and Bart Root
Wed, 11 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P26

In the past few decades, Mars-oriented orbiters and landers have allowed to unravel valuable knowledge about Mars’ surface and interior. With the InSight mission, seismic waves have indicated the presence of more frequent Marsquakes than assumed before the mission (Banerdt et al. 2020). Moreover, active mantle plume is considered below the Elysium Region (Broquet and Andrews-Hanna, 2023). This raises questions regarding the planet's formation and whether Mars is more geologically active than was considered.

An important milestone in studying the interior of Mars is the recovery of static gravity field models. These models have been accomplished using data from the three recent Mars orbiting missions, namely, Mars Global Surveyor (MGS), Mars Odyssey (ODY), and the Mars Reconnaissance Orbiter (MRO). In addition to the static gravity field, seasonal variations of Mars’ gravity field have been observed, providing information regarding the periodic behavior of the polar ice caps (Konopliv et al. 2016, Genova et al. 2016). However, the secular variation of the gravity field and its link to the solid deformation of the planet has been limited studied.

In general, the estimation of the time variations of the gravity field in the very long wavelength can provide insights into activity of the mantle (Wörner et al. 2023). Le Maistre et al., (2023) have studied the spin rate of Mars and its connection to interior mantle flow or atmospheric changes. By analyzing measurements from the Viking and InSight landers, they estimated a long-term change of the rotation rate of Mars and its moment of inertia.  The obtained rotation rate change, along with the J2 coefficient variation over one Martian year, suggests factors such as atmospheric changes, glacial rebound of the polar ice caps (GIA, Glacial Isostatic Adjustment), or substantial deep mantle flow. Therefore, decoupling atmospheric signal from solid Mars deformations in the gravity signal is essential.  

In this study, we focus on a new way of estimating secular variations of the gravity field of Mars from the available tracking data with an open-source orbit estimation tool: TUDAT (TU Delft Astrodynamics Toolbox). First, we review the state-of-the-art literature on studying the plume-lithosphere interaction and model the gravity-rate signal that would come from mantle flow. Then, we perform a sensitivity analysis for decoupling the secular variations from other signals, such as, the atmospheric density variations and ongoing GIA of the polar ice caps. We do this by simulating one-way and two-way Doppler observations of a Mars-orbiting satellite. We include all possible dynamic forces impacting the satellite. Some of these forces are the static and temporal gravity field, the third body gravitation, the solar radiation pressure, the atmospheric drag, and other forces. For the atmospheric drag, we use the Mars-DTM atmosphere density model that models the static, daily, and yearly variations that affect the drag of the satellite (Bruinsma and Lemoine, 2002). We determined the sensitivity of the estimation process for different parameters including: the initial state, the atmospheric drag and the solar radiation coefficients, and the global vs. arc-wise time varying coefficients. Finally, we perform a correlation analysis of these parameters to determine in which estimation scenario we are able to separate the atmospheric signal from the solid Mars gravity changes. This sensitivity analysis will help in decoupling the gravity-rate signal in order to answer the unresolved question about the activity of the Martian interior.

References:

Banerdt, W.B., Smrekar, S.E., Banfield, D. et al. Initial results from the InSight mission on Mars. Nat. Geosci. 13, 183–189 (2020). https://doi.org/10.1038/s41561-020-0544-y

Broquet, A., Andrews-Hanna, J.C. Geophysical evidence for an active mantle plume underneath Elysium Planitia on Mars. Nat Astron 7, 160–169 (2023). https://doi.org/10.1038/s41550-022-01836-3

Bruinsma, S., & Lemoine, F. G. (2002). A preliminary semiempirical thermosphere499
model of mars: Dtm-mars. Journal of Geophysical Research: Planets, 107 (E10). doi: https://doi.org/10.1029/2001JE001508

Genova, A., Goossens, S., Lemoine, F.G., Mazarico, E., Neumann, G.A., Smith, D.E., Zuber, M.T., 2016. Seasonal and static gravity field of Mars from MGS, Mars Odyssey and MRO radio science. Icarus 272, 228–245. http://dx.doi.org/10.1016/j.icarus.2016.02.050

Konopliv, A.S., Park, R.S., Folkner, W.M., 2016. An improved JPL Mars gravity field and orientation from Mars orbiter and lander tracking data. Icarus 274, 253–260. http://dx.doi.org/10.1016/j.icarus.2016.02.052.

Le Maistre, S., Rivoldini, A., Caldiero, A. et al. Spin state and deep interior structure of Mars from InSight radio tracking. Nature 619, 733–737 (2023). https://doi.org/10.1038/s41586-023-06150-0

Wörner, L., Root, B. C., Bouyer, P., Braxmaier, C., Dirkx, D., Encarnação, J., Hauber, E., Hussmann, H., Karatekin, Ö., Koch, A., Kumanchik, L., Migliaccio, F., Reguzzoni, M., Ritter, B., Schilling, M., Schubert, C., Thieulot, C., Klitzing, W. v., & Witasse, O. (2023). MaQuIs—Concept for a Mars Quantum Gravity Mission. Planetary and Space Science, 239, 105800. https://doi.org/10.1016/j.pss.2023.105800

How to cite: Alkahal, R. and Root, B.: Satellite gravity-rate observations to uncover Martian plume-lithosphere dynamics, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-952, https://doi.org/10.5194/epsc2024-952, 2024.

TP5 | Ionospheres of unmagnetized or weakly magnetized bodies

EPSC2024-34 | Posters | TP5 | OPC: evaluations required

Gravity Wave-Induced Ionospheric Irregularities in the Martian Atmosphere 

Rong Tian, Chunhua Jiang, and Beatriz Sánchez‐Cano
Tue, 10 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P23
For the past few decades, it has demonstrated that gravity waves (GWs) and neutral winds can drive ionospheric irregularities on Earth. Still, as far as we know, the formation of ionospheric irregularity on Mars due to GWs has not been well studied. In this study, we use data from NASA's Mars Atmosphere and Volatile Evolution (MAVEN) mission to show evidence of an irregularity event in the Martian ionosphere, which is potentially seeded by the break of GWs (GWB). The statistical findings indicate that the observed ratio of GWB-related irregularity events varies from ~0.25 to 0.57 in each year, and the overall correlation for 2015 to 2020 is ~0.37. Numerical simulations provide further insight into the processes behind irregularities formation, which employs neutral wind shear as a source of perturbation in the context of the GWB. The simulations yield results fundamentally aligned with the observed characteristics of ionospheric irregularities observed in the 2018 event by considering the wind shear as the disturbance source. This study provides supplementary insights into the perturbation sources involved in shaping irregularities within the Martian ionosphere and presents valuable information about the coupling between the Martian ionosphere and the lower atmosphere.

How to cite: Tian, R., Jiang, C., and Sánchez‐Cano, B.: Gravity Wave-Induced Ionospheric Irregularities in the Martian Atmosphere, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-34, https://doi.org/10.5194/epsc2024-34, 2024.

EPSC2024-219 | ECP | Posters | TP5

Evolution of the ion dynamics at comet 67P during the escort phase 

Zoe Lewis, Peter Stephenson, Esa Kallio, Marina Galand, and Arnaud Beth
Tue, 10 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P27

Comet 67P/Churyumov-Gerasimenko was escorted by the Rosetta spacecraft through a 2 year section of its 6 year orbit around the Sun. This enabled the observation of a large variation in comet outgassing and the resulting evolution of the plasma environment. The diamagnetic cavity, a region of negligible magnetic field arising from the interaction of the unmagnetised cometary plasma with the solar wind, began to be detected sporadically by the Rosetta Plasma Consortium/ Magnetometer (RPC/MAG) in April 2015 at a heliocentric distance of 1.8 au [1]. The last detections were in February 2016 at 2.4 au. Within this cavity, the flow of cometary ions has been shown to be largely radial [2]; the ions are accelerated above the neutral gas speed by an ambipolar electric field, but many newborn ions still undergo multiple ion-neutral chemical reactions before escaping [3,4]. Outside the diamagnetic cavity boundary, which is itself highly variable, the ion flow is considerably more complex, and the ambipolar electric field plays a more minor role compared to the convective electric field of the solar wind [2].  At large heliocentric distances (>2.5 au), the total plasma density observed from RPC plasma sensors is well explained by a simple flux conservation model that assumes the ions travel radially away from the nucleus at speed close to that of neutrals [5,6]. However, closer to perihelion and once the diamagnetic cavity has formed, such an approach does not hold [7]. We aim to better understand this transition, the driver of ions' acceleration, and the role that the diamagnetic cavity plays.In this study, we explore the varying ion dynamics both in the presence (e.g. during high outgassing activity) and absence (low outgassing activity) of a diamagnetic cavity. Electric and magnetic fields from hybrid simulations of the cometary environment are used to drive a 3D test particle model of the cometary ions for a range of comet activity levels. We model the behaviour of three key ion species, H2O+, H3O+, and NH4+, in order to assess the impact of the ion dynamics on the ionospheric composition and density.

 [1] Goetz et al. MNRAS S 462, S459–S467 (2016)

[2] Koenders et al., Planetary and Space Science, 101-116, 105 (2015)

[3] Lewis et al., MNRAS, 523, 6208–6219 (2023)

[4] Lewis et al 2024, MNRAS, 530, 66–81 (2024)

[5] Galand et al., MNRAS, S331-S351, 462 (2016)

[6] Heritier et al., A&A, 618 (2018)

[7] Vigren et al., ApJ 6, 881(1) (2019)

How to cite: Lewis, Z., Stephenson, P., Kallio, E., Galand, M., and Beth, A.: Evolution of the ion dynamics at comet 67P during the escort phase, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-219, https://doi.org/10.5194/epsc2024-219, 2024.

EPSC2024-1287 | ECP | Posters | TP5 | OPC: evaluations required

Preliminary Results of the Categorizing and Statistical Survey of Martian Ionospheric Irregularities 

Xin Wan, Jiahao Zhong, Chao Xiong, and Jun Cui
Tue, 10 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P26

A climatological survey of Martian ionospheric plasma density irregularities was conducted by exploring the in-situ measurements of the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft. The irregularities were first classified as enhancement, depletion, and oscillation. By checking the simultaneous magnetic field fluctuation, the irregularities have been classified into two types: with or without magnetic signatures. The classified irregularities exhibit diverse global occurrence patterns, as those with magnetic signatures tend to appear near the periphery of the crustal magnetic anomaly (MA), and those without magnetic signatures prefer to appear either inside of the MA or outside of the MA, depending on the type and solar zenith angle. Under most circumstances, the irregularities have a considerable occurrence rate at altitudes above the ionospheric dynamo height (above 200 km), and the magnetization state of the ions seems irrelevant to their occurrence. In addition, the irregularities do not show dependence on magnetic field geometry, except that the enhancement without magnetic signatures favors the vertical field line, implying its equivalence to the localized bulge (Duru et al., 2006). Other similarities and discrepancies exist in reference to previous studies. We believe this global survey complements previous research and provides crucial research clues for future efforts to clarify the nature of the Martian ionospheric irregularities.

How to cite: Wan, X., Zhong, J., Xiong, C., and Cui, J.: Preliminary Results of the Categorizing and Statistical Survey of Martian Ionospheric Irregularities, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1287, https://doi.org/10.5194/epsc2024-1287, 2024.

TP6 | Late accretion and early differentiation of rocky planetary bodies, from planetesimals to super-Earths

EPSC2024-763 | ECP | Posters | TP6

Resolving the origin of lunar high-Ti basalts by petrologic experiments 

Cordula Haupt, Christian J. Renggli, Arno Rohrbach, Jasper Berndt, and Stephan Klemme
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P35

The origin of the most primitive, picritic lunar basalts, sampled as pyroclastic glass beads in the lunar soils [1,2], remains poorly constrained. Especially the petrogenesis of high-Ti glasses (TiO2 > 6 wt%) seems enigmatic. Three hypotheses have been proposed for the origin of these exotic samples: I) Ascent of a primary, undifferentiated melt from the lunar interior and assimilation of clinopyroxene and ilmenite in the upper section of the not overturned lunar mantle [3]. II) Melting of a hybridized lunar mantle after lunar mantel overturn [4,5]. III) Reaction of sunken low-degree high-density melts of the hybrid magma ocean source with high Mg-cumulates in the deep interior of the lunar mantle and subsequent ascent [5,6,7]. This study re-investigates hypothesis II with the aim to assess whether a one stage melting process of a heterogeneous lunar mantle can cause the compositional variabilities of lunar (high-Ti) picritic glasses. Specifically, the effect of modal mineralogy of different cumulate layers in the hybrid lunar mantle is investigated.

The overturn of the lunar mantle stratification due to Rayleigh-Taylor instabilities will have caused the so-called “Ilmenite-bearing cumulate (IBC)” to sink into the underlying harzburgite lunar mantle [7]. Therefore, an IB cumulate and a harzburgite cumulate appear to be the major components of a hybrid lunar mantle [7]. In this study, the first batch of starting material compositions was mixed similar to [5] using a fixed composition the harzburgite mantle. An IBC was designed by mixing ilmenite, clinopyroxene, and small amounts of plagioclase. These minerals are the basic components of the bulk lunar mantle after cumulate overturn [2]. We investigated how the ilm/cpx ratio within an IBC will affect the melt compositions and melting conditions. A such, we assumed that modal proportions of crystallization were not preserved in the cumulate [e.g., 5,9]. In a second batch of experiments, we slightly adjusted olivine and orthopyroxene ratios in the harzburgite layer. A third component of a hybrid lunar mantle in some of our starting material composition was an urKREEP component [10], which has been proposed to participate in the overturn and melting process [7,11].

To investigate the composition and modal amounts of partial melts from several different starting materials, we conducted high-pressure and high-temperature experiments in an endloaded Piston cylinder apparatus at the University of Münster. Most runs were conducted at a pressure of 1.5 GPa, which corresponds to the lower end of the depth range suggested for the source depth of high-Ti lunar picritic basalts [6]. Run temperatures were varied between 1300 and 1450 °C to investigate the effect of changing melting degree on melt compositions [3]. In order to control fO2 and to minimize Fe-loss in the runs, we used graphite-lined Pt capsules [3,4]. The characterization of experimental runs was conducted by the means of scanning electron microscopy (SEM) and energy dispersive x-ray spectroscopy (EDS). Mineral and melt proportions were determined via mass balance calculations using the major element chemistry of the phases present in each experiment.

All experiments contain glass and forsteritic olivine. Some experiments contain olivine and orthopyroxene. The degree of partial melting (Fmelt) is 0.18–0.75. The experiments that contain both olivine and orthopyroxene were run at lower temperatures (< 1400 °C) and low degrees of partial melting (Fmelt < 0.5). 

 

Figure 1: SiO2 vs. TiO2 of lunar Apollo picritic glasses (yellow, orange, red, black – colored accordingly) from [8,9] and our experimental melts (squares and circles). Indivudual symbols correspond to different starting material compositions.Triangles correspond to melting experiments of a hybrid lunar mantle with modal ilm/cpx from fractional crystallization experiments [5] 

Comparing the composition of experimental melts and natural lunar picritic glasses (Fig. 1), it can be stated that the melting of a heterogeneous lunar mantle produced by the overturn of lunar stratification after the solidification of the lunar magma ocean can generate melts in the range similar high-Ti picritic melts. Experimental temperatures and pressures agree with the temperatures and depth of origin predicted by previous experimental studies [3,5]. A partial melting process of a cumulate-bearing mantle, as modeled by our experiments, is shown to be a viable and simple alternative to the currently accepted complex melting model [5,6]. In our experimental setup, a ratio of ilm/cpx of 1/1 or 4/1 in the IBC layer reproduces high-Ti compositions, similar to the picritic lunar basalts. Good matches are achieved in runs with lower temperatures, which correspond to the comparably lower degree of melting.

The most suitable ol/opx is 3/2. We further find that, following the constraints of [12], some plagioclase has to be entrained in the IBC layer, in order to reproduce Al2O3/CaO in the cumulate mantle melting experiments. Additionally, the presence of urKREEP in the cumulates strongly influences Al2O3/CaO, driving it too high in melts originating from cumulates containing that component.

In the light of our experiments, it is possible to shed some new light on the origin of exotic lunar basalt samples, such as the picritic, high-Ti lunar basalts. We explored the feasibility of a simple melting process of a hybrid lunar mantle after overturn. 

[1] Delano (1986) J Geophys Res-Solid 91, B4 201-213 [2] Shearer C. K. and Papike J. J. (1993) Geochim Cosmochim Ac, 57, 19, 4785-4812 [3] Wagner T. P. and Grove T. L. (1997) Geochim Cosmochim Ac, 61, 6, 1315–1327. [4] Krawczynski M. J. and Grove T. L. (2012) Geochim Cosmochim Ac, 79, 1-19. [5] Singletary G. H. and Grove T. L. (2008) Earth Planet Sc Lett, 208, 182–189. [6] Elkins-Tanton L. T., van Orman J. A. et al. (2002) Earth Planet Sc Lett, 196, 239-249. [7] Hess P. C. and Parmentier E. M. (1995) Earth Planet Sc Lett 134, 501-514 [8] Elkins L. T., Fernandes V. A. et al. (2000) Geochim Cosmochim Ac, 64.13, 2339–2350. [9] Brown, S. M., Grove, T. L. (2015). Geochim Cosmochim Ac171, 201-215. [10] Warren, P. H., Wasson, J. T. (1979). Rev Geophys17(1), 73-88. [11] Elardo, S. M., Laneuville, M. et. al (2020). Nat Geosci13(5), 339-343. [12] Charlier, B., Grove, T. L., et al. (2018). Geochim Cosmochim Ac234, 50-69. 

How to cite: Haupt, C., Renggli, C. J., Rohrbach, A., Berndt, J., and Klemme, S.: Resolving the origin of lunar high-Ti basalts by petrologic experiments, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-763, https://doi.org/10.5194/epsc2024-763, 2024.

EPSC2024-1006 | ECP | Posters | TP6 | OPC: evaluations required

1D Modeling of the Magma Ocean Stage of Rocky Planets  

Meiye Wu and Lena Noack
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P31

To better understand the evolution of rocky planet interiors and their redox transformations, we are developing a comprehensive 1D mineralogical and geochemical interior model. This model is designed to simulate the initial conditions and subsequent evolution of rocky planet interiors in the magma ocean stage.

This project creates a 1D grid-based compressible interior-structure model with a magma ocean thermal evolution and solidification model, building upon existing codes and models [1]. Differing from previous magma ocean models, this 1D model incorporates the evolution of core differentiation and couples it with a magma ocean model. We also further enhance previous models by incorporating depth-dependent thermodynamic properties and implementing high-temperature and high-pressure melt Equations of State (EOS). We are developing a composition-dependent melting temperature formulation that aligns with low-pressure melting temperature laws. For the high-pressure environments of planetary interiors, we apply the extended Lindemann-Stacy melting law [2]. Moreover, the early planet accretion stage is considered in this model and different accretion scenarios are investigated.

Our approach is expected to provide thermal and chemical profiles from planet formation till the end of the magma ocean stage. Our results will help determine the efficiency of core formation under different redox states and planetary conditions during the magma ocean stage. In the future, the results of this model could also be used as inputs for 2D convection simulations or planetary atmospheric models.

 

References

  • [1]  L. Noack, D. Höning, A. Rivoldini, C. Heistracher, N. Zimov, B. Journaux, H. Lammer, T. Van Hoolst, and J.H. Bredehöft. Water-rich planets: How habitable is a water layer deeper than on Earth? Icarus, 277:215–236, 2016.

  • [2]  V. Stamenković, D. Breuer, and T. Spohn. Thermal and transport properties of mantle rock at high pressure: Applications to super-earths. Icarus, 216(2):572– 596, 2011.

How to cite: Wu, M. and Noack, L.: 1D Modeling of the Magma Ocean Stage of Rocky Planets , Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1006, https://doi.org/10.5194/epsc2024-1006, 2024.

TP7 | Planetary volcanism, tectonics, and seismicity

EPSC2024-650 | ECP | Posters | TP7 | OPC: evaluations required

Effects of magmatic styles on the thermal evolution of planetary interiors 

Carianna Herrera, Ana-Catalina Plesa, Julia Maia, and Doris Breuer
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P44

ABSTRACT

Previous studies have suggested that extrusive magmatism efficiently cools planetary interiors but the contribution by intrusive magmatism has been little investigated so far. We study the effects of the magmatic style (i.e., intrusive vs. extrusive magmatism) on the thermal evolution of Mercury-, Venus-, Moon-, and Mars-like bodies. Our results suggest that different magmatic styles strongly affect the thermal evolution of planets, where for instance intrusive magmatism allows for thinner lids, cooler melts and more melt production. In addition, for large and/or highly magmatically active planets such as Venus, an intrusive magmatism allows for more efficient cooling of the interior.

INTRODUCTION

Volcanic activity has been evidenced on planetary bodies across the inner solar system [1]. Mercury's volcanic surface expressions are extrusive volcanic vents, pyroclastic deposits, lava flow margins, etc.; and studies have analyzed their history in terms of the effect of the planetary cooling who led to the global contraction [2,3]. Venus surface is dominated by volcanoes, flow fields, volcanic vents, and coronae, among other features [4]. Recent analysis of radar data collected by the Magellan mission suggests that there could be volcanic activity still ongoing on Venus [5].

Volcanism has also shaped the lunar surface, leading to the formation of a variety of volcanic landforms such as lunar maria, lava flows, domes, cones, etc., features that were confirmed after a series of lunar exploration missions and even directly studied thanks to the return of samples from the lunar surface [6,7]. The extensive volcanic products on Mars account a large variety of explosive [8] and sedimentary volcanism [9], and recent studies presented geophysical evidence for active volcanic processes in the Elysium Planitia region [10,11].

These volcanic features are witnesses of magmatic processes that these bodies have experienced, but the amount of intrusive vs. extrusive melt that was produced during the thermal history is difficult to constrain. Extrusive magmatism has been suggested to allow more efficiently cooling of planetary interiors, but the role of intrusive magmatism on the thermal evolution has been little investigated.

In this study, we analyze the effect of the magmatic style (i.e. ‘fully intrusive’ vs. ‘fully extrusive’ magmatism end-members) by modelling the thermal evolution of Mercury, Venus, Moon, and Mars-like bodies. Our aim is to understand the effect of different magmatic styles for different planetary interiors rather than the particular evolution of the inner solar system bodies.

METHODS

We use the mantle convection code GAIA in a 2D spherical annulus geometry [12, 13]. Our models employ a temperature- and pressure-dependent viscosity that follows an Arrhenius law for diffusion creep [14, 15]. The strong temperature-dependence of the viscosity leads to the formation of a stagnant lid (an immobile layer) at the top of the convecting mantle, due to the cold temperature conditions. The thermal conductivity and thermal expansivity in our models are pressure- and temperature-dependent and we use parametrizations derived from ab-initio calculations and laboratory experiments [16]. We assume a homogeneous distribution of the heat sources and account for the decay in time of radioactive elements (i.e. 238U, 235U, 232Th, and 40K) and consider that the core cools with time.

Melting occurs when the mantle temperature exceeds the melting temperature. To keep the models as best as possible comparable with each other, we use the same melting curve parametrization as derived for the Earth’s interior [17]. We compute partial melting and consider two scenarios, i.e., fully intrusive and fully extrusive magmatism (Fig. 1). For all bodies, the depth of magmatic intrusions is set at 50 km depth. For scenarios where partial melting occurs deeper than the density crossover at ~11GPa [18], the melt is not buoyant enough to rise towards the surface and it is thus not considered in our models.

RESULTS

For all studied bodies, the convection pattern is characterized by stronger mantle plumes and more vigorous mantle flow for the fully intrusive cases than for the fully extrusive cases (Fig. 2). Additionally, the intrusive melt depth seems to control the stagnant lid growth with thinner lids for cases with magmatic intrusions.

While the global average temperature of the entire silicate part is higher for Mercury, Mars, and the Moon in the intrusive scenario compared to extrusive models, for Venus the opposite is the case (Fig. 3a). We explain this by Venus’s smaller lid-to-mantle thickness ratio and high melt production rate compared to Mercury, Mars, and the Moon. Normalized to their silicate mantles, the intrusive magmatism models produce more melt than the extrusive cases (Fig. 3d).

The mechanical lid depth of the intrusive cases is always shallower through time (Fig. 3b), diverging from the extrusive cases up to hundreds of kilometers. This is explained by higher thermal gradients and thus warmer lithospheres in the intrusive cases.

Mercury cools very quickly compared to the other planets and stops producing melt after the first half of its evolution (Fig. 3c). After that moment, the differences between the two scenarios are nearly identical. For all bodies, the intrusive magmatism cases melt at shallower depths with cooler melts during their evolution.

SUMMARY       

Throughout the evolution of all studied bodies, the fully intrusive cases present thinner mechanical lids, cooler melt temperatures, more melt production, and shallower melting depths than the extrusive cases. Our results suggest that large and/or highly magmatically active planets such as Venus efficiently cool their interior through intrusive magmatism, while keeping at the same time a warm and thin lithosphere.

REFERENCES

[1] Byrne et al., Nat.Astron, 2020; [2] Thomas & Rothery, Elements, 2019; [3] Wright et al., J. Volcanol., 2021; [4] Ghail et al., Space Science Reviews, 2024; [5] Herrick & Hensley, Science, 2023; [6] Head, Rev Geophys., 1976; [7] Zhao et al, Space: Science & Technology, 2023; [8] Brož et al., J. Volcano., 2021; [9] Brož et al., Earth Surf. Dynam., 2023; [10] Broquet & Andrews-Hanna, Nat.Astron., 2023; [11] Stähler et al. Nat.Astron., 2022, [12] Hüttig et al., PEPI, 2013; [13] Fleury et al., Geochem.Geophys., 2024; [14] Karato et al., JGR, 1986; [15] Karato & Wu, 1993; [16] Tosi et al., PEPI, 2013; [17] Stixrude et al., EPSL, 2009; [18] Ohtani et al., Chem. Geol., 1995.

How to cite: Herrera, C., Plesa, A.-C., Maia, J., and Breuer, D.: Effects of magmatic styles on the thermal evolution of planetary interiors, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-650, https://doi.org/10.5194/epsc2024-650, 2024.

EPSC2024-1106 | ECP | Posters | TP7

Tectonic influence of multi-ring basins: The case of Mercury’s Discovery Quadrangle and the Andal-Coleridge basin. 

Antonio Sepe, Luigi Ferranti, Valentina Galluzzi, Gene Walter Schmidt, Salvatore Buoninfante, and Pasquale Palumbo
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P47

Introduction
Mercury’s Discovery quadrangle is located at southern mid-latitudes in a heavily cratered region roughly antipodal to Caloris Basin [1]. The quadrangle hosts a probable pre-Tolstojan multi-ring impact basin named Andal-Coleridge, surrounded by a three- to five-ring system [2,3].
Here we present a high-resolution structural map of the quadrangle that will contribute to the 1:3M quadrangle geological map series in preparation for the BepiColombo mission [4]. In addition, we carried out a structural analysis in order to study the structural framework of the quadrangle, validate and ascertain the existence of the Andal-Coleridge basin and investigate its influence on the tectonic evolution of this region.

Data and methods
We produced a high-resolution structural map of the quadrangle utilizing MESSENGER end-of-mission products [5]. Structure strikes were then plotted in rose diagrams to recognize preferential trends at the regional scale. We also investigated the structural relationships between three main scarps (Discovery, Adventure and Resolution Rupes – DAR) –including another scarp (here named Discovery-2) that seems to be the westward continuation of the Discovery Rupes (Discovery-1)– by making several profiles across them and measuring their throw. Furthermore, to validate and ascertain the presence of the Andal–Coleridge basin, we employed gravimetric data and estimated the related stress field via beta-analysis. 

Map and trend of structural features
Our structural map (Fig.1) reveals ∼500 segments of contractional structures –including lobate scarps, high-relief ridges and wrinkle ridges– mainly arranged in a circular pattern at the approximate centre of the quadrangle, encircling both a broad topographic low and a mascon, similar to Caloris Basin [6], although identified only in the Bouguer anomaly map of [7]. Strike directions of the segments within rose diagrams generally show a uniform distribution of bins –together with a NW-SE preferential trend– suggesting an impact-related nature for most structures [8].

Figure 1 Structural map of the Discovery quadrangle overlayed to the Digital Elevation Model of [9] (left) and the Bouguer Anomaly Map of [7] (right). The red dashed line approximately delimits the topographic low and the mascon respectively.

Beta-Analysis
Following the work done for Mercury and Mars by [2] and [10], we estimated the stress field related to concentric structures via beta-analysis. This approach reveals a bimodal distribution consisting of a primary "bull's-eye" distribution, whose centre represents the point from which the causative stress field for faults spread, and a secondary, much less dense, distribution (Fig.2). Discovery Rupes aligns concentrically with the primary distribution, while Adventure and Resolution Rupes exhibits concentric alignment with the secondary distribution roughly coincident with another probable impact basin, informally named b78 [11], where another smaller mascon is present.

Figure 2 (left) The stress field resulting from beta-analysis in stereographic projection centred at quadrangle’s centre. (right) The hypothesised location and extent of Andal-Coleridge and b78 basins from [11].

Throw-Height Analysis
For each profile the scarp height was measured and then plotted against its position on the fault trace, corresponding to the fault length, resulting in a throw profile (Fig.3). The analysis of the throw profiles provide evidence that the three structures represent the morphological expression of two different faults, the Discovery fault and the Adventure-Resolution fault [2], grown hard-linking several segments together. These two faults also appear to be kinematically soft-linked, since the cumulative throw falls approximately at the centre of the system, consistent with terrestrial fault growth patterns [12]. Additionally, Discovery Rupes shows evidence of reactivation within Rameau crater, suggesting a process of fault-trace erosion and its subsequent re-emergence possibly due to Mercury's global contraction. Adopting previous dating of Rameau crater and Discovery Rupes [14,15], we derived a preliminary chronology for Discovery Rupes evolution and found a peak in throw-rate during the Tolstojan period, with declining activity towards the Mansurian period which is possibly continuing to the present [15].

Figure 3 (up) The DAR system and the 40 profiles used to produce the throw-height plot (down). Each peak on the plot represents a fault segment.

Conclusions and future work
The present work reveals a complex structural framework within Mercury’s Discovery quadrangle. The concentric alignment of structures to both a topographic low and a mascon together with the general uniform distribution shown in rose diagrams strongly supports the existence of the Andal-Coleridge basin. The concentricity of Discovery Rupes to the primary distribution identified through beta-analysis emphasizes its connection to the basin. In contrast, Adventure and Resolution Rupes display concentric alignment with the smaller b78 impact basin. This observation implies that these two scarps originated from two distinct impacts and were reactivated by Mercury's global contraction as a linked fault system. The Andal-Coleridge basin likely introduced mechanical discontinuities in the crust, influencing the localization and orientation of faults [2]. These weaker zones were eventually exploited during the subsequent global contraction. Indeed, the throw-height analysis shows signs of reactivation where Discovery Rupes cuts Rameau crater, suggesting a peak in throw-rate during the Tolstojan period.
We aim at improving the structural map of the quadrangle to better understand the NW-SE-trending structures’ behaviour and better investigate the fault reactivation. This study may also represent a contribution for evaluating the rate of global contraction and its magnitude throughout Mercury’s evolution.

Acknowledgements
We gratefully acknowledge funding from the Italian Space Agency (ASI) under ASI-INAF agreement 2017-47-H1.

References
[1] Trask and Dzurisin (1984). USGS, IMAP 1658. [2] Watters et al. (2001). Plan. and Sp. Sci., 49(14-15), 1523-1530. [3] Spudis and Strobell (1984). LPSC, 814-815. [4] Galluzzi et al. (2019). JGR: Planets, 124(10), 2543-2562. [5] Denevi et al. (2017). Sp. Sci. Rev., 214(1). [6] Smith et al. (2012). Science, 336. [7] Buoninfante et al. (2023). Sci. Rep., 13, 19854. [8] Delbo et al. (2019). LPI Contrib. No. 2189. [9] Becker et al. (2016). LPSC Contrib. No. 1903. [10] Wise et al. (1979). Icarus, 38, 456-472. [11] Orgel et al. (2020). JGR: Planets, 125(8). [12] Kim and Sanderson (2005). Earth-Sci. Rev., 68(3-4), 317-334. [13] Kinczyk et al. (2020). Icarus, 341, 113637. [14] Clark et al. (2024). LPSC Contrib. No. 3040. [15] Tosi et al. (2013). JGR: Planets, 118(12), 2474-2487.

 

How to cite: Sepe, A., Ferranti, L., Galluzzi, V., Schmidt, G. W., Buoninfante, S., and Palumbo, P.: Tectonic influence of multi-ring basins: The case of Mercury’s Discovery Quadrangle and the Andal-Coleridge basin., Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1106, https://doi.org/10.5194/epsc2024-1106, 2024.

TP9 | Impact Processes in the Solar System

EPSC2024-372 | ECP | Posters | TP9 | OPC: evaluations required

Evaluation of large craters and basins for lunar production functions 

Astrid Oetting, Wajiha Iqbal, Thomas Heyer, Nico Schmedemann, James Head, Gregory Michael, Harald Hiesinger, and Carolyn van der Bogert
Thu, 12 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P9

Introduction: Large craters and basins that are present on the lunar surface are remants of its early impact history. The relative ages of geological surfaces can be determined with crater size-frequency distribution (CSFD) measurements using a production function (PF), since the impact record on the lunar surface is related to time. Widely used PFs were developed by Neukum (1983) and Neukum et al. (2001) [1,2] and are valid for crater diameters between 0.01–300 km. However, to understand the earlier history of the Moon, when larger impacts were more abundant, an extension of the valid crater diameter range to larger diameters would be beneficial. It could provide input for understanding the number of impactor populations [e.g., 1-4] and the stability of the impact rate on the Moon [e.g., 1,2,5,6]. However, the precise determination of the main basin rim diameter, especially for multi-ring basins, is challenging due to their complex and degraded morphology.

 

Method: The diameters of large craters and basins were measured based on topographic (~100 m/pixel) and gravity data. We used the Lunar Reconnaissance Orbiter (LRO) Wide Angle Camera (WAC) image mosaic [7], LRO Lunar Orbiter Laser Altimeter Digital Elevation Model (LOLA DEM) [8], and the WAC DEM color-shaded relief map [9,10]. In addition, the following geophysical data were used: Gravity Recovery and Interior Laboratory (GRAIL) [11] data, a crustal thickness map (Crustal Thickness – Model 1) [12], and a Bouguer anomaly map [11].

In order to have a consistent measurement of the main basin rim diameter, we follow the approach of Neumann et al. (2015) [13]. They [13] compared the Bouguer anomaly, which reflects changes in the subsurface and/or crustal thickness, with well-preserved basin structures. They suggest that double the diameter of the Bouguer anomaly is consistent with the main or reference rim for multi-ring basins.

The ArcGIS CraterTools add-in [14] was used to perform the CSFD measurements, which were then further analyzed with Craterstats2 [15]. The count area is the entire Moon, but due to significant differences between mare and highlands areas, we subdivide the count area "Entire Moon" into "Highlands" and "Mare".

 

Results & Discussion: We identified 311 craters and basins with diameters between 100-1250 km over the entire Moon. We do not include the South Pole Aitken (SPA) basin as it is unclear which topographic ring is its reference rim and it is assumed to have formed in an even earlier era of lunar basins [16]. The CSFD measurements (Figure 1) of the entire Moon and the highlands are consistent up to diameters of 250 km. However, while the highlands distribution progresses relatively smoothly towards large basin diameters, the CSFD of the entire Moon shows a "step" between 400-700 km. The CSFD of the mare areas is significantly different from the two other areas. Craters <200 km are less abundant than on the other two areas, and craters and basins >250 km are more abundant. The CSFD measurements for the entire Moon and the lunar highlands are better represented by the PF of [2] using the valid crater diameter range of 100-300 km. The mare areas could not be fitted with either PF.

Figure 1: Cumulative CSFD plots for large craters and basins compared with existing PFs (black and blue curves) (a) [1] and (b) [2] in the valid crater diameter range of 100-300 km. The gray line represents the equilibrium function of [17].

 

The significant difference between the CSFDs on highlands and mare areas may result from:

(1) lava flooding [18,19], causing an incomplete identification of basins on the mare areas

(2) asymmetries in cratering rate, e.g. due to the rotation and inclination of the Moon [20,21]

(3) larger basin diameters resultant on the lunar nearside due to a different subsurface temperature and crustal thickness [16]

 

Conclusion: The CSFD measurements indicate that large craters and basins are still in production and, therefore, an extension of PFs for craters and basins from 300-1250 km could be possible. This extented crater diameter range would allow to draw conclusions about the number of impactor populations and the stability of the impact rate in the early lunar history. However, the influence of resurfacing processes on mare units compared to the highlands is not yet entirely understood and a more detailed analysis of the mare regions is necessary.

Acknowledgments: This Project is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 263649064 – TRR 170.

References: [1] Neukum G. (1983) Habil. Thesis LMU, Munich. [2] Neukum G. et al. (2001) Space Sci. Rev., 96, 55. [3] Strom R.G. et al. (2005) Science, 309(5742), 1847-1850. [4] Head J.W. et al. (2010) Science, 329(5998), 1504-1507. [5] Guinness E.A. and Arvidson R.E. (1977) Proc. Lunar Sci Conf., 8, 3475-3494. [6] Bottke W.F. et al. (2005) Icarus, 179(1), 63-94. [7] Robinson M.S. et al. (2010) Space Sci. Rev. 150, 81-124. [8] Smith D.E. et al. (2010) Space Sci. Rev., 150, 209-241. [9] Scholten F. et al. (2012) JGR, 117(E12). [10] Smith D.E. et al. (2010) GRL, 37(18). [11] Zuber M.T. et al. (2013) Science, 339(6120), 668-671. [12] Wieczorek M.A. et al. (2013) Science, 339, 671-675. [13] Neumann G.A. et al. (2015) Sci. Adv., 1(9), e1500852. [14] Kneissl T. et al. (2011) PSS, 59(11-12), 1243-1254. [15] Michael G. et al. (2016) Icarus, 277, 279-285. [16] Miljković K. et al. (2013) Science, 342(6159), 724-726. [17] Hartmann W.K. (1984) Icarus, 60(1), 56-74. [18] Evans A.J. (2016) Geophys. Res. Lett., 43, 2445-2455. [19] Evans A.J. (2018) JGR, 123, 1596-1617. [20] Wiesel W. (1973) Icarus, 15(3), 373-383. [21] Le Feuvre M. and Wieczorek M.A. (2011) Icarus, 214(1), 1-20.

How to cite: Oetting, A., Iqbal, W., Heyer, T., Schmedemann, N., Head, J., Michael, G., Hiesinger, H., and van der Bogert, C.: Evaluation of large craters and basins for lunar production functions, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-372, https://doi.org/10.5194/epsc2024-372, 2024.

EPSC2024-437 | ECP | Posters | TP9 | OPC: evaluations required

Efficacy of classical and spectroscopic techniques for strain quantification in weakly shocked rocks: Results from experimentally impacted Taunus quartzite 

Rajit Das, Amar Agarwal, Arun Ojha, Thomas Kenkmann, and Michael Poelchau
Thu, 12 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P13

Traditional strain estimation techniques, such as the Fry and the Rf/Phi method, have been extensively used in tectonically deformed rocks [1]. However, their applicability in impact rocks has not been explored yet. An essential difference between tectonic and impact cratering-led deformation is that the former causes strain rates in the order of 10-15-10-6 s-1 [2], while the latter leads to a much higher strain rate, ~10-2-106 s-1 [3]. Furthermore, peak pressures in tectonic deformation are in the order of hundreds of megapascals [4], whereas, in the case of naturally shocked rocks, the shock pressures are in the range of gigapascals [5].

To address this, an impact cratering experiment was carried out on a quartzite block to investigate the efficacy of these traditional and spectroscopic tools in determining strain in weakly shocked rocks (Fig. 1). A basalt projectile of diameter 6.18mm was accelerated to a velocity of ~5km/s and impacted a quartzite block, Taunus quartzite from Germany.

The block was sawed through the crater centre. Twelve cores were recovered from the target surface and the sawed (sub-) surface each. Each cylinder was bisected in half to produce 24 specimens, each from the target and subsurface. Specimens 5.1 to 5.24 were from the target surface, whereas the subsurface specimens were labelled 6.1 to 6.24. Specimens 5.1 to 5.12 and 6.1 to 6.12 were from the top half of the cylindrical core obtained from the target surface and subsurface, respectively. Similarly, 5.13 to 5.24 and 6.13 to 6.24 were from the lower half of the cylindrical core obtained from the target surface and the subsurface, respectively.

Figure 1: Photographs of the Taunus quartzite block (edge length 20 cm) (A) before and (B) after the impact experiment. The block was split into several pieces after the experiment. The top/target surface is marked with “Top” and “X” inside a green circle. Note the different orientations of the block in the two images.

Strain analysis was carried out using three methods: (1) Fry, (2) Rf/Phi, and (3) X-ray Diffractometry (XRD) in the impacted Taunus quartzite. These results were then compared to simulated strains in the transient stage from hydrocode modelling (Fig. 2). Radial fractures were observed on the target surface but not in the subsurface in the thin section images (Fig. 3). Results demonstrably suggest that the strain measured using the three methods is similar to the results from the hydrocode simulations. Furthermore, these techniques are sensitive enough to capture strain imparted at shock pressures as low as 0.23 GPa (obtained from hydrocode modelling) at the transient crater stage (at 25 µs post impact).

Figure 2: Results from hydrocode simulation. (A-B) show pressure and density distribution, and (C-D) show damage and total plastic strain at 0 ms and 25 ms. Simulation is run till the end of the transient cratering stage and before the arrival of the reflected wave from the free surface.

Figure 3: Thin section images from the experiment. The arrows indicate the direction of the point of impact. (A) Sample T01 is on the target surface and closest to the point of impact, (B) sample S01 is on the subsurface and just below the point of impact, whereas (C) sample S24 is the point farthest away from the point of impact. Note the radial tensile fractures in the target surface in (A) (parallel to the arrow indicating the direction of the impact point) and their absence in the subsurface (B and C).

However, caution is advised while using Fry and Rf/Phi methods, as grain fracturing and pulverisation may affect the grain centre distribution and grain shape, thus leading to anomalous results. The use of the lattice technique, XRD, can be investigated further to understand the accuracy and efficacy in low-strain conditions. 

Acknowledgements:

We thank the Dept. of Earth Sciences at the Indian Institute of Technology-Kanpur for microscopy and XRD facilities. The authors are also grateful to Dr. Auriol Rae for his input and help on the hydrocode simulation. We are also thankful to Dr. Gaurav Joshi for his suggestions. Yashwant Singh helped in the Fry analysis.

References:

[1] Joshi, G., et al. (2017). Microstructures and strain variation: Evidence of multiple splays in the North Almora Thrust Zone, Kumaun Lesser Himalaya, Uttarakhand, India. Tectonophysics 694, 239–248.

[2] Pfiffner, O. A., et al. (1982). Constraints on geological strain rates: Arguments from finite strain states of naturally deformed rocks, J. Geophys. Res., 87(B1), 311–321.

[3] Ramesh, K. T. (2008). High Rates and Impact Experiments. In: Sharpe, W. (eds) Springer Handbook of Experimental Solid Mechanics. Springer Handbooks. Springer, Boston, MA.

[4] Vaughan-Hammon, J. D., et al. (2021). Alpine peak pressure and tectono-metamorphic history of the Monte Rosa nappe: evidence from the cirque du Véraz, upper Ayas valley, Italy. Swiss J Geosci 114, 20.

[5] Engelhardt, W. V., et al. (1969). Shock induced planar deformation structures in quartz from the Ries crater, Germany. Contr. Mineral. and Petrol. 20, 203–234 (1969).

How to cite: Das, R., Agarwal, A., Ojha, A., Kenkmann, T., and Poelchau, M.: Efficacy of classical and spectroscopic techniques for strain quantification in weakly shocked rocks: Results from experimentally impacted Taunus quartzite, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-437, https://doi.org/10.5194/epsc2024-437, 2024.

EPSC2024-560 | ECP | Posters | TP9 | OPC: evaluations required

Modelling Lunar Impact Flashes from Molten Ejecta 

Paige Rice, Robert Luther, and Kai Wünnemann
Thu, 12 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P6

Introduction:

Lunar impact flashes have been observed from Earth since 1999 [1]. To date, over 600 impact flashes have been detected by individuals as well as multiple observation campaigns [2], including the NASA Lunar Impact Monitoring Program, the MIDAS project, and the NELIOTA project. The record of flashes has been analysed statistically to estimate the current near-Earth meteorite flux and size distribution [3], though the results of such estimates rely on some assumptions about the relationship between the radiative energy of the flashes and impactor velocity and kinetic energy.

Luminous efficiency, the fraction of impactor kinetic energy that is transformed into light emission during the impact process, is a poorly constrained value. Common estimates range over three orders of magnitude, from 10-2 to 10-4 [4]. Recent modelling work produced a detailed model of the radiant process in the impact debris cloud and fit well to the observed flash using calculated impactor properties [5]. Their work builds on previous studies considering the physics of the radiant process, indicating that vaporized material in the instant of impact cannot provide enough thermal radiance to account for impact flash magnitude, and most of the light from impact flashes is created by thermal radiance of droplets of melted ejecta [1]. However, they did not fit their model to the observed diameter and depth of an associated crater.

Method:
We seek to improve the constraints on luminous efficiency by modelling simplified analogues of observed impact flashes with known associated craters. Here, we use the iSALE-Dellen shock physics code [6,7,8] to model the 17 March 2013 lunar impact flash and crater [Fig. 1], adjusting the impact parameters to best replicate the known factors of the impact and match the resultant crater size. Target porosity, modelled by the ε-α-porosity compaction model [8,9], is set at a homogeneous 42% and a vertical impact in a 2D cylindrical model is used for computational efficiency. We use the Drucker-Prager strength model and the ANEOS equation of state to model the lunar regolith. Then, we determine the volume of melted ejecta produced in each impact as a proxy for radiance, and vary velocity-mass balance and target porosity parameters to study the effect on ejected melt volume. Crater data permitting, we will also study additional, recently identified flash-crater pairs [10,11,12].

 

Figure 1: Crater identified with 17 March 2013 impact flash. a: before, b: after. [13]

 

Results:

We show some preliminary results of our first study models in the plot below. Peak pressure on target material is plotted (in Pa x1e10). Material above the black curve is ejected as the model proceeds. All material is shown in position at t=0. The resolution of this model (10 CPPR) is not sufficient for a detailed analysis of peak pressure and melt volume; 40 CPPR resolution will yield a better result, reducing uncertainty to <10% [14]. Critical shock pressure for melt production in regolith at 42% porosity is estimated at no higher than 24 GPa [14]. Our preliminary model shows only a small sliver of overlap between melt volume and ejecta volume, much lower than the 20-30% of melt ejected in [15], suggesting a low luminous efficiency in this impact.

Figure 2: Peak pressure on target material for modelled 17 March 2013 impact. Material above the black curve is ejected. Material above 24 GPa peak pressure is assumed melted.

 

Conclusion:

Further study of the details of thermal radiance of the ejected melt will be left to a future study, enabling a connection from impact parameters to luminous efficiency via melt volume. This work could also be expanded through use of 3D models with realistic replication of the estimated impact angles of known impact flash craters, as well as more detailed, layered target rock properties.

 

References

[1] Yanagisawa, Masahisa et al. (2002), Icarus 159 (1), pp. 31–38. DOI: 10.1006/icar.2002.6931.

[2] Sheward, D. et al. (2024), MNRAS 529 (4), pp. 3828–3837. DOI: 10.1093/mnras/stad2707.

[3] Avdellidou, Chrysa et al. (2021), Planetary and Space Science 200, p. 105201. DOI: 10.1016/j.pss.2021.105201.

[4] Ortiz, J. L. et al. (2000), Nature 405 (6789), pp. 921–923. DOI: 10.1038/35016015.

[5] King, Patrick K. et al. (2022), 16th Hypervelocity Impact Symposium. DOI: 10.1115/HVIS2022-28.

[6] Amsden, A. et al. (1980), LANL, LA8095:101.

[7] Collins, G.S. et al. (2004), Meteoritics & Planetary Science 39:217–231. DOI: 10.1111/j.1945-5100.2004.tb00337.x.

[8] Wünnemann, K. et al. (2006), Icarus 180, pp. 514-527. DOI: 10.1016/j.icarus.2005.10.013.

[9] Collins, G. et al. (2011), Int. J. Imp. Eng., 38:434-439. DOI: 10.1016/j.ijimpeng.2010.10.013,

[10] Sheward, Daniel et al. (2021), Europlanet Science Congress 2021. DOI: 10.5194/epsc2021-590.

[11] Sheward, D. et al. (2022), MNRAS 514 (3), pp. 4320–4328. DOI: 10.1093/mnras/stac1495.

[12] Sheward, Daniel et al. (2024), EGU General Assembly 2024, Abstract EGU24-1216. DOI: 10.5194/egusphere-egu24-1216.

[13] Robinson, Mark S. et al. (2015), Icarus 252, pp. 229–235. DOI: 10.1016/j.icarus.2015.01.019.

[14] Liu, T. et al. (2022), JGR Planets 127 (8), Article e2022JE007264, e2022JE007264. DOI: 10.1029/2022JE007264.

[15] Luther, R. et al. (2017), 48th Annual Lunar and Planetary Science Conference (1964), p. 3012. Bibcode: 2017LPI....48.3012L.

How to cite: Rice, P., Luther, R., and Wünnemann, K.: Modelling Lunar Impact Flashes from Molten Ejecta, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-560, https://doi.org/10.5194/epsc2024-560, 2024.

EPSC2024-956 | ECP | Posters | TP9 | OPC: evaluations required

From Debris to Resource: Simulating High-Velocity Impacts of Space Debris on the Moon 

Nick Langer, Robert Luther, Kai Wünnemann, Frank Koch, and Stefan Linke
Thu, 12 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P7

1. Introduction

The quantity of space debris—defunct human-made objects in orbit—continues to increase, including old satellites, spent rocket stages, and fragments resulting from disintegration, erosion, and collisions. Space debris poses a significant risk for space travel and satellite operations, as even small pieces, traveling at high velocities, can cause substantial damage or destruction upon colliding with satellites or spacecraft. This growing accumulation increases the risk of collisions, poten tially triggering a cascade effect known as the Kessler Syndrome, which could render certain orbits unusable for future satellites.

Despite the impending threat to space infrastructure, measures to mitigate space debris are insufficiently implemented, primarily due to high associated costs. Furthermore, existing regulations are often more akin to recommendations than mandatory directives. A solution with potential economic benefit is to collect larger objects (e.g. a Ariane 5 upper stage in the geostationary transfer orbit), and to transport this material to the lunar surface, where the metal can be used in future projects [1]. To facilitate recycling, the debris is planned to impact the Moon at velocities ranging from 0.5 − 3.0 km/s. Upon impact, the debris would break down, and the resulting fragments would be collected by moon rovers. In this study, we analyze the effect of impact velocity and target parameters on the fate of the projectile and the resulting crater, using numerical simulations.


2. Methods

In this study, we apply iSALE-2D shock physics code [2-4] to simulate the impact of an Ariane 5 ESC-A upper stage (4540 kg, 5.4 m in diameter, 4.711 m in height [5]) onto lunar regolith at velocities from 0.5 − 3 km/s. A realistic representation with thin walls is computationally expensive as it requires a large model resolution. Here, we consider four representations of the projectile: We represent the Ariane upper stage as 1) a small homogeneous cylinder with the density of aluminum alloy 2219 (ρ=2840 kg/m³), 2) a porous cylinder of the same mass, and a correct volume of the projectile, using the epsilon-alpha porosity compaction model [4], 3) a homogeneous cylinder of the same mass and volume but instead of using the porosity model, we modify the density specified in the equation of state, and 4) a hollow cylinder with resolved walls and end caps. We use a grid spacing of 0.2 and a resolution of 14 cells per projectile radius in vertical and horizontal direction.

The projectile and target are simulated by a Tillotson equations of state [6]. It is designed for high-velocity impact computations and hence capable of handling the relevant pressure ranges for this study. A granular behavior is assumed for the regolith target. Therefore the Drucker-Prager strength model with the strength Y=min(Y0p,Ym) is used, where Y0 is the cohesion (yield strength at zero pressure), μ is the coefficient of internal friction for material, Ym is the limiting strength at high pressure and p is pressure. The parameters for the regolith are taken from [7], including porosity model parameters. We apply the Von Mises strength model for the aluminum of the upper stage: Y=Y0. For the porous projectile representation, we derive an elastic threshold based on the Youngs modulus and the yield strength as: ε=315 MPa / 73 GPa=0.004. Both values for this equation are taken from [8].


3. Results

Here, we show results generated using a porous projectile. We observe, as expected, that with increasing kinetic energy, the radius and depth of the crater increase as well (Figure 1). For the slowest impact of 0.5 km/s, the crater grows to ∼ 18 m in diameter, while it reaches ∼ 32 m for an impact at 3 km/s. The faster impact also generates a more heterogeneous pressure distribution (Figure 2).

We find that in the case of the slower impact, the resulting crater is relatively shallow (Figure 2 top panel), with a depth-diameter ratio of 0.17, while it reaches 0.25 for the fastest impact (Figure 2 bottom panel).


Figure 1: Simulated crater radii and depths from the instant of impact until 4 seconds after impact.

 


(2a) Impact velocity vimpact = 0.5 km/s .


(2b) Impact velocity vimpact = 3.0 km/s .
Figure 2: Crater snapshot at 4 seconds after impact. The pressure and density are plotted on the left and right halves, respectively.


4. Discussion & Conclusion

The plotted depth in Figure 1 (lower panel) occasionally deviates from the final value, which is likely due to numerical artefacts or residual material close to the symmetry axis. An adjustment of the analysis procedure can reduce this noise.

The study investigates lunar surface recycling of space debris, specifically through the high-velocity impact simulation of an Ariane 5 upper stage on lunar regolith. The method's potential in repurposing space debris for lunar construction is promising, contributing to sustainable lunar exploration efforts. However, anomalies in pressure distribution and crater rim identification suggest simulation refinement.


References

[1] F. Koch. 8th European Conference on Space Debris, 8(1), 2021.
[2] A.A. Amsden, H.M. Ruppel, and C.W. Hirt. Technical Report LA-8095, 5176006, June 1980.
[3] G.S. Collins, H.J. Melosh, and B.A. Ivanov. Meteoritics & Planetary Science, 39(2):217–231, February 2004.
[4] K. Wünnemann, G.S. Collins, and H.J. Melosh. Icarus, 180(2):514–527, February 2006.
[5] R.Lagier et al. User’s Manual, (1), 2021.
[6] J. H. Tillotson. Technical report, July 1962.
[7] R. Luther et al. The Planetary Science Journal, 3(10):227, October 2022.
[8] C.Y. Ho, J.M. Holt, and H. Mindlin. Number Bd. 2. Cindas/Purdue Univ., 1997.

How to cite: Langer, N., Luther, R., Wünnemann, K., Koch, F., and Linke, S.: From Debris to Resource: Simulating High-Velocity Impacts of Space Debris on the Moon, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-956, https://doi.org/10.5194/epsc2024-956, 2024.

EPSC2024-1111 | ECP | Posters | TP9 | OPC: evaluations required

Quantifying the effect of asteroid structure on hypervelocity impact outcomes with the Material Point Method 

Xiaoran Yan, Wenhan Zhou, Yan Liu, Patrick Michel, and Junfeng Li
Thu, 12 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P8

Introduction Recent advancements in computational power and model sophistication have revolutionized the study of hypervelocity impacts on asteroids, enabling the exploration of broader parameter spaces and more intricate asteroid structures [1, 2]. The outcome of such collisions, including crater formation and catastrophic disruption, is intricately linked to the internal structure of asteroids. However, the nature of this relationship remains an open question, necessitating further investigation. The Material Point Method (MPM), renowned for its ability to accurately track interfaces and resolve contact problems, emerges as a powerful tool for numerical simulations in this context [3]. By harnessing the capabilities of MPM, we aim to shed light on the correlation between an asteroid's structural composition and its response to hypervelocity impacts, ultimately contributing to a deeper understanding of the underlying physical processes.

Method Developed from the particle-in-cell numerical method, the Material Point Method (MPM) combines Eulerian and Lagrangian descriptions, making it well-suited for handling boundary conditions and interface problems. By utilizing a dynamic background grid that only creates nodes with mapped material points for efficient computation, we further implement a contact algorithm in our MPM framework that could accurately calculate the real contact and friction force on each interface, and novelly extend the contact correction from two-object to multi-object scenarios.

When the contact algorithm is activated, if a background grid node is shared by multiple objects, separate node information will be created for each object, along with the associated material point indices, as illustrated in Fig. 1. The execution of contact correction is triggered by assessing both the distance between objects and their relative motion based on non-penetration constraints. It is implemented by applying contact forces to eliminate penetration and enforce suitable friction conditions. Multi-body contact is resolved by solving a linear system of equations for the unknown contact forces between each pair of interacting objects. This algorithmic enhancement bolsters the capability of our MPM code to precisely simulate hypervelocity impacts on asteroids with intricate internal structures.

Simulation and Discussion In this study, we conduct numerical simulations of artificial hypervelocity impacts on asteroids to quantitatively investigate the distinct dynamic responses of rubble-pile asteroids compared to monolithic, homogeneous ones. We employ material parameters of typical basaltic rocks in our simulations [4]. The impact scenarios and results for rubble-pile asteroids are illustrated in Fig. 2. Considering the crucial role of the catastrophic disruption threshold in understanding the evolutionary history of small celestial bodies and designing asteroid impact defense missions, as well as the availability of numerical values for specific asteroid sizes from previous research [5], we aim to quantify the influence of porosity and internal interfaces on impact disruption mechanisms by calculating the catastrophic disruption threshold for rubble-pile asteroids with various configurations.

Our simulation results demonstrate the significant impedance effect of internal interfaces on shock wave propagation, leading to a substantially higher catastrophic disruption threshold for rubble-pile asteroids compared to monolithic asteroids of the same strength. This finding suggests that the presence of a rubble-pile structure enhances the resistance of asteroids to catastrophic disruption, which has important implications for understanding their fragmentation and reassembly processes throughout their evolutionary history.

Furthermore, the quantitative analysis of the catastrophic disruption threshold for different rubble-pile configurations provides valuable insights into the role of porosity and internal structure in the impact response of asteroids. These findings contribute to the development of more accurate models for asteroid evolution and can inform the design of future asteroid impact defense strategies. Additionally, the results of this study have the potential to facilitate the calibration of asteroid ages based on their structural properties and impact history, enhancing our understanding of the chronology of asteroids in the solar system.

Acknowledgment X.Y. and J.L. acknowledge support from the National Natural Science Foundation of China under Grant 12372047. X.Y. acknowledges support from the National Natural Science Foundation of China under Grant 62227901. P.M. acknowledges support from the French space agency CNES and from the French National Centre for Scientific Research (CNRS) through the exploratory research program of the Mission for Transversal and Interdisciplinary Initiatives.

References [1] Raducan, S. D. et al. (2024). PSJ. 5(3). [2] Jutzi, M. et al. (2022) Nat. Commun. 13(1), 7134. [3] Yan, X. et al. (2023) ACM Conference. [4] Jutzi, M. (2015) P&SS. 107(1), 3–9. [5] Benz, W. et al. (1999) Icar. 142, 5–20.

How to cite: Yan, X., Zhou, W., Liu, Y., Michel, P., and Li, J.: Quantifying the effect of asteroid structure on hypervelocity impact outcomes with the Material Point Method, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1111, https://doi.org/10.5194/epsc2024-1111, 2024.

TP10 | Exploring Mercury and its environment

EPSC2024-247 | ECP | Posters | TP10 | OPC: evaluations required

Shape and Albedo from Shading with Planetary Flyby Images of Mercury and the Moon 

Isabel Krüll, Kay Wohlfarth, Moritz Tenthoff, Christian Wöhler, Valentina Galluzzi, Jack Wright, Johannes Benkhoff, and Joe Zender
Thu, 12 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P26

Introduction
Surface reconstruction of planetary bodies such as the Moon and Mercury is crucial for geomorphological analysis, reflectance normalization, thermal modeling, rover landing site planning, and outreach activities. Stereo algorithms and Shape-and-Albedo-from-Shading (SAfS) are well-established methods for planetary 3D reconstruction. The current state-of-the-art combines both methods. SAfS refines the surface slopes of a stereo digital elevation model (DEM) and typically yields 3D models at image resolution [1,2,3,4,5,6]. This approach is generally well-validated for scientifically calibrated instruments that observe the planetary body under favorable conditions. However, the limits of SAfS still need to be explored. This work applied the SAfS algorithm to planetary flyby images acquired with uncalibrated off-the-shelf cameras. We qualitatively and quantitatively assessed the algorithm's performance and found the method robust even under these challenging conditions.

Methods
We considered two images: First, a flyby image of Mercury (Figure 1, left), which was obtained with a monitoring camera during BepiColombo’s third flyby, and second, a flyby image of the Moon captured by a GoPro during the Artemis I mission (Figure 1, right). In both cases, off-the-shelf cameras without a proper radiometric calibration routine were used instead of scientific instruments. Therefore, it was necessary to calibrate the images before applying the SAfS algorithm. For calibration, we estimated a reflectance image of the region of interest with Hapke parameters from [7] to establish a relationship between the digital number output of the camera and the physical radiances. The resulting Mercury flyby DEM was evaluated qualitatively, compared to MDIS WAC images [8]. The lunar flyby DEM was compared to the SLDEM2015 (60-100 m/pixel) [9], which serves as a ground truth. 

Figure 1. Left: Flyby image from the BepiColombo mission [10]. Right: Flyby image from the Artemis I mission [11].

 

Results
Figure 2 shows the color-coded SAfS DEM for the BepiColombo image. We found that the algorithm successfully reconstructed the surface in the center of the image but struggled with the more extremely illuminated sections at the edges. It is obvious that the algorithm reconstructed some details that are not visible in the input DEM and hence improved the resolution. Figure 3 compares a grey-scale representation of the SAfS DEM with MDIS WAC image EW0251718878F. Small craters in Izquierdo (the crater in the right half of the marked section), a few kilometers in diameter, become especially visible. Figure 4 shows the marked section in more detail. Due to the lack of high-resolution ground truth, a detailed algorithm evaluation was not possible for this image.

Figure 2. Color-coded presentation of the reconstructed SAfS DEM from the BepiColombo image.

 

Figure 3. Left: grey scale reconstructed SAfS DEM from the BepiColombo image. Right: wide-angle camera (WAC) image from MDIS [8]. The marked image section was evaluated in more detail (Fig. 4).

 

Figure 4. Comparison of WAC image, input DEM, BepiColombo (BC) image and SAfS DEM. Below the images are the height and slope of the profile (red dashed line).

 

However, the ground truth evaluation with the Artemis I image gives a quantitative measure under comparable conditions. Figure 5 shows the color-coded SAfS DEM of a region of interest (ROI) reconstructed from the flyby image. The results for the Artemis I image for all ROIs were of high quality. Figure 6 shows the elevation profile indicated by the dashed line in Figure 5. The algorithm refines the low-frequency initial DEM (dashed line), yielding a SAfS DEM (red line), which closely resembles the ground truth DEM (black line).  The vertical RMSE between the reconstructed DEM and the ground truth is 523 m, lower than the pixel size of approximately 1500 m. However, there were inaccuracies near the ROIs’ edges, and a preferred direction aligned with the illumination direction became visible.

Figure 5. Color-coded presentation of the reconstructed SAfS DEM from the Artemis image. The black dashed line marks the profile that was analyzed in detail (see Fig. 6).

 

 

Figure 6. Height profile of a selected terrain profile (see black dashed line in Fig. 5). Red line: DEM generated with the SAfS algorithm. Black line: Ground truth DEM. Dashed line: Initial DEM (input for the SAfS algorithm).

 

Conclusion
In conclusion, it is possible to obtain sharp results by applying our SAfS framework to flyby images. Both results show that, despite the challenging conditions, the SAfS algorithm could reconstruct the surface up to image resolution and increase the level of detail of the input DEM. The quality differences between the two images can mainly be attributed to the (spatial) resolution of the original images and the oblique illumination direction. Usually, the image center is distortion-free, and the illumination geometry is best suited for SAfS. We found that the surface reconstruction at the edge of the image is also possible, but the quality decreases significantly. All in all, our flyby-derived DEMs are accurate, and a previous version has been used for ESA outreach activities, similar to [12]:
https://www.esa.int/Science_Exploration/Space_Science/BepiColombo/BepiColombo_s_third_Mercury_flyby_the_movie

 

References
[1] A. Grumpe, F. Belkhir, C. Wöhler. Advances in Space Research, 53(12):1735–1767, 2014.
[2] C. Jiang, S. Douté, B. Luo, L. Zhang, P&RS,130, 2017 
[3] O. Alexandrov, R. Beyer. Earth and Space Science, 5, 2018
[4] B.  Wu, W. C.  Liu, A. Grumpe, C. Wöhler. P&RS, 140, 2018
[5] M. Tenthoff, K. Wohlfarth, C. Wöhler. Remote Sensing, 12(23), 2020.
[6] M. Hess, M. Tenthoff, K. Wohlfarth, and C. Wöhler. Journal of Imaging, 8(6), 2022.
[7] J. Warell, Icarus, 167, 2,2004
[8] Hawkins, S. Edward, et al. Space Science Reviews 131 (2007): 247-338.
[9] M.K. Barker, E. Mazarico, G.A. Neumann, M.T. Zuber, J. Haruyama, D.E. Smith, Icarus, 273, 2016.
[10] ESA. Planetary science archive.2023.
https://archives.esac.esa.int/psa/#!Image%20View/MCAM=instrument, accessed 10th May 2024
[11] NASA. flickr, 2022,
https://www.flickr.com/photos/nasa2explore/52547180935/in/album-72177720303788800/, accessed 10th May 2024
[12] K. Wohlfarth, M. Tenthoff, J. Wright, V. Galluzzi, C. Wöhler, H. Hiesinger, J. Helbert, J. Zender, J. Beckhoff. MExAG Annual Meeting, 02.2023

How to cite: Krüll, I., Wohlfarth, K., Tenthoff, M., Wöhler, C., Galluzzi, V., Wright, J., Benkhoff, J., and Zender, J.: Shape and Albedo from Shading with Planetary Flyby Images of Mercury and the Moon, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-247, https://doi.org/10.5194/epsc2024-247, 2024.

EPSC2024-646 | ECP | Posters | TP10 | OPC: evaluations required

Spectral properties of pyroclastic deposits on Mercury over space and time 

Mireia Leon Dasi, Sebastien Besse, Lauren M. Jozwiak, Erica R. Jawin, and Alain Doressoundiram
Thu, 12 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P34

Explosive volcanic activity on Mercury extended after the end of the widespread effusive volcanism era (Jozwiak et al., 2018). Understanding the precise timing of explosive eruptions has significant implications for the volatile content and thermal evolution of the planet. However, the age of individual pyroclastic deposits remains largely debated and is difficult to assess using crater counting (e.g. Luchitta and Schmitt, 1974). An individual analysis of a selection of faculae has highlighted the spectral diversity across and within deposits (Barraud et al., 2021, Besse et al., 2020). In this work, we constrain the link between the variability in spectral properties across deposits and deposit age. Additionally, we explore the spatial variability in spectral properties inside deposits, and the relationship of this variability to the timing of eruptions inside individual faculae. A combination of morphologic analyses based on MESSENGER/MDIS images and spectral data from MESSENGER/MASCS is analyzed to this end, utilizing a deep learning approach.

We analyze the relationship between the morphological degradation of the vents (physical depressions) and the spectral changes in the associated deposits (faculae). This study shows a correlation between the deposit spectra and vent degradation, characterized by a rapid initial darkening of the deposit and spectral flattening over time, followed by stabilization. The deposits with heavily degraded vents reach the properties of the local background terrain, rendering old deposits spectrally undetectable. To explain these temporal variations in spectral properties, we propose three potential processes: space weathering, mixing with the underlying terrain, and changes in erupted pyroclast size. Space weathering acts “fast”: spectral changes induced by nanophase iron accumulation produced by space weathering on the Moon saturate after ~1 Ga (Tai Udovicic et al., 2021). If a similar mechanism is responsible for most of the spectral modifications observed over time, then a large part of explosive eruptions detected on Mercury could be significantly younger than previously expected.

To explore the relative timing of vents within the same facula, we examine the variability of spectral properties inside the deposits. Using outlines defining the vent profile and the deposit extent defined by Leon-Dasi et al. (2023), we extract the evolution of spectral properties in the direction locally perpendicular to the vent. Using these data, we study (1) the rate of change of spectral properties, (2) the symmetry of the deposit, and (3) the contribution of each vent to the spectral variability. This analysis benefits from leveraging the deep learning-based data reduction performed by Leon-Dasi et al. (2023), which results in a set of latent dimensions integrating spatial and spectral information. Using such latent dimension has proven to highlight the pyroclastic deposits and reduce the spatial noise more effectively than using single MASCS-derived spectral parameters. From a preliminary analysis, we find a trend between the rate of change of the spectral variations across the deposits and the deposit age. Older deposits appear to present slower-changing spectral properties, which is consistent with deposit erasure over time through various space weathering processes. Within multi-vent faculae, we find spatial variations in the association between spectral signatures and different vents—this implies that eruptions within a single facula were separated in time. Overall, this research provides further insight into Mercury’s volcanic history and the processes that have shaped its surface over time.

How to cite: Leon Dasi, M., Besse, S., Jozwiak, L. M., Jawin, E. R., and Doressoundiram, A.: Spectral properties of pyroclastic deposits on Mercury over space and time, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-646, https://doi.org/10.5194/epsc2024-646, 2024.

EPSC2024-757 | ECP | Posters | TP10 | OPC: evaluations required

Linewidth Measurements of Mercury's Alkali Exosphere 

Patrick Lierle, Carl Schmidt, and Emma Lovett
Thu, 12 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P17

Sunlight can shape elements in Mercury’s thin atmosphere into an escaping cometlike tail. Solar photons are absorbed from one direction—the Sun—and scatter isotropically, imparting a net momentum in the anti-sunward direction. This change in momentum is a force that is strongest on atoms that efficiently interact with sunlight, that is, atoms with strong resonance scattering like sodium and potassium. Sunlight is intense at Mercury, and this creates a long sodium tail that extends more than 1000 planetary radii, making it one of the largest structures in our solar system [1].

Precision radial velocity spectrometers are common tools in the exoplanet community that can offer resolved linewidth measurements of Mercury’s exosphere. Here we present observations of the sodium and potassium exosphere from two such instruments. With R~150,000 resolving power and fast tip-tilt image stabilization, the Extreme Precision Spectrometer (EXPRES) at the 4.3m Lowell Discovery Telescope is ideally suited for measuring line broadening in planetary gases [2][3][4]. Over several nights across April 2023, we sampled from 0 to 4 planetary radii along the cometlike tail with EXPRES. Data were obtained surrounding maximum radiation pressure (TAA = 56–80°) and at a 90° phase angle, where Mercury’s tail is oriented perpendicular to the line of sight. Our downtail pointing from 2023 April 10 is illustrated in Figure 1. In addition to probing the exotail, we also acquired spectra on disk to look for small-scale variations in linewidths. In March 2024, we completed a follow-up campaign with the Keck Planet Finder (KPF) spectrometer on the Keck I telescope, once again sampling several planetary radii downtail and more intentionally targeting on-disk regions of interest (cusps, subsolar point, poles, etc.). This run spanned a different season at Mercury (TAA = 20–40°) with the planet still near 90° phase angle. Results from EXPRES are discussed herein and early results from the follow-up KPF campaign will be presented in-person at EPSC2024.

In order to quantify measured linewidths, we derive effective temperatures. While the collisionless exosphere is not inherently thermal, effective temperature estimates are nonetheless a useful energy metric. We obtain these estimates by convolving forward models of the Doppler-broadened hyperfine structure of the sodium and potassium D lines with the instrumental line spread functions of EXPRES and KPF. These convolved models are fit to the solar-subtracted spectra and best-fit temperatures are extracted via a least-squares algorithm. Pointing is determined by matching spectra to slit-viewer images. As verification of this method, we measure sodium gas above the low-latitude dayside limb to be approximately 1200K, consistent with estimates derived from MESSENGER scale heights [5].

We find that both sodium and potassium line profiles exhibit steep growth with downtail distance until their effective temperatures level off near 8000 and 10,000 K, respectively, around 3 planetary radii. We interpret this to be the furthest extent of gravitationally bound gas, where atomic trajectories reach apex and return toward the surface. Beyond 3 radii, the escaping gas populations show a constant effective temperature with distance. Figure 2 shows sodium D1 and D2 profiles from 2023 April 10. As we point downtail, the gravitationally trapped cold population forming the core of the lines is depleted, causing broadening and producing the nonthermal profiles seen past 1.5 radii. The asymmetry in these sodium profiles is consistent across nights and appears to indicate more particles are moving away from us than toward us along the line of sight. Figure 3 plots the leveling off of effective temperature at two TAAs. At maximum radiation pressure, the effective temperature of the distant tail approaches 8000 K, as compared to 5000 K during a lower g-value season. This is expected, as radiation pressure accelerates particles to higher velocities during the peak g-value season. The potassium exosphere is less extended than sodium at only ~700 K above the dayside limb, but both species exhibit significant escape during peak radiation pressure.

Results from on-disk measurements with EXPRES show that sodium above the low-latitude dayside limb is approximately 1200 K, while gas at high latitudes is 100-200 K more energetic. Despite bright enhancements near the magnetic cusps, linewidths here show no evidence for an ion-sputtered component with energies predicted by theory or laboratory time of flight experiments. KPF measurements will help further constrain the contribution of a hot sputtered component. A cold population with enhancement at the cusps may suggest a plasma-stimulated low-energy exospheric source such photon or electron stimulated desorption.  

Figure 1: EXPRES pointing on 2023 April 10 with triangles representing the spectrometer aperture. We sampled the southern lobe of the tail from 0 to 4 radii downtail.

Figure 2: EXPRES line profiles of sodium D1 and D2 at increasing downtail distance from Mercury. Spectra are color-coded to match pointing locations in Figure 1. Data are overlaid with thermal forward models in grey, illustrating how, as linewidths grow, they also become increasingly nonthermal. Best fit temperatures are given in grey. Downtail distances are zeroed at the nightside limb. Temperatures increase from ~1200K on disk to above 7000K past 3 Mercury radii downtail.

Figure 3: Sodium best-fit temperatures to EXPRES data as a function of downtail distance for two true anomaly angles. Temperature levels off just before 3 Mercury radii downtail for both, though effective temperatures far downtail are much higher for TAA=56°.

References:

[1] Baumgardner, J., Wilson, J., & Mendillo, M, 2008, GRL, 35(3), L03201.

[2] Petersburg R. R., Joel Ong J. M., Zhao L. L. et al. 2020 AJ 159 187.

[3] Jurgenson C., Fischer D., McCracken T. et al. 2016 Proc. SPIE 9908 99086T.

[4] Brewer J. M., Fischer D. A., Blackman R. T. et al. 2020 AJ 160 67.

[5] Cassidy, T. A., Merkel, A. W., Burger, M. H., et al., 2015, Icarus, 248.

How to cite: Lierle, P., Schmidt, C., and Lovett, E.: Linewidth Measurements of Mercury's Alkali Exosphere, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-757, https://doi.org/10.5194/epsc2024-757, 2024.

EPSC2024-828 | Posters | TP10

Preliminary temperature analysis of the Region of Interest using MERTIS onboard BepiColombo for the upcoming Mercury's 5th Flyby. 

Nimisha Verma, Joern Helbert, Aurelie Van den Neucker, Mario D'Amore, Solmaz Adeli, Giulia Alemanno, Oceane Barraud, Alessandro Maturilli, Karin Bauch, and Harald Hiesinger
Thu, 12 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P37

Introduction:

The MErcury Radiometer and Thermal infrared Imaging Spectrometer (MERTIS) is an Infrared spectrometer (TIS) and radiometer (TIR) instrument onboard the BepiColombo spacecraft. It is part of the Mercury Planetary Orbiter payload with a spectral wavelength of 7-14 μm and a resolution of 90 nm. The radiometric wavelength of MERTIS is 7-40 μm (Hiesinger et al., 2008). Currently, the MERTIS team is preparing for the upcoming 5th flyby of Mercury by BepiColombo on the 2nd of December, 2024 during which the instrument will observe the Hermean surface to characterize the spectral emission, and map the surface mineralogy and temperature variation of the planet. The observation will be performed using the space view of the instrument.

Figure 1: MERTIS Instrument (ESA, 2024)

Region of Interest:

The currently planned Region of Interest (ROI) expands between 51.54° & -53.18° latitude and -97.17° & -143.30° longitude. It encompasses the Beethoven basin and craters like Michelangelo, Durer and Vieira da Silva. The altitude in the region ranges from 2841 m to -4453.5 m with a maximum slope of around 30°.

Figure 2: Left: Digital Elevation model of the ROI. Right: Slope (in degree) of the ROI.

Methodology:

Due to proximity to the Sun, the surface of Mercury undergoes a large temperature variation while also showcasing difference in regional temperature identified as hot and cold regions. Various factors influence the temperature on Mercury’s surface like the density of impact craters, topography, surface morphology, distance to the sun, density of the material etc. In preparation for the flyby, a preliminary temperature analysis of the ROI is conducted using python language and equations from the Vasavada model.

The surface temperature of Mercury is mainly dependent on three parameters – the radiation received from the sun, the radiative loss of heat and the thermal conduction of the surface and sub-surface (Bauch et al., 2021; Yan et al., 2005). Solar irradiance on the surface of Mercury varies depending on the distance from the sun and the longitude of observation (Yan et al., 2005). During the flyby, mercury will almost be at its closest to the sun with an average distance of 46,959,173 km for the observation period between longitude -97° and -143°. The incidence angle for the period of observation ranges from 81° to 83°. It has been observed that the energy received at longitude 90° and 270° is almost half of the energy received at 0° and 180° during perihelion (Yan et al., 2005). Taking into consideration the above information, the solar irradiance is calculated based on the incidence angle of the sun, the distance and the surface albedo.

The radiative heat loss from the planet during daytime is defined by the upper boundary condition and the lower boundary condition. The upper boundary condition is the calculation of the radiative heat at the surface while the lower boundary is for the sub-surface radiation. For simplifying the process, we do not consider the sub-surface conditions of the planet. Hence, the lower boundary condition is ignored.

In order to calculate the thermal conductivity, we consider the ratio of two parameters – contact conductivity which refers to the ability of two material to conduct heat or electricity through the point of contact and conduction by radiation into a medium. Aubrite is one of the materials considered for the calculation of the thermal conductivity due to its characteristic proximity to the planet Mercury (Keil, 2010).

The results generated from these calculations will be used to develop a temperature map for the ROI. This map will be used along with the emissivity spectral measurement from Planetary Spectroscopy Laboratories (PSL) of the German Aerospace Center (DLR), to better understand the temperature ranges in ROI and characterize the surface mineralogy for the upcoming BepiColombo flyby.

References:

Bauch, K.E., Hiesinger, H., Greenhagen, B.T., Helbert, J., 2021. Estimation of surface temperatures on Mercury in preparation of the MERTIS experiment onboard BepiColombo - ScienceDirect. Icarus 354. https://doi.org/10.1016/j.icarus.2020.114083

Hiesinger, H., Helbert, J., MERTIS Co-I Team, 2008. The Mercury Radiometer and Thermal Infrared Spectrometer (MERTIS) for the BepiColombo mission - ScienceDirect. Planetary and Space Science 58, 144–165. https://doi.org/10.1016/j.pss.2008.09.019

Keil, K., 2010. Enstatite achondrite meteorites (aubrites) and the histories of their asteroidal parent bodies. Geochemistry 70, 295–317. https://doi.org/10.1016/j.chemer.2010.02.002

Yan, N., Chassefière, E., Leblanc, F., Sarkissian, A., 2005. Thermal model of Mercury’s surface and subsurface: Impact of subsurface physical heterogeneities on the surface temperature - ScienceDirect. Advances in Space Research 38, 583–588. https://doi.org/10.1016/j.asr.2005.11.010

How to cite: Verma, N., Helbert, J., Van den Neucker, A., D'Amore, M., Adeli, S., Alemanno, G., Barraud, O., Maturilli, A., Bauch, K., and Hiesinger, H.: Preliminary temperature analysis of the Region of Interest using MERTIS onboard BepiColombo for the upcoming Mercury's 5th Flyby., Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-828, https://doi.org/10.5194/epsc2024-828, 2024.

EPSC2024-1011 | ECP | Posters | TP10 | OPC: evaluations required

Mineralogy of the mantles in sub- Earths and exo- Mercuries 

Camilla Cioria, Giuseppe Mitri, James Alexander Denis Connolly, Jean-Philippe Perrillat, and Fabrizio Saracino
Thu, 12 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P15

The mantle mineralogy of large exoplanets (e.g. super- Earths, planetary bodies having masses between 1-10 ME) has been widely investigated (Duffy et al., 2015).

However, the minerals constituting the mantle of Mercury-sized objects (here referred to as exo-Mercuries) and sub-Earths (planetary bodies having masses ranging from 1 MM < 1ME, respectively 1 Mercury-mass and 1 Earth-mass), orbiting closer to their stars, are still underexplored.  This modeling work has focused on describing stable mineral associations in those mantles equilibrated under low values of oxygen fugacity (fO2). Such reducing conditions are not uncommon in stellar systems, as evidenced by various materials in the solar system, ranging from undifferentiated ones (carbonaceous chondrites belonging to the CH and CB groups) to enstatite chondrites, to already differentiated materials like aubrites, and even entire planets like Mercury.

Assuming Mercury as a proxy, we employed the open-source software Perple_X (Connolly, 1990) to characterize the mineral assemblages forming the mantle of these reduced planetary bodies.

The thermodynamic approach here adopted offers the advantage of allowing the investigation of those planetary interiors otherwise not explorable. Moreover, the employed thermodynamic inputs were extrapolated from known precursor materials, which share several properties with exoplanets under examination. This methodology has already been discussed in literature (Néri et al., 2020; Cioria and Mitri, 2022), ensuring the validity of this approach. Bulk silicate compositions of aubrite and CH, CB, EN chondrites, have been used as thermodynamic inputs in our simulations. Calculations were conducted within the pressure and temperature ranges suggested for the mantle of Mercury: 1200 K -1700 K and 3 GPa -5 GPa (Tosi et al., 2013).

 We found that orthopyroxene, clinopyroxene, olivine, and accessory minerals constitute the mantles of Mercury and reduced exoplanets. These results indicate that the initial bulk compositions have a first-order constraint on the resulting mantle mineralogy; more specifically, the initial abundance of SiO2 determines the chemical equilibrium shift between pyroxene and olivine, respectively constituting the dominant phases at high and low silica content. The predicted mantle mineralogy, dominated by pyroxenes, is consistent with that outlined by Putirka and Rarick, (2019) and Putirka and Xu, (2021).

The differences with Earth’s mantle mineralogy are significant, necessitating a new and more appropriate classification of these rocks, as already suggested in Putirka and Rarick, (2019).

This outcome holds significant implications for the thermochemical evolution and geodynamics of reduced mantles, also exerting a  substantial influence on the properties of the related  crusts and cores.

Future investigations on Mercury, conducted by the ESA BepiColombo mission, could shed new light on the possible mineralogy of its mantle, helping to further detail the minerals stable in those exoplanets formed in similar geochemical contexts.

Acknowledgments

G.M. and C.C. acknowledge support from the Italian Space Agency (2017-40-H.1-2020).

References

Cioria, C., & Mitri, G. (2022). Model of the mineralogy of the deep interior of Triton. Icarus, 388, 115234.

Connolly, J. A. D. (1990). Multivariable phase diagrams; an algorithm based on generalized thermodynamics. American Journal of Science, 290(6), 666- 718. https://doi.org/10.2475/ajs.290.6.6661315

Duffy, T., Madhusudhan, N., and Lee, K.K.M. (2015) Mineralogy of super-Earth planets. Treatise on Geophysics, 2nd Ed., 2, 149-178.1333

Néri, A., Guyot, F., Reynard, B., & Sotin, C. (2020). A carbonaceous chondrite and cometary origin for icy moons of Jupiter and Saturn. Earth and Planetary Science Letters, 530, 115920.

Putirka, K. D., & Rarick, J. C. (2019). The composition and mineralogy of rocky exoplanets: A survey of > 4000 stars from the Hypatia Catalog. American Mineralogist: Journal of Earth and Planetary Materials, 104(6), 817-829. https://doi.org/10.2138/am-2019-67871593

Putirka, K. D., & Xu, S. (2021). Polluted white dwarfs reveal exotic mantle rock types on  exoplanets in our solar neighborhood. Nature Communications, 12(1), 6168.1595

Tosi, N., Grott, M., Plesa, A. C., & Breuer, D. (2013). Thermochemical evolution of Mercury's interior. Journal of Geophysical Research: Planets, 118(12), 2474-2487.  https://doi.org/10.1002/jgre.201681688

 

 

 

How to cite: Cioria, C., Mitri, G., Connolly, J. A. D., Perrillat, J.-P., and Saracino, F.: Mineralogy of the mantles in sub- Earths and exo- Mercuries, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1011, https://doi.org/10.5194/epsc2024-1011, 2024.

TP11 | Unveiling Venus from atmosphere to core

EPSC2024-836 | ECP | Posters | TP11

Analysis of lava flow features on Venus for radar sounder simulations 

Lisa Molaro and Lorenzo Bruzzone
Thu, 12 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P41

Introduction

Previous missions to Venus depicted an environment dominated by volcanic landforms and hostile atmospheric conditions. The surface was imaged by the Magellan mission, and compositional information were extracted from the Venera, VEGA and Venus Express missions. However, these data show low resolution and uncertain geologic correlation. To improve our knowledge, new missions towards Venus are planned, like ESA’s EnVision. EnVision will study the surface and subsurface of Venus and its relationship with the atmosphere with various instruments, among which a subsurface radar sounder (SRS). SRS will operate at a central frequency of 9 MHz with a bandwidth of 5 MHz, penetrating for few hundred meters through the subsurface, with a vertical resolution of about 20 m [1]. One of SRS targets are lava flow features. This work analyses the existing literature related to lava flow features on Venus, to extract morphometric and compositional information to improve the SRS performance prediction through simulations based on geological analogues. This approach exploits existing radargrams in geologically analogous terrains to produce realistic simulations of the investigated target, using parameters related to the composition and morphometry of the target [2].

Lava flows properties and distribution

Several classifications for Venus flows are based on morphology and radar backscatter (Figure 1). These morphologies are then compared to terrestrial common effusive features: pahoehoe, a'a, and blocky. Pahoehoe interpretation prevails based on Magellan radar backscatter of Venus flows, panoramas of the Venera landing sites and comparison with terrestrial lava flows, with values of rms slopes at Magellan SAR resolution of 75 m ranging from 2.5° to 8° [3]. Estimates of flow thicknesses from Magellan altimetric data and stratigraphic relationships span from 10-30 m to 400 m [4]. The extension of flows ranges from tens up to thousands kilometres.

Figure 1: Classifications of lava flow features (elaboration from [5], [6], [7], [8]).

The composition of lava flows inferred by previous missions and morphological observations suggest a predominantly mafic composition, mostly tholeiitic and alkali basalts (Table 1) [9]. Predicted values of porosity are lower for Venus than Earth, varying from 0.05 to 0.75 (bubble volume fraction). More exotic compositions for channels are considered, such as carbonatite or sulphur, and more evolved compositions are possible. Emissivity measurements of Venus flows range from 0.7 to 0.9, consistent with basaltic samples [3].  

Table 1: Characteristics of Venera and VEGA landing sites from [9].

Landing site

(Lat°, Long°)

Geochemistry

Inferred porosity

Venera 8

(-10.70, 335.24)

Very high K, Th, U

-

Venera 9

(31.01, 291.64)

Low K, U; high Th

-

Venera 10

(15.42, 291.51)

Low K, U, Th

1–7%

Venera 13

(-7.55, 303.69)

High K basalt

50–53 %

Venera 14

(13.05, 310.19)

Tholeiitic basalt

60–62 % (top layer)/ 50–53 % (below)

VEGA 1

(8.10, 175.85)

Low K, U, Th

-

VEGA 2

(-7.14, 177.67)

Tholeiitic basalt; low K, U, Th

13%

 

Lava flows are associated with volcanic edifices, rift zones and coronae. There is a concentration of these features in the Beta-Atla-Themis Regios and in Sedna Planitia (Figure 2) [6]. Comparing the probable composition of lava flows to that of their terrestrial analogues, geodynamic contexts for tholeiitic and alkali basalts would be consistent with Earth analogues like NMORB, islands arcs and hot spots, with melting occurring in the shallow mantle or at much greater depth. Carbonatite volcanism occurs on Earth in intraplate regions, on hotspots, associated with orogenic activity or plate separation, and is mantle-derived [8].

 

Figure 2: Map of Venus with major volcanic environments (elaboration from [10], [11]).

Conclusions

The analysis of the literature on lava flows on Venus is useful for both fine-tuning the expected performance of SRS on these features and defining scenarios to be accurately simulated. Considering lava flows compositional and morphological parameters, SRS is expected to be able to detect changes between flows based on composition, porosity, surface roughness and thus to provide a new stratigraphic perspective of Venus history. 

References

[1]           L. Bruzzone et al., ‘Envision Mission to Venus: Subsurface Radar Sounding’, in IGARSS 2020 - 2020 IEEE International Geoscience and Remote Sensing Symposium, Sep. 2020, pp. 5960–5963. doi: 10.1109/IGARSS39084.2020.9324279.

[2]           S. Thakur and L. Bruzzone, ‘An Approach to the Simulation of Radar Sounder Radargrams Based on Geological Analogs’, IEEE Trans. Geosci. Remote Sens., vol. 57, no. 8, pp. 5266–5284, Aug. 2019, doi: 10.1109/TGRS.2019.2898027.

[3]           B. A. Campbell and D. B. Campbell, ‘Analysis of volcanic surface morphology on Venus from comparison of Arecibo, Magellan, and terrestrial airborne radar data’, J. Geophys. Res. Planets, vol. 97, no. E10, pp. 16293–16314, Oct. 1992, doi: 10.1029/92JE01558.

[4]           K. M. Roberts et al., ‘Mylitta Fluctus, Venus: Rift-related, centralized volcanism and the emplacement of large-volume flow units’, J. Geophys. Res. Planets, vol. 97, no. E10, pp. 15991–16015, 1992, doi: 10.1029/92JE01245.

[5]           M. G. Lancaster et al., ‘Great Lava Flow Fields on Venus’, Icarus, vol. 118, no. 1, pp. 69–86, Nov. 1995, doi: 10.1006/icar.1995.1178.

[6]           J. W. Head et al., ‘Venus volcanism: Classification of volcanic features and structures, associations, and global distribution from Magellan data’, J. Geophys. Res. Planets, vol. 97, no. E8, pp. 13153–13197, 1992, doi: https://doi.org/10.1029/92JE01273.

[7]           E. Stofan, ‘Development of Large Volcanoes on Venus: Constraints from Sif, Gula, and Kunapipi Montes’, Icarus, vol. 152, no. 1, pp. 75–95, Jul. 2001, doi: 10.1006/icar.2001.6633.

[8]           V. R. Baker et al., ‘Channels and valleys on Venus: Preliminary analysis of Magellan data’, J. Geophys. Res. Planets, vol. 97, no. E8, pp. 13421–13444, 1992, doi: https://doi.org/10.1029/92JE00927.

[9]           A. Abdrakhimov and A. Basilevsky, ‘Geology of the Venera and Vega Landing-Site Regions’, Sol. Syst. Res., vol. 36, pp. 136–159, Jan. 2002, doi: 10.1023/A:1015222316518.

[10]         R. M. Hahn and P. K. Byrne, ‘A Morphological and Spatial Analysis of Volcanoes on Venus’, J. Geophys. Res. Planets, vol. 128, no. 4, p. e2023JE007753, 2023, doi: https://doi.org/10.1029/2023JE007753.

[11]         E. R. Stofan et al., ‘Preliminary analysis of an expanded corona database for Venus’, Geophys. Res. Lett., vol. 28, no. 22, pp. 4267–4270, 2001, doi: https://doi.org/10.1029/2001GL013307.

How to cite: Molaro, L. and Bruzzone, L.: Analysis of lava flow features on Venus for radar sounder simulations, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-836, https://doi.org/10.5194/epsc2024-836, 2024.

EPSC2024-1042 | ECP | Posters | TP11

Searching for a near-surface particulate layer using near-IR spacecraft observations 

Shubham Kulkarni, Patrick Irwin, and Colin Wilson
Thu, 12 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P63

The Venera 13 probe recorded spectrophotometric data of Venus’s atmosphere from 62 km down to the surface. While the main dataset was lost, a part of it was reconstructed using the originally published graphic material [1]. The reconstructed downward radiance profiles indicate the presence of a possible near-surface particulate layer (NSPL) in the atmosphere. By fitting these Venera 13 radiance data, NEMESIS [2], a radiative transfer and retrieval tool, is used to retrieve the parameters of the particles forming the NSPL, such as size, abundance, and refractive index. The retrieved layer could have an aeolian or volcanic origin, or it could be formed by volatile transport from the surface. Considering the possible formation mechanisms, it is likely that the NSPL exhibits some form of spatiotemporal variability. On nightside of Venus, it is possible to probe the surface and deep atmosphere (surface up to 15 km altitude) in several near-IR thermal emission windows within the wavelength range of 0.78 to 1.18 μm. Thus, if optical thickness variations of the NSPL exist, then it could be possible to detect them using repeated spacecraft observation in near-IR windows. The instruments EnVision/VenSpec-M and VERITAS/VEM will perform such observations using several surface-observing spectral windows [3, 4]. Motivated by these upcoming missions, simulations are performed to study the effect of the NSPL on near-IR spacecraft observations.

NEMESIS is used to perform simulations of the nightside thermal emission while introducing an assumed variability in the above retrieved NSPL. The sensitivity of the simulated thermal emission spectra to an assumed variability of the NSPL is inspected. However, the main cloud deck (MCD) also shows high variability and affects the surface thermal emission. Thus, the spacecraft observations are simulated by also introducing the variability of MCD in simulations. These simulated spacecraft observations are then corrected for MCD optical thickness variations by following the methodology given by [5] to correct the Venus Express/VIRTIS-M-IR observations. The detectability of the NSPL after correcting the observations for MCD variations is then investigated.

 

References: 

[1] Ignatiev, N. I., Moroz, V. I., Moshkin, B. E., Ekonomov, A. P., Gnedykh, V. I., Grigor’ev, A. V., and Khatuntsev, I. V. Cosmic Research 35(1), 1–14 (1997).

[2] Irwin, P. G., Teanby, N. A., de Kok, R., Fletcher, L. N., Howett, C. J., Tsang, C. C., Wilson, C. F., Calcutt, S. B., Nixon, C. A., and Parrish, P. D. Journal of Quantitative Spectroscopy and Radiative Transfer 109(6), 1136–1150 (2008).

[3] Helbert, J., Vandaele, A. C., Marcq, E., Robert, S., Ryan, C., Guignan, G., Rosas-Ortiz, Y. M., Neefs E., Thomas, I. R., Arnold, G., Peter, G., Widemann, T., and Lara, L. M., In Infrared Remote Sensing and Instrumentation XXVII, 1112804, SPIE, (2019).

[4] Helbert, J., Pertenaïs, M., Walter, I., Peter, G., Säuberlich, T., Cacovean, A., Maturilli, A., Alemanno, G., Zender, B., Arcos Carrasco, C., and others. In Infrared Remote Sensing and Instrumentation XXX, 1223302, SPIE, (2022).

[5] Mueller, N. T., Smrekar, S. E., and Tsang, C. C. Icarus 335(August 2019), 113400 (2020).

How to cite: Kulkarni, S., Irwin, P., and Wilson, C.: Searching for a near-surface particulate layer using near-IR spacecraft observations, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1042, https://doi.org/10.5194/epsc2024-1042, 2024.

EPSC2024-1078 | ECP | Posters | TP11 | OPC: evaluations required

Radiative Transfer Modelling of Venus – A comparison of line databases 

Ankita Das, Nils Mueller, Franz Schreier, David Kappel, John Lee Grenfell, Heike Rauer, and Jörn Helbert
Thu, 12 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P64

Introduction:

The Venusian atmosphere is a fascinating object of interest to planetary scientists. Obtaining in-situ data from Venus’ atmosphere and surface is however challenging. Radiative transfer (RT) modelling is an essential tool to understand planetary atmospheres. In a nutshell, radiative transfer codes model absorption, emission, and scattering of light by various components present in the atmosphere and surface. Accurate modelling of the atmosphere is also essential for decoding surface information from remote sensing data collected by probes. In the coming decade, several missions to Venus are planned that aim to image Venus thermal emission in the NIR spectral windows [1]. In order to process the data from these missions once they are available, radiative transfer modelling of the Venusian atmosphere is a necessary first step. One important aspect of the model is absorption by gases. Modelled absorption cross-sections are governed by the line list chosen for the model. These line lists are provided for wavelength ranges over which absorption occurs. The high-resolution transmission molecular absorption database (HITRAN) is a frequently used line database in radiative transfer modelling [2]. Several Venus atmospheric studies [3], however, have relied on the database of [4] for CO2 lines, referred to as “Hot CO2” from here on. This database is structurally similar to high-temperature molecular spectroscopic database (HITEMP) [5]. Fortunately, due to advances in exoplanetary sciences, newer line databases have been developed for high temperature atmospheres which are yet to be applied to Venusian atmospheric studies. In this work we compare absorption cross sections generated by using different line databases for relevant species present in the Venusian atmosphere: HITRAN 2020, HITEMP 2010, Hot CO2, and ExoMol [2,4,5,6]. Additional comparisons are made between radiance spectra generated by radiative transfer methods using different line lists with measured spectra in the NIR wavelength range.

SPICAV dataset from Venus Express:

The NIR wavelength range of 0.8 – 1.2 micron contains spectral windows where Venus’ surface thermal emission radiation is detectable from space, paving the way for surface studies in these bands [4]. Hence, the NIR region of Venus’ spectra is of particular importance. The Spectroscopy for the Investigation of the Characteristics of the Atmosphere of Venus (SPICAV) suite on board Venus Express made observations of Venus’ nightside in the spectral range of 0.65–1.7 um. The data used in our analysis is detailed in [7].

Materials and Methods:

The following computational tools have been utilized in this work: PYthon for Computational ATmospheric Spectroscopy (Py4CATS) [8] and Planetary Spectrum Generator (PSG) [9]. Both Py4CATS and PSG are equipped to read various line lists to produce radiance spectra from calculated absorption coefficients for a specified atmospheric profile. PSG is an online multi-purpose tool capable of computing radiance and transmission for planets and exoplanets with preconfigured settings. Py4CATS implements several scripts for radiative transfer calculations which can be executed from a Python interpreter or from a console. PSG includes a module capable of modelling scattering in the Venusian clouds while Py4CATS has to be combined with additional tools to do so, e.g. [10].

Preliminary results:

In the Venusian atmosphere, absorption features are majorly dominated by CO2 and H2O lines. Thus, we start by comparing absorption cross sections for CO2 computed by Py4CATS for the following line databases: HITRAN 2020, HITEMP 2010, and Hot CO2. From figure 1, one can note that HITEMP 2010 has missing lines in the 8500 – 9000 cm-1 wavenumber region and that absorption cross sections produced from HITRAN 2020 are much closer to those of Hot CO2. This work aims to further explore implications on radiative transfer modelling using recently updated line databases such as HITRAN 2020 and ExoMol.

Figure 1: Comparison of line databases in NIR region generated using Py4CATS, computed for575 K and 2x106 Pa

To decide which database is best suited it is necessary to compare modelled top of atmosphere radiances to observed spectra. We use the PSG to model nadir radiances of the default Venus atmosphere profile and modules (i.e. Lambertian surface with emissivity 0.8, absorption by gases, Rayleigh scattering, multiple scattering at cloud droplets with cloud optical thickness 15). Figure 2 shows that there are noticeable differences in the radiance spectrum generated by the PSG [9] upon using HITRAN 2020 line database and the ExoMol line database.

Fig 2: Radiance spectra generated by PSG (HITRAN 2020 and ExoMol) compared to SPICAV IR data.

It can be noted from figure 2 that radiances generated using HITRAN 2020 database to have a better agreement with SPICAV data than those generated using ExoMol database. It is likely that further modifications to the line shapes and continuum absorption similar to those that are used in other Venus studies will have to be introduced [3, 11].

Conclusions:

The HITRAN 2020 line list is closer to the often used ‘’Hot CO2’’ line list than previous versions of HITRAN. Line databases such as ExoMol which are intended for high temperature atmospheres may improve radiative transfer models for Venus, although our initial comparison does not show an improvement over HITRAN. Overall our investigations have shown HITRAN 2020 to be more promising for Venus RT studies than HITEMP 2010 and ExoMol. This work aims to further explore applying these new databases to radiative transfer models and compare generated spectra to measured data.

References:

[1] Allen D. A. et al. (1984) Nature, 307, 222–224

[2] Gordon I. E. et al. (2022) J. Quant. Spectrosc. Radiat. Transfer, 277, 107949

[3] Bézard B. et al. (2011) Icarus, 216(1), 173–83

[4] Pollack J. B. et al. (1993) Icarus, 103, 1–42

[5] Rothman L. S. et al. (2010) J. Quant. Spectrosc. & Radiat. Transfer, 111(12-13), 2139–2150

[6] Tennyson J. et al. (2016) J. Mol. Spectrosc., 327, 73 – 94

[7] Korablev O. et al. (2006) J. Geophys. Res. 111(E9)

[8] Schreier F. et al. (2019) Atmosphere, 10(5), 262

[9] Villanueva G. L. et al. (2018) J. Quant. Spectrosc. & Radiat. Transfer, 217, 86 – 104

[10] Efremenko D. et al. (2023) Environmental Sciences Proceedings , 29(1), 20

[11] Kappel D. et al. (2016) Icarus 265, 42–62

How to cite: Das, A., Mueller, N., Schreier, F., Kappel, D., Grenfell, J. L., Rauer, H., and Helbert, J.: Radiative Transfer Modelling of Venus – A comparison of line databases, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1078, https://doi.org/10.5194/epsc2024-1078, 2024.

EPSC2024-1152 | ECP | Posters | TP11

Remote Sensing Investigation of Recent Volcanic Activity on Reykjanes Peninsula, Iceland, as an Analog for Venus 

Akin Domac, Solmaz Adeli, Thomas Kenkmann, Gabriele Arnold, and Jörn Helbert
Thu, 12 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P46

Abstract: Venus and Earth share key similarities yet differ fundamentally. It is uncertain how Venus's evolution diverged from Earth's and if its past was significantly different from its present state. Despite the interest in Venus since the dawn of planetary missions, orbiters could only use radar waves to observe the surface due to its dense and opaque atmosphere. Subsequent research revealed that Venus's atmosphere permits the transmission of electromagnetic waves emitted from its hot surface through six spectral bands centered at 0.86 μm, 0.91 μm, 0.99 μm, 1.02 μm, 1.11 μm and 1.18 μm within five atmospheric windows [1]. Advances in near-infrared imagers and laboratory measurements of the emissivity of rocks have spurred a new era of Venus exploration. In June 2021, two orbiters were announced as part of a fleet of Venus missions: NASA’s VERITAS and ESA’s EnVision. Both orbiters will be equipped with emissivity mappers designed and built by the German Aerospace Center (DLR) and aim to use these six bands to map the surface globally: Venus Emissivity Mapper (VEM) and VenSpec-M, respectively. [2, 3].

Here, we inquire into the potential of near-infrared remote sensing observations to investigate recent or active volcanic activity by employing machine learning methods for classifying surface units and change detection in preparation for future Venus missions.

Iceland serves as a primary analog for Venus in addressing the research questions of this study. This is because Iceland has a young, mostly vegetation-free surface that is formed by ongoing volcanic activity. Reykjanes Peninsula has been experiencing frequent fissure eruptions for the last 4 years. Two of these eruptions, known as the 2021 and 2022 Fagradalsfjall eruptions, resulted in the formation of a 5.01 km2 flow field called Fagradalshraun [4, 5, 6]. The study investigated the Reykjanes Peninsula, focusing on Fagradalshraun, using near-infrared bands in the spectral range of VEM of pre- and post-eruption Sentinel-2A datasets. The imagery dataset was selected based on three essential criteria that were satisfied by ESA’s Sentinel-2 constellation: Open public access, spectral bands within the range of VEM, and temporal coverage of the eruptions.

Supervised classification with the Support Vector Machines method and unsupervised classification with the K-means clustering method were performed to distinguish between geological units based on their age differences. In addition, the contributions of the Texture Analysis and the Principal Components Analysis methods for improving the classification results were tested. Then, the Normalized Difference of Time Series method was used for change detection. The results showed that despite employing a limited number of spectral bands as VEM and VenSpec-M will provide, both classification methods could distinguish between fresh basalt, weathered basalt, and aeolian cover. The accuracy of the classification methods improved with the inclusion of texture measures. Additionally, it was observed that the highest accuracy in unsupervised classification was achieved through the utilization of Principal Component Analysis. Finally, it was seen that the Normalized Difference of Time Series method detected the Fagradalshraun with an accuracy of 95%.

The results presented here underscore the capability of VEM and VenSpec-M, especially in detecting potentially active or recently active volcanism on Venus. In addition, they emphasize the contribution of studying the terrestrial analogs to enhancing the science return from future Venus missions.

References: [1] Helbert, J., et al. (2016). SPIE. [2] Helbert, J., et al. (2020). SPIE. [3] Widemann, T., et al. (2023). Space Science Reviews. [4] Adeli, S., et al. (2023). LPSC. [5] Adeli, S., et al. (2024) LPSC. [6] Pedersen, G., et al. (2022). Geophysical Research Letters.

How to cite: Domac, A., Adeli, S., Kenkmann, T., Arnold, G., and Helbert, J.: Remote Sensing Investigation of Recent Volcanic Activity on Reykjanes Peninsula, Iceland, as an Analog for Venus, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1152, https://doi.org/10.5194/epsc2024-1152, 2024.

TP12 | Lunar Science and Exploration Open Session

EPSC2024-87 | ECP | Posters | TP12 | OPC: evaluations required

In-situ extraction and purification of water on the moon – the LUWEX project 

Noria Brecher, Christopher Kreuzig, Gerwin Meier, Christian Schuckart, and Jürgen Blum
Tue, 10 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P30

For the future of sustainable space exploration, In-Situ Resource Utilization (ISRU) plays a key role. Thus, the investigation and development of new technologies focusing on utilizing local resources is crucial to cost-effective robotic and human space missions with less logistical effort. Water does not only serve as consumable for astronauts and plants but also as rocket propellant by electrolysing it into hydrogen and oxygen.

Therefore, the LUWEX (Validation of Lunar Water Extraction and Purification Technologies for In-Situ Propellant and Consumables Production) project aims to develop and validate a technology demonstrator for the extraction of water from icy-lunar-regolith simulant, the following purification and quality monitoring. Several experiments with different mixing ratios and material compositions will be conducted. The Technology Readiness Level (TRL) is increased for the whole process chain from TRL 2 and 3 to TRL 4 based on the planned experiments. The extraction, capturing and liquefaction subsystems are placed inside the thermal-vacuum chamber (TVAC) of the Comet Physics Laboratory (CoPhyLab) at TU Braunschweig, detailed in Kreuzig et al. 2021, to simulate the relevant lunar environment. The TVAC, including the extraction, capturing and liquefaction subsystems inside it, is shown in Fig.1. This image also shows the purification and quality monitoring subsystems placed on the tray on the right side of the chamber. Two liquid nitrogen cooling systems are utilized to cool the TVAC to an ambient temperature of around 100 K. The working pressure inside the TVAC is around 10-5 mbar. To realize the control of all subsystems inside the TVAC, a complex valve system, including 20 electric-pneumatic valves, is installed around it. Further, an electric-pneumatic slider is mounted between the capturing and liquefaction subsystems inside the TVAC (see Fig.1c). This is necessary as a cold trap consisting of two 3D-printed, actively cooled copper-cones is utilized to capture water vapour extracted from the regolith by heating with the liquefaction chamber mounted below it. Once water ice has deposited onto the cold trap, it is heated to release the ice into the liquefaction chamber. There, the ice is liquefied by heating the liquefaction chamber. To liquefy the ice, the pressure inside the liquefaction subsystem needs to increase to the vapour saturation pressure of water. The chamber is, therefore, separated from the cold trap utilizing the above-mentioned slider.

Besides providing an artificial lunar environment, the CoPhyLab is responsible for the production of icy simulant. Therefore, granular water ice with a mean grain size of 2.4 microns, which is produced with the ice machine outlined in Kreuzig et al. 2021, is mixed with lunar Highland simulant by Lunex Technologies. To prevent any thermal evolution of the water ice, the mixing is performed by adding liquid nitrogen to the materials. In the resulting mixture, the water ice particles are expected to be evenly distributed in the regolith. During the water ice production, impurities such as CO2 or methanol can be added.

The experiments with the LUWEX system will be performed in Summer 2024, and the first results will be presented in September. Firstly, an experiment run with an icy-lunar-regolith sample with around 50% water ice in mass is planned to characterize and test the sample preparation, infilling and the extraction, capturing and liquefaction subsystems. Then, icy-regolith samples with 5%, 10% and 15% water ice in mass will be realised, and the complete process chain, including purification and water quality monitoring, will be tested.

More information on the LUWEX project is available at https://luwex.space.

Figure 1: The CoPhyLab TVAC with the extraction (a), capturing (b) and liquefaction (d) subsystems mounted inside. A slider is installed between the capturing subsystem and the liquefaction (c). On the right side of the TVAC, the purification subsystem is placed with two storage tanks (e). Below the purification sits the quality monitoring subsystem (f). Image credit: Luca Kiewiet

References:

Kreuzig, C. et al. (2021). “The CoPhyLab comet-simulation chamber”. Review of Scientific Instruments 92.11, S.115102.

How to cite: Brecher, N., Kreuzig, C., Meier, G., Schuckart, C., and Blum, J.: In-situ extraction and purification of water on the moon – the LUWEX project, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-87, https://doi.org/10.5194/epsc2024-87, 2024.

EPSC2024-196 | Posters | TP12

Astronaut experiences on the slopes along Apollo EVAs 

Wajiha Iqbal, James W. Head, Carolyn H. van der Bogert, Thomas Frueh, Megan Henriksen, Valentin Bickel, David Kring, Harald Hiesinger, David R. Scott, and Thomas Heyer
Tue, 10 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P45

Introduction: The topography of a landing region is a critical factor in the operational safety and planning of human extravehicular activities (EVAs). Consequently, digital terrain models (DTMs) and derived slope maps are employed to impose constraints on proposed EVA paths with minimal slopes and elevation discrepancies. It is clear that there are technical constraints associated with the use and operation of equipment and tools, such as suits and instrument carts, and astronaut walking trafficibility. While these constraints have been improved over the years, it is also important to consider the observational, physiological, and psychological experiences and constraints of astronauts when planning successful EVAs. This study examined Apollo astronaut reports as well as audio and video recordings to ascertain the experience and performance of the astronauts in relation to the topography and slopes encountered during their EVAs.

Data and Methods: We conducted an analysis of the elevations and slopes along the traverses of the Apollo landing sites, utilizing the Lunar Reconnaissance Orbiter (LRO) Narrow Angle Camera (NAC) derived Digital Terrain Model (DTM) with a 2-m pixel resolution and the related slope map with a 6-m baseline [1,2] (Fig. 1). We extracted topographic profiles of each Apollo landing site traverse in ArcGIS Pro and imported them into Matlab for evaluation. The profiles are simplified and vertically exaggerated by a factor of two (Fig. 1).


Furthermore, we examined NASA's extensive archive of Apollo mission journals, images (Fig. 2) and video libraries [3], from which we extracted discussions among the astronauts about their experiences dealing with challenging terrain.

Distance Perception: The Apollo astronauts reported difficulty perceiving the distance and size of objects [3]. In addition, incorrect estimates of size and distance have been reported by travelers on Earth, which are attributed to “the absence of objects of comparison” [4]. Astronauts' perceptions are highly influenced by lighting conditions. In conditions with low sun angles, which are prevalent during landing due to the lunar morning landing time, shadows appear longer (Fig. 2) and the terrain appears more rugged. [5].

 Slopes Perception: The Center of Lunar Science and Exploration report [6] classifies slopes with inclinations of less than 25° as highly accessible. However, more detailed, quantitative data on the performance of Apollo and Lunokhod provides specific guidelines for use under both Earth-based and lunar conditions [7]. The perception of slope steepness can also be affected by lighting conditions, with boulders on slopes casting particularly long shadows which exaggerate the actual steepness of the slopes. During the Apollo 16 mission, the astronauts standing on the North Ray crater rim (Fig. 2) were unable to see the crater floor. Consequently, they perceived the slope of the wall to be approximately 60°. However, this was in fact less than 30° [3] (Fig. 1).

The lunar surface is characterised by a soft regolith, which presents a challenge for walking on uphill slopes. However, the Moon's gravity is ~1/6th of the Earth, which reduces the risk of astronauts sliding downhill [3]. In their own accounts, astronauts have noted that they relied on their center of gravity while walking downhill [3]. Additionally, Apollo astronauts frequently reported the challenge of traversing slopes due to the limited mobility afforded by carrying their tools and samples. This restriction in movement resulted in a longer time needed to complete a traverse [3]. The lunar roving vehicle was introduced to the Apollo 15, 16, and 17 missions to facilitate the exploration of greater distances and steeper slopes by the astronauts.

Future Missions Planning: The 13 Artemis landing sites were selected on ridges and large crater rims at the South Pole in well-illuminated areas with good Earth visibility, and with potential access to the presumably volatile-rich permanently shadowed regions (PSRs) [7]. However, the landing sites are situated in more complex geological terrains, which may necessitate a greater duration of time at individual stations for the completion of tasks. The astronauts climbed up to 150 m during the Apollo missions and walked on slopes of <15°. In contrast, the Artemis landing sites have elevation differences of ~3400 m. The Apollo astronauts traversed near-equatorial regions on relatively smooth surfaces, but perceived them as rough at low sun angles. Near the South Pole, this effect can be more extreme, limiting visibility and mobility during EVA exploration. It is of the utmost importance to prioritize safety measures [6] and to learn from past successful missions in order to optimize success.

 [1] Scholten F., et. al. (2012). J. Geophy. Res. Planet 117 (E12), 2156–2202.

[2] Henriksen M. R., et. al. (2017). Icarus 283, 122–137.

[3] Jones E. M., Glover K. (2017). Apollo lunar surface journal. NASA.

[4] Darwin, C, (1835) ch. XV, p 347, ISBN-13 : 978-1787377424

[5] van der Bogert, C. H., et. al. (2022). LPSC 53., # 1347.

[6] Kring, D. A., Durda, D. (2012). LPI. #1694.

[7] Basilevsky, A.T., et al. (2019). Solar System Research, 53, 383-398.

[8] McKay, D.S., et al. (1991). Lunar sourcebook, 285-356.

 

How to cite: Iqbal, W., Head, J. W., van der Bogert, C. H., Frueh, T., Henriksen, M., Bickel, V., Kring, D., Hiesinger, H., Scott, D. R., and Heyer, T.: Astronaut experiences on the slopes along Apollo EVAs, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-196, https://doi.org/10.5194/epsc2024-196, 2024.

EPSC2024-299 | ECP | Posters | TP12

A Reverse Monte Carlo Method to Investigate the Topography Interaction with the Lunar Exosphere 

Alexander Smolka, Carlos Saenz Reguero, and Philipp Reiss
Tue, 10 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P32
Introduction
Many Monte Carlo-based models of exospheres of our solar system have been developed in the past, for example, by Hodges et al. (1973), Crider et al. (2002), Grava et al. (2014), Hurley et al. (2017), Killen et al. (2019), Kegerreis (2020), Schörghofer et al. (2021), and Smolka et al. (2023). Their main setup includes tracking individual particles from a defined source of the exosphere to the eventual loss of the particles in a forward manner. Due to the size of the particle trajectories through the exosphere, which can span over more than half the Moon’s circumference, the simulations must geometrically include the entire Moon. While this is an excellent tool for investigating global effects and overarching exospheric dynamics, its accuracy on a significantly smaller scale is lower. The smaller the region of interest and, thus, the relevant domain, the less likely it is for any particles of the global simulation to travel through this volume, which negatively impacts the resolution of the output.
 
There are several such scenarios where a traditional Monte Carlo simulation of an exosphere is not ideally suited for the problem, including investigating local and regional topography and its influence on the exosphere. To tackle this problem, we are implementing a reverse Monte Carlo (RMC) method into existing exosphere simulation frameworks, where, instead of tracking a particle from its source to its loss, the simulation starts the particle inside the relevant domain and evaluates the likelihood of this assumption by backtracking its trajectory. This ensures that every simulated particle contributes to the investigation of the region of interest, which greatly reduces the computational effort and increases the resolution of the output. While this method is already used in other applications like the radiative transfer in exospheres (Gratiy et al., 2010), it has not yet been applied to particle tracking in exosphere models. The following section describes the numerical model and its internal mathematical structure.
 
Numerical Model
The main problem of reverse tracking particle trajectories in exosphere models is the unknown kinetic energy of the particle. A backwards trajectory leads from the known landing position, x1, to the unknown starting position, x0. Generally, the shape of the trajectory, including the angular and speed distribution of the launch velocity, is based on the processes that lead to the release of the particle at x0. Since this information is inaccessible a priori, the reverse model must include a probabilistic approach to the particle’s kinetic energy with a weighing factor that is calculated posteriorly.
 
Mathematically speaking, let the particle start its reverse journey at x1, where it draws a landing velocity vector v1 based on a velocity distribution V1 with probability density p1. The reverse trajectory uses the vector −v1 as its starting velocity and backtracks the particle to its actual launch position x0, where it lands with a velocity of −v0. Based on the conditions and processes occurring at x0, the probability density p0 for drawing v0 from V0 can be calculated.
 
Using Bayes’ theorem with the events A that the particle launched at x0 and B that the trajectory in question occurred, p(B) = p1 is the evidence of the trajectory and p(BA) = p0 is the conditional probability of the trajectory given that a particle was launched at x0, leaving only the prior probability p(A) of a particle being launched at x0 to be determined. Together, one can evaluate the weight w of the reverse trajectory as
If multiple hops of a particle are tracked using the RMC method, each i-th trajectory gets assigned its respective weight wi, where the prior probability of launch p(A)i is equal to the weight of the next reverse trajectory, wi+1. Thus, the i-th weight can be recursively calculated as
with a total of N trajectories and an initial launching probability of p(A)N.
 
The principle is shown in Fig. 1, where from an initial position x1, a particle is tracked backwards using the RMC method. The drawn velocities and probabilities are shown in blue. This process can be repeated for multiple hops of the particle until it leaves the computational domain or has come from a known source with a known probability distribution. While p0 and p1 can be directly calculated based on the underlying process, for example, using Maxwellian velocity distributions for thermally desorbed exospheric particles, the prior probability p(A)N has to be estimated either based the conditions at the first x0 outside of the computational domain. A first-order approximation would assume that p(A)N is linearly dependent on the particle number density at x0.
 
Model Scenarios
The main scenario investigated is the interaction of the lunar exosphere with local topography. First, to validate the numerical method, the model is tested with a perfectly flat surface, and the results are compared with the global results of an exosphere simulation. We use Helium particles to further reduce the complexity of the entire model and compare the results with our global model (Smolka et al., 2023). The next step will include a simplified crater geometry where, based on the crater sizes (depth and diameter), we will investigate the exosphere dynamics above and landing probabilities inside the crater. The latter will shed light on whether particles have a preferred landing position inside the crater, which could be used for further studies of permanently shadowed regions and, in extension, the accumulation of volatile species within.
 
Additionally, the method will be used to study the interaction of instruments like mass spectrometers with the lunar exosphere. Particles can be launched from the instrument and tracked backwards to evaluate the likelihood of reaching the instrument. This can be used to optimize the instrument’s position and orientation and estimate the expected particle flux and the required measurement time.
 
Bibliography
  • Crider et al. (2002). Advances in Space Research
  • Gratiy et al. (2010). Icarus
  • Grava et al. (2014). Icarus
  • Hodges et al. (1973). Fourth Lunar Science Conference
  • Hurley et al. (2017). Icarus
  • Kegerreis (2020). Springer Theses Ser.
  • Killen et al. (2019). Icarus
  • Schörghofer et al. (2021). Space Science Reviews
  • Smolka et al. (2023). Icarus

How to cite: Smolka, A., Saenz Reguero, C., and Reiss, P.: A Reverse Monte Carlo Method to Investigate the Topography Interaction with the Lunar Exosphere, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-299, https://doi.org/10.5194/epsc2024-299, 2024.

EPSC2024-1070 | ECP | Posters | TP12 | OPC: evaluations required

Simulating the 6.5 days surface Artemis 3 mission with a first woman and man: EMMPOL20 4Artemis3 EuroMoonMarsPoland May 2024 Analog Astronaut Campaign
(withdrawn after no-show)

Clara Laforet, Matthew Harvey, Bernard Foing, Agata Kołodziejczyk, and Molly Balfe
Tue, 10 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P46

TP14 | Mars Science and Exploration

EPSC2024-877 | ECP | Posters | TP14 | OPC: evaluations required

Study of Mars mesosphere longitudinal temperature variations from NOMAD-SO onboard ESA’s TGO. 

Loïc Trompet, Lori Neary, Ian Thomas, Shohei Aoki, Frank Daerden, Justin Erwin, Arnaud Mahieux, Séverine Robert, Miguel Ángel López-Valverde, Giuliano Liuzzi, Geronimo Villanueva, Adrián Brines, Giancarlo Bellucci, Manish Patel, and Ann Carine Vandaele
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P60

NOMAD [1] is one of the four instruments on board ESA’s Trace Gas Orbiter and consists of three channels: SO, LNO, and UVIS. The SO channel is dedicated to solar occultation measurements and thus probes the Martian terminator. SO is an infrared spectrometer (2.3-4.3 µm) composed of an echelle grating with an acousto-optic tunable filter for the selection of the diffraction orders. SO has been regularly scanning the atmosphere of Mars from the troposphere to the upper thermosphere since the beginning of the science operations of the Trace Gas Orbiter on April 21, 2018. One diffraction order is ~25 cm-1, and six diffraction orders are scanned at each occultation. The spectral resolution is ~0.15 cm-1, and the signal-to-noise ratio is ~2500. The field of view varies between 1.6 km and 1.85 km, and the vertical sampling varies between 0.1 km and 1 km depending on the beta-angle of TGO. The vertical resolution of the profiles is ~2.5 km. Recently, in 2024, SO started to scan 12 diffraction orders per occultation, dividing the vertical sampling by two but improving the coverage of several species. The resulting vertical resolution is then reduced to a maximum of 50 %. The instrument function of SO was described in ref. [2].

The retrieval of CO2 density and temperature was described in ref. [3]. The radiative transfer computations are performed with the ASIMUT software [4]. The regularisation of the profiles is carefully fine-tuned with an iterated Tikhonov method [3, 5]. This fine-tuning of the regularisation helps to better constrain some variabilities in the profiles that are, for instance, produced by tides and gravity waves. The regularisation does not add any a priori information to the retrieved profiles.

The diffraction order 132 (2966 to 2930 cm-1) is used to infer the CO2 density and temperature in the troposphere (altitudes below ~50 km), while the CO2 density and temperature are inferred in the mesosphere (~50 to ~100 km) from diffraction orders 148 (3326 – 3353 cm-1) and 149 (3348 – 3375 cm-1). Previously, the retrievals in the mesosphere were done only with order 148 (3138 measurements from 2018 to 2023) for the thermosphere, but diffraction order 149 (2880 measurements from 2018 to 2023) is now added to the retrieval pipeline. Diffraction order 148 is sensitive to CO2 density but weekly sensitive to temperature, while diffraction order 149 is highly sensitive to temperature in addition to CO2 density. Diffraction order 165 (3708 – 3738 cm-1) is used to retrieve CO2 density and temperature in the upper thermosphere (140 – 190 km). The lower bounds of the diffraction orders are due to the saturation of the lines[3]. Nevertheless, a full profile combining GEM-Mars [6, 7] and the retrieved profiles from the diffraction orders 132 and 148 (altitudes below 100 km) is provided for the retrievals of other species whose lines are dependent on temperature.

We analysed the longitudinal variations of temperature for some profiles with very close solar longitude, local solar time, and latitude around 60°. Only non-migrating tides (non-synchronous with the relative position of the Sun) can be analysed as the set of profiles corresponds to tight ranges in local times: either 0 h, 9 h, or 15 h. The local times close to 0 h cover the Northern hemisphere in the first half year and the Southern hemisphere in the second half year. The local times close to 9 h and 15 h cover the Southern hemisphere in the first half year and the Northern hemisphere in the second half year. Some preliminary results concerning four sets of profiles in Martian year 35 where longitudinal variations could be inferred were presented in [8]. This analysis was now extended to more than fifty of those sets of profiles from Martian years 34 to 37. Amongst the results, we found an important wavenumber-1 structure at 9 h around LS 60° and 110° in the Southern hemisphere with very similar amplitude and phase for Martian years 35 to 37. Still, we found no structure higher than 1% of the background temperature at 15 h around LS 85° in the Southern hemisphere.

Comparisons to simulations from a GCM [6, 7] show some large differences in the amplitude of longitudinal variations in the mesosphere, especially closer to aphelion in the Southern hemisphere, showing that the dynamical processes occurring at that time might still need to be better constrained. Comparisons to the results of measurements from other instruments are ongoing to confirm those results obtained with SO.

How to cite: Trompet, L., Neary, L., Thomas, I., Aoki, S., Daerden, F., Erwin, J., Mahieux, A., Robert, S., López-Valverde, M. Á., Liuzzi, G., Villanueva, G., Brines, A., Bellucci, G., Patel, M., and Vandaele, A. C.: Study of Mars mesosphere longitudinal temperature variations from NOMAD-SO onboard ESA’s TGO., Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-877, https://doi.org/10.5194/epsc2024-877, 2024.

EPSC2024-1334 | ECP | Posters | TP14

Spatio-Temporal Detection, Aggregation and Tracking of Martian Large-Scale Dust Events 

Timoté Lombard and Luca Montabone
Wed, 11 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P67

Introduction: Martian dust storms play a crucial role in the red planet's climate and weather patterns, affecting both atmosphere and surface conditions over various time scales [1]. Orbital data help track these storms, which pose challenges for spacecraft operations and energy production [2], [3], [4] and [5]. Although monitoring is precursor, accurate forecasting is essential for future Mars exploration [6]. Supervised machine learning (SML) has already been used for automatic detection of Martian dust storms from visible images [7], [8]. However, we are not aware of SML and/or unsupervised machine learning (UML) being utilized for this application with other data sources.

Data: Montabone et al. (2015) have already identified a certain degree of interannual repeatability of large-scale dust storms from zonal mean of column dust optical depth normalized to the reference 610 Pa pressure level (CDOD@610) [9]. A novel approach relying exclusively on UML (see Figure 1) has been developed to detect, aggregate and track Martian large-scale dust events both in space and time (ST-DATMADE), from the publicly-available, multiannual (from Martian Year (MY) 24 to 36) CDOD@610 dataset described in Montabone et al. (2015, 2020) [9], [10]. Large-scale dust events are characterized by surface areas ≥1.6×106 km2 and last more than two Martian days [2]. Normalized CDOD is used to remove topographic features (plains, volcanoes etc.) which affect the distribution of dust within the column. The dataset contains gridded observation maps and kriged maps from CDOD satellite retrievals. Being spatially interpolated, regularly kriged maps (see Figure 2, 1st column) are more relevant to use for SML application compared to irregularly gridded maps.

Detection: Detection of dust events involves the identification of spatio-temporal anomalies in CDOD maps. A dust event is defined as a temporal series of dust episodes characterized by foreground dust anomalies. A foreground dust anomaly refers to a significant increase in CDOD observed during a given Sol (Martian day), regardless of spatial dimension, relative to the background dust level. Background dust refers to the typical or baseline level of diffused dust present in the Martian atmosphere. Therefore, to facilitate the recognition of these anomalies, a “background/foreground” segmentation is initially performed to enhance subsequent data partitioning (see Figure 1). The CDOD data distribution is adjusted by mapping the lowest data associated with background dust into a threshold value, , calculated as follows:

The low dust loading season (LDLS) is a multiannual period between LS (Solar Longitude) = 10° and LS = 140° where there is almost no large-scale dust injection [11]. Here, the LDLS background, , is defined as the 95th percentile () of CDOD during LDLS between MY24 and MY36. It represents the multi-annual atmospheric background during the “clear” season. Additionally, the daily background, , is defined as the 33rd percentile () of CDOD at a given Sol and MY. It considers that the background dust may evolve as dust events may increase the surrounding dust opacity.

Furthermore, the adjusted dataset distribution is extremely and positively skewed, in such case an inverse transform may help to reduce the skewness, making the distribution more symmetrical and evenly distributed, particularly when addressing extreme dust events, [12]. This brings extreme values closer to average values. Within an episode of foreground dust anomaly, distinct features are identified as dust core and dust cloud. A dust core is a value-based partition where atmospheric dust injection is likely. A dust cloud is a value-based partition where atmospheric dust diffusion is likely. Grouping between background dust, dust cloud and dust core is realized through CDOD values partitioning in k=3 partitions using a widely-known UML algorithm: k-means [13]. See Figure 1 and Figure 2, 2nd column.

Aggregation: Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm ([14]) is used for spatial clustering of the partition with the highest CDOD values (i.e. dust core). A dust instance is an individual, isolated spatial-based cluster within the dust core partition, produced by DBSCAN (see Figure 1 and Figure 2, 3rd column).

Tracking: A series of space-coherent and time-continuous dust instances is organized into a dust sequence. To track instances from one Sol to another (Sol-to-Sol) the pairwise centroid distances of consecutive Sols is compared to a “sphere of membership” (SOM). As well as the percentage of overlapping (OLP) of instances from consecutive Sols. Thereby, it is possible to track dust event instances Sol-to-Sol and to deduce whether they form sequences (see Figure 1 and Figure 2, 4th column).

Catalog: Previous section of tracking provides Sol-to-Sol cluster assignments (instances organized in sequences) allowing to build a catalog of dust event instances tagged with an “id” with format: MY.._Sol_..._1,2,3 etc. Space-consistent and time-continuous instances are labelled with a sequence “id” as: MY.._A,B,C, etc. (see Table 1).

Statistics: This catalog will allow a comprehensive analysis of large-scale dust events based on CDOD data. In particular, their spatio-temporal distribution through trajectories, region of origination, LS of origination etc. as well as intensity distribution through area and CDOD statistical features (mean, median etc.).

Acknowledgments: The authors would like to acknowledge the use of the publicly available dataset on the LMD webpage: http://www-mars.lmd.jussieu.fr/mars/dust_climatology/index.html and the NASA PDS webpage:

https://pdsatmospheres.nmsu.edu/data_and_services/atmospheres_data/MARS/montabone.html.

References:

[1]  L. Montabone and F. Forget, 2018, http://hdl.handle.net/2346/74226

[2]  B. A. Cantor et al., 2001, https://doi.org/10.1029/2000JE001310

[3]  H. Wang and M. I. Richardson, 2015, https://doi.org/10.1016/j.icarus.2013.10.033

[4]  M. Battalio and H. Wang, 2021, https://doi.org/10.1016/j.icarus.2020.114059

[5]  C. Gebhardt et al., 2022,

https://www-mars.lmd.jussieu.fr/paris2022/abstracts/poster_Gebhardt_Claus.pdf

[6]  L. Montabone et al., 2022,

https://www-mars.lmd.jussieu.fr/paris2022/abstracts/oral_Montabone_Luca.pdf

[7]  R. Alshehhi and C. Gebhardt, 2022, https://doi.org/10.1186/s40645-021-00464-1

[8]  K. Ogohara and R. Gichu, 2022, https://doi.org/10.1016/j.cageo.2022.105043

[9]  L. Montabone et al., 2015, https://doi.org/10.1016/j.icarus.2014.12.034

[10] L. Montabone et al., 2020, https://doi.org/10.1029/2019JE006111

[11] F. Forget and L. Montabone, 2017, http://hdl.handle.net/2346/72982

[12] B. G. Tabachnick and L. S. Fidell, Using Multivariate Statistics. Pearson Education, 2013.

[13] S. Lloyd, 1982, https://doi.org/10.1109/TIT.1982.1056489

[14] M. Ester et al., 1996, https://www2.cs.uh.edu/~ceick/7363/Papers/dbscan.pdf 

How to cite: Lombard, T. and Montabone, L.: Spatio-Temporal Detection, Aggregation and Tracking of Martian Large-Scale Dust Events, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1334, https://doi.org/10.5194/epsc2024-1334, 2024.

OPS1 | Broadening our understanding of Jupiter’s icy moons and their environment

EPSC2024-17 | Posters | OPS1 | OPC: evaluations required

Global morphology of ENA emissions from the atmosphere-magnetosphere interactions at Europa and Callisto 

C. Michael Haynes, Tyler Tippens, Peter Addison, Lucas Liuzzo, Andrew R. Poppe, and Sven Simon
Mon, 09 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P23
We analyze the emission of energetic neutral atom (ENA) flux from charge exchange between Jovian magnetospheric ions and the atmospheres of Callisto and Europa. For this purpose, we combine the draped electromagnetic fields from a hybrid plasma model with a particle tracing tool for the energetic ions. We determine the ENA flux through a spherical detector that encompasses the entirety of each moon's atmosphere, thereby capturing the complete physics imprinted in these emission patterns. In order to constrain the modifications to the ENA emissions that arise from the periodic change of the ambient plasma conditions, we calculate the emission morphology at multiple positions during a Jovian synodic rotation. To isolate the influence of field line draping, we compare to the emission patterns in uniform fields. Our major results are:

(a) At Europa and Callisto, the majority of detectable ENA emissions are concentrated into a band normal to the Jovian magnetospheric field. (b) The fraction of observable ENA flux that contributes to this band depends on the number of complete gyrations that the parent ions can complete within the moon's atmosphere. (c) Field line draping partially deflects impinging parent ions around both moons, thereby attenuating the ENA flux and driving significant morphological changes to the emission patterns. (d) The band of elevated ENA flux contains a local maximum and a local minimum in intensity, on opposite sides of each moon. At Europa, detectable ENA emissions are maximized slightly west of the ramside apex. At Callisto, they maximize near the Jupiter-facing apex.

How to cite: Haynes, C. M., Tippens, T., Addison, P., Liuzzo, L., R. Poppe, A., and Simon, S.: Global morphology of ENA emissions from the atmosphere-magnetosphere interactions at Europa and Callisto, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-17, https://doi.org/10.5194/epsc2024-17, 2024.

EPSC2024-75 | ECP | Posters | OPS1 | OPC: evaluations required

Non-thermal atmospheric escape of sulfur and oxygen on Io driven by photochemistry and atmospheric sputtering 

Xu Huang, Hao Gu, and Jun Cui
Mon, 09 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P31

Introduction

Io, the innermost of Jupiter’s four Galilean satellites, processes a unique SO2-dominated atmosphere comprising sulfides, oxides (such as SO, S, O2, and O), and other minor alkali and chlorine compounds [1]. In general, three mechanisms are responsible for the generation of an atmosphere on Io, including active volcanism, frost sublimation, and surface sputtering, of which the former two are more important [2]. Both the active volcanism arising from Jupiter's powerful tidal forces and the SO2 frost sublimation from Io's surface release large amounts of gases to replenish the tenuous atmosphere of Io, triggering a rich and complicated photochemical network, which may be a significant source of photochemical escape on Io [3-4]. Meanwhile, Io suffers from intense ion bombardment from Jupiter’s magnetosphere [5]. This constant atmospheric erosion by energetic ion precipitation, referred to as atmospheric sputtering, also serves as an important mechanism of Io’s atmospheric escape.

Aims

With the aid of constantly accumulated understandings of Io’s space environment and atmospheric photochemistry, as well as the updated laboratory measurements [6], we evaluate the non-thermal escape of sulfur and oxygen on Io driven by both photochemistry and atmospheric sputtering [7]. A comprehensive review of the atmospheric escape process on Io is also provided.

Methods

The sputtering yield and escape probability are introduced to evaluate the escape intensity driven by the above two mechanisms. A one-dimensional Test Particle Monte Carlo (TPMC) Monte Carlo model is constructed to track the energy degradation of incident energetic ions and atmospheric recoils from which the sputtering yields and escape probabilities of different atmospheric species are determined. Different plasma populations (S++ and O+) and atmospheric conditions are compared, including high-density volcanic and low-density quiet atmospheric states, in which various chemical channels (photodissociation, neutral-neutral, ion-neutral, and dissociative recombination reactions) are considered. The background atmosphere and ionosphere are adapted from previous photochemical models of [3] and [4].

Results and Conclusions

Our calculations suggest a total escape rate of 3×1029 atom s−1 driven by atmospheric sputtering on Io, and SO2 is the dominant sputtered species. The photochemical escape rates are (1.1−2.0) × 1027 s−1 for total O and (1.5−6.7) × 1026 s−1 for total S, occurring mainly in the atomic form. Further investigations reveal that (1) S++ is the most efficient species for atmospheric sputtering on Io and sputtering yields increase substantially with increasing incident ion mass, energy, and incidence angle; (2) The photochemical escape rates vary extensively with the atmospheric conditions, especially in terms of the intensity of volcanic eruption, resulting in the chemical escape rate increases by up to a factor of five. Photochemistry is the most chemical escape channel. (3) By comparing multiple escape mechanisms including thermal escape (Jeans escape) and non-thermal escape, we conclude that atmospheric sputtering is the dominant mechanism driving atmosphere escape at Io. Photochemical escape outweighs Jeans escape for both atomic O and S for the quiet atmosphere scenario, while for the volcanic scenario, it is likely important for atomic S only.

 

Reference:

[1] Giono, G., & Roth, L. (2021). Io's SO2 atmosphere from HST Lyman-α images: 1997 to 2018. Icarus, 359, 114212.

[2] de Pater, I., Goldstein, D., & Lellouch, E. (2023). The Plumes and Atmosphere of Io. Io: A New View of Jupiter’s Moon, 233-290.

[3] Summers, M. E., & Strobel, D. F. (1996). Photochemistry and vertical transport in Io's atmosphere and ionosphere. Icarus, 120(2), 290-316.

[4] Moses, J. I., Zolotov, M. Y., & Fegley Jr, B. (2002). Photochemistry of a volcanically driven atmosphere on Io: Sulfur and oxygen species from a Pele-type eruption. Icarus, 156(1), 76-106.

[5] Crary, F. J., & Bagenal, F. (1997). Coupling the plasma interaction at Io to Jupiter. Geophysical Research Letters, 24(17), 2135-2138.

[6] Bagenal, F., & Dols, V. (2020). The space environment of Io and Europa. Journal of Geophysical Research: Space Physics, 125(5), e2019JA027485.

[7] Huang, X., Gu, H., Cui, J., Sun, M., & Ni, Y. (2023). Non-Thermal Escape of Sulfur and Oxygen on Io Driven by Photochemistry. Journal of Geophysical Research: Planets, 128(9), e2023JE007811.

How to cite: Huang, X., Gu, H., and Cui, J.: Non-thermal atmospheric escape of sulfur and oxygen on Io driven by photochemistry and atmospheric sputtering, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-75, https://doi.org/10.5194/epsc2024-75, 2024.

EPSC2024-323 | ECP | Posters | OPS1 | OPC: evaluations required

Exploring Europa and Ganymede's Internal Structure: A Statistical Perspective  

Terézia Košíková and Marie Běhounková
Mon, 09 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P28

Introduction 

Past missions, such as Galileo, Cassini and Juno, have significantly advanced our understanding of icy moons. Through their measurements, these missions unveiled the presence of subsurface water reservoirs. Building on these discoveries, upcoming missions, JUICE and Europa Clipper, hold the potential to further enhance our understanding of the hydrospheres by delivering new and improved data, including tidal deformation measurements (Cappuccio et al., 2020; Cappuccio et al., 2022; Mazarico et al., 2023).  

Our study aims to characterize the hydrospheres of Europa and Ganymede. Europa, a smaller satellite, possesses a less extensive hydrosphere, where internal pressures are not expected to reach the levels required for the formation of high-pressure ice phases. In contrast, Ganymede, a larger and differentiated satellite, is known to harbour high-pressure ice phases within its hydrosphere. To evaluate these differences, we employed a combination of known thermodynamic properties (Choukroun and Grasset, 2010; Mcdougall and Barker., 2011; Journaux et al., 2020) alongside satellites’ parameters to assess their plausible internal structures. Additionally, we include anticipated Love number measurements from upcoming missions such as JUICE and Europa Clipper into our statistical analysis, aiming to enhance our understanding of the moons' structural characteristics. 

Model 

We employed the PlanetProfile (Styczinski et al., 2023; Vance et al., 2018) to analyse the internal structure, evaluating 1D models of structure based on fundamental planetary properties such as mass (M), moment of inertia (MOI), and radius (R). 
Additionally, we coupled the PlanetProfile with the Markov chain Monte Carlo (MCMC, Foreman-Mackey et al., 2013) method to assess the statistical properties of the interior structure.  Furthermore, to obtain Love numbers, necessary for determining tidal deformation, we will use own library, based on Sabadini and Vermeersen (2004).   

Results 

Europa 

We assume that Europa is fully differentiated into the hydrosphere, mantle, and core to assess the structure. We analyzed the internal structure of two different oceanic compositions: Seawater and a solution containing MgSO4 salt. For Seawater, our statistical analysis includes variables such as M, R, and MOI assuming normal distribution and the temperature at the interface Ih-ocean interface, Tb, between 249K and 272.5K, assuming uniform distribution. For MgSO4 solution, salinity was incorporated as an additional variable, constrained within the range of 1-10 wt%, assuming uniform distribution. 

The results for Seawater is depicted in Figure 1. For MgSO4 solution, findings are presented in Figure 2. As expected, our observations confirm the composition of Europa's hydrosphere consists of Ih ice and ocean.   In the case of Seawater, despite assuming a uniform distribution of Tb, the actual distribution appears to be non-uniform, likely due to the absence of a corresponding model for the given parameters (R, M, MOI, Tb). However, Figure 2 illustrates that the introduction of salinity allows for lower values of Tb. 

Ganymede  

To model internal structure of Ganymede, we used same procedure as in the case of Europe. We also assumed that Ganymede is fully differentiated. However, the hydrosphere is further divided into layers of Ih ice, liquid ocean and layers of high-pressure ice phases. 

In the case of Ganymede, we only worked with the composition of the MgSO4 salt ocean, whose salinity ranged within the same values as in the case of Europa.  The results for Ganymede with MgSO4 solution are presented in Figure 3. By comparing with Figure 2, we observe a greater thickness of the ocean, as was expected. 

Summary 

We modelled the internal structure of the icy moons Europa and Ganymede, for the cases of Seawater composition for Europa, and MgSO4 solution for both Europa and Ganymede. The PlanetProfile was used and modified by implementing Markov chains, using Emcee library. Europa's hydrosphere is composed of a layer of Ih ice and ocean, and the hydrosphere of Ganymede is composed of layers of Ih ice, ocean, and high-pressure phases of ice, which agrees with theoretical assumptions. Furthermore, we integrate anticipated Love number measurements to reduce the uncertainty in the determination of the internal structure. 

Acknowledgement 

This study received funding from project SVV 260709. 

References 

Cappuccio et al. (2020). Ganymede's gravity, tides and rotational state from JUICE's 3 GM experiment simulation. Planetary and Space Science. 187. 104902. 10.1016/j.pss.2020.104902. 

Cappuccio et al. (2022). Callisto and Europa Gravity Measurements from JUICE 3GM Experiment Simulation. Planet. Sci. J. 3 199. 10.3847/PSJ/ac83c4 

Choukroun, Grasset. (2010). Thermodynamic data and modeling of the water and ammonia-water phase diagrams up to 2.2 GPa for planetary geophysics. The Journal of chemical physics. 133. 144502. 10.1063/1.3487520. 

Foreman-Mackey et al. (2013) emcee: The MCMC Hammer, Publications of the Astronomical Society of the Pacific 125(925), p. 306. https://doi.org/10.1086/670067 

Foreman-Mackey (2016). corner.py: Scatterplot matrices in Python. The Journal of Open Source Software, 1(2), 24. doi: 10.21105/joss.00024 

Gomez Casajus et al. (2021). Updated Europa gravity field and interior structure from a reanalysis of Galileo tracking data. Icarus 358, 114187. https://doi.org/10.1016/j.icarus.2020.114187. 

Journaux et al. (2020). Holistic approach for studying planetary hydrospheres: Gibbs representation of ices thermodynamics, elasticity, and the water phase diagram to 2,300 MPa. Journal of Geophysical Research: Planets, 125, e2019JE006176. https://doi.org/10.1029/2019JE006176 

Mazarico et al. The Europa Clipper Gravity and Radio Science Investigation. Space Sci Rev 219, 30 (2023). https://doi.org/10.1007/s11214-023-00972-0 

McDougall and Barker (2011). Getting started with TEOS-10 and the Gibbs seawater (GSW) oceanographic toolbox. SCOR/IAPSO WG, 127, 1–28. 

Sabadini and Vermeersen (2004). Global Dynamics of the Earth: Applications of Normal Mode Relaxation Theory to Solid-Earth Geophysics. Kluwer Academic Publishers. 

Schubert et al., Interior composition, structure and dynamics of the Galilean satellites, in Jupiter. The Planet, Satellites and Magnetosphere, vol. 1, 2004, pp. 281–306. 

Styczinski et al. (2023). PlanetProfile: Self-consistent interior structure modeling for ocean worlds and rocky dwarf planets in Python. Earth and Space Science, 10, e2022EA002748. https://doi.org/10.1029/2022EA002748 

Vance et al. (2018). Geophysical investigations of habitability in ice-covered ocean worlds. Journal of Geophysical Research: Planets Planets, 123, 180–205. https://doi.org/10.1002/2017JE005341 

 

 

 

 

 


 

How to cite: Košíková, T. and Běhounková, M.: Exploring Europa and Ganymede's Internal Structure: A Statistical Perspective , Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-323, https://doi.org/10.5194/epsc2024-323, 2024.

EPSC2024-933 | ECP | Posters | OPS1 | OPC: evaluations required

Strike-slip tectonics as a precursor to crustal spreading in Anshar Sulcus, Ganymede: Implications for grooved terrain formation 

Mafalda Ianiri, Giuseppe Mitri, Davide Sulcanese, Gianluca Chiarolanza, and Camilla Cioria
Mon, 09 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P35

Introduction

The formation of grooved terrains on the icy surface of Ganymede is still debated and it could involve extensive rifting [1, 2] or spreading [3], strike-slip tectonics [4, 5], and minor cryovolcanic resurfacing [6].

To investigate the origin and evolution of the grooved terrains, we performed a geomorphological and structural analysis of Anshar Sulcus, a grooved terrain located in the anti-Jovian hemisphere within the dark terrain of Marius Regio. In addition, we conducted a topographical analysis producing a Digital Elevation Model (DEM) of this region. We selected Anshar Sulcus because it is a terminal portion of a grooved system, and it was imaged at high resolution by the Galileo SSI. These two aspects helped us to perform a palinspastic reconstruction of the area surrounding the sulcus and, consequently, a reconstruction of the different stages of formation of this sulcus.

 

Methods

The geological map was produced based on SSI Galileo images having a spatial resolution of about 152 m/px, at 1:500000 scale, using the differences in tones, textures and patterns and DEM. To import these images to the Geographic Information System (GIS), we calibrated, filtered and geo-referenced them through the Integrated Software for Imagers and Spectrometers (ISIS4) [7]. The DEM was produced using the “shape-from-shading” tool (SfS) provided by the NASA Ames Stereo Pipeline tool suite [8], maintaining the same spatial resolution as the images.

 

Results

Our geological map (Fig. 1) shows that most of the study area is covered by the dark cratered unit (dc) and consists of a heavily cratered surface with several patches of hummocky material and a pervasive fracturing that we have divided into three main sets. The light grooved unit (lg) consists of a prominent lane that crosscuts the dc characterized by sets of sub-parallel linear grooves.  Additionally, topographic profiles traced perpendicular to the sulcus, revealed that the elevation of the grooved terrain increases towards its central part.

From a morphological observation, we identified well preserved rims of cut craters and fractures along the boundary between the lg and dc. These structures interrupt at the boundary with lg, and this allowed us to reconstruct the possible original position of these structures before the formation of the light unit.

Starting from the westernmost crater along the northern boundary of the dc region (Fig. 2a, 2b), we observed morphological coherence and continuity between the fractures within the crater and the set of fractures in the southern region (set a), suggesting a right lateral displacement of about 15 km. Proceeding to east along the boundaries of the two dc regions, we noted the occurrence of two rims of a possible crater (Fig. 2c, 2d) characterized by a lateral displacement of about 13 km. Moreover, in support of morphological evidence, the topographic profiles traced within the analysed structures show a consistency between the elevations of their southern and northern portions (Figs. 3 and 4).

 

Discussion

Through a combination of structural and topographic analyses, we have found evidence that the formation of Anshar Sulcus grooved terrain was a result of spreading and upwelling of new material. This is supported by the presence of well-preserved structures along the boundaries between dc and lg, and by the higher topographic elevation at the centre of the sulcus, which rules out the possibility of rifting as the cause of its formation.

Our reconstruction shows that the first tectonic event was a right-lateral movement in the NW-SE direction (Fig. 5b) that divided the dark terrain of Marius Regio into two distinct regions with conservative margins. The second tectonic stage (Fig. 5c) presents the beginning of the spreading event that caused the separation of the two regions previously formed as consequence of the strike-slip event, corresponding to the formation of light grooved unit. The right lateral strike-slip tectonics continued during the spreading stage, so the separation of the two regions occurred with a transtensional movement toward a WSW- ENE direction.

Furthermore, we have compiled a chronostratigraphic chart (Fig. 6) that identifies three stratigraphic phases for the formation of Anshar Sulcus and surrounding area. The first phase is characterized by the formation of the dark unit due to contamination of the icy shell by impactors. The second phase is characterized by the formation of the three sets of fractures and by the start of the strike-slip event that produced the division of the dark terrain of Marius Regio. Finally, the third phase is characterized by the formation of the lg unit through crustal spreading. The formation of the grooves inside the light unit occurred relatively quickly, due to the brittle faulting and tilting of the new crust during the extension event that have formed the light unit, explaining the absence of deformed or displaced craters inside the sulcus.

 

Conclusions

Our analysis reveals a two-stage tectonic genesis of Anshar Sulcus, with a first strike-slip stage dividing the dark terrain of Marius Regio into two distinct regions, and a second stage characterized by the formation of the grooved terrain of Anshar Sulcus through crustal spreading and upwelling of new, uncontaminated material, with the simultaneous formation of grooves within it due to brittle fracturing.

 

Acknowledgements

M.I. and G.M. acknowledge support from the Italian Space Agency (2023-6- HH.0).

 

References

[1] Collins et al. (1998b)  Icarus, 135, 345–359.

[2] Prockter et al. (2002) J. Geophys. Res.: Planet. 105, 22519-22540.

[3] Pizzi et al. (2017) Icarus, 288, 148-159.

[4] Pappalardo et al. (1998) Icarus, 135, 276-302.

[5] Cameron et al. (2018) Icarus, 315, 92-114.

[6] Showman et al. (2004) Icarus, 172, 625-640.

[7] Houck and Denicola (2000) ADASS IX, 216, 591.

[8] Beyer et al. (2018) Earth and Space Sci. 5(9), 537–548.

How to cite: Ianiri, M., Mitri, G., Sulcanese, D., Chiarolanza, G., and Cioria, C.: Strike-slip tectonics as a precursor to crustal spreading in Anshar Sulcus, Ganymede: Implications for grooved terrain formation, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-933, https://doi.org/10.5194/epsc2024-933, 2024.

EPSC2024-1021 | ECP | Posters | OPS1 | OPC: evaluations required

Rayleigh-Bénard convection in the subsurface ocean of Ganymede.  

Silvia Pagnoscin, Antonello Provenzale, Jost Von Hardenberg, and John Robert Brucato
Mon, 09 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P26

Introduction: Jupiter's icy moon Ganymede is widely known not only for its large dimensions, overcoming those of Mercury, but rather for its possible subsurface ocean that could be the greatest liquid water reservoir of the entire Solar System covering about the 50% of the whole satellite. The presence of liquid water is one of the main necessary conditions for life as we know making Ganymede one of the main targets for Solar System exploration [1,2]. The ocean that seems to be more than 100 km deep is confined by two ice layers with the upper one made by ice I and the lower one made by high pressure ice V or VI. In addition, Ganymede owes several characteristics such as a magnetic field, a complex internal structure, and a geologically interesting surface that make the satellite a primary target for the geosciences. It is thus important to investigate the possible geophysical and fluid dynamical processes ongoing in Ganymede’s putative ocean, such as convective motions that can lead to interactions between different interior layers and ocean mixing.

Methods:  Icy moons’ oceans such that of Ganymede are known to be heated from below facilitating convection that could transport enough heat to melt the upper ice layer. In this perspective, we explored the (expectedly turbulent) convective dynamics of a portion of the hidden ocean. We considered a Newtonian fluid layer, set in a 3D box with thickness D and horizontal sides 2πD, subject to Rayleigh-Bénard convection. Periodic boundary conditions are used on the vertical boundaries at x=0, 2πD, and y=0,2πD while boundary conditions on temperature, salinity, and velocity on the horizontal boundaries are set with proper equations according to the case study. Gravity is aligned with the vertical direction, and it points opposite to the z axis. The model includes rotation that can be set at any direction in the y-z plane. The fluid density is a function of temperature and salinity; here we first use the Boussinesq approximation so that density becomes independent of pressure, and subsequently we consider a situation where the unperturbed state is compressible while density perturbations are kept incompressible. Once the theoretical model has been defined, simulations are carried out using RBSolve [3,4], a 3D Navier-Stokes Fortran code in the Boussinesq approximation.

Aims:  This study will lead to a better understanding of planetary geophysics and comparative oceanography. In addition, investigating subsurface ocean circulation and its effects on the upper and lower layers of ice will lead us to astrobiological constraints on how a possible Ganymede’s environment could work.

References:

  • [1] Vance S. D. et al. (2014).  Ganymede’s internal structure including thermodynamics of magnesium sulfate oceans in contact with ice. Planetary and Space Science, Vol. 96, Iss. 06. https://doi.org/10.1016/j.pss.2014.03.011
  • [2] Vance S. D. et al. (2018). Geophysical investigations of habitability in ice-covered ocean worlds. J. Geophys. Res. Planets 123, 180–205. https://doi.org/10.1002/2017JE005341 
  • [3] Von Hardenberg J. (2008). Large-scale patterns in Rayleigh–Bénard convection, Physics Letters A, Vol. 372, Iss. 13. https://doi.org/10.1016/j.physleta.2007.10.099 
  • [4] Novi L. et al. (2019). Rapidly rotating Rayleigh-Bénard convection with a tilted axis, Physical Review Vol 99, Iss. 5. https://link.aps.org/doi/10.1103/PhysRevE.99.053116

 

How to cite: Pagnoscin, S., Provenzale, A., Von Hardenberg, J., and Brucato, J. R.: Rayleigh-Bénard convection in the subsurface ocean of Ganymede. , Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1021, https://doi.org/10.5194/epsc2024-1021, 2024.

EPSC2024-1100 | ECP | Posters | OPS1 | OPC: evaluations required

Hiding in Plain Sight: Searching for Evidence of Subduction on Europa's Icy Shell 

Hyunseong Kim, Antoniette Grima, and Luke Daly
Mon, 09 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P36

Satellite images from Voyager 1/2 and Galileo indicate that the Jovian moon Europa has a geologically young surface (40-90 Myr) and might host active tectonics within its icy shell. Previous work has interpreted many of the surface features on Europa as evidence of extensional deformation. However, outside the relatively smooth and asymmetric subsumption bands, evidence of compressional topography is very limited. This suggests that compression induced topographic uplift on Europa must be either; a) very diffuse, potentially due to the elastic properties of the ice, b) undetectable in current satellite images due to photoclinometry and resolution limitations, or c) some of the ice mass must subduct below the surface. To investigate this hypothesis, we first calculate the total volume of new ice that is generated at extension bands and rifts for a first order approximation of the expected amount of compressional uplift, assuming icy shell mass conservation and considering isostatic balance. Using the finite element code ASPECT, we will then run visco- elastic-plastic numerical models of subduction to investigate whether any ‘missing’ topographic signal can result from the subduction of ice and its associated (diffuse) compressional deformation at subsumption bands. Our results have the potential to unravel the mystery of Europa’s topography and provide new insights into the tectonics of icy planetary bodies.

How to cite: Kim, H., Grima, A., and Daly, L.: Hiding in Plain Sight: Searching for Evidence of Subduction on Europa's Icy Shell, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-1100, https://doi.org/10.5194/epsc2024-1100, 2024.

OPS2 | Jupiter and Giant Planet Systems: Juno Results

EPSC2024-496 | ECP | Posters | OPS2 | OPC: evaluations required

A long-term study of the Jovian equatorial atmosphere at the upper troposphere-lower stratosphere from HST observations in the 890-nm methane absorption band 

Mikel Sanchez Arregui, Arrate Antuñano, Ricardo Hueso, and Agustin Sanchez-Lavega
Tue, 10 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P52

Despite the obvious differences among Jupiter, Saturn and Earth, there are some remarkably similar phenomena occurring on all three planets. One of them is the observed equatorial stratospheric oscillation of temperatures and zonal winds that on Earth is called the Quasi-Biennial Oscillation (QBO) and that was discovered on wind fields measured from meteorological balloons [1]. On Jupiter, a similar phenomenon was discovered in 1991 from observations of a cycle of temperatures in the equatorial stratosphere that oscillated with a periodicity of about 4 years [2], and was named the quasi-quadrennial oscillation due to its similarity with the QBO. This phenomenon was later shown to be irregularly affected by the development of large convective outbreaks outside the Equator changing the cadence of the regular cycle modifying the period of the oscillation from 3.9 to 5.7 years [3]. In Saturn, a similar oscillation of equatorial temperatures in the stratosphere was discovered with Cassini data [4-5] and was accompanied by the presence of an elevated narrow equatorial jet traced in the motions of the upper atmosphere equatorial hazes [6]. The Jupiter Equatorial Stratospheric Oscillation (JESO) represents an oscillation of temperatures that propagates downwards in time at pressure levels of 0.1-40 mbar and is confined between the North and South Equatorial Belts. Temperatures oscillate from having a local maximum at the equator and local minimum at ±14º latitude to the opposite situation. Because of the thermal wind relation, changes in temperatures should produce changes in stratospheric winds, but because this occurs at the equator, where Coriolis forces becomes negligible, the exact relation between meridional gradients of temperature and vertical wind shears requires the use of modified thermal wind equations that are untested with observational data [7]. Numerical modelling of JESO, and comparisons with Earth meteorology, suggest that gravity waves produced from convection at the troposphere are likely the major contributors to generating the JESO [8-9]. Recently, an intense narrow equatorial jet at stratospheric levels (200-50 mbar), close but below those most affected by the JESO changes in temperatures, has been discovered in analysis of James Webb Space Telescope (JWST) images [10] (data from July 2022 obtained as part of the Early Release Science Program 1373). This jet mimics the behaviour of the elevated equatorial narrow jet in Saturn [5]. The new jovian jet is confined at ±3º of the equator and it could represent a deep counterpart of the JESO phenomena, thus being a key part of the relation between the troposphere and stratosphere.

The lack of further JWST observations equivalent to those obtained in 2022, and the suspected temporal variability of the stratospheric jet, directed our interest to observations obtained by the Hubble Space Telescope (HST) at the strong methane absorption band at 890 nm (filter FQ889N), which are sensitive to the upper aerosols in the atmosphere at levels of around 200-300 mbar below those observed with the JWST. HST images of Jupiter at this wavelength have remained mostly unused in the past to measure winds in the planet due to the lower contrast of the images and the lower image quality than in filters sensitive to deeper levels in the troposphere. We have analysed HST images in this wavelength between 2015 and 2022 to retrieve zonal winds in the equatorial region and study potential variabilities in zonal jets, optical opacities and hazes altitudes. In this talk, we will present a thorough survey of zonal winds and clouds opacity and altitude results together with a close comparison with JWST data in 2022 and the published studies of the thermal aspects of the JESO [e.g. 11] that will enable us to understand in more detail the troposphere-stratosphere connection.

 

References: [1] Baldwin, Gray, Dunkerton et al., The quasi-biennial oscillation. Reviews of Geophysics , 39, 179-229 (2001). [2] Leovy, C., Friedson, A. & Orton, G. The quasiquadrennial oscillation of Jupiter's equatorial stratosphere. Nature 354, 380–382 (1991). [3] Antuñano, A., Cosentino, R.G., Fletcher, L.N. et al. Fluctuations in Jupiter’s equatorial stratospheric oscillation. Nat Astron 5, 71–77 (2021). [4] Fouchet, Guerlet, Strobel et al. An equatorial oscillation in Saturn’s middle atmosphere, Nature, 452, 200-202 (2008). [5] García-Melendo et al. A strong high altitude narrow jet at Saturn’s equator. Geophys. Res. Lett., 37, L22204 (2010). [6] Guerlet et al., Evolution of the equatorial oscillation in Saturn’s stratosphere between 2005 and 2010 from Cassini/CIRS data analysis. Geophys. Res. Lett. 38, L09201 (2011). [7] Marcus, Tollefson, Wong and de Pater. An equatorial thermal wind equation: Applications to Jupiter, Icarus, 324, 198-223 (2019). [8] Cosentino, R. G., Morales-Juberías, R., Greathouse, T., Orton, G., Johnson, P., Fletcher, L. N., & Simon, A. (2017). New observations and modeling of Jupiter's quasi-quadrennial oscillation. Journal of Geophysical Research: Planets, 122, 2719–2744. [9] Cosentino, Greathouse, Simon et al. The Effects of Waves on the Meridional Thermal Structure of Jupiter’s Stratosphere. The Planetary Science Journal. 1. 63 (2020). [10] Hueso, R., Sánchez-Lavega, A., Fouchet, T. et al. An intense narrow equatorial jet in Jupiter’s lower stratosphere observed by JWST. Nat Astron 7, 1454–1462 (2023). [11] Giles, Greathouse, Cosentino et al. Vertically-resolved observations of Jupiter’s quasi-quadriennial oscillation from 2012 to 2019. Icarus, 350,113905 (2020).

How to cite: Sanchez Arregui, M., Antuñano, A., Hueso, R., and Sanchez-Lavega, A.: A long-term study of the Jovian equatorial atmosphere at the upper troposphere-lower stratosphere from HST observations in the 890-nm methane absorption band, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-496, https://doi.org/10.5194/epsc2024-496, 2024.

EPSC2024-669 | ECP | Posters | OPS2 | OPC: evaluations required

Latitudinal Variation in Internal Heat Flux in Jupiter's Atmosphere: Effect on Weather Layer Dynamics 

Xinmiao Hu and Peter Read
Tue, 10 Sep, 10:30–12:00 (CEST) | Poster area Level 2 – Galerie | P53

Conventional weather layer General Circulation Models (GCMs) typically simulate over a height range extending only a short distance beneath the water cloud base, constrained by computational resources. Due to the limited knowledge about the environment at depth, the conditions specified at the bottom boundary of the domain are usually greatly simplified. Consequently, the influence of deeper atmospheric dynamics on cloud-level phenomena remains poorly understood. Recent observations from the Juno mission have provided new insights into the complex conditions prevailing within Jupiter's deep atmosphere. Given these advances, it is timely to re-evaluate the simple assumptions regarding the deep atmosphere currently employed in weather layer GCMs.
In this study, we challenge the conventional approach by introducing latitudinal variations in internal heat flux into a GCM of Jupiter’s atmosphere. Our model incorporates a heat flux profile that decreases from the equator to the poles, with additional complexities such as belt-and-zone contrast and hemispheric asymmetry. Preliminary results show significant deviations in weather layer atmospheric dynamics when compared to constant flux models, particularly in the equatorial regions. We discuss the underlying mechanisms driving these differences, providing insights into the coupling between Jupiter's visible weather layer and its obscured deeper layers. This work represents a step towards developing a more comprehensive GCM for Jupiter, which could also enhance our understanding of other giant planets, by incorporating more realistic conditions at the bottom boundary.

How to cite: Hu, X. and Read, P.: Latitudinal Variation in Internal Heat Flux in Jupiter's Atmosphere: Effect on Weather Layer Dynamics, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-669, https://doi.org/10.5194/epsc2024-669, 2024.

OPS4 | The Mysterious Saturn System

EPSC2024-347 | ECP | Posters | OPS4 | OPC: evaluations required

Pick-up ion distributions in the inner and middle Saturnian Magnetosphere 

Cristian Radulescu, Andrew Coates, Sven Simon, and Geraint Jones
Wed, 11 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P69

Based on the entire dataset collected by the Cassini Plasma Spectrometer, we provide a comprehensive picture of the pitch angle and velocity distributions of pick-up ions (PUIs) in Saturn’s inner and middle magnetosphere. We investigate the dependence of these distributions on Saturnian Local Time and magnetic latitude. We also search for correlations to the signatures of ion cyclotron waves observed by the Cassini magnetometer. Our survey reveals that ion pitch angle distributions have a pancake shape and their full width increases monotonically with magnetic latitude. This increase in angular width is anti-correlated with the observed amplitudes of ion cyclotron waves that are generated during the thermalization of the PUI distribution. We find no evidence of the observed, non-monotonic change of wave amplitudes with magnetic latitude mapping into the width of the pitch angle distributions. This suggests that only a small fraction of the energy deposited into the waves is transferred back to the ions to broaden the distribution. A possible reason for this is wave damping by the Maxwellian core of the distribution, formed by ions that have already been incorporated into the sub-corotating flow. In addition, wave propagation away from the magnetospheric field direction could reduce the efficiency of the energy transfer. When moving away from Saturn’s magnetic equatorial plane, the observed half-width of the velocity distributions does not evolve appreciably with latitude and L shell value. This behavior changes only outside the orbit of Rhea where the observed velocity distributions begin to broaden due to elevated plasma temperatures.

How to cite: Radulescu, C., Coates, A., Simon, S., and Jones, G.: Pick-up ion distributions in the inner and middle Saturnian Magnetosphere, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-347, https://doi.org/10.5194/epsc2024-347, 2024.

OPS5 | Exploration of Titan

EPSC2024-226 | Posters | OPS5 | OPC: evaluations required

On the evolution of the primordial hydrosphere of Titan 

Alizée Amsler Moulanier, Olivier Mousis, Alexis Bouquet, Ngan H. D. Trinh, and Kathleen E. Mandt
Tue, 10 Sep, 14:30–16:00 (CEST) | Poster area Level 2 – Galerie | P69

Titan is the only moon of the solar system harboring a thick N2-rich atmosphere, as well as a subsurface global ocean covered by an ice crust. To explore the origins and evolution of Titan’s present hydrosphere volatile inventory, it is necessary to retrace the influence of the moons’ formation scenarios on this inventory. To do so, we need to explore the post-accretion processes that could impact the distribution of volatiles in the hydrosphere. Especially, we investigate the evolution of the early “open-ocean” phase of Titan, which took place shortly after accretion before the ice crust formation.

Our work focuses on modelling the ocean-atmosphere equilibrium which took place over this period, coupled to a statistical clathrate formation model (fig1). Hence, we compute the vapor-liquid equilibrium between the early atmosphere and ocean, as well as the chemical equilibria happening within the ocean, following the approach described by Marounina et al. (2018) [1], to investigate the primitive hydrosphere of Titan. This approach is coupled to a statistical thermodynamic model introduced in Mousis et al (2013) [2], to assess how clathrates formation impact the early atmosphere and the liquid phase’s composition throughout time. By exploring different building blocks composition and how they affect the primordial hydrosphere composition, we highlight the consequence of the ratio of dissolved CO2 and NH3 on the distribution of partial pressures in the primordial atmosphere of Titan as well as the composition and buoyancy of clathrate hydrates formed in the ocean.

Fig1: Model loop scheme. At each time step, the thermodynamical equilibrium between the primordial atmosphere and ocean is computed. If the conditions to form clathrate hydrates are met, the composition of the liquid phase is updated, and the equilibrium recalculated.

Our coupled model allows for an assessment of the impact of the initial distribution of volatiles on the thermodynamic equilibrium between the early moon’s atmosphere and ocean, with or without clathrate formation throughout time. Based on a range of volatile distributions brought by the accreted material of cometary origins, different formation scenarios of Titan’s primordial hydrosphere are then investigated. Such a model highlights the influence of aqueous species’ chemistry on volatile’s dissolution in the ocean, especially CO2 and NH3, and how clathrate hydrate formation can further impact this distribution.This model will be further improved considering geochemical exchanges taking place between the moon’s rocky core and ocean.

 

References:

[1] N. Marounina, O. Grasset, G. Tobie, and S. Carpy, “Role of the global water ocean on the evolution of Titan’s primitive atmosphere,” Icarus, vol. 310, pp. 127–139, Aug. 2018.

[2] O. Mousis, A. Lakhlifi, S. Picaud, M. Pasek, and E. Chassefière, “On the Abundances of Noble and Biologically Relevant Gases in Lake Vostok, Antarctica,” Astrobiology, vol. 13, pp. 380–390, Apr. 2013

How to cite: Amsler Moulanier, A., Mousis, O., Bouquet, A., Trinh, N. H. D., and Mandt, K. E.: On the evolution of the primordial hydrosphere of Titan, Europlanet Science Congress 2024, Berlin, Germany, 8–13 Sep 2024, EPSC2024-226, https://doi.org/10.5194/epsc2024-226, 2024.

EPSC2024-428 | ECP | Posters | OPS5 | OPC: evaluations required

Seasonal variation of Titan’s stratospheric tilt 

Lucy Wright, Nicholas A. Teanby, Patrick G. J. Irwin, Conor A. Nixon, and Joshua S. Ford
Tue, 10 Sep, 14:30–16:00 (CEST) |