A thor asteroid family as the source of the EL enstatite meteorites

Chrysa Avdellidou, Marco Delbo, Morbidelli Alessandro, K. J Walsh, E.

Munaibari, J. Bourdelle de Micas, M. Devogèle, S. Fornasier, Matthieu

Gounelle, G. van Belle

► To cite this version:

Chrysa Avdellidou, Marco Delbo, Morbidelli Alessandro, K. J Walsh, E. Munaibari, et al. Athor asteroid family as the source of the EL enstatite meteorites. Astronomy and Astrophysics - A&A, 2022, 665, pp.L9. 10.1051/0004-6361/202244590. hal-03803920

HAL Id: hal-03803920

https://hal.science/hal-03803920

Submitted on 7 Oct 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
LETTER TO THE EDITOR

# Athor asteroid family as the source of the EL enstatite meteorites

C. Avdellidou (X. Αβδελλίδου)<fnref n="1" />, M. Delbo<fnref n="1" />, A. Morbidelli<fnref n="1" />, K. J. Walsh<fnref n="2" />, E. Munaibari<fnref n="1" />,<fnref n="3" />, J. Bourdelle de Micas<fnref n="4" />, M. Devogèle<fnref n="5" />, S. Fornasier<fnref n="4" />,<fnref n="6" />, M. Gounelle<fnref n="7" />, and G. van Belle<fnref n="8" />

Received 24 July 2022 / Accepted 30 August 2022

ABSTRACT

The identification of meteorite parent bodies provides the context for understanding planetesimal formation and evolution as well as the key Solar System events they have witnessed. However, identifying such links has proven challenging and some appear ambiguous. Here, we identify that the family of asteroid fragments whose largest member is (161) Athor is the unique source of the rare EL enstatite chondrite meteorites, the closest meteorites to Earth in terms of their isotopic ratios. The Athor family was created by the collisional fragmentation of a parent body 3 Gyr ago in the inner main belt. We calculate that the diameter of the Athor family progenitor was 64 km in diameter, much smaller than the putative size of the EL original planetesimal. Therefore, we deduce that the EL planetesimal that accreted in the terrestrial planet region underwent a first catastrophic collision in that region, and one of its fragments suffered a more recent catastrophic collision in the main belt, generating the current source of the EL meteorites.

**Key words.** minor planets, asteroids: general – meteorites, meteors, meteoroids – techniques: spectroscopic – catalogs

# 1. Introduction

It is established that meteorites were formed inside planetes-
imals that accreted in our protoplanetary disc 4.5 Gyr ago
(Johansen et al. 2015), while current main belt asteroids are
derived from the leftovers of the aforementioned planetesi-
mals. In the main belt, subsequent collisional evolution over
4.5 Gyr further destroyed some planetesimals by breaking them
to create families of smaller asteroid fragments (Nesvorny et al.
2015), and therefore the few surviving planetesimals currently
coexist with millions of asteroid fragments (Delbo et al. 2017,
2019). These small asteroid fragments suffer orbital mobility due
to non-gravitational forces (Vokrouhlický et al. 2015), which
can drive them into zones of orbital instability, where their
orbital eccentricity typically increases, leading to Earth-crossing
orbits (Bottke et al. 2006). These fragments can therefore impact
Earth, producing meteorites. Meteorites, which are the debris of
a primary parent asteroid that broke and produced a family of
smaller fragments (secondary parents; Greenwood et al. 2020),
can come directly from a main belt source. Another pathway is
that whereby the secondary parent asteroids (the family mem-
bers), via the aforementioned orbital instability zones, become
near-Earth asteroids (NEAs), from which debris can detach and
produce meteorites.

According to these scenarios, a meteorite or a meteorite type
is putatively linkable to a family of secondary parent asteroids,

which themselves can be linked to a primary parent asteroid
(family progenitor), and ultimately to an original planetesimal
(Greenwood et al. 2020). Establishing this linkage allows us to
gain insight into planetesimal composition, original size, forma-
tion time, and, indirectly, the heliocentric distance of its accre-
tion zone (Henke et al. 2012; Trieloff et al. 2022). The heliocen-
tric distance can be inferred from the isotopic abundance ratio of
some elements in meteorites, whereas the planetesimal size can
be estimated by fitting size-dependent thermal evolution mod-
els to meteorite cooling curves determined by the isotopic ther-
mochronometers (Henke et al. 2012; Trieloff et al. 2022).

The linkage between meteorites and parent bodies requires
several conditions to be met: (i) the putative parent should be
associated with a family of collisional fragments, and (ii) the
reflectance spectra of family members, their albedo, and colours
should be consistent with those measured for the meteorites,
considering space weathering alteration (Vernazza et al. 2009;
Brunetto et al. 2015). Further supporting information may come
from the pre-atmospheric orbit of the meteoroid that produced
the meteorite (Spurný et al. 2003; Jenniskens et al. 2009) and the
cosmic-ray exposure times (CRE) of similar meteorites, indicat-
ing common origin from their secondary parent bodies.

The most undisputed linkage is between the Howardites,
Eucrites, and Diogenites (HEDs) meteorites, the Vesta asteroid
family (Vestoids), and the asteroid (4) Vesta itself (Russell et al.
2012), with the latter being the primary parent as well as the

<footnote n="1">¹ Université Côte d'Azur, CNRS – Lagrange, Observatoire de la Côte d'Azur, CS 34229, 06304 Nice Cedex 4, France
e-mail: chrysa.avdellidou@oca.eu</footnote>

<footnote n="2">² Southwest Research Institute, 1050 Walnut St. Suite 300, Boulder, CO 80302, USA</footnote>

<footnote n="3">³ Université Côte d'Azur, CNRS – Géoazur, Observatoire de la Côte d'Azur, 250 rue Albert Einstein, Sophia Antipolis, 06560 Valbonne, France</footnote>

<footnote n="4">⁴ LESIA, Université Paris Cité, Observatoire de Paris, Université PSL, Sorbonne Université, CNRS, 92190 Meudon, France</footnote>

<footnote n="5">⁵ Arecibo Observatory, University of Central Florida, HC-3 Box 53995, Arecibo, PR 00612, USA</footnote>

<footnote n="6">⁶ Institut Universitaire de France (IUF), 1 rue Descartes, 75231 Paris Cedex 05, France</footnote>

<footnote n="7">⁷ Muséum National d'Histoire Naturelle, Sorbonne Universités, CNRS, IMPMC – UMR CNRS 7590, 57 rue Cuvier, 75005 Paris, France</footnote>

<footnote n="8">⁸ Lowell Observatory, 1400 West Mars Hill Road, Flagstaff, AZ 86001, USA</footnote>
surviving planetesimal. Moreover, the high-albedo Hungaria family is associated with the enstatite achondrite meteorites (aubrites; Lucas et al. 2019). Additionally, broader associations have been made between the asteroids belonging to the spectroscopic S and C complexes and the ordinary and carbonaceous chondrite meteorites, respectively (Reddy et al. 2015; DeMeo et al. 2022). However, the S- and C-complex main belt asteroid families are numerous (Nesvorny et al. 2015), complicating the search for unique linkages between the different subtypes of ordinary and carbonaceous chondrites and their respectively sources.

On the other hand, the search for meteorites linkable to X-complex families benefits from the advantage that (i) such families are far less numerous, and (ii) the different meteorite types with reflectance spectra within the X-complex, such as the enstatite chondrites (ECs; Vernazza et al. 2009), the aubrites (enstatite achondrites; Gaffey et al. 1992), some types of carbonaceous chondrites (DeMeo et al. 2015), the mesosiderites (Vernazza et al. 2009), and the metallic meteorites (Cloutis et al. 1990), can be distinguished on the basis of weak spectral features and absolute reflectance values (related to albedo). Specifically, the ECs correspond to 2% of all chondrites (Norton & Chitwood 2008) and are divided into the higher and lower iron content EH and EL classes respectively, which contribute equal numbers of meteorites to our collection. Both EHs and ELs have similar oxygen isotopic compositions (Clayton & Mayeda 1996; Herwartz et al. 2014), yet some incompatibilities have been emphasised. In particular, Lin & El Goresy (2002) showed that the chemical and petrography of the EHs and ELs are incompatible with a common parent body, where EHs have a higher Mg/Si ratio than the ELs (Keil 1968). EHs and ELs also have very distinct cooling rates and heating histories, suggesting distinct parents (Zhang et al. 1995; Lin & El Goresy 2002). Each class can be subdivided into different petrologic types that experienced specific degrees of thermal metamorphism: for example, type 6 chondrites were heated to higher metamorphic peak temperatures than those of type 3. ECs are all completely dry as no ECs of petrologic types 1 or 2 have been found. A possible model for the early evolution of their parents is an onion-shell-layered body (Trieloff et al. 2003, 2022). On the other hand, the enstatite achondrites formed from EC-like material on at least two separate parent bodies; one indicated from the peculiar Shallowater meteorite and another one for all the other aubrites. They are nearly monomineralic enstatite pyroxenites, consisting mostly of nearly FeO-free enstatite and other rare minerals that formed under highly reducing conditions (Keil 2010).

Here we study the ELs because the new discovery that thermo-chronometers inside them define the cooling history of the original EL planetesimal allows us to constrain its size (Trieloff et al. 2022). Therefore, identifying the main belt source and parenthood of the ELs offers a unique opportunity to study the formation and evolution of the EL original planetesimal.

## 2. Global search of the EL parent body

There are 15 putative X-complex families (Nesvorny et al. 2015; Delbo et al. 2019) and, as a first criterion, we investigated which of them have geometric albedo values ($p_V$) compatible with those of the ELs. To do so, first we retrieved all the EL spectra from the Reflectance Experiment Laboratory (RELAB) and the Planetary Spectrophotometer Facility (PSF) meteorite spectra databases and using the method of Beck et al. (2021) we calculated their average geometric visible albedo $p_V$ and its standard deviation, finding a value between 0.11 and 0.23 (Fig. 1B).

