Disruption Of Asteroids Research Articles

Asteroid disruption and deflection simulations for multi-modal planetary defense

Planetary defense from asteroids via deflective means alone does not offer viable solutions in terminal scenarios where there is little warning time before impact. The PI method of planetary defense enables operation in terminal interdiction modes where there is little warning time prior to impact, but can also operate in the same extended time scale interdiction modes as made possible by traditional deflection techniques, which results in a versatile, multi-modal planetary defense capability. The method is also practical and cost-effective since it relies solely on launch vehicles and penetrator materials already available today, and thus presents itself as a logical and competitive option for planetary defense. As per the PI method, we investigate the effectiveness of rubble pile asteroid disruption and deflection via hypervelocity impacts with 10:1 aspect ratio cylindrical tungsten penetrators. We present the results of an ongoing simulation campaign dedicated to investigating the PI method, using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian (ALE) hydrodynamics code ALE3D run with the High-End Computing Capability (HECC) at NASA Ames Research Center. We model heterogeneous rubble pile asteroids with a distribution of spherical boulders of varying initial yield strengths set within a weak binder material. We find that rubble pile asteroids of this type in the 20–100 meter-class can be effectively mitigated via 20 km/s impacts with 100–1000 kg penetrators via the coupling of the penetrator kinetic energy into the bulk material of the asteroid.

Optical and acoustic ground effects simulations from terminal defense asteroid disruption via the PI method

Our simulations suggest that Pulverize It (PI), a NASA Phase II NIAC study, is an effective multi-modal approach for planetary defense that can operate in extremely short interdiction modes (with intercepts as short as hours prior to atmospheric entry) as well as long-interdiction time scales with months to years of warning. The basic process is complete disruption of the threat via fragmentation. In scenarios with sufficiently long warning time, the fragment cloud spreads enough to miss Earth, resulting in no ground effects. In “worst-case” scenarios with short warning times (“terminal” scenarios), the fragments (typically <10 m in diameter) will enter Earth’s atmosphere, where their energy is dissipated in a series of ground-level optical pulses and de-correlated shock waves, mitigating any significant damage. We investigate the optical and acoustic ground effects through a set of simulation codes that model the interaction of asteroid fragments with Earth’s atmosphere following terminal threat interception. Even in such cases where fragments enter the atmosphere, our simulations suggest that threats mitigated by PI produce vastly less damage on the ground when compared to an equivalent unfragmented case, yielding optical energy deposition below 200 kJ/m2 and shock wave over-pressures under 3 kPa. Our simulations support the proposition that threats like 2023 PDC, the hypothetical 800 m diameter asteroid from the 2023 Planetary Defense Conference impact exercise, can be effectively mitigated through fragmentation. We find that a terminal defense mitigation scenario that disrupts 2023 PDC into 1 million fragments with an intercept of 60 days before ground impact results in reasonable ground effects with minimal damage.

Acoustic ground effects simulations from asteroid disruption via the “Pulverize It” method

Our simulations show that PI (“Pulverize It”), a NASA Phase II NIAC study, is an effective multi-modal approach for planetary defense that can operate in extremely short interdiction modes (with intercepts as short as hours prior to atmospheric entry) as well as long interdiction time scales with months to years of warning. The basic process is complete disruption of the threat via fragmentation. In long warning time cases, the fragment cloud spreads enough to miss Earth, resulting in no ground effects. In cases where the warning time is short, the fragments (typically &amp;lt;10 m in diameter) will enter Earth’s atmosphere where their energy is dissipated in a series of ground-level optical pulses and de-correlated shock waves, mitigating any significant damage. We investigate the acoustic ground effects through a set of simulation codes that model the interaction of asteroid fragments with the Earth’s atmosphere following threat interception. Even in the short warning time cases where the fragments enter the atmosphere, our simulations show that threats mitigated by the PI method produce vastly less damage on the ground when compared to the same unfragmented case, yielding shock wave over-pressures under 3 kPa. Our simulations support the proposition that threats sized 20 m—1000 m in diameter can be effectively mitigated through fragmentation, resulting in acoustic ground effects that are below estimated damage thresholds and yield short- and long-term non-lethal effects.

JWST Directly Images Giant Planet Candidates Around Two Metal-polluted White Dwarf Stars

