We extend our previous work on the age-chemical abundance structure of the Galactic outer disc to the inner disc (4 < r < 8 kpc) based on the SDSS/APOGEE survey. Different from the outer disc, the inner disc stars exhibit a clear bimodal distribution in the [Mg/Fe]-[Fe/H] plane. While a number of scenarios have been proposed in the literature, it remains challenging to recover this bimodal distribution with theoretical models. To this end, we present a chemical evolution model embedding a complex multi-phase inner disc formation scenario that matches the observed bimodal [Mg/Fe]-[Fe/H] distribution. In this scenario, the formation of the inner disc is dominated by two main starburst episodes 6 Gyr apart with secular, low-level star formation activity in between. In our model, the first starburst occurs at early cosmic times (t~1 Gyr) and the second one 6 Gyr later at a cosmic time of t~7 Gyr. Both these starburst episodes are associated with gas accretion events in our model, and are quenched rapidly. The first starburst leads to the formation of the high-$\alpha$ sequence, and the second starburst leads to the formation of the metal-poor low-$\alpha$ sequence. The metal-rich low-$\alpha$ stars, instead, form during the secular evolution phase between the two bursts. Our model shows that the $\alpha$-dichotomy originates from the rapid suppression of star formation after the first starburst. The two starburst episodes are likely to be responsible for the formation of the geometric thick disc (z >1 kpc), with the old inner thick disc and the young outer thick disc forming during the first and the second starbursts, respectively.

We report the first APOGEE metallicities and alpha-element abundances measured for 3600 red giant stars spanning a large radial range of both the Large (LMC) and Small Magellanic Clouds (SMC), the largest Milky Way dwarf galaxies. Our sample is an order of magnitude larger than that of previous studies, and extends to much larger radial distances. These are the first results presented that make use of the newly installed Southern APOGEE instrument on the du Pont telescope at Las Campanas Observatory. Our unbiased sample of the LMC spans a large range in metallicity, from [Fe/H]=-0.2 to very metal-poor stars with [Fe/H]=-2.5, the most metal-poor Magellanic Clouds (MCs) stars detected to date. The LMC [alpha/Fe]-[Fe/H] distribution is very flat over a large metallicity range, but rises by ~0.1 dex at -1.0<[Fe/H]<-0.5. We interpret this as a sign of the known recent increase in MC star-formation activity, and are able to reproduce the pattern with a chemical evolution model that includes a recent "starburst". At the metal-poor end, we capture the increase of [alpha/Fe] with decreasing [Fe/H], and constrain the "alpha-knee" to [Fe/H]<-2.2 in both MCs, implying a low star-formation efficiency of ~0.01 Gyr^-1. The MC knees are more metal poor than those of less massive Milky Way (MW) dwarf galaxies such as Fornax, Sculptor, or Sagittarius. One possible interpretation is that the MCs formed in a lower-density environment than the MW, a hypothesis that is consistent with the paradigm that the MCs fell into the MW's gravitational potential only recently.

Using 3D positions and kinematics of stars relative to the Sagittarius (Sgr) orbital plane and angular momentum, we identify 166 Sgr stream members observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) that also have Gaia DR2 astrometry. This sample of 63/103 stars in the Sgr trailing/leading arm are combined with an APOGEE sample of 710 members of the Sgr dwarf spheroidal core (385 of them newly presented here) to establish differences of 0.6 dex in median metallicity and 0.1 dex in [$\alpha$/Fe] between our Sgr core and dynamically older stream samples. Mild chemical gradients are found internally along each arm, but these steepen when anchored by core stars. With a model of Sgr tidal disruption providing estimated dynamical ages (i.e., stripping times) for each stream star, we find a mean metallicity gradient of 0.12 +/- 0.03 dex/Gyr for stars stripped from Sgr over time. For the first time, an [$\alpha$/Fe] gradient is also measured within the stream, at 0.02 +/- 0.01 dex/Gyr using magnesium abundances and 0.04 +/- 0.01 dex/Gyr using silicon, which imply that the Sgr progenitor had significant radial abundance gradients. We discuss the magnitude of those inferred gradients and their implication for the nature of the Sgr progenitor within the context of the current family of Milky Way satellite galaxies, and suggest that more sophisticated Sgr models are needed to properly interpret the growing chemodynamical detail we have on the Sgr system.

Using the photometric data from the Next Generation Fornax Survey, we find a significant radial alignment signal among the Fornax dwarf galaxies. For the first time, we report that the radial alignment signal of nucleated dwarfs is stronger than that of non-nucleated ones at 2.4$\sigma$ confidence level, and the dwarfs located in the outer region ($R>R_{\rm{vir}}/3$; $R_{\rm{vir}}$ is the Fornax virial radius) show slightly stronger radial alignment signal than those in the inner region ($R<R_{\rm{vir}}/3$) at $1.5\sigma$ level. We also find that the significance of radial alignment signal is independent of the luminosities or sizes of the dwarfs.

Understanding the evolutionary history of the Magellanic Clouds requires an in-depth exploration and characterization of the stellar content in their outer regions, which ultimately are key to tracing the epochs and nature of past interactions. We present new deep images of a shell-like over-density of stars in the outskirts of the Small Magellanic Cloud (SMC). The shell, also detected in photographic plates dating back to the fifties, is located at ~1.9 degr from the center of the SMC in the north-east direction.The structure and stellar content of this feature were studied with multi-band, optical data from the Survey of the MAgellanic Stellar History (SMASH) carried out with the Dark Energy Camera on the Blanco Telescope at Cerro Tololo Inter-American Observatory. We also investigate the kinematic of the stars in the shell using the Gaia Data Release 2. The shell is composed of a young population with an age ~ 150 Myr, with no contribution from an old population. Thus, it is hard to explain its origin as the remnant of a tidally disrupted stellar system. The spatial distribution of the young main-sequence stars shows a rich sub-structure, with a spiral arm-like feature emanating from the main shell and a separated small arc of young stars close to the globular cluster NGC 362. We find that the absolute g-band magnitude of the shell is M_{g,shell} = -10.78+/- 0.02, with a surface brightness of mu_{g,shell} = 25.81+/- 0.01 mag/arcsec^{2}. We have not found any evidence that this feature is of tidal origin or a bright part of a spiral arm-like structure. Instead, we suggest that the shell formed in a recent star formation event, likely triggered by an interaction with the Large Magellanic Cloud and/or the Milky Way, ~150 Myr ago.

We present a study of the three-dimensional (3D) structure of the Large Magellanic Cloud (LMC) using ~2.2 million red clump (RC) stars selected from the Survey of the MAgellanic Stellar History. To correct for line-of-sight dust extinction, the intrinsic RC color and magnitude and their radial dependence are carefully measured by using internal nearly dust-free regions. These are then used to construct an accurate 2D reddening map (165 square degrees with ~10 arcmin resolution) of the LMC disk and the 3D spatial distribution of RC stars. An inclined disk model is fit to the 2D distance map yielding a best-fit inclination angle i = 25.86(+0.73,-1.39) degrees with random errors of +\-0.19 degrees and line-of-nodes position angle theta = 149.23(+6.43,-8.35) degrees with random errors of +/-0.49 degrees. These angles vary with galactic radius, indicating that the LMC disk is warped and twisted likely due to the repeated tidal interactions with the Small Magellanic Cloud (SMC). For the first time, our data reveal a significant warp in the southwestern part of the outer disk starting at rho ~ 7 degrees that departs from the defined LMC plane up to ~4 kpc toward the SMC, suggesting that it originated from a strong interaction with the SMC. In addition, the inner disk encompassing the off-centered bar appears to be tilted up to 5-15 degrees relative to the rest of the LMC disk. These findings on the outer warp and the tilted bar are consistent with the predictions from the Besla et al. simulation of a recent direct collision with the SMC.

We explore the stellar structure of the Large Magellanic Cloud (LMC) disk using data from the Survey of the MAgellanic Stellar History (SMASH) and the Dark Energy Survey. We detect a ring-like stellar overdensity in the red clump star count map at a radius of ~6 degrees (~5.2 kpc at the LMC distance) that is continuous over ~270 degrees in position angle and is only limited by the current data coverage. The overdensity shows an amplitude up to 2.5 times higher than that of the underlying smooth disk. This structure might be related to the multiple arms found by de Vaucouleurs. We find that the overdensity shows spatial correlation with intermediate-age star clusters, but not with young (< 1 Gyr) main-sequence stars, indicating the stellar populations associated with the overdensity are intermediate in age or older. Our findings on the LMC overdensity can be explained by either of two distinct formation mechanisms of a ring-like overdensity: (1) the overdensity formed out of an asymmetric one-armed spiral wrapping around the LMC main body, which is induced by repeated encounters with the Small Magellanic Cloud (SMC) over the last Gyr, or (2) the overdensity formed very recently as a tidal response to a direct collision with the SMC. Although the measured properties of the overdensity alone cannot distinguish between the two candidate scenarios, the consistency with both scenarios suggests that the ring-like overdensity is likely a product of tidal interaction with the SMC, but not with the Milky Way halo.

