The discovery of the second interstellar object 2I/Borisov on 2019 August 30 raises the question of whether it was ejected recently from a nearby stellar system. Here we compute the asymptotic incoming trajectory of 2I/Borisov, based on both recent and pre-discovery data extending back to December 2018, using a range of force models that account for cometary outgassing. From Gaia DR2 astrometry and radial velocities, we trace back in time the Galactic orbits of 7.4 million stars to look for close encounters with 2I/Borisov. The closest encounter we find took place 910kyr ago with the M0V star Ross 573, at a separation of 0.068pc (90% confidence interval of 0.053-0.09pc) with a relative velocity of 23km/s. This encounter is nine times closer than the closest past encounter identified for the first interstellar object 1I/'Oumuamua. Ejection of 2I/Borisov via a three-body encounter in a binary or planetary system is possible, although such a large ejection velocity is unlikely to be obtained and Ross 573 shows no signs of binarity. We also identify and discuss some other recent close encounters, recognizing that if 2I/Borisov is more than about 10Myr old, our search would be unlikely to find its parent system.

Comet 2I/Borisov, the first unambiguous interstellar comet ever found, was discovered in August 2019 at $\sim3$ au from the Sun on its inbound leg. No pre-discovery detection beyond 3 au has yet been reported, mostly due to the comet's proximity to the Sun as seen from the Earth. Here we present a search for pre-discovery detections of comet Borisov using images taken by the Catalina Sky Survey (CSS), Pan-STARRS and Zwicky Transient Facility (ZTF), with a further comprehensive follow-up campaign being presented in \citet{Bolin2019}. We identified comet Borisov in ZTF images taken in May 2019 and use these data to update its orbit. This allowed us to identify the comet in images acquired as far back as December 2018, when it was 7.8 au from the Sun. The comet was not detected in November 2018 when it was 8.6 au from the Sun, possibly implying an onset of activity around this time. This suggests that the activity of the comet is either driven by a more volatile species other than H$_2$O, such as CO or CO$_2$, or by exothermic crystallization of amorphous ice. We derive the radius of the nucleus to be $<7$ km using the non-detection in November 2018, and estimate an area of $\sim0.5$---$10 \mathrm{km^2}$ has been active between December 2018 and September 2019, though this number is model-dependent and is highly uncertain. The behavior of comet Borisov during its inbound leg is observationally consistent with dynamically new comets observed in our solar system, suggesting some similarities between the two.

The discovery of the first active interstellar object 2I/Borisov provides an unprecedented opportunity to study planetary formation processes in another planetary system. In particular, spectroscopic observations of 2I allow us to constrain the composition of its nuclear ices. We obtained optical spectra of 2I with the 4.2 m William Herschel and 2.5 m Isaac Newton telescopes between 2019 September 30 and October 13, when the comet was between 2.5 au and 2.4 au from the Sun. We also imaged the comet with broadband filters on 15 nights from September 11 to October 17, as well as with a CN narrow-band filter on October 18 and 20, with the TRAPPIST-North telescope. Broadband imaging confirms that the dust coma colours (B-V=0.82$\pm$0.02, V-R=0.46$\pm$0.03, R-I=0.44$\pm$0.03, B-R=1.28$\pm$0.03) are the same as for Solar System comets. We detect CN emission in all spectra and in the TRAPPIST narrow-band images with production rates between 1.6$\times10^{24}$ and 2.1$\times10^{24}$ molec/s. No other species are detected. We determine three-sigma upper limits for C$_2$, C$_3$, and OH production rates of 6$\times10^{23}$ molec/s, 3$\times10^{23}$ molec/s and 2$\times10^{27}$ molec/s, respectively, on October 02. There is no significant increase of the CN production rate or A(0)f$\rho$ during our observing period. Finally, we place a three-sigma upper limit on the Q(C$_2$)/Q(CN) ratio of 0.3 (on October 13). From this, we conclude that 2I is highly depleted in C$_2$, and may have a composition similar to Solar System carbon-chain depleted comets.