We report the discovery of two directly imaged, giant planet candidates orbiting the metal-rich, hydrogen atmosphere white dwarfs WD 1202−232 and WD 2105−82. JWST’s Mid-Infrared Instrument (MIRI) data on these two stars show a nearby resolved source at a projected separation of 11.47 and 34.62 au, respectively. Assuming the planets formed at the same time as their host stars, with total ages of 5.3 and 1.6 Gyr, the MIRI photometry is consistent with giant planets with masses ≈1–7 M Jup. The probability of both candidates being false positives due to red background sources is approximately 1 in 3000. If confirmed, these would be the first directly imaged planets that are similar in both age and separation to the giant planets in our own solar system, and they would demonstrate that widely separated giant planets like Jupiter survive stellar evolution. Giant planet perturbers are widely used to explain the tidal disruption of asteroids around metal-polluted white dwarfs. Confirmation of these two planet candidates with future MIRI imaging would provide evidence that directly links giant planets to metal pollution in white dwarf stars.

Tidal Disruption of Near-Earth Asteroids during Close Encounters with Terrestrial Planets

Numerical modeling has long suggested that gravitationally bound (or so-called rubble-pile) near-Earth asteroids (NEAs) can be destroyed by tidal forces during close and slow encounters with terrestrial planets. However, tidal disruptions of NEAs have never been directly observed nor have they been directly attributed to any families of NEAs. Here we show population-level evidence for the tidal disruption of NEAs during close encounters with Earth and Venus. Debiased model distributions of NEA orbits and absolute magnitudes based on observations by the Catalina Sky Survey during 2005–2012 underpredict the number of NEAs with perihelion distances coinciding with the semimajor axes of Venus and Earth. A detailed analysis of the orbital distributions of the excess NEAs shows that their characteristics agree with the prediction for tidal disruptions, and they cannot be explained by observational selection effects or orbital dynamics. Accounting for tidal disruptions in evolutionary models of the NEA population partly bridges the gap between the predicted rate of impacts by asteroids with diameters of tens of meters and observed statistics of fireballs in the same size range.

Debiased population of very young asteroid families

Context. Asteroid families that are less than one million years old offer a unique possibility to investigate recent asteroid disruption events and test ideas about their dynamical evolution. Observations provided by powerful all-sky surveys have led to an enormous increase in the number of detected asteroids over the past decade. When the known populations are well characterized, they can be used to determine asteroid detection probabilities, including those in young families, as a function of their absolute magnitude. Aims. We use observations from the Catalina Sky Survey (CSS) to determine the bias-corrected population of small members in four young families down to sizes equivalent to several hundred meters. Methods. Using the most recent catalog of known asteroids, we identified members from four young families for which the population has grown appreciably over recent times. A large fraction of these bodies have also been detected by CSS. We used synthetic populations of asteroids, with their magnitude distribution controlled by a small number of parameters, as a template for the bias-corrected model of these families. Applying the known detection probability of the CSS observations, we could adjust these model parameters to match the observed (biased) populations in the young families. Results. In the case of three families, Datura, Adelaide, and Rampo, we find evidence that the magnitude distribution transitions from steep to shallow slopes near 300 to 400 meters. Conversely, the Hobson family population may be represented by a single power-law model. The Lucascavin family has a limited population; no new members have been discovered over the past two decades. We consider a model of parent body rotational fission with the escaping secondary tidally split into two components (thereby providing three members within this family). In support of this idea, we find that no other asteroid with absolute magnitude H ≤ 18.3 accompanies the known three members in the Lucascavin family. A similar result is found for the archetypal asteroid pair Rheinland–Kurpfalz.

NEOMOD: A New Orbital Distribution Model for Near-Earth Objects

Near-Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions—the main belt and trans-Neptunian scattered disk—and ends as bodies impact planets, disintegrate near the Sun, or are ejected from the solar system. Here we develop a new orbital model of NEOs by numerically integrating asteroid orbits from main-belt sources and calibrating the results on observations of the Catalina Sky Survey. The results imply a size-dependent sampling of the main belt with the ν 6 and 3:1 resonances producing ≃30% of NEOs with absolute magnitudes H = 15 and ≃80% of NEOs with H = 25. Hence, the large and small NEOs have different orbital distributions. The inferred flux of H < 18 bodies into the 3:1 resonance can be sustained only if the main-belt asteroids near the resonance drift toward the resonance at the maximal Yarkovsky rate (≃2 × 10−4 au Myr−1 for diameter D = 1 km and semimajor axis a = 2.5 au). This implies obliquities θ ≃ 0° for a < 2.5 au and θ ≃ 180° for a > 2.5 au, both in the immediate neighborhood of the resonance (the same applies to other resonances as well). We confirm the size-dependent disruption of asteroids near the Sun found in previous studies. An interested researcher can use the publicly available NEOMOD Simulator to generate user-defined samples of NEOs from our model.

