Burns et al. (1979) use the parameter beta to describe the ratio of radiation pressure to gravity experienced by a small particle in the solar system. The central potential that these particles experience is effectively reduced by a factor of (1 - beta), lowering the escape velocity. Burns et al. (1979) also provide a simple expression for the value of beta at which particles ejected from a comet follow parabolic orbits and thus leave the solar system; this paper expands on this to calculate critical beta values that take ejection velocity into account. We first derive an expression for the speed at which particles are ejected from a comet that is valid at all particle sizes. We then calculate beta values for both idealized and real materials to determine which particles will be lost from meteor streams. We perform these calculations for cometary ejecta leading, trailing, and following the parent comet's nucleus for 10 major meteor showers. These values bound mass regimes within which meteoroids are ejected from the solar system and therefore cannot contribute to meteor showers.

The number density and flux of a meteoroid stream is enhanced near a massive body due to the phenomenon known as gravitational focusing. The greatest enhancement occurs directly opposite the massive body from the stream radiant: as an observer approaches this location, the degree of focusing is unbound for a perfectly collimated stream. However, real meteoroid streams exhibit some dispersion in radiant and speed that will act to eliminate this singularity. In this paper, we derive an analytic approximation for this smoothing that can be used in meteoroid environment models and is based on real measurements of meteor shower radiant dispersion.

Orbital dissimilarity, or D, criteria are often used to select members of a meteor shower from a set of meteor observations. These criteria provide a quantitative description of the degree to which two orbits differ; if the degree of dissimilarity between a shower's reference orbit and an individual meteor does not exceed a selected threshold, the meteor is considered to be a member of that shower. However, members of a meteor shower tend to disperse in longitude of the ascending node (and thus in solar longitude) while preserving a common Sun-centered ecliptic radiant. Employing dissimilarity criteria to judge shower membership may therefore make the shower appear briefer than it actually is. We demonstrate this effect for two simulated meteor showers and assess the maximum permitted deviation in solar longitude as a function of radiant and velocity measurement error.

Meteoroids pose one of the largest risks to spacecraft outside of low Earth orbit. In order to correctly predict the rate at which meteoroids impact and damage spacecraft, environment models must describe the mass, directionality, velocity, and density distributions of meteoroids. NASA's Meteoroid Engineering Model (MEM) is one such model; MEM 3 is an updated version of the code that better captures the correlation between directionality and velocity and incorporates a bulk density distribution. This paper describes MEM 3 and compares its predictions with the rate of large particle impacts seen on the Long Duration Exposure Facility (LDEF) and the Pegasus II and III satellites.

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.

NASA's Meteoroid Environment Office (MEO) produces an annual meteor shower forecast in order to help spacecraft operators assess the risk posed by meteoroid streams. Previously, this forecast focused on the International Space Station and therefore reported meteoroid fluxes and enhancement factors at an orbital altitude of 400 km. This paper presents an updated forecast algorithm that has an improved calculation of the flux enhancement produced by showers and can calculate fluxes at any selected Earth or lunar orbital altitude. Finally, we discuss and generate forecasted fluxes for the 2018 Draconid meteor shower, which is expected to produce meteoroid flux enhancements near the Sun-Earth L1 and L2 Lagrange points but not at Earth.

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.

We present a new method to detect meteor showers using the Density-Based Spatial Clustering of Applications with Noise algorithm (DBSCAN; Ester et al. 1996). DBSCAN is a modern cluster detection algorithm that is well suited to the problem of extracting meteor showers from all-sky camera data because of its ability to efficiently extract clusters of different shapes and sizes from large datasets. We apply this shower detection algorithm on a dataset that contains 25,885 meteor trajectories and orbits obtained from the NASA All-Sky Fireball Network and the Southern Ontario Meteor Network (SOMN). Using a distance metric based on solar longitude, geocentric velocity, and Sun-centered ecliptic radiant, we find 25 strong cluster detections and 6 weak detections in the data, all of which are good matches to known showers. We include measurement errors in our analysis to quantify the reliability of cluster occurrence and the probability that each meteor belongs to a given cluster. We validate our method through false positive/negative analysis and with a comparison to an established shower detection algorithm.

