Stacking Lake Ice Analog Cryotectonic Dynamics on Dione’s Wispy Terrain

: Wispy Terrain, with its chasmata, is one of the enigmatic regions of Dione. It consists of quasi-parallel graben, and troughs, in parts with horsts, indicating extensional and shear stresses. This study introduces some observations of compression-related features and proposes a new regional formation model. The study of the relationship between impact craters and tectonic features revealed certain “lost” parts of some crosscut craters, indicating additional cryotectonic features, the appearance of accretionary prism-like phenomena, and, theoretically, subsumption-like processes. This study provides new information about the surface renewal processes at one of the youngest and probably still active regions of Dione.


Introduction
Saturn's icy moon Dione and its enigmatic feature, the so-called Wispy Terrain, have been the object of scientific debate for a long time. It was named Wispy Terrain following the analysis of Voyager images, referring to the frequent appearance of wispy streaks, markings, or lines. The early astrogeological mapping of Dione provides a detailed first view (as much it was possible based on the Voyager images) of the geological features of the Wispy Terrain as well, but, beyond a chronostratigraphic proposal, no clear age determination could be made about its formation time [1,2].
Compared to the densely cratered terrains, possibly the oldest surface on the satellite (~3-4.2 Ga), results from the Eurotas Chasmata ( Figure 1) indicate that the cryotectonic processes in the region might be active until~3Ga, or even until~1 Ga ago [3,4]. Almost three decades after the first geological mapping [1,2], the main terrains were re-evaluated in light of the new data [5]. The markings of the Wispy Terrain are defined as a set of quasi-parallel troughs, scarps, and "horst and graben" systems, indicating extensional and shear stress in the region. Please note that since the early studies, mentioned compressional stress-related landforms [2], most of the studies emphasize extension and shear stress in the ice-crust of the satellite. The chasmata systems (e.g., the Eurotas and Palatine Chasmata) are defined as Fractured Cratered Plains and subdivided into three facies types, based on the timing of their formation, which date back to 3.7 Ga (with 100 Ma uncertainties) or alternatively happened between 2.7 Ga and 260 Ma [5].
Following the re-definition of the terrains on Dione, the newly computed ages (4.3 [+0.2/−2.7] Ga-> 45 km crater age, and 2.5 [+2.0/−1.9] Ga-≤ 45 km crater age) for the Faulted Terrain (also known as the Wispy Terrain) showed that these ages do not reflect the timing of cryotectonic activity (i.e., faulting), but the time span over which the larger craters and their ejecta blankets erased the smaller craters in the region [6]. mapping of impact craters, as a part of another ongoing research, (ii) the observation of fragmentary craters, and (iii) a viral video of stacking lake ice, seen during a short break. These three "circumstances" lead to a theory and discovery about the potential of subsumption (a subduction-like) mechanism [14] on Dione. Here, we provide some evidence that supports the existence of compressional stress field and subsumption during the evolution of the Wispy Terrain and its chasmata. Based on the analog geological characteristics between the fault system of chasmata and Earth's regions with the extension-drive divergent tectonic plates, Dione's Faulted Terrain was defined as hemisphere-scale rift zones [7,8].
The possibility of a subsurface ocean under the ice shell [9], similar to Dione's active neighbor, Enceladus [10], brought the promise of a still active surface and a potential for life harbored under the icy shell. The study of the stratigraphic relationship between the craters and faults on the Wispy Terrain suggests that the faulting is a geologically very recent event, dating back to 0.3-0.79 Ga [11]. Some of the newest studies go even further and suggest that the upper limit for the age of studied fault on Wispy Terrain is only 152 Ma, which supports the hypothesis that the cryotectonism might still be active or was active a very short time ago (c.a. 100 Ma) on the satellite [12].
In certain situations, mere coincidence would lead to discovery (or false belief, at the same time) [13]. In the case of this study, such a coincidence has three components: (i) the mapping of impact craters, as a part of another ongoing research, (ii) the observation of fragmentary craters, and (iii) a viral video of stacking lake ice, seen during a short break. These three "circumstances" lead to a theory and discovery about the potential Phys. Sci. Forum 2023, 7, 18 3 of 7 of subsumption (a subduction-like) mechanism [14] on Dione. Here, we provide some evidence that supports the existence of compressional stress field and subsumption during the evolution of the Wispy Terrain and its chasmata.

