The Lyman-$\alpha$ (hereafter Ly-$\alpha$) forest is a probe of large-scale matter density fluctuations at high redshift, $z > 2.1$. It consists of HI absorption spectra along individual lines-of-sight. If the line-of-sight density is large enough, 3D maps of HI absorption can be inferred by tomographic reconstruction. In this article, we investigate the Ly-$\alpha$ forest available in the Stripe 82 field ($220\,\mathrm{deg^{2}}$), based on the quasar spectra from SDSS Data Release DR16. The density of observed quasar spectra is $37\,\mathrm{deg^{-2}}$ with a mean pixel signal-to-noise ratio of two per angstrom. This study provides an intermediate case between the average SDSS density and that of the much denser but smaller CLAMATO survey. We derive a 3D map of large-scale matter fluctuations from these data, using a Wiener filter technique. The total volume of the map is $0.94\,\mathrm{h^{-3} Gpc^{3}}$. Its resolution is $13\,\mathrm{h^{-1} Mpc}$, which is related to the mean transverse distance between nearest lines-of-sight. From this map, we provide a catalog of voids and protocluster candidates in the cosmic web. The map-making and void catalog are compared to simulated eBOSS Stripe 82 observations. A stack over quasar positions provides a visualization of the Ly-$\alpha$ quasar cross-correlation. This tomographic reconstruction constitutes the largest-volume high-redshift 3D map of matter fluctuations.

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 statistical power of Lyman-${\alpha}$ forest Baryon Acoustic Oscillation (BAO) measurements is set to increase significantly in the coming years as new instruments such as the Dark Energy Spectroscopic Instrument deliver progressively more constraining data. Generating mock datasets for such measurements will be important for validating analysis pipelines and evaluating the effects of systematics. With such studies in mind, we present LyaCoLoRe: a package for producing synthetic Lyman-${\alpha}$ forest survey datasets for BAO analyses. LyaCoLoRe transforms initial Gaussian random field skewers into skewers of transmitted flux fraction via a number of fast approximations. In this work we explain the methods of producing mock datasets used in LyaCoLoRe, and then measure correlation functions on a suite of realisations of such data. We demonstrate that we are able to recover the correct BAO signal, as well as large-scale bias parameters similar to literature values. Finally, we briefly describe methods to add further astrophysical effects to our skewers - high column density systems and metal absorbers - which act as potential complications for BAO analyses.

We present a measurement of the 1D Ly$\alpha$ forest flux power spectrum, using the complete Baryon Oscillation Spectroscopic Survey (BOSS) and first extended-BOSS (eBOSS) quasars at $z_{\rm qso}>2.1$, corresponding to the fourteenth data release (DR14) of the Sloan Digital Sky Survey (SDSS). Our results cover thirteen bins in redshift from $z_{\rm Ly\alpha}=2.2$ to 4.6, and scales up to $k=0.02\rm \,(km/s)^{-1}$. From a parent sample of 180,413 visually inspected spectra, we selected the 43,751 quasars with the best quality; this data set improves the previous result from the ninth data release (DR9), both in statistical precision (achieving a reduction by a factor of two) and in redshift coverage. We also present a thorough investigation of identified sources of systematic uncertainties that affect the measurement. The resulting 1D power spectrum of this work is in excellent agreement with the one from the BOSS DR9 data.

We measure the imprint of primordial baryon acoustic oscillations (BAO) in the correlation function of Ly$\alpha$ absorption in quasar spectra from the Baryon Oscillation Spectroscopic Survey (BOSS) and the extended BOSS (eBOSS) in Data Release 14 (DR14) of the Sloan Digital Sky Survey (SDSS)-IV. In addition to 179,965 spectra with absorption in the Lyman-$\alpha$ (Ly$\alpha$) region, we use, for the first time, Ly$\alpha$ absorption in the Lyman-$\beta$ region of 56,154 spectra. We measure the Hubble distance, $D_H$, and the comoving angular diameter distance, $D_M$, relative to the sound horizon at the drag epoch $r_d$ at an effective redshift $z=2.34$. Using a physical model of the correlation function outside the BAO peak, we find $D_H(2.34)/r_d=8.86\pm 0.29$ and $D_M(2.34)/r_d=37.41\pm 1.86$, within 1$\sigma$ from the flat-$\Lambda$CDM model consistent with CMB anisotropy measurements. With the addition of polynomial "broadband" terms, the results remain within one standard deviation of the CMB-inspired model. Combined with the quasar-Ly$\alpha$ cross-correlation measurement presented in a companion paper Blomqvist19, the BAO measurements at $z=2.35$ are within 1.7$\sigma$ of the predictions of this model.

