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.

The angular distribution of galaxies encodes a wealth of information about large scale structure. Ultimately, the Sloan Digital Sky Survey (SDSS) will record the angular positions of order 10^8 galaxies in five bands, adding significantly to the cosmological constraints. This is the first in a series of papers analyzing a rectangular stripe 2.5x90 degrees from early SDSS data. We present the angular correlation function for galaxies in four separate magnitude bins on angular scales ranging from 0.003 degrees to 15 degrees. Much of the focus of this paper is on potential systematic effects. We show that the final galaxy catalog -- with the mask accounting for regions of poor seeing, reddening, bright stars, etc. -- is free from external and internal systematic effects for galaxies brighter than r* = 22. Our estimator of the angular correlation function includes the effects of the integral constraint and the mask. The full covariance matrix of errors in these estimates is derived using mock catalogs with further estimates using a number of other methods.

We present the first measurements of clustering in the Sloan Digital Sky Survey (SDSS) galaxy redshift survey. Our sample consists of 29,300 galaxies with redshifts 5,700 km/s < cz < 39,000 km/s, distributed in several long but narrow (2.5-5 degree) segments, covering 690 square degrees. For the full, flux-limited sample, the redshift-space correlation length is approximately 8 Mpc/h. The two-dimensional correlation function \xi(r_p,\pi) shows clear signatures of both the small-scale, ``fingers-of-God'' distortion caused by velocity dispersions in collapsed objects and the large-scale compression caused by coherent flows, though the latter cannot be measured with high precision in the present sample. The inferred real-space correlation function is well described by a power law, \xi(r)=(r/6.1+/-0.2 Mpc/h)^{-1.75+/-0.03}, for 0.1 Mpc/h < r < 16 Mpc/h. The galaxy pairwise velocity dispersion is \sigma_{12} ~ 600+/-100 km/s for projected separations 0.15 Mpc/h < r_p < 5 Mpc/h. When we divide the sample by color, the red galaxies exhibit a stronger and steeper real-space correlation function and a higher pairwise velocity dispersion than do the blue galaxies. The relative behavior of subsamples defined by high/low profile concentration or high/low surface brightness is qualitatively similar to that of the red/blue subsamples. Our most striking result is a clear measurement of scale-independent luminosity bias at r < 10 Mpc/h: subsamples with absolute magnitude ranges centered on M_*-1.5, M_*, and M_*+1.5 have real-space correlation functions that are parallel power laws of slope ~ -1.8 with correlation lengths of approximately 7.4 Mpc/h, 6.3 Mpc/h, and 4.7 Mpc/h, respectively.

We compute the angular power spectrum C_l from 1.5 million galaxies in early SDSS data on large angular scales, l<600. The data set covers about 160 square degrees, with a characteristic depth of order 1 Gpc/h in the faintest (21<r<22) of our four magnitude bins. Cosmological interpretations of these results are presented in a companion paper by Dodelson et al (2001). The data in all four magnitude bins are consistent with a simple flat ``concordance'' model with nonlinear evolution and linear bias factors of order unity. Nonlinear evolution is particularly evident for the brightest galaxies. A series of tests suggest that systematic errors related to seeing, reddening, etc., are negligible, which bodes well for the sixtyfold larger sample that the SDSS is currently collecting. Uncorrelated error bars and well-behaved window functions make our measurements a convenient starting point for cosmological model fitting.

We present measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