Separating meteor showers from the sporadic meteor background is critical for the study of both showers and the sporadic complex. The linkage of meteors to meteor showers, to parent bodies, and to other meteors is done using measures of orbital similarity. These measures often take the form of so-called D-parameters and are generally paired with some cutoff value within which two orbits are considered related. The appropriate cutoff value can depend on the size of the data-set (Southworth & Hawkins 1963), the sporadic contribution within the observed size range (Jopek 1995), or the inclination of the shower (Galligan 2001). If the goal is to minimize sporadic contamination of the extracted shower, the cutoff value should also reflect the strength of the shower compared to the local sporadic background. In this paper, we present a method for determining, on a per-shower basis, the orbital similarity cutoff value that corresponds to a chosen acceptable false-positive rate. This method also assists us in distinguishing which showers are significant within a set of data. We apply these methods to optical meteor observations from the NASA All-Sky and Southern Ontario Meteor Networks.

Eighty planetary systems of two or more planets are known to orbit stars other than the Sun. For most, the data can be sufficiently explained by non-interacting Keplerian orbits, so the dynamical interactions of these systems have not been observed. Here we present 4 sets of lightcurves from the Kepler spacecraft, which each show multiple planets transiting the same star. Departure of the timing of these transits from strict periodicity indicates the planets are perturbing each other: the observed timing variations match the forcing frequency of the other planet. This confirms that these objects are in the same system. Next we limit their masses to the planetary regime by requiring the system remain stable for astronomical timescales. Finally, we report dynamical fits to the transit times, yielding possible values for the planets' masses and eccentricities. As the timespan of timing data increases, dynamical fits may allow detailed constraints on the systems' architectures, even in cases for which high-precision Doppler follow-up is impractical.

We present a method to confirm the planetary nature of objects in systems with multiple transiting exoplanet candidates. This method involves a Fourier-Domain analysis of the deviations in the transit times from a constant period that result from dynamical interactions within the system. The combination of observed anti-correlations in the transit times and mass constraints from dynamical stability allow us to claim the discovery of four planetary systems Kepler-25, Kepler-26, Kepler-27, and Kepler-28, containing eight planets and one additional planet candidate.

We present a new method for confirming transiting planets based on the combination of transit timingn variations (TTVs) and dynamical stability. Correlated TTVs provide evidence that the pair of bodies are in the same physical system. Orbital stability provides upper limits for the masses of the transiting companions that are in the planetary regime. This paper describes a non-parametric technique for quantifying the statistical significance of TTVs based on the correlation of two TTV data sets. We apply this method to an analysis of the transit timing variations of two stars with multiple transiting planet candidates identified by Kepler. We confirm four transiting planets in two multiple planet systems based on their TTVs and the constraints imposed by dynamical stability. An additional three candidates in these same systems are not confirmed as planets, but are likely to be validated as real planets once further observations and analyses are possible. If all were confirmed, these systems would be near 4:6:9 and 2:4:6:9 period commensurabilities. Our results demonstrate that TTVs provide a powerful tool for confirming transiting planets, including low-mass planets and planets around faint stars for which Doppler follow-up is not practical with existing facilities. Continued Kepler observations will dramatically improve the constraints on the planet masses and orbits and provide sensitivity for detecting additional non-transiting planets. If Kepler observations were extended to eight years, then a similar analysis could likely confirm systems with multiple closely spaced, small transiting planets in or near the habitable zone of solar-type stars.

