Rubble Pile 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.

The Bulk Densities of Small Solar System Bodies as a Probe of Planetesimal Formation

Constraining the formation processes of small solar system bodies is crucial for gaining insights into planetesimal formation. Their bulk densities, determined by their compressive strengths, offer valuable information about their formation history. In this paper, we utilize a formulation of the compressive strength of dust aggregates obtained from dust N-body simulations to establish the relation between the bulk density and diameter. We find that this relation can be effectively approximated by a polytrope with an index of 0.5, coupled with a formulation of the compressive strength of dust aggregates. The lowest-density trans-Neptunian objects (TNOs) and main-belt asteroids (MBAs) are well reproduced by dust aggregates composed of 0.1 μm sized grains. However, most TNOs, MBAs, comets, and near-Earth asteroids (NEAs) exhibit higher densities, suggesting the influence of compaction mechanisms such as collision, dust grain disruption, sintering, or melting, leading to further growth. We speculate that there are two potential formation paths for small solar system bodies. One involves the direct coagulation of primordial dust grains, resulting in the formation of first-generation planetesimals, including the lowest-density TNOs, MBAs, and the parent bodies of comets and NEAs. In this case, comets and NEAs are fragments or rubble piles of first-generation planetesimals, and the objects themselves or the rubble are composed of 0.1 μm sized grains. The other path involves the further potential fragmentation of first-generation planetesimals into the compact dust aggregates observed in protoplanetary disks, resulting in the formation of second-generation planetesimals composed of compact dust aggregates, which may contribute to explaining another formation process of comets and NEAs.

The DECam Ecliptic Exploration Project (DEEP). VII. The Strengths of Three Superfast Rotating Main-belt Asteroids from a Preliminary Search of DEEP Data

Superfast rotators (SFRs) are small solar system objects that rotate faster than generally possible for a cohesionless rubble pile. Their rotational characteristics allow us to make inferences about their interior structure and composition. Here, we present the methods and results from a preliminary search for SFRs in the DECam Ecliptic Exploration Project (DEEP) data set. We find three SFRs from a sample of 686 main-belt asteroids, implying an occurrence rate of 0.4−0.3+0.1 %—a higher incidence rate than has been measured by previous studies. We suggest that this high occurrence rate is due to the small sub-kilometer size regime to which DEEP has access: the objects searched here were as small as ∼500 m. We compute the minimum required cohesive strength for each of these SFRs and discuss the implications of these strengths in the context of likely evolution mechanisms. We find that all three of these SFRs require strengths that are more than that of weak regolith but consistent with many cohesive asteroid strengths reported in the literature. Across the full DEEP data set, we have identified ∼70,000 Main-Belt Asteroids and expect ∼300 SFRs—a result that will be assessed in a future paper.

The geology and evolution of the Near-Earth binary asteroid system (65803) Didymos

Images collected during NASA’s Double Asteroid Redirection Test (DART) mission provide the first resolved views of the Didymos binary asteroid system. These images reveal that the primary asteroid, Didymos, is flattened and has plausible undulations along its equatorial perimeter. At high elevations, its surface is rough and contains large boulders and craters; at low elevations its surface is smooth and possesses fewer large boulders and craters. Didymos’ moon, Dimorphos, possesses an intimate mixture of boulders, several asteroid-wide lineaments, and a handful of craters. The surfaces of both asteroids include boulders that are large relative to their host body, suggesting that both asteroids are rubble piles. Based on these observations, our models indicate that Didymos has a surface cohesion ≤ 1 Pa and an interior cohesion of ∼10 Pa, while Dimorphos has a surface cohesion of <0.9 Pa. Crater size-frequency analyzes indicate the surface age of Didymos is 40–130 times older than Dimorphos, with likely absolute ages of ~\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim$$\\end{document}12.5 Myr and <0.3 Myr, respectively. Solar radiation could have increased Didymos’ spin rate leading to internal deformation and surface mass shedding, which likely created Dimorphos.