For each of these asteroid families, we then retrieved from the Minor Planet Physical Properties Catalogue the $p_V$-values of their members, and computed their average (Table A.1). Twelve families can be excluded as the source of the ELs because their average $p_V$-values are out of the range of the ELs (Fig. 1B). Finally, the San Marcello family has an average $p_V = 0.177$ that is within the range of the ELs (original membership from Nesvorny et al. 2015), but we excluded it as a source of ELs because (i) Broz et al. (2022) proposed that asteroid (22) Kalliope is the potential progenitor of the family, which could be a differentiated body; (ii) the $p_V$-values of its members span a large range between 0.1 and 0.35 (Broz et al. 2022), which could also indicate the presence of a large number of interlopers; and (iii) the position of the San Marcello family in the outer main belt does not favour the delivery of meteorites. The remaining two families with $p_V$-values compatible with the ELs, namely Athor and Zita, are located in the inner main belt, and are therefore dynamically favoured to deliver meteorites (Delbo et al. 2019).

## 3. Spectroscopic observations of potential EL sources

In order to study the potential link between the Athor and Zita families and ELs, we carried out a spectroscopic observational campaign (Appendix B.1, Table B.1), characterised the asteroid reflectances and their spectroscopic types (Appendix B.2, Table B.2), and re-assessed the family membership (Table B.2). Observations were performed mainly in the near-infrared wavelength range using the SpeX instrument at the NASA Infrared Telescope Facility, while in some cases we also performed visible spectroscopic observations using the Copernico Telescope and the Lowell Discovery Telescope (Table B.1). We combined our results with the literature data (Fig. B.1).

In total, we observed 39 members of both families (Appendix B.1) according to the family membership of Delbo et al. (2019), and found a bimodal distribution of their near-infrared spectral reflectance slopes with the two groups clearly separating at $15\% \ \mu\text{m}^{-1}$ (Fig. B.2). These groups have mean and standard deviation of $8.3 \pm 3\% \ \mu\text{m}^{-1}$ and $21 \pm 2.3\% \ \mu\text{m}^{-1}$, respectively (Fig. B.2), being separated by more than $4\sigma$. Even if we adopt the $4.2\% \ \mu\text{m}^{-1}$ SpeX statistical uncertainty of reflectance calibration (Marsset et al. 2020) as standard deviation, the centres of the two groups are separated by more than $3\sigma$. We deduced that the two groups are real, and therefore Athor and Zita families can be distinguished in the near-infrared. The first group predominantly occurs in Athor family and the second in Zita.

Next, we classified each asteroid reflectance in the Bus-DeMeo taxonomy (DeMeo et al. 2009) following established methods (Appendix B.2, Table B.2). The spectrum of the Xc-type is the most recurrent within the Athor-type group. Where applicable, the Xe- and Ch/Cgh-type matches were not preferred as previous dedicated studies showed that Xe-type asteroids have very high $p_V$-values (Lucas et al. 2019), while Ch/Cgh-type have low $p_V$-values, both incompatible with those of Athor family members (Table B.2). In the case of the Zita-type group, the best spectral matches were with Xk-, X-, and T-types.

### 3.1. Re-assessment of family membership

Given our spectroscopic results, we proceeded to a reassessment of the family membership. This was needed because the identification of family members on the basis of clustering in
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Reflectance spectrum of the Athor family and meteorites. (A) Average spectra of Athor and Zita family members, which are clearly distinct, while the curve-matching of EL spectra with the average Athor family spectrum is evident. For comparison, we present three spectra of two EL6 meteorites, Hvittis and Pillistfer, that best-matched with the Athor family. (B) 1.8/0.8 µm reflectance ratio versus the geometric visible albedo of the Athor family and the individual meteorite spectra (CM2 and CV) that achieved a curve-matching with Athor in the first step of the analysis. Another curve-matched meteorite is a bright enstatite achondrite (aubrite) that lies out of the range of the plot (Table B.3). For ELs and EHs, we present the respective averages. The properties of the Athor family members and the ELs overlap, whereas the other meteorite classes are clearly distinct.</figcaption></figure>

orbital elements and the correlations between the orbital proper semi-major axis and the inverse diameter (V-shape criterion, see Bolin et al. 2017; Delbo et al. 2017) is only a rough expression of the true membership (Nesvorny et al. 2015). This means that asteroids unrelated to a family can happen to have values of proper elements within the orbital range of that family, and are therefore grouped together with the true members. This problem could be particularly significant for Athor and Zita families, which are extensively overlapping in the orbital elements, and their members have similar $p_V$-values and visible spectrophotometry (Delbo et al. 2019).

It is natural to attempt to assign Athor members with Zita-type spectra to the Zita family and vice versa. However, we carefully inspected (i) whether or not the Zita-type asteroids that were originally assigned to the Athor family (Delbo et al. 2019) are placed inside the V-shape of Athor family, and (ii) whether or not their orbital elements are overlapping with the Zita family as well, and vice versa. We also visually inspected the spectra of each asteroid in order to assess the possibility that they were spectrally misclassified, that is, members of Athor and Zita are required to be within the spectroscopic X-complex. We find that asteroids (10 309) and (21 475), previously classified within the X-complex on the basis of spectrophotometric data and consequently originally included in Zita family, are in reality S-types, and should therefore be excluded from the Zita family. Asteroid (68 750), originally included in the Athor family (Delbo et al. 2019), is likely a low-albedo C-complex object and should therefore be excluded from the Athor family. Asteroid (757) has been reported in the literature as an Xk-type (Clark et al. 2009), but visual inspection shows that this literature spectrum is intermediate between that of an Xk- and an Xc-type. Moreover, we calculated the near-infrared slope of the Clark et al. (2009) spectrum and find that it is within the Athor-type slope range. We also note that asteroid (757) cannot be assigned to Zita on the basis of its position in the $(a_p, 1/D)$ space because it is outside the V-shape of Zita family, hence we kept it in Athor family.

We find that 17 out of 21 asteroids observed from the core of the Athor family have an Athor-type spectrum and at 95% significance level; this corresponds to 81 ± 15% of the total Athor core

population that counts 110 members. Of these 21 observed asteroids, 3 have a Zita-type spectrum and indeed belong to the Zita family (Table B.2). In the halo of Athor, which is a part of the family that has a far lower member number density than the core, we find that four out of eight asteroids (50%) are Athor-types, the remaining members being Zita-type. The asteroids with a Zita-type spectrum but originally assigned to the Athor halo could indeed belong to the Zita family because they lay inside the V-shape and the proper element range of the latter. The number of objects that we observed is sufficient to also ensure the statistical significance of our findings.

Finally, we calculated the average spectrum for the members of the families (Appendix B.3, Fig. 1A). The average spectrum of the Athor family is in very good agreement with the Xc-type of the Bus-DeMeo taxonomy, while the average spectrum of the Zita family has a slope similar to the X/T-types but also shows the presence of the 0.9 µm band, a characteristic of the Xk-types, placing Zita in an intermediate position between Bus-DeMeo taxonomic types. Furthermore, we calculated the uncertainty-weighed $p_V$ of the Athor and Zita families using the inverse of the uncertainty squared as weights. We selected this method as it combines measurements of the same quantity ($p_V$) obtained with different instruments or phases of space missions, which therefore have different accuracy. For the Athor family, we find that $p_V = 0.19 \pm 0.05$, while for Zita we obtain $p_V = 0.13 \pm 0.03$. The Zita family appears darker than the Athor family, but their average $p_V$s are within $1\sigma$ of uncertainty.

## 4. Spectral match of the Athor family to the ELs

We followed well-established approaches (Appendix B.3), consisting in (i) curve-matching (Popescu et al. 2012) between the average reflectance spectrum of the Athor family and meteorite spectra from RELAB and PSF, and (ii) comparison of the absolute reflectances of meteorites at 0.56 µm and their reflectance ratio at 1.8/0.8 µm (Cloutis et al. 2010; Beck et al. 2021) against the respective quantities of Athor family members (Fig. 1B). We find that Athor family spectroscopically matches with ELs, EHs, some carbonaceous chondrites, and one aubrite (Table B.3). This
is because the reflectance spectral curve of ELs may not be unique amongst meteorites. However, only the ELs have albedo values within the range of the Athor family (Fig. 1), concluding that ELs are the only known meteorites that can be matched to the Athor family. The effects of space weathering are not sufficiently significant to change the spectral slope or albedo of the ELs (Vernazza et al. 2009), implying that the direct match between asteroid and meteorite spectra is satisfactory (Fig. 1).

There is additional evidence of the link between the Athor family and the EL meteorites. First of all, the absence of 3 µm features from the asteroid (161) Athor itself has been reported (Rivkin et al. 2000; Usui et al. 2018), which can be used as an indication that the Athor family progenitor was not hydrated. The uniqueness of the link is also supported by (i) well-established thermal models of EL thermo-chronometers, which show that the meteorites used (all EL6) originate from the same depth of a single parent body (Trieloff et al. 2022), and (ii) the very narrow distribution of the CRE ages of the ELs, which means that they share a common travel history from their source to Earth (Trieloff et al. 2022). The only EL that was observed while entering Earth's atmosphere is the Neuschwanstein. This meteorite, before reaching the ground, had semi-major axis and inclination values (Spurný et al. 2003) inside the range of the Athor family (Delbo et al. 2019), demonstrating that the latter can deliver ELs to Earth. Moreover, EL pieces are found in the Almahata Sitta multi-lithology meteorites, whose parent body was asteroid 2008 TC3 (Zolensky et al. 2010). This body has 82% probability of being delivered to near-Earth space from the inner main belt, where the Athor family is also located. Therefore, it is reasonable to hypothesise that EL pieces coming from the Athor family were implanted via collisions on the 2008 TC3 parent body. We conclude that Athor is the only asteroidal source of the ELs.

## 5. The size of the Athor family progenitor

