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 results from deep and very spatially extended CTIO/DECam $g$ and $r$ photometry (reaching out to $\sim$ 2 mag below the oldest MSTO and covering $\sim$ 20 deg$^2$) around the Sextans dSph. We use this data-set to study the structural properties of Sextans overall stellar population and its different stellar evolutionary phases, as well as to search for signs of tidal disturbance from the MW, which would indicate departure from dynamical equilibrium. We perform the most accurate structural analysis to-date of Sextans' stellar components by applying Bayesian MCMC methods to the individual stars' positions. Surface density maps are built by decontaminating the sample through a matched filter analysis of the CMD, and then analysed for departures from axisymmetry. Sextans is found to be considerably less spatially extended than early studies suggested. No significant distortions or tidal disturbances are found down to a surface brightness of $\sim$ 31.8 mag/arcsec$^{-2}$ in V-band. We identify an overdensity in the central regions that may correspond to previously reported kinematic substructure(s). In agreement with previous findings, old $\&$ metal-poor stars such as BHB stars cover a much larger area than stars in other evolutionary phases, and bright BSs are less spatially extended than faint ones. However, the different spatial distribution of bright and faint BSs appears consistent with the general age/metallicity gradients found in Sextans' stellar component. This is compatible with BSs having formed by evolution of binaries and not necessarily due to the presence of a central disrupted globular cluster, as suggested in the literature. We provide structural parameters for the various populations analyzed and make publicly available the photometric catalogue of point-sources as well as a catalogue of literature spectroscopic measurements with updated membership probabilities.

We present a comprehensive analysis of the structural properties and luminosities of the 23 dwarf spheroidal galaxies that fall within the footprint of the Pan-Andromeda Archaeological Survey (PAndAS). These dwarf galaxies represent the large majority of Andromeda's known satellite dwarf galaxies and cover a wide range in luminosity ($-11.6<M_V<-5.8$ or $10^{4.2}< L <10^{6.5} L_\odot$) and surface brightness ($25.1<\mu_0<29.3$ mag/arcsec$^2$). We confirm most previous measurements, but find And XIX to be significantly larger than before ($r_h=3065^{+1065}_{-935}$ pc, $M_V=-10.1^{+0.8}_{-0.4}$) and cannot derive parameters for And XXVII as it is likely not a bound stellar system. We also significantly revise downward the luminosities of And~XV and And~XVI, which are now $M_V\sim-7.5$ or $L\sim10^5 L_\odot$. Finally, we provide the first detailed analysis of Cas II/And XXX, a fairly faint system ($M_V=-8.0^{+0.4}_{-0.3}$) of typical size ($r_h=270\pm50$ pc), located in close proximity to the two bright elliptical dwarf galaxies NGC 147 & 185. Combined with the set of homogeneous distances published in an earlier contribution, our analysis dutifully tracks all relevant sources of uncertainty in the determination of the properties of the dwarf galaxies from the PAndAS photometric catalogue. We further publish the posterior probability distribution functions of all the parameters we fit for in the form of MCMC chains available online; these inputs should be used in any analysis that aims to remain truthful to the data and properly account for covariance between parameters.

We present a renewed look at M31's Giant Stellar Stream along with the nearby structures Stream C and Stream D, exploiting a new algorithm capable of fitting to the red giant branch (RGB) of a structure in both colour and magnitude space. Using this algorithm, we are able to generate probability distributions in distance, metallicity and RGB width for a series of subfields spanning these structures. Specifically, we confirm a distance gradient of approximately 20 kpc per degree along a 6 degree extension of the Giant Stellar Stream, with the farthest subfields from M31 lying ~ 120 kpc more distant than the inner-most subfields. Further, we find a metallicity that steadily increases from -0.7^{+0.1}_{-0.1} dex to -0.2^{+0.2}_{-0.1} dex along the inner half of the stream before steadily dropping to a value of -1.0^{+0.2}_{-0.2} dex at the farthest reaches of our coverage. The RGB width is found to increase rapidly from 0.4^{+0.1}_{-0.1} dex to 1.1^{+0.2}_{-0.1} dex in the inner portion of the stream before plateauing and decreasing marginally in the outer subfields of the stream. In addition, we estimate Stream C to lie at a distance between 794 and 862 kpc and Stream D between 758 kpc and 868 kpc. We estimate the median metallicity of Stream C to lie in the range -0.7 to -1.6 dex and a metallicity of -1.1^{+0.3}_{-0.2} dex for Stream D. RGB widths for the two structures are estimated to lie in the range 0.4 to 1.2 dex and 0.3 to 0.7 dex respectively. In total, measurements are obtained for 19 subfields along the Giant Stellar Stream, 4 along Stream C, 5 along Stream D and 3 general M31 spheroid fields for comparison. We thus provide a higher resolution coverage of the structures in these parameters than has previously been available in the literature.

NGC147, NGC185 and CassiopeiaII (CassII) have similar positions in the sky, distances and measured line of sight velocities. This proximity in phase space suggests that these three satellites of M31 form a subgroup within the Local Group. Nevertheless, the differences in their star formation history and interstellar medium, and the recent discovery of a stellar stream in NGC~147, combined with the lack of tidal features in the other two satellites, are all indications of complex and diverse interactions between M31 and these three satellites. We use a genetic algorithm to explore the different orbits that these satellites can have and select six sets of orbits that could best explain the observational features of the NGC147, NGC185 and CassII satellites. The parameters of these orbits are then used as a starting point for N-body simulations. We present models for which NGC147, NGC185 and CassII are a bound group for a total time of at least one Gyr but still undergo different interactions with M31 and as a result NGC147 has a clear stellar stream whereas the other two satellites have no significant tidal features. This result shows that it is possible to find solutions that reproduce the contrasting properties of the satellites and for which NGC147-NGC185-CassII have been gravitationally bound.

We present the first detailed analysis of the East Cloud, a highly disrupted diffuse stellar substructure in the outer halo of M31. The core of the substructure lies at a projected distance of $\sim100$ kpc from the centre of M31 in the outer halo, with possible extensions reaching right into the inner halo. Using Pan-Andromeda Archaeological Survey photometry of red giant branch stars, we measure the distance, metallicity and brightness of the cloud. Using Hubble Space Telescope data, we independently measure the distance and metallicity to the two globular clusters coincident with the East Cloud core, PA-57 and PA-58, and find their distances to be consistent with the cloud. Four further globular clusters coincident with the substructure extensions are identified as potentially associated. Combining the analyses, we determine a distance to the cloud of $814^{+20}_{-9}$ kpc, a metallicity of $[Fe/H] = -1.2\pm0.1$, and a brightness of $M_V = -10.7\pm0.4$ mag. Even allowing for the inclusion of the potential extensions, this accounts for less than $20$ per cent of the progenitor luminosity implied by the luminosity-metallicity relation. Using the updated techniques developed for this analysis, we also refine our estimates of the distance and brightness of the South-West Cloud, a separate substructure analyzed in the previous work in this series.

Wrapping around the Milky Way, the Sagittarius stream is the dominant substructure in the halo. Our statistical selection method has allowed us to identify 106 highly likely members of the Sagittarius stream. Spectroscopic analysis of metallicity and kinematics of all members provides us with a new mapping of the Sagittarius stream. We find correspondence between the velocity distribution of stream stars and those computed for a triaxial model of the Milky Way dark matter halo. The Sagittarius trailing arm exhibits a metallicity gradient, ranging from $-0.59$ dex to $-0.97$ dex over 142$^{\circ}$. This is consistent with the scenario of tidal disruption from a progenitor dwarf galaxy that possessed an internal metallicity gradient. We note high metallicity dispersion in the leading arm, causing a lack of detectable gradient and possibly indicating orbital phase mixing. We additionally report on a potential detection of the Sextans dwarf spheroidal in our data.

A central tenet of the current cosmological paradigm is that galaxies grow over time through the accretion of smaller systems. Here, we present new kinematic measurements near the centre of one of the densest pronounced substructures, the South-West Cloud, in the outer halo of our nearest giant neighbour, the Andromeda galaxy. These observations reveal that the kinematic properties of this region of the South-West Cloud are consistent with those of PA-8, a globular cluster previously shown to be co-spatial with the stellar substructure. In this sense the situation is reminiscent of the handful of globular clusters that sit near the heart of the Sagittarius dwarf galaxy, a system that is currently being accreted into the Milky Way, confirming that accretion deposits not only stars but also globular clusters into the halos of large galaxies.

We present deep optical photometry from the DECam imager on the 4m Blanco telescope of over 12 deg$^2$ around the Carina dwarf spheroidal, with complete coverage out to 1 degree and partial coverage extending out to 2.6 degrees. Using a Poisson-based matched filter analysis to identify stars from each of the three main stellar populations, old, intermediate, and young, we confirm the previously identified radial age gradient, distance, tidal radius, stellar radial profiles, relative stellar population sizes, ellipticity, and position angle. We find an angular offset between the three main elliptical populations of Carina, and find only tentative evidence for tidal debris, suggesting that past tidal interactions could not have significantly influenced the Carina dwarf. We detect stars in the vicinity of, but distinct to, the Carina dwarf, and measure their distance to be 46$\pm$2 kpc. We determine this population to be part of the halo of the Large Magellanic Cloud at an angular radius of over 20 degrees. Due to overlap in colour-magnitude space with Magellanic stars, previously detected tidal features in the old population of Carina are likely weaker than previously thought.

We reveal the highly structured nature of the Milky Way stellar halo within the footprint of the PAndAS photometric survey from blue main sequence and main sequence turn-off stars. We map no fewer than five stellar structures within a heliocentric range of ~5 to 30 kpc. Some of these are known (the Monoceros Ring, the Pisces/Triangulum globular cluster stream), but we also uncover three well-defined stellar structures that could be, at least partly, responsible for the so-called Triangulum/Andromeda and Triangulum/Andromeda 2 features. In particular, we trace a new faint stellar stream located at a heliocentric distance of ~17 kpc. With a surface brightness of \Sigma_V ~ 32-32.5 mag/arcsec^2, it follows an orbit that is almost parallel to the Galactic plane north of M31 and has so far eluded surveys of the Milky Way halo as these tend to steer away from regions dominated by the Galactic disk. Investigating our follow-up spectroscopic observations of PAndAS, we serendipitously uncover a radial velocity signature from stars that have colors and magnitudes compatible with the stream. From the velocity of eight likely member stars, we show that this stellar structure is dynamically cold, with an unresolved velocity dispersion that is lower than 7.1 km/s at the 90-percent confidence level. Along with the width of the stream (300-650 pc), its dynamics points to a dwarf-galaxy-accretion origin. The numerous stellar structures we can map in the Milky Way stellar halo between 5 and 30 kpc and their varying morphology is a testament to the complex nature of the stellar halo at these intermediate distances.

In a recent contribution, Bahl \& Baumgardt investigated the incidence of planar alignments of satellite galaxies in the Millennium-II simulation, and concluded that vast thin planes of dwarf galaxies, similar to that observed in the Andromeda galaxy (M31), occur frequently by chance in $\Lambda$-Cold Dark Matter cosmology. However, their analysis did not capture the essential fact that the observed alignment is simultaneously radially extended, yet thin, and kinematically unusual. With the caveat that the Millennium-II simulation may not have sufficient mass resolution to identify confidently simulacra of low-luminosity dwarf galaxies, we re-examine that simulation for planar structures, using the same method as employed by Ibata et al. (2013) on the real M31 satellites. We find that 0.04\% of host galaxies display satellite alignments that are at least as extreme as the observations, when we consider their extent, thickness and number of members rotating in the same sense. We further investigate the angular momentum properties of the co-planar satellites, and find that the median of the specific angular momentum derived from the line of sight velocities in the real M31 structure ($1.3\times10^4$ km/s kpc) is very high compared to systems drawn from the simulations. This analysis confirms that it is highly unlikely that the observed structure around the Andromeda galaxy is due to a chance occurrence. Interestingly, the few extreme systems that are similar to M31 arise from the accretion of a massive sub-halo with its own spatially-concentrated entourage of orphan satellites.

We undertake the first detailed analysis of the stellar population and spatial properties of a diffuse substructure in the outer halo of M31. The South-West Cloud lies at a projected distance of ~100 kpc from the centre of M31, and extends for at least ~50 kpc in projection. We use Pan-Andromeda Archaeological Survey photometry of red giant branch stars to determine a distance to the South-West Cloud of 793 +/- 45 kpc. The metallicity of the cloud is found to be [Fe/H] = -1.3 +/- 0.1. This is consistent with the coincident globular clusters PAndAS-7 and PAndAS-8, which have metallicities determined using an independent technique of [Fe/H] = -1.35 +/- 0.15. We measure a brightness for the Cloud of M_V = -12.1 mag; this is ~75 per cent of the luminosity implied by the luminosity-metallicity relation. Under the assumption that the South-West Cloud is the visible remnant of an accreted dwarf satellite, this suggests that the progenitor object was amongst M31's brightest dwarf galaxies prior to disruption.

We undertake an investigation into the spatial structure of the M31 satellite system utilizing the distance distributions presented in a previous publication. These distances make use of the unique combination of depth and spatial coverage of the Pan-Andromeda Archaeological Survey (PAndAS) to provide a large, homogeneous sample consisting of 27 of M31's satellites, as well as M31 itself. We find that the satellite distribution, when viewed as a whole, is no more planar than one would expect from a random distribution of equal size. A disk consisting of 15 of the satellites is however found to be highly significant, and strikingly thin, with a root-mean-square thickness of just $12.34^{+0.75}_{-0.43}$ kpc. This disk is oriented approximately edge on with respect to the Milky Way and almost perpendicular to the Milky Way disk. It is also roughly orthogonal to the disk like structure regularly reported for the Milky Way satellite system and in close alignment with M31's Giant Stellar Stream. A similar analysis of the asymmetry of the M31 satellite distribution finds that it is also significantly larger than one would expect from a random distribution. In particular, it is remarkable that 20 of the 27 satellites most likely lie on the Milky Way side of the galaxy, with the asymmetry being most pronounced within the satellite subset forming the aforementioned disk. This lopsidedness is all the more intriguing in light of the apparent orthogonality observed between the satellite disk structures of the Milky Way and M31.

Dwarf satellite galaxies are thought to be the remnants of the population of primordial structures that coalesced to form giant galaxies like the Milky Way. An early analysis noted that dwarf galaxies may not be isotropically distributed around our Galaxy, as several are correlated with streams of HI emission, and possibly form co-planar groups. These suspicions are supported by recent analyses, and it has been claimed that the apparently planar distribution of satellites is not predicted within standard cosmology, and cannot simply represent a memory of past coherent accretion. However, other studies dispute this conclusion. Here we report the existence (99.998% significance) of a planar sub-group of satellites in the Andromeda galaxy, comprising approximately 50% of the population. The structure is vast: at least 400 kpc in diameter, but also extremely thin, with a perpendicular scatter <14.1 kpc (99% confidence). Radial velocity measurements reveal that the satellites in this structure have the same sense of rotation about their host. This finding shows conclusively that substantial numbers of dwarf satellite galaxies share the same dynamical orbital properties and direction of angular momentum, a new insight for our understanding of the origin of these most dark matter dominated of galaxies. Intriguingly, the plane we identify is approximately aligned with the pole of the Milky Way's disk and is co-planar with the Milky Way to Andromeda position vector. The existence of such extensive coherent kinematic structures within the halos of massive galaxies is a fact that must be explained within the framework of galaxy formation and cosmology.

In `A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (PART I),' a new technique was introduced for obtaining distances using the TRGB standard candle. Here we describe a useful complement to the technique with the potential to further reduce the uncertainty in our distance measurements by incorporating a matched-filter weighting scheme into the model likelihood calculations. In this scheme, stars are weighted according to their probability of being true object members. We then re-test our modified algorithm using random-realization artificial data to verify the validity of the generated posterior probability distributions (PPDs) and proceed to apply the algorithm to the satellite system of M31, culminating in a 3D view of the system. Further to the distributions thus obtained, we apply a satellite-specific prior on the satellite distances to weight the resulting distance posterior distributions, based on the halo density profile. Thus in a single publication, using a single method, a comprehensive coverage of the distances to the companion galaxies of M31 is presented, encompassing the dwarf spheroidals Andromedas I - III, V, IX-XXVII and XXX along with NGC147, NGC 185, M33 and M31 itself. Of these, the distances to Andromeda XXIV - XXVII and Andromeda XXX have never before been derived using the TRGB. Object distances are determined from high-resolution tip magnitude posterior distributions generated using the Markov Chain Monte Carlo (MCMC) technique and associated sampling of these distributions to take into account uncertainties in foreground extinction and the absolute magnitude of the TRGB as well as photometric errors. The distance PPDs obtained for each object both with, and without the aforementioned prior are made available to the reader in tabular form...

We present a new approach for identifying the Tip of the Red Giant Branch (TRGB) which, as we show, works robustly even on sparsely populated targets. Moreover, the approach is highly adaptable to the available data for the stellar population under study, with prior information readily incorporable into the algorithm. The uncertainty in the derived distances is also made tangible and easily calculable from posterior probability distributions. We provide an outline of the development of the algorithm and present the results of tests designed to characterize its capabilities and limitations. We then apply the new algorithm to three M31 satellites: Andromeda I, Andromeda II and the fainter Andromeda XXIII, using data from the Pan-Andromeda Archaeological Survey (PAndAS), and derive their distances as $731^{(+ 5) + 18}_{(- 4) - 17}$ kpc, $634^{(+ 2) + 15}_{(- 2) - 14}$ kpc and $733^{(+ 13)+ 23}_{(- 11) - 22}$ kpc respectively, where the errors appearing in parentheses are the components intrinsic to the method, while the larger values give the errors after accounting for additional sources of error. These results agree well with the best distance determinations in the literature and provide the smallest uncertainties to date. This paper is an introduction to the workings and capabilities of our new approach in its basic form, while a follow-up paper shall make full use of the method's ability to incorporate priors and use the resulting algorithm to systematically obtain distances to all of M31's satellites identifiable in the PAndAS survey area.