Rubble-pile Structure Research Articles

Expected constraints on Phobos interior from the MMX gravity and rotation observations

The origin of the Martian moons is still uncertain, and knowledge about their interior could provide support to some of its leading theories. In preparation for the JAXA Martian Moons eXploration (MMX) mission, we review our current knowledge on the interior of Phobos, and provide synthetic tests showing how the gravity and rotation determination could allow the detection of specific interior-structure properties. The inversion of the geodetic observables for the retrieval of the internal mass distribution of a set of synthetic interior models is performed via non-linear least-squares, where the interior parameterization is based on the level-set method. We additionally provide simple expressions allowing to relate some of these interior models to the geodetic observables of Phobos. The results, based on realistic measurement resolution and noise scenarios, show good retrievals for most of the models at the data resolutions expected from MMX. Specifically, we find the gravity information is realistically sufficient for the detection of mass anomalies below the Stickney crater, as well as large scale heterogeneous regions within plausible rubble-pile structures. Libration helps retrieve the more degenerate models for gravity, such as those with concentric layers or with density varying linearly with depth. The incremental improvement from further adding an hypothetical mean obliquity measurement is marginal. Finally, we apply the level-set inversion and the analytical formulas to estimate possible interior characteristics of the ‘real’ Phobos from the currently-available scarce geodetic observables. The level-set solutions for the real-data inversion generally converge to a higher mass concentration towards the surface in the equatorial region. Markov chain Monte Carlo estimations of parameters relative to a simple 2-layer model or a radial density distribution similarly hint at a lighter region inside of Phobos.

Experimental and Numerical Simulation of Ejecta Size and Velocity of Hypervelocity Impact Rubble-Pile Asteroid

Rubble-pile asteroids may be the type of near-Earth object most likely to threaten Earth in a future collision event. Small-scale impact experiments and numerical simulations for large-scale impacts were conducted to clarify the size ratio of the boulder/projectile diameter effects on ejecta size–velocity distribution. A series of small-scale impact cratering experiments were performed on porous gypsum–basalt targets at velocities of 2.3 to 5.5 km·s−1. Three successive ejection processes were observed by high-speed and ultra-high-speed cameras. The momentum transfer coefficient and cratering size were measured. A three-dimensional numerical model reflecting the random distribution of the interior boulders of the rubble-pile structure asteroid is established. The size ratio (length to diameter) of the boulder size inside the asteroid to the projectile diameter changed from 0.25 to 1.7. We conducted a smoothed particle hydrodynamics numerical simulation in the AUTODYN software to study the boulder size effect on the ejecta size–velocity distribution. Simulation results suggest that the microscopic porosity on regolith affects the propagation of shock waves and reduces the velocity of ejecta. Experiments and numerical simulation results suggest that both excavation flow and spalling ejection mechanism can eject boulders (0.12–0.72 m) out of the rubble-pile asteroid. These experiments and simulation results help us select the potential impact site in a planetary defense scenario and reduce deflection risk. are comprised primarily of boulders of a range of sizes.

Structural stability of China’s asteroid mission target 2016 HO3 and its possible structure

ABSTRACT Asteroid 2016 HO$_3$, a small asteroid (<60 m) in super fast rotation state ($\sim$28 min), and is the target of China’s Tianwen-2 asteroid sample-return mission. In this work, we investigate its structural stability using an advanced soft-sphere-discrete-element-model code, dembody, which is integrated with bonded-aggregate models to simulate highly irregular boulders. The asteroid body is numerically constructed by tens of thousands particles, and then is slowly spun up until structural failure. Rubble piles with different frictions, cohesions, morphologies, grain size distributions, and structures are investigated. We find a 2016 HO$_3$ shaped granular asteroid would undergo tensile failure at higher strengths as opposed to shear failure in lower strengths, regardless of its shape and constituent grain size ratio. In the tensile failure regime, the critical tensile strength is proportional to the square of the spin rate, but surprisingly, is independent of the internal friction angle. Such relations indicate that the Maximum Tensile Stress criterion emerges as superior paradigm for investigating the failure behaviour of fast-rotating asteroids. We predict that the high-spin rate of asteroid 2016 HO$_3$ requires a surface strength over $\sim$3 Pa and a bulk tensile strength over $\sim$10–30 Pa. Through comparing these strength conditions with the latest data from asteroid missions, we suggest a higher likelihood of a monolithic structure over a typical rubble pile structure. However, the possibility of the latter cannot be completely ruled out. In addition, the asteroid’s surface could retain a loose regolith layer globally or only near its poles, which could be the target for sampling of Tianwen-2 mission.

