# The Appearance of a "Fresh" Surface on 596 Scheila as a Consequence of the 2010 Impact Event

Sunao Hasegawa, Michaël Marsset, Francesca E. Demeo, Schelte J. Bus, Masateru Ishiguro, Daisuke Kuroda, Richard P. Binzel, Josef Hanuš, Akiko M. Nakamura, Bin Yang, et al.

► To cite this version:

Sunao Hasegawa, Michaël Marsset, Francesca E. Demeo, Schelte J. Bus, Masateru Ishiguro, et al. The Appearance of a "Fresh" Surface on 596 Scheila as a Consequence of the 2010 Impact Event. The Astrophysical journal letters, 2022, 924, 10.3847/2041-8213/ac415a. insu-03667485

HAL Id: insu-03667485

https://insu.hal.science/insu-03667485

Submitted on 13 May 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
# The Appearance of a “Fresh” Surface on 596 Scheila as a Consequence of the 2010 Impact Event

Sunao Hasegawa<fnref n="1" />, Michël Marsset<fnref n="2" />,<fnref n="3" />, Francesca E. DeMeo<fnref n="2" />, Schelte J. Bus<fnref n="4" />, Masateru Ishiguro<fnref n="5" />,<fnref n="6" />, Daisuke Kuroda<fnref n="7" />, Richard P. Binzel<fnref n="2" />, Josef Hanuš<fnref n="8" />, Akiko M. Nakamura<fnref n="9" />, Bin Yang<fnref n="3" />,<fnref n="10" />, and Pierre Vernazza<fnref n="11" />

Received 2021 November 1; revised 2021 November 29; accepted 2021 December 7; published 2022 January 4

## Abstract

Dust emission was detected on main-belt asteroid 596 Scheila in 2010 December and was attributed to the collision of a few-tens-of-meters projectile on the surface of the asteroid. In such an impact, the ejected material from the collided body is expected to mainly come from its fresh, unweathered subsurface. Therefore, it is expected that the surface of 596 was partially or entirely refreshed during the 2010 impact. By combining spectra of 596 from the literature and our own observations, we show that the 2010 impact event resulted in a significant slope change in the near-infrared (0.8–2.5 µm) spectrum of the asteroid, from moderately red (T type) before the impact to red (D type) after the impact. This provides evidence that red carbonaceous asteroids become less red with time due to space weathering, in agreement with predictions derived from laboratory experiments on the primitive Tagish Lake meteorite, which is spectrally similar to 596. This discovery provides the very first telescopic confirmation of the expected weathering trend of asteroids spectrally analog to Tagish Lake and/or anhydrous chondritic porous interplanetary dust particles. Our results also suggest that the population of implanted objects from the outer solar system is much larger than previously estimated in the main belt, but many of these objects are hidden below their space-weathered surfaces.

Unified Astronomy Thesaurus concepts: Small Solar System bodies (1469); Asteroids (72); Main belt asteroids (2036); Asteroid surfaces (2209)

Supporting material: data behind figure

## 1. Introduction

596 Scheila, a dark, outer main-belt asteroid with a diameter of 114 km [mean value obtained from the Infrared Astronomical Satellite (Tedesco et al. 2002), AKARI (Usui et al. 2011), and the Wide-field Infrared Survey Explorer (Masiero et al. 2011)] discovered in 1906, was found to exhibit a comet-like appearance on 2010 December 11, about 100 yr after its discovery (Larson et al. 2010). 596 is rather unique in that it is one of the only two main-belt asteroids larger than 100 km in diameter for which cometary-like activity has been detected, the second one being the dwarf planet 1 Ceres (Küppers et al. 2014). All other main-belt comets previously discovered are less than 10 km in diameter (Jewitt 2012). Spectroscopic studies of 596 conducted after it became active appeared to suggest a distinct mechanism for dust emission compared to comets and other main-belt comets. No emission lines originating from volatile components were detected in its UV-visible spectrum from the Swift UV-optical telescope

(Bodewits et al. 2011), the visible spectra by FORS2 onboard the Very Large Telescope (Jehin et al. 2011) and the LRIS-equipped Keck I (Hsieh et al. 2012), the near-infrared (NIR) spectrum by Spex on the InfraRed Telescope Facility (IRTF; Yang & Hsieh 2011), and the microwave spectrum by the Arecibo Observatory (Howell & Lovell 2011). In addition, no coma derived from gas components was detected in the narrow-band images of 596 (Jehin et al. 2011). From spectroscopy in the 3 µm region, the surface of 596 was found to have no more than a few percent of water ice (Yang & Hsieh 2011). These observational results suggested that the cometary activity of 596 is not due to ice sublimation, although they do not themselves establish that the activity of 596 was not due to ice sublimation considering that the production rates of sublimation products could have been below the detection limits of the observations.

Based on analysis of the coma and dust tail distributions of 596, Jewitt et al. (2011), Moreno et al. (2011), Ishiguro et al. (2011a), and Hsieh et al. (2012) confirmed that the cometary activity was not caused by ice sublimation, but by the impact of a projectile. Ishiguro et al. (2011b) further succeeded in perfectly reproducing the observed dust tail morphology of 596 from the model of ejecta distribution during an oblique impact

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

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

<footnote n="1">¹Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, 3-1-1 Yoshinodai, Chuo-ku, Sagamihara, Kanagawa 252-5210, Japan hasehase@isas.jaxa.jp</footnote>

