We present observations of comet 67P/Churyumov-Gerasimenko acquired in support of the $Rosetta$ mission. We obtained usable data on 68 nights from 2014 September until 2016 May, with data acquired regularly whenever the comet was observable. We collected an extensive set of near-IR $J$, $H$, and $Ks$ data throughout the apparition plus visible-light images in $g'$, $r'$, $i'$, and $z'$ when the comet was fainter. We also obtained broadband $R$ and narrowband $CN$ filter observations when the comet was brightest using telescopes at Lowell Observatory. The appearance was dominated by a central condensation and the tail until 2015 June. From 2015 August onwards there were clear asymmetries in the coma, which enhancements revealed to be due to the presence of up to three features (i.e., jets). The features were similar in all broadband filters; $CN$ images did not show these features but were instead broadly enhanced in the southeastern hemisphere. Modeling using the parameters from Vincent et al. (2013) replicated the dust morphology reasonably well, indicating that the pole orientation and locations of active areas have been relatively unchanged over at least the last three apparitions. The dust production, as measured by $A(0^{\circ})f{\rho}$ peaked $\sim$30 days after perihelion and was consistent with predictions from previous apparitions. $A(0^{\circ})f{\rho}$ as a function of heliocentric distance was well fit by a power-law with slope $-$4.2 from 35-120 days post-perihelion. We detected photometric evidence of apparent outbursts on 2015 August 22 and 2015 September 19, although neither was discernible morphologically in this dataset.

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.

Around the time of its perihelion passage the observability of 67P/Churyumov-Gerasimenko from Earth was limited to very short windows each morning from any given site, due to the low solar elongation of the comet. The peak in the comet's activity was therefore difficult to observe with conventionally scheduled telescopes, but was possible where service/queue scheduled mode was possible, and with robotic telescopes. We describe the robotic observations that allowed us to measure the total activity of the comet around perihelion, via photometry (dust) and spectroscopy (gas), and compare these results with the measurements at this time by Rosetta's instruments. The peak of activity occurred approximately two weeks after perihelion. The total brightness (dust) largely followed the predictions from Snodgrass et al. 2013, with no significant change in total activity levels from previous apparitions. The CN gas production rate matched previous orbits near perihelion, but appeared to be relatively low later in the year.

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

As the ESA Rosetta mission approached, orbited, and sent a lander to comet 67P/Churyumov-Gerasimenko in 2014, a large campaign of ground-based observations also followed the comet. We constrain the total activity level of the comet by photometry and spectroscopy to place Rosetta results in context and to understand the large-scale structure of the comet's coma pre-perihelion. We performed observations using a number of telescopes, but concentrate on results from the 8m VLT and Gemini South telescopes in Chile. We use R-band imaging to measure the dust coma contribution to the comet's brightness and UV-visible spectroscopy to search for gas emissions, primarily using VLT/FORS. In addition we imaged the comet in near-infrared wavelengths (JHK) in late 2014 with Gemini-S/Flamingos 2. We find that the comet was already active in early 2014 at heliocentric distances beyond 4 au. The evolution of the total activity (measured by dust) followed previous predictions. No gas emissions were detected despite sensitive searches. The comet maintains a similar level of activity from orbit to orbit, and is in that sense predictable, meaning that Rosetta results correspond to typical behaviour for this comet. The gas production (for CN at least) is highly asymmetric with respect to perihelion, as our upper limits are below the measured production rates for similar distances post-perihelion in previous orbits.

