BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Robert Jedicke

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.IMastro-phastro-ph.GAastro-ph.HEastro-ph.SR

Catalog footprint

What is connected

29works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

29 published item(s)

preprint2024arXiv

Linking tracklets over the years in large datasets

We present a new procedure to identify observations of known objects in large data sets of unlinked detections. It begins with a Keplerian integrals method that allows us to link two tracklets, computing preliminary orbits, even when the tracklets are separated in time by a few years. In the second step, we represent the results in a `graph' where the tracklets are the nodes and the preliminary orbits are the edges. Then, acceptable `3-cycles' are identified and a least squares orbit is computed for each of them. Finally, we construct sequences of $n \geq 4$ tracklets by searching through the orbits of nearby 3-cycles and attempting to attribute the remaining tracklets. We calculate the technique's efficiency at identifying unknown objects using real detections that attempt to mimic key parameters of the Minor Planet Center's Isolated Tracklet File (ITF) and then apply the procedure to the ITF to identify tens of thousands of new objects.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Numerical behavior of the Keplerian Integral methods for initial orbit determination

We investigate the behaviour of two recent methods for the computation of preliminary orbits. These methods are based on the conservation laws of Kepler's problem, and enable the linkage of very short arcs of optical observations even when they are separated in time by a few years. Our analysis is performed using both synthetic and real data of 822 main belt asteroids. The differences between computed and true orbital elements have been analysed for the true linkages, as well as the occurrence of alternative solutions. Some metrics have been introduced to quantify the results, with the aim of discarding as many of the false linkages as possible and keeping the vast majority of true ones. These numerical experiments provide thresholds for the metrics which take advantage of the knowledge of the \emph{ground truth}: the values of these thresholds can be used in normal operation mode, when we do not know the correct values of the orbital elements and whether the linkages are true or false.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

R2-D2: Roman and Rubin -- From Data to Discovery

The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Discovery of A New Retrograde Trans-Neptunian Object: Hint of A Common Orbital Plane for Low Semi-Major Axis, High Inclination TNOs and Centaurs

Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$σ$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The size-frequency distribution of H>13 NEOs and ARM targets detected by Pan-STARRS1