We present a Bayesian surrogate model for the analysis of periodic or quasi-periodic time series data. We describe a computationally efficient implementation that enables Bayesian model comparison. We apply this model to simulated and real exoplanet observations. We discuss the results and demonstrate some of the challenges for applying our surrogate model to realistic exoplanet data sets. In particular, we find that analyses of real world data should pay careful attention to the effects of uneven spacing of observations and the choice of prior for the "jitter" parameter.

We report the distribution of planets as a function of planet radius (R_p), orbital period (P), and stellar effective temperature (Teff) for P < 50 day orbits around GK stars. These results are based on the 1,235 planets (formally "planet candidates") from the Kepler mission that include a nearly complete set of detected planets as small as 2 Earth radii (Re). For each of the 156,000 target stars we assess the detectability of planets as a function of R_p and P. We also correct for the geometric probability of transit, R*/a. We consider first stars within the "solar subset" having Teff = 4100-6100 K, logg = 4.0-4.9, and Kepler magnitude Kp < 15 mag. We include only those stars having noise low enough to permit detection of planets down to 2 Re. We count planets in small domains of R_p and P and divide by the included target stars to calculate planet occurrence in each domain. Occurrence of planets varies by more than three orders of magnitude and increases substantially down to the smallest radius (2 Re) and out to the longest orbital period (50 days, ~0.25 AU) in our study. For P < 50 days, the radius distribution is given by a power law, df/dlogR= k R^\alpha. This rapid increase in planet occurrence with decreasing planet size agrees with core-accretion, but disagrees with population synthesis models. We fit occurrence as a function of P to a power law model with an exponential cutoff below a critical period P_0. For smaller planets, P_0 has larger values, suggesting that the "parking distance" for migrating planets moves outward with decreasing planet size. We also measured planet occurrence over Teff = 3600-7100 K, spanning M0 to F2 dwarfs. The occurrence of 2-4 Re planets in the Kepler field increases with decreasing Teff, making these small planets seven times more abundant around cool stars than the hottest stars in our sample. [abridged]

Doppler planet searches have discovered that giant planets follow orbits with a wide range of orbital eccentricities, revolutionizing theories of planet formation. The discovery of hundreds of exoplanet candidates by NASA's Kepler mission enables astronomers to characterize the eccentricity distribution of small exoplanets. Measuring the eccentricity of individual planets is only practical in favorable cases that are amenable to complementary techniques (e.g., radial velocities, transit timing variations, occultation photometry). Yet even in the absence of individual eccentricities, it is possible to study the distribution of eccentricities based on the distribution of transit durations (relative to the maximum transit duration for a circular orbit). We analyze the transit duration distribution of Kepler planet candidates. We find that for host stars with T_eff > 5100 K we cannot invert this to infer the eccentricity distribution at this time due to uncertainties and possible systematics in the host star densities. With this limitation in mind, we compare the observed transit duration distribution with models to rule out extreme distributions. If we assume a Rayleigh eccentricity distribution for Kepler planet candidates, then we find best-fits with a mean eccentricity of 0.1-0.25 for host stars with T_eff < 5100 K. We compare the transit duration distribution for different subsets of Kepler planet candidates and discuss tentative trends with planetary radius and multiplicity. High-precision spectroscopic follow-up observations for a large sample of host stars will be required to confirm which trends are real and which are the results of systematic errors in stellar radii. Finally, we identify planet candidates that must be eccentric or have a significantly underestimated stellar radius.

We present and discuss five candidate exoplanetary systems identified with the Kepler spacecraft. These five systems show transits from multiple exoplanet candidates. Should these objects prove to be planetary in nature, then these five systems open new opportunities for the field of exoplanets and provide new insights into the formation and dynamical evolution of planetary systems. We discuss the methods used to identify multiple transiting objects from the Kepler photometry as well as the false-positive rejection methods that have been applied to these data. One system shows transits from three distinct objects while the remaining four systems show transits from two objects. Three systems have planet candidates that are near mean motion commensurabilities---two near 2:1 and one just outside 5:2. We discuss the implications that multitransiting systems have on the distribution of orbital inclinations in planetary systems, and hence their dynamical histories; as well as their likely masses and chemical compositions. A Monte Carlo study indicates that, with additional data, most of these systems should exhibit detectable transit timing variations (TTV) due to gravitational interactions---though none are apparent in these data. We also discuss new challenges that arise in TTV analyses due to the presence of more than two planets in a system.