We present a measurement of the baryon acoustic oscillation (BAO) scale at redshift $z=2.35$ from the three-dimensional correlation of Lyman-$\alpha$ (Ly$\alpha$) forest absorption and quasars. The study uses 266,590 quasars in the redshift range $1.77<z<3.5$ from the Sloan Digital Sky Survey (SDSS) Data Release 14 (DR14). The sample includes the first two years of observations by the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS), providing new quasars and re-observations of BOSS quasars for improved statistical precision. Statistics are further improved by including Ly$\alpha$ absorption occurring in the Ly$\beta$ wavelength band of the spectra. From the measured BAO peak position along and across the line of sight, we determined the Hubble distance $D_{H}$ and the comoving angular diameter distance $D_{M}$ relative to the sound horizon at the drag epoch $r_{d}$: $D_{H}(z=2.35)/r_{d}=9.20\pm 0.36$ and $D_{M}(z=2.35)/r_{d}=36.3\pm 1.8$. These results are consistent at $1.5\sigma$ with the prediction of the best-fit spatially-flat cosmological model with the cosmological constant reported for the Planck (2016) analysis of cosmic microwave background anisotropies. Combined with the Ly$\alpha$ auto-correlation measurement presented in a companion paper, the BAO measurements at $z=2.34$ are within $1.7\sigma$ of the predictions of this model.

We present the first attempt at measuring the baryonic acoustic oscillations (BAO) in the large scale cross-correlation between the magnesium-II doublet (MgII) flux transmission field and the position of quasar and galaxy tracers. The MgII flux transmission continuous field at $0.3 < z < 1.6$ is measured from 500,589 quasar spectra obtained in the Baryonic Oscillation Spectroscopic Survey (BOSS) and the extended BOSS (eBOSS). The position of 246,697 quasar tracers and 1,346,776 galaxy tracers are extracted from the Sloan Digital Sky Survey (SDSS) I, II, BOSS and eBOSS catalogs. In addition to measuring the cosmological BAO scale and the biased matter density correlation, this study allows tests and improvements to cosmological Lyman-$\alpha$ analyses. A feature consistent with that of the BAO is detected at a significance of $\Delta \chi^{2} = 7.25$. The measured MgII linear transmission bias parameters are $b_{\mathrm{MgII(2796)}} (z = 0.59) = (-6.82 \pm 0.54) \, \times 10^{-4}$ and $b_{\mathrm{MgII(2804)}} (z = 0.59) = (-5.55 \pm 0.46) \, \times 10^{-4}$, and the MgI bias is $b_{\mathrm{MgI(2853)}} (z = 0.59) = (-1.48 \pm 0.24) \, \times 10^{-4}$. Their redshift evolution is characterized by the power-law index: $\gamma_{\mathrm{Mg}} = 3.36 \pm 0.46$. These measurements open a new window towards using BAO from flux transmission at $z < 1$ in the final eBOSS sample and in the upcoming sample from the Dark Energy Spectroscopic Instrument.

