If a transiting exoplanet has a moon, that moon could be detected directly from the transit itproduces itself, or indirectly via the transit timing variations it produces in its parent planet. There is a range of parameter space where the Kepler Space Telescope is sensitive to the TTVs exomoons might produce, though the moons themselves would be too small to detect photometrically via their own transits. The Earth's Moon, for example, produces TTVs of 2.6 minutes amplitude by causing our planet to move around their mutual center of mass. This ismore than Kepler's short-cadence interval of 1 minute and so nominally detectable (if transit timings are of comparable accuracy), even though the Moon's transit signatureis only 7% that of Earth's, well below Kepler's nominal threshold. Here we explore eight systems from the Kepler data set to examine the exomoon hypothesisas an explanation for their transit timing variations, which we compare with the alternatehypothesis that the TTVs are caused by an non-transiting planet in the system. We find that the TTVs of six of these systems could be plausibly explained by an exomoon, the size of which would not be nominally detectable by Kepler. Though we also find that the TTVsc ould be equally well reproduced by the presence of a non-transiting planet in the system, the observations are nevertheless completely consistent with a existence of a dynamically stablemoon small enough to fall below Kepler's photometric threshold for transit detection, and these systems warrant further observation and analysis.

We present a multi-instrumental, multidecadal analysis of the activity of the Eta-Aquariid and Orionid meteor showers for the purpose of constraining models of 1P/Halley's meteoroid streams. The interannual variability of the showers' peak activity and period of duration is investigated through the compilation of published visual and radar observations prior to 1985 and more recent measurements reported in the International Meteor Organization (IMO) Visual Meteor DataBase, by the IMO Video Meteor Network and by the Canadian Meteor Orbit Radar (CMOR). These techniques probe the range of meteoroid masses from submilligrams to grams. The Eta-Aquariids and Orionids activity duration, shape, maximum zenithal hourly rates (ZHR) values, and the solar longitude of annual peaks since 1985 are analyzed. When available, annual activity profiles recorded by each detection network were measured and are compared. Observations from the three detection methods show generally good agreement in the showers' shape, activity levels, and annual intensity variations. Both showers display several activity peaks of variable location and strength with time. The Eta-Aquariids are usually two to three times stronger than the Orionids, but the two showers display occasional outbursts with peaks two to four times their usual activity level. CMOR observations since 2002 seem to support the existence of an ~12 year cycle in Orionids activity variations; however, additional and longer term radar and optical observations of the shower are required to confirm such periodicity.

It has recently been shown by Egal et al. (2017) that some types of existing meteor in-atmosphere trajectory estimation methods may be less accurate than others, particularly when applied to high precision optical measurements. The comparative performance of trajectory solution methods has previously only been examined for a small number of cases. Besides the radiant, orbital accuracy depends on the estimation of pre-atmosphere velocities, which have both random and systematic biases. Thus it is critical to understand the uncertainty in velocity measurement inherent to each trajectory estimation method. In this first of a series of two papers, we introduce a novel meteor trajectory estimation method which uses the observed dynamics of meteors across stations as a global optimization function and which does not require either a theoretical or empirical flight model to solve for velocity. We also develop a 3D observational meteor trajectory simulator that uses a meteor ablation model to replicate the dynamics of meteoroid flight, as a means to validate different trajectory solvers. We both test this new method and compare it to other methods, using synthetic meteors from three major showers spanning a wide range of velocities and geometries (Draconids, Geminids, Perseids). We determine which meteor trajectory solving algorithm performs better for: all-sky, moderate field of view, and high-precision narrow-field optical meteor detection systems. The results are presented in the second paper in this series. Finally, we give detailed equations for estimating meteor trajectories and analytically computing meteoroid orbits, and provide the Python code of the methodology as open source software.

The DART spacecraft is planned to impact the secondary of the binary asteroid (65803) Didymos in 2022, to assess deflection strategies for planetary defense. The impact will create a crater and release asteroidal material, some of which will escape the Didymos system. Because the closest point of approach of Didymos to Earth's orbit is only 6 million km (about 16 times the Earth-Moon distance), some ejected material will make its way sooner or later to our planet, and the observation of these particles as meteors would increase the scientific payout of the DART mission. The DART project may also represent the first human-generated meteoroids to reach Earth, and a test case for human activity on asteroids and its eventual contribution to the meteoroid environment and spacecraft impact risk. This study examines the amount and timing of the delivery of meteoroids from Didymos to near-Earth space. Very little DART-ejected material will reach our planet, and most only after thousands of years. But some material, the smallest particles and/or those ejected at the highest velocities could be delivered to Earth-crossing trajectories almost immediately, though at very low fluxes. Timing and radiant directions for material reaching the Earth are calculated, though the detection of substantial numbers would indicate more abundant and/or faster ejecta than is expected. The DART impact will also create a new meteoroid stream, though probably not a very dense one. However, larger, more capable asteroid impactors could create meteoroid streams in which the particle flux exceeds that naturally occurring in the Solar System, with implications for spacecraft safety.