We determine the absolute magnitude (H) distribution (or size-frequency distribution, SFD; $N(H) \propto 10^{αH}$ where $α$ is the slope of the distribution) for near-Earth objects (NEO) with $13<H<30$ and Asteroid Retrieval Mission (ARM) targets with $27<H<31$ that were detected by the 1\st\ telescope of the Panoramic Survey Telescope and Rapid Response System - Pan-STARRS1 (e.g. Kaiser et al. 2002, Kaiser 2004, Hodapp et al. 2004). The NEO and ARM target detection efficiencies were calculated using the Greenstreet et al. (2012) NEO orbit distribution. The debiased Pan-STARRS1 NEO absolute magnitude distribution is more complex than a single slope power law - it shows two transitions - at H$\sim$16 from steep to shallow slope, and in the $21<H<23$ interval from a shallow to steep slope, which is consistent with other recent works (e.g. Mainzer et al. 2011c, Brown et al. 2013, Harris and D`Abramo 2015). We fit $α= 0.48\pm0.02$ for NEOs with $13<H<16$, $α= 0.33\pm0.01$ for NEOs with $16<H<22$, and $α= 0.62\pm0.03$ for the smaller objects with $H>22$. There is also another change in slope from steep to shallow around H=27. The three ARM target candidates detected by Pan-STARRS1 in one year of surveying have a corrected SFD with slope $α= 0.40^{+0.33}_{-0.45}$. We also show that the window for follow up observations of small (H$\gtrsim$22) NEOs with the NASA IRTF telescope and Arecibo and Goldstone radars are extremely short - on order of days, and procedures for fast response must be implemented in order to measure physical characteristics of small Earth-approaching objects. CFHT's MegaCam and Pan-STARRS1 have longer observing windows and are capable of following-up more NEOs due to their deeper limiting magnitudes and wider fields of view.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Search for the Return of Activity in Active Asteroid 176P/LINEAR

We present the results of a search for the reactivation of active asteroid 176P/LINEAR during its 2011 perihelion passage using deep optical observations obtained before, during, and after that perihelion passage. Deep composite images of 176P constructed from data obtained between June 2011 and December 2011 show no visible signs of activity, while photometric measurements of the object during this period also show no significant brightness enhancements similar to that observed for 176P between November 2005 and December 2005 when it was previously observed to be active. An azimuthal search for dust emission likewise reveals no evidence for directed emission (i.e., a tail, as was previously observed for 176P), while a one-dimensional surface brightness profile analysis shows no indication of a spherically symmetric coma at any time in 2011. We conclude that 176P did not in fact exhibit activity in 2011, at least not on the level on which it exhibited activity in 2005, and suggest that this could be due to the devolatization or mantling of the active site responsible for its activity in 2005.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The effect of parallax and cadence on asteroid impact probabilities and warning times

We study the time evolution of the impact probability for synthetic but realistic impacting and close approaching asteroids detected in a simulated all-sky survey. We use the impact probability to calculate the impact warning time as the time interval between when an object reaches a Palermo Scale value of -2 and when it impacts Earth. A simple argument shows that warning time ~ D^x with the exponent in the range [1.0,1.5] and our derived value was $x=1.3\pm0.1$. The low-precision astrometry from the single simulated all-sky survey could require many days or weeks to establish an imminent impact for asteroids >100m diameter that are discovered far from Earth. Most close approaching asteroids are quickly identified as not being impactors but a size-dependent percentage, even for those >50m diameter, have a persistent impact probability of >10^{-6} on the day of closest approach. Thus, a single all-sky survey can be of tremendous value in identifying Earth impacting and close approaching asteroids in advance of their closest approach but it can not solve the problem on its own: high-precision astrometry from other optical or radar systems is necessary to rapidly establish an object as an impactor or close approacher. We show that the parallax afforded by surveying the sky from two sites is only of benefit for a small fraction of the smallest objects detected within a couple days before impact: probably not enough to justify the increased operating costs of a 2-site survey. Finally, the survey cadence within a fixed time span is relatively unimportant to the impact probability calculation. We tested three different reasonable cadences and found that one provided ~10x higher (better) value for the impact probability on the discovery night for the smallest (10m diameter) objects but the consequences on the overall impact probability calculation is negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Main-Belt Comets: The Pan-STARRS1 Perspective

We analyze 760475 observations of 333026 main-belt objects obtained by the Pan-STARRS1 (PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discovered two main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 comet detection procedures currently consist of the comparison of the point spread functions (PSFs) of moving objects to those of reference stars, and the flagging of objects that show anomalously large radial PSF widths. Based on the number of missed discovery opportunities among comets discovered by other observers, we estimate an upper limit comet discovery efficiency rate of ~70% for PS1. Additional analyses that could improve comet discovery yields in future surveys include linear PSF analysis, modeling of trailed stellar PSFs for comparison to trailed moving object PSFs, searches for azimuthally localized activity, comparison of point-source-optimized photometry to extended-source-optimized photometry, searches for photometric excesses in objects with known absolute magnitudes, and crowd-sourcing. Analysis of PS1 survey statistics indicates an expected fraction of 59 MBCs per 10^6 outer main-belt asteroids, and a 95% confidence upper limit of 96 MBCs per 10^6 outer main-belt asteroids. We note that more sensitive future surveys could detect many more MBCs than estimated here. We find an excess of high eccentricities (0.1 < e < 0.3) among all known MBCs relative to the background asteroid population. Theoretical calculations show that, given these eccentricities, the sublimation rate for a typical MBC is orders of magnitude larger at perihelion than at aphelion, providing a plausible physical explanation for the observed behavior of MBCs peaking in observed activity strength near perihelion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Main-Belt Comet P/2012 T1 (PANSTARRS)

We present initial results from observations and numerical analyses aimed at characterizing main-belt comet P/2012 T1 (PANSTARRS). Optical monitoring observations were made between October 2012 and February 2013 using the University of Hawaii 2.2 m telescope, the Keck I telescope, the Baade and Clay Magellan telescopes, Faulkes Telescope South, the Perkins Telescope at Lowell Observatory, and the Southern Astrophysical Research (SOAR) telescope. The object's intrinsic brightness approximately doubles from the time of its discovery in early October until mid-November and then decreases by ~60% between late December and early February, similar to photometric behavior exhibited by several other main-belt comets and unlike that exhibited by disrupted asteroid (596) Scheila. We also used Keck to conduct spectroscopic searches for CN emission as well as absorption at 0.7 microns that could indicate the presence of hydrated minerals, finding an upper limit CN production rate of QCN<1.5x10^23 mol/s, from which we infer a water production rate of QH2O<5x10^25 mol/s, and no evidence of the presence of hydrated minerals. Numerical simulations indicate that P/2012 T1 is largely dynamically stable for >100 Myr and is unlikely to be a recently implanted interloper from the outer solar system, while a search for potential asteroid family associations reveal that it is dynamically linked to the ~155 Myr-old Lixiaohua asteroid family.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

NEOKepler: Discovering Near-Earth Objects Using the Kepler Spacecraft

We propose a new Kepler mission, called NEOKepler, that would survey near Earth's orbit to identify potentially hazardous objects (PHOs). To understand its surveying power, Kepler's large field of view produces an 'etendue' (A*Omega) that is 4.5 times larger than the best survey telescope currently in operation. In this paper, we investigate the feasibility of NEOKepler using a double "fence post" survey pattern that efficiently detects PHOs. In a simulated 12-month survey, we estimate that NEOKepler would detect ~150 new NEOs with absolute magnitudes of less than 21.5, ~50 of which would be new PHOs. This would increase the annual PHO discovery rate by at least 50% and improve upon our goal of discovering 90% of PHOs by the end of 2020. Due to its heliocentric orbit, Kepler would also be sensitive to objects inside Earth's orbit, discovering more objects in its first year than are currently known to exist. Understanding this undersampled sub-population of NEOs will reveal new insights into the actual PHO distribution by further constraining current NEO models. As an alternative science goal, NEOKepler could employ a different observing strategy to discover suitable targets for NASA's Asteroid Redirect Mission.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Outgassing Behavior of C/2012 S1 (ISON) From September 2011 to June 2013

We report photometric observations for comet C/2012 S1 (ISON) obtained during the time period immediately after discovery (r=6.28 AU) until it moved into solar conjunction in mid-2013 June using the UH2.2m, and Gemini North 8-m telescopes on Mauna Kea, the Lowell 1.8m in Flagstaff, the Calar Alto 1.2m telescope in Spain, the VYSOS-5 telescopes on Mauna Loa Hawaii and data from the CARA network. Additional pre-discovery data from the Pan STARRS1 survey extends the light curve back to 2011 September 30 (r=9.4 AU). The images showed a similar tail morphology due to small micron sized particles throughout 2013. Observations at sub-mm wavelengths using the JCMT on 15 nights between 2013 March 9 (r=4.52 AU) and June 16 (r=3.35 AU) were used to search for CO and HCN rotation lines. No gas was detected, with upper limits for CO ranging between (3.5-4.5)E27 molec/s. Combined with published water production rate estimates we have generated ice sublimation models consistent with the photometric light curve. The inbound light curve is likely controlled by sublimation of CO2. At these distances water is not a strong contributor to the outgassing. We also infer that there was a long slow outburst of activity beginning in late 2011 peaking in mid-2013 January (r~5 AU) at which point the activity decreased again through 2013 June. We suggest that this outburst was driven by CO injecting large water ice grains into the coma. Observations as the comet came out of solar conjunction seem to confirm our models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing Accuracy and Precision of Existing Photometry Algorithms on Moving Targets

Previous studies determining which astronomical photometry software is best suited for a particular dataset are usually focused on speed, source classification, and/or meeting a sensitivity requirement. For faint objects in particular, the priority is given to maximizing signal-to-noise. Photometry of moving targets offers additional challenges (i) to aperture photometry because background object contamination varies from image to image, and (ii) to routines that build a PSF model from point sources in the image because trailed field stars do not perfectly represent the PSF of the untrailed target. Here, we present the results of testing several photometry algorithms (tphot, DAOPHOT, DoPHOT, APT, and multiple techniques within Source Extractor and IRAF's PHOT) on data for a faint, slow-moving solar system object with a known light curve. We find that the newly-developed tphot software most accurately and precisely reproduces the object's true light curve, with particular advantages in centroiding, exclusion of contaminants from the target's flux, and fitting flux in the wings of the point-spread function.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Identification of known objects in solar system surveys

The discovery of new objects in modern wide-field asteroid and comet surveys can be enhanced by first identifying observations belonging to known solar system objects. The assignation of new observations to a known object is an attribution problem that occurs when a least squares orbit already exists for the object but a separate fit is not possible to just the set of new observations. In this work we explore the strongly asymmetric attribution problem in which the existing least squares orbit is very well constrained and the new data are sparse. We describe an attribution algorithm that introduces new quality control metrics in the presence of strong biases in the astrometric residuals. The main biases arise from the stellar catalogs used in the reduction of asteroid observations and we show that a simple debiasing with measured regional catalog biases significantly improves the results. We tested the attribution algorithm using data from the PS1 survey that used the 2MASS star catalog for the astrometric reduction. We found small but statistically significant biases in the data of up to 0.1 arcsec that are relevant only when the observations reach the level of accuracy made possible by instruments like PS1. The false attribution rate was measured to be < 1/1,000 with a simple additional condition that can reduce it to zero while the attribution efficiency is consistent with 100%.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Improved Asteroid Astrometry and Photometry with Trail Fitting

Asteroid detections in astronomical images may appear as trails due to a combination of their apparent rate of motion and exposure duration. Nearby asteroids in particular typically have high apparent rates of motion and acceleration. Their recovery, especially on their discovery apparition, depends upon obtaining good astrometry from the trailed detections. We present an analytic function describing a trailed detection under the assumption of a Gaussian point spread function (PSF) and constant rate of motion. We have fit the function to both synthetic and real trailed asteroid detections from the Pan-STARRS1 survey telescope to obtain accurate astrometry and photometry. For short trails our trailing function yields the same astrometric and photometry accuracy as a functionally simpler 2-d Gaussian but the latter underestimates the length of the trail - a parameter that can be important for measuring the object's rate of motion and assessing its cometary activity. For trails longer than about 10 pixels (> 3xPSF) our trail fitting provides 3-times better astrometric accuracy and up to 2 magnitudes improvement in the photometry. The trail fitting algorithm can be implemented at the source detection level for all detections to provide trail length and position angle that can be used to reduce the false tracklet rate.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Asteroid age distributions determined by space weathering and collisional evolution models

We provide evidence of consistency between the dynamical evolution of main belt asteroids and their color evolution due to space weathering. The dynamical age of an asteroid's surface \citep{bib.bot05a,bib.nes05} is the time since its last catastrophic disruption event which is a function of the object's diameter. The age of an S-complex asteroid's surface may also be determined from its color using a space weathering model \citep[e.g.][]{bib.wil10,bib.jed04,bib.wil08,bib.mar06}. We used a sample of 95 S-complex asteroids from SMASS and obtained their absolute magnitudes and $u,g,r,i,z$ filter magnitudes from SDSS. The absolute magnitudes yield a size-derived age distribution. The $u,g,r,i,z$ filter magnitudes lead to the principal component color which yields a color-derived age distribution by inverting our color-age relationship, an enhanced version of the `dual $τ$' space weathering model of \citet{bib.wil10}. We fit the size-age distribution to the enhanced dual $τ$ model and found characteristic weathering and gardening times of $τ_w = 2050 \pm 80$ Myr and $τ_g = 4400^{+700}_{-500}$ Myr respectively. The fit also suggests an initial principal component color of $-0.05 \pm 0.01$ for fresh asteroid surface with a maximum possible change of the probable color due to weathering of $ΔPC = 1.34 \pm 0.04$. Our predicted color of fresh asteroid surface matches the color of fresh ordinary chondritic surface of $PC_1 = 0.17 \pm 0.39$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Thousand Asteroid Light Curve Survey

We present the results of our Thousand Asteroid Light Curve Survey (TALCS) conducted with the Canada-France-Hawaii Telescope in September 2006. Our untargeted survey detected 828 Main Belt asteroids to a limiting magnitude of g'~22.5 corresponding to a diameter range of 0.4 km <= D <= 10 km. Of these, 278 objects had photometry of sufficient quality to perform rotation period fits. We debiased the observations and light curve fitting process to determine the true distribution of rotation periods and light curve amplitudes of Main Belt asteroids. We confirm a previously reported excess in the fraction of fast rotators but find a much larger excess of slow rotating asteroids (~15% of our sample). A few percent of objects in the TALCS size range have large light curve amplitudes of ~1 mag. Fits to the debiased distribution of light curve amplitudes indicate that the distribution of triaxial ellipsoid asteroid shapes is proportional to the square of the axis-ratio, (b/a)^2, and may be bi-modal. Finally, we find six objects with rotation periods that may be less than 2 hours with diameters between 400 m and 1.5 km, well above the break-up limit for a gravitationally-bound aggregate. Our debiased data indicate that this population represents <4% of the Main Belt in the 1-10 km size range.

0 citations0 reviewsNo code linkedOpen access