This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

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.

Interesting chemically peculiar field stars may reflect their stellar evolution history and their possible origin in a different environment from where they are found now, which is one of the most important research fields in Galactic archaeology. To explore this further, we have used the CN-CH bands around 4000 A to identify N-rich metal-poor field stars in LAMOST DR3. Here we expand our N-rich metal-poor field star sample to ~100 stars in LAMOST DR5, where 53 of them are newly found in this work. We investigate light elements of the common stars between our sample and APOGEE DR14. While Mg, Al, and Si abundances generally agree with the hypothesis that N-rich metal-poor field stars come from enriched populations in globular clusters, it is still inconclusive for C, N, and O. After integrating the orbits of our N-rich field stars and a control sample of normal metal-poor field stars, we find that N-rich field stars have different orbital parameter distributions compared to the control sample, specifically, apocentric distances, maximum vertical amplitude (Zmax), orbital energy, and z direction angular momentum (Lz). The orbital parameters of N-rich field stars indicate that most of them are inner-halo stars. The kinematics of N-rich field stars support their possible GC origin. The spatial and velocity distributions of our bona fide N-rich field star sample are important observational evidence to constrain simulations of the origin of these interesting objects.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

The large amount of chemical and kinematic information available in large spectroscopic surveys have inspired the search for chemically peculiar stars in the field. Though these metal-poor field stars ([Fe/H$]<-1$) are commonly enriched in nitrogen, their detailed spatial, kinematic, and chemical distributions suggest that various groups may exist, and thus their origin is still a mystery. To study these stars statistically, we increase the sample size by identifying new CN-strong stars with LAMOST DR3 for the first time. We use CN-CH bands around 4000 \AA~to find CN-strong stars, and further separate them into CH-normal stars (44) and CH-strong (or CH) stars (35). The chemical abundances from our data-driven software and APOGEE DR 14 suggest that most CH-normal stars are N-rich, and it cannot be explained by only internal mixing process. The kinematics of our CH-normal stars indicate a substantial fraction of these stars are retrograding, pointing to an extragalactic origin. The chemistry and kinematics of CH-normal stars imply that they may be GC-dissolved stars, or accreted halo stars, or both.

We map the trends of elemental abundance ratios across the Galactic disk, spanning R = 3-15 kpc and midplane distance |Z|= 0-2 kpc, for 15 elements in a sample of 20,485 stars measured by the SDSS/APOGEE survey (O, Na, Mg, Al, Si, P, S, K, Ca, V, Cr, Mn, Fe, Co, Ni). Adopting Mg rather than Fe as our reference element, and separating stars into two populations based on [Fe/Mg], we find that the median trends of [X/Mg] vs. [Mg/H] in each population are nearly independent of location in the Galaxy. The full multi-element cartography can be summarized by combining these nearly universal median sequences with our measured metallicity distribution functions and the relative proportions of the low-[Fe/Mg] (high-alpha) and high-[Fe/Mg] (low-alpha) populations, which depend strongly on R and |Z|. We interpret the median sequences with a semi-empirical "2-process" model that describes both the ratio of core collapse and Type Ia supernova contributions to each element and the metallicity dependence of the supernova yields. These observationally inferred trends can provide strong tests of supernova nucleosynthesis calculations. Our results lead to a relatively simple picture of abundance ratio variations in the Milky Way, in which the trends at any location can be described as the sum of two components with relative contributions that change systematically and smoothly across the Galaxy. Deviations from this picture and future extensions to other elements can provide further insights into the physics of stellar nucleosynthesis and unusual events in the Galaxy's history.