<footnote n="2">²Department of Earth, Atmospheric, and Planetary Sciences, MIT, 77 Massachusetts Avenue, Cambridge, MA 02139, USA</footnote>

<footnote n="3">³European Southern Observatory (ESO), Alonso de Córdova 3107, 1900 Casilla Vitacura, Santiago, Chile</footnote>

<footnote n="4">⁴Institute for Astronomy, University of Hawaii, 2860 Woodlawn Drive, Honolulu, HI 96822-1839, USA</footnote>

<footnote n="5">⁵Department of Physics and Astronomy, Seoul National University, Gwanak-gu, Seoul 08826, Republic of Korea</footnote>

<footnote n="6">⁶SNU Astronomy Research Center, Seoul National University, Gwanak-gu, Seoul 08826, Republic of Korea</footnote>

<footnote n="7">⁷Okayama Observatory, Kyoto University, 3037-5 Honjo, Kamogata-cho, Asakuchi, Okayama 719-0232, Japan</footnote>

<footnote n="8">⁸Institute of Astronomy, Faculty of Mathematics and Physics, Charles University, V Holešovičkách 2, 180 00 Prague 8, Czech Republic</footnote>

<footnote n="9">⁹Department of Planetology, Graduate School of Science, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan</footnote>

<footnote n="10">¹⁰Núcleo de Astronomía, Facultad de Ingeniería Ciencias, Universidad Diego Portales, Chile</footnote>

<footnote n="11">¹¹Aix Marseille Université, CNRS, LAM, Laboratoire d'Astrophysique de Marseille, Marseille, France</footnote>
observed in laboratory experiments. The authors were also able
to determine the impactor size, impact timing and direction, the
size of the resulting crater, and the tensile strength of 596's
surface material. The ejecta discharged from the impact crater
was predicted to completely or partially cover the surface of
596 (Ishiguro et al. 2011a; Bodewits et al. 2014). Licandro
et al. (2011) obtained rotationally resolved NIR spectra after
the impact event and showed that the postimpact spectra in the
0.8–2.4 µm wavelength region are homogeneous, meaning that
the subsurface material ejected from the impact crater likely
entirely resurfaced 596.

Our team is currently conducting a visible and NIR spectroscopic survey of main-belt asteroids with diameters larger than 100 km to investigate the physical properties of planetesimals that existed in the early solar system (see Introduction chapter in Hasegawa et al. 2021). In the course of compiling near-infrared data obtained in the past by the MIT–Hawaii Near-Earth Object Spectroscopic Survey (MITHNEOS),<fnref n="12" /> we found out that spectra of 596 from the 0.8 to 2.5 µm wavelength region had been obtained before and after the 2010 impact as part of the survey. So far, studies of the outcome of the impact on the physical properties of 596 have only been performed by means of visible lightcurve observations (e.g., Ishiguro et al. 2011a; Bodewits et al. 2014), and none of these studies focused on how the reflectance spectrum of 596 changed. The purpose of this study is to investigate whether or not a change happened in the spectrum of 596 from 0.8 to 2.5 µm as a consequence of the impact event, as well as in the visible and 3 µm regions by studying literature. Then, we discuss the factors that may have contributed to this change.

## 2. Observations and Data Analysis

NIR spectroscopic observations of 596 Scheila were
conducted at two distinct epochs: 2002 June 1 and 2011
February 7, with the 3 m IRTF located on Maunakea, Hawaii.
Data from the first epoch of observation were published and
analyzed by DeMeo et al. (2009) while data from the second
epoch were never published. We describe here the observations
and data reduction performed for the data acquired on 2011
February 7.

We used the SpeX NIR spectrograph (Rayner et al. 2003) combined with a 0.′′8 × 15′′ slit in the low-resolution prism mode to measure the spectra over the 0.8–2.5 µm wavelength range. At the time of our observations, 596 had a visual V-band magnitude of 13.4 and was located at an air mass of 1.1. A total of 15 spectral images with 120 s exposure time were recorded in an AB beam pattern to allow efficient removal of the sky background by subtracting pairs of AB images. The following calibration stars, known to be very close spectral analogs to the Sun, were observed on the same night: Hyades 64 and Landolt (1983)'s stars 98–978 and 102–1081. An in-depth analysis of these calibration stars is provided in Marsset et al. (2020).

Data reduction and spectral extraction followed the proce-
dure outlined in Binzel et al. (2019). We summarize it briefly
here. Reduction of the spectral images was performed with the
Image Reduction and Analysis Facility (IRAF; Tody 1993) and
Interactive Data Language (IDL), using the Autospex software
tool to automatically write sets of command files (Rivkin et al.
2005). Reduction steps for the science target and its
corresponding calibration stars included trimming the images,

creating a bad pixel map, flat-fielding the images, sky-
subtracting between AB image pairs, tracing the spectra in
both the wavelength and spatial dimensions, co-adding the
spectral images, extracting the 2D spectra, performing
wavelength calibration, and correcting for air-mass differences
between the asteroid and the solar analogs. Finally, the
resulting asteroid spectrum was divided by the mean stellar
spectra to remove the solar gradient.

3. Spectroscopic Results

We compared the spectra of 596 Scheila obtained before and after the 2010 impact to search for any spectral change induced by this event. Spectral data obtained before the impact event were retrieved from Bus & Binzel (2002), DeMeo et al. (2009), and Licandro et al. (2011). NIR spectra collected after the impact event are from Yang & Hsieh (2011) and Licandro et al. (2011), as well as from this study. Table 1 shows the subobserver coordinates of the asteroid during the pre- and postimpact NIR observations calculated based on the shape model from Hanuš et al. (2021). Since a visible reflectance spectrum is not available after the impact, the average of the multicolor photometry of Trigo-Rodriguez et al. (2011), Betzler et al. (2012), Hsieh et al. (2012), Jewitt (2012), and Shevchenko et al. (2016) was utilized instead.

