We describe an improved technique for using the backscattered phase from meteor radar echo measurements just prior to the specular point ($t_{0}$) to calculate meteor speeds and their uncertainty. Our method, which builds on earlier work of Cervera et al (1997), scans possible speeds in the Fresnel distance - time domain with a dynamic, sliding window and derives a best-speed estimate from the resultant speed distribution. We test the performance of our method, called pre-$t_{0}$ speeds by sliding-slopes technique (PSSST), on transverse scattered meteor echoes observed by the Middle Atmosphere Alomar Radar System (MAARSY) and the Canadian Meteor Orbit Radar (CMOR), and compare the results to time-of-flight and Fresnel transform speed estimates. Our novel technique is shown to produce good results when compared to both model and speed measurements using other techniques. We show that our speed precision is $\pm$5$\%$ at speeds less than 40 km/s and we find that more than 90$\%$ of all CMOR multi-station echoes have PSSST solutions. For CMOR data, PSSST is robust against the selection of critical phase value and poor phase unwrapping. Pick errors of up to $\pm$6 pulses for meteor speeds less than about 50 km/s produce errors of less than $\pm$5$\%$ of the meteoroid speed. In addition, the width of the PSSST speed Kernel density estimate (KDE) is used as a natural measure of uncertainty that captures both noise and $t_0$ pick uncertainties.

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 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.

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.

The Hamburg (H4) meteorite fell on January 17, 2018 at 01:08 UT approximately 10km North of Ann Arbor, Michigan. More than two dozen fragments totaling under one kilogram were recovered, primarily from frozen lake surfaces. The fireball initial velocity was 15.83 $\pm$ 0.05 km/s, based on four independent records showing the fireball above 50 km altitude. The radiant had a zenith angle of 66.1 $\pm$ 0.29 degs and an azimuth of 121.56 $\pm$ 1.2 degs. The resulting low inclination ($<$ 1 deg) Apollo-type orbit has a large aphelion distance and Tisserand value relative to Jupiter (Tj) of $\sim$ 3. Two major flares dominant the energy deposition profile, centred at 24.1 and 21.7 km altitude respectively under dynamic pressures of 5-7 MPa. Our preferred total energy for the Hamburg fireball is 2 to 7 T TNT (8.4 to 28 x 10$^9$ J), which corresponds to a likely initial mass in the range of 60 to 225 kg or diameter between 0.3 to 0.5 m. Based on the model of Granvik et al (2018), the meteorite originated in an escape route from the mid-outer asteroid belt. Hamburg is the 14th known H-chondrite with an instrumentally-derived pre-atmospheric orbit, half of which have small ($<$5 deg) inclinations making connection with (6) Hebe problematic. A definitive parent body consistent with all 14 known H-chondrite orbits remains elusive

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.

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.

Context: Results from previous searches for new meteor showers in the combined Croatian Meteor Network and SonotaCo meteor databases suggested possible parent bodies for several newly identified showers. Aims: We aim to perform an analysis to validate the connection between the identified showers and candidate parent bodies. Methods: Simulated particles were ejected from candidate parent bodies, a dynamical modeling was performed and the results were compared to the real meteor shower observations. Results: From the 13 analysed cases, three were found to be connected with comets, four with asteroids which are possibly dormant comets, four were inconclusive or negative, and two need more observational data before any conclusions can be drawn.