The low velocity of interstellar asteroid 1I/'Oumuamua with respect to our galaxy's Local Standard of Rest implies it is young. Adopting the young age hypothesis, we assess possible origin systems for this interstellar asteroid and for 2I/Borisov, though the latter's higher speed means it is unlikely to be young. First, their past trajectories are modelled under gravitational scattering by galactic components ('disk heating') to assess how far back one can trace them. The stochastic nature of disk heating means that a back-integration can only expect to be accurate to within 15 pc and 2 kms$^{-1}$ at -10 Myr, dropping steeply to 400 pc and 10 kms$^{-1}$ at -100 Myr, sharply limiting our ability to determine a precise origin. Nevertheless, we show `Oumuamua's origin system is likely currently within 1 kpc of Earth, in the local Orion Arm. Second, we back-integrate 'Oumuamua's trajectory to assess source regions, emphasizing young systems and moving groups. Though disk heating allows for only a statistical link to source regions, `Oumuamua passed through a considerable subset of the Carina and Columba moving groups when those groups were forming. This makes them perhaps the most plausible source region, if 'Oumuamua was ejected during planet formation or via intra-cluster interactions. We find three stars in the Ursa Major group, one brown dwarf, and seven other stars to have plausible encounters with 2I/Borisov, within 2 pc and 30 kms$^{-1}$. These encounters' high relative speeds mean none are likely to be the home of 2I/Borisov.

Granvik et al. (2016) report an absence of asteroids on orbits with perihelia near the Sun that they attribute to the 'supercatastrophic disruption' of these bodies. Here we investigate whether there is evidence for this process among other bodies with similarly low perihelia: near-Earth asteroids, SOHO comets, as well as meter-sized and millimeter-sized meteoroids. We determine no known near-Earth asteroids have past (last 10^4 years) histories residing significantly inside the Granvik et al. (2016) limit, indirectly supporting the disruption hypothesis. The exception is asteroid (467372) 2004 LG which spent 2500 years within this limit, and thus presents a challenge to that theory. Phaethon has a perihelion distance hovering just above the limit and may be undergoing slow disruption, which may be the source of its dust complex. We find that the rate at which ungrouped SOHO comets are observed is consistent with expected rates for the injection of small (25 m) class asteroids into the near-Sun region and suggest that this fraction of the SOHO-observed comet population may in fact be asteroidal in origin. We also find that there is an absence of meter-sized bodies with near-Sun perihelia but an excess of millimeter-sized meteoroids. This implies that if near-Sun asteroids disrupt, they do not simply fragment into meter-sized chunks but disintegrate ultimately into millimeter-sized particles. We propose that the disruption of near-Sun asteroids as well as the anomalous brightening and destruction processes that affect SOHO comets occur through meteoroid erosion, that is, the removal of material through impacts by high-speed near-Sun meteoroids.

In the first paper of this series we examined existing methods of optical meteor trajectory estimation and developed a novel method which simultaneously uses both the geometry and the dynamics of meteors to constrain their trajectories. We also developed a simulator which uses an ablation model to generate realistic synthetic meteor trajectories which we use to test meteor trajectory solvers. In this second paper, we perform simulation validation to estimate radiant and velocity accuracy which may be achieved by various meteor observation systems as applied to several meteor showers. For low-resolution all-sky systems, where the meteor deceleration is generally not measurable, the multi-parameter fit method assuming a constant velocity better reproduces the radiant and speed of synthetic meteors. For moderate field of view systems, our novel method performs the best at all convergence angles, while multi-parameter fit methods generally produce larger speed errors. For high-resolution, narrow field of view systems, we find our new method of trajectory estimation reproduces radiant and speed more accurately than all other methods tested. The ablation properties of meteoroids are commonly found to be the limiting factor in velocity accuracy. We show that the true radiant dispersion of meteor showers can be reliably measured with moderate field of view (or more precise) systems provided appropriate methods of meteor trajectory estimation are employed. Finally, we compare estimated and real angular radiant uncertainty and show that for the solvers tested the real radiant error is on average underestimated by a factor of two.

There is a vibrant and effective planetary science community in Canada. We do research in the areas of meteoritics, asteroid and trans-Neptunian object orbits and compositions, and space weather, and are involved in space probe missions to study planetary surfaces and interiors. For Canadian planetary scientists to deliver the highest scientific impact possible, we have several recommendations. Our top recommendation is to join LSST and gain access to the full data releases by hosting a data centre, which could be done by adding to the CADC, which is already highly involved in hosting planetary data and supporting computational modelling for orbital studies. We also support MSE, which can provide spectroscopy and thus compositional information for thousands of small bodies. We support a Canadian-led microsatellite, POEP, which will provide small body sizes by measuring occultations. We support the idea of piggybacking space weather instruments on other astronomical space probes to provide data for the space weather community. Many Canadian planetary scientists are involved in space probe missions, but through haphazard and temporary arrangements like co-appointments at US institutions, so we would like the community to support Canadian researchers to participate in these large, international missions.

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 report on a multi-observatory campaign to examine asteroid 3200 Phaethon during its December 2017 close approach to Earth, in order to improve our measurements of its fundamental parameters, and to search for surface variations, cometary activity and fragmentation. The mean colors of Phaethon are B-V = 0.702 +/- 0.004, V-R = 0.309 +/- 0.003, R-I = 0.266 +/- 0.004, neutral to slightly blue, consistent with previous classifications of Phaethon as a F-type or B-type asteroid. Variations in Phaethon's B-V colors (but not V-R or R-I) with observer sub-latitude are seen and may be associated with craters observed by the Arecibo radar. High cadence photometry over phases from 20 to 100 degrees allows a fit to the values of the HG photometric parameters; H = 14.57 +/- 0.02, 13.63 +/- 0.02, 13.28 +/- 0.02, 13.07 +/- 0.02; G = 0.00 +/- 0.01, -0.09 +/- 0.01, -0.10 +/- 0.01, -0.08 +/- 0.01 in the BVRI filters respectively; the negative G values are consistent with other observations of F type asteroids. Light curve variations were seen that are also consistent with concavities reported by Arecibo, indicative of large craters on Phaethon's surface whose ejecta may be the source of the Geminid meteoroid stream. A search for gas/dust production set an upper limit of 0.06 +/- 0.02 kg/s when Phaethon was 1.449 AU from the Sun, and 0.2 +/- 0.1 kg/s at 1.067 AU. A search for meter-class fragments accompanying Phaethon did not find any whose on-sky motion was not also consistent with background main belt asteroids.

In June 2019 the Earth will approach within 5{\deg} mean anomaly of the centre of the Taurid resonant swarm, its closest post-perihelion encounter with Earth since 1975. This will be the best viewing geometry to detect and place limits on the number of NEOs proposed to reside at the swarm centre until the early 2030s. We present an analysis of the optimal times and pointing locations to image NEOs associated with the swarm

This work presents numerical simulations of meteoroid streams released by comet 21P/Giacobini-Zinner over the period 1850-2030. The initial methodology, based on Vaubaillon et al. (2005), has been updated and modified to account for the evolution of the comet's dust production along its orbit. The peak time, intensity, and duration of the shower were assessed using simulated activity profiles that are calibrated to match observations of historic Draconid outbursts. The characteristics of all the main apparitions of the shower are reproduced, with a peak time accuracy of half an hour and an intensity estimate correct to within a factor of 2 (visual showers) or 3 (radio outbursts). Our model also revealed the existence of a previously unreported strong radio outburst on October 9 1999, that has since been confirmed by archival radar measurements. The first results of the model, presented in Egal et al. (2018), provided one of the best predictions of the recent 2018 outburst. Three future radio outbursts are predicted in the next decade, in 2019, 2025 and 2029. The strongest activity is expected in 2025 when the Earth encounters the young 2012 trail. Because of the dynamical uncertainties associated with comet 21P's orbital evolution between the 1959 and 1965 apparitions, observations of the 2019 radio outburst would be particularly helpful to improve the confidence of subsequent forecasts.

The Kepler Space Telescope has discovered thousands of planets via the transit method. The transit timing variations of these planets allows us not only to infer the existence of other planets, transiting or not, but to characterize a number of parameters of the system. Using the transit timing variations of the planets Kepler-159b and 159c, the transit simulator TTVFast, and the Bayesian Inference tool MultiNest, we predict a new non-transiting planet, Kepler-159d, in a resonant 2:1 orbit with Kepler-159c. This configuration is dynamically stable on at least 10 Myr time scales, though we note that other less stable, higher-order resonances could also produce similar TTVs during the three-year window Kepler was in operation.

Predictions of the 2018 Draconid activity at the Earth and the Sun-Earth L1 and L2 Lagrange points are presented. Numerical simulations of the meteoroids' ejection and evolution from comet 21P/Giacobini-Zinner are performed with a careful implementation of the results analysis and weighting. Model meteoroid fluxes at Earth are derived using as calibration the main peak date, intensity, and shower profiles of previous Draconid outbursts. Good agreement between the model and measurements is found for the 1933, 1946, 1998 and 2011 showers for a meteoroid size distribution index at ejection of about 2.6. A less accurate estimate of the peak time for the 1985, 2005 and 2012 predominantly radio-observed outbursts was found by considering the contribution of individual ejection epochs, while the model peak flux estimate was found to agree with observations to within a factor 3. Despite the promising geometrical configuration in 2018, our simulations predict low Draconid activity is expected on Earth, with a maximum of less than a few tens of meteors per hour around midnight the 9th of October, confirming previous models. At the L1 and L2 Lagrange points, however, the flux estimates suggest a 'meteoroid storm'. The Gaia spacecraft at the L2 region might be able to detect small (~ {\mu}g) Draconid meteoroid impacts centered in a two-hour window around 18h30 UT on the 8th of October, 2018.

Near-Earth asteroid (3200) Phaethon is notable for its association to a strong annual meteor shower, the Geminids, indicative of one or more episodes of mass ejection in the past. The mechanism of Phaethon's past activity is not yet understood. Here we present a Hubble Space Telescope (HST) search of meter-sized fragments in the vicinity of Phaethon, carried out during Phaethon's historic approach to the Earth in mid-December of 2017. Numerical simulation conducted to guide HST's pointing also show that the dynamical evolution of Phaethon-originated particles is quick, as ejected materials take no longer than $\sim250$ yr to spread to the entire orbit of Phaethon. Our search was completed down to 4-meter-class limit (assuming Phaethon-like albedo) and was expected to detect 0.035% particles ejected by Phaethon in the last several decades. The negative result of our search capped the total mass loss of Phaethon over the past few dozen orbits to be $10^{12}$ kg at $3\sigma$ level, taking the best estimates of size power-law from meteor observations and spacecraft data. Our result also implies a millimeter-sized dust flux of $<10^{-12} \mathrm{m^{-2} s^{-1}}$ within 0.1 au of Phaethon, suggesting that any Phaethon-bound mission is unlikely to encounter dense dust clouds.

Observations have revealed rich structures in protoplanetary disks, offering clues about their embedded planets. Due to the complexities introduced by the abundance of gas in these disks, modeling their structure in detail is computationally intensive, requiring complex hydrodynamic codes and substantial computing power. It would be advantageous if computationally simpler models could provide some preliminary information on these disks. Here we apply a particle-only model (that we developed for gas-poor debris disks) to the gas-rich disk, HL Tauri, to address the question of whether such simple models can inform the study of these systems. Assuming three potentially embedded planets, we match HL Tau's radial profile fairly well and derive best-fit planetary masses and orbital radii (0.40, 0.02, 0.21 Jupiter masses for the planets orbiting a 0.55 solar-mass star at 11.22, 29.67, 64.23 AU). Our derived parameters are comparable to those estimated by others, except for the mass of the second planet. Our simulations also reproduce some narrower gaps seen in the ALMA image away from the orbits of the planets. The nature of these gaps is debated but, based on our simulations, we argue they could result from planet-disk interactions via mean-motion resonances, and need not contain planets. Our results suggest that a simple particle-only model can be used as a first step to understanding dynamical structures in gas disks, particularly those formed by planets, and determine some parameters of their hidden planets, serving as useful initial inputs to hydrodynamic models which are needed to investigate disk and planet properties more thoroughly.

Jupiter-family Comet D/1770 L1 (Lexell) was the first discovered Near-Earth Object (NEO), and passed the Earth on 1770 Jul 1 at a recorded distance of 0.015 au. The comet was subsequently lost due to unfavorable observing circumstances during its next apparition followed by a close encounter with Jupiter in 1779. Since then, the fate of D/Lexell has attracted interest from the scientific community, and now we revisit this long-standing question. We investigate the dynamical evolution of D/Lexell based on a set of orbits recalculated using the observations made by Charles Messier, the comet's discoverer, and find that there is a $98\%$ chance that D/Lexell remains in the Solar System by the year of 2000. This finding remains valid even if a moderate non-gravitational effect is imposed. Messier's observations also suggest that the comet is one of the largest known near-Earth comets, with a nucleus of $\gtrsim 10$ km in diameter. This implies that the comet should have been detected by contemporary NEO surveys regardless of its activity level if it has remained in the inner Solar System. We identify asteroid 2010 JL$_{33}$ as a possible descendant of D/Lexell, with a $0.8\%$ probability of chance alignment, but a direct orbital linkage of the two bodies has not been successfully accomplished. We also use the recalculated orbit to investigate the meteors potentially originating from D/Lexell. While no associated meteors have been unambiguously detected, we show that meteor observations can be used to better constrain the orbit of D/Lexell despite the comet being long lost.

Structures observed in debris disks may be caused by gravitational interaction with planetary or stellar companions. These perturbed disks are often thought to indicate the presence of planets and offer insights into the properties of both the disk and the perturbing planets. Gaps in debris disks may indicate a planet physically present within the gap, but such gaps can also occur away from the planet's orbit at mean-motion resonances (MMRs), and this is the focus of our interest here. We extend our study of planet-disk interaction through MMRs, presented in an earlier paper, to systems in which the perturbing planet has moderate orbital eccentricity, a common occurrence in exoplanetary systems. In particular, a new result is that the 3:1 MMR becomes distinct at higher eccentricity, while its effects are absent for circular planetary orbits. We also only consider gravitational interaction with a planetary body of at least 1 Jupiter mass. Our earlier work shows that even a 1 Earth mass planet can theoretically open an MMR gap; however, given the narrow gap that can be opened by a low-mass planet, its observability would be questionable. We find that the widths, locations, and shapes of two prominent structures, the 2:1 and 3:1 MMRs, could be used to determine the mass, semimajor axis, and eccentricity of the planetary perturber and present an algorithm for doing so. These MMR structures can be used to narrow the position and even determine the planetary properties (such as mass) of any inferred but as-yet-unseen planets within a debris disk. We also briefly discuss the implications of eccentric disks on brightness asymmetries and their dependence on the wavelengths with which these disks are observed.

We describe how meteors recorded at the Earth can be used to partly reconstruct the meteoroid environments of the planets if a large sample (i.e. millions of orbits at a minimum) is available. The process involves selecting from the Earth-based sample those meteors which passed near a planet's orbit prior to arriving at the Earth, and so carry information about the planetary meteoroid environment. Indeed this process can be extended to any location in the Solar System, though some regions of space are better sampled than others. From such a reconstruction performed with data from the Canadian Meteor Orbit Radar (CMOR), we reveal that Mars has apex, helion, anti-helion and toroidal sporadic sources, much as Earth does. Such reconstructions, albeit partial, have the potential to provide a wealth of detail about planetary meteoroid environments, and to allow for the ground-truthing of model meteoroid populations without in situ sampling.

We examine the conditions under which material from the martian moons Phobos and Deimos could reach our planet in the form of meteorites. We find that the necessary ejection speeds from these moons (900 and 600 m/s for Phobos and Deimos respectively) are much smaller than from Mars' surface (5000 m/s). These speeds are below typical impact speeds for asteroids and comets (10-40 km/s) at Mars' orbit, and we conclude that the delivery of meteorites from Phobos and Deimos to the Earth can occur.

Centaurs are objects whose orbits are found between those of the giant planets. They are supposed to originate mainly from the TransNeptunian objects, and they are among the sources of NearEarth Objects.TransNeptunian Objects (TNOs) cross Neptune's orbit and produce the Centaurs. We investigate their interactions with main belt asteroids to determine if chaotic scattering caused by close encounters and impacts by these bodies may have played a role in the dynamical evolution of the main belt. We find that Centaurs and TNOs that reach the inner Solar System can modify the orbits of main belt asteroids, though only if their mass is of the order of 1 billion less the mass of the sun for single encounters or, one order less in the case of multiple close encounters. Centaurs and TNOs are unlikely to have significantly dispersed young asteroid families in the main belt, but they can have perturbed some old asteroid families. Current main belt asteroids that originated as Centaurs or Trans-Neptunian Objects may lie in the outer belt with short lifetime less than 4 My, most likely between 2.8 au and 3.2 au at larger eccentricities than typical of main belt asteroids.

Encke-type comet 332P/Ikeya-Murakami is experiencing cascading fragmentation events during its 2016 apparition. It is likely the first splitting Encke-type comet ever being observed. A nongravitational solution to the astrometry reveals a statistical detection of the radial and transverse nongravitational parameters, $A_{1} = \left(1.54 \pm 0.39\right) \times 10^{-8}$ AU day$^{-2}$, and $A_{2} = \left(7.19 \pm 1.92\right) \times 10^{-9}$ AU day$^{-2}$, respectively, which implies a nucleus erosion rate of $\left(0.91 \pm 0.17\right)$% per orbital revolution. The mass-loss rate likely has to be supported by a much larger fraction of an active surface area than known cases of short-period comets; it may be relevant to the ongoing fragmentation. We failed to detect any serendipitous pre-discovery observations of the comet in archival data from major sky surveys, whereby we infer that 332P used to be largely inactive, and is perhaps among the few short-period comets which have been reactivated from weakly active or dormant states. We therefore constrain an upper limit to the nucleus size as $2.0 \pm 0.2$ km in radius. A search for small bodies in similar orbits to that of 332P reveals comet P/2010 B2 (WISE) as the best candidate. From an empirical generalised Jupiter-family (Encke-type included) comet population model, we estimate the likelihood of chance alignment of the 332P--P/2010 B2 pair to be 1 in 33, a small number indicative of a genetic linkage between the two comets on a statistical basis. The pair possibly originated from a common progenitor which underwent a disintegration event well before the twentieth century.

The gravitational influence of a planet on a nearby disk provides a powerful tool for detecting and studying extrasolar planetary systems. Here we demonstrate that gaps can be opened in dynamically cold debris disks at the mean-motion resonances of an orbiting planet. The gaps are opened away from the orbit of the planet itself, revealing that not all disk gaps need contain a planetary body. These gaps are large and deep enough to be detectable in resolved disk images for a wide range of reasonable disk-planet parameters, though we are not aware of any such gaps detected to date. The gap shape and size are diagnostic of the planet location, eccentricity and mass, and allow one to infer the existence of unseen planets, as well as many important parameters of both seen and unseen planets in these systems. We present expressions to allow the planetary mass to be calculated from observed gap width and location.

Temporarily Captured Orbiters (TCOs) are Near-Earth Objects (NEOs) which make a few orbits of Earth before returning to heliocentric orbits. Only one TCO has been observed to date, 2006 RH120, captured by Earth for one year before escaping. Detailed modeling predicts capture should occur from the NEO population predominantly through the Sun-Earth L1 and L2 points, with 1% of TCOs impacting Earth and approximately 0.1% of meteoroids being TCOs. Although thousands of meteoroid orbits have been measured, none until now have conclusively exhibited TCO behaviour, largely due to difficulties in measuring initial meteoroid speed with sufficient precision. We report on a precise meteor observation of January 13, 2014 by a new generation of all-sky fireball digital camera systems operated in the Czech Republic as part of the European Fireball Network, providing the lowest natural object entry speed observed in decades long monitoring by networks world-wide. Modeling atmospheric deceleration and fragmentation yields an initial mass of ~5 kg and diameter of 15 cm, with a maximum Earth-relative velocity just over 11.0 km/s. Spectral observations prove its natural origin. Back-integration across observational uncertainties yields a 92 - 98% probability of TCO behaviour, with close lunar dynamical interaction. The capture duration varies across observational uncertainties from 48 days to 5+ years. We also report on two low-speed impacts recorded by US Government sensors, and we examine Prairie Network event PN39078 from 1965 having an extremely low entry speed of 10.9 km/s. In these cases uncertainties in measurement and origin make TCO designation uncertain.

Previous studies have revealed Jupiter-family comet 252P/LINEAR as a comet that was recently transported into the near-Earth object (NEO) region in $\sim1800$~AD yet only being weakly active. In this Letter, we examine the "formed (almost) dead" hypothesis for 252P/LINEAR using both dynamical and observational approaches. By statistically examining the dynamical evolution of 252P/LINEAR over a period of $10^7$~years, we find the median elapsed residency in the NEO region to be $4\times10^2$~years which highlights the likelihood of 252P/LINEAR as an (almost) first-time NEO. With available cometary and meteor observations, we find the dust production rate of 252P/LINEAR to be at the order of $10^6$~kg per orbit since its entry to the NEO region. These two lines of evidence support the hypothesis that the comet was likely to have formed in a volatile-poor environment. Cometary and meteor observations during the comet's unprecedented close approach to the Earth around 2016 Mar. 21 would be useful for the understanding of the surface and evolutionary properties of this unique comet.

We have analysed the orbits and ablation characteristics in the atmosphere of 59 earth-impacting fireballs, produced by meteoroids one meter in diameter or larger, described here as meter-scale. Using heights at peak luminosity as a proxy for strength, we determine that there is roughly an order of magnitude spread in strengths of the population of meter-scale impactors at the Earth. We use fireballs producing recovered meteorites and well documented fireballs from ground-based camera networks to calibrate our ablation model interpretation of the observed peak height of luminosity as a function of speed. The orbits and physical strength of these objects are consistent with the majority being asteroidal bodies originating from the inner main asteroid belt. We find a lower limit of ~10-15% of our objects have a possible cometary (Jupiter-Family comet and/or Halley-type comet) origin based on orbital characteristics alone. Only half this number, however, also show evidence for weaker than average structure. Two events, Sumava and USG 20131121, have exceptionally high (relative to the remainder of the population) heights of peak brightness. These are physically most consistent with high microporosity objects. We also find three events, including the Oct 8, 2009 airburst near Sulawesi, Indonesia, which display comparatively low heights of peak brightness, consistent with strong monolithic stones or iron meteoroids. Based on orbital similarity, we find a probable connection among several events in our population with the Taurid meteoroid complex; no other major meteoroid streams show probable linkages to the orbits of our meter-scale population. Our impactors cover almost four orders of magnitude in mass, but no trend in height of peak brightness as a function of mass is evident, suggesting no strong trend in strength with size for meter-scale impactors.

It is speculated that some weakly active comets may be transitional objects between active and dormant comets. These objects are at a unique stage of the evolution of cometary nuclei, as they are still identifiable as active comets, in contrast to inactive comets that are observationally indistinguishable from low albedo asteroids. In this paper, we present a synthesis of comet and meteor observations of Jupiter-family comet 209P/LINEAR, one of the most weakly active comets recorded to-date. Images taken by the Xingming 0.35-m telescope and the Gemini Flamingo-2 camera are modeled by a Monte Carlo dust model, which yields a low dust ejection speed ($1/10$ of that of moderately active comets), dominance of large dust grains, and a low dust production of $0.4~\mathrm{kg \cdot s^{-1}}$ at 19~d after the 2014 perihelion passage. We also find a reddish nucleus of 209P/LINEAR that is similar to D-type asteroids and most Trojan asteroids. Meteor observations with the Canadian Meteor Orbit Radar (CMOR), coupled with meteoroid stream modeling, suggest a low dust production of the parent over the past few hundred orbits, although there are hints of a some temporary increase in activity in the 18th century. Dynamical simulations indicate 209P/LINEAR may have resided in a stable near-Earth orbit for $\sim 10^4$~yr, which is significantly longer than typical JFCs. All these lines of evidence imply that 209P/LINEAR as an aging comet quietly exhausting its remaining near surface volatiles. We also compare 209P/LINEAR to other low activity comets, where evidence for a diversity of the origin of low activity is seen.

The Quadrantid meteor shower is among the strongest annual meteor showers, and has drawn the attention of scientists for several decades. The stream is unusual, among others, for several reasons: its very short duration around maximum activity (~12 - 14 hours) as detected by visual, photographic and radar observations, its recent onset (around 1835 AD) and because it had been the only major stream without an obvious parent body until 2003. Ever since, there have been debates as to the age of the stream and the nature of its proposed parent body, asteroid 2003 EH1. In this work, we present results on the most probable age and formation mechanism of the narrow portion of the Quadrantid meteoroid stream. For the first time we use data on eight high precision photographic Quadrantids, equivalent to gram - kilogram size, to constrain the most likely age of the core of the stream. Out of eight high-precision photographic Quadrantids, five pertain directly to the narrow portion of the stream. In addition, we also use data on five high-precision radar Quadrantids, observed within the peak of the shower. We performed backwards numerical integrations of the equations of motion of a large number of 'clones' of both, the eight high-precision photographic and five radar Quadrantid meteors, along with the proposed parent body, 2003 EH1. According to our results, from the backward integrations, the most likely age of the narrow structure of the Quadrantids is between 200 - 300 years. These presumed ejection epochs, corresponding to 1700 - 1800 AD, are then used for forward integrations of large numbers of hypothetical meteoroids, ejected from the parent 2003 EH$_1$, until the present epoch. The aim is to constrain whether the core of the Quadrantid meteoroid stream is consistent with a previously proposed relatively young age (~ 200 years).}

<abridged> NASA's Asteroid Redirect Mission (ARM) has been proposed with the aim to capture a small asteroid a few meters in size and redirect it into an orbit around the Moon. There it can be investigated at leisure by astronauts aboard an Orion or other spacecraft. The target for the mission has not yet been selected, and there are very few potential targets currently known. Though sufficiently small near-Earth asteroids (NEAs) are thought to be numerous, they are also difficult to detect and characterize with current observational facilities. Here we collect the most up-to-date information on the smallest known near-Earth asteroids to outline the properties of these small NEAs as currently understood, in order to examine what the eventual ARM target might be like. Observational biases certainly mean that our sample is not an ideal representation of the true population of small NEAs. However our sample is representative of the eventual target list for the ARM mission, which will be compiled under very similar observational constraints unless dramatic changes are made to the way near-Earth asteroids are searched for and studied. We find that the typical rotation period is 40 minutes. The mean and median axis ratios were 1.43 and 1.29. Rotation rates much faster than the spin barrier are seen, reaching below 30 seconds, and implying that most of these bodies are monoliths. Non-principal axis rotation is uncommon. Axial ratios often reach values as high as two, though no undisputed results reach above three. We find little correlation of axis ratio with size. The most common spectral type in the sample of small NEAs is S-type (> 90%), with only a handful of C and X types known.

The Kreutz family of sungrazing comets contains over 2~000 known members, many of which are believed to be under $\sim 100$~m sizes (\textit{mini~comets}) and have only been studied at small heliocentric distances ($r_\mathrm{H}$) with space-based SOHO/STEREO spacecraft. To understand the brightening process of mini Kreutz comets, we conducted a survey using CFHT/MegaCam at moderate $r_\mathrm{H}$ guided by SOHO/STEREO observations. We identify two comets that should be in our search area but are not detected, indicating that the comets have either followed a steeper brightening rate within the previously-reported rapid brighten stage (the \textit{brightening~burst}), or the brightening burst starts earlier than expected. We present a composite analysis of pre-perihelion light-curves of five Kreutz comets that covered to $\sim 1$~AU. We observe a significant diversity in the light-curves that can be used to grossly classify them into two types: C/Ikeya-Seki and C/SWAN follow the canonical $r_\mathrm{H}^{-4}$ while the others follow $r_\mathrm{H}^{-7}$. In particular, C/SWAN seems to have undergone an outburst ($\Delta m>5$~mag) or a rapid brightening ($n\gtrsim11$) between $r_\mathrm{H}=1.06$~AU to $0.52$~AU, and shows hints of structural/compositional differences compared to other bright Kreutz comets. We also find evidence that the Kreutz comets as a population lose their mass less efficiently than the dynamically new comet, C/ISON, and are relatively devoid of species that drive C/ISON's activity at large $r_\mathrm{H}$. Concurrent observations of C/STEREO in different wavelengths also suggest that a blue-ward species such as CN may be the main driver for brightening burst instead of previously-thought sodium.

<Abridged> The impact of a meteoroid onto an asteroid transfers linear and angular momentum to the larger body, which may affect its orbit and its rotational state. Here we show that the meteoroid environment of our Solar System can have an effect on small asteroids that is comparable to the Yarkovsky and Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effects under certain conditions. The momentum content of the meteoroids themselves is expected to generate an effect much smaller than that of the Yarkovsky effect. However, momentum transport by ejecta may increase the net effective force by two orders of magnitude for impacts into bare rock surfaces. This result is sensitive to the extrapolation of laboratory microcratering experiment results to real meteoroid-asteroid collisions and needs further study. If this extrapolation holds, then meteoroid impacts are more important to the dynamics of small asteroids than had previously been considered.

We acquired Johnson BV photometry of the binary Be disk system delta~Scorpii during its 2009, 2010, 2011, and 2012 observing seasons and used it to probe the innermost regions of the disk. We found that several disk building events have occurred during this time, resulting in an overall brightening in the V-band and reddening of the system. In addition to these long-term trends, we found cyclical variability in each observing season on timescales between 60 and 100 days. We were able to reproduce the changes in the magnitude and colour of delta Sco using our theoretical models and found that variable mass-loss rates in the range 2.5-7.0x10^{-9} M_{sun}/yr over ~35 days can reproduce the observed increase in brightness.

The arrival of solid particles from outside our solar system would present us with an invaluable source of scientific information. Attempts to detect such interstellar particles among the meteors observed in Earth's atmosphere have almost exclusively assumed that those particles moving above the Solar System's escape speed -- particles on orbits hyperbolic with respect to the Sun-- were precisely the extrasolar particles being searched for. Here we show that hyperbolic particles can be generated entirely within the Solar System by gravitational scattering of interplanetary dust and meteoroids by the planets. These particles have necessarily short lifetimes as they quickly escape our star system; nonetheless some may arrive at Earth at speeds comparable to those expected of interstellar meteoroids. Some of these are associated with the encounter of planets with the debris streams of individual comets; however, such encounters are relatively rare. The rates of occurrence of hyperbolically-scattered sporadic meteors are also quite low. Only one of every 10,000 optical meteors observed at Earth is expected to be such a locally generated hyperbolic and its heliocentric velocity is typically only a hundred meters per second above the heliocentric escape velocity at Earth's orbit. Mercury and Venus are predicted to generate weak 'hyperbolic meteor showers': the restrictive geometry of scattering to our planet means that a radiant near the Sun from which hyperbolic meteors arrive at Earth should recur with the planet's synodic period. However, though planetary scattering can produce meteoroids with speeds comparable to interstellar meteors and at fluxes near current upper limits for such events, the majority of this locally-generated component of hyperbolic meteoroids is just above the heliocentric escape velocity and should be easily distinguishable from true interstellar meteoroids.

An unexpected intense outburst of the Draconid meteor shower was detected by the Canadian Meteor Orbit Radar (CMOR) on October 8, 2012. The peak flux occurred at ~16:40 UT on October 8 with a maximum of 2.4 +/- 0.3 hr-1 km-2 (appropriate to meteoroid mass larger than 10-7 kg), equivalent to a ZHRmax = 9000 +/- 1000 using 5-minute intervals, using a mass distribution index of s = 1.88 +/- 0.01 as determined from the amplitude distribution of underdense Draconid echoes. This makes the out- burst among the strongest Draconid returns since 1946 and the highest flux shower since the 1966 Leonid meteor storm, assuming a constant power-law distribution holds from radar to visual meteoroid sizes. The weighted mean geocentric radiant in the time interval of 15-19h UT, Oct 8, 2012 was {\alpha}g = 262.4 +/- 0.1 deg, {\delta}g = 55.7 +/- 0.1 deg (epoch J2000.0). Visual observers also reported increased activity around the peak time, but with a much lower rate (ZHR 200), suggesting that the magnitude-cumulative num- ber relationship is not a simple power-law. Ablation modeling of the observed meteors as a population does not yield a unique solution for the grain size and distribution of Draconid meteoroids, but is consistent with a typical Draconid meteoroid of mtotal between 10-6 to 10-4 kg being composed of 10 - 100 grains. Dynamical simulations indicate that the outburst was caused by dust particles released during the 1966 per- ihelion passage of the parent comet, 21P/Giacobini-Zinner, although there are dis- crepancies between the modelled and observed timing of the encounter, presumably caused by approaches of the comet to Jupiter during 1966-1972. Based on the results of our dynamical simulation, we predict possible increased activity of the Draconid meteor shower in 2018, 2019, 2021 and 2025.

Previous studies have suggested that Comet 209P/LINEAR may produce strong me- teor activity on Earth on 2014 May 24; however, exact timing and activity level is difficult to estimate due to the limited physical observations of the comet. Here we reanalyze the optical observations of 209P/LINEAR obtained during its 2009 appari- tion. We find that the comet is relatively depleted in dust production, with Af {\rho} at 1 cm level within eight months around its perihelion. This feature suggested that this comet may be currently transitioning from typical comet to a dormant comet. Syn- dyne simulation shows that the optical cometary tail is dominated by larger particles with {\beta} ~ 0.003. Numerical simulations of the cometary dust trails confirm the arrival of particles on 2014 May 24 from some of the 1798-1979 trails. The nominal radiant is at RA 122 +/- 1 deg, Dec 79 +/- 1 deg (J2000) in the constellation of Camelopardalis. Given that the comet is found to be depleted in dust production, we concluded that a meteor storm (ZHR>=1000) may be unlikely. However, our simulation also shows that the size distribution of the arrived particles is skewed strongly to larger particles. Coupling with the result of syndyne simulation, we think that the event, if detectable, may be dominated by bright meteors. We encourage observers to monitor the expected meteor event as it will provide us with rare direct information on the dynamical history of 209P/LINEAR which is otherwise irretrievably lost.

The Andromedid meteor shower underwent spectacular outbursts in 1872 and 1885, producing thousands of visual meteors per hour and described as `stars fell like rain' in Chinese records of the time. The shower originates from comet 3D/Biela whose disintegration in the mid-1800's is linked to the outbursts, but the shower has been weak or absent since the late 19th Century. This shower returned in December 2011 with a zenithal hourly rate of approximately 50, the strongest return in over a hundred years. Some 122 probable Andromedid orbits were detected by the Canadian Meteor Orbit Radar. The shower outburst occurred during 2011 Dec 3-5. The radiant at RA +$18\degree$ and Dec +$56\degree$ is typical of the `classical' Andromedids of the early 1800's, whose radiant was actually in Cassiopeia. The orbital elements indicate that the material involved was released before 3D/Biela's breakup prior to 1846. The observed shower in 2011 had a slow geocentric speed (16 km s$^{-1}$) and was comprised of small particles: the mean measured mass from the radar is $\sim5 \times 10^{-7}$ kg corresponding to radii of 0.5 mm at a bulk density of 1000 kg/m$^3$. Numerical simulations of the parent comet indicate that the meteoroids of the 2011 return of the Andromedids shower were primarily ejected during 3D/Biela's 1649 perihelion passage. The orbital characteristics, radiant, timing as well as the absence of large particles in the streamlet are all consistent with simulations. Predictions are made regarding other appearances of the shower in the years 2000-2047 based on our numerical model. We note that the details of the 2011 return can, in principle, be used to better constrain the orbit of 3D/Biela prior to the comets first recorded return in 1772.

The logarithmic potential is of great interest and relevance in the study of the dynamics of galaxies. Some small corrections to the work of Contopoulos & Seimenis (1990) who used the method of Prendergast (1982) to find periodic orbits and bifurcations within such a potential are presented. The solution of the orbital radial equation for the purely radial logarithmic potential is then considered using the p-ellipse (precessing ellipse) method pioneered by Struck (2006). This differential orbital equation is a special case of the generalized Burgers equation. The apsidal angle is also determined, both numerically as well as analytically by means of the Lambert W and the Polylogarithm functions. The use of these functions in computing the gravitational lensing produced by logarithmic potentials is discussed.

We report high resolution multi-station observations of meteors by the Canadian Automated Meteor Observatory (CAMO) recorded from June 2009 to August 2010. Our survey has a limiting detection magnitude of +5 mag in R-band, equivalent to a limiting meteoroid mass of ~2*E-7 kg. The high metric trajectory accuracy (of the order of 30 m perpendicular to the solution and 200 m along-track) allows us to determine velocities with average uncertainty of < 1.5% in speed and ~0.4 degr in radiant direction. A total of 1739 meteors had measured orbits. The data has been searched for meteors in hyperbolic orbits, which are potentially of interstellar origin. We found 22 potential hyperbolic meteors among our sample, with only two of them having a speed at least three sigma above the hyperbolic limit. For our one year survey we find no clear evidence of interstellar meteoroids at mm-sizes in a weighted time-area product of ~1*E4 km^2*h. Backward integrations performed for these 22 potentially hyperbolic meteors to check for close encounters with planets show no considerable changes in their orbits. Detailed examination leads us to conclude that our few identified events are most likely the result of measurement error. We find an upper limit of f_ISP < 2*E-4/(km^2*h) for the flux of interstellar meteoroids at Earth with a limiting mass of m > 2*E-7 kg.

We have carried out a search for Main Belt Asteroids (MBAs) co-orbiting with the large MBA Vesta and the dwarf planet Ceres. Through improving the search criteria used in (Christou, 2000b) and numerical integrations of candidate coorbitals, we have identified approximately 51 (44) objects currently in co-orbital libration with Ceres (Vesta). We show that these form part of a larger population of transient coorbitals; 129 (94) MBAs undergo episodes of co-orbital libration with Ceres (Vesta) within a 2 Myr interval centred on the present. The lifetime in the resonance is typically a few times 10^5 yr butcan exceed 2*10^6 yr. The variational properties of the orbits of several co-orbitals were examined. It was found that their present states with respect to the secondary are well determined but knowledge of it is lost typically after ~2*10^5 years. Objects initially deeper into the coorbital region maintain their coorbital state for longer. Using the model of Namouni et al. (1999) we show that their dynamics are similar to those of temporary coorbital NEAs of the Earth and Venus. As in that case, the lifetime of resonant libration is dictated by planetary secular perturbations, the inherent chaoticity of the orbits and close encounters with massive objects other than the secondary. In particular we present evidence that, while in the coorbital state, close encounters with the secondary are generally avoided and that Ceres affects the stability of tadpole librators of Vesta. Finally we demonstrate the existence of Quasi-satellite orbiters of both Ceres and Vesta and conclude that decametre-sized objects detected in the vicinity of Vesta by the DAWN mission may, in fact, belong to this dynamical class rather than be bona-fide (i.e. keplerian) satellites of Vesta.

Spectral data are presented for comets 2006 VZ13 (LINEAR), 2006 K4 (NEAT), 2006 OF2 (Broughton), 2P/Encke, and 93P/Lovas I, obtained with the Cerro-Tololo Inter-American Observatory 1.5-m telescope in August 2007. Comet 2006 VZ13 is a new Oort cloud comet and shows strong lines of CN (3880 angstroms), the Swan band sequence for C_2 (4740, 5160, and 5630 angstroms), C_3 (4056 angstroms), and other faint species. Lines are also identified in the spectra of the other comets. Flux measurements of the CN, C_2 (Delta v = +1,0), and C_3 lines are recorded for each comet and production rates and ratios are derived. When considering the comets as a group, there is a correlation of C_2 and C_3 production with CN, but there is no conclusive evidence that the production rate ratios depend on heliocentric distance. The continuum is also measured, and the dust production and dust-to-gas ratios are calculated. There is a general trend, for the group of comets, between the dust-to-gas ratio and heliocentric distance, but it does not depend on dynamical age or class. Comet 2006 VZ13 is determined to be in the carbon-depleted (or Tempel 1 type) class.

The Canada-France-Hawaii Telescope Legacy Survey, specifically the Very Wide segment of data, is used to search for possible main-belt comets. In the first data set, 952 separate objects with asteroidal orbits within the main-belt are examined using a three-level technique. First, the full-width-half-maximum of each object is compared to stars of similar magnitude, to look for evidence of a coma. Second, the brightness profiles of each object are compared with three stars of the same magnitude, which are nearby on the image to ensure any extended profile is not due to imaging variations. Finally, the star profiles are subtracted from the asteroid profile and the residuals are compared with the background using an unpaired T-test. No objects in this survey show evidence of cometary activity. The second survey includes 11438 objects in the main-belt, which are examined visually. One object, an unknown comet, is found to show cometary activity. Its motion is consistent with being a main-belt asteroid, but the observed arc is too short for a definitive orbit calculation. No other body in this survey shows evidence of cometary activity. Upper limits of the number of weakly and strongly active main-belt comets are derived to be 630+/-77 and 87+/-28, respectively. These limits are consistent with those expected from asteroid collisions. In addition, data extracted from the Canada-France-Hawaii Telescope image archive of main-belt comet 176P/LINEAR is presented.

Asteroid 3200 Phaethon resembles a comet in some ways, including a highly-eccentric orbit (e=0.89) and a strong associated meteor shower (the Geminids). Yet this object has never been observed to exhibit any cometary activity, i.e., gas production. We observed 3200 Phaethon with the Caltech Submillimeter Observatory on two occasions, once while it was near its closest approach to Earth as it neared perihelion, and another while it was further from Earth post-perihelion. Observations of the J=2-1 and J=3-2 rotational transitions of 12CO, typically strong lines in comets and indicative of gas production, yielded no detection. Upper limits on the 12CO production of 1.8e28 molecules/s and 7.6e28 molecules/s for Phaethon were determined on these two occasions.

We present the results of a project to detect small (~1 km) main-belt asteroids with the 3.6 meter Canada-France-Hawaii Telescope (CFHT). We observed in 2 filters (MegaPrime g' and r') in order to compare the results in each band. Owing to the observational cadence we did not observe the same asteroids through each filter and thus do not have true colour information. However strong differences in the size distributions as seen in the two filters point to a colour-dependence at these sizes, perhaps to be expected in this regime where asteroid cohesiveness begins to be dominated by physical strength and composition rather than by gravity. The best fit slopes of the cumulative size distributions (CSDs) in both filters tend towards lower values for smaller asteroids, consistent with the results of previous studies. In addition to this trend, the size distributions seen in the two filters are distinctly different, with steeper slopes in r' than in g'. Breaking our sample up according to semimajor axis, the difference between the filters in the inner belt is found to be somewhat less pronounced than in the middle and outer belt, but the CSD of those asteroids seen in the r' filter is consistently and significantly steeper than in g' throughout. The CSD slopes also show variations with semimajor axis within a given filter, particularly in r'. We conclude that the size distribution of main belt asteroids is likely to be colour dependent at kilometer sizes and that this dependence may vary across the belt.

In this work we present analytic and numerical treatments of the gravitational wave signal from a pulsar which includes spindown. We consider phase corrections to a received monochromatic signal due to rotational and elliptical orbital motion of the Earth, as well as perturbations due to Jupiter and the Moon. We discuss the Fourier transform of such a signal, which is expressed in terms of well known special functions and lends itself to a tractable numerical analysis.

Observational searches for asteroids orbiting near Earth's triangular Lagrange points face unique obstacles. A population of such asteroids would occupy a large projected area on the sky (possibly hundreds of square degrees) and is not favorably placed with respect to the Sun. Here we examine the properties of synthetic populations of Earth ``Trojans'' in order to aid in the optimization of observational searches for them. We find that the highest on-sky projected number densities are not located at the positions of the L4 and L5 points themselves, but rather a few degrees closer to the Sun. Also, asteroids on orbits about the L4 and L5 points typically brighten as the difference between their ecliptic longitude and that of the Sun increases owing to phase effects, but their number density on the sky concurrently falls rapidly.

A simple question of celestial mechanics is investigated: in what regions of phase space near a binary system can planets persist for long times? The planets are taken to be test particles moving in the field of an eccentric binary system. A range of values of the binary eccentricity and mass ratio is studied, and both the case of planets orbiting close to one of the stars, and that of planets outside the binary orbiting the system's center of mass, are examined. From the results, empirical expressions are developed for both 1) the largest orbit around each of the stars, and 2) the smallest orbit around the binary system as a whole, in which test particles survive the length of the integration (10^4 binary periods). The empirical expressions developed, which are roughly linear in both the mass ratio mu and the binary eccentricity e, are determined for the range 0.0 <= e <= 0.7-0.8 and 0.1 <= mu <= 0.9 in both regions, and can be used to guide searches for planets in binary systems. After considering the case of a single low-mass planet in binary systems, the stability of a mutually-interacting system of planets orbiting one star of a binary system is examined, though in less detail.

We study the evolution of long-period comets by numerical integration of their orbits, following comets from their origin in the Oort cloud until their final escape or destruction, in a model solar system consisting of the Sun, the four giant planets and the Galactic tide. We also examine the effects of non-gravitational forces and the gravitational forces from a hypothetical solar companion or circumsolar disk. We confirm the conclusion of Oort and other investigators that the observed distribution of long-period comet orbits does not match the expected steady-state distribution unless there is fading or some similar process that depletes the population of older comets. We investigate several simple fading laws. We can match the observed orbit distribution if the fraction of comets remaining observable after m apparitions is proportional to m to the power -0.6 +/- 0.1 (close to the fading law originally proposed by Whipple 1962); or if approximately 95% of comets live for only a few (~6) returns and the remainder last indefinitely. Our results also yield statistics such as the expected perihelion distribution, distribution of aphelion directions, frequency of encounters with the giant planets and the rate of production of Halley-type comets.

This paper investigates the long-term orbital stability of small bodies near the central binary of the Alpha Centauri system. Test particles on circular orbits are integrated in the field of this binary for 32000 binary periods or approximately 2.5 Myr. In the region exterior to the binary, particles with semi-major axes less than roughly three times the binary's semi-major axis are unstable. Inside the binary, particles are unstable if further than 0.2 binary semimajor axes from the primary, with stablility closer in a strong function of inclination: orbits inclined near 90 degrees are unstable in as close as 0.01 binary semimajor axes from either star.