Marsset et al. (2020) estimated the uncertainty of spectral slope measurements of asteroids in the NIR wavelength region from 20 yr of MITHNEOS observations with the SpeX spectrograph on the IRTF, and showed that the uncertainty is 4.2% µm⁻¹. Spectra obtained by Yang & Hsieh (2011) after the impact event, which are not from MITHNEOS observations, and the one obtained in this study, which is from MITHNEOS observations, are consistent within that uncertainty range (Figure 1). The total exposure time of the two spectra acquired in Yang & Hsieh (2011) is 400 s each, and the total exposure time of the pre- and postimpact spectra acquired by MITHNEOS observations is 840 s and 1800 s (DeMeo et al. 2009 and this study), respectively: these spectra are data of comparable quality with signal-to-noise ratios of several hundreds. In addition, the spectra collected after the impact event by Licandro et al. (2011) show that the variance of the slope is about 10%, while the spectrum obtained by Yang & Hsieh (2011) is within this range (Figure 1). Considering the dispersion and uncertainty mentioned above, the agreement of the spectra after the impact event is consistent with the point made by Licandro et al. (2011) that the NIR spectra of 596 does not change with rotation phase. The spectrum of Yang & Hsieh (2011) in Figure 1 is a combination of two days of data, and the difference between the two data sets is several percent. This is also an indication of the homogeneity of the surface of the 596 after the impact event.

Comparing the preimpact spectrum of DeMeo et al. (2009) obtained by MITHNEOS with the postimpact spectrum of this study also obtained by MITHNEOS shows that the difference is significantly larger than the indeterminacy obtained by Marsset et al. (2020), which indicates that the spectrum changed before and after the impact event (Figure 1). The slope of 596 at 0.8–2.5 µm before and after the impact event is 25.5 and 46.7% µm⁻¹, respectively. In fact, the spectral type of 596 before the impact event was T type (DeMeo et al. 2009), but after the impact event it changed to D type (Yang & Hsieh 2011). This change is unlikely to be due to varying observing conditions. In particular, the two IRTF spectra

<footnote n="12">¹² http://smass.mit.edu/</footnote>
<table><caption>**Table 1**

List of the Subobserver Coordinates of 596 Scheila at the Time of NIR Observations</caption><thead><tr><th rowspan="2">Observation Time</th><th rowspan="2">Timing of Observation</th><th colspan="4">Sub-Earth Point</th><th rowspan="2">Spectral Slope (% µm<sup>-1</sup>)</th><th rowspan="2">References</th></tr><tr><th>λ<sub>1</sub><sup>a</sup> (°)</th><th>φ<sub>1</sub><sup>a</sup> (°)</th><th>λ<sub>1</sub><sup>b</sup> (°)</th><th>φ<sub>1</sub><sup>b</sup> (°)</th></tr></thead><tbody><tr><td>2002/06/01 09:40</td><td>before impact</td><td>137</td><td>46</td><td>232</td><td>-47</td><td>25.5</td><td>DeMeo et al. (2009)</td></tr><tr><td>2010/12/10 10:58</td><td>after impact</td><td>157</td><td>-40</td><td>292</td><td>41</td><td>36.5<sup>c</sup></td><td>Licandro et al. (2011)</td></tr><tr><td>2010/12/10 14:35</td><td>after impact</td><td>75</td><td>-40</td><td>210</td><td>41</td><td>36.5<sup>c</sup></td><td>Licandro et al. (2011)</td></tr><tr><td>2010/12/10 15:02</td><td>after impact</td><td>64</td><td>-40</td><td>200</td><td>41</td><td>36.5<sup>c</sup></td><td>Licandro et al. (2011)</td></tr><tr><td>2011/01/04 14:50</td><td>after impact</td><td>132</td><td>-39</td><td>263</td><td>43</td><td>40.5<sup>d</sup></td><td>Yang & Hsieh (2011)</td></tr><tr><td>2011/01/05 14:20</td><td>after impact</td><td>319</td><td>-39</td><td>90</td><td>44</td><td>40.5<sup>d</sup></td><td>Yang & Hsieh (2011)</td></tr><tr><td>2011/01/07 02:56</td><td>after impact</td><td>209</td><td>-39</td><td>339</td><td>44</td><td>36.5<sup>c</sup></td><td>Licandro et al. (2011)</td></tr><tr><td>2011/01/08 06:22</td><td>after impact</td><td>306</td><td>-39</td><td>76</td><td>44</td><td>36.5<sup>c</sup></td><td>Licandro et al. (2011)</td></tr><tr><td>2011/01/09 06:50</td><td>after impact</td><td>110</td><td>-39</td><td>241</td><td>44</td><td>36.5<sup>c</sup></td><td>Licandro et al. (2011)</td></tr><tr><td>2011/01/10 01:53</td><td>after impact</td><td>39</td><td>-39</td><td>168</td><td>44</td><td>36.5<sup>c</sup></td><td>Licandro et al. (2011)</td></tr><tr><td>2011/01/10 07:02</td><td>after impact</td><td>282</td><td>-39</td><td>51</td><td>44</td><td>36.5<sup>c</sup></td><td>Licandro et al. (2011)</td></tr><tr><td>2011/02/07 09:31</td><td>after impact</td><td>99</td><td>-40</td><td>225</td><td>50</td><td>46.7</td><td>This study</td></tr></tbody></table>

