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.

Earth's temporarily-captured orbiters (TCOs) are a sub-population of near-Earth objects (NEOs). TCOs can provide constraints for NEO population models in the 1--10-metre-diameter range, and they are outstanding targets for in situ exploration of asteroids due to a low requirement on $\Delta v$. So far there has only been a single serendipitous discovery of a TCO. Here we assess in detail the possibility of their discovery with the upcoming Large Synoptic Survey Telescope (LSST), previously identified as the primary facility for such discoveries. We simulated observations of TCOs by combining a synthetic TCO population with an LSST survey simulation. We then assessed the detection rates, detection linking and orbit computation, and sources for confusion. Typical velocities of detectable TCOs will range from 1$^{\circ}$/day to 50$^{\circ}$/day, and typical apparent $V$ magnitudes from 21 to 23. Potentially-hazardous asteroids have observational characteristics similar to TCOs, but the two populations can be distinguished based on their orbits with LSST data alone. We predict that a TCO can be discovered once every year with the baseline moving-object processing system (MOPS). The rate can be increased to one TCO discovery every two months if tools complementary to the baseline MOPS are developed for the specific purpose of discovering these objects.

(modified from published version) Twelve years ago the Catalina Sky Survey discovered Earth's first known natural geocentric object other than the Moon, a few-meter diameter asteroid designated 2006 RH120. Despite significant improvements in ground-based telescope and detector technology in the past decade the asteroid surveys have not discovered another temporarily-captured orbiter (TCO; colloquially known as minimoons). Within a few years the Large Synoptic Survey Telescope (LSST) will either begin to regularly detect TCOs or force a re-analysis of the creation and dynamical evolution of small asteroids in the inner solar system. The first studies of the provenance, properties, and dynamics of Earth's minimoons suggested that there should be a steady state population with about one 1- to 2-meter diameter captured objects at any time. That model was then improved and extended to include the population of temporarily-captured flybys (TCFs), objects that fail to make an entire revolution around Earth while energetically bound to the Earth-Moon system. Several different techniques for discovering TCOs have been considered but their small diameters, proximity, and rapid motion make them challenging targets for existing ground-based optical, meteor, and radar surveys. We expect that if the TCO population is confirmed, and new objects are frequently discovered, they can provide new opportunities for 1) studying the dynamics of the Earth-Moon system, 2) testing models of the production and dynamical evolution of small asteroids from the asteroid belt, 3) rapid and frequent low delta-v missions to multiple minimoons, and 4) evaluating in-situ resource utilization techniques on asteroidal material. Here we review the past decade of minimoon studies in preparation for capitalizing on the scientific and commercial opportunities of TCOs in the first decade of LSST operations.