What was the size of the Athor family progenitor, before it was collisionally disrupted forming a family? An asteroid family loses members over time due to their collisional disruption and by orbital semi-major axis drift caused by non-gravitational forces that push them into zones of orbital instability (Nesvorny et al. 2015), as detailed earlier. Starting from different possible initial size distributions of the Athor family, we modelled their dynamical and collisional evolution until we found the one best matching with the observed family members (Appendix C, Fig. C.1). By integrating all of the volume from this initial size distribution, we deduce that the Athor family progenitor had a diameter of 64 km before its catastrophic disruption (Appendix C.1), leading to the formation of the Athor family, which is best fit at $3.0_{-0.4}^{+0.5}$ Gyr ago (Delbo et al. 2019). Given the method of Athor family identification (Delbo et al. 2017, 2019), it is likely that those members that contribute most of the volume have been identified.

## 6. Discussion

We show that the ELs originate from a family of asteroid fragments (EL secondary parents), where asteroid (161) Athor is the largest fragment. The link that we find between the Athor family and the ELs is consistent with other works (Vernazza et al. 2009; DeMeo et al. 2022), which indicated that ELs have spectra similar to asteroids classified as Xc-types; indeed, we classified the reflectance spectra of (161) Athor and the Athor family as Xc-type.

How unique is this link? We demonstrated that given the Athor family average reflectance spectrum, only ELs, amongst all meteorites, produce a match with the reflectance spectral curve and the albedo values, and that Athor is the only known collisional family with these characteristics. However, the uniqueness of the EL reflectance spectrum amongst main belt asteroids that are outside the Athor family is another open question. In the search for an answer, we proceeded as follows: First, we calculated the average reflectance spectrum of all the ELs present in the RELAB and PSF databases after their normalisation to 1 at 0.55 µm. We find that the EL reflectance spectra have very small dispersion with a relative standard deviation of <2.5% between 0.7 and 2.5 µm, implying that EL spectra are very similar to each other. Using the method described in Appendix B.2, we then found that only the Xc and the Xe Bus DeMeo spectral types provide a spectral curve match to the normalised reflectance spectra of ELs. However, the absolute reflectance value ($p_V$) of Xe is out of the range of ELs. We therefore deduce that only Xc asteroids with $p_V$ between 0.11 and 0.23 provide a good match to the ELs.

In addition to the Athor family members, we also pointed nine asteroids that appear to be Xc-types with $p_V$ values in the EL range. These are (21) Lutetia, (97) Klotho, (129) Antigone, (224) Oceana, (242) Kriemhild, (1541) Estonia, (2349) Kurchenko, (2567) Elba, and (5632) Ingelehmann (Fig. A.1). Apart from (21), (97), (129), and (224), these asteroids do not have a known near-infrared spectrum to confirm their Xc-type. Of those, only asteroids (21) and (97) were found to spectroscopically match the ELs (Vernazza et al. 2009), but only (21) resides in the inner main belt, which is a condition that dynamically favours the delivery of meteorites. All the rest are located at semi-major axes far from meteorite feeding zones (e.g. v6 and 3:1 MMR). Moreover, none of the asteroids above are associated with a collisional family; we survey all potential X-complex families in section 2, and find it to be unlikely that any of those are the sources of the ELs that we receive on Earth, although they could indeed have an EL composition.

We calculated that the size of the Athor family progenitor (the parent body of the family), finding a value of 64 km in diameter. However, the most likely diameter of the EL original planetesimal is 240–420 km, as is estimated by best-fitting thermal evolution models to EL thermo-chronometric data (Trieloff et al. 2022). A planetesimal of the size of the Athor family progenitor would have been too small and would therefore have had an excessively fast cooling rate to be consistent with the EL thermo-chronometers, which instead indicate a slow cooling rate (Trieloff et al. 2022).

The only possibility to explain the size difference between the Athor family progenitor and the original EL planetesimal is that the Athor family progenitor was a 64 km-diameter collisional fragment of the much larger EL planetesimal (Fig. 2). If the latter had suffered the break up collision at the current location of Athor, a larger collisional family would have been detected by methods that can identify families as old as our Solar System (Bolin et al. 2017; Delbo et al. 2017, 2019; Deienno et al. 2021). However, this large family does not exist, indicating that the breakup of the EL planetesimal did not occur in the main belt. We therefore deduce that the breakup of the original EL planetesimal happened elsewhere, most likely in the terrestrial planetary region. This is because a planetesimal with EL composition most likely accreted in the terrestrial region of our planetary disc (Kallemeyn & Wasson 1986; Zhang et al. 1995), as suggested from its striking isotopic similarity with the Earth and the similar low-oxygen fugacity
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Collisional history of the EL original planetesimal. The EL planetesimal accreted in the terrestrial planet region where it suffered a collisional event. A fragment of it—the Athor family progenitor—was transported to the current location of Athor in the inner main belt by a dynamical process. About 900 Myr later, the Athor family progenitor suffered an additional collisional event that led to the formation of the present Athor family. Meteorites are now transported from the Athor family to Earth via well-established routes.</figcaption></figure>

(Javoy et al. 2010; Dauphas 2017). Differences in chemical composition between ELs and the Earth can be explained via a fractionated condensation sequence starting from the same parent material (Morbidelli et al. 2020); hence, this should not be interpreted as an argument against the Earth-EL relationship.

Additional evidence that the current source of the ELs is a piece of a larger planetesimal comes from the study of the EL metamorphic grades. According to the thermal-evolution models (Trieloff et al. 2022), the EL planetesimal should have formed a large core of EL material of a high metamorphic grade 6 or 7, which was surrounded by thinner layers of lower metamorphic grades towards the surface (from EL5 to EL3). Although the known ELs sample all metamorphic grades, the EL3, EL4, and EL5 are small fractions of the total collection, which instead is dominated by the EL6. This indicates that the ELs come from a parent body (Ather family progenitor) of mainly EL6 composition, and therefore sample only a fraction of the original EL planetesimal, mostly the interior. This is supported by the fact that during the curve matching with the literature meteoritic spectra, the Athor family average spectrum matches better with EL6 (Table B.3). However, this EL Athor family progenitor should have also re-accumulated lower grade EL material from the upper layers of the EL original planetesimal during the initial collision in the terrestrial region, explaining the small presence of EL3-5 in the collections.

How could the EL Athor family progenitor be transported from the terrestrial planet region to the Athor family current location in the main belt? Several dynamical mechanisms have been proposed that could have led to the implantation of planetesimals from the terrestrial region to the main belt at the very early stages following the formation of the Solar System (Bottke et al. 2006; Raymond & Izidoro 2017). Further analysis of the properties of the EL original planetesimal may provide constraints on the time of its initial breakup and subsequent implantation into the main belt. This will clarify whether the existing implantation mechanisms are valid or a new dynamical mechanism is needed. The fact that original planetesimals underwent catastrophic collisions out of the belt and that their pieces were later implanted into the main belt means that some

of the leftover planetesimals from the main belt accreted larger than what is observed. This will have implications for the original size frequency distribution of the planetesimals (Delbo et al. 2017, 2019), which could be conservative at present.

## 7. Conclusions

We report the discovery of a unique source of the rare EL enstatite chondrite meteorites, namely the Athor asteroid family. The EL original planetesimal accreted in the terrestrial region and a comparison of its size to that of the Athor family progenitor suggests that the former broke via an impact, and a fragment of it (the progenitor of Athor) was transported into the inner main belt. The latter suffered a new catastrophic collision, forming the currently observed Athor family.

**Acknowledgements.** We acknowledge support from the ANR ORIGINS (ANR-18-CE31-0014). CA was supported by the French National Research Agency under the project "Investissements d'Avenir" UCA<sup>EDI</sup> (ANR-15-IDEX-01) and by a European Space Agency Research Fellowship. AM acknowledges support from the ERC advanced grant HolyEarth N. 101019380. KJW acknowledges support through Project ESPRESSO, a NASA SSERVI program at SwRI. CA and MD were Visiting Astronomers at the Infrared Telescope Facility, which is operated by the University of Hawaii under contract 80HQTR19D0030 with the National Aeronautics and Space Administration. This work is partially based on observations collected at Copernico telescope (Asiago, Italy) of the INAF – Osservatorio Astronomico di Padova. These results made use of the Lowell Discovery Telescope at Lowell Observatory. Lowell is a private, non-profit institution dedicated to astrophysical research and public appreciation of astronomy and operates the LDT in partnership with Boston University, the University of Maryland, the University of Toledo, Northern Arizona University and Yale University. The upgrade of the DeVeny optical spectrograph has been funded by a generous grant from John and Ginger Giovale and by a grant from the Mt. Cuba Astronomical Foundation. This work is based on data provided by the Minor Planet Physical Properties Catalogue (MP3C) of the Observatoire de la Côte d'Azur (mp3c.oca.eu). Meteorite spectra acquired from the NASA RELAB facility (Brown University) and Planetary Spectrophotometer Facility (University of Winnipeg) databases. We thank S.J. Bus, A.C.A. Boogert and the IRTF telescope operators for their help during the observations. We thank Thomas Burbine for his comments, allowing us to improve the manuscript.

## References

Alf-Lagoa, V., Müller, T. G., Usui, F., & Hasegawa, S. 2018, A&A, 612, A85

Beck, P., Schmitt, B., Potin, S., Pommerol, A., & Brissaud, O. 2021, Icarus, 354, 114066

Bida, T. A., Dunham, E. W., Massey, P., & Roe, H. G. 2014, in Ground-based and Airborne Instrumentation for Astronomy V, eds. S. K. Ramsay, I. S. McLean, & H. Takami, SPIE Conf. Ser., 9147, 91472N

Bolin, B. T., Delbo, M., Morbidelli, A., & Walsh, K. J. 2017, Icarus, 282, 290

Bottke, W. F., Durda, D. D., Nesvorný, D., et al. 2005, Icarus, 179, 63

Bottke, W. F., Nesvorný, D., Grimm, R. E., Morbidelli, A., & O'Brien, D. P. 2006, Nature, 439, 821