Here we present observations of 7 large Kuiper Belt Objects. From these observations, we extract a point source catalog with $\sim0.01"$ precision, and astrometry of our target Kuiper Belt Objects with $0.04-0.08"$ precision within that catalog. We have developed a new technique to predict the future occurrence of stellar occultations by Kuiper Belt Objects. The technique makes use of a maximum likelihood approach which determines the best-fit adjustment to cataloged orbital elements of an object. Using simulations of a theoretical object, we discuss the merits and weaknesses of this technique compared to the commonly adopted ephemeris offset approach. We demonstrate that both methods suffer from separate weaknesses, and thus, together provide a fair assessment of the true uncertainty in a particular prediction. We present occultation predictions made by both methods for the 7 tracked objects, with dates as late as 2015. Finally, we discuss observations of three separate close passages of Quaoar to field stars, which reveal the accuracy of the element adjustment approach, and which also demonstrate the necessity of considering the uncertainty in stellar position when assessing potential occultations.

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.

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 present three new eclipsing white-dwarf / M-dwarf binary systems discovered during a search for transiting planets around M-dwarfs. Unlike most known eclipsing systems of this type, the optical and infrared emission is dominated by the M-dwarf components, and the systems have optical colors and discovery light curves consistent with being Jupiter-radius transiting planets around early M-dwarfs. We detail the PTF/M-dwarf transiting planet survey, part of the Palomar Transient Factory (PTF). We present a Graphics Processing Unit (GPU)-based box-least-squares search for transits that runs approximately 8X faster than similar algorithms implemented on general purpose systems. For the discovered systems, we decompose low-resolution spectra of the systems into white-dwarf and M-dwarf components, and use radial velocity measurements and cooling models to estimate masses and radii for the white dwarfs. The systems are compact, with periods between 0.35 and 0.45 days and semimajor axes of approximately 2 solar radii (0.01 AU). We use the Robo-AO laser guide star adaptive optics system to tentatively identify one of the objects as a triple system. We also use high-cadence photometry to put an upper limit on the white dwarf radius of 0.025 solar radii (95% confidence) in one of the systems. We estimate that 0.08% (90% confidence) of M-dwarfs are in these short-period, post-common-envelope white-dwarf / M-dwarf binaries where the optical light is dominated by the M-dwarf. Similar eclipsing binary systems can have arbitrarily small eclipse depths in red bands and generate plausible small-planet-transit light curves. As such, these systems are a source of false positives for M-dwarf transiting planet searches. We present several ways to rapidly distinguish these binaries from transiting planet systems.

We report the discovery of a new transiting planet in the Southern Hemisphere. It has been found by the WASP-south transit survey and confirmed photometrically and spectroscopically by the 1.2m Swiss Euler telescope, LCOGT 2m Faulkes South Telescope, the 60 cm TRAPPIST telescope and the ESO 3.6m telescope. The orbital period of the planet is 2.94 days. We find it is a gas giant with a mass of 0.88 \pm 0.10 Mj and a radius estimated at 0.96 \pm 0.05 Rj . We have also obtained spectra during transit with the HARPS spectrograph and detect the Rossiter-McLaughlin effect despite its small amplitude. Because of the low signal to noise of the effect and of a small impact parameter we cannot place a constraint on the projected spin-orbit angle. We find two confiicting values for the stellar rotation. Our determination, via spectral line broadening gives v sin I = 2.2 \pm 0.3 km/s, while another method, based on the activity level using the index log R'HK, gives an equatorial rotation velocity of only v = 1.35 \pm 0.20 km/s. Using these as priors in our analysis, the planet could either be misaligned or aligned. This should send strong warnings regarding the use of such priors. There is no evidence for eccentricity nor of any radial velocity drift with time.

PTF/M-dwarfs is a 100,000-target M-dwarf planetary transit survey, a Key Project of the Palomar Transient Factory (PTF) collaboration. The survey is sensitive to Jupiter-radius planets around all of the target stars, and has sufficient precision to reach Neptunes and super-Earths for the best targets. The Palomar Transient Factory is a fully-automated, wide-field survey aimed at a systematic exploration of the optical transient sky. The survey is performed using a new 7.26 square degree camera installed on the 48 inch Samuel Oschin telescope at Palomar Observatory. Each 92-megapixel R-band exposure contains about 3,000 M-dwarfs usable for planet detection. In each PTF observational season PTF/M-dwarfs searches for Jupiter-radius planets around almost 30,000 M-dwarfs, Neptune-radius planets around approximately 500 M-dwarfs, and super-Earths around 100 targets. The full survey is expected to cover more than 100,000 targets over the next several years. Photometric and spectroscopic followup operations are performed on the Palomar 60-inch, LCOGT, Palomar 200-inch, MDM and Keck telescopes. The survey has been running since mid-2009. We detail the survey design, the survey's data analysis pipeline and the performance of the first year of operations.