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

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.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), one of three core programs in the Sloan Digital Sky Survey-IV (SDSS-IV), is an integral-field spectroscopic (IFS) survey of roughly 10,000 nearby galaxies. It employs dithered observations using 17 hexagonal bundles of 2 arcsec fibers to obtain resolved spectroscopy over a wide wavelength range of 3,600-10,300A. To map the internal variations within each galaxy, we need to perform accurate {\it spectral surface photometry}, which is to calibrate the specific intensity at every spatial location sampled by each individual aperture element of the integral field unit. The calibration must correct only for the flux loss due to atmospheric throughput and the instrument response, but not for losses due to the finite geometry of the fiber aperture. This requires the use of standard star measurements to strictly separate these two flux loss factors (throughput versus geometry), a difficult challenge with standard single-fiber spectroscopy techniques due to various practical limitations. Therefore, we developed a technique for spectral surface photometry using multiple small fiber-bundles targeting standard stars simultaneously with galaxy observations. We discuss the principles of our approach and how they compare to previous efforts, and we demonstrate the precision and accuracy achieved. MaNGA's relative calibration between the wavelengths of H$\alpha$ and H$\beta$ has a root-mean-square (RMS) of 1.7%, while that between [NII] $\lambda$6583A and [OII] $\lambda$3727A has an RMS of 4.7%. Using extinction-corrected star formation rates and gas-phase metallicities as an illustration, this level of precision guarantees that flux calibration errors will be sub-dominant when estimating these quantities. The absolute calibration is better than 5% for more than 89% of MaNGA's wavelength range.

We describe the design, manufacture, and performance of bare-fiber integral field units (IFUs) for the SDSS-IV survey MaNGA (Mapping Nearby Galaxies at APO) on the the Sloan 2.5 m telescope at Apache Point Observatory (APO). MaNGA is a luminosity-selected integral-field spectroscopic survey of 10,000 local galaxies covering 360-1030 nm at R ~ 2200. The IFUs have hexagonal dense packing of fibers with packing regularity of 3 um (RMS), and throughput of 96+/-0.5% from 350 nm to 1 um in the lab. Their sizes range from 19 to 127 fibers (3-7 hexagonal layers) using Polymicro FBP 120:132:150 um core:clad:buffer fibers to reach a fill fraction of 56%. High throughput (and low focal-ratio degradation) is achieved by maintaining the fiber cladding and buffer intact, ensuring excellent surface polish, and applying a multi-layer AR coating of the input and output surfaces. In operations on-sky, the IFUs show only an additional 2.3% FRD-related variability in throughput despite repeated mechanical stressing during plate plugging (however other losses are present). The IFUs achieve on-sky throughput 5% above the single-fiber feeds used in SDSS-III/BOSS, attributable to equivalent performance compared to single fibers and additional gains from the AR coating. The manufacturing process is geared toward mass-production of high-multiplex systems. The low-stress process involves a precision ferrule with hexagonal inner shape designed to lead inserted fibers to settle in a dense hexagonal pattern. The ferrule inner diameter is tapered at progressively shallower angles toward its tip and the final 2 mm are straight and only a few um larger than necessary to hold the desired number of fibers. This process scales to accommodate other fiber sizes and to IFUs with substantially larger fiber count. (Abridged)

We use deep, longslit spectra of the nearby edge on galaxy IC 2531 to obtain gas kinematics out to 5 radial scale-lengths (40 kpc) and 4 vertical scale-heights (1.7 kpc). The large vertical range spanned by our data offers unique leverage to constrain three-dimensional models. The shape of the observed emission-line profiles offer insights to line-of-sight density distributions in the disk, and we discuss the possibility that we are seeing disk-flaring in the ionized gas. Finally, we begin to quantify measurements of line shape to allow model galaxies to be compared to data across all radii and heights simultaneously.

We present measurements of how multimode fiber focal-ratio degradation (FRD) and throughput vary with levels of fiber surface polish from 60 to 0.5 micron grit. Measurements used full-beam and laser injection methods at wavelengths between 0.4 and 0.8 microns on 17 meter lengths of Polymicro FBP 300 and 400 micron core fiber. Full-beam injection probed input focal-ratios between f/3 and f/13.5, while laser injection allowed us to isolate FRD at discrete injection angles up to 17 degrees (f/1.6 marginal ray). We find (1) FRD effects decrease as grit size decreases, with the largest gains in beam quality occurring at grit sizes above 5 microns; (2) total throughput increases as grit size decreases, reaching 90% at 790 nm with the finest polishing levels; (3) total throughput is higher at redder wavelengths for coarser polishing grit, indicating surface-scattering as the primary source of loss. We also quantify the angular dependence of FRD as a function of polishing level. Our results indicate that a commonly adopted micro-bending model for FRD is a poor descriptor of the observed phenomenon.

We describe the design, construction, and expected performance of two new fiber integral field units (IFUs) --- HexPak and GradPak --- for the WIYN 3.5m Telescope Nasmyth focus and Bench Spectrograph. These are the first IFUs to provide formatted fiber integral field spectroscopy with simultaneous sampling of varying angular scales. HexPak and GradPak are in a single cable with a dual-head design, permitting easy switching between the two different IFU heads on the telescope without changing the spectrograph feed: the two heads feed a variable-width double-slit. Each IFU head is comprised of a fixed arrangement of fibers with a range of fiber diameters. The layout and diameters of the fibers within each array are scientifically-driven for observations of galaxies: HexPak is designed to observe face-on spiral or spheroidal galaxies while GradPak is optimized for edge-on studies of galaxy disks. HexPak is a hexagonal array of 2.9 arcsec fibers subtending a 40.9 arcsec diameter, with a high-resolution circular core of 0.94 arcsec fibers subtending 6 arcsec diameter. GradPak is a 39 by 55 arcsec rectangular array with rows of fibers of increasing diameter from angular scales of 1.9 arcsec to 5.6 arcsec across the array. The variable pitch of these IFU heads allows for adequate sampling of light profile gradients while maintaining the photon limit at different scales.

We present near infra-red light curves of supernova (SN) 2011fe in M101, including 34 epochs in H band starting fourteen days before maximum brightness in the B-band. The light curve data were obtained with the WIYN High-Resolution Infrared Camera (WHIRC). When the data are calibrated using templates of other Type Ia SNe, we derive an apparent H-band magnitude at the epoch of B-band maximum of 10.85 \pm 0.04. This implies a distance modulus for M101 that ranges from 28.86 to 29.17 mag, depending on which absolute calibration for Type Ia SNe is used.