Broz, M., Ferrais, M., Vernazza, P., Sevecek, P., & Jutzi, M. 2022, A&A, 664, A69

Brunetto, R., Loeffler, M. J., Nesvorný, D., Sasaki, S., & Strazzulla, G. 2015, in Asteroids IV, eds. P. Michel, F. E. DeMeo, W. F. Bottke, et al., 597

Bus, S. J., & Binzel, R. P. 2002, Icarus, 158, 106

Carvano, J. M., Hasselmann, P. H., Lazzaro, D., & Mothé-Diniz, T. 2010, A&A, 510, A43

Clark, B. E., Ockert-Bell, M. E., Cloutis, E. A., et al. 2009, Icarus, 202, 119

Clayton, R. N., & Mayeda, T. K. 1996, Geochim. Cosmochim. Acta., 60, 1999

Cloutis, E. A., Gaffey, M. J., Smith, D. G. W., & Lambert, R. S. J. 1990, J. Geophys. Res., 95, 8323

Cloutis, E. A., Hardersen, P. S., Bish, D. L., et al. 2010, MAPS, 45, 304

Cushing, M. C., Vacca, W. D., & Rayner, J. T. 2004, PASP, 116, 362

Dauphas, N. 2017, Nature, 541, 521

Deienno, R., Walsh, K. J., & Delbo, M. 2021, Icarus, 357, 114218

Delbo, M., Walsh, K., Bolin, B., Avdellidou, C., & Morbidelli, A. 2017, Science, 357, 1026

Delbo, M., Avdellidou, C., & Morbidelli, A. 2019, A&A, 624, A69

DeMeo, F. E., & Carry, B. 2013, Icarus, 226, 723

DeMeo, F. E., Binzel, R. P., Slivan, S. M., & Bus, S. J. 2009, Icarus, 202, 160
DeMeo, F. E., Alexander, C. M. O., Walsh, K. J., Chapman, C. R., Binzel, R. P., et al. 2015, in *Asteroids IV*, eds. P. Michel, F. E. DeMeo, W. F. Bottke, et al.

DeMeo, F. E., Burt, B. J., Marsset, M., et al. 2022, *Icarus*, **380**, 114218

Devogèle, M., Moskovitz, N., Thirouin, A., et al. 2019, *AJ*, **158**, 196

Dohnanyi, J. S. 1969, *J. Geophys. Res.*, **74**, 2531

Fornasier, S., Lazzarin, M., Barbieri, C., & Barucci, M. A. 1999, *A&AS*, **135**, 65

Fornasier, S., Clark, B. E., Dotto, E., et al. 2010, *Icarus*, **210**, 655

Gaffey, M. J., Reed, K. L., & Kelley, M. S. 1992, *Icarus*, **100**, 95

Gayon-Markt, J., Delbo, M., Morbidelli, A., & Marchi, S. 2012, *MNRAS*, **424**, 508

Greenwood, R. C., Burbine, T. H., & Franchi, I. A. 2020, *Geochim. Cosmochim. Acta.*, **277**, 377

Gulbis, A. A. S., Bus, S. J., Elliot, J. L., et al. 2011, *PASP*, **123**, 461

Henke, S., Gail, H. P., Trieloff, M., Schwarz, W. H., & Kleine, T. 2012, *A&A*, **545**, A135

Herwartz, D., Pack, A., Friedrichs, B., & Bischoff, A. 2014, *Science*, **344**, 1146

Ivezić, Ž., Tabachnik, S., Rafikov, R., et al. 2001, *AJ*, **122**, 2749

Javoy, M., Kaminski, E., Guyot, F., et al. 2010, *Earth Planet. Sci. Lett.*, **293**, 259

Jenniskens, P., Shaddad, M. H., Numan, D., et al. 2009, *Nature*, **458**, 485

Johansen, A., Jacquet, E., Cuzzi, J. N., Morbidelli, A., Gounelle, M., et al. 2015, in *Asteroids IV*, eds. P. Michel, F. E. DeMeo, W. F. Bottke, et al., 471

Kallemeyn, G. W., & Wasson, J. T. 1986, *Geochim. Cosmochim. Acta.*, **50**, 2153

Keil, K. 1968, *J. Geophys. Res.*, **73**, 6945

Keil, K. 2010, *Geochemistry*, **70**, 295

Lazzaro, D., Angeli, C. A., Carvano, J. M., et al. 2004, *Icarus*, **172**, 179

Levison, H. F., & Duncan, M. J. 1994, *Icarus*, **108**, 18

Lin, Y., & El Goresy, A. 2002, *MAPS*, **37**, 577

Lucas, M. P., Emery, J. P., MacLennan, E. M., et al. 2019, *Icarus*, **322**, 227

Marsset, M., DeMeo, F. E., Binzel, R. P., et al. 2020, *ApJS*, **247**, 73

Morbidelli, A., Libourel, G., Palme, H., Jacobson, S. A., & Rubie, D. C. 2020, *Earth Planet. Sci. Lett.*, **538**

Nesvorny, D., Broz, M., Carruba, V., et al. 2015, in *Asteroids IV*, eds. P. Michel, F. E. DeMeo, W. F. Bottke, et al., 297

Norton, O. R., & Chitwood, L. A. 2008, *Field Guide to Meteors and Meteorites*

Popescu, M., Birlan, M., & Nedelcu, D. A. 2012, *A&A*, **544**, A130

Popescu, M., Licandro, J., Morate, D., et al. 2016, *A&A*, **591**, A115

Raymond, S. N., & Izidoro, A. 2017, *Sci. Adv.*, **3**, e1701138

Rayner, J. T., Toomey, D. W., Onaka, P. M., et al. 2003, *PASP*, **115**, 362

Reddy, V., Emery, J. P., Gaffey, M. J., et al. 2009, *MAPS*, **44**, 1917

Reddy, V., Dunn, T. L., Thomas, C. A., Moskovitz, N. A., Burbine, T. H., et al. 2015, in *Asteroids IV*, eds. P. Michel, F. E. DeMeo, W. F. Bottke, et al., 43

Rivkin, A. S., Howell, E. S., Lebofsky, L. A., Clark, B. E., & Britt, D. T. 2000, *Icarus*, **145**, 351

Rozehnal, J., Brož, M., Nesvorný, D., et al. 2022, *Icarus*, **383**, 115064

Russell, C. T., Raymond, C. A., Coradini, A., et al. 2012, *Science*, **336**, 684

Spurný, P., Oberst, J., & Heinlein, D. 2003, *Nature*, **423**, 151

Tanga, P., Cellino, A., Michel, P., et al. 1999, *Icarus*, **141**, 65

Trieloff, M., Jessberger, E. K., Herrwerth, I., et al. 2003, *Nature*, **422**, 502

Trieloff, M., Hopp, J., & Gail, H.-P. 2022, *Icarus*, **373**, 114762

Usui, F., Hasegawa, S., Ootsubo, T., & Onaka, T. 2018, *PASJ*, **71**, 1

Vernazza, P., Brunetto, R., Binzel, R. P., et al. 2009, *Icarus*, **202**, 477

Vokrouhlický, D., Bottke, W. F., Chesley, S. R., Scheeres, D. J., Statler, T. S., et al. 2015, in *Asteroids IV*, eds. P. Michel, F. E. DeMeo, W. F. Bottke, et al., 509

Xu, S., Binzel, R. P., Burbine, T. H., & Bus, S. J. 1995, *Icarus*, **115**, 1

Zhang, Y., Benoit, P. H., & Sears, D. W. G. 1995, *J. Geophys. Res.*, **100**, 9417

Zolensky, M., Herrin, J., Mikouchi, T., et al. 2010, *MAPS*, **45**, 1618
# Appendix A: Table of X-complex asteroid families and location of individual Xc-type asteroids

<table><caption>**Table A.1.** Average $p_V$-values of potential X-complex families.</caption><thead><tr><th>Asteroid Family</th><th>Location</th><th>$p_V$</th></tr></thead><tbody><tr><td>87 Sylvia</td><td>OMB</td><td>0.078</td></tr><tr><td>161 Athor</td><td>IMB</td><td>0.184</td></tr><tr><td>293 Brasilia</td><td>OMB</td><td>0.081</td></tr><tr><td>322 Phaeo</td><td>CMB</td><td>0.077</td></tr><tr><td>363 Padua</td><td>CMB</td><td>0.076</td></tr><tr><td>396 Aeolia</td><td>CMB</td><td>0.083</td></tr><tr><td>369 Aeria</td><td>CMB</td><td>0.081</td></tr><tr><td>589 Croatia</td><td>OMB</td><td>0.073</td></tr><tr><td>689 Zita</td><td>IMB</td><td>0.177</td></tr><tr><td>709 Fringilla</td><td>OMB</td><td>0.075</td></tr><tr><td>909 Ulla</td><td>OMB</td><td>0.057</td></tr><tr><td>1222 Tina</td><td>CMB</td><td>0.076</td></tr><tr><td>1303 Luthera</td><td>OMB</td><td>0.072</td></tr><tr><td>5567 Durisen</td><td>OMB</td><td>0.076</td></tr><tr><td>7481 San Marcello</td><td>OMB</td><td>0.177</td></tr></tbody></table>

**Notes.** Reported average $p_V$-values for the Athor and Zita families are the ones that were retrieved using the original family membership by Delbo et al. (2019) and before a re-assessment of the family membership, which was performed after our observations. Column location refers to the region of the main belt where the family is located following Nesvorny et al. (2015), namely IMB = Inner Main Belt, $2.0 < a < 2.5$ au, $i < 17.5^{\circ}$; CMB = Central Main Belt, $2.5 < a < 2.82$ au, $i < 17.5^{\circ}$; OMB = Outer Main Belt, $2.82 < a < 3.7$ au, $i < 17^{\circ}$, where *a* and *i* are the orbital semi-major axis and inclination.

<figure><img src="image_4.png" /><figcaption>**Fig. A.1.** The location in the main belt of nine potential Xc-type asteroids with $p_V$ values within the range of the ELs. No one of these belong to an asteroid family.</figcaption></figure>