Mechanical properties of rubble pile asteroids (Dimorphos, Itokawa, Ryugu, and Bennu) through surface boulder morphological analysis

Planetary defense efforts rely on estimates of the mechanical properties of asteroids, which are difficult to constrain accurately from Earth. The mechanical properties of asteroid material are also important in the interpretation of the Double Asteroid Redirection Test (DART) impact. Here we perform a detailed morphological analysis of the surface boulders on Dimorphos using images, the primary data set available from the DART mission. We estimate the bulk angle of internal friction of the boulders to be 32.7 ± 2. 5° from our measurements of the roundness of the 34 best-resolved boulders ranging in size from 1.67–6.64 m. The elongated nature of the boulders around the DART impact site implies that they were likely formed through impact processing. Finally, we find striking similarities in the morphology of the boulders on Dimorphos with those on other rubble pile asteroids (Itokawa, Ryugu and Bennu). This leads to very similar internal friction angles across the four bodies and suggests that a common formation mechanism has shaped the boulders. Our results provide key inputs for understanding the DART impact and for improving our knowledge about the physical properties, the formation and the evolution of both near-Earth rubble-pile and binary asteroids.

Evidence for multi-fragmentation and mass shedding of boulders on rubble-pile binary asteroid system (65803) Didymos

Asteroids smaller than 10 km are thought to be rubble piles formed from the reaccumulation of fragments produced in the catastrophic disruption of parent bodies. Ground-based observations reveal that some of these asteroids are today binary systems, in which a smaller secondary orbits a larger primary asteroid. However, how these asteroids became binary systems remains unclear. Here, we report the analysis of boulders on the surface of the stony asteroid (65803) Didymos and its moonlet, Dimorphos, from data collected by the NASA DART mission. The size-frequency distribution of boulders larger than 5 m on Dimorphos and larger than 22.8 m on Didymos confirms that both asteroids are piles of fragments produced in the catastrophic disruption of their progenitors. Dimorphos boulders smaller than 5 m have size best-fit by a Weibull distribution, which we attribute to a multi-phase fragmentation process either occurring during coalescence or during surface evolution. The density per km2 of Dimorphos boulders ≥1 m is 2.3x with respect to the one obtained for (101955) Bennu, while it is 3.0x with respect to (162173) Ryugu. Such values increase once Dimorphos boulders ≥5 m are compared with Bennu (3.5x), Ryugu (3.9x) and (25143) Itokawa (5.1x). This is of interest in the context of asteroid studies because it means that contrarily to the single bodies visited so far, binary systems might be affected by subsequential fragmentation processes that largely increase their block density per km2. Direct comparison between the surface distribution and shapes of the boulders on Didymos and Dimorphos suggest that the latter inherited its material from the former. This finding supports the hypothesis that some asteroid binary systems form through the spin up and mass shedding of a fraction of the primary asteroid.

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