**Notes.** There are two shape/pole solutions, which is a common thing in the inversion of disk-integrated data without stellar occultations or disk-resolved data.

<sup>a</sup> Pole solution 1: $\lambda_1 = 110$ (°), $\beta_1 = -24$ (°), Period = 15.8605 (hr).

<sup>b</sup> Pole solution 2: $\lambda_1 = 273$ (°), $\beta_1 = -38$ (°), Period = 15.8609 (hr).

<sup>c</sup> The value of the combined spectrum in Licandro et al. (2011).

<sup>d</sup> The value of the combined spectrum in Yang & Hsieh (2011).

obtained as part of MITHNEOS before and after the impact were taken at slightly different air-mass values (1.4 before the impact versus 1.1 after). Marsset et al. (2020) showed that air mass only contributes to a systematic slope change of $-0.92\% \ \mu m^{-1}$ per 0.1 unit air mass: this effect is insufficient to explain the slope difference between the two spectra.

Possible additional external factors causing slope variability of the asteroid may include varying weather conditions and close encounters of the asteroid with a bright star during the observations. However, the sky was reported to be clean during the night of the preimpact observations and no stellar encounter was confirmed based on a combination of the JPL HORIZONS ephemeris generator<fnref n="13" /> and the SIMBAD Astronomical Database.<fnref n="14" /> In addition, a search for bright (V-band magnitude $< 10.5$) stars using the SIMBAD Astronomical Database (see footnote 14) reveals no such stars at proximity (within $20'$) of 596 during the MITHNEOS observations. Finally, another asteroid observed on the same, preimpact night as 596, the S-type asteroid 57 Mnemosyne, was also observed on 2010 July 11 and exhibited an entirely consistent shape of its spectra between the two observing nights (see Figure A1).

Based on the change in the lightcurve of 596 recorded in the R-band filter before and after the impact event (Ishiguro et al. 2011a), there may be an increase or decrease in visible reflectance of up to 4.5%. However, the slopes of the visible spectra taken before and after the impact event are consistent within the observational error range, and there is no significant change in the visible slope before and after the impact event. Spectral observations of 596 in the 3 µm band were also made before and after the impact event, and the 3 µm spectra are consistent within the observational error range (Figure 2). Yang & Hsieh (2011) found that there was almost no ice on the surface of 596 through 3 µm spectra of 596 after the impact event. Similarly, the 3 µm spectrum of 596 obtained before the impact (Licandro et al. 2011) does not exhibit any water ice signature either.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Comparison of 596's NIR spectra before and after the 2010 impact event. The dotted and solid lines show the spectra before and after the impact event, respectively. The figure shows spectra of 596 in the 0.4-2.5 µm wavelength region, normalized to 1 at 0.55 µm. The black and blue arrows in the figure indicate the slope uncertainty from Marsset et al. (2020), and the dispersion range of the spectra obtained by Licandro et al. (2011), centered at the reflectance value of the spectra at 2.35 µm as shown by the black dashed lines.</figcaption></figure>

(The data used to create this figure are available.)

## 4. Discussion

The slope of the NIR spectrum of 596 Scheila after the impact, with $46.7 \ \mu m^{-1}$, is redder than the spectrum before the impact event, with $25.5 \ \mu m^{-1}$. This section discusses the factors that may have contributed to this.

The first possibility is surface heterogeneity. Table 1 shows the subobserver coordinates of the pre- and postimpact observations for two different solutions of the pole orientation. Since 596 was observed before the impact in June, and after the impact around January, the direction of observation of 596

<footnote n="13">¹³ https://ssd.jpl.nasa.gov/horizons/</footnote>

<footnote n="14">¹⁴ http://simbad.u-strasbg.fr/simbad/</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Comparison of the 3 µm spectra of 596 before and after the 2010 impact event. The dotted and solid lines show the spectra before and after the impact event, respectively. The figure shows the spectra in the 1.8–4.0 µm wavelength regions, normalized to 1 at 2.2 µm.</figcaption></figure>

from the Earth is different by about 180 degrees. Sub-Earth latitude is reversed because the observation period before and after the impact happens to be off by half a year. The homogeneity of the postimpact spectra indicates that the southern hemisphere,<fnref n="15" /> observed after the impact, exhibits a homogeneous surface. Comparing the surface of the southern hemisphere (see footnote 15) observed after the impact and the surface of the northern hemisphere (see footnote 15) observed before the impact may be heterotropic. However, considering that most of the material from the crater was expected to be ejected isotropically (e.g., Schultz & Gault 1985; Shuvalov 2011), and that the ejected dust covered the entire surface of 596 as it became comet-like (e.g., Ishiguro et al. 2011a, 2011b), we expect no difference between the northern and southern hemisphere spectra after the impact.

Note that the argument for partial resurfacing was based on the presence or absence of changes in the visible lightcurve before and after the impact event (Bodewits et al. 2014). However, the present study does not show any change in the visible spectrum before and after the impact, and the maximum albedo change is 4.5% (Ishiguro et al. 2011a). NIR spectra after the impact were obtained at longitudes with both many and little changes in the visible lightcurve before and after the impact event, but they show a similar spectral slope. Therefore, it is appropriate to assume that the spectral change occurred in the whole area observed by NIR wavelength region before and after the impact event. These results indicate that the only partial resurfacing of 596 inferred from visible lightcurve may be an underestimation of the true ejecta distribution. Since this study shows that the NIR region is more sensitive to spectral changes before and after the impact event than the visible wavelength region, it is natural to assume that ejecta were deposited over the entire surface of 596.

Also, the area where the impactor is thought to have collided with 596 (Ishiguro et al. 2011b) is at approximately the same

latitude as the one observed before the impact event (Figure 3). This shows that the surface material in the area observed before the impact must have fallen the area observed after the impact event as ejecta. In the case of spectral heterogeneity between the northern and southern hemispheres of 596, the spectra should have been similar to the spectrum taken before the impact, however, the pre- and postimpact spectra are different, indicating that geographic heterogeneity is not responsible for the observed spectral change of 596.

The second possibility is that the postimpact spectrum is redder due to a change in size of the particles that compose the surface layer of 596. The reddening of the spectrum after the impact suggests that finer particles may be present compared to before the impact (e.g., Johnson & Fanale 1973; Vernazza et al. 2016). Fine particles ejected from craters by the excavation flow can have a wide range of ejection velocities (Singer et al. 2020). However, finer particles have faster ejection velocities on average than larger particles, meaning that large particles should have preferentially been re-accreted at the surface of 596 after the impact. In general, if large sized particles remain, the spectrum should change to blue, but the observed change, a spectral reddening, is the opposite (see also Figure B1). Note that the original surface layer of 596 must also be composed of particles ejected from previous impact events, so the particle size distribution should basically be the same before and after the impact event. Due to thermal fatigue and meteoroid impacts (e.g., Delbo et al. 2014; Basilevsky et al. 2015), the surface particles of asteroids tend to become smaller as time goes by after large impacts such as the 2010 impact. Given that, it is not possible to explain spectral reddening by grain-size effect.

The third possibility is that the spectral slope change for the case of 596 may be due to the presence of the dust coma. The large dust particles that are inside the Hill radius of 596 and are not affected by radiation pressure may be around 596 for a long time. As shown in Ishiguro et al. (2011b), dust more than 140 µm may be mixed with the spectrum of 596 for a long time. (Dust less than 140 µm is stretched by radiation pressure on 2011 February when the observations for this study were made.) However, since about a few weeks after the impact the visible apparent magnitude of 596 approached the predicted value from before the impact (e.g., Jewitt et al. 2011; Trigo-Rodriguez et al. 2011), it is natural to assume that the scattering cross section of such a dust cloud is negligible compared to the entire surface of 596. Also, if the dust coma was affecting the spectrum, the spectrum from the observation closest to the collision event may be the reddest, but Figure 1 and Licandro et al. (2011) show that this is not the case.

Finally, space weathering is another mechanism that may have also played a role in the spectral change of 596 during the impact. The outer layer of asteroids is thought to undergo space weathering, where the spectrum is altered by micrometeorite impacts and solar wind irradiation (e.g., Clark et al. 2002; Brunetto et al. 2015). Most of the particles that cause cometary activity due to impacts or ice sublimation are thought to be ejected from the asteroid's subsurface rather than from the surface (e.g., impacts: Stoeffler et al. 1975; ice sublimation: A'Hearn et al. 2005). If covered by those particles, the surface layer after cometary activity would be by a fresh material that has not undergone space weathering.