We measure the effective optical depth in the Lyman-alpha (Ly$\alpha$) forest using 40,035 quasar spectra from the Twelfth Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS) of SDSS-IV. A rigorous selection based on spectral index and the equivalent width of the C IV emission line is applied to choose seven uniform samples with minimal intrinsic variations across redshifts. Modeling the redshift evolution of the effective optical depth with a power-law, $\tau_\mathrm{eff}=\tau_0 (1 + z)^\gamma$ produces $\tau_0=(5.54 \pm 0.64)\times 10^{-3}$ and $\gamma=3.182 \pm 0.074$. The 2.2% precision estimate on $\gamma$ is dominated by systematic errors, likely arising from the bias and uncertainties in spectral index estimates. Even after incorporating the systematic errors, this work provides the most precise estimates of optical depth parameters to date. Finally, using the reconstructed Ly$\alpha$ forest continuum to directly measure the transmitted flux ratio as a function of redshift, we find deviations of less than 2.5% from the predictions from the global model and no convincing evidence for signal associated with He II reionization.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

Using 16,421 spectra from a sample of 340 quasars ($1.62<z<3.30$) from the SDSS Reverberation Mapping Project, we present an analysis of quasar spectral variability. We confirm the intrinsic Baldwin Effect and brighter-means-bluer trends in which emission line strength and color are associated with changes in luminosity. We then create a composite differential spectrum that captures the wavelength dependence of quasar variability as a function of change in luminosity. When using a bandpass around 1740 \AA\ to describe the luminosity, the differential spectrum follows a power law at longer wavelengths that breaks blueward of 1700 \AA. However, the shape of the differential spectrum, the location of the power law break, and the slope of the intrinsic Baldwin Effect all vary with the choice of bandpass used to define the change in quasar luminosity. We propose that the observed behavior can be explained by inhomogeneous accretion or slim accretion disk models where delays in the reprocessing of incident light on the accretion disk cause local deviations in temperature from the thin disk model. Finally, we quantify the effects on cosmology studies due to the variations in the quasar spectrum in the Lyman-$\alpha$ forest wavelength range. Using the observed spectroscopic signatures to predict the quasar continuum over the interval $1040 < \lambda < 1200$ \AA, we find that the derived spectral templates can reduce the uncertainty of the Lyman-$\alpha$ forest continuum level in individual epochs from 17.2% to 7.7%, near the level where systematic errors in SDSS flux calibration are expected to dominate.

We present measurements of the cross-correlation of the triply-ionized carbon (CIV) forest with quasars using Sloan Digital Sky Survey Data Release 14. The study exploits a large sample of new quasars from the first two years of observations by the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). The CIV forest is a weaker tracer of large-scale structure than the Ly$\alpha$ forest, but benefits from being accessible at redshifts $z<2$ where the quasar number density from eBOSS is high. Our data sample consists of 287,651 CIV forest quasars in the redshift range $1.4<z<3.5$ and 387,315 tracer quasars with $1.2<z<3.5$. We measure large-scale correlations from CIV absorption occuring in three distinct quasar rest-frame wavelength bands of the spectra referred to as the CIV forest, the SiIV forest and the Ly$\alpha$ forest. From the combined fit to the quasar-CIV cross-correlations for the CIV forest and the SiIV forest, the CIV redshift-space distortion parameter is $\beta_{\rm CIV}=0.27_{\ -0.14}^{\ +0.16}$ and its combination with the CIV linear transmission bias parameter is $b_{\rm CIV}(1+\beta_{\rm CIV})=-0.0183_{\ -0.0014}^{\ +0.0013}$ ($1\sigma$ statistical error) at the mean redshift $z=2.00$. Splitting the sample at $z=2.2$ to constrain the bias evolution with redshift yields the power-law exponent $\gamma=0.60\pm0.63$, indicating a significantly weaker redshift-evolution than for the Ly$\alpha$ forest linear transmission bias. We demonstrate that CIV absorption has the potential to be used as a probe of baryon acoustic oscillations (BAO). While the current data set is insufficient for a detection of the BAO peak feature, the final quasar samples for redshifts $1.4<z<2.2$ from eBOSS and the Dark Energy Spectroscopic Instrument (DESI) are expected to provide measurements of the isotropic BAO scale to $\sim7\%$ and $\sim3\%$ precision, respectively, at $z\simeq1.6$.

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.

