The paper presents an analysis of Polish Fireball Network (PFN) observations of enhanced activity of the Southern Taurid meteor shower in 2005 and 2015. In 2005, between October 20 and November 10, seven stations of PFN determined 107 accurate orbits with 37 of them belonging to the Southern Taurid shower. In the same period of 2015, 25 stations of PFN recorded 719 accurate orbits with 215 orbits of the Southern Taurids. Both maxima were rich in fireballs which accounted to 17% of all observed Taurids. The whole sample of Taurid fireballs is quite uniform in the sense of starting and terminal heights of the trajectory. On the other hand a clear decreasing trend in geocentric velocity with increasing solar longitude was observed. Orbital parameters of observed Southern Taurids were compared to orbital elements of Near Earth Objects (NEO) from the NEODYS-2 database. Using the Drummond criterion $D'$ with threshold as low as 0.06, we found over 100 fireballs strikingly similar to the orbit of asteroid 2015 TX24. Several dozens of Southern Taurids have orbits similar to three other asteroids, namely: 2005 TF50, 2005 UR and 2010 TU149. All mentioned NEOs have orbital periods very close to the 7:2 resonance with Jupiter's orbit. It confirms a theory of a "resonant meteoroid swarm" within the Taurid complex that predicts that in specific years, the Earth is hit by a greater number of meteoroids capable of producing fireballs.

On 2010 October 13, the Apollo type 20-meter asteroid 2010 TB54 passed within 6.1 lunar distances from the Earth. On the same date, but 11.4 hours earlier, exactly at 02:52:32 UT, the sky over central Poland was illuminated by -8.6 magnitude PF131010 Ciechanow fireball. The trajectory and orbit of the fireball was computed using multi-station data of Polish Fireball Network (PFN). The results indicate that the orbit of the meteoroid which caused the PF131010 fireball is similar to the orbit of 2010 TB54 asteroid and both bodies may be related. Moreover, two days before appearance of Ciechanow fireball another small asteroid denoted as 2010 SX11 passed close to the Earth-Moon system. Its orbit is even more similar to the orbit of Ciechanow fireball parent body than in case of 2010 TB54. The PF131010 Ciechanow entered Earth's atmosphere with the velocity of 12.9 +\- 0.2 km/s and started to shine at height of 82.5 +\- 0.3 km. Clear deceleration started after first three seconds of flight, and the terminal velocity of the meteor was only 5.8 +\- 0.2 km/s at height of 29.3 +\- 0.1 km. Such a low value of terminal velocity indicates that fragments with total mass of around 2 kg could survive the atmospheric passage and cause fall of the meteorites. The predicted area of possible meteorite impact is computed and it is located near Grabowo village south of Ostroleka city.

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

On the night of Oct 18/19, 2012 at 00:23 UT a -14.7 mag Orionid fireball occurred over northeastern Poland. The precise orbit and atmospheric trajectory of the event is presented, based on the data collected by five video and one photographic Polish Fireball Network (PFN) stations. The beginning height of the meteor is 168.4 +\- 0.6 km which makes the PF191012 Myszyniec fireball the highest ever observed, well documented meteor not belonging to the Leonid shower. The ablation became the dominant source of light of the meteor at a height of around 115 km. The thermalization of sputtered particles is suggested to be the source of radiation above that value. The transition height of 115 km is 10-15 km below the transition heights derived for Leonids and it might suggest that the material of Leonids should be more fragile and have probably smaller bulk density than in case of Orionids.

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

We analyse the coarse-grained phase-space structure of the six Galaxy-scale dark matter haloes of the Aquarius Project using a state-of-the-art 6D substructure finder. Within r_50, we find that about 35% of the mass is in identifiable substructures, predominantly tidal streams, but including about 14% in self-bound subhaloes. The slope of the differential substructure mass function is close to -2, which should be compared to around -1.9 for the population of self-bound subhaloes. Near r_50 about 60% of the mass is in substructures, with about 30% in self-bound subhaloes. The inner 35 kpc of the highest resolution simulation has only 0.5% of its mass in self-bound subhaloes, but 3.3% in detected substructure, again primarily tidal streams. The densest tidal streams near the solar position have a 3-D mass density about 1% of the local mean, and populate the high velocity tail of the velocity distribution.

A new multi-dimensional Hierarchical Structure Finder (HSF) to study the phase-space structure of dark matter in N-body cosmological simulations is presented. The algorithm depends mainly on two parameters, which control the level of connectivity of the detected structures and their significance compared to Poisson noise. By working in 6D phase-space, where contrasts are much more pronounced than in 3D position space, our HSF algorithm is capable of detecting subhaloes including their tidal tails, and can recognise other phase-space structures such as pure streams and candidate caustics. If an additional unbinding criterion is added, the algorithm can be used as a self-consistent halo and subhalo finder. As a test, we apply it to a large halo of the Millennium Simulation, where 19 % of the halo mass are found to belong to bound substructures, which is more than what is detected with conventional 3D substructure finders, and an additional 23-36 % of the total mass belongs to unbound HSF structures. The distribution of identified phase-space density peaks is clearly bimodal: high peaks are dominated by the bound structures and low peaks belong mostly to tidal streams. In order to better understand what HSF provides, we examine the time evolution of structures, based on the merger tree history. Bound structures typically make only up to 6 orbits inside the main halo. Still, HSF can identify at the present time at least 80 % of the original content of structures with a redshift of infall as high as z <= 0.3, which illustrates the significant power of this tool to perform dynamical analyses in phase-space.

This paper reviews and analyses methods used to identify neighbours in 6D space and estimate the corresponding phase-space density. It compares SPH methods to 6D Delaunay tessellation on statical and dynamical realisation of single halo profiles, paying attention to the unknown scaling, S_G, used to relate the spatial dimensions to the velocity dimensions. The methods with local adaptive metric provide the best phase-space estimators. They make use of a Shannon entropy criterion combined with a binary tree partitioning and with SPH interpolation using 10-40 neighbours. Local scaling implemented by such methods, which enforces local isotropy of the distribution function, can vary by about one order of magnitude in different regions within the system. It presents a bimodal distribution, in which one component is dominated by the main part of the halo and the other one is dominated by the substructures. While potentially better than SPH techniques, since it yields an optimal estimate of the local softening volume (and the local number of neighbours required to perform the interpolation), the Delaunay tessellation in fact poorly estimates the phase-space distribution function. Indeed, it requires, the choice of a global scaling S_G. We propose two methods to estimate S_G that yield a good global compromise. However, the Delaunay interpolation still remains quite sensitive to local anisotropies in the distribution. We also compare 6D phase-space density estimation with the proxy, Q=rho/sigma^3, where rho is the local density and 3 sigma^2 is the local 3D velocity dispersion. We show that Q only corresponds to a rough approximation of the true phase-space density, and is not able to capture all the details of the distribution in phase-space, ignoring, in particular, filamentation and tidal streams.