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

The fundamental metallicity relation (FMR) states that galaxies of the same stellar mass but larger star formation rate (SFR) tend to have smaller gas-phase metallicity (<Zg>). It is thought to be fundamental because it naturally arises from the stochastic feeding of star-formation from external metal-poor gas accretion, a process extremely elusive to observe but essential according the cosmological simulations of galaxy formation. In this letter, we show how the FMR emerges from the local anti-correlation between SFR surface density and Zg recently observed to exist in disk galaxies. We analytically derive the global FMR from the local law, and then show that both relations agree quantitatively when considering the star-forming galaxies of the MaNGA survey. Thus, understanding the FMR becomes equivalent to understanding the origin of the anti-correlation between SFR and metallicity followed by the set of star-forming regions of any typical galaxy. The correspondence between local and global laws is not specific of the FMR, so that a number of local relations should exist associated with known global relations.

We report the detection of a highly collimated linear emission-line structure in the spiral galaxy NGC\,232 through the use of integral field spectroscopy data from the All-weather MUse Supernova Integral field Nearby Galaxies (AMUSING) survey. This jet--like feature extends radially from the nucleus and is primarily detected in [oiii]$\lambda$5007 without clear evidence of an optical continuum counterpart. The length of the radial structure projected on sky reaches $\sim 3$ kpc, which makes NGC\,232 the second longest emission-line jet reported. The ionized gas presents extreme [Oiii]/H$\beta$ and [Nii]/H$\alpha$ line ratios, increasing along the jet-like structure. We discuss three possible scenarios to explain the observed structure: (i) direct ionization of in-falling material from the intergalactic medium by the AGN; (ii) photo-ionization by an un-detected optical counter-part of the radio jet and (iii) fast shocks ionization due to the lateral expansion of the radio jet across the ISM. Our analysis favors in-situ ionization.

We study a sample of 28 S0 galaxies extracted from the integral-field spectroscopic (IFS) survey CALIFA. We combine an accurate two-dimensional (2D) multi-component photometric decomposition with the IFS kinematic properties of their bulges to understand their formation scenario. Our final sample is representative of S0s with high stellar masses ($M_{star}/M_{sun} > 10^{10}$). They lay mainly on the red sequence and live in relatively isolated environments similar to that of the field and loose groups. We use our 2D photometric decomposition to define the size and photometric properties of the bulges, as well as their location within the galaxies. We perform mock spectroscopic simulations mimicking our observed galaxies to quantify the impact of the underlying disc on our bulge kinematic measurements ($\lambda$ and $v/\sigma$). We compare our bulge corrected kinematic measurements with the results from Schwarzschild dynamical modelling. The good agreement confirms the robustness of our results and allows us to use bulge reprojected values of $\lambda$ and $v/\sigma$. We find that the photometric ($n$ and $B/T$) and kinematic ($v/\sigma$ and $\lambda$) properties of our field S0 bulges are not correlated. We demonstrate that this morpho-kinematic decoupling is intrinsic to the bulges and it is not due to projection effects. We conclude that photometric diagnostics to separate different types of bulges (disc-like vs classical) might not be useful for S0 galaxies. The morpho-kinematics properties of S0 bulges derived in this paper suggest that they are mainly formed by dissipation processes happening at high redshift, but dedicated high-resolution simulations are necessary to better identify their origin.

We present a two-dimensional multi-component photometric decomposition of 404 galaxies from the CALIFA Data Release 3. 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 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(Mstar/Msun)>11), the galaxy population is dominated by galaxies modelled with a single Sersic or a bulge+disc with a bulge-to-total (B/T) luminosity ratio B/T>0.2. At intermediate masses (9.5<log(Mstar/Msun)<11), galaxies described with bulge+disc but B/T < 0.2 are preponderant, whereas, at the low mass end (log(Mstar/Msun)<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(Mstar/Msun)<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, Type II, and Type III 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.