Data and Methods
The studied region, the Wispy Terrains, is spreading approximately between latitude 40 • and −60 • and longitude 260 • to 350 • (westward) (Figure 1). The focus of this research is between Palatine and the "southern" section of Eurotas Chasmata. The map is based on Cassini-Voyager Global Mosaic 154m v1 map, which can be found at the internet site Astropedia-Lunar and Planetary Cartographic Catalog [15][16][17][18]. The applied nomenclature follows the recommendation of the Gazetteer of Planetary Nomenclature [19].
The GIS research and the mapping of craters and other geological features were performed by QGIS 3.22 software. In previous studies, the mapping of craters in Dione was limited to craters with ≥ 4 km crater diameter due to problems during the identification of smaller craters in image mosaics with lower resolution [6]. The crater mapping also had some limitations in certain areas due to the distortion of the map.

Crater Preservation in Various Cryotectonic Settings-Theory behind the Research
In some cases of short-distance horizontal allocation of the crater fragments (strike-slip faults), the cryotectonic activity "only" causes the distortion of the crater rim, resulting in visible sharp edges in the supposedly curved-shaped rim. The whole crater is still recognizable (Figure 2a). In certain cases of horizontal allocation (strike-slip faults), the parts of the crater may move further distance, but even in such cases, the component parts are recognizable along the fault trace, and the crater is re-constructible ( Figure 2b).

Data and Methods
The studied region, the Wispy Terrains, is spreading approximately between latitude 40° and −60° and longitude 260° to 350° (westward) (Figure 1). The focus of this research is between Palatine and the "southern" section of Eurotas Chasmata. The map is based on Cassini-Voyager Global Mosaic 154m v1 map, which can be found at the internet site Astropedia-Lunar and Planetary Cartographic Catalog [15][16][17][18]. The applied nomenclature follows the recommendation of the Gazetteer of Planetary Nomenclature [19].
The GIS research and the mapping of craters and other geological features were performed by QGIS 3.22 software. In previous studies, the mapping of craters in Dione was limited to craters with ⌀ ≥ 4 km crater diameter due to problems during the identification of smaller craters in image mosaics with lower resolution [6]. The crater mapping also had some limitations in certain areas due to the distortion of the map.

Crater Preservation in Various Cryotectonic Settings-Theory behind the Research
In some cases of short-distance horizontal allocation of the crater fragments (strikeslip faults), the cryotectonic activity "only" causes the distortion of the crater rim, resulting in visible sharp edges in the supposedly curved-shaped rim. The whole crater is still recognizable (Figure 2a). In certain cases of horizontal allocation (strike-slip faults), the parts of the crater may move further distance, but even in such cases, the component parts are recognizable along the fault trace, and the crater is re-constructible (Figure 2b).  Normal (and in some case reverse) faulting may cause the allocation of the crater fragments along a fault plane with a certain dip to the horizontal plane. Based on the angle of the dip (deep to shallow) and the orientation of the fault planes of individual blocks in a fault system, various fault types may form on the crust of icy satellites, such as normal faults, horst, and graben structures, domino-style fault blocks, and listric normal faults [20]. Among many similarities, there is one that is important regarding the preservation of the crater fragments, separated by deep or shallow angle normal faulting: despite the allocation of the crater parts, they are recognizable on the original surface of the footwall and hanging wall block as well, and the "whole" crater is re-constructible (Figure 2c).
The commonality in the suggested processes above is that in every case, the parts of the crater are recognizable, and the crater itself is re-constructible. One of the exceptional cases would be the formation of thrust faults, a reverse fault with a shallow dip, with a hanging wall that moves up and over the footwall. In such cases, the hanging wall may "cover" a certain part of the footwall top, resulting in fragmentary or incomplete craters ("half-craters") without any recognizable and completing other connecting crater parts in their surroundings (Figure 2d).
The key (and simplified) difference between the various faults is that various types of normal faults are formed in extensional-and thrust faults (low-angle reverse faults) in a compressional stress field. Strike-slip faults are caused by horizontal compression, but compared to thrust faults, their displacement is quasi-parallel to the compressional force (shear stress). In this study, the focus is on the fragmentary craters, which may indicate thrust fault cryotectonic settings, formed in a compressional stress field (Figure 2d).