We present reverberation mapping results from the first year of combined spectroscopic and photometric observations of the Sloan Digital Sky Survey Reverberation Mapping Project. We successfully recover reverberation time delays between the $g+i$-band emission and the broad H$\beta$ emission line for a total of 44 quasars, and for the broad H$\alpha$ emission line in 18 quasars. Time delays are computed using the JAVELIN and CREAM software and the traditional interpolated cross-correlation function (ICCF): Using well defined criteria, we report measurements of 32 H$\beta$ and 13 H$\alpha$ lags with JAVELIN, 42 H$\beta$ and 17 H$\alpha$ lags with CREAM, and 16 H$\beta$ and 8 H$\alpha$ lags with the ICCF. Lag values are generally consistent among the three methods, though we typically measure smaller uncertainties with JAVELIN and CREAM than with the ICCF, given the more physically motivated light curve interpolation and more robust statistical modeling of the former two methods. The median redshift of our H$\beta$-detected sample of quasars is 0.53, significantly higher than that of the previous reverberation-mapping sample. We find that in most objects, the time delay of the H$\alpha$ emission is consistent with or slightly longer than that of H$\beta$. We measure black hole masses using our measured time delays and line widths for these quasars. These black hole mass measurements are mostly consistent with expectations based on the local M-sigma relationship, and are also consistent with single-epoch black hole mass measurements. This work increases the current sample size of reverberation-mapped active galaxies by about two-thirds and represents the first large sample of reverberation mapping observations beyond the local universe (z < 0.3).

Data and analysis methodology used for the SDSS/APOGEE Data Releases 13 and 14 are described, highlighting differences from the DR12 analysis presented in Holtzman (2015). Some improvement in the handling of telluric absorption and persistence is demonstrated. The derivation and calibration of stellar parameters, chemical abundances, and respective uncertainties are described, along with the ranges over which calibration was performed. Some known issues with the public data related to the calibration of the effective temperatures (DR13), surface gravity (DR13 and DR14), and C and N abundances for dwarfs (DR13 and DR14) are highlighted. We discuss how results from a data-driven technique, The Cannon (Casey 2016), are included in DR14, and compare those with results from the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP). We describe how using The Cannon in a mode that restricts the abundance analysis of each element to regions of the spectrum with known features from that element leads to Cannon abundances can lead to significantly different results for some elements than when all regions of the spectrum are used to derive abundances.

We use deep multi-epoch near-IR images of the VISTA Variables in the Via Lactea (VVV) Survey to measure proper motions (PMs) of stars in the Milky Way globular cluster FSR1716 = VVV-GC05. The color-magnitude diagram of this object, made using PM selected members, shows an extended horizontal branch, nine confirmed RR Lyrae members in the instability strip, and possibly several hotter stars extending to the blue. Based on the fundamental-mode (ab-type) RR Lyrae stars that move coherently with the cluster, we confirmed that FSR1716 is an Oosterhoff I globular cluster with a mean period Pab = 0.574 days. Intriguingly, we detect tidal extensions to both sides of this cluster in the spatial distribution of PM selected member stars. Also, one of the confirmed RRabs is located 11 arcmin in projection from the cluster center, suggesting that FSR1716 may be losing stars due to the gravitational interaction with the Galaxy. We also measure radial velocities (RVs) for five cluster red giants selected using the PMs. The combination of RVs and PMs allow us to compute for the first time the orbit of this globular cluster, using an updated Galactic potential. The orbit results to be confined within |Zmax | < 2.0 kpc, and has eccentricity 0.4 < e < 0.6, with perigalactic distance 1.5 < Rperi (kpc) < 2.3, and apogalactic distance 5.3 < Rapo (kpc) < 6.4. We conclude that, in agreement with its relatively low metallicity ([Fe/H]= -1.4 dex), this is an inner halo globular cluster plunging into the disk of the Galaxy. As such, this is a unique object to test the dynamical processes that contribute to the disruption of Galactic globular clusters.