We describe a deep, systematic imaging study of satellites in the outer halo of the Milky Way. Our sample consists of 58 stellar overdensities --- i.e., substructures classified as either globular clusters, classical dwarf galaxies, or ultra-faint dwarf galaxies --- that are located at Galactocentric distances of R$_{\rm GC}$ > 25 kpc (outer halo) and out to ~400 kpc. This includes 44 objects for which we have acquired deep, wide-field, $g-$ and $r-$band imaging with the MegaCam mosaic cameras on the 3.6m Canada-France-Hawaii Telescope and the 6.5m Magellan-Clay telescope. These data are supplemented by archival imaging, or published $gr$ photometry, for an additional 14 objects, most of which were discovered recently in the Dark Energy Survey (DES). We describe the scientific motivation for our survey, including sample selection, observing strategy, data reduction pipeline, calibration procedures, and the depth and precision of the photometry. The typical 5$\sigma$ point-source limiting magnitudes for our MegaCam imaging --- which collectively covers an area of ~52 deg$^{2}$ --- are $g_{\rm lim}$ ~25.6 and $r_{\rm lim}$ ~25.3 AB mag. These limits are comparable to those from the coadded DES images and are roughly a half-magnitude deeper than will be reached in a single visit with LSST. Our photometric catalog thus provides the deepest and most uniform photometric database of Milky Way satellites available for the foreseeable future. In other papers in this series, we have used these data to explore the blue straggler populations in these objects, their density distributions, star formation histories, scaling relations and possible foreground structures.

We present structural parameters from a wide-field homogeneous imaging survey of Milky Way satellites carried out with the MegaCam imagers on the 3.6m Canada-France-Hawaii Telescope (CFHT) and 6.5m Magellan-Clay telescope. Our survey targets an unbiased sample of "outer halo" satellites (i.e., substructures having Galactocentric distances greater than 25 kpc) and includes classical dSph galaxies, ultra-faint dwarfs, and remote globular clusters. We combine deep, panoramic $gr$ imaging for 44 satellites and archival $gr$ imaging for 14 additional objects (primarily obtained with the DECam instrument as part of the Dark Energy Survey), to measure photometric and structural parameters for 58 outer halo satellites. This is the largest and most uniform analysis of Milky Way satellites undertaken to date and represents roughly three quarters ($58/81 \simeq$72\%) of all known outer halo satellites. We use a maximum-likelihood method to fit four density laws to each object in our survey: exponential, Plummer, King and Sersic models. We examine systematically the isodensity contour maps and color magnitude diagrams for each of our program objects, present a comparison with previous results, and tabulate our best-fit photometric and structural parameters, including ellipticities, position angles, effective radii, Sersic indices, absolute magnitudes, and surface brightness measurements. We investigate the distribution of outer halo satellites in the size-magnitude diagram, and show that the current sample of outer halo substructures spans a wide range in effective radius, luminosity and surface brightness, with little evidence for a clean separation into star cluster and galaxy populations at the faintest luminosities and surface brightnesses.

We present the detection of very extended stellar populations around the Large Magellanic Cloud (LMC) out to R~21 degrees, or ~18.5 kpc at the LMC distance of 50 kpc, as detected in the Survey of the MAgellanic Stellar History (SMASH) performed with the Dark Energy Camera on the NOAO Blanco 4m Telescope. The deep (g~24) SMASH color magnitude diagrams (CMDs) clearly reveal old (~9 Gyr), metal-poor ([Fe/H]=-0.8 dex) main-sequence stars at a distance of 50 kpc. The surface brightness of these detections is extremely low with our most distant detection having 34 mag per arcsec squared in g-band. The SMASH radial density profile breaks from the inner LMC exponential decline at ~13-15 degrees and a second component at larger radii has a shallower slope with power-law index of -2.2 that contributes ~0.4% of the LMC's total stellar mass. In addition, the SMASH densities exhibit large scatter around our best-fit model of ~70% indicating that the envelope of stellar material in the LMC periphery is highly disturbed. We also use data from the NOAO Source catalog to map the LMC main-sequence populations at intermediate radii and detect a steep dropoff in density on the eastern side of the LMC (at R~8 deg) as well as an extended structure to the far northeast. These combined results confirm the existence of a very extended, low-density envelope of stellar material with disturbed shape around the LMC. The exact origin of this structure remains unclear but the leading options include a classical accreted halo or tidally stripped outer disk material.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present Keck/DEIMOS spectroscopy of globular clusters (GCs) around the ultra-diffuse galaxies (UDGs) VLSB-B, VLSB-D, and VCC615 located in the central regions of the Virgo cluster. We spectroscopically identify 4, 12, and 7 GC satellites of these UDGs, respectively. We find that the three UDGs have systemic velocities ($V_{sys}$) consistent with being in the Virgo cluster, and that they span a wide range of velocity dispersions, from $\sim 16$ to $\sim 47$ km/s, and high dynamical mass-to-light ratios within the radius that contains half the number of GCs ($ 407^{+916}_{-407}$, $21^{+15}_{-11}$, $60^{+65}_{-38}$, respectively). VLSB-D shows possible evidence for rotation along the stellar major axis and its $V_{sys}$ is consistent with that of the massive galaxy M84 and the center of the Virgo cluster itself. These findings, in addition to having a dynamically and spatially ($\sim 1$ kpc) off-centered nucleus and being extremely elongated, suggest that VLSB-D could be tidally perturbed. On the contrary, VLSB-B and VCC615 show no signals of tidal deformation. Whereas the dynamics of VLSB-D suggest that it has a less massive dark matter halo than expected for its stellar mass, VLSB-B and VCC615 are consistent with a $\sim 10^{12}$ M$_{\odot}$ dark matter halo. Although our samples of galaxies and GCs are small, these results suggest that UDGs may be a diverse population, with their low surface brightnesses being the result of very early formation, tidal disruption, or a combination of the two.

We derive stellar population parameters for a representative sample of ultracompact dwarfs (UCDs) and a large sample of massive globular clusters (GCs) with stellar masses $\gtrsim$ 10$^{6}$ $M_{\odot}$ in the central galaxy M87 of the Virgo galaxy cluster, based on model fitting to the Lick-index measurements from both the literature and new observations. After necessary spectral stacking of the relatively faint objects in our initial sample of 40 UCDs and 118 GCs, we obtain 30 sets of Lick-index measurements for UCDs and 80 for GCs. The M87 UCDs have ages $\gtrsim$ 8 Gyr and [$\alpha$/Fe] $\simeq$ 0.4 dex, in agreement with previous studies based on smaller samples. The literature UCDs, located in lower-density environments than M87, extend to younger ages and smaller [$\alpha$/Fe] (at given metallicities) than M87 UCDs, resembling the environmental dependence of the Virgo dE nuclei. The UCDs exhibit a positive mass-metallicity relation (MZR), which flattens and connects compact ellipticals at stellar masses $\gtrsim$ 10$^{8}$ $M_{\odot}$. The Virgo dE nuclei largely follow the average MZR of UCDs, whereas most of the M87 GCs are offset towards higher metallicities for given stellar masses. The difference between the mass-metallicity distributions of UCDs and GCs may be qualitatively understood as a result of their different physical sizes at birth in a self-enrichment scenario or of galactic nuclear cluster star formation efficiency being relatively low in a tidal stripping scenario for UCD formation. The existing observations provide the necessary but not sufficient evidence for tidally stripped dE nuclei being the dominant contributors to the M87 UCDs.

We present the results of a wide spectroscopic survey aimed at detecting extragalactic globular clusters (GCs) in the core of the Fornax cluster. About 4500 low resolution spectra (from 4800 to 10000 \AA ) were observed in 25 VLT/VIMOS masks covering the central 1 deg$^{2}$ around the dominant galaxy NGC 1399 corresponding to $\sim$175 kpc galactocentric radius. We describe the methodology used for data reduction and data analysis. We found a total of 387 unique physical objects (372 GCs and 15 ultra compact dwarfs) in the field covered by our observations. Most of these objects are associated with NGC 1399, with only 10% likely belonging to other giant galaxies. The new VIMOS dataset is complementary to the many GC catalogues already present in the literature and it brings the total number of tracer particles around NGC 1399 to more than 1130 objects. With this comprehensive radial velocity sample we have found that the velocity dispersion of the GC population (equally for red and blue GC populations) shows a relatively sharp increase from low velocity dispersion ($\sim250$-$350$ kms$^{-1}$) to high velocity dispersion ($\sim300$-$400$ kms$^{-1}$) at projected radius of $\approx10$ arcmin ($\sim60$ kpc) from the galaxy centre. This suggests that at a projected radius of $\approx60$ kpc both blue and red GC populations begin to be governed by the dominating Fornax cluster potential, rather than by the central NGC 1399 galaxy potential. This kinematic evidence corroborates similar results found using surface brightness analysis and planetary nebulae kinematics.

