Source author record

Mikael Granvik

Mikael Granvik appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2026arXiv

Instantaneous thermally-driven erosion can explain dearth of dark near-Sun asteroids

Recent models of the near-Earth asteroid population show that asteroids must be super-catastrophically destroyed when they evolve to orbits with perihelion passages well inside of Mercury's orbit. The heliocentric distances at which the disruptions typically occur are tens of solar radii, which is too far from the Sun for asteroids to be destroyed by sublimation and tidal disruption. The typical disruption distance also appears to be larger for darker asteroids. Here, by carrying out irradiance experiments in vacuum that replicate the conditions in the near-Sun environment, we show that CI meteorite simulants are destroyed within minutes when exposed to the level of solar irradiance encountered at heliocentric distances of up to about 0.2 au. Our results provide an explanation for the scarcity of dark, carbonaceous asteroids with perihelion distances less than 0.2 au, and for the observed mass-loss rate of the asteroid-like object 322P/SOHO~1 assuming its composition is similar to CI carbonaceous chondrites.

preprint2023arXiv

(523599) 2003 RM: The asteroid that wanted to be a comet

We report a statistically significant detection of nongravitational acceleration on the sub-kilometer near-Earth asteroid (523599) 2003 RM. Due to its orbit, 2003 RM experiences favorable observing apparitions every 5 years. Thus, since its discovery, 2003 RM has been extensively tracked with ground-based optical facilities in 2003, 2008, 2013, and 2018. We find that the observed plane-of-sky positions cannot be explained with a purely gravity-driven trajectory. Including a transverse nongravitational acceleration allows us to match all observational data, but its magnitude is inconsistent with perturbations typical of asteroids such as the Yarkovsky effect or solar radiation pressure. After ruling out that the orbital deviations are due to a close approach or collision with another asteroid, we hypothesize that this anomalous acceleration is caused by unseen cometary outgassing. A detailed search for evidence of cometary activity with archival and deep observations from Pan-STARRS and the VLT does not reveal any detectable dust production. However, the best-fitting H$_2$O sublimation model allows for brightening due to activity consistent with the scatter of the data. We estimate the production rate required for H$_2$O outgassing to power the acceleration, and find that, assuming a diameter of 300 m, 2003 RM would require Q(H$_2$O)$\sim10^{23}$ molec s$^{-1}$ at perihelion. We investigate the recent dynamical history of 2003 RM and find that the object most likely originated in the mid-to-outer main belt ($\sim86\%$) as opposed to from the Jupiter-family comet region ($\sim11\%$). Further observations, especially in the infrared, could shed light on the nature of this anomalous acceleration.

preprint2023arXiv

Tidal disruption of near-Earth asteroids during close encounters with terrestrial planets

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

preprint2022arXiv

(3200) Phaethon Polarimetry in the Negative Branch: New Evidence for the Anhydrous Nature of the DESTINY+ Target Asteroid

We report on the first polarimetric study of (3200) Phaethon, the target of JAXA's DESTINY$^+$ mission, in the negative branch to ensure its anhydrous nature and to derive an accurate geometric albedo. We conducted observations at low phase angles (Sun-target-observer angle, alpha = 8.8-32.4 deg) from 2021 October to 2022 January and found that Phaethon has a minimum polarization degree $P_{min}$ = -1.3 +- 0.1 %, a polarimetric slope h = 0.22 +- 0.02 % deg$^{-1}$, and an inversion angle alpha$_0$ = 19.9 +- 0.3 deg. The derived geometric albedo is $p_V$ = 0.11 (in the range of 0.08-0.13). These polarimetric properties are consistent with anhydrous chondrites, and contradict hydrous chondrites and typical cometary nuclei.

preprint2022arXiv

Evidence of surface heterogeneity on active asteroid (3200) Phaethon

