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.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We report on the diversity in quasar spectra from the Baryon Oscillation Spectroscopic Survey. After filtering the spectra to mitigate selection effects and Malmquist bias associated with a nearly flux-limited sample, we create high signal-to-noise ratio composite spectra from 58,656 quasars (2.1 \le z \le 3.5), binned by luminosity, spectral index, and redshift. With these composite spectra, we confirm the traditional Baldwin effect (BE, i.e., the anticorrelation of C IV equivalent width (EW) and luminosity) that follows the relation W_\lambda \propto L^{\beta_w} with slope \beta_w = -0.35 \pm 0.004, -0.35 \pm 0.005, and -0.41 \pm 0.005 for z = 2.25, 2.46, and 2.84, respectively. In addition to the redshift evolution in the slope of the BE, we find redshift evolution in average quasar spectral features at fixed luminosity. The spectroscopic signature of the redshift evolution is correlated at 98% with the signature of varying luminosity, indicating that they arise from the same physical mechanism. At a fixed luminosity, the average C IV FWHM decreases with increasing redshift and is anti-correlated with C IV EW. The spectroscopic signature associated with C IV FWHM suggests that the trends in luminosity and redshift are likely caused by a superposition of effects that are related to black hole mass and Eddington ratio. The redshift evolution is the consequence of a changing balance between these two quantities as quasars evolve toward a population with lower typical accretion rates at a given black hole mass.

The Baryon Oscillation Spectroscopic Survey (BOSS) has collected more than 150,000 $2.1 \leq z \leq 3.5$ quasar spectra since 2009. Using this unprecedented sample, we create a composite spectrum in the rest-frame of 102,150 quasar spectra from 800 \AA\ to 3300 \AA\ at a signal-to-noise ratio close to 1000 per pixel ($\Delta v$ of 69 km~s$^{-1}$). Included in this analysis is a correction to account for flux calibration residuals in the BOSS spectrophotometry. We determine the spectral index as a function of redshift of the full sample, warp the composite spectrum to match the median spectral index, and compare the resulting spectrum to SDSS photometry used in target selection. The quasar composite matches the color of the quasar population to within 0.02 magnitudes in $g-r$, 0.03 magnitudes in $r-i$, and 0.01 magnitudes in $i-z$ over the redshift range $2.2<z<2.6$. The composite spectrum deviates from the imaging photometry by 0.05 magnitudes around $z = 2.7$, likely due to differences in target selection as the quasar colors become similar to the stellar locus at this redshift. Finally, we characterize the line features in the high signal-to-noise composite and identify nine faint lines not found in the previous composite spectrum from SDSS.