Binary stars make up a significant portion of all stellar systems. Consequently, an understanding of the bulk properties of binary stars is necessary for a full picture of star formation. Binary surveys indicate that both multiplicity fraction and typical orbital separation increase as functions of primary mass. Correlations with higher order architectural parameters such as mass ratio are less well constrained. We seek to identify and characterize double-lined spectroscopic binaries (SB2s) among the 1350 M dwarf ancillary science targets with APOGEE spectra in the SDSS-III Data Release 13. We measure the degree of asymmetry in the APOGEE pipeline cross-correlation functions (CCFs), and use those metrics to identify a sample of 44 high-likelihood candidate SB2s. At least 11 of these SB2s are known, having been previously identified by Deshapnde et al, and/or El Badry et al. We are able to extract radial velocities (RVs) for the components of 36 of these systems from their CCFs. With these RVs, we measure mass ratios for 29 SB2s and 5 SB3s. We use Bayesian techniques to fit maximum likelihood (but still preliminary) orbits for 4 SB2s with 8 or more distinct APOGEE observations. The observed (but incomplete) mass ratio distribution of this sample rises quickly towards unity. Two-sided Kolmogorov-Smirnov tests and probabilities of 18.3% and 18.7%, demonstrating that the mass ratio distribution of our sample is consistent with those measured by Pourbaix et al. and Fernandez et al., respectively.

Bars are common galactic structures in the local universe that play an important role in the secular evolution of galaxies, including the Milky Way. In particular, the velocity distribution of individual stars in our galaxy is useful to shed light on stellar dynamics, and provides information complementary to that inferred from the integrated light of external galaxies. However, since a wide variety of models reproduce the distribution of velocity and the velocity dispersion observed in the Milky Way, we look for signatures of the bar on higher-order moments of the line-of-sight velocity ($V_{los}$) distribution. We make use of two different numerical simulations --one that has developed a bar and one that remains nearly axisymmetric-- to compare them with observations in the latest APOGEE data release (SDSS DR14). This comparison reveals three interesting structures that support the notion that the Milky Way is a barred galaxy. A high skewness region found at positive longitudes constrains the orientation angle of the bar, and is incompatible with the orientation of the bar at $\ell=0^\circ$ proposed in previous studies. We also analyse the $V_{los}$ distributions in three regions, and introduce the Hellinger distance to quantify the differences among them. Our results show a strong non-Gaussian distribution both in the data and in the barred model, confirming the qualitative conclusions drawn from the velocity maps. In contrast to earlier work, we conclude it is possible to infer the presence of the bar from the kurtosis distribution.

The nature of the Triangulum-Andromeda (TriAnd) system has been debated since the discovery of this distant, low-latitude Milky Way (MW) overdensity more than a decade ago. Explanations for its origin are either as a halo substructure from the disruption of a dwarf galaxy or a distant extension of the Galactic disk. We test these hypotheses using chemical abundances of a dozen TriAnd members from the Sloan Digital Sky Survey's 14th Data Release of Apache Point Observatory Galactic Evolution Experiment (APOGEE) data to compare to APOGEE abundances of stars with similar metallicity from both the Sagittarius (Sgr) dSph, and the outer MW disk. We find that TriAnd stars are chemically distinct from Sgr across a variety of elements, (C+N), Mg, K, Ca, Mn, and Ni, with a separation in [X/Fe] of about 0.1 to 0.4 dex depending on the element. Instead, the TriAnd stars, with a median metallicity of about -0.8, exhibit chemical abundance ratios similar to those of the lowest metallicity ([Fe/H] ~ -0.7) stars in the outer Galactic disk, and are consistent with expectations of extrapolated chemical gradients in the outer disk of the MW. These results suggest that TriAnd is associated with the MW disk, and, therefore, that the disk extends to this overdensity --- i.e., past a Galactocentric radius of 24 kpc --- albeit vertically perturbed about 7 kpc below the nominal disk midplane in this region of the Galaxy.

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.