Experiments simulating space weathering have been conducted on the C2-un grouped Tagish Lake meteorite, which has a similar spectrum to that of 596. In these experiments, the

<footnote n="15">¹⁵ This applies to the first shape/pole solution. For the second shape/pole solution, the northern and southern hemispheres are reversed.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Orbit diagrams of 596 at the time of preimpact NIR observation (left panel) and the 2010 impact event (right panel). Assuming an impact angle of 45°, which is the highest impact probability, the small asteroid would have collided with 596 from the direction (α, δ) = (60°, -40°) (Ishiguro et al. 2011b). The values in blue and red fonts are the values of the two pole solutions of 596 (Hanusš et al. 2021).</figcaption></figure>

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

spectrum tends to become bluer after space weathering<fnref n="16" /> (Hiroi et al. 2013; Vernazza et al. 2013; Lantz et al. 2017). Along these lines, Lantz et al. (2018) also stated that Tagish-Lake-like parent bodies are hidden in X- and T-type objects based on the results of laboratory space-weathering simulations. Since the spectrum of 596 before the 2010 impact event is bluer (T type; Bus & Binzel 2002; DeMeo et al. 2009), the fact that the surface of 596 before the impact was space weathered is in agreement with the claim of Lantz et al. (2018).

The maximum increase or decrease in absolute reflectance of 596 after the impact in the R-band is about 4.5% (Ishiguro et al. 2011a), which is very small compared to the spectral change obtained at minimum irradiance in the space-weathering experiment of Lantz et al. (2017). At minimum irradiance, the visible spectrum does not change at the 0.4–0.9 μm wavelength region (see the yellow area in Figure B1) and becomes bluer at 0.9–2.5 μm wavelength region (see the blue area in Figure B1). This trend in this laboratory experiment is highly consistent with the spectral change of 596 in the same wavelength range before and after the impact event. Note that the lightcurves acquired in R band were measured as relative values, and it is not possible to determine whether the absolute reflectance has increased or decreased after the impact only by comparing the lightcurves collected pre- and postimpact.

