We present a novel method for determining the total matter surface density of the Galactic disk by analysing the kinematics of a dynamically cold stellar stream that passes through or close to the Galactic plane. The method relies on the fact that the vertical component of energy for such stream stars is approximately constant, such that their vertical positions and vertical velocities are interrelated via the matter density of the Galactic disk. By testing our method on mock data stellar streams, with realistic phase-space dispersions and Gaia uncertainties, we demonstrate that it is applicable to small streams out to a distance of a few kilo-parsec, and that the surface density of the disk can be determined to a precision of 6 %. This method is complementary to other mass measurements. In particular, it does not rely on any equilibrium assumption for stars in the Galactic disk, and also makes it possible to measure the surface density to good precision at large distances from the Sun. Such measurements would inform us of the matter composition of the Galactic disk and its spatial variation, place stronger constraints on dark disk sub-structure, and even diagnose possible non-equilibrium effects that bias other types of dynamical mass measurements.

The steepness of the central density profiles of dark matter (DM) in low-mass galaxy halos (e.g. dwarf galaxies) is a powerful probe of the nature of DM. We propose a novel scheme to probe the inner profiles of galaxy subhalos using stellar streams. We show that the present day morphological and dynamical properties of accreted globular cluster (GC) streams - those produced from tidal stripping of GCs that initially evolved within satellite galaxies and later merged with the Milky Way (MW) - are sensitive to the central DM density profile and mass of their parent satellites. GCs that accrete within cuspy CDM subhalos produce streams that are physically wider and dynamically hotter than streams that accrete inside cored subhalos. A first comparison of MW streams "GD-1" and "Jhelum" (likely of accreted GC origin) with our simulations indicates a preference for cored subhalos. If these results hold up in future data, the implication is that either the DM cusps were erased by baryonic feedback, or their subhalos naturally possessed cored density profiles implying DM models beyond CDM. Moreover, accreted GC streams are highly structured and exhibit complex morphological features (e.g., parallel structures and "spurs"). This implies that the accretion scenario can naturally explain the recently observed peculiarities in some of the MW streams. We also propose a novel mechanism for forming "gaps" in streams when the remnant of the parent subhalo later passes through the stream. This encounter can last a longer time (and have more of an impact) than the random encounters with DM subhalos previously considered, because the GC stream and its parent subhalo are on similar orbits with small relative velocities. Current and future surveys of the MW halo will uncover numerous faint stellar streams and provide the data needed to substantiate our preliminary tests with this new probe of DM.

We present a new spectroscopic study of the faint Milky Way satellite Sagittarius II. Using multi-object spectroscopy from the Fibre Large Array Multi Element Spectrograph, we supplement the dataset of Longeard et al. (2020) with 47 newly observed stars, 19 of which are identified as members of the satellite. These additional member stars are used to put tighter constraints on the dynamics and the metallicity properties of the system. We find a low velocity dispersion of SgrII v = 1.7 +/- 0.5 km s-1, in agreement with the dispersion of Milky Way globular clusters of similar luminosity. We confirm the very metal-poor nature of the satellite ([Fe/H]_SgrII = -2.23 +/- 0.07) and find that the metallicity dispersion of Sgr II is not resolved, reaching only 0.20 at the 95% confidence limit. No star with a metallicity below -2.5 is confidently detected. Therefore, despite the unusually large size of the system (rh = 35.5 +1.4-1.2 pc), we conclude that Sgr II is an old and metal-poor globular cluster of the Milky Way.

We reconsider the case for the association of Galactic globular clusters (GCs) to the tidal stream of the Sagittarius dwarf spheroidal galaxy (Sgr dSph), using Gaia DR2 data. We use RR Lyrae to trace the stream in 6D and we select clusters matching the observed stream in position and velocity. In addition to the clusters residing in the main body of the galaxy (M 54, Ter 8, Ter 7, Arp 2) we confirm the membership of Pal 12 and Whiting 1 to the portion of the trailing arm populated by stars lost during recent perigalactic passages. NGC 2419, NGC 5634 and NGC 4147 are very interesting candidates, possibly associated to more ancient wraps of the stream. We note that all these clusters, with the exception of M 54, that lies within the stellar nucleus of the galaxy, are found in the trailing arm of the stream. The selected clusters are fully consistent with the [Fe/H] vs. [Mg/Fe], [Ca/Fe] patterns and the age-metallicity relation displayed by field stars in the main body of Sgr dSph.