We report the results of a systematic photometric survey of the peripheral regions of a sample of fourteen globular clusters in the outer halo of the Milky Way at distances d_GC>25 kpc from the Galactic centre. The survey is aimed at searching for the remnants of the host satellite galaxies where these clusters could originally have been formed before being accreted onto the Galactic halo. The limiting surface brightness varies within our sample, but reaches muV_lim=30-32 mag arcsec^-2. For only two globular clusters (NGC 7492 and Whiting 1; already suggested to be associated with the Sagittarius galaxy) we detect extended stellar populations that cannot be associated with either the clusters themselves or with the surrounding Galactic field population. We show that the lack of substructures around globular clusters at these Galactocentric distances is still compatible with the predictions of cosmological simulations whereby in the outer halo the Galactic globular cluster system is built up through hierarchical accretion at early epochs.

We present deep, wide-field Subaru Hyper Suprime-Cam photometry of two recently discovered satellites of the Milky Way (MW): Columba I and Triangulum II. The color magnitude diagrams of both objects point to exclusively old and metal-poor stellar populations. We re-derive structural parameters and luminosities of these satellites, and find $M_{\rm V, Col~I} = -4.2\pm0.2$ for Col I and $M_{\rm V, Tri~II} = -1.2\pm0.4$ for Tri II, with corresponding half-light radii of $r_{\rm h, Col~I} = 117\pm17$ pc and $r_{\rm h, Tri~II} = 21\pm4$ pc. The properties of both systems are consistent with observed scaling relations for MW dwarf galaxies. Based on archival data, we derive upper limits on the neutral gas content of these dwarfs, and find that they lack HI, as do the majority of observed satellites within the MW virial radius. Neither satellite shows evidence of tidal stripping in the form of extensions or distortions in matched-filter stellar density maps or surface density profiles. However, the smaller Tri II system is relatively metal-rich for its luminosity (compared to other MW satellites), possibly because it has been tidally stripped. Through a suite of orbit simulations, we show that Tri II is approaching pericenter of its eccentric orbit, a stage at which tidal debris is unlikely to be seen. In addition, we find that Tri II may be on its first infall into the MW, which helps explain its unique properties among MW dwarfs. Further evidence that Tri II is likely an ultra-faint dwarf comes from its stellar mass function, which is similar to those of other MW dwarfs.

APOGEE-2 is a high-resolution, near-infrared spectroscopic survey observing roughly 300,000 stars across the entire sky. It is the successor to APOGEE and is part of the Sloan Digital Sky Survey IV (SDSS-IV). APOGEE-2 is expanding upon APOGEE's goals of addressing critical questions of stellar astrophysics, stellar populations, and Galactic chemodynamical evolution using (1) an enhanced set of target types and (2) a second spectrograph at Las Campanas Observatory in Chile. APOGEE-2 is targeting red giant branch (RGB) and red clump (RC) stars, RR Lyrae, low-mass dwarf stars, young stellar objects, and numerous other Milky Way and Local Group sources across the entire sky from both hemispheres. In this paper, we describe the APOGEE-2 observational design, target selection catalogs and algorithms, and the targeting-related documentation included in the SDSS data releases.

We present the survey strategy and early results of the "Satellites Around Galactic Analogs" (SAGA) Survey. The SAGA Survey's goal is to measure the distribution of satellite galaxies around 100 systems analogous to the Milky Way down to the luminosity of the Leo I dwarf galaxy ($ M_r < -12.3 $). We define a Milky Way analog based on $K$-band luminosity and local environment. Here, we present satellite luminosity functions for 8 Milky Way analog galaxies between 20 to 40 Mpc. These systems have nearly complete spectroscopic coverage of candidate satellites within the projected host virial radius down to $ r_o < 20.75 $ using low redshift $gri$ color criteria. We have discovered a total of 25 new satellite galaxies: 14 new satellite galaxies meet our formal criteria around our complete host systems, plus 11 additional satellites in either incompletely surveyed hosts or below our formal magnitude limit. Combined with 13 previously known satellites, there are a total of 27 satellites around 8 complete Milky Way analog hosts. We find a wide distribution in the number of satellites per host, from 1 to 9, in the luminosity range for which there are five Milky Way satellites. Standard abundance matching extrapolated from higher luminosities predicts less scatter between hosts and a steeper luminosity function slope than observed. We find that the majority of satellites (26 of 27) are star-forming. These early results indicate that the Milky Way has a different satellite population than typical in our sample, potentially changing the physical interpretation of measurements based only on the Milky Way's satellite galaxies.

The Large and Small Magellanic Clouds (LMC and SMC) are unique local laboratories for studying the formation and evolution of small galaxies in exquisite detail. The Survey of the MAgellanic Stellar History (SMASH) is an NOAO community DECam survey of the Clouds mapping 480 square degrees (distributed over ~2400 square degrees at ~20% filling factor) to ~24th mag in ugriz with the goal of identifying broadly distributed, low surface brightness stellar populations associated with the stellar halos and tidal debris of the Clouds. SMASH will also derive spatially-resolved star formation histories covering all ages out to large radii from the MCs that will further complement our understanding of their formation. Here, we present a summary of the survey, its data reduction, and a description of the first public Data Release (DR1). The SMASH DECam data have been reduced with a combination of the NOAO Community Pipeline, PHOTRED, an automated PSF photometry pipeline based mainly on the DAOPHOT suite, and custom calibration software. The attained astrometric precision is ~15 mas and the accuracy is ~2 mas with respect to the Gaia DR1 astrometric reference frame. The photometric precision is ~0.5-0.7% in griz and ~1% in u with a calibration accuracy of ~1.3% in all bands. The median 5 sigma point source depths in ugriz bands are 23.9, 24.8, 24.5, 24.2, 23.5 mag. The SMASH data already have been used to discover the Hydra II Milky Way satellite, the SMASH 1 old globular cluster likely associated with the LMC, and very extended stellar populations around the LMC out to R~18.4 kpc. SMASH DR1 contains measurements of ~100 million objects distributed in 61 fields. A prototype version of the NOAO Data Lab provides data access, including a data discovery tool, SMASH database access, an image cutout service, and a Jupyter notebook server with example notebooks for exploratory analysis.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We build a background cluster candidate catalog from the Next Generation Virgo Cluster Survey, using our detection algorithm RedGOLD. The NGVS covers 104$deg^2$ of the Virgo cluster in the $u^*,g,r,i,z$-bandpasses to a depth of $ g \sim 25.7$~mag (5$\sigma$). Part of the survey was not covered or has shallow observations in the $r$--band. We build two cluster catalogs: one using all bandpasses, for the fields with deep $r$--band observations ($\sim 20 \ deg^2$), and the other using four bandpasses ($u^*,g,i,z$) for the entire NGVS area. Based on our previous CFHT-LS W1 studies, we estimate that both of our catalogs are $\sim100\%$($\sim70\%$) complete and $\sim80\%$ pure, at $z\le 0.6$($z\lesssim1$), for galaxy clusters with masses of $M\gtrsim10^{14}\ M_{\odot}$. We show that when using four bandpasses, though the photometric redshift accuracy is lower, RedGOLD detects massive galaxy clusters up to $z\sim 1$ with completeness and purity similar to the five-band case. This is achieved when taking into account the bias in the richness estimation, which is $\sim40\%$ lower at $0.5\le z<0.6$ and $\sim20\%$ higher at $0.6<z< 0.8$, with respect to the five-band case. RedGOLD recovers all the X-ray clusters in the area with mass $M_{500} > 1.4 \times 10^{14} \rm M_{\odot}$ and $0.08<z<0.5$. Because of our different cluster richness limits and the NGVS depth, our catalogs reach to lower masses than the published redMaPPer cluster catalog over the area, and we recover $\sim 90-100\%$ of its detections.

We present a study of the kinematics of 101 stars observed with VIMOS around Whiting1, a globular cluster embedded in the Sagittarius tidal stream. The obtained velocity distribution shows the presence of two wraps of that halo substructure at the same heliocentric distance as that of the cluster and with well differentiated mean radial velocities. The most prominent velocity component seems to be associated with the trailing arm of Sagittarius with $<v_{\rm r}>$ ~ -130 km/s, which is consistent with the velocity of Whiting1. This result supports that this globular cluster was formed in Sagittarius and recently accreted by the Milky Way. The second component with $<v_{\rm r}>$ ~ 120 km/s might correspond to the leading arm of Sagittarius, which has been predicted by numerical simulations but with no conclusive observational evidence of its existence presented so far. This detection of the old leading wrap of Sagittarius in the southern hemisphere may be used to confirm and further constrain the models for its orbit and evolution.

We report chemical abundances obtained by SDSS-III/APOGEE for giant stars in five globular clusters located within 2.2 kpc of the Galactic centre. We detect the presence of multiple stellar populations in four of those clusters (NGC 6553, NGC 6528, Terzan 5, and Palomar 6) and find strong evidence for their presence in NGC 6522. All clusters present a significant spread in the abundances of N, C, Na, and Al, with the usual correlations and anti-correlations between various abundances seen in other globular clusters. Our results provide important quantitative constraints on theoretical models for self-enrichment of globular clusters, by testing their predictions for the dependence of yields of elements such as Na, N, C, and Al on metallicity. They also confirm that, under the assumption that field N-rich stars originate from globular cluster destruction, they can be used as tracers of their parental systems in the high- metallicity regime.