Detailed chemical abundance distributions for fourteen elements are derived for eight high-probability stellar members of the solar metallicity old open cluster M67 with an age of $\sim$4 Gyr. The eight stars consist of four pairs, with each pair occupying a distinct phase of stellar evolution: two G-dwarfs, two turnoff stars, two G-subgiants, and two red clump K-giants. The abundance analysis uses near-IR high-resolution spectra ($\lambda$1.5 -- 1.7$\mu$m) from the APOGEE survey and derives abundances for C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, and Fe. Our derived stellar parameters and metallicity for 2M08510076+113115 suggest that this star is a solar-twin, exhibiting abundance differences relative to the Sun of $\leq$ 0.04 dex for all elements. Chemical homogeneity is found within each class of stars ($\sim$0.02 dex), while significant abundance variations ($\sim$0.05 -- 0.20 dex) are found across the different evolutionary phases; the turnoff stars typically have the lowest abundances, while the red clump tend to have the largest. Non-LTE corrections to the LTE-derived abundances are unlikely to explain the differences. A detailed comparison of the derived Fe, Mg, Si, and Ca abundances with recently published surface abundances from stellar models that include chemical diffusion, provides a good match between the observed and predicted abundances as a function of stellar mass. Such agreement would indicate the detection of chemical diffusion processes in the stellar members of M67.

We study the radial acceleration relation (RAR) for early-type galaxies (ETGs) in the SDSS MaNGA MPL5 dataset. The complete ETG sample show a slightly offset RAR from the relation reported by McGaugh et al. (2016) at the low-acceleration end; we find that the deviation is due to the fact that the slow rotators show a systematically higher acceleration relation than the McGaugh's RAR, while the fast rotators show a consistent acceleration relation to McGaugh's RAR. There is a 1\sigma significant difference between the acceleration relations of the fast and slow rotators, suggesting that the acceleration relation correlates with the galactic spins, and that the slow rotators may have a different mass distribution compared with fast rotators and late-type galaxies. We suspect that the acceleration relation deviation of slow rotators may be attributed to more galaxy merger events, which would disrupt the original spins and correlated distributions of baryons and dark matter orbits in galaxies.

In our previous work, we found that only two scenarios are capable of reproducing the observed integrated mass-metallicity relations for the gas and stellar components of local star-forming galaxies simultaneously. One scenario invokes a time-dependent metal outflow loading factor with stronger outflows at early times. The other scenario uses a time-dependent IMF slope with a steeper IMF at early times. In this work, we extend our study to investigate the radial profile of gas and stellar metallicity in local star-forming galaxies using spatially resolved spectroscopic data from the SDSS-IV MaNGA survey. We find that most galaxies show negative gradients in both gas and stellar metallicity with steeper gradients in stellar metallicity. The stellar metallicity gradient tend to be mass dependent with steeper gradients in more massive galaxies while no clear mass dependence is found for the gas metallicity gradient. Then we compare the observations with the predictions from a chemical evolution model of the radial profiles of gas and stellar metallicities. We confirm that the two scenarios proposed in our previous work are also required to explain the metallicity gradients. Based on these two scenarios we successfully reproduce the radial profiles of gas metallicity, stellar metallicity, stellar mass surface density, and star formation rate (SFR) surface density simultaneously. The origin of the negative gradient in stellar metallicity turns out to be driven by either radially dependent metal outflow or IMF slope. In contrast, the radial dependence of the gas metallicity is less constrained because of the degeneracy in model parameters.

In this work we make use of DR14 APOGEE spectroscopic data to study a sample of 92 known OB stars. We developed a near-infrared semi-empirical spectral classification method that was successfully used in case of four new exemplars, previously classified as later B-type stars. Our results agree well with those determined independently from ECHELLE optical spectra, being in line with the spectral types derived from the "canonical" MK blue optical system. This confirms that the APOGEE spectrograph can also be used as a powerful tool in surveys aiming to unveil and study large number of moderately and highly obscured OB stars still hidden in the Galaxy.