We present the first full six-dimensional panoramic portrait of the Sagittarius stream, obtained by searching for wide stellar streams in the Gaia DR2 dataset with the STREAMFINDER algorithm. We use the kinematic behavior of the sample to devise a selection of Gaia RR Lyrae, providing excellent distance measurements along the stream. The proper motion data are complemented with radial velocities from public surveys. We find that the global morphological and kinematic properties of the Sagittarius stream are still reasonably well reproduced by the simple Law & Majewski (2010) model (LM10), although the model overestimates the leading arm and trailing arm distances by up to $\sim 15$%. The sample newly reveals the leading arm of the Sagittarius stream as it passes into very crowded regions of the Galactic disk towards the Galactic Anticenter direction. Fortuitously, this part of the stream is almost exactly at the diametrically opposite location from the Galactic Center to the progenitor, which should allow an assessment of the influence of dynamical friction and self-gravity in a way that is nearly independent of the underlying Galactic potential model.

The density variations in thin stellar streams may encode important information on the nature of the dark matter. For instance, if dark matter aggregates into massive sub-halos, these perturbers are expected to scatter stars out of dynamically cold stellar streams, possibly leading to detectable gaps in those structures. Here we re-examine the density variations in the GD-1 stream, using Gaia DR2 astrometry, Pan-STARRS photometry, together with high precision radial velocities measured with the CFHT/ESPaDOnS and VLT/UVES instruments and complemented with public radial velocity catalogs. We show that after correcting for projection effects, the density profile exhibits high contrast periodic peaks, separated by $2.64\pm0.18$ kpc. An N-body simulation is presented that reproduces this striking morphology with simple epicyclic motion in a smooth Galactic potential. We also discuss the reliability of measuring density variations using ground-based photometric surveys, and for the particular case of GD-1 we highlight some of the artifacts present in the Gaia DR2 catalog along its track. Massive dark subhalos do not appear to be required to explain the density clumping along GD-1.

The orbits of the least chemically enriched stars open a window on the formation of our Galaxy when it was still in its infancy. The common picture is that these low-metallicity stars are distributed as an isotropic, pressure-supported component since these stars were either accreted from the early building blocks of the assembling Milky Way, or were later brought by the accretion of faint dwarf galaxies. Combining the metallicities and radial velocities from the Pristine and LAMOST surveys and Gaia DR2 parallaxes and proper motions for an unprecedented large and unbiased sample of very metal-poor stars at $[Fe/H]\leq-2.5$ we show that this picture is incomplete. This sample shows strong statistical evidence (at the $5.0\sigma$ level) of asymmetry in their kinematics, favouring prograde motion. Moreover, we find that $31\%$ of the stars that currently reside in the disk do not venture outside of the disk plane throughout their orbit. The discovery of this population implies that a significant fraction of stars with iron abundances $[Fe/H]\leq-2.5$ formed within or concurrently with the Milky Way disk and that the history of the disk was quiet enough to allow them to retain their disk-like orbital properties.

