Metal-poor stars are important tools for tracing the early history of the Milky Way, and for learning about the first generations of stars. Simulations suggest that the oldest metal-poor stars are to be found in the inner Galaxy. Typical bulge surveys, however, lack low metallicity ([Fe/H] < -1.0) stars because the inner Galaxy is predominantly metal-rich. The aim of the Pristine Inner Galaxy Survey (PIGS) is to study the metal-poor and very metal-poor (VMP, [Fe/H] < -2.0) stars in this region. In PIGS, metal-poor targets for spectroscopic follow-up are selected from metallicity-sensitive CaHK photometry from the CFHT. This work presents the ~250 deg^2 photometric survey as well as intermediate-resolution spectroscopic follow-up observations for ~8000 stars using AAOmega on the AAT. The spectra are analysed using two independent tools: ULySS with an empirical spectral library, and FERRE with a library of synthetic spectra. The comparison between the two methods enables a robust determination of the stellar parameters and their uncertainties. We present a sample of 1300 VMP stars -- the largest sample of VMP stars in the inner Galaxy to date. Additionally, our spectroscopic dataset includes ~1700 horizontal branch stars, which are useful metal-poor standard candles. We furthermore show that PIGS photometry selects VMP stars with unprecedented efficiency: 86%/80% (lower/higher extinction) of the best candidates satisfy [Fe/H] < -2.0, as do 80%/63% of a larger, less strictly selected sample. We discuss future applications of this unique dataset that will further our understanding of the chemical and dynamical evolution of the innermost regions of our Galaxy.

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 Pristine survey uses narrow-band photometry to derive precise metallicities down to the extremely metal-poor regime ([Fe/H] < -3), and currently consists of over 4 million FGK-type stars over a sky area of $\sim 2~500\, \mathrm{deg}^2$. We focus our analysis on a subsample of $\sim 80~000$ main sequence turnoff stars with heliocentric distances between 6 and 20 kpc, which we take to be a representative sample of the inner halo. The resulting metallicity distribution function (MDF) has a peak at [Fe/H] = -1.6, and a slope of $\Delta$(LogN)/$\Delta[Fe/H] = 1.0 \pm 0.1$ in the metallicity range of -3.4 < [Fe/H] < -2.5. This agrees well with a simple closed-box chemical enrichment model in this range, but is shallower than previous spectroscopic MDFs presented in the literature, suggesting that there may be a larger proportion of metal-poor stars in the inner halo than previously reported. We identify the Monoceros/TriAnd/ACS/EBS/A13 structure in metallicity space in a low latitude field in the anticenter direction, and also discuss the possibility that the inner halo is dominated by a single, large merger event, but cannot strongly support or refute this idea with the current data. Finally, based on the MDF of field stars, we estimate the number of expected metal-poor globular clusters in the Milky Way halo to be 5.4 for [Fe/H] < -2.5 and 1.5 for [Fe/H] < -3, suggesting that the lack of low metallicity globular clusters in the Milky Way is not due simply to statistical undersampling.

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.

