This paper documents the 16th data release (DR16) from the Sloan Digital Sky Surveys (SDSS), 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 and new data from the SPectroscopic IDentification of ERosita Survey 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 measured the gas abundance profiles in a sample of 122 face-on spiral galaxies observed by the CALIFA survey and included all spaxels whose line emission was consistent with star formation. This type of analysis allowed us to improve the statistics with respect to previous studies, and to properly estimate the oxygen distribution across the entire disc to a distance of up to 3−4 disc effective radii (r e ). We confirm the results obtained from classical H  region analysis. In addition to the general negative gradient, an outer flattening can be observed in the oxygen abundance radial profile. An inner drop is also found in some cases. There is a common abundance gradient between 0.5 and 2.0 r e of α O/H = − 0.075 dex/r e with a scatter of σ = 0.016 dex/r e when normalising the distances to the disc effective radius. By performing a set of Kolmogorov-Smirnov tests, we determined that this slope is independent of other galaxy properties, such as morphology, absolute magnitude, and the presence or absence of bars. In particular, barred galaxies do not seem to display shallower gradients, as predicted by numerical simulations. Interestingly, we find that most of the galaxies in the sample with reliable oxygen abundance values beyond ∼2 effective radii (57 galaxies) present a flattening of the abundance gradient in these outer regions. This flattening is not associated with any morphological feature, which suggests that it is a common property of disc galaxies. Finally, we detect a drop or truncation of the abundance in the inner regions of 27 galaxies in the sample; this is only visible for the most massive galaxies.

We present a two-dimensional multi-component photometric decomposition of 404 galaxies from the Calar Alto Legacy Integral Field Area Data Release 3 (CALIFA-DR3). They represent all possible galaxies with no clear signs of interaction and not strongly inclined in the final CALIFA data release. Galaxies are modelled in the g, r, and i Sloan Digital Sky Survey (SDSS) images including, when appropriate, a nuclear point source, bulge, bar, and an exponential or broken disc component. We use a human-supervised approach to determine the optimal number of structures to be included in the fit. The dataset, including the photometric parameters of the CALIFA sample, is released together with statistical errors and a visual analysis of the quality of each fit. The analysis of the photometric components reveals a clear segregation of the structural composition of galaxies with stellar mass. At high masses (log(M ⋆ /M ⊙ ) > 11), the galaxy population is dominated by galaxies modelled with a single Sérsic or a bulge+disc with a bulge-tototal (B/T ) luminosity ratio B/T > 0.2. At intermediate masses (9.5 < log(M ⋆ /M ⊙ ) < 11), galaxies described with bulge+disc but B/T < 0.2 are preponderant, whereas, at the low mass end (log(M ⋆ /M ⊙ ) < 9.5), the prevailing population is constituted by galaxies modelled with either pure discs or nuclear point sources+discs (i.e., no discernible bulge). We obtain that 57% of the volume corrected sample of disc galaxies in the CALIFA sample host a bar. This bar fraction shows a significant drop with increasing galaxy mass in the range 9.5 < log(M ⋆ /M ⊙ ) < 11.5. The analyses of the extended multi-component radial profile result in a volume-corrected distribution of 62%, 28%, and 10% for the so-called Type I (pure exponential), Type II (down-bending), and Type III (up-bending) disc profiles, respectively. These fractions are in discordance with previous findings. We argue that the different methodologies used to detect the breaks are the main cause for these differences.