We introduce a new technique to estimate the comet nuclear size frequency distribution (SFD) that combines a cometary activity model with a survey simulation and apply it to 150 long period comets (LPC) detected by the Pan-STARRS1 near-Earth object survey. The debiased LPC size-frequency distribution is in agreement with previous estimates for large comets with nuclear diameter $>\sim 1$~km but we measure a significant drop in the SFD slope for small objects with diameters $<1$~km and approaching only $100$~m diameter. Large objects have a slope $\alpha_{big} = 0.72 \pm 0.09 (stat.) \pm 0.15 (sys.)$ while small objects behave as $\alpha_{small} = 0.07 \pm 0.03 (stat.) \pm 0.09 (sys.)$ where the SFD is $\propto 10^{\alpha H_N}$ and $H_N$ represents the cometary nuclear absolute magnitude. The total number of LPCs that are $>1$~km diameter and have perihelia $q<10$~au is $0.46 \pm 0.15 \times 10^9$ while there are only $2.4 \pm 0.5 (stat.) \pm 2 (sys.) \times 10^9$ objects with diameters $>100$~m due to the shallow slope of the SFD for diameters $<1$~m. We estimate that the total number of `potentially active' objects with diameters $\ge 1$~km in the Oort cloud, objects that would be defined as LPCs if their perihelia evolved to $<10$~au, is $(1.5\pm1)\times10^{12}$ with a combined mass of $1.3\pm0.9 \, M_{Earth}$. The debiased LPC orbit distribution is broadly in agreement with expectations from contemporary dynamical models but there are discrepancies that could point towards a future ability to disentangle the relative importance of stellar perturbations and galactic tides in producing the LPC population.

We describe the digest2 software package, a fast, short-arc orbit classifier for small Solar System bodies. The digest2 algorithm has been serving the community for more than 13 years. The code provides a score, D2, which represents a pseudo-probability that a tracklet belongs to a given Solar System orbit type. digest2 is primarily used as a classifier for Near-Earth Object (NEO) candidates, to identify those to be prioritized for follow-up observation. We describe the historical development of digest2 and demonstrate its use on real and synthetic data. We find that digest2 can accurately and precisely distinguish NEOs from non- NEOs. At the time of detection, 14% of NEO tracklets and 98.5% of non-NEOs tracklets have D2 below the critical value of D2 = 65. 94% of our simulated NEOs achieved the maximum D2 = 100 and 99.6% of NEOs achieved $D2 \ge 65$ at least once during the simulated 10-year timeframe. We demonstrate that D2 varies as a function of time, rate of motion, magnitude and sky-plane location, and show that NEOs tend to have lower D2 at low Solar elongations close to the ecliptic. We use our findings to recommend future development directions for the digest2 code.

As part of our multi-observatory, multi-filter campaign, we present \rmi color observations of 82 Near-Earth Objects (NEOs) obtained with the RATIR instrument on the 1.5m robotic telescope at the San Pedro Martir's National Observatory in Mexico. Our project is particularly focused on rapid response observations of small ($\lesssim 850$ m) NEOs. The rapid response and the use of spectrophotometry allows us to constrain the taxonomic classification of NEOs with high efficiency. Here we present the methodology of our observations and our result, suggesting that the ratio of C-type to S-type asteroids in a size range of $\sim$30-850m is 1.1, which is in accordance with our previous results. We also find that 10$\%$ of all NEOs in our sample are neither C- nor S-type asteroids

This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3\pi$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema, and how the database tables are related. Examples of queries for a range of science goals are included. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI.

The Pan-STARRS1 Science Consortium have carried out a set of imaging surveys using the 1.4 giga-pixel GPC1 camera on the PS1 telescope. As this camera is composed of many individual electronic readouts, and covers a very large field of view, great care was taken to ensure that the many instrumental effects were corrected to produce the most uniform detector response possible. We present the image detrending steps used as part of the processing of the data contained within the public release of the Pan-STARRS1 Data Release 1 (DR1). In addition to the single image processing, the methods used to transform the 375,573 individual exposures into a common sky-oriented grid are discussed, as well as those used to produce both the image stack and difference combination products.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

The Pan-STARRS Data Processing System is responsible for the steps needed to downloaded, archive, and process all images obtained by the Pan-STARRS telescopes, including real-time detection of transient sources such as supernovae and moving objects including potentially hazardous asteroids. With a nightly data volume of up to 4 terabytes and an archive of over 4 petabytes of raw imagery, Pan-STARRS is solidly in the realm of Big Data astronomy. The full data processing system consists of several subsystems covering the wide range of necessary capabilities. This article describes the Image Processing Pipeline and its connections to both the summit data systems and the outward-facing systems downstream. The latter include the Moving Object Processing System (MOPS) & the public database: the Published Science Products Subsystem (PSPS).

We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3\pi$ Survey. The photometric goals were to reduce the systematic effects introduced by the camera and detectors, and to place all of the observations onto a photometric system with consistent zero points over the entire area surveyed, the ~30,000 square degrees north of $\delta$ = -30 degrees. The astrometric calibration compensates for similar systematic effects so that positions, proper motions, and parallaxes are reliable as well. The Pan-STARRS Data Release 2 (DR2) astrometry is tied to the Gaia DR1 release.

Over 3 billion astronomical objects have been detected in the more than 22 million orthogonal transfer CCD images obtained as part of the Pan-STARRS1 $3\pi$ survey. Over 85 billion instances of those objects have been automatically detected and characterized by the Pan-STARRS Image Processing Pipeline photometry software, psphot. This fast, automatic, and reliable software was developed for the Pan-STARRS project, but is easily adaptable to images from other telescopes. We describe the analysis of the astronomical objects by psphot in general as well as for the specific case of the 3rd processing version used for the first two public releases of the Pan-STARRS $3\pi$ survey data, DR1 & DR2.

Using the first 18 months of the Pan-STARRS 1 survey we have identified 33 candidate high-amplitude objects for follow-up observations and carried out observations of 22 asteroids. 4 of the observed objects were found to have observed amplitude $A_{obs}\geq 1.0$ mag. We find that these high amplitude objects are most simply explained by single rubble pile objects with some density-dependent internal strength, allowing them to resist mass shedding even at their highly elongated shapes. 3 further objects although below the cut-off for 'high-amplitude' had a combination of elongation and rotation period which also may require internal cohesive strength, depending on the density of the body. We find that none of the 'high-amplitude asteroids' identified here require any unusual cohesive strengths to resist rotational fission. 3 asteroids were sufficiently observed to allow for shape and spin pole models to be determined through light curve inversion. 45864 was determined to have retrograde rotation with spin pole axes $\lambda=218\pm 10^{\circ}, \beta=-82\pm 5^{\circ}$ and asteroid 206167 was found to have best fit spin pole axes $\lambda= 57 \pm 5^{\circ}$, $\beta=-67 \pm 5^{\circ}$. An additional object not initially measured with $A_{obs}>1.0$ mag, 49257, was determined to have a shape model which does suggest a high-amplitude object. Its spin pole axes were best fit for values $\lambda=112\pm 6^{\circ}, \beta=6\pm 5^{\circ}$. In the course of this project to date no large super-fast rotators ($P_{rot} < 2.2$ h) have been identified.

The debiased absolute-magnitude and orbit distributions as well as source regions for near-Earth objects (NEOs) provide a fundamental frame of reference for studies of individual NEOs and more complex population-level questions. We present a new four-dimensional model of the NEO population that describes debiased steady-state distributions of semimajor axis, eccentricity, inclination, and absolute magnitude $H$ in the range $17<H<25$. The modeling approach improves upon the methodology originally developed by Bottke et al. (2000; Science 288, 2190-2194) in that it is, for example, based on more realistic orbit distributions and uses source-specific absolute-magnitude distributions that allow for a power-law slope that varies with $H$. We divide the main asteroid belt into six different entrance routes or regions (ER) to the NEO region: the $\nu_6$, 3:1J, 5:2J and 2:1J resonance complexes as well as Hungarias and Phocaeas. In addition we include the Jupiter-family comets as the primary cometary source of NEOs. We calibrate the model against NEO detections by Catalina Sky Surveys' stations 703 and G96 during 2005-2012, and utilize the complementary nature of these two systems to quantify the systematic uncertainties associated to the resulting model. The most important ERs are the $\nu_6$ and 3:1J resonance complexes with JFCs contributing a few percent of NEOs on average. We predict that there are $962^{+52}_{-56}$ ($802^{+48}_{-42}\times10^3$) NEOs with $H<17.75$ ($H<25$) and these numbers are in agreement with the most recent estimates found in the literature (the uncertainty estimates only account for the random component).

The majority of known asteroid diameters are derived from thermal-infrared observations. Diameters are derived using asteroid thermal models that approximate their surface temperature distributions and compare the measured thermal-infrared flux with model-dependent predictions. The most commonly used thermal model is the Near-Earth Asteroid Thermal Model (NEATM), which is usually perceived as superior to other models like the Fast-Rotating Model (FRM). We investigate the applicability of the NEATM and the FRM to thermal-infrared observations of Near-Earth Objects using synthetic asteroids with properties based on the real Near-Earth Asteroid (NEA) population. We find the NEATM to provide more accurate diameters and albedos than the FRM in most cases, with a few exceptions. The modeling results are barely affected by the physical properties of the objects, but we find a large impact of the solar phase angle on the modeling results. We conclude that the NEATM provides statistically more robust diameter estimates for NEAs observed at solar phase angles less than ~65{\deg}, while the FRM provides more robust diameter estimates for solar phase angles greater than ~65{\deg}. We estimate that <5% of all NEA diameters and albedos derived up to date are affected by systematic effects that are of the same order of magnitude as the typical thermal model uncertainties. We provide statistical correction functions for diameters and albedos derived using the NEATM and FRM as a function of solar phase angle.

During the formation and evolution of the Solar System, significant numbers of cometary and asteroidal bodies were ejected into interstellar space$^{1,2}$. It can be reasonably expected that the same happened for planetary systems other than our own. Detection of such Inter- stellar Objects (ISOs) would allow us to probe the planetesimal formation processes around other stars, possibly together with the effects of long-term exposure to the interstellar medium. 1I/2017 U1 'Oumuamua is the first known ISO, discovered by the Pan-STARRS1 telescope in October 2017$^3$.The discovery epoch photometry implies a highly elongated body with radii of $\sim 200 \times 20$ m when a comet-like geometric albedo of 0.04 is assumed. Here we report spectroscopic characterisation of 'Oumuamua, finding it to be variable with time but similar to organically rich surfaces found in the outer Solar System. The observable ISO population is expected to be dominated by comet-like bodies in agreement with our spectra, yet the reported inactivity implies a lack of surface ice. We show this is consistent with predictions of an insulating mantle produced by long-term cosmic ray exposure. An internal icy composition cannot therefore be ruled out by the lack of activity, even though 'Oumuamua passed within 0.25 au of the Sun.

Comet C/2017 K2 (PANSTARRS) was discovered by the Pan-STARRS1 (PS1) Survey on 2017 May 21 at a distance 16.09 au from the Sun, the second most distant discovery of an active comet. Pre-discovery images in the PS1 archive back to 2014 and additional deep CFHT images between 2013 May 10-13 showed the comet to be active at 23.75 au. We derive an upper limit to the nucleus radius of $R_N$=80 km, assuming a 4\% albedo. The spectral reflectivity of the comet surface is similar to "fresh" regions seen on comet 67P/Churyumov-Gerasimenko using the $Rosetta$ OSIRIS camera. Pre-discovery photometry combined with new data obtained with Megacam on the CFHT show that the activity is consistent with CO-ice sublimation and inconsistent with CO$_2$-ice sublimation. The ice sublimation models were run out to perihelion in 2022 at 1.8 au to predict the CO production rates, assuming that the outgassing area does not change. Assuming a canonical 4\% active surface area for water-ice sublimation, we present production rate ratios, $Q_{\rm CO}$/$Q_{\rm H2O}$, for a range of nucleus sizes. Comparing these results with other CO-rich comets we derive a lower limit to the nucleus radius of $\sim$14 km. We present predictions for $Q_{\rm CO}$ at a range of distances that will be useful for planning observations with JWST and large ground-based facilities.

Context. A lot of photometric data is produced by surveys such as Pan-STARRS, LONEOS, WISE or Catalina. These data are a rich source of information about the physical properties of asteroids. There are several possible approaches for utilizing these data. Lightcurve inversion is a typical method that works with individual asteroids. Our approach in this paper is statistical when we focused on large groups of asteroids like dynamical families and taxonomic classes, and the data were not sufficient for individual models. Aims. Our aim was to study the distributions of shape elongation $b/a$ and the spin axis latitude $\beta$ for various subpopulations of asteroids and to compare our results, based on Pan-STARRS1 survey, with statistics previously done using different photometric data (Lowell database, WISE data). Methods. We use the LEADER algorithm to compare the $b/a$ and $\beta$ distributions for different subpopulations of asteroids. The algorithm creates a cumulative distributive function (CDF) of observed brightness variations, and computes the $b/a$ and $\beta$ distributions using analytical basis functions that yield the observed CDF. A variant of LEADER is used to solve the joint distributions for synthetic populations to test the validity of the method. Results. When comparing distributions of shape elongation for groups of asteroids with different diameters $D$, we found that there are no differences for $D < 25$ km. We also constructed distributions for asteroids with different rotation periods and revealed that the fastest rotators with $P = 0 - 4$ h are more spheroidal than the population with $P = 4 - 8$ h.

We report the discovery of a $H_r = 3.4\pm0.1$ dwarf planet candidate by the Pan-STARRS Outer Solar System Survey. 2010 JO$_{179}$ is red with $(g-r)=0.88 \pm 0.21$, roughly round, and slowly rotating, with a period of $30.6$ hr. Estimates of its albedo imply a diameter of 600--900~km. Observations sampling the span between 2005--2016 provide an exceptionally well-determined orbit for 2010 JO$_{179}$, with a semi-major axis of $78.307\pm0.009$ au, distant orbits known to this precision are rare. We find that 2010 JO$_{179}$ librates securely within the 21:5 mean-motion resonance with Neptune on hundred-megayear time scales, joining the small but growing set of known distant dwarf planets on metastable resonant orbits. These imply a substantial trans-Neptunian population that shifts between stability in high-order resonances, the detached population, and the eroding population of the scattering disk.

Eos family was created during a catastrophic impact about 1.3 Gyr ago. Rotation states of individual family members contain information about the history of the whole population. We aim to increase the number of asteroid shape models and rotation states within the Eos collision family, as well as to revise previously published shape models from the literature. Such results can be used to constrain theoretical collisional and evolution models of the family, or to estimate other physical parameters by a thermophysical modeling of the thermal infrared data. We use all available disk-integrated optical data (i.e., classical dense-in-time photometry obtained from public databases and through a large collaboration network as well as sparse-in-time individual measurements from a few sky surveys) as input for the convex inversion method, and derive 3D shape models of asteroids together with their rotation periods and orientations of rotation axes. We present updated shape models for 15 asteroids and new shape model determinations for 16 asteroids. Together with the already published models from the publicly available DAMIT database, we compiled a sample of 56 Eos family members with known shape models that we used in our analysis of physical properties within the family. Rotation states of asteroids smaller than ~20 km are heavily influenced by the YORP effect, whilst the large objects more or less retained their rotation state properties since the family creation. Moreover, we also present a shape model and bulk density of asteroid (423) Diotima, an interloper in the Eos family, based on the disk-resolved data obtained by the Near InfraRed Camera (Nirc2) mounted on the W.M. Keck II telescope.

We present deep imaging observations, orbital dynamics, and dust tail model analyses of the double-component asteroid P/2016 J1 (J1-A and J1-B). The observations were acquired at the Gran Telescopio Canarias (GTC) and the Canada-France-Hawaii Telescope (CFHT) from mid March to late July, 2016. A statistical analysis of backward-in-time integrations of the orbits of a large sample of clone objects of P/2016 J1-A and J1-B shows that the minimum separation between them occurred most likely $\sim$2300 days prior to the current perihelion passage, i.e., during the previous orbit near perihelion. This closest approach was probably linked to a fragmentation event of their parent body. Monte Carlo dust tail models show that those two components became active simultaneously $\sim$250 days before the current perihelion, with comparable maximum loss rates of $\sim$0.7 kg s$^{-1}$ and $\sim$0.5 kg s$^{-1}$, and total ejected masses of 8$\times$10$^{6}$ kg and 6$\times$10$^{6}$ kg for fragments J1-A and J1-B, respectively. In consequence, the fragmentation event and the present dust activity are unrelated. The simultaneous activation times of the two components and the fact that the activity lasted 6 to 9 months or longer, strongly indicate ice sublimation as the most likely mechanism involved in the dust emission process.

We derived 90% confidence limits (CL) on the interstellar number density ($\rho_{IS}^{CL}$) of interstellar objects (ISO; comets and asteroids) as a function of the slope of their size-frequency distribution and limiting absolute magnitude. To account for gravitational focusing, we first generated a quasi-realistic ISO population to ~750 au from the Sun and propagated it forward in time to generate a steady state population of ISOs with heliocentric distance <50 au. We then simulated the detection of the synthetic ISOs using pointing data for each image and average detection efficiencies for each of three contemporary solar system surveys --- PS1, the Mt. Lemmon Survey, and the Catalina Sky Survey. These simulations allowed us to determine the surveys' combined ISO detection efficiency under several different but realistic modes of identifying ISOs in the survey data. Some of the synthetic detected ISOs had eccentricities as small as 1.01 --- in the range of the largest eccentricities of several known comets. Our best CL of $\rho_{IS}^{CL} = 1.4 \times 10^{-4}$ au$^{-3}$ implies that the expectation that extra-solar systems form like our solar system, eject planetesimals in the same way, and then distribute them throughout the galaxy, is too simplistic, or that the SFD or behavior of ISOs as they pass through our solar system is far from expectations.

We determine the absolute magnitude (H) distribution (or size-frequency distribution, SFD; $N(H) \propto 10^{\alpha H}$ where $\alpha$ 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 $\alpha = 0.48\pm0.02$ for NEOs with $13<H<16$, $\alpha = 0.33\pm0.01$ for NEOs with $16<H<22$, and $\alpha = 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 $\alpha = 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.

In this work we report the detection of seven Neptune Trojans (NTs) in the Pan-STARRS 1 (PS1) survey. Five of these are new discoveries, consisting of four L4 Trojans and one L5 Trojan. Our orbital simulations show that the L5 Trojan stably librates for only several million years. This suggests that the L5 Trojan must be of recent capture origin. On the other hand, all four new L4 Trojans stably occupy the 1:1 resonance with Neptune for more than 1 Gyr. They can, therefore, be of primordial origin. Our survey simulation results show that the inclination width of the Neptune Trojan population should be between $7^{\circ}$ and $27^{\circ}$ at $>$ 95% confidence, and most likely $\sim 11^{\circ}$. In this paper, we describe the PS1 survey, the Outer Solar System pipeline, the confirming observations, and the orbital/physical properties of the new Neptune Trojans.

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-$\sigma$. 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.

We present a method to identify distant solar system objects in long-term wide-field asteroid survey data, and conduct a search for them in the Pan-STARRS1 (PS1) image data acquired from 2010 to mid-2015. We demonstrate that our method is able to find multi-opposition orbital links, and we present the resulting orbital distributions which consist of 154 Centaurs, 255 classical Trans-Neptunian Objects (TNOs), 121 resonant TNOs, 89 Scattered Disc Objects (SDOs) and 10 comets. Our results show more than half of these are new discoveries, including a newly discovered 19th magnitude TNO. Our identified objects do not show clustering in their argument of perihelia, which if present, might support the existence of a large unknown planetary-sized object in the outer solar system.

The near-Earth asteroid (3200) Phaethon is an intriguing object: its perihelion is at only 0.14 au and is associated with the Geminid meteor stream. We aim to use all available disk-integrated optical data to derive a reliable convex shape model of Phaethon. By interpreting the available space- and ground-based thermal infrared data and Spitzer spectra using a thermophysical model, we also aim to further constrain its size, thermal inertia, and visible geometric albedo. We applied the convex inversion method to the new optical data obtained by six instruments and to previous observations. The convex shape model was then used as input for the thermophysical modeling. We also studied the long-term stability of Phaethon's orbit and spin axis with a numerical orbital and rotation-state integrator. We present a new convex shape model and rotational state of Phaethon: a sidereal rotation period of 3.603958(2) h and ecliptic coordinates of the preferred pole orientation of (319$^{\circ}$, $-$39$^{\circ}$) with a 5$^{\circ}$ uncertainty. Moreover, we derive its size ($D$=5.1$\pm$0.2 km), thermal inertia ($\Gamma$=600$\pm$200 J m$^{-2}$ s$^{-1/2}$ K$^{-1}$), geometric visible albedo ($p_{\mathrm{V}}$=0.122$\pm$0.008), and estimate the macroscopic surface roughness. We also find that the Sun illumination at the perihelion passage during the past several thousand years is not connected to a specific area on the surface, which implies non-preferential heating.

The rotational state of asteroids is controlled by various physical mechanisms including collisions, internal damping and the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. We have analysed the changes in magnitude between consecutive detections of approximately 60,000 asteroids measured by the PanSTARRS 1 survey during its first 18 months of operations. We have attempted to explain the derived brightness changes physically and through the application of a simple model. We have found a tendency toward smaller magnitude variations with decreasing diameter for objects of 1 < D < 8 km. Assuming the shape distribution of objects in this size range to be independent of size and composition our model suggests a population with average axial ratios 1 : 0.85 \pm 0.13 : 0.71 \pm 0.13, with larger objects more likely to have spin axes perpendicular to the orbital plane.

Using the Wide Field Camera for the United Kingdom Infrared Telescope, we measure the near-infrared colors of near-Earth objects (NEOs) in order to put constraints on their taxonomic classifications. The rapid-response character of our observations allows us to observe NEOs when they are close to the Earth and bright. Here we present near-infrared color measurements of 86 NEOs, most of which were observed within a few days of their discovery, allowing us to characterize NEOs with diameters of only a few meters. Using machine-learning methods, we compare our measurements to existing asteroid spectral data and provide probabilistic taxonomic classifications for our targets. Our observations allow us to distinguish between S-complex, C/X-complex, D-type, and V-type asteroids. Our results suggest that the fraction of S-complex asteroids in the whole NEO population is lower than the fraction of ordinary chondrites in the meteorite fall statistics. Future data obtained with UKIRT will be used to investigate the significance of this discrepancy.

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).

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.

Efficient identification and follow-up of astronomical transients is hindered by the need for humans to manually select promising candidates from data streams that contain many false positives. These artefacts arise in the difference images that are produced by most major ground-based time domain surveys with large format CCD cameras. This dependence on humans to reject bogus detections is unsustainable for next generation all-sky surveys and significant effort is now being invested to solve the problem computationally. In this paper we explore a simple machine learning approach to real-bogus classification by constructing a training set from the image data of ~32000 real astrophysical transients and bogus detections from the Pan-STARRS1 Medium Deep Survey. We derive our feature representation from the pixel intensity values of a 20x20 pixel stamp around the centre of the candidates. This differs from previous work in that it works directly on the pixels rather than catalogued domain knowledge for feature design or selection. Three machine learning algorithms are trained (artificial neural networks, support vector machines and random forests) and their performances are tested on a held-out subset of 25% of the training data. We find the best results from the random forest classifier and demonstrate that by accepting a false positive rate of 1%, the classifier initially suggests a missed detection rate of around 10%. However we also find that a combination of bright star variability, nuclear transients and uncertainty in human labelling means that our best estimate of the missed detection rate is approximately 6%.

The Pan-STARRS1 (PS1) survey has obtained imaging in 5 bands (grizy_P1) over 10 Medium Deep Survey (MDS) fields covering a total of 70 square degrees. This paper describes the search for apparently hostless supernovae (SNe) within the first year of PS1 MDS data with an aim of discovering new superluminous supernovae (SLSNe). A total of 249 hostless transients were discovered down to a limiting magnitude of M_AB ~ 23.5, of which 76 were classified as Type Ia SNe. There were 57 SNe with complete light curves that are likely core-collapse SNe (CCSNe) or SLSNe and 12 of these have had spectra taken. Of these 12 hostless, non-Type Ia SNe, 7 were SLSNe of Type Ic at redshifts between 0.5-1.4. This illustrates that the discovery rate of Type Ic SLSNe can be maximised by concentrating on hostless transients and removing normal SNe Ia. We present data for two new possible SLSNe; PS1-10pm (z = 1.206) and PS1-10ahf (z = 1.1), and estimate the rate of SLSNe-Ic to be between 3^{+3}_{-2} * 10^{-5} and 8^{+2}_{-1} * 10^{-5} of the CCSNe rate within 0.3 <= z <= 1.4 by applying a Monte-Carlo technique. The rate of slowly evolving, SN2007bi-like explosions is estimated as a factor of 10 lower than this range.

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.

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.

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.

In the past decade, several rapidly-evolving transients have been discovered whose timescales and luminosities are not easily explained by traditional supernovae (SN) models. The sample size of these objects has remained small due, at least in part, to the challenge of detecting short timescale transients with traditional survey cadences. Here we present the results from a search within the Pan-STARRS1 Medium Deep Survey (PS1-MDS) for rapidly-evolving and luminous transients. We identify 10 new transients with a time above half-maximum of less than 12 days and -16.5 > M > -20 mag. This increases the number of known events in this region of SN phase space by roughly a factor of three. The median redshift of the PS1-MDS sample is z=0.275 and they all exploded in star forming galaxies. In general, the transients possess faster rise than decline timescale and blue colors at maximum light (g - r < -0.2). Best fit blackbodies reveal photospheric temperatures/radii that expand/cool with time and explosion spectra taken near maximum light are dominated by a blue continuum, consistent with a hot, optically thick, ejecta. We find it difficult to reconcile the short timescale, high peak luminosity (L > 10^43 erg/s), and lack of UV line blanketing observed in many of these transients with an explosion powered mainly by the radioactive decay of Ni-56. Rather, we find that many are consistent with either (1) cooling envelope emission from the explosion of a star with a low-mass extended envelope which ejected very little (<0.03 M_sun) radioactive material, or (2) a shock breakout within a dense, optically thick, wind surrounding the progenitor star. After calculating the detection efficiency for objects with rapid timescales in the PS1-MDS we find a volumetric rate of 4800 - 8000 events/yr/Gpc^3 (4-7% of the core-collapse SN rate at z=0.2).

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.

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\tau_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\tau_{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.

P/2011 S1 (Gibbs) is an outer solar system comet or active Centaur with a similar orbit to that of the famous 29P/Schwassmann-Wachmann 1. P/2011 S1 (Gibbs) has been observed by the Pan-STARRS 1 (PS1) sky survey from 2010 to 2012. The resulting data allow us to perform multi-color studies of the nucleus and coma of the comet. Analysis of PS1 images reveals that P/2011 S1 (Gibbs) has a small nucleus $< 4$ km radius, with colors $g_{P1}-r_{P1} = 0.5 \pm 0.02$, $r_{P1}-i_{P1} = 0.12 \pm 0.02$ and $i_{P1}-z_{P1} = 0.46 \pm 0.03$. The comet remained active from 2010 to 2012, with a model-dependent mass-loss rate of $\sim100$ kg s$^{-1}$. The mass-loss rate per unit surface area of P/2011 S1 (Gibbs) is as high as that of 29P/Schwassmann-Wachmann 1, making it one of the most active Centaurs. The mass-loss rate also varies with time from $\sim 40$ kg s$^{-1}$ to 150 kg s$^{-1}$. Due to its rather circular orbit, we propose that P/2011 S1 (Gibbs) has 29P/Schwassmann-Wachmann 1-like outbursts that control the outgassing rate. The results indicate that it may have a similar surface composition to that of 29P/Schwassmann-Wachmann 1. Our numerical simulations show that the future orbital evolution of P/2011 S1 (Gibbs) is more similar to that of the main population of Centaurs than to that of 29P/Schwassmann-Wachmann 1. The results also demonstrate that P/2011 S1 (Gibbs) is dynamically unstable and can only remain near its current orbit for roughly a thousand years.

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.

Membership identification is the first step to determine the properties of a star cluster. Low-mass members in particular could be used to trace the dynamical history, such as mass segregation, stellar evaporation, or tidal stripping, of a star cluster in its Galactic environment. We identified member candidates with stellar masses $\sim$0.11--2.4 Msun of the intermediate-age Praesepe cluster (M44), by using Pan-STARRS and 2MASS photometry, and PPMXL proper motions. Within a sky area of 3 deg radius, 1040 candidates are identified, of which 96 are new inclusions. Using the same set of selection criteria on field stars led to an estimate of a false positive rate 16%, suggesting 872 of the candidates being true members. This most complete and reliable membership list allows us to favor the BT-Settl model in comparison with other stellar models. The cluster shows a distinct binary track above the main sequence, with a binary frequency of 20--40%, and a high occurrence rate of similar mass pairs. The mass function is consistent with that of the disk population but shows a deficit of members below 0.3 solar masses. A clear mass segregation is evidenced, with the lowest-mass members in our sample being evaporated from this disintegrating cluster.

We propose a low-cost robotic optical survey aimed at $1-300$ m Near Earth Objects (NEO) based on four state-of-the-art telescopes having extremely wide field of view. The small Near-Earth Asteroids (NEA) represent a potential risk but also easily accessible space resources for future robotic or human space in-situ exploration, or commercial activities. The survey system will be optimized for the detection of fast moving - trailed - asteroids, space debris and will provide real-time alert notifications. The expected cost of the system including 1-year development and 2-year operation is 1,000,000 EUR. The successful demonstration of the system will promote cost-efficient ADAM-WFS (Automatic Detection of Asteroids and Meteoroids - A Wide Field Survey) systems to be built around the world.

We have used Pan-STARRS1 to discover an extremely red late-L dwarf, which has (J-K)_MKO = 2.84 and (J-K)_2MASS = 2.78, making it the reddest known field dwarf and second only to 2MASS J1207-39b among substellar companions. Near-IR spectroscopy shows a spectral type of L7 and reveals a triangular H-band continuum and weak alkali (K I and Na I) lines, hallmarks of low surface gravity. Near-IR astrometry from the Hawaii Infrared Parallax Program gives a distance of 24.6+/-1.4 pc and indicates a much fainter J-band absolute magnitude than field L dwarfs. The position and kinematics of PSO J318-22 point to membership in the beta Pictoris moving group. Evolutionary models give a temperature of 1160 (-40,+30) K and a mass of 6.5 (-1.0, +1.3) Mjup, making PSO J318-22 one of the lowest mass free-floating objects in the solar neighborhood. This object adds to the growing list of low-gravity field L dwarfs and is the first to be strongly deficient in methane relative to its estimated temperature. Comparing their spectra suggests that young L dwarfs with similar ages and temperatures can have different spectral signatures of youth. For the two objects with well constrained ages (PSO J318-22 and 2MASS J0355+11), we find their temperatures are ~400 K cooler than field objects of similar spectral type but their luminosities are comparable, i.e., these young L dwarfs are very red and unusually cool but not "underluminous." Altogether, PSO J318-22 is the first free-floating object with the colors, magnitudes, spectrum, luminosity, and mass that overlap the young dusty planets around HR 8799 and 2MASS J1207-39.

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.

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.

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.

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.

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.

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.

We have used numerical routines to model the evolution of a simulated Baptistina family to constrain its age in light of new measurements of the diameters and albedos of family members from the Wide-field Infrared Survey Explorer. We also investigate the effect of varying the assumed physical and orbital parameters on the best-fitting age. We find that the physically allowed range of assumed values for the density and thermal conductivity induces a large uncertainty in the rate of evolution. When realistic uncertainties in the family members' physical parameters are taken into account we find the best-fitting age can fall anywhere in the range of 140-320 Myr. Without more information on the physical properties of the family members it is difficult to place a more firm constraint on Baptistina's age.

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{\sigma}. The most statistically significant cluster at the \sim 2{\sigma} 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.

We present the preliminary results of an analysis of the sub-populations within the near-Earth asteroids, including the Atens, Apollos, Amors, and those that are considered potentially hazardous using data from the Wide-field Infrared Survey Explorer (WISE). In order to extrapolate the sample of objects detected by WISE to the greater population, we determined the survey biases for asteroids detected by the project's automated moving object processing system (known as NEOWISE) as a function of diameter, visible albedo, and orbital elements. Using this technique, we are able to place constraints on the number of potentially hazardous asteroids (PHAs) larger than 100 m and find that there are $\sim4700\pm1450$ such objects. As expected, the Atens, Apollos, and Amors are revealed by WISE to have somewhat different albedo distributions, with the Atens being brighter than the Amors. The cumulative size distributions of the various near-Earth object (NEO) subgroups vary slightly between 100 m and 1 km. A comparison of the observed orbital elements of the various sub-populations of the NEOs with the current best model is shown.

We report on our serendipitous pre-discovery detection and detailed follow-up of the broad-lined Type Ic supernova (SN) 2010ay at z = 0.067 imaged by the Pan-STARRS1 3pi survey just ~4 days after explosion. The SN had a peak luminosity, M_R ~ -20.2 mag, significantly more luminous than known GRB-SNe and one of the most luminous SNe Ib/c ever discovered. The absorption velocity of SN 2010ay is v_Si ~ 19,000 km/s at ~40 days after explosion, 2-5 times higher than other broad-lined SNe and similar to the GRB-SN 2010bh at comparable epochs. Moreover, the velocity declines ~2 times slower than other SNe Ic-BL and GRB-SNe. Assuming that the optical emission is powered by radioactive decay, the peak magnitude implies the synthesis of an unusually large mass of 56 Ni, M_Ni = 0.9 M_solar. Modeling of the light-curve points to a total ejecta mass, M_ej ~ 4.7 M_sol, and total kinetic energy, E_K ~ 11x10^51 ergs. The ratio of M_Ni to M_ej is ~2 times as large for SN 2010ay as typical GRB-SNe and may suggest an additional energy reservoir. The metallicity (log(O/H)_PP04 + 12 = 8.19) of the explosion site within the host galaxy places SN 2010ay in the low-metallicity regime populated by GRB-SNe, and ~0.5(0.2) dex lower than that typically measured for the host environments of normal (broad-lined) Ic supernovae. We constrain any gamma-ray emission with E_gamma < 6x10^{48} erg (25-150 keV) and our deep radio follow-up observations with the Expanded Very Large Array rule out relativistic ejecta with energy, E > 10^48 erg. We therefore rule out the association of a relativistic outflow like those which accompanied SN 1998bw and traditional long-duration GRBs, but place less-stringent constraints on a weak afterglow like that seen from XRF 060218. These observations challenge the importance of progenitor metallicity for the production of a GRB, and suggest that other parameters also play a key role.

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%.

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.

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.

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.

With the NEOWISE portion of the \emph{Wide-field Infrared Survey Explorer} (WISE) project, we have carried out a highly uniform survey of the near-Earth object (NEO) population at thermal infrared wavelengths ranging from 3 to 22 $\mu$m, allowing us to refine estimates of their numbers, sizes, and albedos. The NEOWISE survey detected NEOs the same way whether they were previously known or not, subject to the availability of ground-based follow-up observations, resulting in the discovery of more than 130 new NEOs. The survey's uniformity in sensitivity, observing cadence, and image quality have permitted extrapolation of the 428 near-Earth asteroids (NEAs) detected by NEOWISE during the fully cryogenic portion of the WISE mission to the larger population. We find that there are 981$\pm$19 NEAs larger than 1 km and 20,500$\pm$3000 NEAs larger than 100 m. We show that the Spaceguard goal of detecting 90% of all 1 km NEAs has been met, and that the cumulative size distribution is best represented by a broken power law with a slope of 1.32$\pm$0.14 below 1.5 km. This power law slope produces $\sim13,200\pm$1,900 NEAs with $D>$140 m. Although previous studies predict another break in the cumulative size distribution below $D\sim$50-100 m, resulting in an increase in the number of NEOs in this size range and smaller, we did not detect enough objects to comment on this increase. The overall number for the NEA population between 100-1000 m are lower than previous estimates. The numbers of near-Earth comets will be the subject of future work.

The seven known main belt comets (MBCs) have orbital characteristics of main belt asteroids yet exhibit dust ejection like comets. In order to constrain their physical and orbital properties we searched the Thousand Asteroid Light Curve Survey (TALCS; Masiero et al. 2009) for additional candidates using two diagnostics: tail and coma detection. This was the most sensitive MBC survey effort to date, extending the search from MBCs with H~18 (D~1 km) to H~21 (D~150 m). We fit each of the 924 TALCS objects to a PSF model incorporating both a coma and nuclear component to measure the fractional contribution of the coma to the total surface brightness. We determined the significance of the coma detection using the same algorithm on a sample of comparable null detections. We did not identify any MBC candidates with this technique to a sensitivity limit on the order of cometary mass loss rate of about 0.1 kg/s. Our tail detection algorithm identified statistically significant flux in a segmented annulus around the candidate object. We show that the technique can detect tail activity throughout the asteroid belt to the level of the currently known MBCs. Although we did not identify any MBC candidates with this technique, we find a statistically significant detection of faint activity in the entire ensemble of TALCS asteroids. This suggests that many main belt asteroids are active at very low levels. We set 90% upper confidence limits on the number distribution of MBCs as a function of absolute magnitude, semimajor axis, eccentricity, and inclination. There are <~ 400000 MBCs in the main belt brighter than H_V=21 (~150 m) and the MBC:MBA ratio is <~ 1:400. We further comment on the ability of observations to meaningfully constrain the snow line's location. Under some reasonable and simple assumptions we claim 85% confidence that the contemporary snow line lies beyond 2.5 AU.

The \emph{Wide-field Infrared Survey Explorer} has surveyed the entire sky at four infrared wavelengths with greatly improved sensitivity and spatial resolution compared to its predecessors, the \emph{Infrared Astronomical Satellite} and the \emph{Cosmic Background Explorer}. NASA's Planetary Science Division has funded an enhancement to the \WISE\ data processing system called "NEOWISE" that allows detection and archiving of moving objects found in the \WISE\ data. NEOWISE has mined the \WISE\ images for a wide array of small bodies in our Solar System, including Near-Earth Objects (NEOs), Main Belt asteroids, comets, Trojans, and Centaurs. By the end of survey operations in February 2011, NEOWISE identified over 157,000 asteroids, including more than 500 NEOs and $\sim$120 comets. The NEOWISE dataset will enable a panoply of new scientific investigations.

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 $\tau$' space weathering model of \citet{bib.wil10}. We fit the size-age distribution to the enhanced dual $\tau$ model and found characteristic weathering and gardening times of $\tau_w = 2050 \pm 80$ Myr and $\tau_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 $\Delta 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$.

We present the early UV and optical light curve of Type IIP supernova (SN) 2010aq at z=0.0862, and compare it to analytical models for thermal emission following SN shock breakout in a red supergiant star. SN 2010aq was discovered in joint monitoring between the Galaxy Evolution Explorer (GALEX) Time Domain Survey (TDS) in the NUV and the Pan-STARRS1 Medium Deep Survey (PS1 MDS) in the g, r, i, and z bands. The GALEX and Pan-STARRS1 observations detect the SN less than 1 day after shock breakout, measure a diluted blackbody temperature of 31,000 +/- 6,000 K 1 day later, and follow the rise in the UV/optical light curve over the next 2 days caused by the expansion and cooling of the SN ejecta. The high signal-to-noise ratio of the simultaneous UV and optical photometry allows us to fit for a progenitor star radius of 700 +/- 200 R_sun, the size of a red supergiant star. An excess in UV emission two weeks after shock breakout compared to SNe well fitted by model atmosphere-code synthetic spectra with solar metallicity, is best explained by suppressed line blanketing due to a lower metallicity progenitor star in SN 2010aq. Continued monitoring of PS1 MDS fields by the GALEX TDS will increase the sample of early UV detections of Type II SNe by an order of magnitude, and probe the diversity of SN progenitor star properties.

We have extended our earlier work on space weathering of the youngest S-complex asteroid families to include results from asteroid clusters with ages <10^6 years and to newly identified asteroid pairs with ages <5x10^5 years. We have identified three S-complex asteroid clusters with ages in the range 10^{5-6} years. The average color of the objects in these clusters agree with the prediction of Willman et al., 2008. SDSS photometry of the members of very young asteroid pairs with ages <10^5 years was used to determine their taxonomy. The average color of the S-complex pairs is PC_1=0.49+/-0.03, over 5-sigma redder than predicted by Willman et al., 2008. Therefore, the most likely pair formation mechanism is gentle separation due to YORP spin-up leaving much of the aged and reddened surface undisturbed. In this case our color measurement allows us to set an upper limit of ~64% on the disturbed surface portion. Using pre-existing color data and our new results for the youngest S-complex asteroid clusters we have extended our space weather model to explicitly include the effects of regolith gardening and fit separate weathering and gardening characteristic timescales of tau_w=960+/-160My and tau_g=2000+/-290My respectively. The first principal component color for fresh S-complex material is 0.37+/-0.01 while the maximum amount of local reddening is 0.33+/-0.06. Our first-ever determination of the gardening time is in stark contrast to our calculated gardening time of tau_g~270My based on main belt impact rates and reasonable assumptions about crater and ejecta blanket sizes. A possible resolution for the discrepancy is through a `honeycomb' mechanism in which the surface regolith structure absorbs small impactors without producing significant ejecta. This mechanism could also account for the paucity of small craters on (433) Eros.

Rewriting a shorter abstract.

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.

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.

We present the observational results of a survey designed to target and detect asteroids whose colors are similar to those of Vesta family members and thus may be considered as candidates for having a basaltic composition. Fifty basaltic candidates were selected with orbital elements that lie outside of the Vesta dynamical family. Optical and near-infrared spectra were used to assign a taxonomic type to 11 of the 50 candidates. Ten of these were spectroscopically confirmed as V-type asteroids, suggesting that most of the candidates are basaltic and can be used to constrain the distribution of basaltic material in the Main Belt. Using our catalog of V-type candidates and the success rate of the survey, we calculate unbiased size-frequency and semi-major axis distributions of V-type asteroids. These distributions, in addition to an estimate for the total mass of basaltic material, suggest that Vesta was the predominant contributor to the basaltic asteroid inventory of the Main Belt, however scattered planetesimals from the inner Solar System (a < 2.0 AU) and other partially/fully differentiated bodies likely contributed to this inventory. In particular, we infer the presence of basaltic fragments in the vicinity of asteroid 15 Eunomia, which may be derived from a differentiated parent body in the middle Main Belt (2.5 < a < 2.8). We find no asteroidal evidence for a large number of previously undiscovered basaltic asteroids, which agrees with previous theories suggesting that basaltic fragments from the ~100 differentiated parent bodies represented in meteorite collections have been "battered to bits" [Burbine, T.H., Meibom, A., Binzel, R.P., 1996. Mantle material in the Main Belt: Battered to bits? Met. & Planet. Sci. 31, 607].

We present visible and near-infrared reflectance spectra and interpreted surface mineralogy for asteroid 10537 (1991 RY16). The spectrum of this object is without precedent amongst the Main Belt asteroids. A unique absorption band centered at 0.63 microns could be attributed to one of several mineralogies. Pronounced 1- and 2-micron absorption bands suggest that the composition of 10537 is a mixture of pyroxenes and olivine and that it originated from a parent body that was partially or fully differentiated. The closest available analog is the large Main Belt asteroid 349 Dembowska but 10537 may be an isolated fragment from a completely eroded parent body.

Given a set of astrometric observations of the same object, the problem of orbit determination is to compute the orbit and to assess its uncertainty and reliability. For the next generation surveys, with much larger number density of observed objects, new algorithms or substantial revisions of the classical ones are needed. The problem has three main steps, preliminary orbit, least squares orbit, and quality control. The classical theory of preliminary orbits was incomplete: the consequences of the topocentric correction had not been fully studied. We show that it is possible to account for this correction, possibly with an increase in the number of preliminary solutions, without impairing the overall orbit determination performance. We have developed modified least squares orbit determination algorithms that can be used to improve the reliability of the procedure. We have tested the complete procedure on two simulations with number densities comparable to that expected from the next generation surveys such as Pan-STARRS and LSST. To control the problem of false identifications we have introduced a quality control on the fit residuals based on an array of metrics and a procedure to remove duplications and contradictions in the output. The results confirm that large sets of discoveries can be obtained with good quality orbits and very high success rate losing only 0.6 to 1.3% of objects and a false identification rate in the range 0.02 to 0.06%.