Thermal infrared emission and thermophysical modeling techniques are powerful tools in deciphering the surface properties of asteroids. The near-Earth asteroid (3200) Phaethon is an active asteroid with a very small perihelion distance and is likely the source of the Geminid meteor shower. We estimate and interpret the thermal inertia of this extraordinary asteroid using observations that span ten distinct sightings. The variation in thermal inertia over these sightings is inconsistent with the expected temperature-dependent thermal inertia theorized from radiative heat transfer within the regolith. Thus, we test whether the variation in thermal inertia can be explained by modeling a regolith layer over bedrock and two spatially heterogeneous scenarios. We find that the model in which Phaethon's north and south hemispheres have distinctly different thermophysical properties can sufficiently explain the thermal-inertias determined herein. In particular, we find that a boundary located between latitudes -30 deg and +10 deg separates fine-grained southern latitudes from a northern hemisphere that is dominated by coarse-grained regolith and/or a high coverage of porous boulders. We discuss the implications related to Phaethon's activity and potential association with 2005 UD.

preprint2022arXiv

Physical Characterization of 2015 JD1: A Possibly Inhomogeneous Near-Earth Asteroid

The surfaces of airless bodies such as asteroids are exposed to many phenomena that can alter their physical properties. Bennu, the target of the OSIRIS-REx mission, has demonstrated how complex the surface of a small body can be. In 2019 November, the potentially hazardous asteroid 2015 JD1 experienced a close approach of 0.0331 au from the Earth. We present results of the physical characterization of 2015 JD1 based on ground-based radar, spectroscopy, and photometric observations acquired during 2019 November. Radar polarimetry measurements from the Arecibo Observatory indicate a morphologically complex surface. The delay-Doppler images reveal a contact-binary asteroid with an estimated visible extent of ~150 m. Our observations suggest that 2015 JD1 is an E-type asteroid with a surface composition similar to aubrites, a class of differentiated enstatite meteorites. The dynamical properties of 2015 JD1 suggest it came from the $ν_6$ resonance with Jupiter, and spectral comparison with major E-type bodies suggest that it may have been derived from a parental body similar to the progenitor of the E-type (64) Angelina. Significantly, we find rotational spectral variation across the surface of 2015 JD1 from red to blue spectral slope. Our compositional analysis suggests that the spectral slope variation could be due to the lack of iron and sulfides in one area of the 2015 JD1 and/or differences in grain sizes.

preprint2022arXiv

SHINeS: Space and High-Irradiance Near-Sun Simulator

We present SHINeS, a space simulator which can be used to replicate the thermal environment in the immediate neighborhood of the Sun down to a heliocentric distance r~0.06 au. The system consists of three main parts: the solar simulator which was designed and constructed in-house, a vacuum chamber, and the probing and recording equipment needed to monitor the experimental procedures. Our motivation for building this experimental system was to study the effect of intense solar radiation on the surfaces of asteroids when their perihelion distances become smaller than the semi-major axis of the orbit of Mercury. Comparisons between observational data and recent orbit and size-frequency models of the population of near-Earth asteroids suggest that asteroids are super-catastrophically destroyed when they approach the Sun. Whereas the current models are agnostic about the disruption mechanism, SHINeS was developed to study the mechanism or mechanisms responsible. The system can, however, be used for other applications that need to study the effects of high solar radiation on other natural or artificial objects.

preprint2022arXiv

The debiased compositional distribution of MITHNEOS: Global match between the near-Earth and main-belt asteroid populations and excess of D-type Near-Earth Objects

We report 491 new near-infrared spectroscopic measurements of 420 near-Earth objects (NEOs) collected on the NASA InfraRed Telescope Facility (IRTF) as part of the MIT-Hawaii NEO Spectroscopic Survey (MITHNEOS). These measurements were combined with previously published data (Binzel et al. 2019) and bias-corrected to derive the intrinsic compositional distribution of the overall NEO population, as well as of subpopulations coming from various escape routes (ERs) in the asteroid belt and beyond. The resulting distributions reflect well the overall compositional gradient of the asteroid belt, with decreasing fractions of silicate-rich (S- and Q-type) bodies and increasing fractions of carbonaceous (B-, C-, D- and P-type) bodies as a function of increasing ER distance from the Sun. The close compositional match between NEOs and their predicted source populations validates dynamical models used to identify ERs and argues against any strong composition change with size in the asteroid belt between ~5 km down to ~100 m. A notable exception comes from the over-abundance of D-type NEOs from the 5:2J and, to a lesser extend, the 3:1J and nu6 ERs, hinting at the presence of a large population of small D-type asteroids in the main belt. Alternatively, this excess may indicate preferential spectral evolution from D-type surfaces to C- and P-types as a consequence of space weathering, or to the fact that D-type objects fragment more often than other spectral types in the NEO space. No further evidence for the existence of collisional families in the main belt, below the detection limit of current main-belt surveys, was found in this work.