Exsolution in alkali feldspar in ordinary chondrites: Ubiquitous evidence for rapid cooling at high temperatures

Thermal metamorphism in undisrupted ordinary chondrite (OC) parent bodies is thought to occur through the radioactive decay of 26Al, producing an onion-shell-like structure with higher peak metamorphic temperatures corresponding to increasing depth. During retrograde metamorphism, the onion-shell model predicts slower cooling rates with increasing petrologic type. However, cooling rates determined by pyroxene diffusion, metallographic, and other methods are inconsistent with onion-shell-like cooling, leading to a model of asteroid disruption and reaccretion into a rubble pile, after peak metamorphism. Potassium-feldspar exsolution in albite, in a perthite texture, has been noted in OCs and can be used as another method for determining cooling rates. We conducted a survey of K-feldspar occurrences and textures, within chondrules, in petrologic type 3.6–6 H, L, and LL OCs. Potassium-feldspar is present as a secondary feature, in primary and secondary albite, as fine-scale exsolution lamellae, 0.1–1.5 μm wide, as well as in larger patches up to 50 μm in size. Exsolution is present in all OC groups and is most common in petrologic type 4.In the H4 chondrite Avanhandava, we estimate the cooling rate from perthite to be 3–17 °C/yr over a temperature interval of 765–670 °C. Peristerite is also present in Avanhandava for which we estimate a cooling rate of 0.2–2.4 × 10−3 °C/yr from 570-540 °C. In general, the relatively fast, high-temperature cooling rate determined by perthite is similar to cooling rates recovered from two-pyroxene speedometry. The peristerite cooling rate is closer to the slow, lower temperature metallographic cooling rates. Because K-feldspar exsolution is present in similar fine-scale lamellae in all OC groups, we suggest that all OC parent bodies experienced the same cooling history at high temperatures. These results are inconsistent with predictions of OC asteroid cooling from undisturbed onion-shell metamorphism but are consistent with models involving disruption after peak metamorphism followed by reassembly.

The young Hobson family: Possible binary parent body and low-velocity dispersal

Context.Asteroid families with ages younger than 1 Myr offer an interesting possibility of studying the outcomes of asteroid disruptions that are little modified by subsequent evolutionary processes.Aims.We analyze a very young asteroid family associated with (18777) Hobson in the central part of the main belt. We aim at (i) understanding its peculiar size distribution, and (ii) setting an upper limit on the characteristic dispersal velocity at subkilometer sizes corresponding to the smallest visible Hobson members.Methods.We identified the Hobson family using an up-to-date asteroid catalog. A significant increase in the number of its known members allowed us to study their size distribution and compare it with computer simulations of catastrophic disruptions. Backward orbital integrations of the heliocentric orbits allowed us to confirm the previously suggested age of Hobson and helped to estimate limits of the ejection speed.Results.The Hobson family has an unusual size distribution: two nearly equal-size bodies, followed by a population of smaller asteroids, whose distribution takes a characteristic power law. There are two possibilities to explain these data. Either a canonical impact onto a single parent body, requiring fine-tuned impact conditions that have not been studied so far, or an unconventional model for the parent object of the Hobson family, namely a binary with ≃7−9 km primary and a ≃2.5 km secondary. In the latter case, the primary was disrupted, leaving behind the largest remnant (18777) Hobson and a suite of subkilometer asteroids. The second largest asteroid, (57738) 2001 UZ160, is the nearly intact satellite of the parent binary. The excellent convergence of nominal orbits of Hobson members sets an upper limit of ≃(10−20) m s−1for the initial dispersal velocity of the known members, which is consistent with both formation models. The Hobson family provides a so far rare opportunity of studying disruptions of small asteroids in a situation in which both the material strength and reaccumulation efficiency play an important role.