The detection of Interstellar Objects passing through the Solar System offers the promise of constraining the physical and chemical processes involved in planetary formation in other extrasolar systems. While the effect of outgassing by 1I/2017 U1 ('Oumuamua) was dynamically observed, no direct detection of the ejected material was made. The discovery of the active interstellar comet 2I/Borisov means spectroscopic investigations of the sublimated ices is possible for this object. We report the first detection of gas emitted by an interstellar comet via the near-UV emission of CN from 2I/Borisov at a heliocentric distance of $r$ = 2.7 au on 2019 September 20. The production rate was found to be Q(CN) = $(3.7\pm0.4)\times10^{24}$ s$^{-1}$, using a simple Haser model with an outflow velocity of 0.5 km s$^{-1}$. No other emission was detected, with an upper limit to the production rate of C$_2$ of $4\times10^{24}$ s$^{-1}$. The spectral reflectance slope of the dust coma over $3900$ \AA\ $< \lambda< 6000$ \AA \ is steeper than at longer wavelengths, as found for other comets. Broad band $R_c$ photometry on 2019 September 19 gave a dust production rate of $Af\rho=143\pm10$ cm. Modelling of the observed gas and dust production rates constrains the nuclear radius to $0.7-3.3$ km assuming reasonable nuclear properties. Overall, we find the gas, dust and nuclear properties for the first active Interstellar Object are similar to normal Solar System comets.

Near-Earth asteroid (410777) 2009 FD is a potentially hazardous asteroid with potential impacts on Earth at the end of the 22nd century. The astrometry collected during the 2019 apparition provides information on the trajectory of (410777) by constraining the Yarkovsky effect, which is the main source of uncertainty for future predictions, and informing the impact hazard assessment. We included the Yarkovsky effect in the force model and estimated its magnitude from the fit to the (410777) optical and radar astrometric data. We performed the (410777) hazard assessment over 200 years by using two independent approaches: the NEODyS group adopted a generalisation of the Line Of Variations method in a 7-dimensional space, whereas the JPL team resorted to the Multilayer Clustered Sampling technique. We obtain a 4-sigma detection of the Yarkovsky effect acting on (410777), which corresponds to a semimajor axis drift of (3.6 +/- 0.9) x 10^{-3} au/Myr. As for the hazard results of both teams, the main impact possibility in 2185 is ruled out and the only remaining one is in 2190, but with a probability smaller than 10^{-8}.

We report on the catastrophic disintegration of P/2016 G1 (PANSTARRS), an active asteroid, in April 2016. Deep images over three months show that object is constituted by a central concentration of fragments surrounded by an elongated coma, and presents previously unreported sharp arc-like and a narrow linear features. The morphology and evolution of these characteristics independently point toward a brief event on 2016 March 6. The arc and the linear feature can be reproduced by large particles on a ring, moving at ~2.5 m/s. The expansion of the ring defines a cone with a ~40deg half-opening. We propose that G1 was hit by a small object which caused its (partial or total) disruption, and that the ring corresponds to large fragments ejected during the final stages of the crater formation.}

The Large Synoptic Survey Telescope (LSST) is expected to increase known small solar system object populations by an order of magnitude or more over the next decade, enabling a broad array of transformative solar system science investigations to be performed. In this white paper, we discuss software tools and infrastructure that we anticipate will be needed to conduct these investigations and outline possible approaches for implementing them. Feedback from the community or contributions to future updates of this work are welcome. Our aim is for this white paper to encourage further consideration of the software development needs of the LSST solar system science community, and also to be a call to action for working to meet those needs in advance of the expected start of the survey in late 2022.

On 2019 January 5 a streamer associated with the 4--10 km main-belt asteroid (6478)~Gault was detected by the ATLAS sky survey, a rare discovery of activity around a main-belt asteroid. Archival data from ATLAS and Pan-STARRS1 show the trail in early December 2018, but not between 2010 and January 2018. The feature has significantly changed over one month, perfectly matching predictions of pure dust dynamical evolution and changes in observing geometry for a short release of dust around 2018 October 28. Follow-up observations with HST show a second narrow trail corresponding to a brief release of dust on 2018 December 30. Both releases occurred with negligible velocity. We find the dust grains to be fairly large, with power-law size distributions in the $10^{-5} - 10^{-3}$~m range and power-law indices of $\sim -1.5$. Three runs of ground-based data find a signature of $\sim 2\,\rm h$ rotation, close to the rotational limit, suggesting that the activity is the result of landslides or reconfigurations after YORP spin-up.