preprint2021arXiv

Masses, bulk densities, and macroporosities of asteroids (15) Eunomia, (29) Amphitrite, (52) Europa, and (445) Edna based on Gaia astrometry

Gaia Data Release 2 (DR2) includes milliarcsecond-accuracy astrometry for 14,099 asteroids. We explore the practical impact of this data for the purpose of asteroid mass and orbit estimation by estimating the masses individually for four large asteroids. We use various combinations of Gaia astrometry and/or Earth-based astrometry so as to determine the impact of Gaia on the estimated masses. By utilizing published information about estimated volumes and meteorite analogs, we also derive estimates for bulk densities and macroporosities. We apply a Markov chain Monte Carlo (MCMC) algorithm for asteroid mass and orbit estimation by modeling asteroid-asteroid close encounters to four separate large asteroids in an attempt to estimate their masses based on multiple simultaneously studied close encounters with multiple test asteroids. In order to validate our algorithm and data treatment, we apply the MCMC algorithm to pure orbit determination for the main-belt asteroid (367) Amicitia and compare the residuals to previously published ones. In addition, we attempt to estimate a mass for (445) Edna with Gaia astrometry alone based on its close encounter with (1764) Cogshall. In the case of the orbit of (367) Amicitia, we find a solution that improves on the previously published solution. The study of (445) Edna reveals that mass estimation with DR2 astrometry alone is unfeasible and that it must be combined with astrometry from other sources to achieve meaningful results. We find that a combination of DR2 and Earth-based astrometry results in dramatically reduced uncertainties and, by extension, significantly improved results in comparison to those computed based on less accurate Earth-based astrometry alone.

preprint2021arXiv

The Sariçiçek howardite fall in Turkey: Source crater of HED meteorites on Vesta and impact risk of Vestoids

The Sariçiçek howardite meteorite shower consisting of 343 documented stones occurred on 2 September 2015 in Turkey and is the first documented howardite fall. Cosmogenic isotopes show that Sariçiçek experienced a complex cosmic ray exposure history, exposed during ~12-14 Ma in a regolith near the surface of a parent asteroid, and that an ca.1 m sized meteoroid was launched by an impact 22 +/- 2 Ma ago to Earth (as did one third of all HED meteorites). SIMS dating of zircon and baddeleyite yielded 4550.4 +/- 2.5 Ma and 4553 +/- 8.8 Ma crystallization ages for the basaltic magma clasts. The apatite U-Pb age of 4525 +/- 17 Ma, K-Ar age of ~3.9 Ga, and the U,Th-He ages of 1.8 +/- 0.7 and 2.6 +/- 0.3 Ga are interpreted to represent thermal metamorphic and impact-related resetting ages, respectively. Petrographic, geochemical and O-, Cr- and Ti- isotopic studies confirm that Sariçiçek belongs to the normal clan of HED meteorites. Petrographic observations and analysis of organic material indicate a small portion of carbonaceous chondrite material in the Sariçiçek regolith and organic contamination of the meteorite after a few days on soil. Video observations of the fall show an atmospheric entry at 17.3 +/- 0.8 kms-1 from NW, fragmentations at 37, 33, 31 and 27 km altitude, and provide a pre-atmospheric orbit that is the first dynamical link between the normal HED meteorite clan and the inner Main Belt. Spectral data indicate the similarity of Sariçiçek with the Vesta asteroid family spectra, a group of asteroids stretching to delivery resonances, which includes (4) Vesta. Dynamical modeling of meteoroid delivery to Earth shows that the disruption of a ca.1 km sized Vesta family asteroid or a ~10 km sized impact crater on Vesta is required to provide sufficient meteoroids <4 m in size to account for the influx of meteorites from this HED clan.

preprint2020arXiv

Debiased albedo distribution for Near Earth Objects

