Consequences of large impacts on Enceladus' core shape

International audienceThe intense activity on Enceladus suggests a differentiated interior consisting of a rocky core, an internal ocean and an icy mantle. However, topography and gravity data suggests large heterogeneity in the interior, possibly including significant core topography. In the present study, we investigated the consequences of collisions with large impactors on the core shape. We performed impact simulations using the code iSALE2D considering large differentiated impactors with radius ranging between 25 and 100 km and impact velocities ranging between 0.24 and 2.4 km/s. Our simulations showed that the main controlling parameters for the post-impact shape of Enceladus’ rock core are the impactor radius and velocity and to a lesser extent the presence of an internal water ocean and the porosity and strength of the rock core. For low energy impacts, the impactors do not pass completely through the icy mantle. Subsequent sinking and spreading of the impactor rock core lead to a positive core topographic anomaly. For moderately energetic impacts, the impactors completely penetrate through the icy mantle, inducing a negative core topography surrounded by a positive anomaly of smaller amplitude. The depth and lateral extent of the excavated area is mostly determined by the impactor radius and velocity. For highly energetic impacts, the rocky core is strongly deformed, and the full body is likely to be disrupted. Explaining the long-wavelength irregular shape of Enceladus’ core by impacts would imply multiple low velocity (<2.4 km/s) collisions with deca-kilometric differentiated impactors, which is possible only after the LHB period

Enhanced visualisation of complex thermofluid data: Vertical and horizontal combined convection and microscale heat transfer cases

Copyright @ 2000 UITIn general, convective heat transfer is an 'n-dimensional' problem where n is well in excess of 3 for steady flows. Traditionally, the method of dimensional analysis results in a small number of dimensionless groups. In the case of steady forced convection these can reduce to three, namely the Nusselt (Nu), Reynolds (Re) and Prandtl (Pr) numbers, for heat transfer, fluid flow regime and fluid properties respectively. Again, traditionally, data are presented on log-log graphs, say of Nu versus Re, with Pr being a possible third parameter. For natural convection, the Grashof number (Gr) expresses buoyancy effects in place of Re, while for combined (natural and forced) convection Gr becomes additional to, rather than replacing, Re. Using sets of data for: (a) vertical combined convection in nuclear safety, (b) horizontal combined convection review material, and (c) microchannel heat transfer, in the first part of this paper we survey this problem. We reach the following conclusions: that heat transfer data are presented in either 'holistic' or 'reductive' modes, and that other thermodynamic performance data are related to the generic scientific cases of (a) 3-dimensional space and (b) multi-dimensional space. In the second part of the paper we present a first attempt at applying design-type procedures to specifying this problem. Visualisation priorities are suggested from which particular solutions will be developed in future

Complex crater formation by oblique impacts on the Earth and Moon

Almost all meteorite impacts occur at oblique incidence angles, but the effect of impact angle on crater size is not well understood, especially for large craters. To improve oblique impact crater scaling, we present a suite of simulations of complex crater formation on Earth and the Moon over a range of impact angles, velocities and impactor sizes. We show that crater diameter is larger than predicted by existing scaling relationships for oblique impacts; there is little dependence on obliquity for impacts steeper than 45° from the horizontal. Crater depth, volume and diameter depend on impact angle in different ways—relatively shallower craters are formed by more oblique impacts. Our simulation results have implications for how crater populations are determined from impactor populations and vice-versa. They suggest that existing approaches to account for impact obliquity may underestimate the number of complex craters larger than a given size by as much as one-third

Spud 1.0: generalising and automating the user interfaces of scientific computer models

The interfaces by which users specify the scenarios to be simulated by scientific computer models are frequently primitive, under-documented and ad-hoc text files which make using the model in question difficult and error-prone and significantly increase the development cost of the model. In this paper, we present a model-independent system, Spud, which formalises the specification of model input formats in terms of formal grammars. This is combined with an automated graphical user interface which guides users to create valid model inputs based on the grammar provided, and a generic options reading module, libspud, which minimises the development cost of adding model options. &lt;br&gt;&lt;br&gt; Together, this provides a user friendly, well documented, self validating user interface which is applicable to a wide range of scientific models and which minimises the developer input required to maintain and extend the model interface

Combining shock barometry with numerical modeling: insights into complex crater formation – The example of the Siljan impact structure (Sweden)

