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).

The Open Cluster Chemical Abundances and Mapping (OCCAM) survey aims to constrain key Galactic dynamical and chemical evolution parameters by the construction of a large, comprehensive, uniform, infrared-based spectroscopic data set of hundreds of open clusters. This fourth contribution from the OCCAM survey presents analysis using SDSS/APOGEE DR16 of a sample of 128 open clusters, 71 of which we designate to be "high quality" based on the appearance of their color-magnitude diagram. We find the APOGEE DR16 derived [Fe/H] abundances to be in good agreement with previous high resolution spectroscopic open cluster abundance studies. Using the high quality sample, we measure Galactic abundance gradients in 16 elements, and find evolution of some of the [X/Fe] gradients as a function of age. We find an overall Galactic [Fe/H] vs R_GC gradient of $-0.068 \pm 0.001$ dex kpc$^{-1}$ over the range of $6 <$ R_GC $< 13.9$ kpc; however, we note that this result is sensitive to the distance catalog used, varying as much as 15%. We formally derive the location a break in the [Fe/H] abundance gradient as a free parameter in the gradient fit for the first time. We also measure significant Galactic gradients in O, Mg, S, Ca, Mn, Cr, Cu, Na, Al, and K, some of which are measured for the first time. Our large sample allows us to explore four well-populated age bins to explore the time evolution of gradients for a large number of elements and comment on possible implications for Galactic chemical evolution and radial migration.

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.

We present spectroscopic determinations of the effective temperatures, surface gravities and metallicities for 21 M-dwarfs observed at high-resolution (R $\sim$ 22,500) in the \textit{H}-band as part of the SDSS-IV APOGEE survey. The atmospheric parameters and metallicities are derived from spectral syntheses with 1-D LTE plane parallel MARCS models and the APOGEE atomic/molecular line list, together with up-to-date H$_{2}$O and FeH molecular line lists. Our sample range in $T_{\rm eff}$ from $\sim$ 3200 to 3800K, where eleven stars are in binary systems with a warmer (FGK) primary, while the other 10 M-dwarfs have interferometric radii in the literature. We define an $M_{K_{S}}$--Radius calibration based on our M-dwarf radii derived from the detailed analysis of APOGEE spectra and Gaia DR2 distances, as well as a mass-radius relation using the spectroscopically-derived surface gravities. A comparison of the derived radii with interferometric values from the literature finds that the spectroscopic radii are slightly offset towards smaller values, with $\Delta$ = -0.01 $\pm$ 0.02 $R{\star}$/$R_{\odot}$. In addition, the derived M-dwarf masses based upon the radii and surface gravities tend to be slightly smaller (by $\sim$5-10\%) than masses derived for M-dwarf members of eclipsing binary systems for a given stellar radius. The metallicities derived for the 11 M-dwarfs in binary systems, compared to metallicities obtained for their hotter FGK main-sequence primary stars from the literature, shows excellent agreement, with a mean difference of [Fe/H](M-dwarf - FGK primary) = +0.04 $\pm$ 0.18 dex, confirming the APOGEE metallicity scale derived here for M-dwarfs.

Chemical abundances for 15 elements (C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni) are presented for 83 stellar members of the 4 Gyr old solar-metallicity open cluster M67. The sample contains stars spanning a wide range of evolutionary phases, from G dwarfs to red clump stars. The abundances were derived from near-IR ($\lambda$1.5 -- 1.7$\mu$m) high-resolution spectra ($R$ = 22,500) from the SDSS-IV/APOGEE survey. A 1-D LTE abundance analysis was carried out using the APOGEE synthetic spectral libraries, via chi-square minimization of the synthetic and observed spectra with the qASPCAP code. We found significant abundance differences ($\sim$0.05 -- 0.30 dex) between the M67 member stars as a function of the stellar mass (or position on the HR diagram), where the abundance patterns exhibit a general depletion (in [X/H]) in stars at the main-sequence turnoff. The amount of the depletion is different for different elements. We find that atomic diffusion models provide, in general, good agreement with the abundance trends for most chemical species, supporting recent studies indicating that measurable atomic diffusion operates in M67 stars.

Internal mixing on the giant branch is an important process which affects the evolution of stars and the chemical evolution of the galaxy. While several mechanisms have been proposed to explain this mixing, better empirical constraints are necessary. Here, we use [C/N] abundances in 26097 evolved stars from the SDSS-IV/APOGEE-2 Data Release 14 to trace mixing and extra mixing in old field giants with -1.7< [Fe/H] < 0.1. We show that the APOGEE [C/N] ratios before any dredge-up occurs are metallicity dependent, but that the change in [C/N] at the first dredge-up is metallicity independent for stars above [Fe/H] ~ -1. We identify the position of the red giant branch (RGB) bump as a function of metallicity, note that a metallicity-dependent extra mixing episode takes place for low-metallicity stars ([Fe/H] <-0.4) 0.14 dex in log g above the bump, and confirm that this extra mixing is stronger at low metallicity, reaching $\Delta$ [C/N] = 0.58 dex at [Fe/H] = -1.4. We show evidence for further extra mixing on the upper giant branch, well above the bump, among the stars with [Fe/H] < -1.0. This upper giant branch mixing is stronger in the more metal-poor stars, reaching 0.38 dex in [C/N] for each 1.0 dex in log g. The APOGEE [C/N] ratios for red clump (RC) stars are significantly higher than for stars at the tip of the RGB, suggesting additional mixing processes occur during the helium flash or that unknown abundance zero points for C and N may exist among the red clump RC sample. Finally, because of extra mixing, we note that current empirical calibrations between [C/N] ratios and ages cannot be naively extrapolated for use in low-metallicity stars specifically for those above the bump in the luminosity function.

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.