# Appendix B: Spectroscopic Analysis

## B.1. Literature data and observations

Considering the Athor and Zita family membership, we initially searched for literature spectral and spectro-photometric data. Spectra exist in the near-infrared (NIR) only for (161) Athor (Bus & Binzel 2002; Fornasier et al. 2010) and (757) Portlandia

(Clark et al. 2009) and in the visible (VIS) for other 22 asteroids (Table B.1). The remaining family members have at least spectrophotometric data from the Sloan Digital Sky Survey (SDSS) (Ivezić et al. 2001) and/or near-infrared colours from VISTA-VHS survey (MOVIS) (Popescu et al. 2016).

We performed new compositionally diagnostic spectroscopic observations (Table B.1) mostly in the NIR, but also in the VIS of Athor and Zita family members. We used the SpeX (Rayner et al. 2003) and MORIS instruments (Gulbis et al. 2011) at the NASA Infrared Telescope Facility (IRTF), the Near-Infrared High-Throughput Spectrograph (NIHTS; Bida et al. 2014) and the DeVeny optical spectrograph (Bida et al. 2014) at the Lowell Discovery Telescope (LDT), as well as the Asiago Faint Object Spectrograph and Camera (AFOSC) at the Coper-nico Telescope (Table B.1). SpeX was used with a slit 0.8×15" in PRISM mode, covering wavelengths between 0.7 and 2.5 µm with a spectral resolution ~200 in a single configuration. NIHTS covers the 0.84 to 2.4 µm range and was used with the 1.34"-wide slit. For this slit the spectral resolution is 2194 pix⁻¹ at 1 µm, 667 pix⁻¹ at 1.5 µm and 535 pix⁻¹ at 2 µm. For NIR observations we used the commonly adopted ABBA observation procedure which consists in alternating our targets between two positions along the slit (A and B), therefore allowing us to remove the background contribution. DeVeny was used with the R150 grating providing a dispersion of 0.43 nm pix⁻¹ covering a spectral range between 0.32 to 1 µm. AFOSC was used with the VPH#6 grism covering the 0.4–1.0 µm range and the 4.22"-wide slit resulting in a spectral resolution of 0.295 nm pix⁻¹.

In the NIR, following well-established procedures (Reddy et al. 2009), we calculated the asteroid's reflectance $R(\lambda)$ as a function of the wavelength $\lambda$ using Eq. B.1:

$$R(\lambda) = \frac{A(\lambda)}{S_L(\lambda)} \times \text{Poly}\left(\frac{S_L(\lambda)}{S_T(\lambda)}\right), \qquad (B.1)$$

where $A(\lambda)$, $S_L(\lambda)$, $S_T(\lambda)$ are the wavelength-calibrated raw spectra respectively of the asteroid, a local G2-type star observed within ~300" of the asteroid, and a well-studied (trusted) solar analogue star that was observed at similar airmass when possible respectively (table S2). The function *Poly*() represents a polynomial fit of the stars ratio, excluding the regions affected by the telluric water vapour absorption (1.3 < $\lambda$ < 1.5, 1.78 < $\lambda$ < 2.1, and $\lambda$ > 2.4 µm). Asteroid spectra were shifted to sub-pixel accuracy to align with the calibration star spectra. In general, the local star ensures accurate removal of the telluric features, but may require the slope correction $Poly(\frac{S_L(\lambda)}{S_T(\lambda)})$ due to difference between the spectrum of the local star and that of the Sun. When the local solar analogue was not observed, the correction of the telluric features was performed directly by dividing the asteroid spectrum by that of the trusted solar analogue, i.e: $R(\lambda) = \frac{A(\lambda)}{S_T(\lambda)}$. This latter procedure is also used for visible-light spectroscopy. In the NIR we collected data with the slit oriented to the parallactic angle so that any atmospheric dispersion is along the slit, which minimises the light loss. In the VIS observations from Asiago, the rotator was not used because the 4.22" slit was large enough to compensate for differential refraction.

SpeX observations were reduced using the interactive data language (IDL)-based spectral reduction tool Spextool (v4.1 Cushing et al. 2004), while NITHS observations with a version of the Spextool adapted to this instrument. DeVeny data were reduced using standard reduction techniques for spectroscopy as detailed in Devogèle et al. (2019). AFOSC data were analysed using the ESO-MIDAS with standard reduction methods as detailed in Fornasier et al. (1999). In all cases we performed flat
<figure><img src="image_5.png" /><figcaption>**Fig. B.1.** NIR spectra (black points) were obtained from our survey and are presented after we applied a $2\sigma$ clipping in order to eliminate the scattered datapoints inside the telluric bands (grey areas) by as much as possible. We also present their existing VIS counterpart (grey points) obtained also from our survey or retrieved from the literature, as well as the SDSS data (red points).</figcaption></figure>

fielding and wavelength calibration, while for the VIS we also
performed bias and atmospheric extinction corrections.