Minimum perihelion distances and associated dwell times for near-Earth asteroids

ABSTRACT The observed near-Earth asteroid (NEA) population contains very few objects with small perihelion distances, say, $q \lesssim 0.2\, \mathrm{au}$. NEAs that currently have orbits with larger q might be hiding a past evolution during which they have approached closer to the Sun. We present a probabilistic assessment of the minimum q that an asteroid has reached during its orbital history. At the same time, we offer an estimate of the dwell time, that is, the time q has been in a specific range. We have re-analysed orbital integrations of test asteroids from the moment they enter the near-Earth region until they either collide with a major body or are thrown out from the inner Solar system. We considered a total disruption of asteroids at certain q as a function of absolute magnitude (H). We calculated the probability that an asteroid with given orbital elements and H has reached a q smaller than a given threshold value and its respective dwell time in that range. We have constructed a look-up table that can be used to study the past orbital and thermal evolution of asteroids as well as meteorite falls and their possible parent bodies. An application to 25 meteorite falls shows that carbonaceous chondrites typically have short dwell times at small q, whereas for ordinary chondrites it ranges from $10\, 000$ to $500\, 000$ yr. A dearth of meteorite falls with long dwell times and small minimum q supports a supercatastrophic disruption of asteroids at small q.

On the Mass Distribution of Fragments of an Asteroid Disrupted in the Earth’s Atmosphere

To model the interaction with the atmosphere of fragments of a disrupted asteroid, which move independently of each other, it is necessary to know their mass distribution. In this regard, an analogy is drawn with fragmentation in high-speed impact experiments performed to simulate the disruption of asteroids at their collisions in outer space. Based on the results of impact experiments and assuming a power law for the mass distribution in a differential form, we obtained the cumulative number of fragments as a function of the fragment mass m normalized to the total mass of fragments, the mass fraction of the largest fragment(s), the number of the largest fragments, and the power index. The formula for the cumulative number of fragments of a disrupted body is used to describe the results of impact experiments for different fragmentation types. The proposed fragment mass distribution is also tested by comparison with the mass distributions of recovered meteorites in the cases of Mbale, Bassikounou, Almahata Sitta, Košice, and Chelyabinsk meteorite falls.

Accurate 3D fireball trajectory and orbit calculation using the 3D-firetoc automatic Python code

ABSTRACT The disruption of asteroids and comets produces cm-sized meteoroids that end up impacting the Earth’s atmosphere and producing bright fireballs that might have associated shock waves or, in geometrically favourable occasions excavate craters that put them into unexpected hazardous scenarios. The astrometric reduction of meteors and fireballs to infer their atmospheric trajectories and heliocentric orbits involves a complex and tedious process that generally requires many manual tasks. To streamline the process, we present a software package called SPMN 3D Fireball Trajectory and Orbit Calculator (3D-firetoc), an automatic Python code for detection, trajectory reconstruction of meteors, and heliocentric orbit computation from video recordings. The automatic 3D-firetoc package comprises of a user interface and a graphic engine that generates a realistic 3D representation model, which allows users to easily check the geometric consistency of the results and facilitates scientific content production for dissemination. The software automatically detects meteors from digital systems, completes the astrometric measurements, performs photometry, computes the meteor atmospheric trajectory, calculates the velocity curve, and obtains the radiant and the heliocentric orbit, all in all quantifying the error measurements in each step. The software applies corrections such as light aberration, refraction, zenith attraction, diurnal aberration, and atmospheric extinction. It also characterizes the atmospheric flight and consequently determines fireball fates by using the α − β criterion that analyses the ability of a fireball to penetrate deep into the atmosphere and produce meteorites. We demonstrate the performance of the software by analysing two bright fireballs recorded by the Spanish Fireball and Meteorite Network (SPMN).

Supercatastrophic Disruption of Asteroids in the Context of SOHO Comet, Fireball, and Meteor Observations

