We explore the vicinity of the Milky Way through the use of spectro-photometric data from the Sloan Digital Sky Survey and high-quality proper motions derived from multi-epoch positions extracted from the Guide Star Catalogue II database. In order to identify and characterise streams as relics of the Milky Way formation, we start with classifying, select, and study $2417$ subdwarfs with $rm{[Fe/H] < -1.5}$ up to $3$ kpc away from the Sun as tracers of the local halo system. Then, through phase-space analysis, we find statistical evidence of five discrete kinematic overdensities among $67$ of the fastest-moving stars, and compare them to high-resolution N-body simulations of the interaction between a Milky-Way like galaxy and orbiting dwarf galaxies with four representative cases of merging histories. The observed overdensities can be interpreted as fossil substructures consisting of streamers torn from their progenitors, such progenitors appear to be satellites on prograde and retrograde orbits on different inclinations. In particular, of the five detected overdensities, two appear to be associated, yelding twenty-one additional main-sequence members, with the stream of Helmi et al. (1999) that our analysis confirms on a high inclination prograde orbit. The three newly identified kinematic groups could be associated with the retrograde streams detected by Dinescu (2002) and Kepley et al. (2007), whatever their origin, the progenitor(s) would be on retrograde orbit(s) and inclination(s) within the range $10^{circ} div 60^{circ}$. Finally, we use our simulations to investigate the impact of observational errors and compare the current picture to the promising prospect of highly improved data expected from the Gaia mission.
Download