We use deep optical photometry from the Next Generation Virgo Cluster Survey [NGVS] to investigate the color-magnitude diagram for the galaxies inhabiting the core of this cluster. The sensitivity of the NGVS imaging allows us to continuously probe galaxy colors over a factor of $\sim 2 \times 10^5$ in luminosity, from brightest cluster galaxies to scales overlapping classical satellites of the Milky Way [$M_{g^{\prime}}$ $\sim$ $-$9; $M_{*}$ $\sim 10^6$ M$_{\odot}$], within a single environment. Remarkably, we find the first evidence that the RS flattens in all colors at the faint-magnitude end [starting between $-$14 $\le$ $M_{g^{\prime}}$ $\le$ $-$13, around $M_{*}$ $\sim 4 \times 10^7$ M$_{\odot}$], with the slope decreasing to $\sim$60% or less of its value at brighter magnitudes. This could indicate that the stellar populations of faint dwarfs in Virgo's core share similar characteristics [e.g. constant mean age] over $\sim$3 mags in luminosity, suggesting that these galaxies were quenched coevally, likely via pre-processing in smaller hosts. We also compare our results to galaxy formation models, finding that the RS in model clusters have slopes at intermediate magnitudes that are too shallow, and in the case of semi-analytic models, do not reproduce the flattening seen at both extremes [bright/faint] of the Virgo RS. Deficiencies in the chemical evolution of model galaxies likely contribute to the model-data discrepancies at all masses, while overly efficient quenching may also be a factor at dwarf scales. Deep UV and near-IR photometry are required to unambiguously diagnose the cause of the faint-end flattening.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Deep observations of the dwarf elliptical (dE) galaxy NGC 1396 (M$_V = -16.60$, Mass $\sim 4\times10^8$ M$_\odot$), located in the Fornax cluster, have been performed with the VLT/ MUSE spectrograph in the wavelength region from $4750-9350$ \AA{}. In this paper we present a stellar population analysis studying chemical abundances, the star formation history (SFH) and the stellar initial mass function (IMF) as a function of galacto-centric distance. Different, independent ways to analyse the stellar populations result in a luminosity-weighted age of $\sim$ 6 Gyr and a metallicity [Fe/H]$\sim$ $-0.4$, similar to other dEs of similar mass. We find unusually overabundant values of [Ca/Fe] $\sim +0.1$, and under-abundant Sodium, with [Na/Fe] values around $-0.1$, while [Mg/Fe] is overabundant at all radii, increasing from $\sim+0.1$ in the centre to $\sim +0.2$ dex. We notice a significant metallicity and age gradient within this dwarf galaxy. To constrain the stellar IMF of NGC 1396, we find that the IMF of NGC 1396 is consistent with either a Kroupa-like or a top-heavy distribution, while a bottom-heavy IMF is firmly ruled out. An analysis of the abundance ratios, and a comparison with galaxies in the Local Group, shows that the chemical enrichment history of NGC 1396 is similar to the Galactic disc, with an extended star formation history. This would be the case if the galaxy originated from a LMC-sized dwarf galaxy progenitor, which would lose its gas while falling into the Fornax cluster.

We report the discovery of 158 previously undetected dwarf galaxies in the Fornax cluster central regions using a deep coadded $u, g$ and $i$-band image obtained with the DECam wide-field camera mounted on the 4-meter Blanco telescope at the Cerro Tololo Interamerican Observatory as part of the {\it Next Generation Fornax Survey} (NGFS). The new dwarf galaxies have quasi-exponential light profiles, effective radii $0.1\!<\!r_e\!<\!2.8$ kpc and average effective surface brightness values $22.0\!<\!\mu_i\!<\!28.0$ mag arcsec$^{-2}$. We confirm the existence of ultra-diffuse galaxies (UDGs) in the Fornax core regions that resemble counterparts recently discovered in the Virgo and Coma galaxy clusters.~We also find extremely low surface brightness NGFS dwarfs, which are several magnitudes fainter than the classical UDGs. The faintest dwarf candidate in our NGFS sample has an absolute magnitude of $M_i\!=\!-8.0$\,mag. The nucleation fraction of the NGFS dwarf galaxy sample appears to decrease as a function of their total luminosity, reaching from a nucleation fraction of $>\!75\%$ at luminosities brighter than $M_i\!\simeq\!-15.0$ mag to $0\%$ at luminosities fainter than $M_i\!\simeq\!-10.0$ mag. The two-point correlation function analysis of the NGFS dwarf sample shows an excess on length scales below $\sim\!100$ kpc, pointing to the clustering of dwarf galaxies in the Fornax cluster core.

We use imaging from the Next Generation Virgo cluster Survey (NGVS) to present a comparative study of ultra-compact dwarf (UCD) galaxies associated with three prominent Virgo sub-clusters: those centered on the massive, red-sequence galaxies M87, M49 and M60. We show how UCDs can be selected with high completeness using a combination of half-light radius and location in color-color diagrams ($u^*iK_s$ or $u^*gz$). Although the central galaxies in each of these sub-clusters have nearly identical luminosities and stellar masses, we find large differences in the sizes of their UCD populations, with M87 containing ~3.5 and 7.8 times more UCDs than M49 and M60, respectively. The relative abundance of UCDs in the three regions scales in proportion to sub-cluster mass, as traced by X-ray gas mass, total gravitating mass, number of globular clusters, and number of nearby galaxies. We find that the UCDs are predominantly blue in color, with ~85% of the UCDs having colors similar to blue GCs and stellar nuclei of dwarf galaxies. We present evidence that UCDs surrounding M87 and M49 may follow a morphological sequence ordered by the prominence of their outer, low surface brightness envelope, ultimately merging with the sequence of nucleated low-mass galaxies, and that envelope prominence correlates with distance from either galaxy. Our analysis provides evidence that tidal stripping of nucleated galaxies is an important process in the formation of UCDs.

The spiral structure of the Milky Way is highly uncertain and is the subject of much discussion nowadays. We present the first result from a program that determines the properties of the Local spiral arm (LOA), together with a full description of the program. In this context we have made a comprehensive study of the young LOA open cluster NGC 2302, which includes a UBVRI photometric analysis and determination of its kinematic properties - proper motion and radial velocity - and of its orbital parameters. We determined the mean PM of NGC 2302 relative to the local field of disk stars, and, through a comparison with the UCAC4 catalog, we transformed this relative PM into an absolute one. Using medium-resolution spectroscopy of 26 stars in the field of NGC 2302, we derived its mean RV. Isochrone fits to the photometric diagrams allowed us to determine the fundamental parameters of NGC 2302, including reddening, distance, and age. The kinematic data and derived distance allowed us to determine the space motion of NGC 2302. This was done by adopting a time-independent, axisymmetric, and fully analytic gravitational potential for the MW. We obtained an absolute PM for NGC 2302 of ($\mu_{\alpha} \cos\delta,\mu_{\delta}) = (-2.09,-2.11)$ mas/yr, with standard errors of 0.410 and 0.400 mas/yr. The mean RV of NGC 2302 turned out to be 31.2 km/sec with a standard error of 0.7 km/sec. Isochrone fits displaced for this reddening and for a distance modulus of (m-M)o = 10.69 indicate an age of log(t) = 7.90-8.00 with a slight tendency toward the younger age. Inspection of the shape of the orbit of NGC 2302 and the resulting orbital parameters indicate that it is a typical population I object.

The origin of the Galactic halo stellar structure known as the Monoceros ring is still under debate. In this work, we study that halo substructure using deep CFHT wide-field photometry obtained for the globular clusters NGC2419 and Koposov2, where the presence of Monoceros becomes significant because of their coincident projected position. Using Sloan Digital Sky Survey photometry and spectroscopy in the area surrounding these globulars and beyond, where the same Monoceros population is detected, we conclude that a second feature, not likely to be associated with Milky Way disk stars along the line-of-sight, is present as foreground population. Our analysis suggests that the Monoceros ring might be composed of an old stellar population of age t ~ 9Gyr and a new component ~ 4Gyr younger at the same heliocentric distance. Alternatively, this detection might be associated with a second wrap of Monoceros in that direction of the sky and also indicate a metallicity spread in the ring. The detection of such a low-density feature in other sections of this halo substructure will shed light on its nature.