The Sloan Digital Sky Survey is one of the first multicolor photometric and spectroscopic surveys designed to measure the statistical properties of galaxies within the local Universe. In this Letter we present some of the initial results on the angular 2-point correlation function measured from the early SDSS galaxy data. The form of the correlation function, over the magnitude interval 18<r*<22, is shown to be consistent with results from existing wide-field, photographic-based surveys and narrower CCD galaxy surveys. On scales between 1 arcminute and 1 degree the correlation function is well described by a power-law with an exponent of ~ -0.7. The amplitude of the correlation function, within this angular interval, decreases with fainter magnitudes in good agreement with analyses from existing galaxy surveys. There is a characteristic break in the correlation function on scales of approximately 1-2 degrees. On small scales, < 1', the SDSS correlation function does not appear to be consistent with the power-law form fitted to the 1'< theta <0.5 deg data. With a data set that is less than 2% of the full SDSS survey area, we have obtained high precision measurements of the power-law angular correlation function on angular scales 1' < theta < 1 deg, which are robust to systematic uncertainties. Because of the limited area and the highly correlated nature of the error covariance matrix, these initial results do not yet provide a definitive characterization of departures from the power-law form at smaller and larger angles. In the near future, however, the area of the SDSS imaging survey will be sufficient to allow detailed analysis of the small and large scale regimes, measurements of higher-order correlations, and studies of angular clustering as a function of redshift and galaxy type.

We have created a variety of composite quasar spectra using a homogeneous data set of over 2200 spectra from the Sloan Digital Sky Survey (SDSS). The quasar sample spans a redshift range of 0.044 <= z <= 4.789, and an absolute r' magnitude range of -18.0 to -26.5. The median composite covers a rest wavelength range from 800 - 8555A, and reaches a peak signal-to-noise ratio of over 300 per 1A resolution element in the rest frame. We have identified over 80 emission line features in the spectrum. Emission line peak velocity shifts of the broad permitted and semi-forbidden lines are strongly correlated with ionization energy, as previously suggested, but we find that the narrow forbidden lines are also shifted by amounts which are strongly correlated with ionization energy. The magnitude of the forbidden line shifts is < 100 km/s, compared to shifts of up to 550 km/s for some of the permitted and semi-forbidden lines. At wavelengths longer than the Ly-a emission, the continuum of the geometric mean composite is well-fit by two power-laws, with a break at approximately 5000A. The frequency power law index, alpha_nu, is -0.44 from 1300 - 5000A, and -2.45 redward of about 5000A. The abrupt change in slope can be accounted for partly by host galaxy contamination at low redshift. Stellar absorption lines, including higher-order Balmer lines, seen in the composites suggest that young or intermediate age stars make a significant contribution to the light of the host galaxies. An electronic table of the median quasar template is available.

During commissioning observations, the Sloan Digital Sky Survey (SDSS) has produced one of the largest existing galaxy redshift samples selected from CCD images. Using 11,275 galaxies complete to r^* = 17.6 over 140 square degrees, we compute the luminosity function of galaxies in the r^* band over a range -23 < M < -16 (for h=1). The result is well-described by a Schechter function with parameters phi_* = 0.0146 +/- 0.0012 h^3 Mpc^{-3}, M_* = -20.83 +/- 0.03, and alpha = -1.20 +/- 0.03. The implied luminosity density in r^* is j = (2.6 +/- 0.3) x 10^8 h L_sun Mpc^{-3}. The surface brightness selection threshold has a negligible impact for M < -18. We measure the luminosity function in the u^*, g^*, i^*, and z^* bands as well; the slope at low luminosities ranges from alpha=-1.35 to alpha=-1.2. We measure the bivariate distribution of r^* luminosity with half-light surface brightness, intrinsic color, and morphology. High surface brightness, red, highly concentrated galaxies are on average more luminous than low surface brightness, blue, less concentrated galaxies. If we synthesize results for R-band or b_j-band using the Petrosian magnitudes with which the SDSS measures galaxy fluxes, we obtain luminosity densities 2.0 times that found by the Las Campanas Redshift Survey in R and 1.4 times that found by the Two-degree Field Galaxy Redshift Survey in b_j. We are able to reproduce the luminosity functions obtained by these surveys if we also mimic their isophotal limits for defining galaxy magnitudes, which are shallower and more redshift dependent than the Petrosian magnitudes used by the SDSS. (Abridged)

