We identify shell structures in the Milky Way for the first time. We find 2 shells in the Virgo Overdensity (VOD) region and 2 shells in the Hercules Aquila Cloud (HAC) region using Sloan Digital Sky Survey and Gaia data. These shell stars are a subset of the substructure previously identified as the Virgo Radial Merger (VRM). Timing arguments for these shells indicate that their progenitor dwarf galaxy passed through the Galactic center 2.4 $\pm$ 0.2 Gyr ago. Based on the time of collision, it is also possible that the VRM is related to the phenomenon that created phase-space spirals in the vertical motion of the disk and the Splash, and could have caused a burst of star formation in the inner disk. We analyze phase mixing in a collection of radial merger N-body simulations, and find that shell structure similar to that observed in Milky Way data disappears by 5 Gyr after collision with the Galactic center. The method used to calculate the merger time of the VRM was able to reliably recover the correct merger times for these simulations. Previous work supports the idea that the VRM and the Gaia Sausage/Gaia-Enceladus Merger are the same. However, the Gaia Sausage is widely believed to be 8--11 Gyr old. The disparate ages could be reconciled if the larger age is associated with an infall time when the progenitor crossed the virial radius; we do not constrain the time at which the progenitor became bound to the Milky Way. Alternatively, the Gaia Sausage could be younger than previously thought.

We suggest that the Virgo Overdensity (VOD) of stars in the stellar halo is the result of a radial dwarf galaxy merger that we call the Virgo Radial Merger. Because the dwarf galaxy passed very near to the Galactic center, the debris has a large range of energies but nearly zero $L_z$ angular momentum. The debris appears to extend from 5 to 50 kpc from the Sun in the Virgo region. We connect different parts of this merger debris to the Perpendicular and Parallel Streams (the Virgo Stellar Stream is associated with either or both of these streams), the Hercules-Aquila Cloud (HAC), and possibly the Eridanus Phoenix Overdensity (EriPhe). This radial merger can explain the majority of the observed moving groups of RR Lyrae and blue horizontal branch stars that have previously been identified in Virgo. This merger also produces debris in the Solar neighborhood similar to that identified as the Gaia-Enceladus or it Gaia-sausage merger. Orbits are provided for components of the Virgo Radial Merger progenitor and for debris that appears to be related to the Cocytos Stream, which was also recovered in the Virgo region.

Based on the second Gaia data release (Gaia DR2) and spectroscopy from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data, we identified 20,089 halo stars kinematically and chemically. The halo streams in the solar neighborhood could be detected in the space of energy and angular momentum. We reshuffle the velocities of these stars to determine the significance of substructure. Finally, we identify 4 statistically significant substructures that are labeled GL-1 through 4. Among these substructures, GL-1 is previously known stream ("N2" stream) and the rest 3 substructures are new. These substructures may be the debris of dwarf galaxies accretion event, their dynamical and chemical information can help to understand the history of the Milky Way.

Based on the first Gaia data release and spectroscopy from the LAMOST Data Release 4, we study the kinematics and chemistry of the local halo stars. The halo stars are identified kinematically with a relative speed of at least 220 km s$^{-1}$ with respect to the local standard of rest. In total, 436 halo stars are identified. From this halo sample, 16 high velocity (HiVel) stars are identified. We studied the metallicity and [$\alpha$/Fe] distribution of these HiVel stars. Though most of HiVel stars are metal-poor, there are several stars that have metallicity above $-0.5$ dex. To understand the origin of high velocity stars, we evolve the trajectory of the star backwards along the orbit in our adopted Galaxy potential model to determine the orbital parameters and assess whether the star could have originated in the Galactic center. We found that some high velocity stars could have originated from the Galactic center, but other stars were probably kicked up from the Galactic disk.