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

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.

A large population of ultra-diffuse galaxies (UDGs) was recently discovered in the Coma cluster. Here we present optical spectra of three such UDGs, DF7, DF44 and DF17, which have central surface brightnesses of $\mu_g \approx 24.4-25.1$ mag arcsec$^{-2}$. The spectra were acquired as part of an ancillary program within the SDSS-IV MaNGA Survey. We stacked 19 fibers in the central regions from larger integral field units (IFUs) per source. With over 13.5 hours of on-source integration we achieved a mean signal-to-noise ratio (S/N) in the optical of $9.5$\AA$^{-1}$, $7.9$\AA$^{-1}$ and $5.0$\AA$^{-1}$, respectively, for DF7, DF44 and DF17. Stellar population models applied to these spectra enable measurements of recession velocities, ages and metallicities. The recession velocities of DF7, DF44 and DF17 are $6599^{+40}_{-25}$km/s, $6402^{+41}_{-39}$km/s and $8315^{+43}_{-43}$km/s, spectroscopically confirming that all of them reside in the Coma cluster. The stellar populations of these three galaxies are old and metal-poor, with ages of $7.9^{+3.6}_{-2.5}$Gyr, $8.9^{+4.3}_{-3.3}$Gyr and $9.1^{+3.9}_{-5.5}$Gyr, and iron abundances of $\mathrm{[Fe/H]}$ $-1.0^{+0.3}_{-0.4}$, $-1.3^{+0.4}_{-0.4}$ and $-0.8^{+0.5}_{-0.5}$, respectively. Their stellar masses are $3$-$6\times10^8 M_\odot$. The UDGs in our sample are as old or older than galaxies at similar stellar mass or velocity dispersion (only DF44 has an independently measured dispersion). They all follow the well-established stellar mass$-$stellar metallicity relation, while DF44 lies below the velocity dispersion-metallicity relation. These results, combined with the fact that UDGs are unusually large for their stellar mass, suggest that stellar mass plays a more important role in setting stellar population properties for these galaxies than either size or surface brightness.

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 measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

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.

Recent observations and a photometric search for variable stars in the Ursa Minor dwarf spheroidal galaxy (UMi dSph) are presented. Our observations were taken at Apache Point Observatory in 2014 and 2016 using the 0.5m ARCSAT telescope and the West Mountain Observatory 0.9m telescope of Brigham Young University in 2016. Previously known RR Lyrae stars in our field of view of the UMi dSph are identified, and we also catalog new variable star candidates. Tentative classifications are given for some of the new variable stars. We have conducted period searches with the data collected with the WMO telescope. Our ultimate goal is to create an updated catalog of variable stars in the UMi dSph and to compare the RR Lyrae stellar characteristics to other RR Lyrae stars found in the Local Group dSph galaxies.

Using kinematic maps from the Sloan Digital Sky Survey (SDSS) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we reveal that the majority of low-mass quenched galaxies exhibit coherent rotation in their stellar kinematics. Our sample includes all 39 quenched low-mass galaxies observed in the first year of MaNGA. The galaxies are selected with $M_{r} > -19.1$, stellar masses $10^{9}$ M$_{\odot} < M_{\star} < 5\times10^{9}$ M$_{\odot}$, EW$_{H\alpha} <2$ \AA, and all have red colours $(u-r)>1.9$. They lie on the size-magnitude and $\sigma$-luminosity relations for previously studied dwarf galaxies. Just six ($15\pm5.7$ per cent) are found to have rotation speeds $v_{e,rot} < 15$ km s$^{-1}$ at $\sim1$ $R_{e}$, and may be dominated by pressure support at all radii. Two galaxies in our sample have kinematically distinct cores in their stellar component, likely the result of accretion. Six contain ionised gas despite not hosting ongoing star formation, and this gas is typically kinematically misaligned from their stellar component. This is the first large-scale Integral Field Unit (IFU) study of low mass galaxies selected without bias against low-density environments. Nevertheless, we find the majority of these galaxies are within $\sim1.5$ Mpc of a bright neighbour ($M_{K} < -23$; or M$_{\star} > 5\times10^{10}$ M$_{\odot}$), supporting the hypothesis that galaxy-galaxy or galaxy-group interactions quench star formation in low-mass galaxies. The local bright galaxy density for our sample is $\rho_{proj} = 8.2\pm2.0$ Mpc$^{-2}$, compared to $\rho_{proj} = 2.1\pm0.4$ Mpc$^{-2}$ for a star forming comparison sample, confirming that the quenched low mass galaxies are preferentially found in higher density environments.

We quantify and correct systematic errors in PPMXL proper motions using extragalactic sources from the first two LAMOST data releases and the Veron-Cetty & Veron Catalog of Quasars. Although the majority of the sources are from the Veron catalog, LAMOST makes important contributions in regions that are not well-sampled by previous catalogs, particularly at low Galactic latitudes and in the south Galactic cap. We show that quasars in PPMXL have measureable and significant proper motions, which reflect the systematic zero-point offsets present in the catalog. We confirm the global proper motion shifts seen by Wu, Ma, & Zhou (2011), and additionally find smaller-scale fluctuations of the QSO-derived corrections to an absolute frame. We average the proper motions of 158,106 extragalactic objects in bins of 3x3 degrees and present a table of proper motion corrections.

We find that Galactic disk stars near the anticenter exhibit velocity asymmetries in both the Galactocentric radial and vertical components across the mid-plane as well as azimuthally. These findings are based on LAMOST spectroscopic velocities for a sample of ~400,000 F-type stars, combined with proper motions from the PPMXL catalog for which we have derived corrections to the zero points based in part on spectroscopically discovered galaxies and QSOs from LAMOST. In the region within 2 kpc outside the Sun's radius and +/-2 kpc from the Galactic midplane, we show that stars above the plane exhibit net outward radial motions with downward vertical velocities, while stars below the plane have roughly the opposite behavior. We discuss this in the context of other recent findings, and conclude that we are likely seeing the signature of vertical disturbances to the disk due to an external perturbation.