On 26 May 1999, one of the Sloan Digital Sky Survey (SDSS) fiber-fed spectrographs saw astronomical first light. This was followed by the first spectroscopic commissioning run during the dark period of June 1999. We present here the first hour of extra-galactic spectroscopy taken during these early commissioning stages: an observation of the Coma cluster of galaxies. Our data samples the Southern part of this cluster, out to a radius of 1.5degrees and thus fully covers the NGC 4839 group. We outline in this paper the main characteristics of the SDSS spectroscopic systems and provide redshifts and spectral classifications for 196 Coma galaxies, of which 45 redshifts are new. For the 151 galaxies in common with the literature, we find excellent agreement between our redshift determinations and the published values. As part of our analysis, we have investigated four different spectral classification algorithms: spectral line strengths, a principal component decomposition, a wavelet analysis and the fitting of spectral synthesis models to the data. We find that a significant fraction (25%) of our observed Coma galaxies show signs of recent star-formation activity and that the velocity dispersion of these active galaxies (emission-line and post-starburst galaxies) is 30% larger than the absorption-line galaxies. We also find no active galaxies within the central (projected) 200 h-1 Kpc of the cluster. The spatial distribution of our Coma active galaxies is consistent with that found at higher redshift for the CNOC1 cluster survey. Beyond the core region, the fraction of bright active galaxies appears to rise slowly out to the virial radius and are randomly distributed within the cluster with no apparent correlation with the potential merger of the NGC 4839 group. [ABRIDGED]

(abridged) We present measurements of galaxy-galaxy lensing from early commissioning imaging data from the Sloan Digital Sky Survey (SDSS). We measure a mean tangential shear around a stacked sample of foreground galaxies in three bandpasses out to angular radii of 600'', detecting the shear signal at very high statistical significance. The shear profile is well described by a power-law. A variety of rigorous tests demonstrate the reality of the gravitational lensing signal and confirm the uncertainty estimates. We interpret our results by modeling the mass distributions of the foreground galaxies as approximately isothermal spheres characterized by a velocity dispersion and a truncation radius. The velocity dispersion is constrained to be 150-190 km/s at 95% confidence (145-195 km/s including systematic uncertainties), consistent with previous determinations but with smaller error bars. Our detection of shear at large angular radii sets a 95% confidence lower limit $s>140^{\prime\prime}$, corresponding to a physical radius of $260h^{-1}$ kpc, implying that galaxy halos extend to very large radii. However, it is likely that this is being biased high by diffuse matter in the halos of groups and clusters. We also present a preliminary determination of the galaxy-mass correlation function finding a correlation length similar to the galaxy autocorrelation function and consistency with a low matter density universe with modest bias. The full SDSS will cover an area 44 times larger and provide spectroscopic redshifts for the foreground galaxies, making it possible to greatly improve the precision of these constraints, measure additional parameters such as halo shape, and measure the properties of dark matter halos separately for many different classes of galaxies.

We report the discovery of a second field methane brown dwarf from the commissioning data of the Sloan Digital Sky Survey (SDSS). The object, SDSS J134646.45-003150.4 (SDSS 1346-00), was selected because of its very red color and stellar appearance. Its spectrum between 0.8-2.5 mic is dominated by strong absorption bands of H_2O and CH_4 and closely mimics those of Gliese 229B and SDSS 162414.37+002915.6 (SDSS 1624+00), two other known methane brown dwarfs. SDSS 1346-00 is approximately 1.5 mag fainter than Gliese 229B, suggesting that it lies about 11 pc from the sun. The ratio of flux at 2.1 mic to that at 1.27 mic is larger for SDSS 1346-00 than for Gliese 229B and SDSS 1624+00, which suggests that SDSS 1346-00 has a slightly higher effective temperature than the others. Based on a search area of 130 sq. deg. and a detection limit of z* = 19.8, we estimate a space density of 0.05 pc^-3 for methane brown dwarfs with T_eff ~ 1000 K in the 40 pc^3 volume of our search. This estimate is based on small-sample statistics and should be treated with appropriate caution.

This is the second paper in a series aimed at finding high-redshift quasars from five-color (u'g'r'i'z') imaging data taken along the Celestial Equator by the Sloan Digital Sky Survey (SDSS) during its commissioning phase. In this paper, we present 22 high-redshift quasars (z>3.6) discovered from ~250 deg^2 of data in the spring Equatorial Stripe, plus photometry for two previously known high-redshift quasars in the same region of sky. Our success rate of identifying high-redshift quasars is 68%. Five of the newly discovered quasars have redshifts higher than 4.6 (z=4.62, 4.69, 4.70, 4.92 and 5.03). All the quasars have i* < 20.2 with absolute magnitude -28.8 < M_B < -26.1 (h=0.5, q_0=0.5). Several of the quasars show unusual emission and absorption features in their spectra, including an object at z=4.62 without detectable emission lines, and a Broad Absorption Line (BAL) quasar at z=4.92.

We report observations of a luminous unresolved object at redshift z = 4.62, with a featureless optical spectrum redward of the Lyman alpha forest region, discovered from Sloan Digital Sky Survey (SDSS) commissioning data. The redshift is determined by the onset of the Lyman-alpha forest at lambda ~ 6800 A, and a Lyman Limit System at lambda = 5120 A. A strong Ly_alpha absorption system with weak metal absorption lines at z=4.58 is also identified in the spectrum. The object has a continuum absolute magnitude of -26.6 at 1450 A in the rest-frame (h_0=0.5, q_0=0.5), and therefore cannot be an ordinary galaxy. It shows no radio emission (the 3-sigma upper limit of its flux at 6 cm is 60 micron Jy), indicating an radio-to-optical flux ratio at least as small as that of the radio-weakest known BL Lacs. It is also not linearly polarized in the observed I band to a 3-sigma upper limit of 4%. Therefore, it is either the most distant BL Lac object known to date, with very weak radio emission, or a new type of unbeamed quasar, whose broad emission-line region is very weak or absent.

This paper describes the discovery of seven dwarf objects of spectral type `L' (objects cooler than the latest M dwarfs) in commissioning imaging data taken by the Sloan Digital Sky Survey (SDSS). Low-resolution spectroscopy shows that these objects have spectral types from L0 to L8. Comparison of the SDSS and 2MASS photometry for several of these objects indicates the presence of significant opacity at optical wavelengths, perhaps due to atmospheric dust. This comparison also demonstrates the high astrometric accuracy (better than 1'' for these faint sources) of both surveys. The L dwarfs are shown to occupy a distinctive region of color-color space as measured in the SDSS filters, which should enable their identification in a straightforward way. This should lead eventually to a complete sample of many hundreds of these low mass objects, or about one per 15 square degrees to i'~20, in the complete SDSS data set.

We report the discovery of the coolest field dwarf yet known, selected as a stellar object with extremely red colors from commissioning imaging data of the Sloan Digital Sky Survey. Its spectrum from 0.8 to 2.5 microns is dominated by strong bands of H_2 O and CH_4. Its spectrum and colors over this range are very similar to those of Gliese 229B, the only other known example of a methane dwarf. It is roughly 1.2 mag fainter than Gliese 229B, implying that it lies at a distance of roughly 10 pc. Such a cool object must have a mass well below the hydrogen-burning limit of 0.08 solar masses, and therefore is a genuine brown dwarf, with a probable mass in the range 0.015-0.06 solar masses for an age range of 0.3-5 Gyr.

We present photometric and spectroscopic observations of 15 high-redshift quasars (z > 3.6) discovered from ~140 deg^2 of five-color (u'g'r'i'z') imaging data taken by the Sloan Digital Sky Survey (SDSS) during its commissioning phase. The quasars are selected by their distinctive colors in SDSS multicolor space. Four of the quasars have redshifts higher than 4.6 (z=4.63, 4.75, 4.90 and 5.00, the latter being the highest redshift quasar yet known). In addition, two previously known z > 4 objects were recovered from the data. The quasars all have i* < 20 and have luminosities comparable to that of 3C273. The spectra of the quasars have similar features (strong, broad emission lines and substantial absorption blueward of the Ly alpha emission line) seen in previously known high-redshift quasars. Although the photometric accuracy and image quality fail to meet the final survey requirements, our success rate for identifying high-redshift quasars (17 quasars from 27 candidates) is much higher than that of previous multicolor surveys. However, the numbers of high-redshift quasars found is in close accord with the number density inferred from previous surveys.