Modelling the impact of DART on the asteroid Dimorphos reveals its rubble-pile structure

Modelling the impact of DART on the asteroid Dimorphos reveals its rubble-pile structure

Lessons Learned from NASA’s DART Impact about Disrupting Rubble-pile Asteroids

We present a series of numerical simulations using a shock physics smoothed particle hydrodynamics code, investigating energetic impacts on small celestial bodies characterized by diverse internal structures, ranging from weak and homogeneous compositions to rubble-pile structures with varying boulder volume packing. Our findings reveal that the internal structure of these rubble-pile bodies significantly influences the impact outcomes. Specifically, we observe that the same impact energy can either catastrophically disrupt a target with a low boulder packing (≲30 vol%), or result in the ejection of only a small fraction of material from a target with the same mass but high boulder packing (≳40 vol%). This finding highlights the pivotal role played by the rubble-pile structure, effectively acting as a bulk shear strength, which governs the size and behavior of the resulting impact. Consequently, understanding and characterizing the internal structure of asteroids will be of paramount importance for any future efforts to deflect or disrupt an asteroid on a collision course with Earth.

Exploring density and strength variations in asteroid 16 Psyche’s composition with 3D hydrocode modeling of its deepest impact structure

Asteroid 16 Psyche is the largest metallic Main Belt Asteroid and is the subject of a forthcoming NASA mission. The composition of Psyche is still unknown and subject of recent debate. In particular, how much porosity is within Psyche, along with how much of Psyche consists of non-metallic versus metallic materials, are central questions to the issue of Psyche’s composition. If Psyche is indeed predominantly composed of metallic materials, it would need to have considerable porosity (∼ 30%–50%) for a composition consistent with its expected bulk density (∼ 3.7–4.1 g/cm3). In this work, we vary the density and strength of Psyche by including uniform and layered fields of pseudo-microporosity, in addition to investigating the presence of macroscopic voids, i.e., spaces larger than the size of the simulation’s mesh cells, in rubble-pile configurations. All configurations result in bulk densities within the uncertainties of measured values, however the strength of Psyche and the distribution of pseudo-pores are varied. Through 3D computational models of Psyche’s deepest impact structure, we show that Psyche’s composition is unlikely to contain only pseudo-microporosity. Rather, rubble pile structures, which include macroscopic voids, are shown to match the crater’s measured aspect ratio better than simulations of structures that included only pseudo-microporosity.

The Formation of a Rubble Pile Asteroid: Insights from the Asteroid Ryugu

The Hayabusa2 mission returned primitive samples from the C-type asteroid Ryugu to Earth. The C-type asteroids hold clues to the origin of Earth’s water and the building blocks of life. The rubble pile structure of C-type asteroids is a crucial physical feature relating to their origin and evolution. A rubble pile asteroid is hypothesized to be bound primarily by self-gravity with a significant void space among irregularly shaped materials after catastrophic impacts between larger asteroids. However, the geological observations from Hayabusa2 and the analyses of the returned sample from Ryugu revealed that the high microporosity was common to various >10 m- to mm-sized materials of Ryugu, which suggests that the asteroid Ryugu is not just a loosely bound agglomeration of massive rocky debris from shattered asteroids. For a better understanding of the origin and evolution of the rubble pile asteroid, the current most accepted hypothesis should be verified by observations and laboratory analyses and improved upon based on this information. Here, the previous models are examined using Hayabusa2’s geological observations of the asteroid and the analytical data from the samples returned from Ryugu’s surface and subsurface material. Incorporating the new findings, a hypothesis for the evolution of the rubble pile asteroid Ryugu from a cometary nucleus through sublimation and subsequent dynamic resurfacing is proposed. The proposed hypothesis is applicable to other rubble-pile asteroids and would provide perspectives for near-Earth objects in general.