Previous examples of observations of fresh surfaces for small solar system bodies include: (1): P/2010 A2, which was shattered by a catastrophic disruption (Kim et al. 2012, 2017); (2): 6478 Gault, which is believed to have blown away its surface dust due to increased rotation (Marsset et al. 2019; Purdum et al. 2021), and (3): 162173 Ryugu, on the surface of which an artificial crater was formed by the Small Carry-on Impactor of the Hayabusa2 spacecraft (Arakawa et al. 2020; Morota et al. 2020). Due to the lack of spectroscopic observations acquired before the destruction of P/2010 A2, it is not known whether the spectrum of that object changed as a consequence of its disruption. In the case of 6478, important

spectral bluing occurred and was attributed to the loss of surface regolith (Marsset et al. 2019). Finally, the artificial crater experiment on the surface of the C-type asteroid Ryugu resulted in a slightly bluer visible light spectrum of that region of the asteroid (Morota et al. 2020). The different spectral trend observed for these objects compared to 596 is not entirely surprising considering their distinct spectral type (hence surface composition), environment (heliocentric distance, hence surface and subsurface thermal evolution), and resurfacing mechanisms. Our results provide the first spectroscopic evidence that the space-weathered surface layer of an asteroid was covered by fresh material from the subsurface during a natural impact event. Meanwhile, this is also the first time that actual astronomical observations show that D-type asteroids with spectra similar to the Tagish Lake meteorite become less red as a consequence of space weathering. Our results are also the first to provide astronomical evidence that Tagish-Lake parent bodies may be found among asteroids that are spectrally bluer than that meteorite, a prediction previously made by Lantz et al. (2018) based on laboratory experiments.

It is generally accepted that the Tagish Lake meteorite, which has been collected as a meteorite with a size large enough to conduct space weathering experiments, is of D-type asteroid origin. Vernazza et al. (2015) stated that most surfaces of dark C- and X-complex and D-type (including T-type) asteroids have not been collected as meteorites on Earth, but only as anhydrous chondritic porous interplanetary dust particles (IDPs). Hasegawa et al. (2017) also pointed out that in addition to their spectroscopic properties, the radar and polarization properties of these bodies indicate that they are related to chondritic porous IDPs. Most dark X-complex and D-type asteroid surfaces have been shown to have anhydrous surface compositions by spectroscopic observations in the 3 μm wavelength region (e.g., Takir & Emery 2012; Usui et al. 2019), and 596 also has anhydrous surface layers, which indicates an association with chondritic porous IDPs (Figure 2). Meanwhile, the Tagish Lake meteorite, which is aqueously altered, exhibits a 3 μm absorption feature (Hiroi et al. 2001), meaning that it may not be exactly analog in composition to 596. Vernazza et al. (2021) showed that large dark X-complex

<footnote n="16">¹⁶ The spectral slope of dark carbonaceous chondrites (CI and CM) including the Tagish lake meteorite becomes bluer due to space weathering, while that of bright carbonaceous chondrites (CO and CV) becomes redder.</footnote>
and D-type asteroids that were implanted inward of 3.3 au would have lost outer shell due to water-ice sublimation and impact erosion. However, since the erosion rate at the location of 596 is only 1 km over the age of the solar system (Vernazza et al. 2021), it is thought that 596 does not expose any of aqueously altered materials expected to be contained and abundant in its interior (Vernazza et al. 2017).

Our finding of a collision-induced change in the spectral slope and the spectral type of 596 supports the idea that more D types were originally implanted in the main belt than previously thought. This agrees with Vernazza et al. (2021)'s proposal that dark X-complex and D-type and some C-type asteroids may share the same origin in the transplanetary and/or trans-Neptunian planetesimal region. Most of the heavily weathered D-type asteroids may be very red objects, such as those discovered by Hasegawa et al. (2021), that have become bluer in their spectra as a result of space weathering. Our discovery therefore relaxes the observational constraints on dynamical models simulating the early implantation of outer solar system bodies in the asteroid belt. For instance, Levison et al. (2009) used the low observed number of dark X- and D-type asteroids as a stringent constraint to their simulations. Our results that some implanted objects may in fact be hidden among moderately red asteroids loosens that constraint.

Assuming that a collision event of the magnitude of the 2010 impact on 596 would entirely refresh the surface of an asteroid, the probability of an asteroid collision at 596's location in the main belt, the number of asteroids with similar size to the impactor of 596, and the cross sectional area of 596 yields a resurfacing timescale of $10^{3.5}-10^{4.5}$ yr (Jewitt 2012). Although it is not clear how much smaller an impact can be while still being able to resurface 596, this means that the typical timescale for surface renewal is shorter than $10^{3.5}-10^{4.5}$ yr.

## 5. Conclusions

We compared the spectra of 596 Scheila obtained before and after the 2010 impact event on that asteroid and showed that the NIR spectrum of its postimpact fresh surface was redder than the NIR spectrum of its space-weathered surface before the impact. This is consistent with the results of experiments simulating space weathering of the Tagish Lake meteorite, which has a spectrum similar to that of 596 after the impact. Our results are the first spectroscopic evidence of asteroid resurfacing by fresh material from its subsurface during a naturally occurring impact, and the first observational evidence that the spectra of D-type asteroids become less red owing to space weathering. Our results also imply that under space-weathering process, carbon-free, silicate-rich assemblages will redden whereas carbon-rich assemblages will be bluish. Finally, our results suggest that the number of implanted outer solar system objects in the main belt is much higher than previously estimated, but many of these objects are hidden below a space-weathered crust.

