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.

Context. Several authors predicted an outburst of the Draconid meteor shower in 2018, but with an uncertain level of activity. Aims. Optical meteor observations were used to derive the population and mass indices, flux, and radiant positions of Draconid meteors. Methods. 90 minutes of multi-station observations after the predicted peak of activity were performed using highly sensitive Electron Multiplying Charge Coupled Device (EMCCD) cameras. The data calibration is discussed in detail. A novel maximum likelihood estimation method of computing the population and mass index with robust error estimation was developed. We apply the method to observed Draconids and use the values to derive the flux. Meteor trajectories are computed and compared to predicted radiant positions from meteoroid ejection models. Results. We found that the mass index was $1.74 \pm 0.18$ in the 30 minute bin after the predicted peak, and $2.32 \pm 0.27$ in the next 60 minutes. The location and the dispersion of the radiant matches well to modeled values, but there is an offset of $0.4^{\circ}$ in solar longitude.

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.

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 physical composition and structure of meteoroids gives us insight into the formation processes of their parent asteroids and comets. The strength of and fundamental grain sizes in meteoroids tell us about the environment in which small solar system bodies formed, and the processes which built up these basic planetary building blocks. The structure of meteorites can be studied directly, but the set of objects which survive entry through the atmosphere is biased toward large, strong objects with slow encounter speeds. Fragile objects, small objects and objects with high relative speeds are very unlikely to survive impact with the atmosphere. Objects between 100 microns and 1 meter, which are not strong enough to survive the ablation process, must be studied by radar or optical methods. Large meteoroids, which produce bright fireballs, are generally studied by investigating their compressive strength when they fragment, and their strength can also be inferred indirectly from their end heights. Fragmentation in faint meteors can be inferred from interference in radar observations, or observed directly with high-resolution optical systems. Meteor light curves, begin heights and time-evolving spectra can also be used to infer meteoroid structure. This paper presents a review of the meteor observation methods currently used to infer the structure and fragmentation of meteoroids in the millimeter to meter size range, and the current state of understanding these observations have given us.

Sixteen years of meteor radar data from the Canadian Meteor Orbit Radar (CMOR) were used to investigate the link between observed meteor rates and both solar and geomagnetic activity. Meteor rates were corrected for transmitter power and receiver noise, and seasonal effects were removed. A strong negative correlation is seen between solar activity, as measured with the 10.7 cm flux, and observed meteor rates. This lends support to the idea that heating in the atmosphere at times of elevated solar activity changes the scale height and therefore the length and maximum brightness of meteors; a larger scale height near solar maximum leads to longer, fainter meteors and therefore lower rates. A weaker negative correlation was observed with geomagnetic activity as measured with the $K$ index; this correlation was still present when solar activity effects were removed. Meteor activity at solar maximum is as much as 30\% lower than at solar minimum, strictly due to observing biases; geomagnetic activity usually affects meteor rates by less than 10 percent.

Twenty-one meteors showing double peaked light curves were analysed with observations collected with the Canadian Automated Meteor Observatory tracking system. Each event has orbital information, photometry, and at least one high-resolution observation. Two distinct light curve shapes were found: sudden double peaked curves, and smooth double peaked curves. The sudden peaked curves were produced by objects on asteroidal orbits and mostly showed noticeable fragmentation, while the smooth peaked curves were produced by cometary meteoroids and predominantly showed little to no visible fragmentation. An attempt to model these meteors as single bodies with two chemical components was unsuccessful, implying that fragmentation must be included in meteoroid ablation models.

Many existing optical meteor trajectory estimation methods use the approximation that the velocity of the meteor at the beginning of its luminous phase is equivalent to its velocity before atmospheric entry. Meteoroid kinetic energy loss prior to the luminous phase cannot be measured, but for some masses and entry geometries neglecting this loss may lead to non-negligible deceleration prior to thermal ablation. Using a numerical meteoroid ablation model, we simulate the kinematics of meteoroids beginning at 180 km with initial velocities ranging from 11 km/s to 71 km/s, and compare model velocities at the moment of detection to measurements. We validate the simulations by comparing the simulated luminous beginning heights with observed beginning heights of different populations of meteors detected with different optical systems. We find that most low-velocity meteoroids have a significant velocity difference of 100 m/s to 750 m/s (depending on meteoroid type, mass, and observation system). This systematic underestimate of meteoroid speeds also results in systematically lower semi-major axes for meteoroid orbits.

Luminous efficiency is a necessary parameter for determining meteoroid mass from optical emission. Despite this importance, it is very poorly known, with previous results varying by up to two orders of magnitude for a given speed. We present the most recent study of luminous efficiency values determined with modern high-resolution instruments, by directly comparing dynamic and photometric meteoroid masses. Fifteen non-fragmenting meteoroids were used, with a further five clearly fragmenting events for comparison. Twelve of the fifteen non-fragmenting meteoroids had luminous efficiencies less than 1%, while the fragmenting meteoroids had upper limits of a few tens of percent. No clear trend with speed was seen, but there was a weak negative trend of luminous efficiency on meteoroid mass, implying that smaller meteoroids radiate more efficiently.