Simulating hypervelocity impacts into rubble pile structures for planetary defense

Kinetic impactors are one of the most mature technologies for planetary defense and this technique was tested for the first time by NASA's Double Asteroid Redirection Test (DART) mission. The material and physical properties of the asteroid are very important in understanding the efficacy of a kinetic impactor. To constrain how the local topography of an impactor alters the response to the hypervelocity impact of a kinetic impactor, we used the adaptive smoothed particle hydrodynamics code, Spheral++, to impact plates into a rubble pile asteroid. By varying both the impact location and the strength between the boulders and regolith of the formed asteroid, we find that variations in the rubble pile contributes approximately 10% uncertainty to the resulting momentum enhancement in regolith dominated impacts. We further demonstrate that the craters produced from the impacts into the rubble pile are wide, shallow, and show deviations from symmetry due to the presence of boulders on the surface near the impact site. These simulations help us constrain the potential outcomes that may occur from a hypervelocity kinetic impact in a planetary defense scenario.

Rubble pile asteroids are forever

Rubble piles asteroids consist of reassembled fragments from shattered monolithic asteroids and are much more abundant than previously thought in the solar system. Although monolithic asteroids that are a kilometer in diameter have been predicted to have a lifespan of few 100 million years, it is currently not known how durable rubble pile asteroids are. Here, we show that rubble pile asteroids can survive ambient solar system bombardment processes for extremely long periods and potentially 10times longer than their monolith counterparts. We studied three regolith dust particles recovered by the Hayabusa space probe from the rubble pile asteroid 25143 Itokawa using electron backscatter diffraction, time-of-flight secondary ion mass spectrometry, atom probe tomography, and 40Ar/39Ar dating techniques. Our results show that the particles have only been affected by shock pressure of ca. 5 to 15 GPa. Two particles have 40Ar/39Ar ages of 4,219 ± 35 and 4,149 ± 41 My and when combined with thermal and diffusion models; these results constrain the formation age of the rubble pile structure to ≥4.2 billion years ago. Such a long survival time for an asteroid is attributed to the shock-absorbent nature of rubble pile material and suggests that rubble piles are hard to destroy once they are created. Our results suggest that rubble piles are probably more abundant in the asteroid belt than previously thought and provide constrain to help develop mitigation strategies to prevent asteroid collisions with Earth.

Is (3599) Basov a large C-type super-fast rotator?

The rubble-pile model of asteroids deduced from the existence of a spin cutoff is an important paradigm in the interior structure of asteroids as a consequence of their collisional evolution. Many objects with size less than 200 ​m have been found to be rotating with a spin period less than 2.2 ​h but seldom so for bigger objects. Thus, asteroids with rotation periods significantly less than this spin cutoff limit are called super-fast rotators. Most often this picture of spin cutoff and hence the rubble-pile model has been applied to the S-type asteroids. In this study, we report the interesting case of (3599) Basov which is a member of the Themis family of C-class composition. The spin barrier of the Themis family asteroids has been shown to be P ​= ​3 ​h. What is unusual about (3599) Basov is that its rotation period has been confirmed to be P ​= ​2.157 ​± ​0.037 ​h which is significantly shorter than the critical value for centrifugal breakup. Furthermore, its size of 18 ​km is much larger than those of the S-class super-fast rotators. The physical property of (3599) Basov might therefore be quite unique. • C-type asteroids of rubble-pile structure should have a spin-cutoff at about 3 h • The rotation period of (3599) Basov is confirmed to be P=2.157±0.037 h and a C-type asteroid by photometry and spectroscopy observations. • If (3599) Basov has a rubble-pile structure, it should not have such a short rotation period unless it has a very large cohesive force.