We thank the referee for a careful review and constructive suggestions, which helped us to improve the manuscript significantly. We are grateful to Dr. Cateline Lantz for sharing her Tagish Lake spectral data. We express our gratitude to Dr. Takahiro Hiroi for useful comments about the Tagish Lake data that he acquired at RELAB. We greatly appreciate Dr. Masahiko Arakawa and Dr. Tatsuhiro Michikami for useful comments about impact. This work is based on observations

<figure><img src="image_7.png" /><figcaption>Figure A1. Spectra of 57 Mnemosyne.</figcaption></figure>

collected at the Infrared Telescope Facility, which is operated by the University of Hawaii under contract 80HQTR19D0030 with the National Aeronautics and Space Administration. The authors acknowledge the sacred nature of Maunakea and appreciate the opportunity to observe from the mountain. The spectra of the Tagish Lake meteorite used in this study are from the RELAB database of Brown University. This study has utilized the SIMBAD database, operated at CDS, Strasbourg, France, and the JPL HORIZONS ephemeris generator system, operated at JPL, Pasadena, USA. M.M. and F.D. were supported by the National Aeronautics and Space Administration under grant No. 80NSSC18K0849 and 80NSSC18K1004 issued through the Planetary Astronomy Program. M.I. was supported by the NRF grant No. 2018R1D1A1A09084105. The work of J.H. has been supported by the Czech Science Foundation through grant 20-08218S. This study was supported by JSPS KAKENHI (grant Nos. JP18K03723, JP19H00719, JP20K04055, JP21H01140, and JP21H01148) and by the Hypervelocity Impact Facility (former facility name: the Space Plasma Laboratory), ISAS, JAXA.

Facility: IRTF:3.0m.

Software: IRAF (Tody 1993), IDL.

## Appendix A
Spectra of 57 Mnemosyne

The NIR spectrum of 57 Mnemosyne was obtained on the same day that the spectrum of 596 was obtained before the impact event (DeMeo et al. 2009). The spectra were acquired on different days, and both spectra are consistent (Figure A1). This shows the homogeneity of surface of 57, and also indicate that the observing conditions were stable and good on the day the two spectra were acquired.

## Appendix B
Spectra of the Tagish Lake Meteorite

Spectroscopic data of the Tagish Lake meteorite is provided for reference. Powder data and pellet data of the Tagish Lake meteorite were obtained from the Reflectance Experiment
<figure><img src="image_8.png" /><figcaption>**Figure B1.** Spectra of the Tagish Lake meteorite. Yellow and cyan areas indicate regions where the slope of the spectrum did not change and regions where it did change after minimum irradiation, respectively.</figcaption></figure>

Laboratory (RELAB)<fnref n="17" /> at Brown University and Lantz et al. (2017), respectively (Figure B1). Pellet samples exhibit a flat spectral slope and are likely not representative of asteroid regolith. In general, the spectra of pellet samples are much bluer than those of powder samples (e.g., Hasegawa et al. 2019).

## ORCID iDs

Sunao Hasegawa https://orcid.org/0000-0001-6366-2608

Michaël Marsset https://orcid.org/0000-0001-8617-2425

Francesca E. DeMeo https://orcid.org/0000-0002-8397-4219

Schelte J. Bus https://orcid.org/0000-0003-4191-6536

Masateru Ishiguro https://orcid.org/0000-0002-7332-2479

Daisuke Kuroda https://orcid.org/0000-0002-7363-187X

Richard P. Binzel https://orcid.org/0000-0002-9995-7341

Josef Hanuš https://orcid.org/0000-0002-2934-3723

Akiko M. Nakamura https://orcid.org/0000-0001-6990-8496

Bin Yang https://orcid.org/0000-0002-5033-9593

Pierre Vernazza https://orcid.org/0000-0002-2564-6743

## References

* A'Hearn, M. F., Belton, M. J. S., Delamere, W. A., et al. 2005, *Sci*, **310**, 258

* Arakawa, M., Saiki, T., Wada, K., et al. 2020, *Sci*, **368**, 67

* Basilevsky, A. T., Head, J. W., Horz, F., & Ramsley, K. 2015, *P&SS*, **117**, 312

* Betzler, A. S., Novaes, A. B., Santos, A. C. P., et al. 2012, *MPBu*, **39**, 5

* Binzel, R. P., DeMeo, F. E., Turtelboom, E. V., et al. 2019, *Icar*, **324**, 41

* Bodewits, D., Kelley, M. S., Li, J. Y., et al. 2011, *ApJL*, **733**, L3

* Bodewits, D., Vincent, J. B., & Kelley, M. S. P. 2014, *Icar*, **229**, 190

* Brunetto, R., Loeffler, M. J., Nesvorný, D., Sasaki, S., & Strazzulla, G. 2015, in *Asteroid Surface Alteration by Space Weathering Processes*, Asteroids

* IV, ed. P. Michel, F. E. DeMeo, & W. F. Bottke (Tucson, AZ: Univ. Arizona Press), 597

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

* Clark, B. E., Hapke, B., Pieters, C., & Britt, D. 2002, in *Asteroid Space Weathering and Regolith Evolution*, Asteroids III, ed. W. F. Bottke et al. (Tucson, AZ: Univ. Arizona Press), 585

* Delbo, M., Libourel, G., Wilkerson, J., et al. 2014, *Natur*, **508**, 233

* DeMeo, F. E., Binzel, R. P., Slivan, S. M., & Bus, S. J. 2009, *Icar*, **202**, 160