Siljan, central Sweden, is the largest known impact structure in Europe. It was formed at about 380 Ma, in the late Devonian period. The structure has been heavily eroded to a level originally located underneath the crater floor, and to date, important questions about the original size and morphology of Siljan remain unanswered. Here we present the results of a shock barometry study of quartz-bearing surface and drill core samples combined with numerical modeling using iSALE. The investigated 13 bedrock granitoid samples show that the recorded shock pressure decreases with increasing depth from 15 to 20 GPa near the (present) surface, to 10–15 GPa at 600 m depth. A best-fit model that is consistent with observational constraints relating to the present size of the structure, the location of the downfaulted sediments, and the observed surface and vertical shock barometry profiles is presented. The best-fit model results in a final crater (rim-to-rim) diameter of ~65 km. According to our simulations, the original Siljan impact structure would have been a peak-ring crater. Siljan was formed in a mixed target of Paleozoic sedimentary rocks overlaying crystalline basement. Our modeling suggests that, at the time of impact, the sedimentary sequence was approximately 3 km thick. Since then, there has been around 4 km of erosion of the structure

The effects of asteroid layering on ejecta mass-velocity distribution and implications for impact momentum transfer

Most bodies in the Solar System do not have a homogeneous structure. Understanding the outcome of an impact into regolith layers of different properties is especially important for NASA’s Double Asteroid Redirection Test (DART) and ESA’s Hera missions. Here we used the iSALE shock physics code to simulate the DART impact into three different target scenarios in the strength regime: a homogeneous porous half-space; layered targets with a porous weak layer overlying a stronger bedrock; and targets with exponentially decreasing porosity with depth. For each scenario we determined the sensitivity of crater morphology, ejecta mass-velocity distribution and momentum transferred from the impact for deflection, , to target properties and structure. We found that for a homogeneous porous half-space, cohesion and porosity play a significant role and the DART impact is expected to produce a between 1 and 3. In a two-layer target scenario, the presence of a less porous, stronger lower layer close to the surface can cause both amplification and reduction of ejected mass and momentum relative to the homogeneous upper-layer case. For the case of DART, the momentum enhancement can change by up to 90%. Impacts into targets with an exponentially decreasing porosity with depth only produced an enhancement in the ejected mass and momentum for sharp decreases in porosity that occur within 6 m of the asteroid surface. Together with measurements of the DART crater by the Hera mission, these results can be used to test the predictive capabilities of numerical models of asteroid deflection

Clastic polygonal networks around Lyot crater, Mars: Possible formation mechanisms from morphometric analysis

Polygonal networks of patterned ground are a common feature in cold-climate environments. They can form through the thermal contraction of ice-cemented sediment (i.e. formed from fractures), or the freezing and thawing of ground ice (i.e. formed by patterns of clasts, or ground deformation). The characteristics of these landforms provide information about environmental conditions. Analogous polygonal forms have been observed on Mars leading to inferences about environmental conditions. We have identified clastic polygonal features located around Lyot crater, Mars (50°N, 30°E). These polygons are unusually large ( > 100 m diameter) compared to terrestrial clastic polygons, and contain very large clasts, some of which are up to 15 metres in diameter. The polygons are distributed in a wide arc around the eastern side of Lyot crater, at a consistent distance from the crater rim. Using high-resolution imaging data, we digitised these features to extract morphological information. These data are compared to existing terrestrial and Martian polygon data to look for similarities and differences and to inform hypotheses concerning possible formation mechanisms. Our results show the clastic polygons do not have any morphometric features that indicate they are similar to terrestrial sorted, clastic polygons formed by freeze-thaw processes. They are too large, do not show the expected variation in form with slope, and have clasts that do not scale in size with polygon diameter. However, the clastic networks are similar in network morphology to thermal contraction cracks, and there is a potential direct Martian analogue in a sub-type of thermal contraction polygons located in Utopia Planitia. Based upon our observations, we reject the hypothesis that polygons located around Lyot formed as freeze-thaw polygons and instead an alternative mechanism is put forward: they result from the infilling of earlier thermal contraction cracks by wind-blown material, which then became compressed and/or cemented resulting in a resistant fill. Erosion then leads to preservation of these polygons in positive relief, while later weathering results in the fracturing of the fill material to form angular clasts. These results suggest that there was an extensive area of ice-rich terrain, the extent of which is linked to ejecta from Lyot crater