Reshaping and ejection processes on rubble-pile asteroids from impacts

Context. Most small asteroids (< 50 km in diameter) are the result of the breakup of a larger parent body and are often considered to be rubble-pile objects. Similar structures are expected for the secondaries of small asteroid binaries, including Dimorphos, the smaller component of the 65 803 Didymos binary system and the target of NASA’s Double Asteroid Redirection Test (DART) and ESA’s Hera mission. The DART impact will occur on September 26, 2022, and will alter the orbital period of Dimorphos around Didymos. Aims. In this work we assume Dimorphos-like bodies with a rubble-pile structure and quantify the effects of boulder packing in its interior on the post-impact morphology, degree of shape change, and material ejection processes. Methods. We used the Bern smoothed particle hydrodynamics shock physics code to numerically model hypervelocity impacts on small, 160 m in diameter, rubble-pile asteroids with a variety of boulder distributions. Results. We find that the post-impact target morphology is most sensitive to the mass fraction of boulders comprising the target, while the asteroid deflection efficiency depends on both the mass fraction of boulders on the target and on the boulder size distribution close to the impact point. Our results may also have important implications for the structure of small asteroids.

The Coupling Orbit–Attitude–Structure Evolution of Rubble-Pile Asteroid with Earth Flyby in the Restricted Three-Body Problem

Some asteroids flying close to Earth may pose a threat of impact. Among them, the structural and dynamical characteristics of rubble-pile asteroids can be changed because of the tidal force of the Earth in this process. This can provide key information for predicting the dynamical evolution of potentially hazardous asteroids. In this study, the long-term evolution of the coupling orbit–attitude–structure of these small bodies is presented numerically based on the integration of two models. One is the 3D discrete element method, which models the structure and irregular shape of the rubble-pile asteroid. The other is the dynamical model of the circular restricted three-body problem (CRTBP). This provides a more precise dynamical environment of the asteroid orbital deflection, morphological modification, and attitude angles analysis compared to the frequently adopted two-body problem. Parametric studies on the asteroid evolution were performed focusing on its flyby distance and the bulk porosity. Numerical results indicate that the Earth flyby can form different patterns of modification of asteroids, where the rubble-pile structure can be destructed by considering the bulk porosity. The asteroid orbital deflection and attitude variational trends are also summarized based on the simulations of multi-orbital revolutions.

The Large Superfast Rotators Discovered by the Zwicky Transient Facility

Two high-cadence surveys aiming for rotation period measurements of asteroids have been conducted in 2019 January and October using the Zwicky Transient Facility. From the surveys, 25 large superfast rotators (SFRs) were discovered and they are all main-belt asteroids (MBAs), except for one Mars crosser. These large SFRs have a diameter ranging from 0.43 to 7.87 km and a rotation period between 0.48 and 1.95 hr. Considering their diameters and fast rotations, they cannot be explained by rubble-pile structure unless using extraordinary high bulk densities. Cohesion, if available, can conserve these large SFRs. The estimated cohesion for these large SFRs could be up to thousands of pascals, much higher than the currently estimated cohesion for asteroids and that generated by the regolith of Moon and Mars. Such high-level cohesion can be produced from fine-grain regolith, like clay. However, the availability of such fine-grain regolith for asteroids is still unknown. Although the possibility of these large SFRs being large monolithic objects cannot be ruled out, this scenario is very unlikely given that the timescale of disruptive impact for MBAs in a similar diameter range is 107–108 yr.

FOSSIL. II. The Rotation Periods of Small-sized Hilda Asteroids