Metal-poor globular clusters (GCs) exhibit intriguing Al-Mg anti-correlations and possible Si-Al correlations, which are important clues to decipher the multiple-population phenomenon. NGC 5053 is one of the most metal-poor GCs in the nearby Universe, and has been suggested to be associated with the Sagittarius (Sgr) dwarf galaxy, due to its similarity in location and radial velocity with one of the Sgr arms. In this work, we simulate the orbit of NGC 5053, and argue against a physical connection between Sgr and NGC 5053. On the other hand, the Mg, Al, and Si spectral lines, which are difficult to detect in the optical spectra of NGC 5053 stars, have been detected in the near-infrared APOGEE spectra. We use three different sets of stellar parameters and codes to derive the Mg, Al, and Si abundances. Regardless of which method is adopted, we see a large Al variation, and a substantial Si spread. Along with NGC 5053, metal-poor GCs exhibit different Mg, Al, and Si variations. Moreover, NGC 5053 has the lowest cluster mass among the GCs that have been identified to exhibit an observable Si spread until now.

The Hercules stream is a group of co-moving stars in the Solar neighbourhood, which can potentially be explained as a signature of either the outer Lindblad resonance (OLR) of a fast Galactic bar or the corotation resonance of a slower bar. In either case, the feature should be present over a large area of the disc. With the recent commissioning of the APOGEE-2 Southern spectrograph we can search for the Hercules stream at $(l,b)=(270^\circ,0)$, a direction in which the Hercules stream, if caused by the bar's OLR, would be strong enough to be detected using only the line-of-sight velocities. We clearly detect a narrow, Hercules-like feature in the data that can be traced from the solar neighbourhood to a distance of about 4 kpc. The detected feature matches well the line-of-sight velocity distribution from the fast-bar (OLR) model. Confronting the data with a model where the Hercules stream is caused by the corotation resonance of a slower bar leads to a poorer match, as the corotation model does not predict clearly separated modes, possibly because the slow-bar model is too hot.

A cold high-velocity (HV, $\sim$ 200 km/s) peak was first reported in several Galactic bulge fields based on the APOGEE commissioning observations. Both the existence and the nature of the high-velocity peak are still under debate. Here we revisit this feature with the latest APOGEE DR13 data. We find that most of the low latitude bulge fields display a skewed Gaussian distribution with a HV shoulder. However, only 3 out of 53 fields show distinct high-velocity peaks around 200 km/s. The velocity distribution can be well described by Gauss-Hermite polynomials, except the three fields showing clear HV peaks. We find that the correlation between the skewness parameter ($h_{3}$) and the mean velocity ($\bar{v}$), instead of a distinctive HV peak, is a strong indicator of the bar. It was recently suggested that the HV peak is composed of preferentially young stars. We choose three fields showing clear HV peaks to test this hypothesis using the metallicity, [$\alpha$/M] and [C/N] as age proxies. We find that both young and old stars show HV features. The similarity between the chemical abundances of stars in the HV peaks and the main component indicates that they are not systematically different in terms of chemical abundance or age. In contrast, there are clear differences in chemical space between stars in the Sagittarius dwarf and the bulge stars. The strong HV peaks off-plane are still to be explained properly, and could be different in nature.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

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 take advantage of the Gaia-ESO Survey iDR4 bulge data to search for abundance anomalies that could shed light on the composite nature of the Milky Way bulge. The alpha-elements (Mg, Si, and whenever available, Ca) abundances, and their trends with Fe abundances have been analysed for a total of 776 bulge stars. In addition, the aluminum abundances and their ratio to Fe and Mg have also been examined. Our analysis reveals the existence of low-alpha element abundance stars with respect to the standard bulge sequence in the [alpha/Fe] vs. [Fe/H] plane. 18 objects present deviations in [alpha/Fe] ranging from 2.1 to 5.3 sigma with respect to the median standard value. Those stars do not show Mg-Al anti-correlation patterns. Incidentally, this sign of the existence of multiple stellar populations is reported firmly for the first time for the bulge globular cluster NGC 6522. The identified low-alpha abundance stars have chemical patterns compatible with those of the thin disc. Their link with massive dwarf galaxies accretion seems unlikely, as larger deviations in alpha abundance and Al would be expected. The vision of a bulge composite nature and a complex formation process is reinforced by our results. The used approach, a multi-method and model-driven analysis of high resolution data seems crucial to reveal this complexity.

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.