Abstract Granvik et al. reported an absence of asteroids on orbits with perihelia near the Sun that they attribute to the “supercatastrophic disruption” of these bodies. Here we investigate whether there is evidence for this process among other bodies with similarly low perihelia: near-Earth asteroids, Solar and Heliospheric Observatory (SOHO) comets, and meter- and millimeter-sized meteoroids. We determine that no known near-Earth asteroids have past (last 104 yr) histories residing significantly inside the Granvik et al. limit, indirectly supporting the disruption hypothesis. The exception is asteroid (467372) 2004 LG, which spent 2500 yr within this limit and thus presents a challenge to that theory. Phaethon has a perihelion distance hovering just above the limit and may be undergoing slow disruption, which may be the source of its dust complex. We find that the rate at which ungrouped SOHO comets are observed is consistent with expected rates for the injection of small (25 m) class asteroids into the near-Sun region and suggest that this fraction of the SOHO-observed comet population may in fact be asteroidal in origin. We also find that there is an absence of meter-sized bodies with near-Sun perihelia but an excess of millimeter-sized meteoroids. This implies that if near-Sun asteroids disrupt, they do not simply fragment into meter-sized chunks but ultimately disintegrate into millimeter-sized particles. We propose that the disruption of near-Sun asteroids, as well as the anomalous brightening and destruction processes that affect SOHO comets, occur through meteoroid erosion, that is, the removal of material through impacts by high-speed near-Sun meteoroids.

Constraining the origin of the planetary debris surrounding ZTF J0139+5245 through rotational fission of a triaxial asteroid

ABSTRACT White dwarfs containing orbiting planetesimals or their debris represent crucial benchmarks by which theoretical investigations of post-main-sequence planetary systems may be calibrated. The photometric transit signatures of likely planetary debris in the ZTF J0139+5245 white dwarf system have an orbital period of about 110 d. An asteroid which breaks up to produce this debris may spin itself to destruction through repeated close encounters with the star without entering its Roche radius and without influence from the white dwarf’s luminosity. Here, we place coupled constraints on the orbital pericentre (q) and the ratio (β) of the middle to longest semiaxes of a triaxial asteroid which disrupts outside of this white dwarf’s Roche radius (rRoche) soon after attaining its 110-d orbit. We find that disruption within tens of years is likely when β ≲ 0.6 and q ≈ 1.0–2.0rRoche, and when β ≲ 0.2 out to q ≈ 2.5rRoche. Analysing the longer time-scale disruption of triaxial asteroids around ZTF J0139+5245 is desirable but may require either an analytical approach relying on ergodic theory or novel numerical techniques.

Disassociation energies for the finite-density N-body problem

This paper considers the energy required for collections of finite density bodies to undergo escape under internal gravitational interactions alone. As the level of the system energy is increased there are different combinations of components that can escape, until the total energy becomes positive, when the entire system can undergo mutual disruption. The results are also defined for bodies modeled as a continuum. These results provide rigorous constraints for the disruption of rubble pile asteroids when only considering gravitational interaction effects, with the energy provided by rotation of an initial collection of the system. These issues are considered for discrete particles in the N body problem and for size distributions of discrete particles in the continuum limit.

Thermal and impact histories of 25143 Itokawa recorded in Hayabusa particles

Understanding the origin and evolution of near-Earth asteroids (NEAs) is an issue of scientific interest and practical importance because NEAs are potentially hazardous to the Earth. However, when and how NEAs formed and their evolutionary history remain enigmas. Here, we report the U-Pb systematics of Itokawa particles for the first time. Ion microprobe analyses of seven phosphate grains from a single particle provide an isochron age of 4.64 ± 0.18 billion years (1σ). This ancient phosphate age is thought to represent the thermal metamorphism of Itokawa’s parent body, which is identical to that of typical LL chondrites. In addition, the incorporation of other particles suggests that a significant shock event might have occurred 1.51 ± 0.85 billion years ago (1σ), which is significantly different from the shock ages of 4.2 billion years of the majority of shocked LL chondrites and similar to that of the Chelyabinsk meteorite. Combining these data with recent Ar-Ar studies on particles from a different landing site, we conclude that a globally intense impact, possibly a catastrophic event, occurred ca. 1.4 Ga ago. This conclusion enables us to establish constraints on the timescale of asteroid disruption frequency, the validity of the crater chronology and the mean lifetime of small NEAs.

Spinup and Disruption of Interstellar Asteroids by Mechanical Torques, and Implications for 1I/2017 U1 (‘Oumuamua)

