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 use 80922 galaxies in the Galaxy And Mass Assembly (GAMA) survey to measure the galaxy luminosity function (LF) in different environments over the redshift range 0.04<z<0.26. The depth and size of GAMA allows us to define samples split by colour and redshift to measure the dependence of the LF on environment, redshift and colour. We find that the LF varies smoothly with overdensity, consistent with previous results, with little environmental dependent evolution over the last 3 Gyrs. The modified GALFORM model predictions agree remarkably well with our LFs split by environment, particularly in the most overdense environments. The LFs predicted by the model for both blue and red galaxies are consistent with GAMA for the environments and luminosities at which such galaxies dominate. Discrepancies between the model and the data seen in the faint end of the LF suggest too many faint red galaxies are predicted, which is likely to be due to the over-quenching of satellite galaxies. The excess of bright blue galaxies predicted in underdense regions could be due to the implementation of AGN feedback not being sufficiently effective in the lower mass halos.

We use the Galaxy And Mass Assembly (GAMA) survey to measure the local Universe mass dependent merger fraction and merger rate using galaxy pairs and the CAS structural method, which identifies highly asymmetric merger candidate galaxies. Our goals are to determine which types of mergers produce highly asymmetrical galaxies, and to provide a new measurement of the local galaxy major merger rate. We examine galaxy pairs at stellar mass limits down to $M_{*} = 10^{8}M_{\odot}$ with mass ratios of $<$100:1 and line of sight velocity differences of $\Delta V<500$ km s$^{-1}$. We find a significant increase in mean asymmetries for projected separations less than the sum of the individual galaxy's Petrosian 90 radii. For systems in major merger pairs with mass ratios of $<$4:1 both galaxies in the pair show a strong increase in asymmetry, while in minor merger systems (with mass ratios of $>$4:1) the lower mass companion becomes highly asymmetric, while the larger galaxy is much less affected. The fraction of highly asymmetric paired galaxies which have a major merger companion is highest for the most massive galaxies and drops progressively with decreasing mass. We calculate that the mass dependent major merger fraction is fairly constant at $\sim1.3-2\%$ between $10^{9.5}<M_{*}<10^{11.5} M_{\odot}$, and increases to $\sim4\%$ at lower masses. When the observability time scales are taken into consideration, the major merger rate is found to approximately triple over the mass range we consider. The total co-moving volume major merger rate over the range $10^{8.0}<M_{*}<10^{11.5} M_{\odot}$ is $(1.2 \pm 0.5) \times 10^{-3}$ $h^{3}_{70}$ Mpc$^{-3}$ Gyr$^{-1}$.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.