Abstract Using the high-cadence lightcurves collected from the FOSSIL survey, rotation periods of 17 small (diameter 1 km < D < 3 km) Hilda asteroids (hereinafter Hildas) were obtained. Combined with the previously measured rotation periods of Hildas, a spin-rate limit appears at around 3 hr. Assuming rubble-pile structures for the Hildas, a bulk density of ∼1.5 g cm−3 is required to withstand this spin-rate limit. This value is similar to that of the C-type asteroids (1.33 g cm−3) and higher than the ∼1 g cm−3 bulk density of the Jupiter Trojans. This suggests that the Hildas population may contain more C-type asteroids than expected, and the limit at 3 hr simply reflects the spin-rate limit for C-type asteroids. In addition, a Hilda superfast rotator was found, which has a rotation period of 1.633 hr and an estimated diameter of 0.7 km. This object is unlikely to be explained by a rubble-pile or monolithic structure.

Interior of top-shaped asteroids with cohesionless surface

Multiple sources of evidence suggest that asteroids ranging from hundreds meters to few kilometers in size are rubble piles, i.e. gravitational aggregates of loosely consolidated material. However, no direct data on their internal structure is available to date. Cohesion between rubble-pile building blocks has been invoked in the past to explain the stability of top-shaped asteroids, which in most cases would not be capable of maintaining their large-scale shape features (low flattening, and a pronounced equatorial ridge) otherwise. However, the physical origin of cohesion is unclear and there is no direct evidence of it. Recent close-range imaging and local sampling of the surfaces of top-shaped Near Earth Asteroids (NEA) suggest the presence of very porous surface structure with minimal strength and nearly no cohesion. This raises new questions about the internal structure of such objects, with important implications on their origin and evolution. Here we show by numerically simulating the dynamics of irregular rocky fragments, that the presence of a rigid core within the asteroid’s rubble-pile structure can explain the top shape and surface features observed recently on Bennu and Ryugu, without the need of cohesion between building blocks. Also, we find that the rigid core model produces more easily equatorial mass shedding, which is thought to be responsible for satellite formation. The presence of a rigid core has never been revealed so far, but is consistent with the accretion history of those objects, and with recent estimates of their internal mass distribution. Our findings will be tested directly by ESA’s Hera mission, which will scan the interior of Dimorphos, the small moon of Didymos binary system, providing for the first time direct data on the interior of a NEA.

The Asteroid 162173 Ryugu: a Cometary Origin

Abstract The Japanese Hayabusa2 mission has revealed in detail the physical characteristics of the C-type asteroid 162173 Ryugu, in particular, its spinning top-shaped rubble-pile structure and potentially high organic content. A widely accepted formation scenario for Ryugu is catastrophic collision between larger asteroids and the subsequent slow gravitational accumulation of collisional debris. An alternative scenario is that Ryugu is an extinct comet that lost its icy components. Here, we numerically simulated the sublimation of water ice from a porous cometary nucleus until the refractory components, such as silicate rocks and organic matter, were left behind as evaporative residues. Such a process represents the transformation from a comet to an asteroid. The spin-up related to the shrinking nucleus, associated with water ice sublimation, was also calculated. The result of the calculation indicates that the cometary origin scenario can account for all the features of Ryugu discussed above. We conclude that organic-rich spinning top-shaped rubble-pile asteroids, such as Ryugu, are comet–asteroid transition objects or extinct comets.

Size Distributions of Bluish and Reddish Small Main-belt Asteroids Obtained by Subaru/Hyper Suprime-Cam*

Abstract We performed a wide-field survey observation of small asteroids using the Hyper Suprime-Cam installed on the 8.2 m Subaru Telescope. We detected more than 3000 main-belt asteroids with a detection limit of 24.2 mag in the r-band, which were classified into two groups (bluish C-like and reddish S-like) by the g–r color of each asteroid and obtained size distributions of each group. We found that the shapes of the size distributions of asteroids with C-like and S-like colors agree with each other in the size range of 0.4–5 km in diameter. Assuming the asteroid population in this size range is under collision equilibrium, our results indicate that compositional difference hardly affects the size dependence of impact strength, at least for the size range between several hundred meters and several kilometers. This size range corresponds to the size range of “spin barrier,” an upper limit observed in the rotation rate distribution. Our results are consistent with the view that most asteroids in this size range have a rubble-pile structure.