(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

We use the Pristine survey CaHK narrow-band photometry, combined with the SDSS ugr photometry, to provide a cleaner sample of blue horizontal branch stars in the Galactic halo out to large distances. We demonstrate a completeness of 91% and a purity of 93% with respect to available spectroscopic classifications. We subsequently use our new clean sample of these standard candles to investigate the substructure in the Galactic halo over the Pristine footprint. Among other features, this allows for a careful tracing of multiple parts of the Sagittarius stream, providing a measurement independent from other tracers used and reaching larger distances. Moreover, we demonstrate with this clean and complete sample that the halo follows a density profile with a negative power-law slope of 3.5 - 4.0. As the relatively shallow SDSS u-band is the limiting factor in this technique, we foresee large potential for combining Pristine survey photometry with the much deeper u-band photometry from the Canada-France-Imaging Survey.

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 a new fully data-driven algorithm that uses photometric data from the Canada-France-Imaging-Survey (CFIS; $u$), Pan-STARRS 1 (PS1; $griz$), and Gaia ($G$) to discriminate between dwarf and giant stars and to estimate their distances and metallicities. The algorithm is trained and tested using the SDSS/SEGUE spectroscopic dataset and Gaia photometric/astrometric dataset. At [Fe/H]$<-1.2$, the algorithm succeeds in identifying more than 70% of the giants in the training/test set, with a dwarf contamination fraction below 30% (with respect to the SDSS/SEGUE dataset). The photometric metallicity estimates have uncertainties better than 0.2 dex when compared with the spectroscopic measurements. The distances estimated by the algorithm are valid out to a distance of at least $\sim 80$ kpc without requiring any prior on the stellar distribution, and have fully independent uncertainities that take into account both random and systematic errors. These advances allow us to estimate these stellar parameters for approximately 12 million stars in the photometric dataset. This will enable studies involving the chemical mapping of the distant outer disc and the stellar halo, including their kinematics using the Gaia proper motions. This type of algorithm can be applied in the Southern hemisphere to the first release of LSST data, thus providing an almost complete view of the external components of our Galaxy out to at least $\sim 80$ kpc. Critical to the success of these efforts will be ensuring well-defined spectroscopic training sets that sample a broad range of stellar parameters with minimal biases. A catalogue containing the training/test set and all relevant parameters within the public footprint of CFIS is available online.

We present the results of a 3-year long, medium-resolution spectroscopic campaign aimed at identifying very metal-poor stars from candidates selected with the CaHK, metallicity-sensitive Pristine survey. The catalogue consists of a total of 1007 stars, and includes 146 rediscoveries of metal-poor stars already presented in previous surveys, 707 new very metal-poor stars with [Fe/H] < -2.0, and 95 new extremely metal-poor stars with [Fe/H] < -3.0. We provide a spectroscopic [Fe/H] for every star in the catalogue, and [C/Fe] measurements for a subset of the stars (10% with [Fe/H] < -3 and 24% with -3 < [Fe/H] < -2) for which a carbon determination is possible, contingent mainly on the carbon abundance, effective temperature and S/N of the stellar spectra. We find an average carbon enhancement fraction ([C/Fe] >= +0.7) of 41 +- 4% for stars with -3 < [Fe/H] < -2 and 58 +- 14% for stars with [Fe/H] < -3, and report updated success rates for the Pristine survey of 56 % and 23 % to recover stars with [Fe/H] < -2.5 and [Fe/H] < -3, respectively. Finally, we discuss the current status of the survey and its preparation for providing targets to upcoming multi-object spectroscopic surveys such as WEAVE.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics.

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.

The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products.

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.

Kinematically cold tidal streams of globular clusters (GC) are excellent tracers of the Galactic gravitational potential at moderate Galactocentric distances, and can also be used as probes of the law of gravity on Galactic scales. Here, we compare for the first time the generation of such streams in Newtonian and Milgromian gravity (MOND). We first compute analytical results to investigate the expected shape of the GC gravitational potential in both frameworks, and we then run N-body simulations with the Phantom of Ramses code. We find that the GCs tend to become lopsided in MOND. This is a consequence of the external field effect which breaks the strong equivalence principle. When the GC is filling its tidal radius the lopsidedness generates a strongly asymmetric tidal stream. In Newtonian dynamics, such markedly asymmetric streams can in general only be the consequence of interactions with dark matter subhalos, giant molecular clouds, or interaction with the Galactic bar. In these Newtonian cases, the asymmetry is the consequence of a very large gap in the stream, whilst in MOND it is a true asymmetry. This should thus allow us in the future to distinguish these different scenarios by making deep observations of the environment of the asymmetric stellar stream of Palomar 5. Moreover, our simulations indicate that the high internal velocity dispersion of Palomar 5 for its small stellar mass would be natural in MOND.

We present two spectroscopic surveys of the tidal stellar stream of the Palomar 5 globular cluster, undertaken with the VLT/FLAMES and AAT/AAOmega instruments. We use these data in conjunction with photometric data presented in the previous contribution in this series to classify the survey stars in terms of their probability of belonging to the Palomar 5 stellar stream. We find that high-probability candidates are only found in a very narrow spatial interval surrounding the locus of the stream on the sky. PanSTARRS RRLyrae stars in this region of sky are also distributed in a similar manner. The absence of significant "fanning" of this stellar stream confirms that Palomar 5 does not follow a chaotic orbit. Previous studies have found that Palomar 5 is largely devoid of low-mass stars, and we show that this is true also of the stellar populations along the trailing arm out to $6\deg$. Within this region, which contains 73\% of the detected stars, the population is statistically identical to the core, implying that the ejection of the low-mass stars occurred before the formation of the stream. We also present an updated structural model fit to the bound remnant, which yields a total mass of $4297\pm98 {\rm\,M_\odot}$ and a tidal radius $0.145\pm0.009$ kpc. We estimate the mass of the observed system including the stream to be $12200\pm400 {\rm\,M_\odot}$, and the initial mass to have been $\sim47000\pm1500 {\rm\,M_\odot}$. These observational constraints will be employed in our next study to model the dynamics of the system in detail.

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.

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.

Tidal streams of disrupting dwarf galaxies orbiting around their host galaxy offer a unique way to constrain the shape of galactic gravitational potentials. Such streams can be used as leaning tower gravitational experiments on galactic scales. The most well motivated modification of gravity proposed as an alternative to dark matter on galactic scales is Milgromian dynamics (MOND), and we present here the first ever N-body simulations of the dynamical evolution of the disrupting Sagittarius dwarf galaxy in this framework. Using a realistic baryonic mass model for the Milky Way, we attempt to reproduce the present-day spatial and kinematic structure of the Sagittarius dwarf and its immense tidal stream that wraps around the Milky Way. With very little freedom on the original structure of the progenitor, constrained by the total luminosity of the Sagittarius structure and by the observed stellar mass-size relation for isolated dwarf galaxies, we find reasonable agreement between our simulations and observations of this system. The observed stellar velocities in the leading arm can be reproduced if we include a massive hot gas corona around the Milky Way that is flattened in the direction of the principal plane of its satellites. This is the first time that tidal dissolution in MOND has been tested rigorously at these mass and acceleration scales.

We present an analysis of the presence of substructures in the stellar stream of the Palomar 5 globular cluster, as derived from Sloan Digital Sky Survey data. Using a matched filter technique, we recover the positions and sizes of overdensities reported in previous studies. To explore the reality of these structures, we also create an artificial model of the stream, in which we construct a realistic background on top of which we add a perfectly smooth stream structure, taking into account the effects of photometric completeness and interstellar extinction. We find that the smooth artificial stream then shows similarly-pronounced substructures as the real structure. Interestingly, our best-fit N-body simulation does display real projected density variations linked to stellar epicyclic motions, but these become less significant when taking into account the SDSS star-count constraints. The substructures found when applying our matched filter technique to the N-body particles converted into observable stars are thus mostly unrelated to these epicyclic motions. This analysis suggests that the majority of the previously-detected substructures along the tidal tail of Palomar 5 are artefacts of observational inhomogeneities.

(abridged) We report eclipse mapping of time resolved spectroscopy of the dwarf nova IP Pegasi on the late decline from the May 1993 outburst. The spiral arms are still visible some 8 days after the onset of the outburst. Their fractional contribution to the continuum emission is similar to that measured close to outburst maximum, whereas their orientation is rotated by 58 degrees with respect to the spirals seen in the eclipse map at outburst maximum. Velocity-resolved light curves across the Halpha and the HeI lines show the classical rotational disturbance, with the blue side of the line being eclipsed earlier than the red side. The differences between the Halpha and the HeI maps are significant. The spiral arms are clearly seen in the HeI maps, with the receding arm being stronger in the red side while the approaching arm is stronger in the blue side of the line. The analysis of the Halpha maps suggests that this emission arises mainly from a large and vertically-extended region which we interpret as an outflowing (and spiraling) disc wind. The spectrum of the uneclipsed light is dominated by a strong, blueshifted and narrow Halpha emission line superimposed on a red continuum and can be understood as a combination of emission from an M5V secondary star plus optically thin emission from the outer parts of the vertically-extended disc wind. The inner disc regions show an emission line spectrum with a strong and broad Halpha component superimposed on a flat continuum. This is in marked contrast with the results from the spectral mapping of nova-like variables of comparable inclination and mass ratio and suggests that intrinsically different physical conditions hold in the inner disc regions of outbursting dwarf novae and nova-like systems.