Results and Discussion-Putative Thrust Faults of the Wispy Terrain
The preliminary geological mapping of the area resulted in the identification of close to 450 main tectonic lines and 2500 craters, including the ones with missing, "untraceable" parts ( Figure 1c). Despite the preliminary nature of the map, the relationship between the tectonically active zones and the increasing number of fragmentary craters along the tectonic lines seems clear, which may indicate thrust faulting (Figure 2d) as one of the key cryotectonic processes in the region. Compared to the previous studies [7,8,20], thrust faults indicate a compressional stress field and convergence of plates. In a terrestrial environment, thrust fault may develop when two blocks collide, and it may lead to the formation of so-called splay and décollement faults, one of the most characteristic components of subduction zones (e.g., found in the accretionary wedges) [21].
To verify the results of the crater mapping, the walls of larger craters found in a region and crosscut by tectonic lines were investigated. Despite the potential bias due to the resolution of the images, certain structural geological characteristics might be recognizable (Figure 3). Without any knowledge about the appearance of fragmentary craters in the surroundings, putative faults, observed in Figure 3a,b, might be identified as domino-style normal faults [20,22] or thrust faults, but supported by various pieces of evidence, the possibility of the former ones was already excluded [20]. As an alternative to thrust faults, the features in Figure 3c may be recognized as listric normal faults. Unfortunately, in this resolution, it is hard to recognize key features such as the thrust sheets (the "imbricating splays" of the splay and décollement faults) or the tilted face of the blocks (a listric normal fault), which may help to decide between the two features. As a comparison of the putative thrust faults, the features in Figure 3d most likely indicate a normal fault system of a "graben" structure between two uplifted blocks. Despite the uncertainty about the faults in Figure 3a-c, they may be identified as thrust/splay and décollement faults. Phys. Sci. Forum 2023, 7, x 5 of 7 -orientation of (quasi-)vertical block movement, 9-orientation of (quasi-)horizontal block movement; 10-exposed crater wall.

Conclusions-Subsumption or Chaotic Ice Stacking… or Both?
Extensional stress fields seem common on Dione`s crust, resulting in the development of various types of normal fault systems [20]. The scale of the dilatational features spans from step faults, horst and graben structures to hemisphere-scale rift zones [7,8]. The discovery of potential thrust/splay and décollement faults may make the surface renewal models more complex than thought before. The appearance of compression-related tectonic features, especially the ones whose morphology is reminiscent of the structure of an accretionary prism in a terrestrial environment, may suggest the existence of subsumption (subduction-like) processes on Dione (Figure 4). Ideally, it means that Wilson cyclelike tectonic cycles may appear in icy planetary bodies [23], and Dione`s Wispy Terrain, which consists of divergent and convergent sections as well, may turn out to be the model region for it. Rifting may be triggered by endogenic (e.g., phase change within the satellite, solid-state convection in the crust, and thermal plumes) and exogenic (e.g., diurnal tides, (quasi-)vertical block movement, 9-orientation of (quasi-)horizontal block movement; 10-exposed crater wall.

Conclusions-Subsumption or Chaotic Ice Stacking . . . or Both?
Extensional stress fields seem common on Dione's crust, resulting in the development of various types of normal fault systems [20]. The scale of the dilatational features spans from step faults, horst and graben structures to hemisphere-scale rift zones [7,8]. The discovery of potential thrust/splay and décollement faults may make the surface renewal models more complex than thought before. The appearance of compression-related tectonic features, especially the ones whose morphology is reminiscent of the structure of an accretionary prism in a terrestrial environment, may suggest the existence of subsumption (subduction-like) processes on Dione (Figure 4). Ideally, it means that Wilson cycle-like tectonic cycles may appear in icy planetary bodies [23], and Dione's Wispy Terrain, which consists of divergent and convergent sections as well, may turn out to be the model region for it. Rifting may be triggered by endogenic (e.g., phase change within the satellite, solidstate convection in the crust, and thermal plumes) and exogenic (e.g., diurnal tides, tidal forces, orbital forcing, and non-synchronous rotation of the ice shell) processes [7,8,24,25]. This section of the ice crust may spread due to the continuous material accretion via cryovolcanic activity [12] and at a certain geological moment it may collide with thicker, more stable, and possibly older terrain of the crust (Figure 1c); e.g., some evolving part of the Faulted Terrain collides the much older Intermediate Cratered Terrain [6].
Phys. Sci. Forum 2023, 7, x 6 of 7 tidal forces, orbital forcing, and non-synchronous rotation of the ice shell) processes [7,8,24,25]. This section of the ice crust may spread due to the continuous material accretion via cryovolcanic activity [12] and at a certain geological moment it may collide with thicker, more stable, and possibly older terrain of the crust (Figure 1c); e.g., some evolving part of the Faulted Terrain collides the much older Intermediate Cratered Terrain [6]. Only one question remains, would such collision develop Earth analog subduction connected features (e.g., the identified thrust faults; Figure 4), or would result in a more chaotic cryotetconic setting, similar to the stacking lake ice, observed at Lake Superior [26]? The next goal of our research is to map various crater-fault relations and reconstruct the stress field that appears in the wispy Terrain, especially in the region of the Palatine and Eurotas Chasmata, and may be able to answer such questions.