As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $\sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $\pm$ 0.5) Gyr ago, with a peak rate of (8.8 $\pm$ 1.4) M$_{\odot}$yr$^{-1}$ at (9.8 $\pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $\pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs.

We report the discovery of a 70-degree long stellar stream in the Milky Way halo, which criss-crosses the well known "GD-1" stream. We show that this new stellar structure ("Kshir") and GD-1 lie at similar distance, and are remarkably correlated in kinematics. We propose several explanations for the nature of this new structure and its possible association with GD-1. However, a scenario in which these two streams were accreted onto the Milky Way within the same dark matter sub-halo seems to provide a natural explanation for their phase-space entanglement, and other complexities of this coupled-system.

The rotation curve of the Milky Way is commonly used to estimate the local dark matter density $\rho_{{\rm DM},\odot}$. However, the estimates are subject to the choice of the distribution of baryons needed in this type of studies. In this work we explore several Galactic mass models that differ in the distribution of baryons and dark matter, in order to determine $\rho_{{\rm DM},\odot}$. For this purpose we analyze the precise circular velocity curve measurement of the Milky Way up to $\sim 25$ kpc from the Galactic centre obtained from Gaia DR2 [1]. We find that the estimated value of $\rho_{{\rm DM},\odot}$ stays robust to reasonable changes in the spherical dark matter halo. However, we show that $\rho_{{\rm DM},\odot}$ is affected by the choice of the model for the underlying baryonic components. In particular, we find that $\rho_{{\rm DM},\odot}$ is mostly sensitive to uncertainties in the disk components of the Galaxy. We also show that, when choosing one particular baryonic model, the estimate of $\rho_{{\rm DM},\odot}$ has an uncertainty of only about $10\%$ of its best-fit value, but this uncertainty gets much bigger when we also consider the variation of the baryonic model. In particular, the rotation curve method does not allow to exclude the presence of an additional very thin component, that can increase $\rho_{{\rm DM},\odot}$ by more than a factor of 8 (the thin disk could even be made of dark matter). Therefore, we conclude that exclusively using the rotation curve of the Galaxy is not enough to provide a robust estimate of $\rho_{{\rm DM},\odot}$. For all the models that we study without the presence of an additional thin component, our resulting estimates of the local dark matter density take values in the range $\rho_{{\rm DM},\odot} \simeq \text{0.3--0.4}\,\mathrm{GeV/cm^3}$, consistent with many of the estimates in the literature.

A chemo-dynamical analysis of 115 metal-poor candidate stars selected from the narrow-band Pristine photometric survey is presented based on CFHT high-resolution ESPaDoNS spectroscopy. We have discover 28 new bright (V < 15) stars with [Fe/H]<-2.5 and 5 with [Fe/H]<-3.0 for success rates of 40% (28/70) and 19% (5/27), respectively. A detailed model atmospheres analysis is carried out for the 28 new metal-poor stars. Stellar parameters were determined from SDSS photometric colours, Gaia DR2 parallaxes, MESA/MIST stellar isochrones, and the initial Pristine survey metallicities, following a Bayesian inference method. Chemical abundances are determined for 10 elements (Na, Mg, Ca, Sc, Ti, Cr, Fe, Ni, Y, Ba). Most stars show chemical abundance patterns that are similar to the normal metal-poor stars in the Galactic halo; however, we also report the discoveries of a new r-process rich star, a new CEMP-s candidate with [Y/Ba]>0, and a metal-poor star with very low [Mg/Fe]. The kinematics and orbits for all of the highly probable metal-poor candidates are determined by combining our precision radial velocities with Gaia DR2 proper motions. Some stars show unusual kinematics for their chemistries, including planar orbits, unbound orbits, and highly elliptical orbits that plunge deeply into the Galactic bulge (Rperi < 0.5 kpc); also, eight stars have orbital energies and actions consistent with the Gaia-Enceladus accretion event. This paper contributes to our understanding of the complex chemo-dynamics of the metal-poor Galaxy, and increases the number of known bright metal-poor stars available for detailed nucleosynthetic studies.

Our Galaxy is known to contain a central boxy/peanut-shaped bulge, yet the importance of a classical, pressure-supported component within the central part of the Milky Way is still being debated. It should be most visible at low metallicity, a regime that has not yet been studied in detail. Using metallicity-sensitive narrow-band photometry, the Pristine Inner Galaxy Survey (PIGS) has collected a large sample of metal-poor ([Fe/H] < -1.0) stars in the inner Galaxy to address this open question. We use PIGS to trace the metal-poor inner Galaxy kinematics as function of metallicity for the first time. We find that the rotational signal decreases with decreasing [Fe/H], until it becomes negligible for the most metal-poor stars. Additionally, the velocity dispersion increases with decreasing metallicity for -3.0 < [Fe/H] < -0.5, with a gradient of -44 $\pm$ 4 km$\,$s$^{-1}\,$dex$^{-1}$. These observations may signal a transition between Galactic components of different metallicities and kinematics, a different mapping onto the boxy/peanut-shaped bulge for former disk stars of different metallicities and/or the secular dynamical and gravitational influence of the bar on the pressure-supported component. Our results provide strong constraints on models that attempt to explain the properties of the inner Galaxy.

We present radial velocities for five member stars of the recently discovered young (age= 100-150 Myr) stellar system Price-Whelan 1 (PW 1), that is located far away in the Galactic Halo (D~ 29 kpc, Z~ 15 kpc), and that is probably associated to the Leading Arm (LA) of the Magellanic Stream. We measure the systemic radial velocity of PW 1, V_r=275 +/- 10 km/s, significantly larger than the velocity of the LA gas in the same direction. We re-discuss the main properties and the origin of this system in the light of these new observations, computing the orbit of the system and comparing its velocity with that of the HI in its surroundings. We show that the bulk of the gas at the velocity of the stars is more than 10 deg (5~kpc) away from PW 1 and the velocity difference between the gas and the stars become larger as gas closer to the stars is considered. We discuss the possibilities that (a) the parent gas cloud was dissolved by the interaction with the Galactic gas, and (b) that the parent cloud is the high velocity cloud HVC287.5+22.5+240, lagging behind the stellar system by ~25 km/s and ~10 deg.

Tidally disrupted globular cluster streams are usually observed, and therefore perceived, as narrow, linear and one-dimensional structures in the 6D phase-space. Here we show that the GD-1 stellar stream ($\approx$ 30 pc wide), which is the tidal debris of a disrupted globular cluster, possesses a secondary diffuse and extended stellar component ($\approx$ 100 pc wide) around it, detected at >5$\sigma$ confidence level. Similar morphological properties are seen in synthetic streams that are produced from star clusters that are formed within dark matter sub-halos and then accrete onto a massive host galaxy. This lends credence to the idea that the progenitor of the highly retrograde GD-1 stream was originally formed outside of the Milky Way in a now defunct dark satellite galaxy. We deem that in future studies, this newly found $cocoon$ component may serve as a structural hallmark to distinguish between the in-situ and ex-situ (accreted) formed globular cluster streams.

We use ESA/Gaia astrometry together with SEGUE and LAMOST measurements of the GD-1 stellar stream to explore the improvement on the Galactic gravitational potential that these new data provide. Assuming a realistic universal model for the dark matter halo together with reasonable models of the baryonic components, we find that the orbital solutions for GD-1 require the circular velocity at the Solar radius to be $V_{\rm circ}(R_\odot) =244 \pm 4\,{\rm km\, s^{-1}}$, and also that the density flattening of the dark halo is $q_{\rho}=0.82^{+0.25}_{-0.13}$. The corresponding Galactic mass within 20kpc was estimated to be $M_{\rm MW}(<20kpc)=2.5 \pm 0.2 \times 10^{11}\, M_\odot$. Moreover, Gaia's excellent proper motions also allowed us to constrain the velocity dispersion of the GD-1 stream in the direction tangential to the line of sight to be $<2.30\,{\rm km\, s^{-1}}$ (95 % confidence limit), confirming the extremely cold dynamical nature of this system.

Omega Centauri is the most massive globular cluster of the Milky Way, and possesses many peculiar properties. In particular, the cluster contains distinct multiple stellar populations, with a large spread in metallicity and different kinematics as a function of light elements abundance, implying that it formed over an extended period of time. This has lead to the suggestion that $\omega$ Cen is the remnant core of an accreted dwarf galaxy. If this scenario is correct, $\omega$ Cen should be tidally limited, and one should expect to find tidal debris spread along its orbit. Here we use N-body simulations to show that the recently-discovered `Fimbulthul' structure, identified in the second data release (DR2) of the Gaia mission, is the long sought-for tidal stream of $\omega$ Cen, extending up to $28\deg$ from the cluster. Follow-up high-resolution spectroscopy of 5 stars in the stream show that they are closely-grouped in velocity, and have metallicities consistent with having originated in that cluster. Guided by our N-body simulations, we devise a selection filter that we apply to Gaia data to also uncover the portion of the stream in the highly-contaminated and crowded field within $10\deg$ of $\omega$ Cen. Further modelling of the stream may help to constrain the dynamical history of the dwarf galaxy progenitor of this disrupting system and guide future searches for its remnant stars in the Milky Way.

We present the discovery of a large population of stellar streams that surround the inner Galaxy, found in the Gaia DR2 catalog using the new STREAMFINDER algorithm. Here we focus on the properties of eight new high-significance structures found at Heliocentric distances between 1 and 10 kpc and at Galactic latitudes $|b|>20\deg$, named Slidr, Sylgr, Ylgr, Fimbulthul, Sv\"ol, Fj\"orm, Gj\"oll and Leiptr. Spectroscopic measurements of seven of the streams confirm the detections, which are based on Gaia astrometry and photometry alone, and show that these streams are predominantly metal-poor. The sample possesses diverse orbital properties, although most of the streams appear to be debris of inner-halo globular clusters. Many more candidate streams are visible in our maps, but require follow-up spectroscopy to confirm their nature. We also explain in detail the workings of the algorithm, and gauge the incidence of false detections by running the algorithm on a smooth model of the Gaia catalog.

We present the kinematics of Blue Straggler (BS) stars identified in the Canada-France-Imaging-Survey (CFIS), covering 4000 deg$^2$ on the sky in the $u$-band. The BS sample, characterised through CFIS and Pan-STARRS photometry, has been kinematically decomposed into putative halo and disc populations after cross-matching with Gaia astrometry and SDSS/SEGUE/LAMOST spectroscopy. This decomposition clearly reveals the strong flaring of the outer Milky Way disc. In particular, we show that we can detect this flaring up to a vertical height of $|Z| \simeq 8$ kpc at a Galactocentric distance of $R\sim $27 kpc. While some small level of flaring is expected for extended discs built up by radial migration, we demonstrate that this observed very strong flaring of the Milky Way disc is more likely consistent with it being dynamically heated by the repeated passage of the Sagittarius dwarf spheroidal galaxy through the midplane.

We present the stellar density profile of the outer halo of the Galaxy traced over a range of Galactocentric radii from $15< R_{GC} < 220$ kpc by blue horizontal branch (BHB) stars. These stars are identified photometrically using deep $u-$band imaging from the new Canada-France-Imaging-Survey (CFIS) that reaches 24.5 mag. This is combined with $griz$ bands from Pan-STARRS 1 and covers a total of $\sim4000$ deg$^2$ of the northern sky. We present a new method to select BHB stars that has low contamination from blue stragglers and high completeness. We use this sample to measure and parameterize the three dimensional density profile of the outer stellar halo. We fit the profile using (i) a simple power-law with a constant flattening (ii) a flattening that varies as a function of Galactocentric radius (iii) a broken power law profile. We find that outer stellar halo traced by the BHB is well modelled by a broken power law with a constant flattening of $q=0.86 \pm 0.02$, with an inner slope of $\gamma=4.24 \pm 0.08$. This is much steeper than the preferred outer profile that has a slope of $\beta=3.21\pm 0.07$ after a break radius of $r_b=41.4^{+2.5}_{-2.4}$ kpc. The outer profile of the stellar halo trace by BHB stars is shallower than that recently measured using RR Lyrae, a surprising result given the broad similarity of the ages of these stellar populations.

We report the discovery of a $75\deg$ long stellar stream in Gaia DR2 catalog, found using the new STREAMFINDER algorithm. The structure is probably the remnant of a now fully disrupted globular cluster, lies $\approx 3.8 \, {\rm kpc}$ away from the Sun in the direction of the Galactic bulge, and possesses highly retrograde motion. We find that the system orbits close to the Galactic plane at Galactocentric distances between $4.9$ and $19.8 \, {\rm kpc}$. The discovery of this extended and extremely low surface brightness stream ($\Sigma_G\sim 34.3 \, {\rm mag \, arcsec^{-2}}$) with a mass of only $2580\pm140 \, {\rm\,M_\odot}$, demonstrates the power of the STREAMFINDER algorithm to detect even very nearby and ultra-faint structures. Due to its proximity and length we expect that Phlegethon will be a very useful probe of the Galactic acceleration field.

We present a panoramic map of the stellar streams of the Milky Way based upon astrometric and photometric measurements from the Gaia DR2 catalogue. In this first contribution, we concentrate on the halo at heliocentric distances beyond 5 kpc, and at Galactic latitudes $|b|>30\deg$, using the STREAMFINDER algorithm to detect structures along plausible orbits that are consistent with the Gaia proper motion measurements. We find a rich network of criss-crossing streams in the halo. Some of these structures were previously-known, several are new discoveries, but others are potentially artefacts of the Gaia scanning law and will require confirmation. With these initial discoveries, we are starting to unravel the complex formation of the halo of our Galaxy.

STREAMFINDER is a new algorithm that we have built to detect stellar streams in an automated and systematic way in astrophysical datasets that possess any combination of positional and kinematic information. In Paper I, we introduced the methodology and the workings of our algorithm and showed that it is capable of detecting ultra-faint and distant halo stream structures containing as few as $\sim 15$ members ($\Sigma_{\rm G} \sim 33.6\, {\rm mag \, arcsec^{-2}}$) in the Gaia dataset. Here, we test the method with real proper motion data from the Pan-STARRS1 survey, and by selecting targets down to $r_{0}=18.5$ mag we show that it is able to detect the GD-1 stellar stream, whereas the structure remains below a useful detection limit when using a Matched Filter technique. The radial velocity solutions provided by STREAMFINDER for GD-1 candidate members are found to be in good agreement with observations. Furthermore, our algorithm detects a $\sim 40^{\deg}$ long structure approximately parallel to GD-1, and which fans out from it, possibly a sign of stream-fanning due to the triaxiality of the Galactic potential. This analysis shows the promise of this method for detecting and analysing stellar streams in the upcoming Gaia DR2 catalogue.

We have designed a powerful new algorithm to detect stellar streams in an automated and systematic way. The algorithm, which we call the STREAMFINDER, is well suited for finding dynamically cold and thin stream structures that may lie along any simple or complex orbits in Galactic stellar surveys containing any combination of positional and kinematic information. In the present contribution we introduce the algorithm, lay out the ideas behind it, explain the methodology adopted to detect streams and detail its workings by running it on a suite of simulations of mock Galactic survey data of similar quality to that expected from the ESA/Gaia mission. We show that our algorithm is able to detect even ultra-faint stream features lying well below previous detection limits. Tests show that our algorithm will be able to detect distant halo stream structures $>10^{\circ}$ long containing as few as $\sim 15$ members ($\Sigma_{\rm G} \sim 33.6\, {\rm mag \, arcsec^{-2}}$) in the Gaia dataset.

The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies.

We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30\deg < |b| < 70\deg$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger.

We present a method for measuring the Sun's motion using the proper motions of Galactic halo star streams. The method relies on the fact that the motion of the stars perpendicular to a stream from a low-mass progenitor is close to zero when viewed from a non-rotating frame at rest with respect to the Galaxy, and that the deviation from zero is due to the reflex motion of the observer. The procedure we implement here has the advantage of being independent of the Galactic mass distribution. We run a suite of simulations to test the algorithm we have developed, and find that we can recover the input Solar motion to good accuracy with data of the quality that will soon become available from the ESA/Gaia mission.