We present an overview of the distributions of 11 elemental abundances in the Milky Way's inner regions, as traced by APOGEE stars released as part of SDSS Data Release 14/15 (DR14/DR15), including O, Mg, Si, Ca, Cr, Mn, Co, Ni, Na, Al, and K. This sample spans ~4000 stars with R_GC<4 kpc, enabling the most comprehensive study to date of these abundances and their variations within the innermost few kiloparsecs of the Milky Way. We describe the observed abundance patterns ([X/Fe]-[Fe/H]), compare to previous literature results and to patterns in stars at the solar Galactic radius, and discuss possible trends with DR14/DR15 effective temperatures. We find that the position of the [Mg/Fe]-[Fe/H] "knee" is nearly constant with R_GC, indicating a well-mixed star-forming medium or high levels of radial migration in the early inner Galaxy. We quantify the linear correlation between pairs of elements in different subsamples of stars and find that these relationships vary; some abundance correlations are very similar between the alpha-rich and alpha-poor stars, but others differ significantly, suggesting variations in the metallicity dependencies of certain supernova yields. These empirical trends will form the basis for more detailed future explorations and for the refinement of model comparison metrics. That the inner Milky Way abundances appear dominated by a single chemical evolutionary track and that they extend to such high metallicities underscore the unique importance of this part of the Galaxy for constraining the ingredients of chemical evolution modeling and for improving our understanding of the evolution of the Galaxy as a whole.

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.

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.

The first detailed chemical abundance analysis of the M dwarf (M4.0) exoplanet-hosting star Ross 128 is presented here, based upon near-infrared (1.5--1.7 \micron) high-resolution ($R$$\sim$22,500) spectra from the SDSS-APOGEE survey. We determined precise atmospheric parameters $T_{\rm eff}$=3231$\pm$100K, log$g$=4.96$\pm$0.11 dex and chemical abundances of eight elements (C, O, Mg, Al, K, Ca, Ti, and Fe), finding Ross 128 to have near solar metallicity ([Fe/H] = +0.03$\pm$0.09 dex). The derived results were obtained via spectral synthesis (1-D LTE) adopting both MARCS and PHOENIX model atmospheres; stellar parameters and chemical abundances derived from the different adopted models do not show significant offsets. Mass-radius modeling of Ross 128b indicate that it lies below the pure rock composition curve, suggesting that it contains a mixture of rock and iron, with the relative amounts of each set by the ratio of Fe/Mg. If Ross 128b formed with a sub-solar Si/Mg ratio, and assuming the planet's composition matches that of the host-star, it likely has a larger core size relative to the Earth. The derived planetary parameters -- insolation flux (S$_{\rm Earth}$=1.79$\pm$0.26) and equilibrium temperature ($T_{\rm eq}$=294$\pm$10K) -- support previous findings that Ross 128b is a temperate exoplanet in the inner edge of the habitable zone.

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 report the first detailed chemical abundance analysis of the exoplanet-hosting M-dwarf stars Kepler-138 and Kepler-186 from the analysis of high-resolution ($R$ $\sim$ 22,500) $H$-band spectra from the SDSS IV - APOGEE survey. Chemical abundances of thirteen elements - C, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, and Fe - are extracted from the APOGEE spectra of these early M-dwarfs via spectrum syntheses computed with an improved line list that takes into account H$_{2}$O and FeH lines. This paper demonstrates that APOGEE spectra can be analyzed to determine detailed chemical compositions of M-dwarfs. Both exoplanet-hosting M-dwarfs display modest sub-solar metallicities: [Fe/H]$_{Kepler-138}$ = -0.09 $\pm$ 0.09 dex and [Fe/H]$_{Kepler-186}$ = -0.08 $\pm$ 0.10 dex. The measured metallicities resulting from this high-resolution analysis are found to be higher by $\sim$0.1-0.2 dex than previous estimates from lower-resolution spectra. The C/O ratios obtained for the two planet-hosting stars are near-solar, with values of 0.55 $\pm$ 0.10 for Kepler-138 and 0.52 $\pm$ 0.12 for Kepler-186. Kepler-186 exhibits a marginally enhanced [Si/Fe] ratio.

The Gaia-ESO Survey is obtaining high-quality spectroscopic data for about 10^5 stars using FLAMES at the VLT. UVES high-resolution spectra are being collected for about 5000 FGK-type stars. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO 2nd internal release and will be part of its 1st public release of advanced data products. The final parameter scale is tied to the one defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. A set of open and globular clusters is used to evaluate the physical soundness of the results. Each methodology is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted-medians of those from the individual methods. The recommended results successfully reproduce the benchmark stars atmospheric parameters and the expected Teff-log g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log g, and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50-100 K for Teff, 0.10-0.25 dex for log g, and 0.05-0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex.