A large number of new members ($\sim$150) of the Cetus Stream (CS) were identified from their clustering features in dynamical space using 6D kinematic data by combining LAMOST DR5 and Gaia DR2 surveys. They map a diffuse structure that extends over at least 100 degrees in the northern and southern Galactic hemispheres, at heliocentric distances between 20 to 50 kpc. Taking advantage of this expanded dataset, we model the stream with a suite of tailored N-body simulations. Our findings exclude the possibility that the NGC 5824 globular cluster is the core of the progenitor of the stream, as postulated by previous studies. Our best models, which successfully reproduce the features of the CS indicate that the progenitor is likely a dwarf galaxy of $\sim$ 2$\times$10$^9$M$_{\odot}$, with a diffuse disc morphology. The merger occured $\sim$ 5 Gyr ago and since then it has experienced approximately eight apo-center passages. Our results suggest that NGC 5824 was either a globular cluster situated off-centre in the dwarf progenitor or, alternatively, it was the nuclear star cluster of another dwarf galaxy that has very similar orbit as the progenitor of the CS. In both scenarios, the progenitor systems would leave streams around NGC 5824, but with distinct distance distributions. To discriminate between these scenarios, the detection and accurate distance measurements of the predicted stream around the GC are crucial, which will be possible in the upcoming LSST era. Our simulations also predict that part of the Southern Cetus stream is very likely the newly discovered Palca stream, and possibly related to another, more diffuse Southern substructure, the Eridanus-Pheonix overdensity.

We use a sample of ~350 RR Lyrae stars with radial velocities and Gaia DR2 proper motions to study orbital properties of the Hercules-Aquila Cloud (HAC) and Virgo Over-density (VOD). We demonstrate that both structures are dominated by stars on highly eccentric orbits, with peri-centres around ~1 kpc and apo-centres between 15 and 25 kpc from the Galactic centre. Given that the stars in the HAC and the VOD occupy very similar regions in the space spanned by integrals of motion, we conclude that these diffuse debris clouds are part of the same accretion event. More precisely, these inner halo sub-structures likely represent two complementary not-fully-mixed portions of an ancient massive merger, also known as the "sausage" event.

We designed a follow-up program to find the spectroscopic properties of the Hercules-Aquila Cloud (HAC) and test scenarios for its formation. We measured the radial velocities (RVs) of 45 RR Lyrae in the southern portion of the HAC using the facilities at the MDM observatory, producing the first large sample of velocities in the HAC. We found a double-peaked distribution in RVs, skewed slightly to negative velocities. We compared both the morphology of HAC projected onto the plane of the sky and the distribution of velocities in this structure outlined by RR Lyrae and other tracer populations at different distances to N-body simulations. We found that the behaviour is characteristic of an old, well-mixed accretion event with small apo-galactic radius. We cannot yet rule out other formation mechanisms for the HAC. However, if our interpretation is correct, HAC represents just a small portion of a much larger debris structure spread throughout the inner Galaxy whose distinct kinematic structure should be apparent in RV studies along many lines of sight.

We study the structure of the inner Milky Way using the latest data release of the Vista Variables in Via Lactea (VVV) survey. The VVV is a deep near-infrared, multi-colour photometric survey with a coverage of 300 square degrees towards the Bulge/Bar. We use Red Clump (RC) stars to produce a high-resolution dust map of the VVV's field of view. From de-reddened colour-magnitude diagrams we select Red Giant Branch stars to investigate their 3D density distribution within the central 4 kpc. We demonstrate that our best-fit parametric model of the Bulge density provides a good description of the VVV data, with a median percentage residual of 5$\%$ over the fitted region. The strongest of the otherwise low-level residuals are overdensities associated with a low-latitude structure as well as the so-called X-shape previously identified using the split RC. These additional components contribute only $\sim5\%$ and $\sim7\%$ respectively to the Bulge mass budget. The best-fit Bulge is `boxy' with an axis ratio of [1:0.44:0.31] and is rotated with respect to the Sun-Galactic Centre line by at least $20^{\circ}$. We provide an estimate of the total, full sky, mass of the Bulge of $M_\mathrm{Bulge}^{\mathrm{Chabrier}} = 2.36 \times 10^{10} M_{\odot}$ for a Chabrier initial mass function. We show there exists a strong degeneracy between the viewing angle and the dispersion of the Red Clump absolute magnitude distribution. The value of the latter is strongly dependent on the assumptions made about the intrinsic luminosity function of the Bulge.

We map the large-scale sub-structure in the Galactic stellar halo using accurate 3D positions of ~14,000 RR Lyrae reported by the Catalina Sky Survey. In the heliocentric distance range of 10-25 kpc, in the region of the sky approximately bounded by 30{\deg} < l < 55{\deg} and -45{\deg} < b < -25{\deg}, there appears to be a strong excess of RRab stars. This overdensity, peaking at 18 kpc, is most likely associated with the so-called Hercules-Aquila Cloud, previously detected using Main Sequence tracers at similar distances in the Sloan Digital Sky Survey data. Our analysis of the period-amplitude distribution of RR Lyrae in this region indicates that the HAC is dominated by the Oosterhoff I type population. By comparing the measured RR Lyrae number density to models of a smooth stellar halo, we estimate the significance of the observed excess and provide an updated estimate of the total luminosity of the Cloud's progenitor.