Context. Generating a synthetic dataset of meteoroid orbits is a crucial step in analysing the probabilities of random grouping of meteoroid orbits in automated meteor shower surveys. Recent works have shown the importance of choosing a low similarity threshold value of meteoroid orbits, some pointing out that the recent meteor shower surveys produced false positives due to similarity thresholds which were too high. On the other hand, the methods of synthetic meteoroid orbit generation introduce additional biases into the data, thus making the final decision on an appropriate threshold value uncertain. Aims. As a part of the ongoing effort to determine the nature of meteor showers and improve automated methods, it was decided to tackle the problem of synthetic meteoroid orbit generation, the main goal being to reproduce the underlying structure and the statistics of the observed data in the synthetic orbits. Methods. A new method of generating synthetic meteoroid orbits using the Kernel Density Estimation method is presented. Several types of approaches are recommended, depending on whether one strives to preserve the data structure, the data statistics or to have a compromise between the two. Results. The improvements over the existing methods of synthetic orbit generation are demonstrated. The comparison between the previous and newly developed methods are given, as well as the visualization tools one can use to estimate the influence of different input parameters on the final data.

The luminous efficiency of meteors is poorly known, but critical for determining the meteoroid mass. We present an uncertainty analysis of the luminous efficiency as determined by the classical ablation equations, and suggest a possible method for determining the luminous efficiency of real meteor events. We find that a two-term exponential fit to simulated lag data is able to reproduce simulated luminous efficiencies reasonably well.

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 Canadian Automated Meteor Observatory (CAMO) detects occasional meteors with two maxima in the image intensified CCD based light curves. We report early results from an analysis of 21 of these events. Most of these events show qualitatively similar light curves, with a rounded luminous peak, followed by an almost linear sharp rise in the second peak, and a relatively rapid curved decay of the second peak. While a number of mechanisms could explain two maxima in the light curves, numerical modelling shows that most of these events can be matched by a simple dustball model in which some grains have been released well before intensive ablation begins, followed by a later release of core grains at a single time. Best fits to observations are obtained with the core grains being larger than the pre-released outer grains, with the core grains typically $10^{-6}$ kg while the early release grains are of the order of $10^{-9}$ kg.

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.

A strong outburst of the October Draconid meteor shower was predicted for October 8, 2011. Here we present the observations obtained by the Canadian Meteor Orbit Radar (CMOR) during the 2011 outburst. CMOR recorded 61 multi-station Draconid echoes and 179 single-station overdense Draconid echoes (covering the magnitude range of +3 <= MV <= +7) between 16-20h UT on October 8, 2011. The mean radiant for the outburst was determined to be a_g = 261.9o +/- 0.3o, d_g = +55.3o +/- 0.3o (J2000) from observations of the underdense multi-station echoes. This radiant location agrees with model predictions to ~1o . The determined geocentric velocity was found to be ~10-15% lower than the model value (17.0 - 19.1 km s^-1 versus 20.4 km s^-1), a discrepancy we attribute to undercorrection for atmospheric deceleration of low density Draconid meteoroids as well as to poor radar radiant geometry during the outburst peak. The mass index at the time of the outburst was determined to be ~ 1.75 using the amplitude distribution of underdense echoes...

The distribution of dust in the ecliptic plane between 0.96 and 1.04 AU has been inferred from impacts on the two STEREO spacecraft through observation of secondary particle trails and unexpected off-points in the Heliospheric Imager (HI) cameras. This study made use of analysis carried out by members of a distributed web-based project, Solar Stormwatch. A comparison between observations of the brightest particle trails and a survey of fainter trails shows consistent distributions. While there is no obvious correlation between this distribution and the occurrence of individual meteor streams at Earth, there are some broad longitudinal features in these distributions that are also observed in sources of the sporadic meteor population. The asymmetry in the number of trails seen by each spacecraft and the fact that there are many more unexpected off-points in the HI-B than in HI-A, indicates that the majority of impacts are coming from the apex direction. For impacts causing off-points in the HI-B camera these dust particles are estimated to have masses in excess of 10-17 kg with radii exceeding 0.1 {\mu}m. For off-points observed in the HI-A images, which can only have been caused by particles travelling from the anti-apex direction, the distribution is consistent with that of secondary 'storm' trails observed by HI-B, providing evidence that these trails also result from impacts with primary particles from an anti-apex source. It is apparent that the differential mass index of particles from the apex direction is consistently above 2. This indicates that the majority of the mass is within the smaller particles of this population. In contrast, the differential mass index of particles from the anti-apex direction (causing off-points in HI-A) is consistently below 2, indicating that the majority of the mass is to be found in larger particles of this distribution.

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.