We present the clustering measurements of quasars in configuration space based on the Data Release 14 (DR14) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey. This dataset includes 148,659 quasars spread over the redshift range $0.8\leq z \leq 2.2$ and spanning 2112.9 square degrees. We use the Convolution Lagrangian Perturbation Theory (CLPT) approach with a Gaussian Streaming (GS) model for the redshift space distortions of the correlation function and demonstrate its applicability for dark matter halos hosting eBOSS quasar tracers. At the effective redshift $z_{\rm eff} = 1.52$, we measure the linear growth rate of structure $f\sigma_{8}(z_{\rm eff})= 0.426 \pm 0.077$, the expansion rate $H(z_{\rm eff})= 159^{+12}_{-13}(r_{s}^{\rm fid}/r_s){\rm km.s}^{-1}.{\rm Mpc}^{-1}$, and the angular diameter distance $D_{A}(z_{\rm eff})=1850^{+90}_{-115}\,(r_s/r_{s}^{\rm fid}){\rm Mpc}$, where $r_{s}$ is the sound horizon at the end of the baryon drag epoch and $r_{s}^{\rm fid}$ is its value in the fiducial cosmology. The quoted errors include both systematic and statistical contributions. The results on the evolution of distances are consistent with the predictions of flat $\Lambda$-Cold Dark Matter ($\Lambda$-CDM) cosmology with Planck parameters, and the measurement of $f\sigma_{8}$ extends the validity of General Relativity (GR) to higher redshifts($z>1$) This paper is released with companion papers using the same sample. The results on the cosmological parameters of the studies are found to be in very good agreement, providing clear evidence of the complementarity and of the robustness of the first full-shape clustering measurements with the eBOSS DR14 quasar sample.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

We present a measurement of baryon acoustic oscillations (BAO) in the cross-correlation of quasars with the Ly$\alpha$-forest flux-transmission at a mean redshift $z=2.40$. The measurement uses the complete SDSS-III data sample: 168,889 forests and 234,367 quasars from the SDSS Data Release DR12. In addition to the statistical improvement on our previous study using DR11, we have implemented numerous improvements at the analysis level allowing a more accurate measurement of this cross-correlation. We also developed the first simulations of the cross-correlation allowing us to test different aspects of our data analysis and to search for potential systematic errors in the determination of the BAO peak position. We measure the two ratios $D_{H}(z=2.40)/r_{d} = 9.01 \pm 0.36$ and $D_{M}(z=2.40)/r_{d} = 35.7 \pm 1.7$, where the errors include marginalization over the non-linear velocity of quasars and the metal - quasar cross-correlation contribution, among other effects. These results are within $1.8\sigma$ of the prediction of the flat-$\Lambda$CDM model describing the observed CMB anisotropies. We combine this study with the Ly$\alpha$-forest auto-correlation function [2017A&A...603A..12B], yielding $D_{H}(z=2.40)/r_{d} = 8.94 \pm 0.22$ and $D_{M}(z=2.40)/r_{d} = 36.6 \pm 1.2$, within $2.3\sigma$ of the same flat-$\Lambda$CDM model.