We extend the most recent orbital and absolute magnitude Near Earth Object (NEO) model (Granvik et al., 2018) to provide a statistical description of NEO geometric albedos. Our model is calibrated on NEOWISE albedo data for the NEO population and reproduces these data very well once a simple model for the NEOWISE observational biases is applied. The results are consistent with previous estimates. There are about 1,000 NEOs with diameter D>1km and the mean albedo to convert absolute magnitude into diameter is 0.147. We do not find any statistically significant evidence that the albedo distribution of NEOs depends on NEO size. Instead, we find evidence that the disruption of NEOs at small perihelion distances found in Granvik et al. (2016) occurs preferentially for dark NEOs. The interval between km-sized bodies striking the Earth should occur on average once every 750,000 years. Low and high albedo NEOs are well mixed in orbital space, but a trend remains with higher albedo objects being at smaller semimajor axes and lower albedo objects more likely found at larger semimajor axes.

preprint2016arXiv

OSSOS: IV. Discovery of a dwarf planet candidate in the 9:2 resonance with Neptune

We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has a diameter of $\sim670$ km. Based on astrometric measurements from OSSOS and Pan-STARRS1, we find that 2015 RR$_{245}$ is securely trapped on ten-Myr timescales in the 9:2 mean-motion resonance with Neptune. It is the first TNO identified in this resonance. On hundred-Myr timescales, particles in 2015 RR$_{245}$-like orbits depart and sometimes return to the resonance, indicating that 2015 RR$_{245}$ likely forms part of the long-lived metastable population of distant TNOs that drift between resonance sticking and actively scattering via gravitational encounters with Neptune. The discovery of a 9:2 TNO stresses the role of resonances in the long-term evolution of objects in the scattering disk, and reinforces the view that distant resonances are heavily populated in the current Solar System. This object further motivates detailed modelling of the transient sticking population.

preprint2016arXiv

Realistic Detectability of Close Interstellar Comets

During the planet formation process, billions of comets are created and ejected into interstellar space. The detection and characterization of such interstellar comets (also known as extra-solar planetesimals or extra-solar comets) would give us in situ information about the efficiency and properties of planet formation throughout the galaxy. However, no interstellar comets have ever been detected, despite the fact that their hyperbolic orbits would make them readily identifiable as unrelated to the solar system. Moro-Martín et al. 2009 have made a detailed and reasonable estimate of the properties of the interstellar comet population. We extend their estimates of detectability with a numerical model that allows us to consider "close" interstellar comets, e.g., those that come within the orbit of Jupiter. We include several constraints on a "detectable" object that allow for realistic estimates of the frequency of detections expected from the Large Synoptic Survey Telescope (LSST) and other surveys. The influence of several of the assumed model parameters on the frequency of detections is explored in detail. Based on the expectation from Moro-Martín et al. 2009, we expect that LSST will detect 0.001-10 interstellar comets during its nominal 10-year lifetime, with most of the uncertainty from the unknown number density of small (nuclei of $\sim$0.1-1 km) interstellar comets. Using simulated LSST-like astrometric data, we study the problem of orbit determination for these bodies, finding that LSST could identify their orbits as hyperbolic and determine an ephemeris sufficiently accurate for follow-up in about 4-7 days. We give the hyperbolic orbital parameters of the most detectable interstellar comets. Taking the results into consideration, we give recommendations to future searches for interstellar comets.

preprint2016arXiv

The Outer Solar System Origins Survey: I. Design and First-Quarter Discoveries

We report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis uncertainty $<0.1\%$. We achieve this precision in just two oppositions, as compared to the normal 3-5 oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper belt, and infer the existence of an extension of the "stirred" cold classical Kuiper belt to at least several AU beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. (2011) remains a plausible representation of the Kuiper belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the Solar System.

preprint2015arXiv

Absolute magnitudes and slope parameters for 250,000 asteroids observed by Pan-STARRS PS1 - preliminary results