We also included in our analysis spectrophotometry from the
SDSS MOC4 (Ivezić et al. 2001), which contains magnitudes
at five spectral bands (u', g', r', i', z') centred at 0.354, 0.477,

0.623, 0.763, and 0.913 µm, respectively Carvano et al. (2010),
Gayon-Markt et al. (2012). We transformed SDSS magnitudes
to reflectance values following the procedure of DeMeo & Carry
(2013) and we did not use the u' magnitudes because of their
very large errors.
<figure><img src="image_6.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>

## B.2. Classification of spectra

First, we normalised the NIR reflectances of asteroids at 1 µm and calculated their slopes between 0.8 and 2.4 µm avoiding the telluric absorption regions, as it is commonly done (Marsset et al. 2020; DeMeo et al. 2009). Next, we classified each asteroid reflectance in the Bus-DeMeo taxonomy (DeMeo et al. 2009). To do so, we performed a reflectance curve-matching by minimising value of the $\chi^2$-figure of merit (Eq. B.2):

$$ \chi_M^2 = \sum_{\lambda} \frac{[O(\lambda) - \beta M(\lambda)]^2}{\sigma(\lambda)^2} \quad \text{with} \quad \beta = \frac{\sum_{\lambda} O(\lambda) M(\lambda)}{\sum_{\lambda} M(\lambda)^2}, \qquad (B.2) $$

where $O(\lambda)$ is the observed asteroid reflectance (i.e. $O(\lambda) = R(\lambda)$), $\sigma$ its uncertainty, and $M(\lambda)$ the reflectance model of each Bus-DeMeo class DeMeo et al. (2009), downloaded from the Small Main-Belt Asteroid Spectroscopic Survey<fnref n="1" />. The parameter $\beta$ is required to take into account that reflectances are arbitrarily normalised. For each asteroid, we calculated independently the $\chi_M^2$ for the NIR and visible spectra as well as the SDSS spectrophotometry where applicable, and sum the results. We present the two best-matching classes in Table B.3, that is, those with the

lowest and second-to-lowest $\chi_M^2$. The definition of the best spectral class is not based only on the output of the best fits, but also on the object's $p_V$ together with the visual inspection of the spectrum. Classes such as Xe, Xk, and Ch/Cgh cannot be proposed as the best ones because the absorption bands at 0.49 µm, 0.9 µm, and 0.7 µm, respectively, are missing.

## B.3. Meteorite spectral matching

For each asteroid, we combined NIR, VIS, and SDSS reflectances by normalising them at 0.85 µm. The normalisation was performed by fitting a natural spline (Python csaps module) with a smoothing coefficient value of 0.9 to each reflectance independently, the values of which were divided by the values of the spline at 0.85 µm. Next, we calculated the median and the absolute median deviation (MAD) of the reflectances of the Athor members in 0.0125 µm-wide wavelength bins between 0.5 and 2.5 µm. We removed in each bin those data-points which distance was > 2xMAD from the median; next we calculated the mean and standard deviation in each wavelength bin. The result of this computation represents the average spectra of the Athor and Zita families (Fig. 1 A).

In order to match the Athor family average spectrum we used the databases of the meteorite reflectance spectra from RELAB (Brown University) and PSF (University of Winnipeg).

<footnote n="1">¹ http://smass.mit.edu/_documents/busdemeo-meanspectra.xlsx</footnote>
<table><caption>**Table B.1.** VIS and NIR observations for asteroids and calibrating stars are reported along with the telescope and instrument that was used each time. The phase angle at the time of the observations is also reported.</caption><thead><tr><th>Asteroid Number</th><th>Name</th><th>Date (UT)</th><th>Phase angle(°)</th><th>Telescope</th><th>Instrument</th><th>SA trusted</th><th>SA local</th></tr></thead><tbody><tr><td>161</td><td>Athor</td><td>03-10-2017</td><td>23.5</td><td>IRTF</td><td>SpeX</td><td>SA 93-101</td><td>-</td></tr><tr><td>689</td><td>Zita</td><td>23-12-2019</td><td>17.4</td><td>Copernico</td><td>AFOSC</td><td>Hyades 64</td><td>-</td></tr><tr><td>689</td><td>Zita</td><td>04-11-2019</td><td>23.2</td><td>IRTF</td><td>SpeX</td><td>Hyades 64</td><td>HD 76446</td></tr><tr><td>757</td><td>Portlandia</td><td>01-10-2017</td><td>23.7</td><td>IRTF</td><td>SpeX</td><td>SA 112-1333</td><td>-</td></tr><tr><td>1063 1063 1490</td><td>Limpopo</td><td>13-12-2020</td><td>4.8</td><td>Copernico</td><td>AFOSC</td><td>Hyades 64</td><td>-</td></tr><tr><td>1998</td><td>Titius</td><td>23-12-2019</td><td>11.7</td><td>Copernico</td><td>AFOSC</td><td>Hyades 64</td><td>-</td></tr><tr><td>1998</td><td>Titius</td><td>04-11-2019</td><td>12.0</td><td>IRTF</td><td>SpeX</td><td>Hyades 64</td><td>HD 280114</td></tr><tr><td>2015</td><td>Kachuevskaya</td><td>23-12-2019</td><td>15.2</td><td>Copernico</td><td>AFOSC</td><td>Hyades 64</td><td>-</td></tr><tr><td>2015</td><td>Kachuevskaya</td><td>27-09-2019</td><td>20.7</td><td>IRTF</td><td>SpeX</td><td>Hyades 64</td><td>HD 275951</td></tr><tr><td>2419</td><td>Moldavia</td><td>27-09-2019</td><td>19.4</td><td>IRTF</td><td>SpeX</td><td>SA 93-101</td><td>HD 198259</td></tr><tr><td>3007</td><td>Reaves</td><td>27-09-2019</td><td>15.2</td><td>IRTF</td><td>SpeX</td><td>Hyades 64</td><td>HD 16640</td></tr><tr><td>3427</td><td>Szentmartoni</td><td>04-11-2019</td><td>24.1</td><td>IRTF</td><td>SpeX</td><td>SA 115-271</td><td>-</td></tr><tr><td>3665</td><td>Fitzgerald</td><td>22-01-2020</td><td>6.3</td><td>IRTF</td><td>SpeX</td><td>SA 102-1081</td><td>HD 77730</td></tr><tr><td>3704</td><td>Gaoshiqi</td><td>03-01-2020</td><td>15.5</td><td>IRTF</td><td>SpeX</td><td>SA 102-1081</td><td>HD 82606</td></tr><tr><td>4256</td><td>Kagamigawa</td><td>23-12-2019</td><td>2.2</td><td>Copernico</td><td>AFOSC</td><td>Hyades 64</td><td>-</td></tr><tr><td>4256</td><td>Kagamigawa</td><td>03-01-2020</td><td>3.5</td><td>IRTF</td><td>SpeX</td><td>SA 98-978</td><td>TYC 1881-556-1</td></tr><tr><td>4353</td><td>Onizaki</td><td>01-10-2017</td><td>17.0</td><td>IRTF</td><td>SpeX</td><td>SA 110-361</td><td>-</td></tr><tr><td>4548</td><td>Wielen</td><td>04-11-2019</td><td>2.0</td><td>IRTF</td><td>SpeX</td><td>Hyades 64</td><td>TYC 1230-630-1</td></tr><tr><td>5171</td><td>Augustesen</td><td>03-10-2017</td><td>25.1</td><td>IRTF</td><td>SpeX</td><td>SA 113-276</td><td>-</td></tr><tr><td>5236</td><td>Yoko</td><td>23-01-2020</td><td>6.0</td><td>IRTF</td><td>SpeX</td><td>SA 98-978</td><td>HD 60696</td></tr><tr><td>5343</td><td>Ryzhov</td><td>22-01-2020</td><td>20.9</td><td>IRTF</td><td>SpeX</td><td>SA 102-1081</td><td>TYC 872-431-1</td></tr><tr><td>6123</td><td>Aristoteles</td><td>04-11-2019</td><td>19.0</td><td>IRTF</td><td>SpeX</td><td>SA 115-271</td><td>HD 224651</td></tr><tr><td>6245</td><td>Ikufumi</td><td>21-08-2021</td><td>21.7</td><td>IRTF</td><td>SpeX</td><td>SA 110361</td><td>HD 172404</td></tr><tr><td>6247</td><td>Amanogawa</td><td>27-09-2019</td><td>14.6</td><td>IRTF</td><td>SpeX</td><td>SA 115-271</td><td>HD 212231</td></tr><tr><td>6855</td><td>Armellini</td><td>23-01-2020</td><td>6.9</td><td>IRTF</td><td>SpeX</td><td>SA 98-978</td><td>HD 56513</td></tr><tr><td>7111</td><td>1985QA1</td><td>27-09-2019</td><td>5.0</td><td>IRTF</td><td>SpeX</td><td>SA 93-101</td><td>HD 1267</td></tr><tr><td>7116</td><td>Mentall</td><td>03-08-2021</td><td>24.8</td><td>LDT</td><td>DeVeny</td><td>SA 93-101</td><td>-</td></tr><tr><td>7116</td><td>Mentall</td><td>01-10-2017</td><td>12.5</td><td>IRTF</td><td>SpeX</td><td>SA 115-271</td><td>-</td></tr><tr><td>10309</td><td>1990QC6</td><td>03-01-2020</td><td>7.3</td><td>IRTF</td><td>SpeX</td><td>SA 98-978</td><td>U Gem star 9</td></tr><tr><td>13101</td><td>Fransson</td><td>22-08-2021</td><td>5.9</td><td>IRTF</td><td>SpeX</td><td>SA 115-271</td><td>HD 201529</td></tr><tr><td>14950</td><td>1996BE2</td><td>03-01-2020</td><td>26.1</td><td>IRTF</td><td>SpeX</td><td>SA 102-1081</td><td>HD 99392</td></tr><tr><td>15996</td><td>1998YC12</td><td>23-01-2020</td><td>4.9</td><td>IRTF</td><td>SpeX</td><td>SA 102-1081</td><td>HD 79078</td></tr><tr><td>16007</td><td>Kaasalainen</td><td>03-08-2021</td><td>4.8</td><td>LDT</td><td>DeVeny</td><td>SA 110-361</td><td>-</td></tr><tr><td>16007</td><td>Kaasalainen</td><td>06-02-2020</td><td>0.8</td><td>LDT</td><td>NITHS</td><td>Hyades 64</td><td>HD 79930</td></tr><tr><td>17635</td><td>1996OC1</td><td>04-11-2019</td><td>7.7</td><td>IRTF</td><td>SpeX</td><td>Hyades 64</td><td>TYC 72-921-1</td></tr><tr><td>19207</td><td>1992QS1</td><td>21-08-2021</td><td>8.6</td><td>IRTF</td><td>SpeX</td><td>SA 115-271</td><td>HD 217254</td></tr><tr><td>20566</td><td>Laurielee</td><td>24-08-2021</td><td>12.6</td><td>IRTF</td><td>SpeX</td><td>SA 110361</td><td>HD 196164</td></tr><tr><td>21475</td><td>Jasonclain</td><td>22-01-2020</td><td>4.0</td><td>IRTF</td><td>SpeX</td><td>SA 98-978</td><td>HD 73708</td></tr><tr><td>33110</td><td>1998AM10</td><td>27-09-2019</td><td>4.0</td><td>IRTF</td><td>SpeX</td><td>SA 93-101</td><td>HD 1267</td></tr><tr><td>34120</td><td>2000PL28</td><td>23-12-2019</td><td>6.1</td><td>Copernico</td><td>AFOSC</td><td>Hyades 64</td><td>-</td></tr><tr><td>34120</td><td>2000PL28</td><td>22-01-2020</td><td>11.9</td><td>IRTF</td><td>SpeX</td><td>SA 98-978</td><td>HD 257880</td></tr><tr><td>37983</td><td>1998HB136</td><td>27-09-2019</td><td>4.9</td><td>IRTF</td><td>SpeX</td><td>SA 93-101</td><td>HD 1267</td></tr><tr><td>38233</td><td>1999NS57</td><td>01-10-2017</td><td>2.7</td><td>IRTF</td><td>SpeX</td><td>SA 115-271</td><td>-</td></tr><tr><td>54355</td><td>2000KJ33</td><td>23-01-2020</td><td>11.3</td><td>IRTF</td><td>SpeX</td><td>SA 102-1081</td><td>HD 88371</td></tr><tr><td>68750</td><td>2002EM68</td><td>02-10-2017</td><td>3.8</td><td>IRTF</td><td>SpeX</td><td>SA 112-1333</td><td>-</td></tr><tr><td>74119</td><td>1998QH52</td><td>03-10-2017</td><td>16.8</td><td>IRTF</td><td>SpeX</td><td>SA 112-1333</td><td>-</td></tr></tbody></table>

For curve matching, we performed a Monte Carlo simulation, where at each iteration, the Athor family average reflectance $\bar{R}(\lambda)$ is varied within the standard deviation of the reflectance distribution of the family members $\bar{\sigma}(\lambda)$ and compared to all meteorite reflectances $R(\lambda)_M$. For this, we again use Eq. B.2 assigning $O = \bar{R}, M = \bar{R}_M$, and $\sigma = \bar{\sigma}$. The meteorite spectrum scoring the best comparison (minimum $\chi^2_M$) was recorded. After $10^5$ Monte Carlo iterations, the distribution of scoring is a good approximation of the probability that a meteorite spectrum

matches the average reflectance of Athor family members. There are 17 meteorite spectra that carry 99% of the cumulative probability to match, corresponding to nine individual meteorites of which five are enstatite chondrites (EL, EH groups), three are carbonaceous chondrites (CM2 and CV2 groups), and one is an enstatite achondrite (aubrite).

Curve-matching of the reflectances is a condition necessary but not always sufficient to establish an asteroid-meteorite link (especially in the case of spectra with weak or no absorption
<table><caption>**Table B.2.** Summary of the asteroid members of Athor and Zita families for which either we performed spectroscopic observations or retrieved literature spectra.</caption><thead><tr><th>Asteroid Number</th><th>Original Membership</th><th>NIR</th><th>VIS</th><th>SDSS</th><th>NIR Slope %µm<sup>-1</sup></th><th>Spectral Class</th><th>p<sub>V</sub></th><th>Revised Membership</th></tr></thead><tbody><tr><td>161</td><td>Athor core</td><td>This work</td><td>(a)</td><td>n/a</td><td>7.22 ± 0.1</td><td>Xc, Xn</td><td>0.23</td><td>Athor</td></tr><tr><td>689</td><td>Zita</td><td>This work</td><td>This work</td><td>n/a</td><td>25.35 ± 0.33</td><td>T, X</td><td>0.10</td><td>Zita</td></tr><tr><td>757</td><td>Athor core</td><td>This work</td><td>(a)</td><td>yes</td><td>6.19 ± 0.15</td><td>Xe, Xc</td><td>0.22</td><td>Athor</td></tr><tr><td>1490</td><td>Athor halo</td><td>n/a</td><td>This work</td><td>n/a</td><td>-</td><td>Xe, Xc</td><td>0.17</td><td>Athor</td></tr><tr><td>1998</td><td>Athor halo</td><td>This work</td><td>This work</td><td>yes</td><td>8.80 ± 0.14</td><td>Xc, Cgh</td><td>0.10</td><td>Athor</td></tr><tr><td>1998</td><td>Athor halo</td><td>This work</td><td>(a)</td><td>yes</td><td>8.80 ± 0.14</td><td>Cgh, Xc</td><td>0.10</td><td>Athor</td></tr><tr><td>2015</td><td>Athor halo</td><td>This work</td><td>This work</td><td>yes</td><td>22.84 ± 0.29</td><td>X, Xk</td><td>0.14</td><td>Zita</td></tr><tr><td>2419</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>8.26 ± 0.30</td><td>Xc, Cgh</td><td>0.24</td><td>Athor</td></tr><tr><td>3007</td><td>Athor halo</td><td>This work</td><td>(a)</td><td>n/a</td><td>17.91 ± 0.22</td><td>Xk, X</td><td>0.15</td><td>Zita</td></tr><tr><td>3427</td><td>Zita</td><td>This work</td><td>n/a</td><td>n/a</td><td>7.20 ± 0.18</td><td>Ch, Xc</td><td>0.22</td><td>Athor</td></tr><tr><td>3665</td><td>Athor halo</td><td>This work</td><td>(b)</td><td>yes</td><td>20.06 ± 0.21</td><td>X, Xk</td><td>0.20</td><td>Zita</td></tr><tr><td>3704</td><td>Athor halo</td><td>This work</td><td>(a)</td><td>n/a</td><td>10.82 ± 0.25</td><td>Xc, Xe</td><td>0.21</td><td>Athor</td></tr><tr><td>4256</td><td>Zita</td><td>This work</td><td>(a)</td><td>n/a</td><td>22.23 ± 0.47</td><td>Sq, Xk</td><td>0.23</td><td>Zita</td></tr><tr><td>4256</td><td>Zita</td><td>This work</td><td>This work</td><td>n/a</td><td>22.23 ± 0.47</td><td>Sq, X</td><td>0.23</td><td>Zita</td></tr><tr><td>4353</td><td>Athor core</td><td>This work</td><td>(a)</td><td>yes</td><td>4.33 ± 0.17</td><td>Xe, Xc</td><td>0.15</td><td>Athor</td></tr><tr><td>4353</td><td>Athor core</td><td>This work</td><td>(b)</td><td>yes</td><td>4.33 ± 0.17</td><td>Xe, Xc</td><td>0.15</td><td>Athor</td></tr><tr><td>4548</td><td>Athor halo</td><td>This work</td><td>(a)</td><td>yes</td><td>14.11 ± 0.24</td><td>Xk, Cb</td><td>0.18</td><td>Zita</td></tr><tr><td>5171</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>17.27 ± 0.38</td><td>T, Xk</td><td>0.11</td><td>Zita</td></tr><tr><td>5236</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>18.43 ± 0.23</td><td>T, X</td><td>0.15</td><td>Zita</td></tr><tr><td>5343</td><td>Athor core</td><td>This work</td><td>(c)</td><td>yes</td><td>9.91 ± 0.39</td><td>Xc, Cgh</td><td>0.29</td><td>Athor</td></tr><tr><td>6123</td><td>Athor halo</td><td>This work</td><td>n/a</td><td>yes</td><td>9.52 ± 0.42</td><td>Xe, Xc</td><td>0.16</td><td>Athor</td></tr><tr><td>6245</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>6.19 ± 0.21</td><td>Xe, Xc</td><td>0.13</td><td>Athor</td></tr><tr><td>6247</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>9.21 ± 0.22</td><td>Xc, Xe</td><td>0.17</td><td>Athor</td></tr><tr><td>6855</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>9.91 ± 0.39</td><td>Xe, Xc</td><td>0.13</td><td>Athor</td></tr><tr><td>7111</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>7.29 ± 0.33</td><td>Xc, Xe</td><td>0.22</td><td>Athor</td></tr><tr><td>7116</td><td>Athor core</td><td>This work</td><td>This work</td><td>yes</td><td>3.63 ± 0.23</td><td>Xe, Xc</td><td>0.15</td><td>Athor</td></tr><tr><td>10309</td><td>Zita</td><td>This work</td><td>n/a</td><td>n/a</td><td>18.18 ± 0.61</td><td>Sr, S</td><td>0.27</td><td>Misclassification</td></tr><tr><td>13101</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>10.81 ± 0.47</td><td>Cg, Xc</td><td>0.28</td><td>Athor</td></tr><tr><td>14950</td><td>Zita</td><td>This work</td><td>n/a</td><td>yes</td><td>23.17 ± 0.34</td><td>X, Xk</td><td>0.19</td><td>Zita</td></tr><tr><td>15996</td><td>Zita</td><td>This work</td><td>n/a</td><td>yes</td><td>23.39 ± 0.32</td><td>T, X</td><td>0.14</td><td>Zita</td></tr><tr><td>16007</td><td>Athor core</td><td>This work</td><td>This work</td><td>yes</td><td>19.45 ± 0.95</td><td>T, X</td><td>0.11</td><td>Zita</td></tr><tr><td>17635</td><td>Zita</td><td>This work</td><td>n/a</td><td>yes</td><td>19.54 ± 0.24</td><td>X, Xk</td><td>0.11</td><td>Zita</td></tr><tr><td>19207</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>5.53 ± 0.17</td><td>Xc, Xe</td><td>0.16</td><td>Athor</td></tr><tr><td>20566</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>8.35 ± 0.40</td><td>Xc, Cgh</td><td>0.19</td><td>Athor</td></tr><tr><td>21475</td><td>Zita</td><td>This work</td><td>n/a</td><td>yes</td><td>23.89 ± 0.57</td><td>Sq, Q</td><td>0.13</td><td>Misclassification</td></tr><tr><td>33110</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>3.29 ± 0.14</td><td>Xn, Xc</td><td>0.25</td><td>Athor</td></tr><tr><td>34120</td><td>Zita</td><td>This work</td><td>This work</td><td>yes</td><td>18.26 ± 0.24</td><td>Xk, Cb</td><td>0.24</td><td>Zita</td></tr><tr><td>37983</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>10.04 ± 0.18</td><td>Xe, Xc</td><td>0.14</td><td>Athor</td></tr><tr><td>38233</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>4.79 ± 0.15</td><td>Xc, Xe</td><td>0.12</td><td>Athor</td></tr><tr><td>54355</td><td>Zita</td><td>This work</td><td>n/a</td><td>yes</td><td>19.73 ± 0.54</td><td>Xk, X</td><td>0.14</td><td>Zita</td></tr><tr><td>68750</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>13.74 ± 0.67</td><td>C, Cg</td><td>0.10</td><td>Misclassification</td></tr><tr><td>74119</td><td>Athor core</td><td>This work</td><td>n/a</td><td>yes</td><td>11.18 ± 0.89</td><td>Xe, Xc</td><td>0.27</td><td>Athor</td></tr><tr><td colspan="9" style="text-align: left; font-style: italic;">Literature</td></tr><tr><td>161</td><td>Athor core</td><td>(d)</td><td>(d)</td><td>n/a</td><td>7.94 ± 0.07</td><td>Xc</td><td>0.23</td><td>Athor</td></tr><tr><td>757</td><td>Athor core</td><td>(e)</td><td>(a)</td><td>yes</td><td>10.79 ± 0.26</td><td>Xk</td><td>0.22</td><td>Athor</td></tr><tr><td>1490</td><td>Athor halo</td><td>n/a</td><td>(a)</td><td>n/a</td><td>-</td><td>Xc</td><td>0.17</td><td>Athor</td></tr><tr><td>1697</td><td>Athor halo</td><td>n/a</td><td>(b)</td><td>n/a</td><td>-</td><td>TX</td><td>0.15</td><td>Zita</td></tr><tr><td>2194</td><td>Athor core</td><td>n/a</td><td>(a)</td><td>yes</td><td>-</td><td>Xc</td><td>0.21</td><td>Athor</td></tr><tr><td>2328</td><td>Zita</td><td>n/a</td><td>(a)</td><td>yes</td><td>-</td><td>C</td><td>0.11</td><td>Misclassification</td></tr><tr><td>2991</td><td>Zita</td><td>n/a</td><td>(c)</td><td>yes</td><td>-</td><td>C</td><td>0.16</td><td>Misclassification</td></tr><tr><td>3137</td><td>Zita</td><td>n/a</td><td>(a)</td><td>yes</td><td>-</td><td>C</td><td>0.21</td><td>Misclassification</td></tr><tr><td>3375</td><td>Zita</td><td>n/a</td><td>(a)</td><td>yes</td><td>-</td><td>C</td><td>0.12</td><td>Misclassification</td></tr><tr><td>3865</td><td>Athor halo</td><td>n/a</td><td>(a)</td><td>n/a</td><td>-</td><td>Xc</td><td>0.18</td><td>Athor</td></tr><tr><td>4838</td><td>Zita</td><td>n/a</td><td>(a)</td><td>n/a</td><td>-</td><td>Xc</td><td>0.14</td><td>Athor</td></tr><tr><td>4839</td><td>Athor halo</td><td>n/a</td><td>(a)</td><td>n/a</td><td>-</td><td>Xc</td><td>0.20</td><td>Athor</td></tr><tr><td>4845</td><td>Athor halo</td><td>n/a</td><td>(a)</td><td>n/a</td><td>-</td><td>Xc</td><td>0.18</td><td>Athor</td></tr></tbody></table>

**Notes.** The spectral classes are derived after the classification of each combined spectrum in the Bus-DeMeo taxonomy; in bold we indicate the preferred class. The original family membership and the $p_V$ s are from Delbo et al. (2019), while the revised membership is derived after the new spectroscopic observations in combination with the orbital elements inspection. Refs.: (a) Bus & Binzel (2002), (b) Xu et al. (1995), (c) Lazzaro et al. (2004), (d) Fornasier et al. (2010), (e) Clark et al. (2009)

bands). A further matching parameter is the absolute value of the reflectance, which is typically calculated at 0.56 µm. We addressed this matching at a second step. Initially, using an established procedure (Eq. 2 from Beck et al. 2021), we corrected the meteorite absolute reflectances, which are typically measured with an emission angle of 30° and incident illuminating

angle 0° (phase angle 30°) to 0° of phase in order to compare them with the geometric visible albedo $p_V$ of the Athor family members. In addition, we computed the 1.8/0.8 µm reflectance ratio for the meteorites and the Athor family members in order to construct compositionally diagnostic plots (Fig. 1 B); following a well established method (Cloutis et al. 2010), we built
<figure><img src="image_7.png" /><figcaption>**Fig. B.2.** NIR slopes of the asteroid members of the Athor and Zita families appear to have a bimodal distribution. The separation at around 15% µm⁻¹ distinguishes the members of these two overlapping families, with the Zita members having the redder slopes.</figcaption></figure>

<table><caption>**Table B.3.** Meteorites that satisfied the curve matching with the average spectrum of the Athor family.</caption><thead><tr><th>File name</th><th>Database</th><th>Group</th></tr></thead><tbody><tr><td>ccmt40.txt</td><td>RELAB</td><td>EH</td></tr><tr><td>c1mt122c.txt</td><td>RELAB</td><td>EH</td></tr><tr><td>90814c.011.txt</td><td>PSF</td><td>CV</td></tr><tr><td>camt40.txt</td><td>RELAB</td><td>EH</td></tr><tr><td>camp226.txt</td><td>RELAB</td><td>CM</td></tr><tr><td>chmt40.txt</td><td>RELAB</td><td>EH</td></tr><tr><td>clmb64.txt</td><td>RELAB</td><td>CM</td></tr><tr><td>091104a.016.txt</td><td>PSF</td><td>EH</td></tr><tr><td>90814c.012.txt</td><td>PSF</td><td>CV</td></tr><tr><td>s1tb46.txt</td><td>RELAB</td><td>Aubrites</td></tr><tr><td>cgn021.txt</td><td>RELAB</td><td>EL</td></tr><tr><td>c1mp226a2.txt</td><td>RELAB</td><td>CM</td></tr><tr><td>mgn021.txt</td><td>RELAB</td><td>EL</td></tr><tr><td>c1mt122d.txt</td><td>RELAB</td><td>EH</td></tr><tr><td>cgn017.txt</td><td>RELAB</td><td>EL</td></tr><tr><td>cgn019.txt</td><td>RELAB</td><td>EL</td></tr><tr><td>mgn017.txt</td><td>RELAB</td><td>EL</td></tr></tbody></table>

the 1.8/0.8 µm reflectance ratio vs. 0.56 µm absolute reflectance plot for the reflectance spectra of those carbonaceous chondrites passing the curve-matching step and for all EL and EH reflectance spectra in the RELAB and PSF databases including those passing the curve-matching step. The filenames of the curve-matched meteorites are given in Table B.3.

## Appendix C: Dynamical and collisional loss of Athor family members

Firstly, we estimated the size-dependent probability that an Athor family member is dynamically removed from main belt. To do so, we integrated the orbital evolution of a randomised synthetic Athor family for 3 Gyr, monitoring how many of these bodies leave the inner main belt. The code used is a modified version of SWIFT RMVS3 (Levison & Duncan 1994). Our orbital integrations included all planets except Mercury on their current orbits and also the orbital drift due to the Yarkovsky effect for

all the asteroids. The large asteroids, such as Ceres and Vesta, were not included in the integrations because it has been shown (Delbo et al. 2017) that their effect in the dynamical depletion is negligible for our purposes. We modelled 375 asteroids whose diameters (D) were selected randomly from a cumulative size frequency distribution (SFD) given by Eq. C.1

$$N(> D) = N_0 D^{\alpha}, \qquad (C.1)$$

where N is the number of asteroids larger than D, $N_0 = D_{Max}^{-\alpha}$, α the cumulative slope set equal to -3, and $3 < D < 50$ km (implying that $D_{Max} = 50$ km). The cumulative slopes of asteroid families have been measured and modelled to be steeper (Tanga et al. 1999; Rozehnal et al. 2022) than the -2.5, meaning steeper than the collisional equilibrium of Dohnanyi (1969). On the other hand, cumulative family size slopes that are steeper than -3 lead to a mass function that cannot be integrated down to dust sizes. Therefore, we chose to run the simulation with $\alpha=-3$. The selection of $D_{Max} = 50$ km was done in order to account for the current largest fragment of the family; asteroid (161) Athor with $D = 42 \pm 4$ km (Alí-Lagoa et al. 2018). Each body was given an orbital drift rate (da/dt, Vokrouhlický et al. 2015) depending on its diameter and on random obliquity (ε) with a rate of $1.93 \times 10^{-4} \times \cos(\epsilon) \times (1 \text{ km/D}^{-1}) \text{ au Myr}^{-1}$. A 10-day time step was used for the integrations. Asteroids were eliminated from the simulation if their distance from the Sun reached less than 0.5 au or greater than 30 au or if they passed within a Hill Sphere of a planet (the area where the planet's gravity dominates). All the asteroids had initial averaged orbital elements (averaged over 20 Myr to approximate proper elements) within the range of eccentricity $0.055 < e < 0.1$, and inclination $7.15^{\circ} < i < 9.2^{\circ}$, typical of the observed Athor family. The initial semi-major axes were centred around $a = 2.379$ au (the semi-major axis of Athor), and an initial range of semi-major axes was calculated for an initial family velocity dispersion of 65 m s⁻¹. After the 3 Gyr simulation, 254 of the initial 375 bodies were still on orbits in the inner main belt allowing us to estimate the size-dependent probability of the dynamical loss (1 minus the ratio between the initial and the final number of asteroids) in different diameter bins.

Given the above, we were able to estimate the dynamical and collisional evolution of a synthetic Athor family using a probabilistic approach: We began by creating an initial family of 1000 asteroids whose diameters were randomly drawn from a cumulative SFD whose form is given by Eq. C.1. We use $\alpha_{ini}$ to refer to the initial slope of the SFD. To apply the dynamical loss we performed the following. For each asteroid we extracted a random number from a uniform distribution between 1 and 0. If this number was smaller than the probability of the dynamical loss of said asteroid (calculated at the previous step), then the latter is eliminated. This procedure resulted in an erosion of the initial SFD. Next, we applied a second SFD erosional step to simulate the elimination of asteroids due to collisions. For this, we adopted a classical method that consists in determining the probability that an asteroid undergoes a catastrophic collision during a time $\Delta t$. This probability is equal to the size-dependent inverse collisional lifetime from Bottke et al. (2005), multiplied by the time $\Delta t$. Subsequently, we ran a Monte Carlo simulation for a number of steps equal to $3 \text{ Gyr}/\Delta t$ (where $\Delta t=10$ Myr). At each step, an asteroid is eliminated from the simulation if a random number, uniformly distributed between 0 and 1, was found smaller than that probability. Next, we measured the dynamical and collisional loss of family members by adjusting the parameters $\alpha$ and $N_0$ to obtain the best fit of Eq. C.1 to the post-evolutionary SFD and we call $\alpha_{act}$ this best fit value of $\alpha$. We repeated the
<figure><img src="image_8.png" /><figcaption>**Fig. C.1.** Cumulative size distribution approximated by a power law, $N(>D) = N_0D^\alpha$, where the $\alpha$-value is the slope of a straight line in the $\log_{10}N - \log_{10}D$ space. Here we show the actual slope, $\alpha_{act}$, and its standard error (grey area), of a synthetic Athor family resulting from the dynamical and collisional evolution of an initial distribution of family members with slope $\alpha_{ini}$ for 3 Gyr. The horizontal dashed line denotes the actual observed slope of the real Athor family (-1.67). The two curves cross at $\alpha_{ini}=-2.37$, which indicates the initial slope of the cumulative size distribution of the family.</figcaption></figure>

Monte Carlo simulations and fitting of Eq. C.1 1,000 times for discrete values of $\alpha_{ini}$ between -3.7 and -1.7 with a step of 0.1;

for each value of $\alpha_{ini}$ we recorded the mean value of $\alpha_{act}$ and its standard error. The result is presented in Fig. C.1.

The $\alpha$-value of the real Athor family was derived using the diameters of its members from Delbo et al. (2019) after having removed the interlopers spectroscopically identified here (Table B.2). Using the $\alpha$-value of the real Athor family as the $\alpha_{act}$ (-1.67) we indicated the most likely $\alpha_{ini}$-slope of the SFD of the family to be -2.37.

### C.1. Initial size of the parent body of the Athor family

We determined the size of the precursor to the Athor family by integrating the initial SFD of this latter given by Eq. C.1 with $\alpha$ equal to the best fit $\alpha_{ini}$ (Fig. C1), namely:

$$V_p = \frac{\pi}{6} \alpha_{ini} N_0 \int_{D_{\text{Max}}}^{0} D^{\alpha_{ini}+2} dD, \qquad (C.2)$$

where $V_p$ is the volume of Athor family precursor body and $D_{\text{Max}}$ is the diameter of the largest family member, which is (161) Athor of ~42 km in diameter (Alí-Lagoa et al. 2018). Equation C.2 can be solved analytically for $\alpha_{ini} > -3$, which is our case, yielding the diameter of the Athor family progenitor (~64 km):

$$D_p = \left(-\alpha_{ini} N_0 \frac{1}{\alpha_{ini} + 3} D^{\alpha_{ini} + 3}\right)^{\frac{1}{3}}. \qquad (C.3)$$