* Hanuš, J., Pejcha, O., Shappee, B. J., et al. 2021, *A&A*, **654**, A48

* Hasegawa, S., Hiroi, T., Ohtsuka, K., et al. 2019, *PASJ*, **71**, 103

* Hasegawa, S., Kuroda, D., Yanagisawa, K., & Usui, F. 2017, *PASJ*, **69**, 99

* Hasegawa, S., Marsset, M., DeMeo, F. E., et al. 2021, *ApJL*, **916**, L6

* Hiroi, T., Sasaki, S., Misu, T., & Nakamura, T. 2013, *LPSC*, **44**, 1276

* Hiroi, T., Zolensky, M. E., & Pieters, C. M. 2001, *Sci*, **293**, 2234

* Howell, E. S., & Lovell, A. J. 2011, *IAUC*, **9191**, 2

* Hsieh, H. H., Yang, B., & Haghighipour, N. 2012, *ApJ*, **744**, 9

* Ishiguro, M., Hanayama, H., Hasegawa, S., et al. 2011a, *ApJL*, **740**, L11

* Ishiguro, M., Hanayama, H., Hasegawa, S., et al. 2011b, *ApJL*, **741**, L24

* Jehin, E., Manfroid, J., Hutsemekers, D., Gillon, M., & Magain, P. 2011, *CBET*, **2632**, 2

* Jewitt, D. 2012, *AJ*, **143**, 66

* Jewitt, D., Weaver, H., Mutchler, M., Larson, S., & Agarwal, J. 2011, *ApJL*, **733**, L4

* Johnson, T. V., & Fanale, F. P. 1973, *JGR*, **78**, 8507

* Kim, J., Ishiguro, M., Hanayama, H., et al. 2012, *ApJL*, **746**, L11

* Kim, Y., Ishiguro, M., Michikami, T., & Nakamura, A. M. 2017, *AJ*, **153**, 228

* Küppers, M., O'Rourke, L., Bockelée-Morvan, D., et al. 2014, *Natur*, **505**, 525

* Landolt, A. U. 1983, *AJ*, **88**, 439

* Lantz, C., Binzel, R. P., & DeMeo, F. E. 2018, *Icar*, **302**, 10

* Lantz, C., Brunetto, R., Barucci, M. A., et al. 2017, *Icar*, **285**, 43

* Larson, S., Kowalski, R., Sato, H., Yoshimoto, K., & Nakano, S. 2010, *CBET*, **2583**, 1

* Levison, H. F., Bottke, W. F., Gounelle, M., et al. 2009, *Natur*, **460**, 364

* Licandro, J., de León, J., Kelley, M. S., et al. 2011, in *EPSC-DPS Joint Meeting, Multi-wavelength study of activated asteroid* (596) Scheila, 2011, 1109

* Marsset, M., DeMeo, F., Sonka, A., et al. 2019, *ApJL*, **882**, L2

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

* Masiero, J. R., Mainzer, A. K., Grav, T., et al. 2011, *ApJ*, **741**, 68

* Moreno, F., Licandro, J., Ortiz, J. L., et al. 2011, *ApJ*, **738**, 130

* Morota, T., Sugita, S., Cho, Y., et al. 2020, *Sci*, **368**, 654

* Purdum, J. N., Lin, Z.-Y., Bolin, B. T., et al. 2021, *ApJL*, **911**, L35

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

* Rivkin, A. S., Binzel, R. P., & Bus, S. J. 2005, *Icar*, **175**, 175

* Schultz, P. H., & Gault, D. E. 1985, *JGR*, **90**, 3701

* Shevchenko, V. G., Belskaya, I. N., Muinonen, K., et al. 2016, *P&SS*, **123**, 101

* Shuvalov, V. 2011, *M&PS*, **46**, 1713

* Singer, K. N., Jolliff, B. L., & McKinnon, W. B. 2020, *JGRE*, **125**, e06313

* Stoeffler, D., Gault, D. E., Wedekind, J., & Polkowski, G. 1975, *JGR*, **80**, 4062

* Takir, D., & Emery, J. P. 2012, *Icar*, **219**, 641

* Tedesco, E. F., Noah, P. V., Noah, M., & Price, S. D. 2002, *AJ*, **123**, 1056

* Tody, D. 1993, in *ASP Conf. Ser.* 52, Astronomical Data Analysis Software and Systems II, ed. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes (San Francisco, CA: ASP), 173

* Trigo-Rodriguez, J. M., García-Hernández, D. A., Montañés-Rodríguez, P., et al. 2011, in *EPSC-DPS Joint Meeting 2011, Multiband photometry of asteroid* 596 Scheila after its December 2011 outburst, 2011, 264

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

* Usui, F., Kuroda, D., Müller, T. G., et al. 2011, *PASJ*, **63**, 1117

* Vernazza, P., Castillo-Rogez, J., Beck, P., et al. 2017, *AJ*, **153**, 72

* Vernazza, P., Ferrais, M., Jorda, L., et al. 2021, *A&A*, **654**, A56

* Vernazza, P., Fulvio, D., Brunetto, R., et al. 2013, *Icar*, **225**, 517

* Vernazza, P., Marsset, M., Beck, P., et al. 2015, *ApJ*, **806**, 204

* Vernazza, P., Marsset, M., Beck, P., et al. 2016, *AJ*, **152**, 54

* Yang, B., & Hsieh, H. 2011, *ApJL*, **737**, L39

<footnote n="17">17 http://www.planetary.brown.edu/relab/</footnote>