The observed extrasolar planets possess both large masses (with a median M sin i of 1.65 MJ) and a wide range in orbital eccentricity (0 < e < 0.94). As planets are thought to form in circumstellar disks, one important question in planet formation is determining whether, and to what degree, a gaseous disk affects an eccentric planet's orbit. Recent studies have probed the interaction between a disk and a terrestrial planet on an eccentric orbit, and the interaction between a disk and a gas giant on a nearly circular orbit, but little is known about the interaction between a disk and an eccentric gas giant. Such a scenario could arise due to scattering while the disk is still present, or perhaps through planet formation via gravitational instability. We fill this gap with simulations of eccentric, massive (gap-forming) planets in disks using the hydrodynamical code FARGO. Although the long-term orbital evolution of the planet depends on disk properties, including the boundary conditions used, the disk always acts to damp eccentricity when the planet is released into the disk. This eccentricity damping takes place on a timescale of 40 years, 15 times faster than the migration timescale.

The extrasolar planets discovered to date possess unexpected orbital elements. Most orbit their host stars with larger eccentricities and smaller semi-major axes than similarly sized planets in our own solar system do. It is generally agreed that the interaction between giant planets and circumstellar disks (Type II migration) drives these planets inward to small radii, but the effect of these same disks on orbital eccentricity, e, is controversial. Several recent analytic calculations suggest that disk-planet interactions can excite eccentricity, while numerical studies generally produce eccentricity damping. This paper addresses this controversy using a quasi-analytic approach, drawing on several preceding analytic studies. This work refines the current treatment of eccentricity evolution by removing several approximations from the calculation of disk torques. We encounter neither uniform damping nor uniform excitation of orbital eccentricity, but rather a function de/dt that varies in both sign and magnitude depending on eccentricity and other solar system properties. Most significantly, we find that for every combination of disk and planet properties investigated herein, corotation torques produce negative values of de/dt for some range in e within the interval [0.1, 0.5]. If corotation torques are saturated, this region of eccentricity damping disappears, and excitation occurs on a short timescale of less than 0.08 Myr. Thus, our study does not produce eccentricity excitation on a timescale of a few Myr -- we obtain either eccentricity excitation on a short time scale, or eccentricity damping on a longer time scale. Finally, we discuss the implications of this result for producing the observed range in extrasolar planet eccentricity.

This paper presents a parametric study of giant planet migration through the combined action of disk torques and planet-planet scattering. The torques exerted on planets during Type II migration in circumstellar disks readily decrease the semi-major axes, whereas scattering between planets increases the orbital eccentricities. This paper presents a parametric exploration of the possible parameter space for this migration scenario using two (initial) planetary mass distributions and a range of values for the time scale of eccentricity damping (due to the disk). For each class of systems, many realizations of the simulations are performed in order to determine the distributions of the resulting orbital elements of the surviving planets; this paper presents the results of 8500 numerical experiments. Our goal is to study the physics of this particular migration mechanism and to test it against observations of extrasolar planets. The action of disk torques and planet-planet scattering results in a distribution of final orbital elements that fills the a-e plane, in rough agreement with the orbital elements of observed extrasolar planets. In addition to specifying the orbital elements, we characterize this migration mechanism by finding the percentages of ejected and accreted planets, the number of collisions, the dependence of outcomes on planetary masses, the time spent in 2:1 and 3:1 resonances, and the effects of the planetary IMF. We also determine the distribution of inclination angles of surviving planets and the distribution of ejection speeds for exiled planets.