The origin of ultra-compact dwarfs (UCDs)--objects larger and more massive than typical globular clusters (GCs), but more compact than typical dwarf galaxies--has been hotly debated in the 15 years since their discovery. Even whether UCDs should be considered galactic in origin, or simply the most extreme GCs, is not yet settled. We present the dynamical properties of 97 spectroscopically confirmed UCDs (rh >~10 pc) and 911 GCs associated with central cD galaxy of the Virgo cluster, M87. Our UCDs, of which 89% have M_star > ~2X10^6 M_sun and 92% are as blue as the classic blue GCs, nearly triple the sample of previous confirmed Virgo UCDs, providing by far the best opportunity for studying the global dynamics of a UCD system. We found that (1) UCDs have a surface number density profile that is shallower than that of the blue GCs in the inner ~ 70 kpc and as steep as that of the red GCs at larger radii; (2) UCDs exhibit a significantly stronger rotation than the GCs, and the blue GCs seem to have a velocity field that is more consistent with that of the surrounding dwarf ellipticals than with that of UCDs; (3) UCDs have a radially increasing orbital anisotropy profile, and are tangentially-biased at radii < ~ 40 kpc and radially-biased further out. In contrast, the blue GCs become more tangentially-biased at larger radii beyond ~ 40 kpc; (4) GCs with M_star > 2X10^6 M_sun have rotational properties indistinguishable from the less massive ones, suggesting that it is the size, instead of mass, that differentiates UCDs from GCs as kinematically distinct populations. We conclude that most UCDs in M87 are not consistent with being merely the most luminous and extended examples of otherwise normal GCs. The radially-biased orbital structure of UCDs at large radii is in general agreement with the "tidally threshed dwarf galaxy" scenario.

We present neutral hydrogen (HI) imaging observations with the Westerbork Synthesis Radio Telescope of AGC198606, an HI cloud discovered in the ALFALFA 21cm survey. This object is of particular note as it is located 16 km/s and 1.2 degrees from the gas-bearing ultra-faint dwarf galaxy Leo T while having a similar HI linewidth and approximately twice the flux density. The HI imaging observations reveal a smooth, undisturbed HI morphology with a full extent of 23'x16' at the 5x10^18 atoms cm^-2 level. The velocity field of AGC198606 shows ordered motion with a gradient of ~25 km/s across ~20'. The global velocity dispersion is 9.3 km/s with no evidence for a narrow spectral component. No optical counterpart to AGC198606 is detected. The distance to AGC198606 is unknown, and we consider several different scenarios: physical association with Leo T, a minihalo at a distance of ~150 kpc based on the models of Faerman et al. (2013), and a cloud in the Galactic halo. At a distance of 420 kpc, AGC198606 would have an HI mass of 6.2x10^5 Msun, an HI radius of 1.4 kpc, and a dynamical mass within the HI extent of 1.5x10^8 Msun.

We present new constraints on the star formation histories of six ultra-faint dwarf galaxies: Bootes I, Canes Venatici II, Coma Berenices, Hercules, Leo IV, and Ursa Major I. Our analysis employs a combination of high-precision photometry obtained with the Advanced Camera for Surveys on the Hubble Space Telescope, medium-resolution spectroscopy obtained with the DEep Imaging Multi-Object Spectrograph on the W.M. Keck Observatory, and updated Victoria-Regina isochrones tailored to the abundance patterns appropriate for these galaxies. The data for five of these Milky Way satellites are best fit by a star formation history where at least 75% of the stars formed by z~10 (13.3 Gyr ago). All of the galaxies are consistent with 80% of the stars forming by z~6 (12.8 Gyr ago) and 100% of the stars forming by z~3 (11.6 Gyr ago). The similarly ancient populations of these galaxies support the hypothesis that star formation in the smallest dark matter sub-halos was suppressed by a global outside influence, such as the reionization of the universe.

In the hierarchical formation scenario in which the outer halo of the Milky Way is the result of the continuous accretion of low-mass galaxies, a fraction of the Galactic globular cluster system might have originated in and been accreted with already extinct dwarf galaxies. In this context, we expect that the remnants of these progenitor galaxies might be still populating the surroundings of those accreted globulars. In this work, we present wide-field photometry of a sample of 23 globular clusters in the Galactocentric distance range 10 < Rg < 40kpc, which we use to search for remnants of their hypothetical progenitor systems. Our deep photometry reveals the presence of underlying stellar populations along the line-of-sight of about half of the globulars included in our sample. Among the detections lying in the footprint of the Sagittarius tidal stream, which we identify via the comparison with its orbit derived from numerical simulations, only Whiting1 and NGC7492 seem to be inmersed in that remnant at a compatible heliocentric distance. We also confirm the existence of a subjacent Main-Sequence feature in the surroundings of NGC1851. A tentative detection of the vast Hercules-Aquila cloud is unveiled in the background of NGC7006.

We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalogue of 922 globular cluster line-of- sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. 263 globular clusters, mainly located beyond 40 kpc, are newly observed by the Next Generation Virgo Survey (NGVS). For the M2M modelling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modelled as a logarithmic potential. Our best dynamical model returns a stellar mass-to-light ratio in the I band of M/LI = 6.0(+ -0.3) M_sun/L_sun with a dark matter potential scale velocity of 591(+ -50) km/s and scale radius of 42(+ -10) kpc. We determine the total mass of M87 within 180 kpc to be (1.5 + - 0.2) 10^13 M_sun. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond 45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derived is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy beta parameter for the globular clusters in M87 is small, varying from -0.2 at the centre to 0.2 at 40 kpc, and gradually decreasing to zero at 120 kpc.

Deep photometry of the Small Magellanic Cloud (SMC) stellar periphery (R=4 deg, 4.2 kpc) is used to study its line-of-sight depth with red clump (RC) stars. The RC luminosity function is affected little by young (<1 Gyr) blue-loop stars in these regions because their main-sequence counterparts are not observed in the color magnitude diagrams. The SMC's eastern side is found to have a large line-of-sight depth (~23 kpc) while the western side has a much shallower depth (~10 kpc), consistent with previous photographic plate photometry results. We use a model SMC RC luminosity function to deconvolve the observed RC magnitudes and construct the density function in distance for our fields. Three of the eastern fields show a distance bimodality with one component at the "systemic" ~67 kpc SMC distance and a second component at ~55 kpc. Our data are not reproduced well by the various extant Magellanic Cloud and Stream simulations. However, the models predict that the known HI Magellanic Bridge (stretching from the SMC eastward towards the LMC) has a decreasing distance with angle from the SMC and should be seen in both the gaseous and stellar components. From comparison with these models we conclude that the most likely explanation for our newly identified ~55 kpc stellar structure in the eastern SMC is a stellar counterpart of the HI Magellanic Bridge that was tidally stripped from the SMC ~200 Myr ago during a close encounter with the LMC. This discovery has important implications for microlensing surveys of the SMC.

We present early results from a Hubble Space Telescope survey of the ultra-faint dwarf galaxies. These Milky Way satellites were discovered in the Sloan Digital Sky Survey, and appear to be an extension of the classical dwarf spheroidals to low luminosities, offering a new front in the efforts to understand the missing satellite problem. Because they are the least luminous, most dark matter dominated, and least chemically evolved galaxies known, the ultra-faint dwarfs are the best candidate fossils from the early universe. The primary goal of the survey is to measure the star-formation histories of these galaxies and discern any synchronization due to the reionization of the universe. We find that the six galaxies of our survey have very similar star-formation histories, and that each is dominated by stars older than 12 Gyr.

We present Keck/DEIMOS spectroscopy and CFHT/MegaCam photometry for the Milky Way globular cluster Palomar 13. We triple the number of spectroscopically confirmed members, including many repeat velocity measurements. Palomar 13 is the only known globular cluster with possible evidence for dark matter, based on a Keck/HIRES 21 star velocity dispersion of sigma=2.2+/-0.4 km/s. We reproduce this measurement, but demonstrate that it is inflated by unresolved binary stars. For our sample of 61 stars, the velocity dispersion is sigma=0.7(+0.6/-0.5) km/s. Combining our DEIMOS data with literature values, our final velocity dispersion is sigma=0.4(+0.4/-0.3) km/s. We determine a spectroscopic metallicity of [Fe/H]=-1.6+/-0.1 dex, placing a 1-sigma upper limit of sigma_[Fe/H]~0.2 dex on any internal metallicity spread. We determine Palomar 13's total luminosity to be M_V=-2.8+/-0.4, making it among the least luminous known globular clusters. The photometric isophotes are regular out to the half-light radius and mildly irregular outside this radius. The outer surface brightness profile slope is shallower than typical globular clusters (eta=-2.8+/-0.3). Thus at large radius, tidal debris is likely affecting the appearance of Palomar 13. Combining our luminosity with the intrinsic velocity dispersion, we find a dynamical mass of of M_1/2=1.3(+2.7/-1.3)x10^3 M_sun and a mass-to-light ratio of M/L_V=2.4(+5.0/-2.4) M_sun/L_sun. Within our measurement errors, the mass-to-light ratio agrees with the theoretical predictions for a single stellar population. We conclude that, while there is some evidence for tidal stripping at large radius, the dynamical mass of Palomar 13 is consistent with its stellar mass and neither significant dark matter, nor extreme tidal heating, is required to explain the cluster dynamics.