When new objects are detected in the sky, an orbit determination needs to be performed immediately to find out their origin, to determine the probability of an Earth impact and possibly also to estimate the impact region on Earth. ESA's Flyeye telescope is expected to revolutionize the effort of predicting potential asteroid or deep space debris impact hazards due to the expected increase of near-Earth object discoveries. As the observed orbit arc for such an object is very short, classical Gaussian orbit determination cannot be used. We adopt the systematic ranging technique to overcome the lack of information and predict a region of the sky where the body can most likely be found. We also provide a detection probability for follow-up observations and investigate potential follow-up telescopes for the Flyeye telescope.

We present a photometric and dynamical study of comet C/2010 U3 (Boattini), which was seen active in prediscovery data as early as 2005 November at a new inbound record heliocentric distance $r_{\rm H} = 25.8$ au. Two outburst events around 2009 and 2017 were observed. The coma and tail of the comet consist of dust grains of $\sim$10 $\mu$m in radius, ejected protractedly at speeds $\lesssim$50 m s$^{-1}$ near the subsolar point, and are subjected to the Lorentz force, solar gravitation and radiation pressure force altogether. The prolonged activity indicates that sublimation of supervolatiles (e.g., CO, CO$_2$) is at play, causing a net mass-loss rate $\gtrsim$1 kg s$^{-1}$. To sustain the mass loss, the nucleus radius has to be $\gtrsim$0.1 km. The color of the cometary dust, similar to other long-period comets, is redder than the solar colors, but we also observed potential color variations when the comet was at $10 < r_{\rm H} < 15$ au, concurrent with the onset of crystallisation of amorphous water ice, if at all. Using publicly available and our refined astrometric measurements, we estimated the precise trajectory of the comet and propagated it backward to its previous perihelion. We found that the comet has visited the planetary region $1.96 \pm 0.04$ Myr ago, with barycentric perihelion distance $q = 8.364 \pm 0.004$ au. Thus, C/2010 U3 (Boattini) is almost certainly a dynamically old comet from the Oort cloud, and the observed activity cannot be caused by retained heat from the previous apparition.

1I/`Oumuamua is the first confirmed interstellar body in our Solar System. Here we report on observations of `Oumuamua made with the Spitzer Space Telescope on 2017 November 21--22 (UT). We integrated for 30.2~hours at 4.5 micron (IRAC channel 2). We did not detect the object and place an upper limit on the flux of 0.3 uJy (3sigma). This implies an effective spherical diameter less than [98, 140, 440] meters and albedo greater than [0.2, 0.1, 0.01] under the assumption of low, middle, or high thermal beaming parameter eta, respectively. With an aspect ratio for `Oumuamua of 6:1, these results correspond to dimensions of [240:40, 341:57, 1080:180] meters, respectively. We place upper limits on the amount of dust, CO, and CO2 coming from this object that are lower than previous results; we are unable to constrain the production of other gas species. Both our size and outgassing limits are important because `Oumuamua's trajectory shows non-gravitational accelerations that are sensitive to size and mass and presumably caused by gas emission. We suggest that `Oumuamua may have experienced low-level post-perihelion volatile emission that produced a fresh, bright, icy mantle. This model is consistent with the expected eta value and implied high albedo value for this solution, but, given our strict limits on CO and CO2, requires another gas species --- probably H2O --- to explain the observed non-gravitational acceleration. Our results extend the mystery of `Oumuamua's origin and evolution.

We report observations of the reactivations of main-belt comets 238P/Read and 288P/(300163) 2006 VW139, that also track the evolution of each object's activity over several months in 2016 and 2017. We additionally identify and analyze archival SDSS data showing 288P to be active in 2000, meaning that both 238P and 288P have now each been confirmed to be active near perihelion on three separate occasions. From data obtained of 288P from 2012-2015 when it appeared inactive, we find best-fit R-band H,G phase function parameters of H_R=16.80+/-0.12 mag and G_R=0.18+/-0.11, corresponding to effective component radii of r_c=0.80+/-0.04 km, assuming a binary system with equally-sized components. Fitting linear functions to ejected dust masses inferred for 238P and 288P soon after their observed reactivations in 2016, we find an initial average net dust production rate of 0.7+/-0.3 kg/s and a best-fit start date of 2016 March 11 (when the object was at a true anomaly of -63 deg) for 238P, and an initial average net dust production rate of 5.6+/-0.7 kg/s and a best-fit start date of 2016 August 5 (when the object was at a true anomaly of -27 deg) for 288P. Applying similar analyses to archival data, we find similar start points for previous active episodes for both objects, suggesting that minimal mantle growth or ice recession occurred between the active episodes in question. Some changes in dust production rates between active episodes are detected, however. More detailed dust modeling is suggested to further clarify the process of activity evolution in main-belt comets.

The first detected interstellar object 'Oumuamua that passed within 0.25au of the Sun on 2017 September 9 was presumably ejected from a stellar system. We use its newly determined non-Keplerian trajectory together with the reconstructed Galactic orbits of 7 million stars from Gaia DR2 to identify past close encounters. Such an "encounter" could reveal the home system from which 'Oumuamua was ejected. The closest encounter, at 0.60pc (0.53-0.67pc, 90% confidence interval), was with the M2.5 dwarf HIP 3757 at a relative velocity of 24.7km/s, 1Myr ago. A more distant encounter (1.6pc) but with a lower encounter (ejection) velocity of 10.7km/s was with the G5 dwarf HD 292249, 3.8Myr ago. Two more stars have encounter distances and velocities intermediate to these. The encounter parameters are similar across six different non-gravitational trajectories for 'Oumuamua. Ejection of 'Oumuamua by scattering from a giant planet in one of the systems is plausible, but requires a rather unlikely configuration to achieve the high velocities found. A binary star system is more likely to produce the observed velocities. None of the four home candidates have published exoplanets or are known to be binaries. Given that the 7 million stars in Gaia DR2 with 6D phase space information is just a small fraction of all stars for which we can eventually reconstruct orbits, it is a priori unlikely that our current search would find 'Oumuamua's home star system. As 'Oumuamua is expected to pass within 1pc of about 20 stars and brown dwarfs every Myr, the plausibility of a home system depends also on an appropriate (low) encounter velocity.

The near-Earth object (NEO) population is a window into the original conditions of the protosolar nebula, and has the potential to provide a key pathway for the delivery of water and organics to the early Earth. In addition to delivering the crucial ingredients for life, NEOs can pose a serious hazard to humanity since they can impact the Earth. To properly quantify the impact risk, physical properties of the NEO population need to be studied. Unfortunately, NEOs have a great variation in terms of mitigation-relevant quantities (size, albedo, composition, etc.) and less than 15% of them have been characterized to date. There is an urgent need to undertake a comprehensive characterization of smaller NEOs (D<300m) given that there are many more of them than larger objects. One of the main aims of the NEOShield-2 project (2015--2017), financed by the European Community in the framework of the Horizon 2020 program, is therefore to retrieve physical properties of a wide number of NEOs in order to design impact mitigation missions and assess the consequences of an impact on Earth. We carried out visible photometry of NEOs, making use of the DOLORES instrument at the Telescopio Nazionale Galileo (TNG, La Palma, Spain) in order to derive visible color indexes and the taxonomic classification for each target in our sample. We attributed for the first time the taxonomical complex of 67 objects obtained during the first year of the project. While the majority of our sample belong to the S-complex, carbonaceous C-complex NEOs deserve particular attention. These NEOs can be located in orbits that are challenging from a mitigation point of view, with high inclination and low minimum orbit intersection distance (MOID). In addition, the lack of carbonaceous material we see in the small NEO population might not be due to an observational bias alone.

We present observations of main-belt comet 358P/PANSTARRS (P/2012 T1) obtained using the Gemini South telescope from 2017 July to 2017 December, as the object approached perihelion for the first time since its discovery. We find best-fit IAU phase function parameters of H_R=19.5+/-0.2 mag and G_R=-0.22+/-0.13 for the nucleus, corresponding to an effective radius of r_N=0.32+/-0.03 km (assuming an albedo of p_R=0.05). The object appears significantly brighter (by >1 mag) than expected starting in 2017 November, while a faint dust tail oriented approximately in the antisolar direction is also observed on 2017 December 18. We conclude that 358P has become active again for the first time since its previously observed active period in 2012-2013. These observations make 358P the seventh main-belt comet candidate confirmed to exhibit recurrent activity near perihelion with intervening inactivity away from perihelion, strongly indicating that its activity is sublimation-driven. Fitting a linear function to the ejected dust masses inferred for 358P in 2017 when it is apparently active, we find an average net dust production rate of 2.0+/-0.6 kg/s (assuming a mean effective particle radius of 1 mm) and an estimated activity start date of 2017 November 8+/-4 when the object was at a true anomaly of 316+/-1 deg and a heliocentric distance of R=2.54 AU. Insufficient data is currently available to ascertain whether activity strength has changed between the object's 2012-2013 and 2017 active periods. Further observations are therefore highly encouraged during the object's upcoming observing window (2018 August through 2019 May).

We show that `Oumuamua's excited spin could be in a high energy LAM state, which implies that its shape could be far from the highly elongated shape found in previous studies. CLEAN and ANOVA algorithms are used to analyze `Oumuamua's lightcurve using 818 observations over 29.3~days. Two fundamental periodicities are found at frequencies (2.77$\pm$0.11) and (6.42$\pm$0.18)~cycles/day, corresponding to (8.67$\pm$0.34)~h and (3.74$\pm$0.11)~h, respectively. The phased data show that the lightcurve does not repeat in a simple manner, but approximately shows a double minimum at 2.77~cycles/day and a single minimum at 6.42~cycles/day. This is characteristic of an excited spin state. `Oumuamua could be spinning in either the long (LAM) or short (SAM) axis mode. For both, the long axis precesses around the total angular momentum vector with an average period of (8.67$\pm$0.34)~h. For the three LAMs we have found, the possible rotation periods around the long axis are 6.58, 13.15, or 54.48~h, with 54.48~h being the most likely. `Oumuamua may also be nutating with respective periods of half of these values. We have also found two possible SAM states where `Oumuamua oscillates around the long axis with possible periods at 13.15 and 54.48~h, the latter as the most likely. In this case any nutation will occur with the same periods. Determination of the spin state, the amplitude of the nutation, the direction of the TAMV, and the average total spin period may be possible with a direct model fit to the lightcurve. We find that `Oumuamua is "cigar-shaped"', if close to its lowest rotational energy, and an extremely oblate spheroid if close to its highest energy state for its total angular momentum.

On 2015 November 13, the small artificial object designated WT1190F entered the Earth atmosphere above the Indian Ocean offshore Sri Lanka after being discovered as a possible new asteroid only a few weeks earlier. At ESA's SSA-NEO Coordination Centre we took advantage of this opportunity to organize a ground-based observational campaign, using WT1190F as a test case for a possible similar future event involving a natural asteroidal body.

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.

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.

We imaged Comet 252P/2000 G1 (LINEAR) (hereafter 252P) with the Hubble Space Telescope and both 252P and P/2016 BA$_{14}$ (PanSTARRS) (hereafter BA$_{14}$) with the Discovery Channel Telescope in March and April 2016, surrounding its close encounter to Earth. The r'-band $Af\rho$ of 252P in a 0.2"-radius aperture were $16.8\pm0.3$ and $57\pm1$ cm on March 14 and April 4, respectively, and its gas production rates were: $Q$(OH) = $(5.8\pm0.1)\times10^{27}$ s$^{-1}$, and $Q$(CN) = $(1.25\pm0.01)\times10^{25}$ s$^{-1}$ on April 17. The r'-band upper limit $Af\rho$ of BA1$_{14}$ was $0.19\pm0.01$ cm in a 19.2"-radius aperture, and $Q$(CN) = $(1.4\pm0.1)10^{22}$ s$^{-1}$ on April 17, 2017. 252P shows a bright and narrow jet of a few hundred kilometers long in the sunward direction, changing its projected position angle in the sky with a periodicity consistent with 7.24 hours. However, its photometric lightcurve is consistent with a periodicity of 5.41 hours. We suggest that the nucleus of 252P is likely in a non-principal axis rotation. The nucleus radius of 252P is estimated to be about $0.3\pm0.03$ km, indicating an active fraction of 40% to >100% in its 2016 apparition. Evidence implies a possible cloud of slow-moving grains surrounding the nucleus. The activity level of 252P in the 2016 apparition increased by two orders of magnitude from its previous apparitions, making this apparition unusual. On the other hand, the activity level of BA14 appears to be at least three orders of magnitude lower than that of 252P, despite its ten times or larger surface area.

We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies.

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 present initial time-resolved observations of the split comet 332P/Ikeya-Murakami taken using the Hubble Space Telescope. Our images reveal a dust-bathed cluster of fragments receding from their parent nucleus at projected speeds in the range 0.06 to 3.5 m s$^{-1}$ from which we estimate ejection times from October to December 2015. The number of fragments with effective radii $\gtrsim$20 m follows a differential power law with index $\gamma$ = -3.6$\pm$0.6, while smaller fragments are less abundant than expected from an extrapolation of this power-law. We argue that, in addition to losses due to observational selection, torques from anisotropic outgassing are capable of destroying the small fragments by driving them quickly to rotational instability. Specifically, the spin-up times of fragments $\lesssim$20 m in radius are shorter than the time elapsed since ejection from the parent nucleus. The effective radius of the parent nucleus is $r_e \le$ 275 m (geometric albedo 0.04 assumed). This is about seven times smaller than previous estimates and results in a nucleus mass at least 300 times smaller than previously thought. The mass in solid pieces, $2\times10^9$ kg, is about 4% of the mass of the parent nucleus. As a result of its small size, the parent nucleus also has a short spin-up time. Brightness variations in time-resolved nucleus photometry are consistent with rotational instability playing a role in the release of fragments.

We describe 2016 January to April observations of the fragments of 332P/Ikeya-Murakami, a comet earlier observed in a 2010 October outburst (Ishiguro et al 2014). We present photometry of the fragments, and perform simulations to infer the time of breakup. We argue that the eastern-most rapidly brightening fragment ($F4$) best corresponds to the original nucleus, rather than the initial bright fragment $F1$. We compute radial and tangential non-gravitational parameters, $A_1 = (1.5 \pm 0.4) \times 10^{-8}$ AU day$^{-2}$ and $(7.2 \pm 1.9) \times 10^{-9}$ AU day$^{-2}$; both are consistent with zero at the $4\sigma$ level. Monte Carlo simulations indicate that the fragments were emitted on the outbound journey well after the 2010 outburst, with bright fragment $F1$ splitting in mid--2013 and the fainter fragments within months of the 2016 January recovery. Western fragment $F7$ is the oldest, dating from 2011. We suggest that the delayed onset of the splitting is consistent with a self-propagating crystallization of water ice.

In this work we propose the Amor-type asteroid 2009 WN25 as the likely progenitor of the November i-Draconids (NID, IAU#392), a recently detected weak annual meteoroid stream. We first describe our recovery and follow-up effort to obtain timely ground based astrometry with large aperture telescopes, and ensure that 2009 WN25 would not become lost. We then discuss the possible parent-stream association, using its updated orbit to model the ejection of dust particles from the surface of the parent body and match the observed properties of the stream.

The Mars flyby of C/2013 A1 (Siding Spring) represented a unique opportunity for imaging a long-period comet and resolving its nucleus and rotation period. Because of the small encounter distance and the high relative velocity, the goal of successfully observing C/2013 A1 from the Mars orbiting spacecrafts posed strict accuracy requirements on the comet's ephemerides. These requirements were hard to meet, as comets are known for being highly unpredictable: astrometric observations can be significantly biased and nongravitational perturbations affect comet trajectories. Therefore, even prior to the encounter, we remeasured a couple of hundred astrometric images obtained with ground-based and Earth-orbiting telescopes. We also observed the comet with the Mars Reconnaissance Orbiter's High Resolution Imaging Science Experiment (HiRISE) camera on 2014 October 7. In particular, these HiRISE observations were decisive in securing the trajectory and revealed that out-of-plane nongravitational perturbations were larger than previously assumed. Though the resulting ephemeris predictions for the Mars encounter allowed observations of the comet from the Mars orbiting spacecrafts, post-encounter observations show a discrepancy with the pre-encounter trajectory. We reconcile this discrepancy by employing the Rotating Jet Model, which is a higher fidelity model for nongravitational perturbations and provides an estimate of C/2013 A1's spin pole.

We describe systematic ranging, an orbit determination technique especially suitable to assess the near-term Earth impact hazard posed by newly discovered asteroids. For these late warning cases, the time interval covered by the observations is generally short, perhaps a few hours or even less, which leads to severe degeneracies in the orbit estimation process. The systematic ranging approach gets around these degeneracies by performing a raster scan in the poorly-constrained space of topocentric range and range rate, while the plane of sky position and motion are directly tied to the recorded observations. This scan allows us to identify regions corresponding to collision solutions, as well as potential impact times and locations. From the probability distribution of the observation errors, we obtain a probability distribution in the orbital space and then estimate the probability of an Earth impact. We show how this technique is effective for a number of examples, including 2008 TC3 and 2014 AA, the only two asteroids to date discovered prior to impact.

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.

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. We find that the object is clearly active both in observations obtained in 2014 and in precovery observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggesting that its activity is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in the total brightness of the comet over the 2014 observing period, and dust modeling results showing that the dust emission persists over at least three months during both active periods, where we find start dates for emission no later than 2003 July 24+/-10 for the 2003 active period and 2014 July 28+/-10 for the 2014 active period. From serendipitous observations by the Subaru Telescope in 2004 when the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude of H_R=17.1+/-0.3, corresponding to an effective nucleus radius of r_e~1.00+/-0.15 km. The object's faintness at that time means we cannot rule out the presence of activity, and so this computed radius should be considered an upper limit. We find that 313P's orbit is intrinsically chaotic, having a Lyapunov time of T_l=12000 yr and being located near two 3-body mean-motion resonances with Jupiter and Saturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 T1 (PANSTARRS), to belong to the ~155 Myr old Lixiaohua asteroid family.

Asteroid 2009 FD could impact Earth between 2185 and 2196. The long term propagation to the possible impacts and the intervening planetary encounters make 2009 FD one of the most challenging asteroids in terms of hazard assessment. To compute accurate impact probabilities we model the Yarkovsky effect by using the available physical characterization of 2009 FD and general properties of the Near Earth Asteroid population. We perform the hazard assessment with two independent methods: the first method is a generalization of the standard impact monitoring algorithms in use by NEODyS and Sentry, while the second one is based on a Monte Carlo approach. Both methods generate orbital samples in a 7 dimensional space that includes orbital elements and the parameter characterizing the Yarkovsky effect. The highest impact probability is $2.7 \times 10^{-3}$ for an impact during the 2185 Earth encounter. Impacts after 2185 corresponding to resonant returns are possible, the most relevant being in 2190 with a probability of $3 \times 10^{-4}$. Both numerical methods can be used in the future to handle similar cases. The structure of resonant returns and the list of the possible keyholes on the Target Plane of the scattering encounter in 2185 can be predicted by an analytic theory.

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 our astrometric observations of the small near-Earth object 2011~MD ($H \sim 28.0$), obtained after its very close fly-by to Earth in June 2011. Our set of observations extends the observational arc to $73$ days, and together with the published astrometry obtained around the Earth fly-by allows a direct detection of the effect of radiation pressure on the object, with a confidence of $5\sigma$. The detection can be used to put constraints on the density of the object, pointing to either an unexpectedly low value of $\rho = (640 \pm 330) \mbox{ kg} / \mbox{m} ^3$ ($68\%$ confidence interval) if we assume a typical probability distribution for the unknown albedo, or to an unusually high reflectivity of its surface. This result may have important implications both in terms of impact hazard from small objects and in light of a possible retrieval of this target.

P/2013 P5 PANSTARRS was discovered in Aug. 2013, displaying a cometary tail, but with orbital elements typical for a member of the inner asteroid Main Belt. We monitored the object from 2013 Aug. 30 until Oct. 05 using the CFHT, NTT, CA 1.23m, Perkins 1.8m (Lowell), and the 0.6m TRAPPIST telescopes. We measured its nuclear radius to be r < 0.25-0.29km, and its colours g-r = 0.58+/-0.05 and r-i = 0.23+/-0.06, typical for an S-class asteroid. We failed to detect any rotational light curve, with an amplitude < 0.05mag and a double-peaked rotation period < 20h. A detailed Finson-Probstein analysis of deep NTT and CFHT images indicated that the object was active since at least late January 2013 until the time of the latest observations in 2013 September, with at least two peaks of activity around 2013 June 14+/-10d and 2013 July 22+/-3d. The changes of activity level and the activity peaks were extremely sharp and short, shorter than the temporal resolution of our observations (about 1d). The dust distribution was similar during these two events, with dust grains covering at least the 1-1000{\mu}m range. The total mass ejected in grains <1mm was estimated to be 3.0 10$^6$kg and 2.6 10$^7$kg around the two activity peaks. Rotational disruption cannot be ruled out as the cause of the dust ejection. We also propose that the components of a contact binary might gently rub and produce the observed emission. Volatile sublimation might also explain what appears as cometary activity over a period of 8 months. However, while Main Belt comets best explained by ice sublimation are found in the outskirts of the Main Belt, where water ice is believed to be able to survive buried in moderately large objects for the age of the solar system deeply, the presence of volatiles in an object smaller than 300m in radius would be very surprising in the inner asteroid belt.

We present the results of our search for a dynamical family around the active asteroid P/2012F5 (Gibbs). By applying the hierarchical clustering method, we discover an extremely compact 9-body cluster associated with P/2012F5. The statistical significance of this newly discovered Gibbs cluster is estimated to be >99.9%, strongly suggesting that its members share a common origin. The cluster is located in a dynamically cold region of the outer main-belt at a proper semi-major axis of about 3.005 AU, and all members are found to be dynamically stable over very long time-scales. Backward numerical orbital integrations show that the age of the cluster is only 1.5 $\pm$ 0.1 Myr. Taxonomic classifications are unavailable for most of the cluster members, but SDSS spectrophotometry available for two cluster members indicate that both appear to be $Q$-type objects. We also estimate a lower limit of the size of the parent body to be about 10 km, and find that the impact event which produced the Gibbs cluster is intermediate between a cratering and a catastrophic collision. In addition, we search for new main-belt comets in the region of the Gibbs cluster by observing seven asteroids either belonging to the cluster, or being very close in the space of orbital proper elements. However, we do not detect any convincing evidence of the presence of a tail or coma in any our targets. Finally, we obtain optical images of P/2012F5, and find absolute R-band and V-band magnitudes of $H_R$ = 17.0 $\pm$ 0.1 mag and $H_V$ = 17.4 $\pm$ 0.1 mag, respectively, corresponding to an upper limit on the diameter of the P/2012F5 nucleus of about 2 km.

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.

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.

We assess the risk of an Earth impact for asteroid (99942) Apophis by means of a statistical analysis accounting for the uncertainty of both the orbital solution and the Yarkovsky effect. We select those observations with either rigorous uncertainty information provided by the observer or a high established accuracy. For the Yarkovsky effect we perform a Monte Carlo simulation that fully accounts for the uncertainty in the physical characterization, especially for the unknown spin orientation. By mapping the uncertainty information onto the 2029 b-plane and identifying the keyholes corresponding to subsequent impacts we assess the impact risk for future encounters. In particular, we find an impact probability greater than 10^-6 for an impact in 2068. We analyze the stability of the impact probability with respect to the assumptions on Apophis' physical characterization and consider the possible effect of the early 2013 radar apparition.

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 report the direct detection of radiation pressure on the asteroid 2009 BD, one of the smallest multi-opposition near-Earth objects currently known, with H ~ 28.4. Under the purely gravitational model of NEODyS the object is currently considered a possible future impactor, with impact solutions starting in 2071. The detection of a radiation-related acceleration allows us to estimate an Area to Mass Ratio (AMR) for the object, that can be converted (under some assumptions) into a range of possible values for its average density. Our result AMR = (2.97 \pm 0.33) x 10^(-4) m^2 kg^(-1) is compatible with the object being of natural origin, and it is narrow enough to exclude a man-made nature. The possible origin of this object, its future observability, and the importance of radiation pressure in the impact monitoring process, are also discussed.

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.

The asteroid 2003 EH1, proposed as the parent body of the Quadrantid meteor shower, is thought to be the remnant of a past cometary object, tentatively identified with the historical comets C/1490 Y1 and C/1385 U1. In the present work we use recovery astrometry to extend the observed arc of 2003 EH1 from 10 months to about 5 years, enough to exclude the proposed direct relationship of the asteroid with both of the comets.