-
Modeling the formation of Selk impact crater on Titan: Implications for Dragonfly
Authors:
Shigeru Wakita,
Brandon C. Johnson,
Jason M. Soderblom,
Jahnavi Shah,
Catherine D. Neish,
Jordan K. Steckloff
Abstract:
Selk crater is an $\sim$ 80 km diameter impact crater on the Saturnian icy satellite, Titan. Melt pools associated with impact craters like Selk provide environments where liquid water and organics can mix and produce biomolecules like amino acids. It is partly for this reason that the Selk region has been selected as the area that NASA's Dragonfly mission will explore and address one of its prima…
▽ More
Selk crater is an $\sim$ 80 km diameter impact crater on the Saturnian icy satellite, Titan. Melt pools associated with impact craters like Selk provide environments where liquid water and organics can mix and produce biomolecules like amino acids. It is partly for this reason that the Selk region has been selected as the area that NASA's Dragonfly mission will explore and address one of its primary goals: to search for biological signatures on Titan. Here we simulate Selk-sized impact craters on Titan to better understand the formation of Selk and its melt pool. We consider several structures for the icy target material by changing the thickness of the methane clathrate layer, which has a substantial effect on the target thermal structure and crater formation. Our numerical results show that a 4 km-diameter-impactor produces a Selk-sized crater when 5-15 km thick methane clathrate layers are considered. We confirm the production of melt pools in these cases and find that the melt volumes are similar regardless of methane clathrate layer thickness. The distribution of the melted material, however, is sensitive to the thickness of the methane clathrate layer. The melt pool appears as a torus-like shape with a few km depth in the case of 10-15 km thick methane clathrate layer, and as a shallower layer in the case of a 5 km thick clathrate layer. Melt pools of this thickness may take tens of thousands of years to freeze, allowing more time for complex organics to form.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
The Effects of Early Collisional Evolution on Amorphous Water Ice Bodies
Authors:
Jordan K. Steckloff,
Gal Sarid,
Brandon C. Johnson
Abstract:
Conditions in the outer protoplanetary disk during Solar System formation were thought to be favorable for the formation of amorphous water ice (AWI),a glassy phase of water ice. However, subsequent collisional processing could have shock crystallized any AWI present. Here we use the iSALE shock physics hydrocode to simulate impacts between large icy bodies at impact velocities relevant to these c…
▽ More
Conditions in the outer protoplanetary disk during Solar System formation were thought to be favorable for the formation of amorphous water ice (AWI),a glassy phase of water ice. However, subsequent collisional processing could have shock crystallized any AWI present. Here we use the iSALE shock physics hydrocode to simulate impacts between large icy bodies at impact velocities relevant to these collisional environments, and then feed these results into a custom-built AWI crystallization script, to compute how much AWI crystallizes/survives these impact events. We find that impact speeds between icy bodies post-planet migration (i.e., between trans-Neptunian Objects or TNOs) are too slow to crystallize any meaningful fraction of AWI. During planet migration, however, the amount of AWI that crystallizes is highly stochastic: relatively little AWI crystallizes at lower impact velocities (less than ~2 km/s), yet most AWI present in the bodies (if equal sized) or impactor and impact site (if different sizes) crystallizes at higher impact velocities (greater than ~4 km/s). Given that suspected impact speeds during planet migration were ~2-4 km/s, this suggests that primordial AWI's ability to survive planet migration is highly stochastic. However, if proto-EKB objects and their fragments experienced multiple impact events, nearly all primordial AWI could have crystallized; such a highly collisional proto-EKB during planet migration is consistent with the lack of any unambiguous direct detection of AWI on any icy body. Ultimately, primordial AWI's survival to the present day depends sensitively on the proto-EKB's size-frequency distribution, which is currently poorly understood.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Effect of impact velocity and angle on deformational heating and post-impact temperature
Authors:
Shigeru Wakita,
Hidenori Genda,
Kosuke Kurosawa,
Thomas M. Davison,
Brandon C. Johnson
Abstract:
The record of impact induced shock-heating in meteorites is an important key for understanding the collisional history of the solar system. Material strength is important for impact heating, but the effect of impact angle and impact velocity on shear heating remains poorly understood. Here, we report three-dimensional oblique impact simulations, which confirm the enhanced heating due to material s…
▽ More
The record of impact induced shock-heating in meteorites is an important key for understanding the collisional history of the solar system. Material strength is important for impact heating, but the effect of impact angle and impact velocity on shear heating remains poorly understood. Here, we report three-dimensional oblique impact simulations, which confirm the enhanced heating due to material strength and explore the effects of impact angle and impact velocity. We find that oblique impacts with an impact angle that is steeper than 45 degree produce a similar amount of heated mass as vertical impacts. On the other hand, grazing impacts produce less heated mass and smaller heated regions compared to impacts at steeper angles. We derive an empirical formula of the heated mass, as a function of the impact angle and velocity. This formula can be used to estimate the impact conditions (velocity and angle) that had occurred and caused Ar loss in the meteoritic parent bodies. Furthermore, our results indicate that grazing impacts at higher impact velocities could generate a similar amount of heated material as vertical impacts at lower velocities. As the heated material produced by grazing impacts has experienced lower pressure than the material heated by vertical impacts, our results imply that grazing impacts may produce weakly shock-heated meteorites.
△ Less
Submitted 16 August, 2022;
originally announced August 2022.
-
Methane-saturated layers limit the observability of impact craters on Titan
Authors:
Shigeru Wakita,
Brandon C. Johnson,
Jason M. Soderblom,
Jahnavi Shah,
Catherine D. Neish
Abstract:
As the only icy satellite with a thick atmosphere and liquids on its surface, Titan represents a unique end-member to study the impact cratering process. Unlike craters on other Saturnian satellites, Titan's craters are preferentially located in high-elevation regions near the equator. This led to the hypothesis that the presence of liquid methane in Titan's lowlands affects crater morphology, mak…
▽ More
As the only icy satellite with a thick atmosphere and liquids on its surface, Titan represents a unique end-member to study the impact cratering process. Unlike craters on other Saturnian satellites, Titan's craters are preferentially located in high-elevation regions near the equator. This led to the hypothesis that the presence of liquid methane in Titan's lowlands affects crater morphology, making them difficult to identify. This is because surfaces covered by weak fluid-saturated sediment limit the topographic expression of impact craters, as sediment moves into the crater cavity shortly after formation. Here we simulate crater-forming impacts on Titan's surface, exploring how a methane-saturated layer overlying a methane-clathrate layer affects crater formation. Our numerical results show that impacts form smaller craters in a methane-clathrate basement than a water-ice basement, due to the differences in strength. We find that the addition of a methane-saturated layer atop this basement reduces crater depths and influences crater morphology. The morphology of impact craters formed in a thin methane-saturated layer are similar to those in a "dry" target, but a thick saturated layer produces an impact structure with little to no topography. A thick methane-saturated layer (thicker than 40% of the impactor diameter) could explain the dearth of craters in the low-elevation regions on Titan.
△ Less
Submitted 24 January, 2022;
originally announced January 2022.
-
Effect of ice sheet thickness on formation of the Hiawatha impact crater
Authors:
Elizabeth A. Silber,
Brandon C. Johnson,
Evan Bjonnes,
Joseph A. MacGregor,
Nicolaj K. Larsen,
Sean E. Wiggins
Abstract:
The discovery of a large putative impact crater buried beneath Hiawatha Glacier along the margin of the northwestern Greenland Ice Sheet has reinvigorated interest into the nature of large impacts into thick ice masses. This circular structure is relatively shallow and exhibits a small central uplift, whereas a peak-ring morphology is expected. This discrepancy may be due to long-term and ongoing…
▽ More
The discovery of a large putative impact crater buried beneath Hiawatha Glacier along the margin of the northwestern Greenland Ice Sheet has reinvigorated interest into the nature of large impacts into thick ice masses. This circular structure is relatively shallow and exhibits a small central uplift, whereas a peak-ring morphology is expected. This discrepancy may be due to long-term and ongoing subglacial erosion but may also be explained by a relatively recent impact through the Greenland Ice Sheet, which is expected to alter the final crater morphology. Here we model crater formation using hydrocode simulations, varying pre-impact ice thickness and impactor composition over crystalline target rock. We find that an ice-sheet thickness of 1.5 or 2 km results in a crater morphology that is consistent with the present morphology of this structure. Further, an ice sheet that thick substantially inhibits ejection of rocky material, which might explain the absence of rocky ejecta in most existing Greenland deep ice cores if the impact occurred during the late Pleistocene. From the present morphology of the putative Hiawatha impact crater alone, we cannot distinguish between an older crater formed by a pre-Pleistocene impact into ice-free bedrock or a younger, Pleistocene impact into locally thick ice, but based on our modeling we conclude that latter scenario is possible.
△ Less
Submitted 16 April, 2021;
originally announced April 2021.
-
HD 145263: Spectral Observations of Silica Debris Disk Formation via Extreme Space Weathering?
Authors:
C. M. Lisse,
H. Y. A. Meng,
M. L. Sitko,
A. Morlok,
B. C. Johnson,
A. P. Jackson,
R. J. Vervack Jr.,
C. H. Chen,
S. J. Wolk,
M. D. Lucas,
M. Marengo,
D. T. Britt
Abstract:
We report here time domain infrared spectroscopy and optical photometry of the HD145263 silica-rich circumstellar disk system taken from 2003 through 2014. We find an F4V host star surrounded by a stable, massive 1e22 - 1e23 kg (M_Moon to M_Mars) dust disk. No disk gas was detected, and the primary star was seen rotating with a rapid ~1.75 day period. After resolving a problem with previously repo…
▽ More
We report here time domain infrared spectroscopy and optical photometry of the HD145263 silica-rich circumstellar disk system taken from 2003 through 2014. We find an F4V host star surrounded by a stable, massive 1e22 - 1e23 kg (M_Moon to M_Mars) dust disk. No disk gas was detected, and the primary star was seen rotating with a rapid ~1.75 day period. After resolving a problem with previously reported observations, we find the silica, Mg-olivine, and Fe-pyroxene mineralogy of the dust disk to be stable throughout, and very unusual compared to the ferromagnesian silicates typically found in primordial and debris disks. By comparison with mid-infrared spectral features of primitive solar system dust, we explore the possibility that HD 145263's circumstellar dust mineralogy occurred with preferential destruction of Fe-bearing olivines, metal sulfides, and water ice in an initially comet-like mineral mix and their replacement by Fe-bearing pyroxenes, amorphous pyroxene, and silica. We reject models based on vaporizing optical stellar megaflares, aqueous alteration, or giant hypervelocity impacts as unable to produce the observed mineralogy. Scenarios involving unusually high Si abundances are at odds with the normal stellar absorption near-infrared feature strengths for Mg, Fe, and Si. Models involving intense space weathering of a thin surface patina via moderate (T < 1300 K) heating and energetic ion sputtering due to a stellar superflare from the F4V primary are consistent with the observations. The space weathered patina should be reddened, contain copious amounts of nanophase Fe, and should be transient on timescales of decades unless replenished.
△ Less
Submitted 15 March, 2020;
originally announced March 2020.
-
Ferrovolcanism on metal worlds and the origin of pallasites
Authors:
Brandon C. Johnson,
Michael M. Sori,
Alexander J. Evans
Abstract:
As differentiated planetesimals cool, their cores can solidify from the outside-in, as evidenced by paleomagnetic measurements and cooling rate estimates of iron meteorites. The details of outside-in solidification and fate of residual core melt are poorly understood. For a core primarily composed of Fe and Ni alloyed with lighter constituent elements, like sulfur, such inward core growth would li…
▽ More
As differentiated planetesimals cool, their cores can solidify from the outside-in, as evidenced by paleomagnetic measurements and cooling rate estimates of iron meteorites. The details of outside-in solidification and fate of residual core melt are poorly understood. For a core primarily composed of Fe and Ni alloyed with lighter constituent elements, like sulfur, such inward core growth would likely be achieved by growth of solid FeNi dendrites. Growth of FeNi dendrites results in interconnected pockets of residual melt that become progressively enriched in sulfur up to a eutectic composition of 31 wt percent sulfur as FeNi continues to solidify. Here we show that regions of residual sulfur-enriched FeNi melt in the core attain sufficient excess pressures to propagate via dikes into the mantle. Thus, core material will intrude into the overlying rocky mantle or possibly even erupt onto the plantesimals surface. We refer to these processes collectively as ferrovolcanism. Our calculation show that ferrovolcanic surface eruptions are more likely on bodies with mantles less than 50 km thick. We show that intrusive ferromagmatism can produce pallasites, an enigmatic class of meteorites composed of olivine crystals entrained in a matrix of FeNi metal. Ferrovolcanic eruptions may explain the observations that Psyche has a bulk density inconsistent with iron metorites yet shows evidence of a metallic surface composition.
△ Less
Submitted 16 September, 2019;
originally announced September 2019.
-
Impact crater morphology and the structure of Europa's ice shell
Authors:
Elizabeth A. Silber,
Brandon C. Johnson
Abstract:
We performed numerical simulations of impact crater formation on Europa to infer the thickness and structure of its ice shell. The simulations were performed using iSALE to test both the conductive ice shell over ocean and the conductive lid over warm convective ice scenarios for a variety of conditions. The modeled crater depth-diameter is strongly dependent on thermal gradient and temperature of…
▽ More
We performed numerical simulations of impact crater formation on Europa to infer the thickness and structure of its ice shell. The simulations were performed using iSALE to test both the conductive ice shell over ocean and the conductive lid over warm convective ice scenarios for a variety of conditions. The modeled crater depth-diameter is strongly dependent on thermal gradient and temperature of the warm convective ice. Our results indicate that both a fully conductive (thin) shell and a conductive-convective (thick) shell can reproduce the observed crater depth-diameter and morphologies. For the conductive ice shell over ocean, the best fit is an approximately 8 km thick conductive ice shell. Depending on the temperature (255 - 265 K) and therefore strength of warm convective ice, the thickness of the conductive ice lid is estimated at 5 - 7 km. If central features within the crater, such as pits and domes, form during crater collapse, our simulations are in better agreement with the fully conductive shell (thin shell). If central features form well after the impact, however, our simulations suggest a conductive-convective shell (thick shell) is more likely. Although our study does not provide firm conclusion regarding the thickness of Europa's ice shell, our work indicates that Valhalla-class multiring basins on Europa may provide robust constraints on the thickness of Europa's ice shell.
△ Less
Submitted 24 November, 2017;
originally announced November 2017.
-
Effect of impact velocity and acoustic fluidization on the simple-to-complex transition of lunar craters
Authors:
Elizabeth A. Silber,
Gordon R. Osinski,
Brandon C. Johnson,
Richard A. F. Grieve
Abstract:
We use numerical modeling to investigate the combined effects of impact velocity and acoustic fluidization on lunar craters in the simple-to-complex transition regime. To investigate the full scope of the problem, we employed the two widely adopted Block-Model of acoustic fluidization scaling assumptions (scaling block size by impactor size and scaling by coupling parameter) and compared their out…
▽ More
We use numerical modeling to investigate the combined effects of impact velocity and acoustic fluidization on lunar craters in the simple-to-complex transition regime. To investigate the full scope of the problem, we employed the two widely adopted Block-Model of acoustic fluidization scaling assumptions (scaling block size by impactor size and scaling by coupling parameter) and compared their outcomes. Impactor size and velocity were varied, such that large/slow and small/fast impactors would produce craters of the same diameter within a suite of simulations, ranging in diameter from 10-26 km, which straddles the simple-to-complex crater transition on Moon. Our study suggests that the transition from simple to complex structures is highly sensitive to the choice of the time decay and viscosity constants in the Block-Model of acoustic fluidization. Moreover, the combination of impactor size and velocity plays a greater role than previously thought in the morphology of craters in the simple-to-complex size range. We propose that scaling of block size by impactor size is an appropriate choice for modeling simple-to-complex craters on planetary surfaces, including both varying and constant impact velocities, as the modeling results are more consistent with the observed morphology of lunar craters. This scaling suggests that the simple-to-complex transition occurs at a larger crater size, if higher impact velocities are considered, and is consistent with the observation that the simple-to-complex transition occurs at larger sizes on Mercury than Mars.
△ Less
Submitted 13 April, 2017;
originally announced April 2017.
-
Dynamic Sublimation Pressure and the Catastrophic Breakup of Comet ISON
Authors:
Jordan K. Steckloff,
Brandon C. Johnson,
Timothy Bowling,
H. Jay Melosh,
David Minton,
Carey M. Lisse,
Karl Battams
Abstract:
Previously proposed mechanisms have difficulty explaining the disruption of Comet C/2012 S1 (ISON) as it approached the Sun. We describe a novel cometary disruption mechanism whereby comet nuclei fragment and disperse through dynamic sublimation pressure, which induces differential stresses within the interior of the nucleus. When these differential stresses exceed its material strength, the nucle…
▽ More
Previously proposed mechanisms have difficulty explaining the disruption of Comet C/2012 S1 (ISON) as it approached the Sun. We describe a novel cometary disruption mechanism whereby comet nuclei fragment and disperse through dynamic sublimation pressure, which induces differential stresses within the interior of the nucleus. When these differential stresses exceed its material strength, the nucleus breaks into fragments. We model the sublimation process thermodynamically and propose that it is responsible for the disruption of Comet ISON. We estimate the bulk unconfined crushing strength of Comet ISON's nucleus and the resulting fragments to be 0.5 Pa and 1-9 Pa respectively, assuming typical Jupiter Family Comet (JFC) albedos. However, if Comet ISON has an albedo similar to Pluto, this strength estimate drops to 0.2 Pa for the intact nucleus and 0.6-4 Pa for its fragments. Regardless of assumed albedo, these are similar to previous strength estimates of JFCs. This suggests that, if Comet ISON is representative of dynamically new comets, then low bulk strength is a primordial property of some comet nuclei, and not due to thermal processing during migration into the Jupiter Family.
△ Less
Submitted 9 July, 2015;
originally announced July 2015.
-
A Self-Consistent Model of the Circumstellar Debris Created by a Giant Hypervelocity Impact in the HD172555 System
Authors:
B. C. Johnson,
C. M. Lisse,
C. H. Chen,
H. J. Melosh,
M. C. Wyatt,
P. Thebault,
W. G. Henning,
E. Gaidos,
L. T. Elkins-Tanton,
J. C. Bridges,
A. Morlok
Abstract:
Spectral modeling of the large infrared excess in the Spitzer IRS spectra of HD 172555 suggests that there is more than 10^19 kg of sub-micron dust in the system. Using physical arguments and constraints from observations, we rule out the possibility of the infrared excess being created by a magma ocean planet or a circumplanetary disk or torus. We show that the infrared excess is consistent with…
▽ More
Spectral modeling of the large infrared excess in the Spitzer IRS spectra of HD 172555 suggests that there is more than 10^19 kg of sub-micron dust in the system. Using physical arguments and constraints from observations, we rule out the possibility of the infrared excess being created by a magma ocean planet or a circumplanetary disk or torus. We show that the infrared excess is consistent with a circumstellar debris disk or torus, located at approximately 6 AU, that was created by a planetary scale hypervelocity impact. We find that radiation pressure should remove submicron dust from the debris disk in less than one year. However, the system's mid-infrared photometric flux, dominated by submicron grains, has been stable within 4 percent over the last 27 years, from IRAS (1983) to WISE (2010). Our new spectral modeling work and calculations of the radiation pressure on fine dust in HD 172555 provide a self-consistent explanation for this apparent contradiction. We also explore the unconfirmed claim that 10^47 molecules of SiO vapor are needed to explain an emission feature at 8 um in the Spitzer IRS spectrum of HD 172555. We find that unless there are 10^48 atoms or 0.05 Earth masses of atomic Si and O vapor in the system, SiO vapor should be destroyed by photo-dissociation in less than 0.2 years. We argue that a second plausible explanation for the 8 um feature can be emission from solid SiO, which naturally occurs in submicron silicate "smokes" created by quickly condensing vaporized silicate.
△ Less
Submitted 23 October, 2012;
originally announced October 2012.