Abstract The discovery of the first interstellar asteroid, 1I/2017 U1 (‘Oumuamua), has opened a new era for research on interstellar objects. In this paper, we study the rotational dynamics of interstellar asteroids (ISAs) of irregular shapes moving through the interstellar gas. We find that regular mechanical torques resulting from the bombardment of gas flow on the irregular body could be important for the dynamics and destruction of ISAs. Mechanical torques can spin up the ISA, resulting in the breakup of the original ISA into small binary asteroids when the rotation rate exceeds the critical frequency. We find that the breakup timescale is short for ISAs of highly irregular shapes and low tensile strength. We apply our results to the first observed ISA, ‘Oumuamua, and suggest that its extreme elongated shape may originate from a reassembly of the binary fragments, due to gravity along its journey in the interstellar medium. The tumbling of ‘Oumuamua could have been induced by rotational disruption due to mechanical torques. Finally, we discuss the survival possibility of high-velocity asteroids presumably formed by tidal disruption of planetary systems by the black hole at the Galactic center.

Internal gravity, self-energy, and disruption of comets and asteroids

The internal gravity and self-gravitational energy of a comet, asteroid, or small moon have applications to their geophysics, including their formation, evolution, cratering, and disruption, the stresses and strains inside such objects, sample return, eventual asteroid mining, and planetary defense strategies for potentially hazardous objects. This paper describes the relation of an object’s self-energy to its collisional disruption energy, and shows how to determine an object’s self-energy from its internal gravitational potential.Any solid object can be approximated to any desired accuracy by a polyhedron of sufficient complexity. An analytic formula is known for the gravitational potential of any homogeneous polyhedron, but it is widely believed that this formula applies only on the surface or outside of the object. Here we show instead that this formula applies equally well inside the object.We have used these formulae to develop a numerical code which evaluates the self-energy of any homogeneous polyhedron, along with the gravitational potential and attraction both inside and outside of the object, as well as the slope of its surface. Then we use our code to find the internal, external, and surface gravitational fields of the Platonic solids, asteroid (216) Kleopatra, and comet 67P/Churyumov–Gerasimenko, as well as their surface slopes and their self-gravitational energies. We also present simple spherical, ellipsoidal, cuboidal, and duplex models of Kleopatra and comet 67P, and show how to generalize our methods to inhomogeneous objects and magnetic fields.At present, only the self-energies of spheres, ellipsoids, and cuboids (boxes) are known analytically (or semi-analytically). The Supplementary Material contours the central potential and self-energy of homogeneous ellipsoids and cuboids of all aspect ratios, and also analytically the self-gravitational energy of a “duplex” consisting of two coupled spheres. The duplex is a good model for “contact binary” comets and asteroids; in fact, most comets seem to be bilobate, and might be described better as “dirty snowmen” than as “dirty snowballs”.

Rocks, Shocks and Asteroids, and Some Interesting Research Directions in Mechanics

Mechanics lies at the heart of many of the underpinnings of modern technological civilization: materials, infrastructure, transportation, health and security. The mechanics of dynamic failure processes also has a major bearing on the potential catastrophes that threaten civilization, including airbursts and major asteroid impacts. Recent events (such as the Chelyabinsk meteoroid) have demonstrated the need to understand major impact and fragmentation events. Many of the fundamental problems of current interest in national security also involve impact and fragmentation, typically studied through large-scale computational simulations. In the Murray lecture, these issues were addressed by describing fundamental high-strain-rate experiments, high-speed visualization, theoretical and computational modeling of failure processes, and simulations of asteroid damage and disruption. This paper focuses on experimental results on meteorites and a basalt.

GPU-accelerated computational tool for studying the effectiveness of asteroid disruption techniques

This paper presents the development of a new Graphics Processing Unit (GPU) accelerated computational tool for asteroid disruption techniques. Numerical simulations are completed using the high-order spectral difference (SD) method. Due to the compact nature of the SD method, it is well suited for implementation with the GPU architecture, hence solutions are generated at orders of magnitude faster than the Central Processing Unit (CPU) counterpart. A multiphase model integrated with the SD method is introduced, and several asteroid disruption simulations are conducted, including kinetic-energy impactors, multi-kinetic energy impactor systems, and nuclear options. Results illustrate the benefits of using multi-kinetic energy impactor systems when compared to a single impactor system. In addition, the effectiveness of nuclear options is observed.