The Panoramic Survey Telescope And Rapid Response System (Pan-STARRS) under development at the University of Hawaii's Institute for Astronomy is creating the first fully automated end-to-end Moving Object Processing System (MOPS) in the world. It will be capable of identifying detections of moving objects in our solar system and linking those detections within and between nights, attributing those detections to known objects, calculating initial and differentially-corrected orbits for linked detections, precovering detections when they exist, and orbit identification. Here we describe new kd-tree and variable-tree algorithms that allow fast, efficient, scalable linking of intra and inter-night detections. Using a pseudo-realistic simulation of the Pan-STARRS survey strategy incorporating weather, astrometric accuracy and false detections we have achieved nearly 100% efficiency and accuracy for intra-night linking and nearly 100% efficiency for inter-night linking within a lunation. At realistic sky-plane densities for both real and false detections the intra-night linking of detections into `tracks' currently has an accuracy of 0.3%. Successful tests of the MOPS on real source detections from the Spacewatch asteroid survey indicate that the MOPS is capable of identifying asteroids in real data.

We estimate the total number and the slope of the size frequency distribution (SFD) of dormant Jupiter Family Comets (JFCs) by fitting a one-parameter model to the known population. We first select 61 Near Earth Objects (NEOs) that are likely to be dormant JFCs because their orbits are dynamically coupled to Jupiter (Bottke et al, 2002). Then, from the numerical simulations of Levison & Duncan (1997), we construct an orbit distribution model for JFCs in the NEO orbital element space. We assume an orbit independent SFD for all JFCs, the slope of which is our unique free parameter. Finally, we compute observational biases for dormant JFCs using a calibrated NEO survey simulator (Jedicke et al. 2003). By fitting the biased model to the data, we estimate that there are ~75 dormant JFCs with H<18 in the NEO region and that the slope of their cumulative SFD is -1.5 +/- 0.3. Our slope for the SFD of dormant JFCs is very close to that of active JFCs as determined by Weissman and Lowry (2003). Thus, we argue that when JFCs fade they are likely to become dormant rather than to disrupt and that the fate of faded comets is size independent. Our results imply that the size distribution of the JFC progenitors - the scattered disk trans-Neptunian population - either (i) has a similar and shallow SFD or (i') is slightly steeper and physical processes acting on the comets in a size-dependent manner creates the shallower active comet SFD. Our measured slope, typical of collisionally evolved populations with a size dependent impact strength (Benz and Asphaug 1999), suggests that scattered disk bodies reached collisional equilibrium inside the proto-planetary disk prior to their removal from the planetary region.

The study of asteroid families has provided tremendous insight into the forces that sculpted the main belt and continue to drive the collisional and dynamical evolution of asteroids. The identification of asteroid families within the NEO population could provide a similar boon to studies of their formation and interiors. In this study we examine the purported identification of NEO families by Drummond (2000) and conclude that it is unlikely that they are anything more than random fluctuations in the distribution of NEO osculating orbital elements. We arrive at this conclusion after examining the expected formation rate of NEO families, the identification of NEO groups in synthetic populations that contain no genetically related NEOs, the orbital evolution of the largest association identified by Drummond (2000), and the decoherence of synthetic NEO families intended to reproduce the observed members of the same association. These studies allowed us to identify a new criterion that can be used to select real NEO families for further study in future analyses, based on the ratio of the number of pairs and the size of strings to the number of objects in an identified association.

We present a detailed observational strategy for finding Near Earth Objects (NEOs) with the Sloan Digital Sky Survey (SDSS) telescope. We investigate strategies in normal, unbinned mode as well as binning the CCDs 2x2 or 3x3, which affects the sky coverage rate and the limiting apparent magnitude. We present results from 1 month, 3 year and 10 year simulations of such surveys. For each cadence and binning mode, we evaluate the possibility of achieving the Spaceguard goal of detecting 90% of 1 km NEOs (absolute magnitude H <= 18 for an albedo of 0.1). We find that an unbinned survey is most effective at detecting H <= 20 NEOs in our sample. However, a 3x3 binned survey reaches the Spaceguard Goal after only seven years of operation. As the proposed large survey telescopes (PanStarss; LSST) are at least 5-10 years from operation, an SDSS NEO survey could make a significant contribution to the detection and photometric characterization of the NEO population.

We have developed a model-independent analytical method for debiasing the four-dimensional (a,e,i,H) distribution obtained in any asteroid observation program and have applied the technique to results obtained with the 0.9m Spacewatch Telescope. From 1992 to 1995 Spacewatch observed ~3740 deg^2 near the ecliptic and made observations of more than 60,000 asteroids to a limiting magnitude of V~21. The debiased semi-major axis and inclination distributions of Main Belt asteroids in this sample with 11.5<= H <16 match the distributions of the known asteroids with H <11.5. The absolute magnitude distribution was studied in the range 8< H <17.5. We have found that the set of known asteroids is complete to about absolute magnitudes 12.75, 12.25 and 11.25 in the inner, middle and outer regions of the belt respectively. The number distribution as a function of absolute magnitude cannot be represented by a single power-law (10^{alpha H}) in any region. We were able to define broad ranges in H in each part of the belt where alpha was nearly constant. Within these ranges of H the slope does not correspond to the value of 0.5 expected for an equilibrium cascade in self-similar collisions (Dohnanyi, 1971). The value of alpha varies with absolute magnitude and shows a `kink' in all regions of the belt for H~13. This absolute magnitude corresponds to a diameter ranging from about 8.5 to 12.5 km depending on the albedo or region of the belt.