ABSTRACT Asteroid 2016 HO$_3$, a small asteroid (&amp;lt;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.

Rapid formation of binary asteroid systems post rotational failure: A recipe for making atypically shaped satellites

Binary asteroid formation is a highly complex process, which has been highlighted with recent observations of satellites with unexpected shapes, such as the oblate Dimorphos by the NASA DART mission and the contact binary Selam by NASA’s Lucy mission. There is no clear consensus on which dynamical mechanisms determine the final shape of these objects. In turn, we explore a formation pathway where spin-up and rotational failure of a rubble pile asteroid lead to mass-shedding and a wide circumasteroidal debris disk in which the satellite forms. Using a combination of smooth-particle hydrodynamical and N-body simulations, we study the dynamical evolution in detail. We find that a debris disk containing matter corresponding to a few percent of the primary asteroid mass extending beyond the fluid Roche limit can consistently form both oblate and bilobate satellites via a series of tidal encounters with the primary body and mergers with other gravitational aggregates. Principally, satellites end up prolate (elongated) and on synchronous orbits, accreting mainly in a radial direction while tides from the primary asteroid keep the shape intact. However, close encounters and mergers can break the orbital state, leading to orbital migration and deformation. Satellite–satellite impacts occurring in this regime have lower impact velocities than merger-driven moon formation in e.g. planetary rings, leading to soft impacts between differently sized, non-spherical bodies. The resulting post-merger shape of the satellite is highly dependent on the impact geometry. Only moons having experienced a prior mild or catastrophic tidal disruption during a close encounter with the primary asteroid can become oblate spheroids, which is consistent with the predominantly prolate observed population of binary asteroid satellites.

The dynamical origins of the dark comets and a proposed evolutionary track

So-called ‘dark comets’ are small, morphologically inactive near-Earth objects (NEOs) that exhibit nongravitational accelerations inconsistent with radiative effects. These objects exhibit short rotational periods (minutes to hours), where measured. We find that the strengths required to prevent catastrophic disintegration are consistent with those measured in cometary nuclei and expected in rubble pile objects. We hypothesize that these dark comets are the end result of a rotational fragmentation cascade, which is consistent with their measured physical properties. We calculate the predicted size-frequency distribution for objects evolving under this model. Using dynamical simulations, we further demonstrate that the majority of these bodies originated from the ν6 resonance, implying the existence of volatiles in the current inner main belt. Moreover, one of the dark comets, (523599) 2003 RM, likely originated from the outer main belt, although a JFC origin is also plausible. These results provide strong evidence that volatiles from a reservoir in the inner main belt are present in the near-Earth environment.

Numerical simulations suggest asteroids (101955) Bennu and (162173) Ryugu are likely second or later generation rubble piles

Rubble pile asteroids are widely understood to be composed of reaccumulated debris following a catastrophic collision between asteroids in the main asteroid belt, where each disruption can make a family of new asteroids. Near-Earth asteroids Ryugu and Bennu have been linked to collisional families in the main asteroid belt, but surface age analyses of each asteroid suggest these bodies are substantially younger than their putative families. Here we show, through a coupled collisional and dynamical evolution of members of these families, that neither asteroid was likely to have been created at the same time as the original family breakups, but rather are likely remnants of later disruptions of original family members, making them second, or later, generation remnants. Our model finds about 80% and 60% of asteroids currently being delivered to near-Earth orbits from the respective families of New Polana and Eulalia are second or later generation. These asteroids delivered today in the 0.5-1 km size range have median ages since their last disruption that are substantially younger than the family age, reconciling their measured crater retention ages with membership in these families.

BYORP and Dissipation in Binary Asteroids: Lessons from DART

The near-Earth binary asteroid Didymos was the target of the planetary defense demonstration mission DART in 2022 September. The smaller binary component, Dimorphos, was impacted by the spacecraft in order to measure momentum transfer in kinetic impacts into rubble piles. DART and associated Earth-based observation campaigns have provided a wealth of scientific data on the Didymos–Dimorphos binary. DART revealed the largely oblate and ellipsoidal shape of Dimorphos before the impact, while the postimpact observations suggest that Dimorphos now has a prolate shape. Here we add those data points to the known properties of small binary asteroids and propose new paradigms of the radiative binary Yarkovsky–O’Keefe–Radzievskii–Paddack (BYORP) effect as well as tidal dissipation in small binaries. We find that relatively spheroidal bodies like Dimorphos made of small debris may experience a weaker and more size-dependent BYORP effect than previously thought. This could explain the observed values of period drift in several well-characterized binaries. We also propose that energy dissipation in small binaries is dominated by relatively brief episodes of large-scale movement of (likely surface) materials, rather than long-term steady-state tidal dissipation. We propose that one such episode was triggered on Dimorphos by the DART impact. Depending on the longevity of this high-dissipation regime, it is possible that Dimorphos will be more dynamically relaxed in time for the Hera mission than it was in the weeks following the impact.

DART Impact Ejecta Plume Evolution: Implications for Dimorphos

The NASA Double Asteroid Redirection Test (DART) spacecraft impacted the moon Dimorphos of the [65803] Didymos binary system and changed the binary orbit period, demonstrating asteroid deflection by a kinetic impact and indicating that more momentum was transferred to Dimorphos by escaping impact ejecta than was incident with DART. Images of the DART impact ejecta plume were obtained by the Light Italian cubesat for Imaging of Asteroids (LICIACube) in the first few minutes after the DART impact. The ejecta plume imaged by LICIACube 158 s after the DART impact prior to closest approach shows no evidence for plume clearing at low altitude. The ejecta plume imaged 175 s after the DART impact is optically thick up to projected altitudes of 200 m above the surface of Dimorphos. These observations are compared with models of the impact ejecta plume optical depth, structure, and evolution, which are developed from point-source scaling models fitted to numerical simulations of the DART impact into a rubble pile Dimorphos with different material strengths. The observations of the impact plume optical depth and the high momentum transfer from the DART impact are not consistent with impact and ejecta plume models assuming the Dimorphos cohesive strength to be as high as 5000 Pa. Models with 5 and 50 Pa Dimorphos cohesive strength provide the overall best consistency with plume opacity observations and high momentum transfer.

Physical properties of asteroid Dimorphos as derived from the DART impact

AbstractOn 26 September 2022, NASA’s Double Asteroid Redirection Test (DART) mission successfully impacted Dimorphos, the natural satellite of the binary near-Earth asteroid (65803) Didymos. Numerical simulations of the impact provide a means to find the surface material properties and structures of the target that are consistent with the observed momentum deflection efficiency, ejecta cone geometry and ejected mass. Our simulation that best matches the observations indicates that Dimorphos is weak, with a cohesive strength of less than a few pascals, like asteroids (162173) Ryugu and (101955) Bennu. We find that the bulk density of Dimorphos ρB is lower than ~2,400 kg m−3 and that it has a low volume fraction of boulders (≲40 vol%) on the surface and in the shallow subsurface, which are consistent with data measured by the DART experiment. These findings suggest that Dimorphos is a rubble pile that might have formed through rotational mass shedding and reaccumulation from Didymos. Our simulations indicate that the DART impact caused global deformation and resurfacing of Dimorphos. ESA’s upcoming Hera mission may find a reshaped asteroid rather than a well-defined crater.

The Strength and Shapes of Contact Binary Objects

While contact binary objects are common in the solar system, their formation mechanism is unclear. In this work we examine several contact binaries and calculate the necessary strength parameters that allow the two lobes to merge without the smaller of the two being gravitationally destroyed by the larger. We find a small but nonzero amount of cohesion or a large friction angle is required for the smaller lobe to survive the merging process, consistent with observations. This means it is possible for two previously separated rubble piles to experience a collapse of their mutual orbit and form a contact binary. The necessary strength required to survive this merger depends on the relative size, shape, and density of the body, with prolate shapes requiring more cohesion than oblate shapes.

Asteroid regolith strength: Role of fine-fractions

Most smaller asteroids (<1 km diameter) are granular material loosely bound together primarily by self-gravity known as rubble piles. In an effort to better understand the evolution of rubble-pile asteroids, we performed bulk measurements using granular simulant to study the effects of the presence of fine grains on the strength of coarse grains. Our laboratory samples consisted of fine–coarse mixtures of varying percentages of fine grains by volume of the sample. We measured the material’s angle of repose, Young’s Modulus, angle of internal friction, cohesion, and tensile strength by subjecting the samples to compressive and shear stresses. The coarse grains comprising the fine–coarse mixtures ranged from 1 mm to 20 mm (2 cm) and the fines were sieved to sub-millimeter sizes (<1 mm). The measured angles of repose varied between 32°–45° which increased with increasing fine percentage. In compression, samples generally increased in strength with increasing fine percentage for both confined and unconfined environments. In all cases, the peak strengths were not for purely fine grains but for a mixture of fine and coarse grains. Shear stress measurements yielded angles of internal friction ranging between 25° and 45° with a trend opposite that of the angle of repose, 300–550 Pa for bulk cohesion, and 0.5–1.1 kPa for tensile strength. Using other published works that include data from telescopic and in-situ observations as well as numerical simulations, we discussed the implications of our findings regarding rubble-pile formation, composition, evolution, and disruption. We find that the presence of fine grains in subsurface layers of regolith on an asteroid (confined environment) aids the avoidance of disruption due to impact. However these same fines increase an asteroid’s chance to disrupt or deform from high rotation speeds due to reduced grain interlocking. In surface layers (unconfined environments), we find that the presence of fine grains between coarse ones generates stronger cohesion and aids in the prevention of mass loss and surface shedding.

Tidal dissipation of binaries in asteroid pairs

Tidal dissipation in a celestial body can be used to probe its internal structure. Tides govern the orbital evolution of binary systems and therefore constraints on the interior of binary system members can be derived by knowing the age and tidal state of the binary system. For asteroids, age estimates are challenging due to a lack of direct observation of their surface. However, the age of asteroid pairs formed by rotational fission of a parent body can be derived from dynamical modeling, and as such can be used to constrain the age of binary systems existing within asteroid pairs. We study 13 binary asteroid systems existing in asteroid pairs by modeling their tidal locking and eccentricity damping timescales from tidal dissipation in the primaries and secondaries. We consider the impact of thermal torques on these timescales from the YORP and BYORP effects. The resulting constraints on the tidal dissipation ratio Q/k2 are compared to monolithic and rubble pile asteroid theories, showing that all secondaries are consistent with rubble piles with regolith layers greater than 3 m and suggest that Q/k2 for rubble piles increases with radius. A particular case is the first bound secondary of asteroid (3749) Balam, whose Q/k2 is constrained to be between 2.7 ×104 and 1.4 ×106, consistent with a rubble-pile with a regolith thickness between 30 m and 100 m.

An Efficient Numerical Approach to Modeling the Effects of Particle Shape on Rubble-pile Dynamics

We present an approach for the inclusion of nonspherical constituents in high-resolution N-body discrete element method (DEM) simulations. We use aggregates composed of bonded spheres to model nonspherical components. Though the method may be applied more generally, we detail our implementation in the existing N-body code pkdgrav. It has long been acknowledged that nonspherical grains confer additional shear strength and resistance to flow when compared with spheres. As a result, we expect that rubble-pile asteroids will also exhibit these properties and may behave differently than comparable rubble piles composed of idealized spheres. Since spherical particles avoid some significant technical challenges, most DEM gravity codes have used only spherical particles or have been confined to relatively low resolutions. We also discuss the work that has gone into improving performance with nonspherical grains, building on pkdgrav's existing leading-edge computational efficiency among DEM gravity codes. This allows for the addition of nonspherical shapes while maintaining the efficiencies afforded by pkdgrav's tree implementation and parallelization. As a test, we simulated the gravitational collapse of 25,000 nonspherical bodies in parallel. In this case, the efficiency improvements allowed for an increase in speed by nearly a factor of 3 when compared with the naive implementation. Without these enhancements, large runs with nonspherical components would remain prohibitively expensive. Finally, we present the results of several small-scale tests: spin-up due to the YORP effect, tidal encounters, and the Brazil nut effect. In all cases, we find that the inclusion of nonspherical constituents has a measurable impact on simulation outcomes.

A slowly cooled deep crust on asteroid 4 Vesta and the recent impact history of rubble pile vestoids recorded by diogenites

In this study, we investigate the 40Ar/39Ar systematics of nineteen diogenites thought to come from deep crustal levels of asteroid 4 Vesta. We applied both Electron Backscattered Diffraction (EBSD) and 40Ar/39Ar and methods to the unbrecciated diogenite LAP 031381. We obtained three plateau ages resulting in a combined weighted mean age of 4441 ± 15 Ma (P = 0.16). The EBSD analyses suggest that LAP 031381 displays minimal evidence of shock and, when combined with petrography observations, diffusion modelling and 40Ar/39Ar data, these results suggest that the crustal volume that initially contained this diogenite, reached a temperature of ca. 630 °C at ∼ 4.44 Ga. This corresponds to a linear cooling rate of ∼ 5 °C / Ma for a crystallization age of 4550 Ma. Independent thermal models suggest that these conditions were present at a depth of 60 to 65 km at 4.44 Ga.The other eighteen diogenites yielded 40Ar/39Ar results that indicate that they have been variously shocked by impact events and seven of them yielded plateau ages ranging from 2413 ± 189 Ma to 84 ± 162 Ma. We combined these results with 40Ar/39Ar ages from eucrites and howardites and propose that the HED (Howardite, Eucrite, Diogenite) meteorites recorded impact events at the surface of Vesta until ∼ 3.4 Ga when they were then ejected during a large collision. The eucrites, diogenites and howardites were then recombined into small rubble pile asteroids which probably make up a large part of the Vestoid family. After ejection, the K/Ar system in plagioclase crystals ceased in most cases to be fully reset by impact events as the temperature spikes reached during small impacts lack enough energy to trigger significant 40Ar* diffusion. On the other hand, ultra-transient and high-temperature – sensitive pyroxene crystals kept a more systematic record of small impacts until recent time. 38Arc cosmochron cosmogenic exposure ages on diogenites mostly range from 51 ± 7 Ma to 0 ± 1 Ma and when combined with other HED cosmochron ages, suggest that almost all the HED meteorites were continuously ejected from secondary rubble pile asteroids mostly between 50 Ma and present.

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.

Investigation of the porosity of L/LL4 ordinary chondrite Bjurböle using synchrotron radiation microtomography and scanning electron microscopy: Implications for parent body evolution

Porosity is an essential property of chondritic meteorites and is closely related to the genesis, thermal evolution, metamorphism, and thermal properties of the meteorite parent bodies. We study porosity, its texture, and shapes at sub-micron resolution in 3D and 2D within a 0.35 cm3 sample of the L/LL4 ordinary chondrite (OC) Bjurböle using two techniques, synchrotron radiation microtomography (SRμCT) and scanning electron microscopy (SEM). We employ automated segmentation tools that can be applied to both SRμCT and SEM data. Successful segmentation results can be achieved by combining visual qualitative examination and machine learning algorithms.We report novel measurement results of three-dimensional porosity properties of Bjurböle, such as aspect ratio and connectivity of void spaces, and compare the results of 2D and 3D porosity analysis. The Bjurböle sample in this study is a complex, highly porous, and friable medium and the dominant type of porosity is intergranular, continuous porosity, which contains almost all porosity volume. The shapes of the void volumes have an important effect on the connectivity of the porosity and thermal transport properties. Positive correlations between void diameter and aspect ratio as well as void volume and connectivity are present in Bjurböle, which indicate that smaller voids have lower aspect ratios and lower connectivity. In Bjurböle, small, near-spherical voids with few connections have the highest relative frequency, whereas larger void spaces with higher aspect ratios and connectivity are significantly fewer. Completely isolated pores, i.e., voids surrounded by solid material, also have a high relative frequency, and they exist within the chondrules and the matrix. However, the volume percentage of these pores is negligible compared to that of the continuous porosity.Our results support the previously measured high porosities of Bjurböle. The volume percentage of intergranular void spaces, in particular in the matrix, and the measured high porosities are not in line with the results of thermal evolution and sintering models of chondritic parent bodies regarding petrologic type 4, which implies that Bjurböle originates from a parent body with an initial onion shell structure that fragmented during or after its metamorphic peak and quickly reaccreted into a rubble pile.