Context. M-type subdwarfs are metal-poor low-mass stars and probe for the old populations in our Galaxy. Accurate knowledge of their atmospheric parameters and especially the composition is essential for understanding the chemical history of our Galaxy. Aims. The purpose of this work is to perform a detailed study of M-subdwarf spectra covering the full wavelength range from the optical to the near-infrared. It allows to do a more detailed analysis of the atmospheric composition in order to determine the stellar parameters, and to constrain the atmospheric models. The study will allow us to further understand physical and chemical processes such as increasing condensation of gas into dust, to point out the missing continuum opacities and see how the main band features are reproduced by the models. The spectral resolution and the large wavelength coverage used is a unique combination to constrain the process that occur in cool atmosphere. Methods. We obtained medium-resolution (R = 5000-7000) spectra over the wavelength range 0.3-2.5 {\mu}m of ten M-type subdwarfs with X-SHOOTER at VLT. These data constitute an unique atlas of M-subdwarfs from optical to near-infrared. We performed spectral synthesis analysis using a full grid of synthetic spectra computed from BT-Settl models and obtained consistent stellar parameters such as effective temperature, surface gravity and metallicity. Results. We show that state-of the-art atmospheric models correctly represent the overall shape of their spectral energy distribution, as well as atomic and molecular line profiles both in the optical and near-infrared. We found that the actual fitted gravities of almost all our sample are consistent with old objects, except for LHS 73 where it is found surprisingly low.

We present preliminary results of the CIDA Equatorial Variability Survey (CEVS), looking for quasar (hereafter QSO) candidates near the Galactic plane. The CEVS contains photometric data from extended and adjacent regions of the Milky Way disk ($\sim$ 500 sq. deg.). In this work 2.5 square degrees with moderately high temporal sampling in the CEVS were analyzed. The selection of QSO candidates was based on the study of intrinsic optical photometric variability of 14,719 light curves. We studied samples defined by cuts in the variability index (Vindex $>$ 66.5), periodicity index (Q $>$ 2), and the distribution of these sources in the plane (AT , ${\gamma}$), using a slight modification of the first-order of the structure function for the temporal sampling of the survey. Finally, 288 sources were selected as QSO candidates. The results shown in this work are a first attempt to develop a robust method to detect QSO towards the Galactic plane in the era of massive surveys such as VISTA and Gaia.

We estimate the expected distribution of displacements between the two dominant dark matter (DM) peaks (DM-DM displacements) and between DM and gaseous baryon peak (DM-gas displacements) in dark matter halos with masses larger than $10^{13}$ Msun/h. We use as a benchmark the observation of SL2S J08544-0121, which is the lowest mass system ($1.0\times 10^{14}$ Msun/h) observed so far featuring a bi-modal dark matter distribution with a dislocated gas component. We find that $(50 \pm 10)$% of the dark matter halos with circular velocities in the range 300 km/s to 700 km/s (groups) show DM-DM displacements equal or larger than $186 \pm 30$ kpc/h as observed in SL2S J08544-0121. For dark matter halos with circular velocities larger than 700 km/s (clusters) this fraction rises to 70 $\pm$ 10%. Using the same simulation we estimate the DM-gas displacements and find that 0.1 to 1.0% of the groups should present separations equal or larger than $87\pm 14$kpc/h corresponding to our observational benchmark; for clusters this fraction rises to (7 $\pm$ 3)%, consistent with previous studies of dark matter to baryon separations. Considering both constraints on the DM-DM and DM-gas displacements we find that the number density of groups similar to SL2S J08544-0121 is $\sim 6.0\times 10^{-7}$ Mpc$^{-3}$, three times larger than the estimated value for clusters. These results open up the possibility for a new statistical test of LCDM by looking for DM-gas displacements in low mass clusters and groups.