We present the results of a Monte Carlo technique to calculate the absolute magnitudes (H) and slope parameters (G) of about 240000 asteroids observed by the Pan-STARRS1 telescope during the first 15 months of its 3-year all-sky survey mission. The system's exquisite photometry with photometric errors < 0.04 mag, and well-defined filter and photometric system, allowed us to derive accurate H and G even with a limited number of observations and restricted range in phase angles. Our Monte Carlo method simulates each asteroid's rotation period, amplitude and color to derive the most-likely H and G, but its major advantage is in estimating realistic statistical+systematic uncertainties and errors on each parameter. The method was tested by comparison with the well-established and accurate results for about 500 asteroids provided by Pravec et al. (2012) and then applied to determining H and G for the Pan-STARRS1 asteroids using both the Muinonen et al. (2010) and Bowell et al. (1989) phase functions. Our results confirm the bias in MPC photometry discovered by (Juric et al. 2002).

preprint2015arXiv

Surveys, Astrometric Follow-up & Population Statistics

Asteroid surveys are the backbone of asteroid science, and with this in mind we begin with a broad review of the impact of asteroid surveys on our field. We then provide a brief history of asteroid discoveries so as to place contemporary and future surveys in perspective. Surveys in the United States have discovered the vast majority of the asteroids and this dominance has been consolidated since the publication of Asteroids III. Our descriptions of the asteroid surveys that have been operational since that time are focussed upon those that have contributed the vast majority of asteroid observations and discoveries. We also provide some insight into upcoming next-generation surveys that are sure to alter our understanding of the small bodies in the inner solar system and provide evidence to untangle their complicated dynamical and physical histories. The Minor Planet Center, the nerve center of the asteroid discovery effort, has improved its operations significantly in the past decade so that it can manage the increasing discovery rate, and ensure that it is well-placed to handle the data rates expected in the next decade. We also consider the difficulties associated with astrometric follow-up of newly identified objects. It seems clear that both of these efforts must operate in new modes in order to keep pace with expected discovery rates of next-generation ground- and space-based surveys.

preprint2015arXiv

The daily processing of asteroid observations by Gaia

The Gaia mission started its regular observing program in the summer of 2014, and since then it is regularly obtaining observations of asteroids. This paper draws the outline of the data processing for Solar System objects, and in particular on the daily "short-term" processing, from the on-board data acquisition to the ground-based processing. We illustrate the tools developed to compute predictions of asteroid observations, we discuss the procedures implemented by the daily processing, and we illustrate some tests and validations of the processing of the asteroid observations. Our findings are overall consistent with the expectations concerning the performances of Gaia and the effectiveness of the developed software for data reduction.

preprint2014arXiv

Detecting Earth's Temporarily-Captured Natural Satellites - Minimoons

We present a study on the discoverability of temporarily captured orbiters (TCOs) by present day or near-term anticipated ground-based and space-based facilities. TCOs (Granvik et al. 2012) are potential targets for spacecraft ren- dezvous or human exploration (Chyba et al. 2014) and provide an opportunity to study the population of the smallest asteroids in the solar system. We find that present day ground-based optical surveys such as Pan-STARRS and ATLAS can discover the largest TCOs over years of operation. A targeted survey conducted with the Subaru telescope can discover TCOs in the 0.5 m to 1.0 m diameter size range in about 5 nights of observing. Furthermore, we discuss the application of space-based infrared surveys, such as NEOWISE, and ground-based meteor detection systems such as CAMS, CAMO and ASGARD in discovering TCOs. These systems can detect TCOs but at a uninteresting rate. Finally, we discuss the application of bi-static radar at Arecibo and Green Bank to discover TCOs. Our radar simulations are strongly dependent on the rotation rate distribution of the smallest asteroids but with an optimistic distribution we find that these systems have > 80% chance of detecting a > 10 cm diameter TCO in about 40 h of operation.

preprint2014arXiv

Observational Constraints on the Catastrophic Disruption Rate of Small Main Belt Asteroids