FOSSIL. I. The Spin Rate Limit of Jupiter Trojans

Abstract Rotation periods of 53 small (diameters 2 km < D < 40 km) Jupiter Trojans (JTs) were derived using the high-cadence lightcurves obtained by the FOSSIL phase I survey, a Subaru/Hyper Suprime-Cam intensive program. These are the first reported periods measured for JTs with D < 10 km. We found a lower limit of the rotation period near 4 hr, instead of the previously published result of 5 hr found for larger JTs. Assuming a rubble-pile structure for JTs, a bulk density of ≈0.9 g cm−3 is required to withstand this spin rate limit, consistent with the value ∼0.8–1.0 g cm−3 derived from the binary JT system, (617) Patroclus–Menoetius system.

(216) Kleopatra, a low density critically rotating M-type asteroid

Context. The recent estimates of the 3D shape of the M/Xe-type triple asteroid system (216) Kleopatra indicated a density of ~5 g cm−3, which is by far the highest for a small Solar System body. Such a high density implies a high metal content as well as a low porosity which is not easy to reconcile with its peculiar “dumbbell” shape. Aims. Given the unprecedented angular resolution of the VLT/SPHERE/ZIMPOL camera, here, we aim to constrain the mass (via the characterization of the orbits of the moons) and the shape of (216) Kleopatra with high accuracy, hence its density. Methods. We combined our new VLT/SPHERE observations of (216) Kleopatra recorded during two apparitions in 2017 and 2018 with archival data from the W. M. Keck Observatory, as well as lightcurve, occultation, and delay-Doppler images, to derive a model of its 3D shape using two different algorithms (ADAM, MPCD). Furthermore, an N-body dynamical model allowed us to retrieve the orbital elements of the two moons as explained in the accompanying paper. Results. The shape of (216) Kleopatra is very close to an equilibrium dumbbell figure with two lobes and a thick neck. Its volume equivalent diameter (118.75 ± 1.40) km and mass (2.97 ± 0.32) × 1018 kg (i.e., 56% lower than previously reported) imply a bulk density of (3.38 ± 0.50) g cm−3. Such a low density for a supposedly metal-rich body indicates a substantial porosity within the primary. This porous structure along with its near equilibrium shape is compatible with a formation scenario including a giant impact followed by reaccumulation. (216) Kleopatra’s current rotation period and dumbbell shape imply that it is in a critically rotating state. The low effective gravity along the equator of the body, together with the equatorial orbits of the moons and possibly rubble-pile structure, opens the possibility that the moons formed via mass shedding. Conclusions. (216) Kleopatra is a puzzling multiple system due to the unique characteristics of the primary. This system certainly deserves particular attention in the future, with the Extremely Large Telescopes and possibly a dedicated space mission, to decipher its entire formation history.

A lighter core for Phobos?

Context.The origin of the Martian moons Phobos and Deimos is still poorly understood, and is the focus of intense debate.Aims.We demonstrate that a stratified internal structure of Phobos is compatible with the observed gravity coefficients.Methods.We fit previously derivedC20andC22Phobos gravity coefficients derived from the combined MEX Doppler-tracking data from the close flybys in +2010 and 2013 with respect to the corresponding coefficients of a core–mantle stratification model of Phobos, with two opposite cases: a core denser than the mantle, and a core lighter than the mantle.Results.Only the case with a core lighter than the mantle fits at the 3σlevel the previously reported observed second degree and order coefficientC20, but a homogeneous Phobos cannot be strictly ruled out at the 3σlevel.Conclusions.This possible loosening of the core density might be the result of a displacement of material toward the surface, may be caused by centrifugal forces acting on a loosely packed rubble-pile structure, and/or by a hot-then-cold in-orbit accretion process. These two hypotheses are by no means exhaustive.