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 present the first release of the MaNGA Stellar Library (MaStar), which is a large, well-calibrated, high-quality empirical library covering the wavelength range of 3,622-10,354A at a resolving power of R~1800. The spectra were obtained using the same instrument as used by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, by piggybacking on the SDSS-IV/APOGEE-2N observations. Compared to previous empirical libraries, the MaStar library will have a higher number of stars and a more comprehensive stellar-parameter coverage, especially of cool dwarfs, low-metallicity stars, and stars with different [alpha/Fe], achieved by a sophisticated target selection strategy that takes advantage of stellar-parameter catalogs from the literature. This empirical library will provide a new basis for stellar population synthesis, and is particularly well-suited for stellar-population analysis of MaNGA galaxies. The first version of the library contains 8646 high-quality per-visit spectra for 3321 unique stars. Compared to photometry, the relative flux calibration of the library is accurate to 3.9% in g-r, 2.7% in r-i, and 2.2% in i-z. The data are released as part of Sloan Digital Sky Survey Data Release 15. We expect the final release of the library to contain more than 10,000 stars.

We present an elemental-abundance analysis, in the near-ultraviolet (NUV) spectral range, for the bright carbon-enhanced metal-poor (CEMP) stars HD196944 (V = 8.40, [Fe/H] = -2.41) and HD201626 (V = 8.16, [Fe/H] = -1.51), based on data acquired with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope. Both of these stars belong to the sub-class CEMP-s, and exhibit clear over-abundances of heavy elements associated with production by the slow neutron-capture process. HD196944 has been well-studied in the optical region, but we are able to add abundance results for six species (Ge, Nb, Mo, Lu, Pt, and Au) that are only accessible in the NUV. In addition, we provide the first determination of its orbital period, P=1325 days. HD201626 has only a limited number of abundance results based on previous optical work -- here we add five new species from the NUV, including Pb. We compare these results with models of binary-system evolution and s-process element production in stars on the asymptotic giant branch, aiming to explain their origin and evolution. Our best-fitting models for HD 196944 (M1,i = 0.9Mo, M2,i = 0.86Mo, for [Fe/H]=-2.2), and HD 201626 (M1,i = 0.9Mo , M2,i = 0.76Mo , for [Fe/H]=-2.2; M1,i = 1.6Mo , M2,i = 0.59Mo, for [Fe/H]=-1.5) are consistent with the current accepted scenario for the formation of CEMP-s stars.