We have calculated 90% confidence limits on the steady-state rate of catastrophic disruptions of main belt asteroids in terms of the absolute magnitude at which one catastrophic disruption occurs per year (HCL) as a function of the post-disruption increase in brightness (delta m) and subsequent brightness decay rate (tau). The confidence limits were calculated using the brightest unknown main belt asteroid (V = 18.5) detected with the Pan-STARRS1 (Pan-STARRS1) telescope. We measured the Pan-STARRS1's catastrophic disruption detection efficiency over a 453-day interval using the Pan-STARRS moving object processing system (MOPS) and a simple model for the catastrophic disruption event's photometric behavior in a small aperture centered on the catastrophic disruption event. Our simplistic catastrophic disruption model suggests that delta m = 20 mag and 0.01 mag d-1 < tau < 0.1 mag d-1 which would imply that H0 = 28 -- strongly inconsistent with H0,B2005 = 23.26 +/- 0.02 predicted by Bottke et al. (2005) using purely collisional models. We postulate that the solution to the discrepancy is that > 99% of main belt catastrophic disruptions in the size range to which this study was sensitive (100 m) are not impact-generated, but are instead due to fainter rotational breakups, of which the recent discoveries of disrupted asteroids P/2013 P5 and P/2013 R3 are probable examples. We estimate that current and upcoming asteroid surveys may discover up to 10 catastrophic disruptions/year brighter than V = 18.5.

preprint2014arXiv

Properties and evolution of NEO families created by tidal disruption at Earth

We have calculated the coherence and detectable lifetimes of synthetic near-Earth object (NEO) families created by catastrophic disruption of a progenitor as it suffers a very close Earth approach. The closest or slowest approaches yield the most violent `s-class' disruption events. We found that the average slope of the absolute magnitude (H) distribution, $N(H)\propto10^{(0.55\pm0.04)\,H}$, for the fragments in the s-class families is steeper than the slope of the NEO population \citep{mainzer2011} in the same size range. The families remain coherent as statistically significant clusters of orbits within the NEO population for an average of $\barτ_c = (14.7\pm0.6)\times10^3$ years after disruption. The s-class families are detectable with the techniques developed by \citet{fu2005} and \citet{Schunova2012} for an average duration ($\barτ_{det}$) ranging from about 2,000 to about 12,000 years for progenitors in the absolute magnitude ($H_p$) range from 20 to 13 corresponding to diameters in the range from about 0.5 to 10$\km$ respectively. The short detectability lifetime explains why zero NEO families have been discovered to-date. Nonetheless, every tidal disruption event of a progenitor with D$>0.5\km$ is capable of producing several million fragments in the $1\meter$ to $10\meter$ diameter range that can contribute to temporary local density enhancements of small NEOs in Earth's vicinity. We expect that there are about 1,200 objects in the steady state NEO population in this size range due to tidal disruption assuming that one $1\km$ diameter NEO tidally disrupts at Earth every 2,500 years. These objects may be suitable targets for asteroid retrieval missions due to their Earth-like orbits with corresponding low $v_{\infty}$. The fragments from the tidal disruptions at Earth have $\sim5\times$ the collision probability with Earth compared to the background NEO population.

preprint2013arXiv

The Pan-STARRS Moving Object Processing System

We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development.

preprint2012arXiv

Discovery of Main-Belt Comet P/2006 VW139 by Pan-STARRS1

Main belt asteroid (300163) 2006 VW139 (later designated P/2006 VW139) was discovered to exhibit comet-like activity by the Pan-STARRS1 survey telescope using automated point-spread-function analyses performed by PS1's Moving Object Processing System. Deep follow-up observations show both a short (\sim 10") antisolar dust tail and a longer (\sim 60") dust trail aligned with the object's orbit plane, similar to the morphology observed for another main-belt comet, P/2010 R2 (La Sagra), and other well-established comets, implying the action of a long-lived, sublimation-driven emission event. Photometry showing the brightness of the near-nucleus coma remaining constant over \sim 30 days provides further evidence for this object's cometary nature, suggesting it is in fact a main-belt comet, and not a disrupted asteroid. A spectroscopic search for CN emission was unsuccessful, though we find an upper limit CN production rate of Q_CN < 1.3x10^24 mol/s, from which we infer a water production rate of Q_H2O < 10^26 mol/s. We also find an approximately linear optical spectral slope of 7.2%/1000A, similar to other cometary dust comae. Numerical simulations indicate that P/2006 VW139 is dynamically stable for > 100 Myr, while a search for a potential asteroid family around the object reveals a cluster of 24 asteroids within a cutoff distance of 68 m/s. At 70 m/s, this cluster merges with the Themis family, suggesting that it could be similar to the Beagle family to which another main-belt comet, 133P/Elst-Pizarro, belongs.

preprint2012arXiv

Modeling collision probability for Earth-impactor 2008 TC3

We study the evolution of the Earth collision probability of asteroid 2008 TC3 using a short observational arc and small numbers of observations. To assess impact probability, we use techniques that rely on the orbital-element probability density function characterized using both Markov-chain Monte-Carlo orbital ranging and Monte-Carlo ranging. First, we evaluate the orbital uncertainties for the object from the night of discovery onwards and examine the collapse of the orbital-element distributions in time. Second, we examine the sensitivity of the results to the assumed astrometric noise. Each of the orbits obtained from the MCMC ranging method is propagated into the future (within chosen time bounds of the expected impact), and the collision probability is calculated as a weighted fraction of the orbits leading to a collision from the Earth. We compare the results obtained with both methods.

preprint2012arXiv

Observational and Dynamical Characterization of Main-Belt Comet P/2010 R2 (La Sagra)

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out.

preprint2012arXiv

Searching for the first Near-Earth Object family

We report on our search for genetically related asteroids amongst the near-Earth object (NEO) population - families of NEOs akin to the well known main belt asteroid families. We used the technique proposed by Fu et al. (2005) supplemented with a detailed analysis of the statistical significance of the detected clusters. Their significance was assessed by comparison to identical searches performed on 1,000 'fuzzy-real' NEO orbit distribution models that we developed for this purpose. The family-free 'fuzzy-real' NEO models maintain both the micro and macro distribution of 5 orbital elements (ignoring the mean anomaly). Three clusters were identified that contain four or more NEOs but none of them are statistically significant at \geq 3σ. The most statistically significant cluster at the \sim 2σ level contains 4 objects with H < 20 and all members have long observational arcs and concomitant good orbital elements. Despite the low statistical significance we performed several other tests on the cluster to determine if it is likely a genetic family. The tests included examining the cluster's taxonomy, size-frequency distribution, consistency with a family-forming event during tidal disruption in a close approach to Mars, and whether it is detectable in a proper element cluster search. None of these tests exclude the possibility that the cluster is a family but neither do they confirm the hypothesis. We conclude that we have not identified any NEO families.

preprint2011arXiv

The population of natural Earth satellites

We have for the first time calculated the population characteristics of the Earth's irregular natural satellites (NES) that are temporarily captured from the near-Earth-object (NEO) population. The steady-state NES size-frequency and residence-time distributions were determined under the dynamical influence of all the massive bodies in the solar system (but mainly the Sun, Earth, and Moon) for NEOs of negligible mass. To this end, we compute the NES capture probability from the NEO population as a function of the latter's heliocentric orbital elements and combine those results with the current best estimates for the NEO size-frequency and orbital distribution. At any given time there should be at least one NES of 1-meter diameter orbiting the Earth. The average temporarily-captured orbiter (TCO; an object that makes at least one revolution around the Earth in a co-rotating coordinate system) completes $(2.88\pm0.82)\rev$ around the Earth during a capture event that lasts $(286\pm18)\days$. We find a small preference for capture events starting in either January or July. Our results are consistent with the single known natural TCO, 2006 RH$_{120}$, a few-meter diameter object that was captured for about a year starting in June 2006. We estimate that about 0.1% of all meteors impacting the Earth were TCOs.

preprint2009arXiv

Detection of Earth-impacting asteroids with the next generation all-sky surveys

We have performed a simulation of a next generation sky survey's (Pan-STARRS 1) efficiency for detecting Earth-impacting asteroids. The steady-state sky-plane distribution of the impactors long before impact is concentrated towards small solar elongations (Chesley and Spahr, 2004) but we find that there is interesting and potentially exploitable behavior in the sky-plane distribution in the months leading up to impact. The next generation surveys will find most of the dangerous impactors (>140m diameter) during their decade-long survey missions though there is the potential to miss difficult objects with long synodic periods appearing in the direction of the Sun, as well as objects with long orbital periods that spend much of their time far from the Sun and Earth. A space-based platform that can observe close to the Sun may be needed to identify many of the potential impactors that spend much of their time interior to the Earth's orbit. The next generation surveys have a good chance of imaging a bolide like 2008TC3 before it enters the atmosphere but the difficulty will lie in obtaining enough images in advance of impact to allow an accurate pre-impact orbit to be computed.