We present a homogeneous study of blue straggler stars across ten outer halo globular clusters, three classical dwarf spheroidal and nine ultra-faint galaxies based on deep and wide-field photometric data taken with MegaCam on the Canada-France-Hawaii Telescope. We find blue straggler stars to be ubiquitous among these Milky Way satellites. Based on these data, we can test the importance of primordial binaries or multiple systems on blue straggler star formation in low density environments. For the outer halo globular clusters we find an anti-correlation between the specific frequency of blue straggler and absolute magnitude, similar to that previously observed for inner halo clusters. When plotted against density and encounter rate, the frequency of blue stragglers are well fitted by single trends with smooth transitions between dwarf galaxies and globular clusters, which points to a common origin for their blue stragglers. The fraction of blue stragglers stays constant and high in the low encounter rate regime spanned by our dwarf galaxies, and decreases with density and encounter rate in the range spanned by our globular clusters. We find that young stars can mimic blue stragglers in dwarf galaxies only if their ages are 2.5+/-0.5 Gyr and they represent ~1-7% of the total number of stars, which we deem highly unlikely. These results point to mass-transfer or mergers of primordial binaries or multiple systems as the dominant blue straggler formation mechanism in low density systems.

We present constraints on the stellar initial mass function (IMF) in two ultra-faint dwarf (UFD) galaxies, Hercules and Leo IV, based on deep HST/ACS imaging. The Hercules and Leo IV galaxies are extremely low luminosity (M_V = -6.2, -5.5), metal-poor (<[Fe/H]>= -2.4, -2.5) systems that have old stellar populations (> 11 Gyr). Because they have long relaxation times, we can directly measure the low-mass stellar IMF by counting stars below the main-sequence turnoff without correcting for dynamical evolution. Over the stellar mass range probed by our data, 0.52 - 0.77 Msun, the IMF is best fit by a power-law slope of alpha = 1.2^{+0.4}_{-0.5} for Hercules and alpha = 1.3 +/- 0.8 for Leo IV. For Hercules, the IMF slope is more shallow than a Salpeter IMF (alpha=2.35) at the 5.8-sigma level, and a Kroupa IMF (alpha=2.3 above 0.5 Msun) at 5.4-sigma level. We simultaneously fit for the binary fraction, finding f_binary = 0.47^{+0.16}_{-0.14} for Hercules, and 0.47^{+0.37}_{-0.17} for Leo IV. The UFD binary fractions are consistent with that inferred for Milky Way stars in the same mass range, despite very different metallicities. In contrast, the IMF slopes in the UFDs are shallower than other galactic environments. In the mass range 0.5 - 0.8 Msun, we see a trend across the handful of galaxies with directly measured IMFs such that the power-law slopes become shallower (more bottom-light) with decreasing galactic velocity dispersion and metallicity. This trend is qualitatively consistent with results in elliptical galaxies inferred via indirect methods and is direct evidence for IMF variations with galactic environment.