We present a measurement of the DLA mean bias from the cross-correlation of DLA and the Ly$\alpha$ forest, updating earlier results of Font-Ribera et al. 2012 with the final BOSS Data Release and an improved method to address continuum fitting corrections. Our cross-correlation is well fitted by linear theory with the standard $\Lambda CDM$ model, with a DLA bias of $b_{\rm DLA} = 1.99\pm 0.11$; a more conservative analysis, which removes DLA in the Ly$\beta$ forest and uses only the cross-correlation at $r> 10{\rm h^{-1}\,Mpc}$, yields $b_{\rm DLA} = 2.00\pm 0.19$. This assumes the cosmological model from \cite{Planck2015} and the Ly$\alpha$ forest bias factors of Bautista et al. 2017, and includes only statistical errors obtained from bootstrap analysis. The main systematic errors arise from possible impurities and selection effects in the DLA catalogue, and from uncertainties in the determination of the Ly$\alpha$ forest bias factors and a correction for effects of high column density absorbers. We find no dependence of the DLA bias on column density or redshift. The measured bias value corresponds to a host halo mass $\sim 4\cdot10^{11} {\rm M_{\odot}}$ if all DLA were hosted in halos of a similar mass. In a realistic model where host halos over a broad mass range have a DLA cross section $\Sigma(M_h) \propto M_h^{\alpha}$ down to $M_h > M_{\rm min} =10^{8.5} {\rm M_{\odot}}$, we find that $\alpha > 1$ is required to have $b_{\rm DLA}> 1.7$, implying a steeper relation or higher value of $M_{\rm min}$ than is generally predicted in numerical simulations of galaxy formation.

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 present constraints on masses of active and sterile neutrinos. We use the one-dimensional Ly$\alpha$-forest power spectrum from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey (SDSS-III) and from the VLT/XSHOOTER legacy survey (XQ-100). In this paper, we present our own measurement of the power spectrum with the publicly released XQ-100 quasar spectra. Fitting Ly$\alpha$ data alone leads to cosmological parameters in excellent agreement with the values derived independently from Planck 2015 Cosmic Microwave Background (CMB) data. Combining BOSS and XQ-100 Ly$\alpha$ power spectra, we constrain the sum of neutrino masses to $\sum m_\nu < 0.8$ eV (95\% C.L). With the addition of CMB data, this bound is tightened to $\sum m_\nu < 0.14$ eV (95\% C.L.). With their sensitivity to small scales, Ly$\alpha$ data are ideal to constrain $\Lambda$WDM models. Using XQ-100 alone, we issue lower bounds on pure dark matter particles: $m_X \gtrsim 2.08 \: \rm{keV}$ (95\% C.L.) for early decoupled thermal relics, and $m_s \gtrsim 10.2 \: \rm{keV}$ (95\% C.L.) for non-resonantly produced right-handed neutrinos. Combining the 1D Ly$\alpha$ forest power spectrum measured by BOSS and XQ-100, we improve the two bounds to $m_X \gtrsim 4.17 \: \rm{keV}$ and $m_s \gtrsim 25.0 \: \rm{keV}$ (95\% C.L.). The $3~\sigma$ bound shows a more significant improvement, increasing from $m_X \gtrsim 2.74 \: \rm{keV}$ for BOSS alone to $m_X \gtrsim 3.10 \: \rm{keV}$ for the combined BOSS+XQ-100 data set. Finally, we include in our analysis the first two redshift bins ($z=4.2$ and $z=4.6$) of the power spectrum measured with the high-resolution HIRES/MIKE spectrographs. The addition of HIRES/MIKE power spectrum allows us to further improve the two limits to $m_X \gtrsim 4.65 \: \rm{keV}$ and $m_s \gtrsim 28.8 \: \rm{keV}$ (95\% C.L.).

In this study, we probe the transition to cosmic homogeneity in the Large Scale Structure (LSS) of the Universe using the CMASS galaxy sample of BOSS spectroscopic survey which covers the largest effective volume to date, $3\ h^{-3}\ \mathrm{Gpc}^3$ at $0.43 \leq z \leq 0.7$. We study the scaled counts-in-spheres, $\mathcal{N}(<r)$, and the fractal correlation dimension, $\mathcal{D}_2(r)$, to assess the homogeneity scale of the universe using a $Landy\ \&\ Szalay$ inspired estimator. Defining the scale of transition to homogeneity as the scale at which $\mathcal{D}_2(r)$ reaches 3 within $1\%$, i.e. $\mathcal{D}_2(r)>2.97$ for $r>\mathcal{R}_H$, we find $\mathcal{R}_H = (63.3\pm0.7) \ h^{-1}\ \mathrm{Mpc}$, in agreement at the percentage level with the predictions of the $\Lambda$CDM model $\mathcal{R}_H=62.0\ h^{-1}\ \mathrm{Mpc}$. Thanks to the large cosmic depth of the survey, we investigate the redshift evolution of the transition to homogeneity scale and find agreement with the $\Lambda$CDM prediction. Finally, we find that $\mathcal{D}_2$ is compatible with $3$ at scales larger than $300\ h^{-1}\ $Mpc in all redshift bins. These results consolidate the Cosmological Principle and represent a precise consistency test of the $\Lambda CDM$ model.

