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 analyze the anisotropic clustering of the Sloan Digital Sky Survey-IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy Data Release 14 (DR14) sample combined with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of galaxies in the redshift range 0.6$<z<$1.0, which consists of 80,118 galaxies from eBOSS and 46,439 galaxies from the BOSS-CMASS sample. The eBOSS-CMASS Luminous Red Galaxy sample has a sky coverage of 1,844 deg$^2$, with an effective volume of 0.9 Gpc$^3$. The analysis was made in configuration space using a Legendre multipole expansion. The Redshift Space Distortion signal is modeled as a combination of the Convolution Lagrangian Perturbation Model and the Gaussian Streaming Model. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f (z_{\rm eff})\sigma_8(z_{\rm eff})=0.454 \pm0.139 $, and the Alcock-Paczynski dilation scales which constraints the angular diameter distance $D_A(z_{eff})=1466.5 \pm 136.6 (r_s/r_s^{\rm fid})$ and $H(z_{\rm eff})=105.8 \pm 16 (r_s^{\rm fid}/r_s) \mathrm{km\,s^{-1}\,Mpc^{-1}}$, where $r_s$ is the sound horizon at the end of the baryon drag epoch and $r_s^{\rm fid}$ is its value in the fiducial cosmology at an effective redshift $z_{\rm eff}=0.72$. These results are in full agreement with the current $\Lambda$-Cold Dark Matter ($\Lambda$-CDM) cosmological model inferred from Planck measurements. This study is the first eBOSS LRG full-shape analysis i.e. including Redshift-Space Distortions (RSD) simultaneously with the Alcock-Paczynski (AP) effect and the Baryon Acoustic Oscillation (BAO) scale.