We present kinematical profiles and metallicity for the M31 dwarf spheroidal (dSph) satellite galaxy Andromeda II (And II) based on Keck DEIMOS spectroscopy of 531 red giant branch stars. Our kinematical sample is among the largest for any M31 satellite and extends out to two effective radii (r_eff = 5.3' = 1.1 kpc). We find a mean systemic velocity of -192.4+-0.5 km/s and an average velocity dispersion of sigma_v = 7.8+-1.1 km/s. While the rotation velocity along the major axis of And II is nearly zero (<1 km/s), the rotation along the minor axis is significant with a maximum rotational velocity of v_max=8.6+-1.8 km/s. We find a kinematical major axis, with a maximum rotational velocity of v_max=10.9+-2.4 km/s, misaligned by 67 degrees to the isophotal major axis. And II is thus the first dwarf galaxy with evidence for nearly prolate rotation with a v_max/sigma_v = 1.1, although given its ellipticity of epsilon = 0.10, this object may be triaxial. We measured metallicities for a subsample of our data, finding a mean metallicity of [Fe/H] = -1.39+- 0.03 dex and an internal metallicity dispersion of 0.72+-0.03 dex. We find a radial metallicity gradient with metal-rich stars more centrally concentrated, but do not observe a significant difference in the dynamics of two metallicity populations. And II is the only known dwarf galaxy to show minor axis rotation making it a unique system whose existence offers important clues on the processes responsible for the formation of dSphs.

We have assembled a large-area spectroscopic survey of giant stars in the Sagittarius (Sgr) dwarf galaxy core. Using medium resolution (R ~15,000), multifiber spectroscopy we have measured velocities of these stars, which extend up to 12 degrees from the galaxy's center (3.7 core radii or 0.4 times the King limiting radius). From these high quality spectra we identify 1310 Sgr members out of 2296 stars surveyed distributed across 24 different fields across the Sgr core. Additional slit spectra were obtained of stars bridging from the Sgr core to its trailing tail. Our systematic, large area sample shows no evidence for significant rotation, a result at odds with the ~20 km/s rotation required as an explanation for the bifurcation seen in the Sgr tidal stream; the observed small (<= 4 km/s) velocity trend along primarily the major axis is consistent with models of the projected motion of an extended body on the sky with no need for intrinsic rotation. The Sgr core is found to have a flat velocity dispersion (except for a kinematically colder center point) across its surveyed extent and into its tidal tails, a property that matches the velocity dispersion profiles measured for other Milky Way dwarf spheroidal (dSph) galaxies. We comment on the possible significance of this observed kinematical similarity for the dynamical state of the other classical Milky Way dSphs in light of the fact that Sgr is clearly a strongly tidally disrupted system.

We present new constraints on the star formation histories of the ultra-faint dwarf (UFD) galaxies, using deep photometry obtained with the Hubble Space Telescope (HST). A galaxy class recently discovered in the Sloan Digital Sky Survey, the UFDs appear to be an extension of the classical dwarf spheroidals to low luminosities, offering a new front in efforts to understand the missing satellite problem. They are the least luminous, most dark-matter dominated, and least chemically-evolved galaxies known. Our HST survey of six UFDs seeks to determine if these galaxies are true fossils from the early universe. We present here the preliminary analysis of three UFD galaxies: Hercules, Leo IV, and Ursa Major I. Classical dwarf spheroidals of the Local Group exhibit extended star formation histories, but these three Milky Way satellites are at least as old as the ancient globular cluster M92, with no evidence for intermediate-age populations. Their ages also appear to be synchronized to within ~1 Gyr of each other, as might be expected if their star formation was truncated by a global event, such as reionization.

We report the discovery of a new ultra-faint globular cluster in the constellation of Ursa Minor, based on stellar photometry from the MegaCam imager at the Canada-France-Hawaii Telescope (CFHT). We find that this cluster, Munoz 1, is located at a distance of 45 +/- 5 kpc and at a projected distance of only 45 arcmin from the center of the Ursa Minor dSph galaxy. Using a Maximum Likelihood technique we measure a half-light radius of 0.5 arcmin, or equivalently 7 pc and an ellipticity consistent with being zero. We estimate its absolute magnitude to be M_V=-0.4 +/- 0.9, which corresponds to L_V=120 (+160, -65) L_sun and we measure a heliocentric radial velocity of -137 +/- 4 km/s based on Keck/DEIMOS spectroscopy. This new satellite is separate from Ursa Minor by ~30 kpc and 110 km/s suggesting the cluster is not obviously associated with the dSph, despite the very close angular separation. Based on its photometric properties and structural parameters we conclude that Munoz 1 is a new ultra-faint stellar cluster. Along with Segue 3 this is one of the faintest stellar clusters known to date.

The discovery of Ultra-Faint Dwarf (UFD) galaxies in the halo of the Milky Way extends the faint end of the galaxy luminosity function to a few hundred solar luminosities. This extremely low luminosity regime poses a significant challenge for the photometric characterization of these systems. We present a suite of simulations aimed at understanding how different observational choices related to the properties of a low luminosity system impact our ability to determine its true structural parameters such as half-light radius and central surface brightness. We focus on estimating half-light radii (on which mass estimates depend linearly) and find that these numbers can have up to 100% uncertainties when relatively shallow photometric surveys, such as SDSS, are used. Our simulations suggest that to recover structural parameters within 10% or better of their true values: (a) the ratio of the field-of-view to the half-light radius of the satellite must be greater than three, (b) the total number of stars, including background objects should be larger than 1000, and (c) the central to background stellar density ratio must be higher than 20. If one or more of these criteria are not met, the accuracy of the resulting structural parameters can be significantly compromised. In the context of future surveys such as LSST, the latter condition will be closely tied to our ability to remove unresolved background galaxies. Assessing the reliability of measured structural parameters will become increasingly critical as the next generation of deep wide-field surveys detects UFDs beyond the reach of current spectroscopic limits.

We investigate the kinematic and photometric properties of the Segue 3 Milky Way companion using Keck/DEIMOS spectroscopy and Magellan/IMACS g and r-band imaging. Using maximum likelihood methods to analyze the photometry, we study the structure and stellar population of Segue 3. We find the half-light radius of Segue 3 is 26" +/- 5" (2.1 +/- 0.4 pc, for a distance of 17 kpc) and the absolute magnitude is a mere M_V = 0.0 +/- 0.8 mag, making Segue 3 the least luminous old stellar system known. We find Segue 3 to be consistent with a single stellar population, with an age of 12.0 +1.5/-0.4 Gyr and an [Fe/H] of -1.7 +0.07/-0.27. Line-of-sight velocities from the spectra are combined with the photometry to determine a sample of 32 stars which are likely associated with Segue 3. The member stars within three half-light radii have a velocity dispersion of 1.2 +/- 2.6 km/s. Photometry of the members indicates the stellar population has a spread in [Fe/H] of <0.3 dex. These facts, together with the small physical size of Segue 3, imply the object is likely an old, faint stellar cluster which contains no significant dark matter. We find tentative evidence for stellar mass loss in Segue 3 through the eleven candidate member stars outside of three half-light radii, as expected from dynamical arguments. Interpretation of the data outside of three half-light radii, is complicated by the object's spatial coincidence with a previously known halo substructure, which may enhance contamination of our member sample.

The Magellanic Clouds are a local laboratory for understanding the evolution and properties of dwarf irregular galaxies. To reveal the extended structure and interaction history of the Magellanic Clouds we have undertaken a large-scale photometric and spectroscopic study of their stellar periphery (the MAgellanic Periphery Survey, MAPS). We present first MAPS results for the Small Magellanic Cloud (SMC): Washington M, T2 + DDO51 photometry reveals metal-poor red giant branch stars in the SMC that extend to large radii (~11 kpc), are distributed nearly azimuthally symmetrically (ellipticity=0.1), and are well-fitted by an exponential profile (out to R~7.5 deg). An ~6 Gyr old, [Fe/H] -1.3 main-sequence turnoff is also evident to at least R=7.3 deg, and as far as 8.4 deg in some directions. We find evidence for a "break" population beyond ~8 radial scalelengths having a very shallow radial density profile that could be either a bound stellar halo or a population of extratidal stars. The distribution of the intermediate stellar component (3<R<7.5 deg) contrasts with that of the inner stellar component (R<3 deg), which is both more elliptical (ellipticity~0.3) and offset from the center of the intermediate component by 0.59 deg, although both components share a similar radial exponential scale length. This offset is likely due to a perspective effect because stars on the eastern side of the SMC are closer on average than stars on the western side. This mapping of its outer stellar structures indicates that the SMC is more complex than previously thought.

We derive an accurate mass estimator for dispersion-supported stellar systems and demonstrate its validity by analyzing resolved line-of-sight velocity data for globular clusters, dwarf galaxies, and elliptical galaxies. Specifically, by manipulating the spherical Jeans equation we show that the dynamical mass enclosed within the 3D deprojected half-light radius r_1/2 can be determined with only mild assumptions about the spatial variation of the stellar velocity dispersion anisotropy. We find M_1/2 = 3 \sigma_los^2 r_1/2 / G ~ 4 \sigma_los^2 R_eff / G, where \sigma_los^2 is the luminosity-weighted square of the line-of-sight velocity dispersion and R_eff is the 2D projected half-light radius. While deceptively familiar in form, this formula is not the virial theorem, which cannot be used to determine accurate masses unless the radial profile of the total mass is known a priori. We utilize this finding to show that all of the Milky Way dwarf spheroidal galaxies (MW dSphs) are consistent with having formed within a halo of mass approximately 3 x 10^9 M_sun in Lambda CDM cosmology. The faintest MW dSphs seem to have formed in dark matter halos that are at least as massive as those of the brightest MW dSphs, despite the almost five orders of magnitude spread in luminosity. We expand our analysis to the full range of observed dispersion-supported stellar systems and examine their I-band mass-to-light ratios (M/L). The M/L vs. M_1/2 relation for dispersion-supported galaxies follows a U-shape, with a broad minimum near M/L ~ 3 that spans dwarf elliptical galaxies to normal ellipticals, a steep rise to M/L ~ 3,200 for ultra-faint dSphs, and a more shallow rise to M/L ~ 800 for galaxy cluster spheroids.

CARMENES, Calar Alto high-Resolution search for M dwarfs with Exo-earths with a Near-infrared Echelle Spectrograph, is a study for a next-generation instrument for the 3.5m Calar Alto Telescope to be designed, built, integrated, and operated by a consortium of nine German and Spanish institutions. Our main objective is finding habitable exoplanets around M dwarfs, which will be achieved by radial velocity measurements on the m/s level in the near-infrared, where low-mass stars emit the bulk of their radiation.

We present deep CFHT/MegaCam photometry of the ultra-faint Milky Way satellite galaxies Coma Berenices (ComBer) and Ursa Major II (UMa II). These data extend to r~25, corresponding to three magnitudes below the main sequence turn-offs in these galaxies. We robustly calculate a total luminosity of M_V=-3.8 +/- 0.6 for ComBer and M_V=-3.9 +/- 0.5 for UMa II, in agreement with previous results. ComBer shows a fairly regular morphology with no signs of active tidal stripping down to a surface brightness limit of 32.4 magarcsec^-2. Using a maximum likelihood analysis, we calculate the half-light radius of ComBer to be r_half=74 +/- 4 pc (5.8 +/- 0.3 arcmin) and its ellipticity e=0.36 +/- 0.04. In contrast, UMa II shows signs of on-going disruption. We map its morphology down to mu_V=32.6 mag arcsec^-2 and found that UMa II is larger than previously determined, extending at least ~700 pc (1.2 deg on the sky) and it is also quite elongated with an ellipticity of e=0.50 +/- 0.2. However, our estimate for the half-light radius, 123 +/- 3 pc (14.1 +/- 0.3 arcmin) is similar to previous results. We discuss the implications of these findings in the context of potential indirect dark matter detections and galaxy formation. We conclude that while ComBer appears to be a stable dwarf galaxy, UMa II shows signs of on-going tidal interaction.

We report the results of a spectroscopic study of the Bootes III (BooIII) stellar overdensity carried out with the Hectospec multifiber spectrograph on the MMT telescope. Radial velocities have been measured for 193 BooIII candidate stars selected to have magnitudes and colors consistent with its upper main sequence and lower red giant branch, as well as a number of horizontal branch candidates. From 20 identified candidate BooIII members, we measure a systemic velocity of V_sun=197.5+-3.8 km/s and a velocity dispersion of sigma_o=14.0+-3.2 km/s. We use the somewhat large velocity dispersion and the implied highly radial orbit, along with morphological evidence from Grillmair (2009) and stellar abundances, to argue that BooIII is likely the first known object observed in a transitional state between being a bound dwarf galaxy and a completely unbound tidal stream.

We report the results of N-body simulations of disrupting satellites aimed at exploring whether the observed features of dSphs can be accounted for with simple, mass-follows-light (MFL) models including tidal disruption. As a test case, we focus on the Carina dwarf Spheroidal (dSph), which presently has the most extensive data at large radius. We find that previous N-body, MFL simulations of dSphs did not sufficiently explore the parameter space of satellite mass, density and orbital shape to find adequate matches to Galactic dSph systems, whereas with a systematic survey of parameter space we are able to find tidally disrupting, MFL satellite models that rather faithfully reproduce Carina's velocity profile, velocity dispersion profile and projected density distribution over its entire sampled radius. The successful MFL model satellites have very eccentric orbits, currently favored by Cold Dark Matter (CDM) models, and central velocity dispersions that still yield an accurate representation of the bound mass and observed central M/L~40 of Carina, despite inflation of the velocity dispersion outside the dSph core by unbound debris. Our survey of parameter space also allows us to address a number of commonly held misperceptions of tidal disruption and its observable effects on dSph structure and dynamics. The simulations suggest that even modest tidal disruption can have a profound effect on the observed dynamics of dSph stars at large radii. Satellites that are well-described by tidally disrupting MFL models could still be fully compatible with L-CDM if for example they represent a later stage in the evolution of luminous subhalos.

We present reliable measurements of the metallicity distribution function (MDF) at different points along the tidal stream of the Sagittarius (Sgr) dwarf spheroidal (dSph) galaxy, based on high resolution, echelle spectroscopy of candidate M giant members of the Sgr system. The Sgr MDF is found to evolve significantly from a median [Fe/H] ~-0.4 in the core to ~-1.1 dex over a Sgr leading arm length representing ~2.5-3.0 Gyr of dynamical (i.e. tidal stripping) age. This is direct evidence that there can be significant chemical differences between current dSph satellites and the bulk of the stars they have contributed to the halo. Our results suggest that Sgr experienced a significant change in binding energy over the past several Gyr, which has substantially decreased its tidal boundary across a radial range over which there must have been a significant metallicity gradient in the progenitor galaxy. By accounting for MDF variation along the debris arms, we approximate the MDF Sgr would have had several Gyr ago. We also analyze the MDF of a moving group of M giants we previously discovered towards the North Galactic Cap having opposite radial velocities to the infalling Sgr leading arm stars there and propose that most of these represent Sgr trailing arm stars overlapping the Sgr leading arm in this part of the sky. If so, these trailing arm stars further demonstrate the strong MDF evolution within the Sgr stream.

We report the results of a spectroscopic study of the Bootes (Boo) dwarf spheroidal (dSph) galaxy carried out with the WIYN telescope and the Hydra multifiber spectrograph. Radial velocities have been measured for 58 Boo candidate stars selected to have magnitudes and colors consistent with its red and asymptotic giant branches. Within the 13' half-light radius, seven members of Boo yield a systemic velocity of V_r=95.6+-3.4 km/s and a velocity dispersion of 6.6+-2.3 km/s. This implies a mass on the order of 1 x 10^7 M_sun, similar to the inferred masses of other Galactic dSphs. Adopting a total Boo luminosity of L=1.8 x 10^4 L_sun to 8.6 x 10^4 L_sun implies M/L ~ 610 to 130, making Boo, the most distorted known Milky Way dwarf galaxy, potentially also the darkest. From the spectra of Boo member stars we estimate its metallicity to be [Fe/H] ~ -2.5, which would make it the most metal poor dSph known to date.

Overdensities in the distribution of low latitude, 2MASS giant stars are revealed by systematically peeling away from sky maps the bulk of the giant stars conforming to ``isotropic'' density laws generally accounting for known Milky Way components. This procedure, combined with a higher resolution treatment of the sky density of both giants and dust allows us to probe to lower Galactic latitudes than previous 2MASS giant star studies. While the results show the swath of excess giants previously associated with the Monoceros ring system in the second and third Galactic quadrants at distances of 6-20 kpc, we also find a several times larger overdensity of giants in the same distance range concentrated in the direction of the ancient constellation Argo. Isodensity contours of the large structure suggest that it is highly elongated and inclined by about 3 deg to the disk, although details of the structure -- including the actual location of highest density, overall extent, true shape -- and its origin, remain unknown because only a fraction of it lies outside highly dust-obscured, low latitude regions. Nevertheless, our results suggest that the 2MASS M giant overdensity previously claimed to represent the core of a dwarf galaxy in Canis Major (l ~ 240 deg) is an artifact of a dust extinction window opening to the overall density rise to the more significant Argo structure centered at larger longitude (l ~ 290 +- 10 deg, b ~ -4 +- 2 deg).

A new large-area Washington M,T_2+DDO51 filter survey of more than 10 deg^2 around the Carina dSph galaxy reveals a spectroscopically confirmed power law radial density "break" population of Carina giant stars extending several degrees beyond the central King profile. Magellan telescope MIKE spectroscopy establishes the existence of Carina stars to at least 4.5 times its central King limiting radius, r_lim and primarily along Carina's major axis. To keep these stars bound to the dSph would require a global Carina mass-to-light ratio of M/L > 6,300 M/L_sun. The MIKE velocities, supplemented with ~950 additional Carina field velocities from archived VLT+GIRAFFE spectra with r<=r_lim, demonstrate a nearly constant Carina velocity dispersion to just beyond r = r_lim, and both a rising velocity dispersion and a velocity shear at still larger radii. Together, the observational evidence suggests that the discovered extended Carina population represents tidal debris from the dSph. Of 65 giant candidates at large angular radii from the Carina center for which MIKE spectra have been obtained 94% are associated either with Carina or a second, newly discovered diffuse, but strongly radial velocity-coherent (velocity dispersion of 9.8 km s^-1), foreground halo system. The fifteen stars in this second, retrograde velocity population have (1) a mean metallicity ~1 dex higher than that of Carina, and (2) colors and magnitudes consistent with the red clump of the Large Magellanic Cloud (LMC). Additional spectroscopy of giant star candidates in fields linking Carina and the LMC shows a smooth velocity gradient between the LMC and the retrograde Carina moving group. We conclude that we have found Magellanic stars almost twice as far (22 deg) from the LMC center than previously known.

We analyze velocity dispersion profiles for the Draco and Ursa Minor (UMi) dwarf spheroidal (dSph) galaxies based on published and new Keck HIRES spectra for stars in the outer UMi field. Washington+DDO51 filter photometric catalogs provide additional leverage on membership of individual stars, and beyond 0.5 King limiting radii (R_lim) identify bona fide dSph members up to 4.5 times more efficiently than simple color-magnitude diagram selections. Previously reported ``cold populations'' R_lim are not obvious in the data and appear only with particular binning; more or less constant and platykurtic dispersion profiles are characteristic of these dSphs to large radii. We report discovery of UMi stars to at least 2.7 R_lim (i.e.,210 arcmin or 4 kpc). Even with conservative assumptions, a UMi mass of M > 4.9 x 10^8 M_(sun) is required to bind these stars, implying an unlikely global mass-to-light ratio of M/L > 900 (M/L)_(sun). We conclude that we have found stars tidally stripped from UMi.

It has been previously proposed that some distant open clusters in the Milky Way may have been accreted during a dwarf galaxy merger, perhaps associated with the same event that led to the formation of the Galactic anticenter stellar structure (GASS), also known as the ``Monoceros Ring''. We have obtained VI and Washington+DDO51 photometric and medium resolution (R ~ 8000) multi-fiber spectroscopic data for the three distant old open clusters Berkeley 29, Saurer 1, and vdB-Hagen 176 (BH 176). These clusters are spatially coincident with GASS, but radial velocities and spectroscopic metallicities had not been available during previous studies of the GASS candidate cluster system. Similar data for the clusters Berkeley 20 and Berkeley 39 have been obtained for calibration purposes. We provide the first {\it reliable} radial velocity for BH 176 (V_{helio} = 11.2 +/- 5.3 km/s). We also find that V_{helio} = +95.4 +/- 3.6 and +28.4 +/- 3.6 km/s, for Saurer 1(A) and Berkeley 29, respectively. We show that alpha-enhanced isochrones, while spectroscopically motivated, provide a poor fit to Be29 in contrast to previous findings. We find that the clusters Berkeley 29 and Saurer 1 are consistent with the previously reported characteristics for GASS candidate clustersand the GASS stellar stream as derived from M-giant observations. However, the radial velocity and photometric metallicity ([Fe/H] ~ 0.0 dex) for BH 176 suggests that a connection of this cluster with the putative GASS cluster system is unlikely. We reassess the age-metallicity relation for the most likely members of the GASS clusters system for which spectroscopic metallicities are now available.

The question of the existence of active tidal disruption around various dSph galaxies remains controversial. That debate often centers on the nature (bound vs. unbound) of extended populations of stars. However, the more fundamental issue of the very existence of the extended populations is still contentious. We present an evaluation of the debate centering on one particular dSph, Carina, for which claims both for and against the existence of stars beyond the King radius have been made. Our review includes an examination of all previous studies bearing on the Carina radial profile and shows that the survey method which achieves the highest detected dSph signal-to-background in the outer parts of the galaxy is the Washington M, T2 + DDO51 (MTD) filter approach from Paper II in this series. We then address statistical methods used to evaluate the reliability of MTD surveys in the presence of photometric errors and for which a new, a posteriori statistical analysis methodology is provided. Finally, these statistical methods are tested by new spectroscopy of stars in the MTD-selected Carina candidate sample. Of 74 candidate giants with follow-up spectroscopy, the MTD technique identified 61 new Carina members, including 8 stars outside the King radius. From a sample of 29 stars not initially identified as candidate Carina giants but that lie just outside of our selection criteria, 12 have radial velocities consistent with membership, including 5 extratidal stars. Carina is shown to have an extended population of giant stars extending to a major axis radius of 40' (1.44x the nominal King radius).

The Galactic Anticenter Stellar Structure (GASS) has been identified with excess surface densities of field stars in several large area sky surveys, and with an unusual, string-like grouping of star clusters. Some members of the cluster grouping have radial velocities (RVs) consistent with the observed GASS velocity-longitude trend. We provide new RV measurements of stars in six clusters that have been suggested to be associated with the GASS. We show that the RVs of at least four clusters are consistent with the previously measured RV trend for GASS. We also derive spectroscopic metallicities for four clusters, and provide an improved age-metallicity relation for the clusters apparently associated with GASS.

The Galactic Anticenter Stellar Structure (GASS) has been identified with excess surface densities of field stars in several large area sky surveys, and with an unusual, string-like grouping of five globular clusters. At least two of these are diffuse, young ``transitional'' clusters between open and globular types. Here we call attention to the fact that four {\it younger} open or transitional clusters extend the previously identified, string-like cluster grouping, with at least one having a radial velocity consistent with the previously found GASS velocity-longitude trend. All nine clusters lie close to a plane tipped 17$\arcdeg$ to the Galactic plane. This planar orientation is used to forage for additional potential cluster members in the inner Galaxy, and a number are found along the same plane and string-like sequence, {\it including almost all fifteen known outer, old open clusters}. Tidal accretion of a dwarf satellite galaxy on a low inclination orbit -- perhaps the GASS system -- appears to be a plausible explanation for the origin of the outer, old open and transitional clusters of the Milky Way. We use these clusters to explore the age-metallicity relation of the putative accreted GASS progenitor. Finally, we provide the first radial velocity of a star in the cluster BH 176 and discuss its implications.