We study the first year of the eBOSS quasar sample in the redshift range $0.9<z<2.2$ which includes 68,772 homogeneously selected quasars. We show that the main source of systematics in the evaluation of the correlation function arises from inhomogeneities in the quasar target selection, particularly related to the extinction and depth of the imaging data used for targeting. We propose a weighting scheme that mitigates these systematics. We measure the quasar correlation function and provide the most accurate measurement to date of the quasar bias in this redshift range, $b_Q = 2.45 \pm 0.05$ at $\bar z=1.55$, together with its evolution with redshift. We use this information to determine the minimum mass of the halo hosting the quasars and the characteristic halo mass, which we find to be both independent of redshift within statistical error. Using a recently-measured quasar-luminosity-function we also determine the quasar duty cycle. The size of this first year sample is insufficient to detect any luminosity dependence to quasar clustering and this issue should be further studied with the final $\sim$500,000 eBOSS quasar sample.

We use flux-transmission correlations in \Lya forests to measure the imprint of baryon acoustic oscillations (BAO). The study uses spectra of 157,783 quasars in the redshift range $2.1\le z \le 3.5$ from the Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12). Besides the statistical improvements on our previous studies using SDSS DR9 and DR11, we have implemented numerous improvements in the analysis procedure, allowing us to construct a physical model of the correlation function and to investigate potential systematic errors in the determination of the BAO peak position. The Hubble distance, $\DHub=c/H(z)$, relative to the sound horizon is $\DHub(z=2.33)/r_d=9.07 \pm 0.31$. The best-determined combination of comoving angular-diameter distance, $\DM$, and the Hubble distance is found to be $\DHub^{0.7}\DM^{0.3}/r_d=13.94\pm0.35$. This value is $1.028\pm0.026$ times the prediction of the flat-\lcdm model consistent with the cosmic microwave background (CMB) anisotropy spectrum. The errors include marginalization over the effects of unidentified high-density absorption systems and fluctuations in ultraviolet ionizing radiation. Independently of the CMB measurements, the combination of our results and other BAO observations determine the open-\lcdm density parameters to be $\om=0.296 \pm 0.029$, $\ol=0.699 \pm 0.100$ and $\Omega_k = -0.002 \pm 0.119$.

The BOSS quasar sample is used to study cosmic homogeneity with a 3D survey in the redshift range $2.2<z<2.8$. We measure the count-in-sphere, $N(<\! r)$, i.e. the average number of objects around a given object, and its logarithmic derivative, the fractal correlation dimension, $D_2(r)$. For a homogeneous distribution $N(<\! r) \propto r^3$ and $D_2(r)=3$. Due to the uncertainty on tracer density evolution, 3D surveys can only probe homogeneity up to a redshift dependence, i.e. they probe so-called "spatial isotropy". Our data demonstrate spatial isotropy of the quasar distribution in the redshift range $2.2<z<2.8$ in a model-independent way, independent of any FLRW fiducial cosmology, resulting in $3-\langle D_2 \rangle < 1.7 \times 10^{-3}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc for the quasar distribution. If we assume that quasars do not have a bias much less than unity, this implies spatial isotropy of the matter distribution on large scales. Then, combining with the Copernican principle, we finally get homogeneity of the matter distribution on large scales. Alternatively, using a flat $\Lambda$CDM fiducial cosmology with CMB-derived parameters, and measuring the quasar bias relative to this $\Lambda$CDM model, our data provide a consistency check of the model, in terms of how homogeneous the Universe is on different scales. $D_2(r)$ is found to be compatible with our $\Lambda$CDM model on the whole $10<r<1200 \, h^{-1}$Mpc range. For the matter distribution we obtain $3-\langle D_2 \rangle < 5 \times 10^{-5}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc, consistent with homogeneity on large scales.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.