Draft version September 26, 2017 Preprint typeset using LATEX style emulateapj v. 01/23/15

THE THIRTEENTH DATA RELEASE OF THE SLOAN DIGITAL SKY SURVEY: FIRST SPECTROSCOPIC DATA FROM THE SDSS-IV SURVEY MAPPING NEARBY GALAXIES AT APACHE POINT OBSERVATORY Franco D. Albareti1,2,3, Carlos Allende Prieto4,5, Andres Almeida6, Friedrich Anders7, Scott Anderson8, Brett H. Andrews9, Alfonso Aragon-Salamanca´ 10, Maria Argudo-Fernandez´ 11,12, Eric Armengaud13, Eric Aubourg14, Vladimir Avila-Reese15, Carles Badenes9, Stephen Bailey16, Beatriz Barbuy17,18, Kat Barger19, Jorge Barrera-Ballesteros20, Curtis Bartosz8, Sarbani Basu21, Dominic Bates22, Giuseppina Battaglia4,5, Falk Baumgarten7,23, Julien Baur13, Julian Bautista24, Timothy C. Beers25, Francesco Belfiore26,27, Matthew Bershady28, Sara Bertran de Lis4,5, Jonathan C. Bird29, Dmitry Bizyaev30,31,32, Guillermo A. Blanc33,34,35, Michael Blanton36, Michael Blomqvist37, Adam S. Bolton38,24, J. Borissova39,40, Jo Bovy41,42, William Nielsen Brandt43,44,45, Jonathan Brinkmann30, Joel R. Brownstein24, Kevin Bundy46, Etienne Burtin13, Nicolas´ G. Busca14, Hugo Orlando Camacho Chavez47,18, M. Cano D´ıaz15, Michele Cappellari48, Ricardo Carrera4,5, Yanping Chen49, Brian Cherinka20, Edmond Cheung46, Cristina Chiappini7, Drew Chojnowski31, Chia-Hsun Chuang7, Haeun Chung50, Rafael Fernando Cirolini51,18, Nicolas Clerc52, Roger E. Cohen53, Julia M. Comerford54, Johan Comparat1,2, Janaina Correa do Nascimento55,18, Marie-Claude Cousinou56, Kevin Covey57, Jeffrey D. Crane33, Rupert Croft58, Katia Cunha59, Jeremy Darling54, James W. Davidson Jr.60, Kyle Dawson24, Luiz Da Costa18,59, Gabriele Da Silva Ilha51,18, Alice Deconto Machado51,18, Timothee´ Delubac61, Nathan De Lee62,29, Axel De la Macorra15, Sylvain De la Torre37, Aleksandar M. Diamond-Stanic28, John Donor19, Juan Jose Downes15,63, Niv Drory64, Cheng Du65,Helion´ Du Mas des Bourboux13, Tom Dwelly52, Garrett Ebelke60, Arthur Eigenbrot28, Daniel J. Eisenstein66, Yvonne P. Elsworth67,68, Eric Emsellem69,70, Michael Eracleous43,44, Stephanie Escoffier56, Michael L. Evans8, Jesus´ Falcon-Barroso´ 4,5, Xiaohui Fan71, Ginevra Favole1,2, Emma Fernandez-Alvar15, J. G. Fernandez-Trincado72, Diane Feuillet31, Scott W. Fleming73,74, Andreu Font-Ribera16,46, Gordon Freischlad30,31, Peter Frinchaboy19, Hai Fu75, Yang Gao (高扬)12, Rafael A. Garcia76, R. Garcia-Dias4,5, D. A. Garcia-Hernandez´ 4,5, Ana E. Garcia Perez´ 4,5, Patrick Gaulme30,31, Junqiang Ge77, Douglas Geisler53, Bruce Gillespie20, Hector Gil Marin78,79,Leo´ Girardi80,18, Daniel Goddard81, Yilen Gomez Maqueo Chew15, Violeta Gonzalez-Perez81, Kathleen Grabowski30,31, Paul Green66, Catherine J. Grier43,44, Thomas Grier19,82, Hong Guo12, Julien Guy83, Alex Hagen43,44, Matt Hall60, Paul Harding84, R. E. Harley57, Sten Hasselquist31, Suzanne Hawley8, Christian R. Hayes60, Fred Hearty43, Saskia Hekker85, Hector Hernandez Toledo15, Shirley Ho 58,16, David W. Hogg36, Kelly Holley-Bockelmann29, Jon A. Holtzman31, Parker H. Holzer24, Jian Hu (胡剑)64, Daniel Huber86,87,68, Timothy Alan Hutchinson24, Ho Seong Hwang50,Hector´ J. Ibarra-Medel15, Inese I. Ivans24, KeShawn Ivory19,88, Kurt Jaehnig28, Trey W. Jensen24, Jennifer A. Johnson89,90, Amy Jones91, Eric Jullo37, T. Kallinger92, Karen Kinemuchi30,31, David Kirkby93, Mark Klaene30, Jean-Paul Kneib61, Juna A. Kollmeier33, Ivan Lacerna94, Richard R. Lane94, Dustin Lang42, Pierre Laurent13, David R. Law73, Alexie Leauthaud46, Jean-Marc Le Goff13, Chen Li12, Cheng Li12,65, Niu Li65, Ran Li77, Fu-Heng Liang (梁赋珩)64, Yu Liang65, Marcos Lima47,18, Lihwai Lin (林俐暉)94, Lin Lin (林琳)12, Yen-Ting Lin (林彥廷)94, Chao Liu77, Dan Long30, Sara Lucatello80, Nicholas MacDonald8, Chelsea L. MacLeod66, J. Ted Mackereth96, Suvrath Mahadevan43, Marcio Antonio Geimba Maia18,59, Roberto Maiolino26,27, Steven R. Majewski60, Olena Malanushenko30,31, Viktor Malanushenko30,31,N´ıcolas Dullius Mallmann55,18, Arturo Manchado4,5,97, Claudia Maraston81, Rui Marques-Chaves4,5, Inma Martinez Valpuesta4,5, Karen L. Masters81, Savita Mathur98, Ian D. McGreer71, Andrea Merloni52, Michael R. Merrifield10, Szabolcs Meszaros´ 99, Andres Meza100, Andrea Miglio67, Ivan Minchev7, Karan Molaverdikhani54,101, Antonio D. Montero-Dorta24, Benoit Mosser102, Demitri Muna89,90, Adam Myers103, Preethi Nair104, Kirpal Nandra52, Melissa Ness101, Jeffrey A. Newman9, Robert C. Nichol81, David L. Nidever71, Christian Nitschelm11, Julia O’Connell19, Audrey Oravetz30,31, Daniel J. Oravetz30,31, Zachary Pace28, Nelson Padilla94, Nathalie Palanque-Delabrouille13, Kaike Pan30,31, John Parejko8, Isabelle Paris37, Changbom Park50, John A. Peacock105, Sebastien Peirani106,46, Marcos Pellejero-Ibanez4,5, Samantha Penny81, Will J. Percival81, Jeffrey W. Percival28, Ismael Perez-Fournon4,5, Patrick Petitjean106, Matthew Pieri37, Marc H. Pinsonneault89, Alice Pisani56,106, Francisco Prada1,2,107, Abhishek Prakash9, Natalie Price-Jones42, M. Jordan Raddick20, Mubdi Rahman20, Anand Raichoor13, arXiv:1608.02013v2 [astro-ph.GA] 25 Sep 2017 Sandro Barboza Rembold51,18, A. M. Reyna57, James Rich13, Hannah Richstein19, Jethro Ridl52, Rogemar A. Riffel51,18, Rogerio´ Riffel55,18, Hans-Walter Rix101, Annie C. Robin72, Constance M. Rockosi108, Sergio Rodr´ıguez-Torres1, Tha´ıse S. Rodrigues80,109,18, Natalie Roe16, A. Roman Lopes110, Carlos Roman-Z´ u´niga˜ 111, Ashley J. Ross90, Graziano Rossi112, John Ruan8, Rossana Ruggeri81, Jessie C. Runnoe43,44, Salvador Salazar-Albornoz52, Mara Salvato52, Sebastian F. Sanchez15, Ariel G. Sanchez52, Jose´ R. Sanchez-Gallego8, Bas´ılio Xavier Santiago55,18, Ricardo Schiavon96, Jaderson S. Schimoia55,18, Eddie Schlafly16,113, David J. Schlegel16, Donald P. Schneider43,44, Ralph Schonrich¨ 48, Mathias Schultheis114, Axel Schwope7, Hee-Jong Seo115, Aldo Serenelli116, Branimir Sesar101, Zhengyi Shao12, Matthew Shetrone117, Michael Shull54, Victor Silva Aguirre68, M. F. Skrutskie60, Anˇze Slosar118, Michael Smith28, Verne V. Smith119, Jennifer Sobeck60, Garrett Somers89,29, Diogo Souto59, David V. Stark46, Keivan G. Stassun29, Matthias Steinmetz7, Dennis Stello86,120, Thaisa Storchi Bergmann55,18, Michael A. Strauss121, Alina Streblyanska4,5, Guy S. Stringfellow54, Genaro Suarez111, Jing Sun19, Manuchehr Taghizadeh-Popp20, Baitian Tang53, Charling Tao65,56, Jamie Tayar89, Mita Tembe60, Daniel Thomas81, Jeremy Tinker36, Rita Tojeiro22, Christy Tremonti28, Nicholas Troup60, Jonathan R. Trump43,113, Eduardo Unda-Sanzana11, O. Valenzuela15, Remco Van den Bosch101, 2 SDSS Collaboration

Mariana Vargas-Magana˜ 122, Jose Alberto Vazquez118, Sandro Villanova53, M. Vivek24, Nicole Vogt31, David Wake123,124, Rene Walterbos31, Yuting Wang77,81, Enci Wang12, Benjamin Alan Weaver36, Anne-Marie Weijmans22, David H. Weinberg89,90, Kyle B. Westfall81, David G. Whelan125, Eric Wilcots28, Vivienne Wild22, Rob A. Williams96, John Wilson60, W. M. Wood-Vasey9, Dominika Wylezalek20, Ting Xiao (肖 婷)12, Renbin Yan126, Meng Yang65, Jason E. Ybarra111,127, Christophe Yeche13, Fang-Ting Yuan12, Nadia Zakamska20, Olga Zamora4,5, Gail Zasowski20, Kai Zhang126, Cheng Zhao65, Gong-Bo Zhao77,81, Zheng Zheng77, Zheng Zheng24, Zhi-Min Zhou77, Guangtun Zhu20,113, Joel C. Zinn89, Hu Zou77 Draft version September 26, 2017

Abstract The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2), Mapping Nearby Galaxies at APO (MaNGA), and the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). As well as its core program, eBOSS contains two major subprograms: the Time Domain Spectroscopic Survey (TDSS) and the SPectroscopic IDentification of ERosita Sources (SPIDERS). This paper describes the first data release from SDSS-IV, Data Release 13 (DR13). DR13 makes publicly available the first 1390 spatially resolved integral field unit obser- vations of nearby galaxies from MaNGA. It includes new observations from eBOSS, completing the Sloan Extended QUasar, Emission-line galaxy, Luminous red galaxy Survey (SEQUELS), which also targeted variability-selected objects and X-ray selected objects. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification, and new reductions of the SDSS-III APOGEE-1 data, improving stellar parameters for dwarf stars and cooler stars. DR13 provides more robust and precise photometric calibrations. Value-added target catalogs relevant for eBOSS, TDSS, and SPIDERS and an updated red-clump catalog for APOGEE are also available. This paper describes the location and format of the data and provides references to important technical papers. The SDSS website, www.sdss.org, provides links to the data, tutori- als, 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. Subject headings: Atlases — Catalogs — Surveys

1 Instituto de F´ısica Te´orica, (UAM/CSIC), Universidad Aut´onomade Madrid, Cantoblanco, E-28049 Madrid, Spain University, Fort Worth, TX 76129, USA 2 Campus of International Excellence UAM+CSIC, Canto- 20 Center for Astrophysical Sciences, Department of Physics blanco, E-28049 Madrid, Spain and Astronomy, Johns Hopkins University, 3400 North Charles 3 la Caixa-Severo Ochoa Scholar Street, Baltimore, MD 21218, USA 4 Instituto de Astrof´ısica de Canarias, E-38205 La Laguna, 21 Department of Astronomy, Yale University, 52 Hillhouse Tenerife, Spain Avenue, New Haven, CT 06511, USA 5 Departamento de Astrof´ısica, Universidad de La Laguna 22 School of Physics and Astronomy, University of St An- (ULL), E-38206 La Laguna, Tenerife, Spain drews, North Haugh, St Andrews KY16 9SS, UK 6 Instituto de Investigaci´on Multidisciplinario en Ciencia 23 Humboldt-Universitat zu Berlin, Institut f¨ur Physik, y Tecnolog´ıa, Universidad de La Serena, Benavente 980, La Newtonstrasse 15,D-12589, Berlin, Germany Serena, Chile 24 Department of Physics and Astronomy, University of Utah, 7 Leibniz-Institut fuer Astrophysik Potsdam (AIP), An der 115 S. 1400 E., Salt Lake City, UT 84112, USA Sternwarte 16, D-14482 Potsdam, Germany 25 Department of Physics and JINA Center for the Evolution 8 Department of Astronomy, University of Washington, Box of the Elements, University of Notre Dame, Notre Dame, IN 351580, Seattle, WA 98195, USA 46556, USA 9 PITT PACC, Department of Physics and Astronomy, 26 Cavendish Laboratory, University of Cambridge, 19 J. J. University of Pittsburgh, Pittsburgh, PA 15260, USA Thomson Avenue, Cambridge CB3 0HE, UK 10 School of Physics and Astronomy, University of Notting- 27 University of Cambridge, Kavli Institute for Cosmology, ham, University Park, Nottingham, NG7 2RD, UK Cambridge, CB3 0HE, UK 11 Unidad de Astronom´ıa, Universidad de Antofagasta, 28 Department of Astronomy, University of Wisconsin- Avenida Angamos 601, Antofagasta 1270300, Chile Madison, 475 N. Charter St., Madison WI 53706, USA 12 Shanghai Astronomical Observatory, Chinese Academy of 29 Department of Physics and Astronomy, Vanderbilt Univer- Science, 80 Nandan Road, Shanghai 200030, P. R. China sity, 6301 Stevenson Center, Nashville, TN, 37235, USA 13 IRFU, CEA, Centre d’Etudes Saclay, 91191 Gif-Sur-Yvette 30 Apache Point Observatory, P.O. Box 59, Sunspot, NM Cedex, France 88349, USA 14 APC, University of Paris Diderot, CNRS/IN2P3, 31 Department of Astronomy, New Mexico State University, CEA/IRFU, Observatoire de Paris, Sorbonne Paris Cite, Las Cruces, NM 88003, USA France 32 Sternberg Astronomical Institute, Moscow State Univer- 15 Instituto de Astronom´ıa,Universidad Nacional Aut´onoma sity, Moscow 119992, Russia de M´exico,Apartado Postal 70-264, M´exicoD.F., 04510 Mexico 33 Observatories of the Carnegie Institution for Science, 813 16 Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Santa Barbara St, Pasadena, CA, 91101, USA Berkeley, CA 94720, USA 34 Departamento de Astronom´ıa, Universidad de Chile, 17 Universidade de S˜aoPaulo, IAG, Rua do Mat˜ao1226, S˜ao Camino del Observatorio 1515, Las Condes, Santiago, Chile Paulo 05508-900, Brazil 35 Centro de Astrof´ısica y Tecnolog´ıas Afines (CATA), 18 Laborat´orio Interinstitucional de e-Astronomia - LIneA, Camino del Observatorio 1515, Las Condes, Santiago, Chile Rua Gal. Jos´eCristino 77, Rio de Janeiro, RJ - 20921-400, 36 Center for Cosmology and Particle Physics, Department of Brazil Physics, New York University, New York, NY 10003, USA 19 Department of Physics & Astronomy, Texas Christian 37 Aix Marseille Universit´e, CNRS, LAM (Laboratoire Data Release 13 3

d’Astrophysique de Marseille) UMR 7326, F-13388, Marseille, 66 Harvard-Smithsonian Center for Astrophysics, 60 Garden France St., MS 20, Cambridge, MA 02138, USA 38 National Optical Astronomy Observatory, 950 N Cherry 67 School of Physics and Astronomy, University of Birming- Ave, Tucson, AZ 85719, USA ham, Edgbaston, Birmingham B15 2TT, UK 39 Instituto de F´ısica y Astronom´ıa, Universidad de Val- 68 Stellar Astrophysics Centre, Department of Physics and para´ıso, Av. Gran Breta˜na1111, Playa, Ancha, Casilla 5030, Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Chile Aarhus C, Denmark 40 Millennium Institute of Astrophysics (MAS),Monse˜nor 69 European Southern Observatory, Karl-Schwarzschild-Str. Sotero Sanz 100, oficina 104, Providencia, Santiago, Chile 2, 85748 Garching, Germany 41 Department of Astronomy and Astrophysics, University of 70 Universit´e Lyon 1, Observatoire de Lyon, Centre de Toronto, 50 St. George Street, Toronto, ON, M55 3H4, Canada ´ 42 Recherche Astrophysique de Lyon and Ecole Normale Sup´erieure Dunlap Institute for Astronomy & Astrophysics, University de Lyon, 9 avenue Charles Andr´e,F-69230 Saint-Genis Laval, of Toronto, 50 St. George Street, Toronto, Ontario, M5S 3H4, France Canada 71 43 Steward Observatory, University of Arizona, 933 North Department of Astronomy and Astrophysics, Eberly Col- Cherry Avenue, Tucson, AZ 85721, USA lege of Science, The Pennsylvania State University, 525 Davey 72 Institut UTINAM, CNRS-UMR6213, OSU THETA, Laboratory, University Park, PA 16802, USA 44 Universit´e Bourgogne-Franche-Comt´e, 41bis avenue de Institute for Gravitation and the Cosmos, Pennsylvania l’Observatoire, 25010 Besan¸conCedex, France State University, University Park, PA 16802, USA 73 45 Space Telescope Science Institute, 3700 San Martin Drive, Department of Physics, The Pennsylvania State University, Baltimore, MD 21218, USA University Park, PA 16802, USA 74 46 CSRA, Inc., 3700 San Martin Drive, Baltimore, MD 21218, Kavli IPMU (WPI), UTIAS, The University of Tokyo, USA Kashiwa, Chiba 277-8583, Japan 75 47 Department of Physics & Astronomy, University of Iowa, Departamento de F´ısica Matem´atica,Instituto de F´ısica, Iowa City, IA 52245 Universidade de S˜aoPaulo, CP 66318, CEP 05314-970, S˜ao 76 Laboratoire AIM, CEA/DRF – CNRS - Univ. Paris Paulo, SP, Brazil 48 Diderot – IRFU/SAp, Centre de Saclay, 91191 Gif-sur-Yvette Sub-department of Astrophysics, Department of Physics, Cedex, France University of Oxford, Denys Wilkinson Building, Keble Road, 77 National Astronomy Observatories, Chinese Academy of Oxford OX1 3RH, UK 49 Science, 20A Datun Road, Chaoyang District, Beijing, 100012, New York University Abu Dhabi, P.O. Box 129188, Abu P. R. China Dhabi, United Arab Emirates 78 50 Sorbonne Universit´es,UPMC Univ Paris 06, UMR 7095, Korea Institute for Advanced Study, 85 Hoegiro, Institut d’Astrophysique de Paris, 98 bis boulevard Arago, Dongdaemun-gu, Seoul 02455, Republic of Korea 51 F-75014, Paris, France Departamento de F´ısica, Centro de Ciˆencias Naturais e 79 Laboratoire de Physique Nucl´eaireet de Hautes Energies, Exatas, Universidade Federal de Santa Maria, 97105-900, Santa Universit´ePierre et Marie Curie, 4 Place Jussieu,75005 Paris, Maria, RS, Brazil 52 France Max-Planck-Institut f¨urextraterrestrische Physik, Postfach 80 Osservatorio Astronomico di Padova – INAF, Vicolo 1312, Giessenbachstraße, 85741 Garching, Germany. 53 dell’Osservatorio 5, I-35122, Padova, Italy Departamento de Astronom´ıa,Universidad de Concepci´on, 81 Institute of Cosmology & Gravitation, University of Casilla 160-C, Concepci´on, Chile 54 Portsmouth, Dennis Sciama Building, Portsmouth, PO1 3FX, Center for Astrophysics and Space Astronomy, Department UK of Astrophysical and Planetary Sciences, University of Colorado, 82 DePauw University, Greencastle, IN 46135, USA 389 UCB, Boulder, CO 80309, USA 83 55 LPNHE, CNRS/IN2P3, Universit´ePierre et Marie Curie Instituto de F´ısica,Universidade Federal do Rio Grande do Paris 6, Universit´eDenis Diderot Paris 7, 4 place Jussieu, 75252 Sul, Campus do Vale, Porto Alegre, RS, , 91501-907, Brazil 56 Paris CEDEX, France Aix Marseille Universit´e, CNRS/IN2P3, Centre de 84 Department of Astronomy, Case Western Reserve Univer- Physique des Particules de Marseille, UMR 7346, 13288, Mar- sity, Cleveland, OH 44106, USA seille, France 85 57 Max Planck Institute for Solar System Research,Justus- Department of Physics and Astronomy, Western Washing- von-Liebig-Weg 3, 37077 Goettingen, Germany ton University, 516 High Street, Bellingham, WA 98225, USA 86 58 Sydney Institute for Astronomy (SIfA), School of Physics, Department of Physics and McWilliams Center for Cos- University of Sydney, NSW 2006, Australia mology, Carnegie Mellon University, 5000 Forbes Avenue, 87 SETI Institute, 189 Bernardo Avenue, Mountain View, CA Pittsburgh, PA 15213, USA 59 94043, USA Observat´orio Nacional, 77 Rua General Jos´eCristino, Rio 88 Rice University, 6100 Main St, Houston, TX 77005, USA de Janeiro, 20921-400, Brazil 89 60 Department of Astronomy, The Ohio State University, 140 Department of Astronomy, University of Virginia, Char- W. 18th Ave., Columbus, OH 43210, USA lottesville, VA 22904-4325, USA 90 Center for Cosmology and AstroParticle Physics, The Ohio 61 Laboratoire d’Astrophysique, Ecole´ Polytechnique F´ed´erale State University, 191 W. Woodruff Ave., Columbus, OH 43210, de Lausanne, 1015 Lausanne, Switzerland USA 62 Department of Physics, Geology, and Engineering Tech, 91 Max-Planck-Institut fuer Astrophysik, Karl-Schwarzschild- Northern Kentucky University, Highland Heights, KY 41099, Str. 1, D-85748 Garching, Germany USA 92 Institute for Astronomy, University of Vienna, T¨urken- 63 Centro de Investigaciones de Astronom´ıa,AP 264, M´erida schanzstrasse 17, 1180 Vienna, Austria 5101-A, Venezuela 93 Department of Physics and Astronomy, University of 64 McDonald Observatory, The University of Texas at Austin, California, Irvine, Irvine, CA 92697, USA 1 University Station, Austin, TX 78712, USA 94 Instituto de Astrof´ısica,Pontificia Universidad Cat´olicade 65 Tsinghua Center for Astrophysics and Department of Chile, Av. Vicuna Mackenna 4860, 782-0436 Macul, Santiago, Physics, Tsinghua University, Beijing 100084, P. R. China Chile 4 SDSS Collaboration

1. INTRODUCTION drift-scan camera (Gunn et al. 1998) and two fiber-fed The Sloan Digital Sky Survey (SDSS) has been ob- optical R∼1800 spectrographs (Smee et al. 2013), each serving the Universe using the 2.5-meter Sloan Foun- with 320 fibers. dation Telescope (Gunn et al. 2006) at Apache Point The imaging and spectroscopy goals were not entirely Observatory (APO) for over 15 years. The goal of the fulfilled in the initial five-year period, and thus SDSS- original survey (2000–2005; York et al. 2000) was to map I was followed by SDSS-II (2005–2008; Abazajian et al. large-scale structure with five-band imaging over ∼ π 2009). Its first goal was to complete the planned initial steradians of the sky and spectra of ∼ 106 galaxies and large scale structure redshift survey as the Legacy pro- ∼105 quasars. This program was accomplished using a gram. It added SEGUE (Sloan Extension for Galactic Understanding and Exploration; Yanny et al. 2009), a spectroscopic survey focused on stars, and imaged an 95 Academia Sinica Institute of Astronomy and Astrophysics, P.O. Box 23-141, Taipei 10617, Taiwan average of once every five days a ∼ 200 sq. deg area 96 Astrophysics Research Institute, Liverpool John Moores along the celestial equator with repeated scans in SDSS- University, IC2, Liverpool Science Park, 146 Brownlow Hill, I (“Stripe 82”), to search for Type Ia supernovae and Liverpool L3 5RF, UK other transients (Frieman et al. 2008). 97 CSIC, Serrano, 117 - 28006, Madrid, Spain 98 Space Science Institute, 4750 Walnut Street, Suite 205, The success of SDSS as a cosmological probe, partic- Boulder, CO 80301, USA ularly the detection of the clustering of luminous red 99 ELTE Gothard Astrophysical Observatory, H-9704 Szom- galaxies (LRG) on the 100 h−1 Mpc scale expected bathely, Szent Imre herceg st. 112, Hungary from baryon acoustic oscillations (BAO; Eisenstein et al. 100 Departamento de Ciencias Fisicas, Universidad Andres Bello, Sazie 2212, Santiago, Chile 2005), led to the conception and implementation of BOSS 101 Max-Planck-Institut f¨ur Astronomie, Konigstuhl 17, (Baryon Oscillation Spectroscopic Survey; Dawson et al. D-69117 Heidelberg, Germany 2013) as the flagship program in the third version of the 102 LESIA, UMR 8109, Universit´e Pierre et Marie Curie, Universit´e Denis Diderot, Observatoire de Paris, F-92195 survey, SDSS-III (2008–2014; Eisenstein et al. 2011). As Meudon Cedex, France part of BOSS, SDSS-III imaged additional areas in the 103 Department of Physics and Astronomy, University of part of the south Galactic cap visible from the Northern Wyoming, Laramie, WY 82071, USA hemisphere. At the conclusion of these observations, the 104 Department of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487-0324, USA SDSS imaging camera was retired and is now part of the 105 Institute for Astronomy, University of Edinburgh, Royal permanent collection of the Smithsonian National Air Observatory, Edinburgh EH9 3HJ, UK and Space Museum1. During the summer shutdown in 106 Universit´eParis 6 et CNRS, Institut d’Astrophysique de 2009, the original SDSS spectrographs were replaced by Paris, 98bis blvd. Arago, 75014 Paris, France 107 Instituto de Astrof´ısicade Andaluc´ıa(CSIC), Glorieta de new, more efficient, spectrographs to be used by BOSS. la Astronom´ıa,E-18080 Granada, Spain The BOSS spectrographs featured expanded wavelength 108 Department of Astrnonomy and Astrophysics, University coverage (3560A˚< λ <10400 A),˚ new CCD detectors of California, Santa Cruz and UC Observatories, Santa Cruz, CA, 95064, USA with improved read noise, smaller pixels (15µm), and 109 Dipartimento di Fisica e Astronomia, Universit`a di improved quantum efficiency, and VPH gratings instead Padova, Vicolo dell’Osservatorio 2, I-35122 Padova, Italy of the original replicated surface relief gratings (Smee 110 Departamento de Fisica y Astronomia, Universidad de La et al. 2013). The two spectrographs were now fed by 500 Serena, Cisternas 1200, La Serena, 0000-0002-1379-4204, Chile 111 Instituto de Astronom´ıa,Universidad Nacional Aut´onoma fibers each so that the desired number of redshifts could de M´exico, Unidad Acad´emica en Ensenada, Ensenada BC be reached in the planned survey lifetime. 22860, Mexico During the first year of SDSS-III (2008–2009), the 112 Department of Astronomy and Space Science, Sejong SEGUE-2 survey (Rockosi et al., in preparation) used the University, Seoul 143-747, Korea 113 Hubble Fellow original SDSS spectrographs to observe additional Milky 114 Observatoire de la Cˆoted’Azur, Laboratoire Lagrange, Way halo fields to target distant halo samplers and trace 06304 Nice Cedex 4, France substructure. In SDSS-III all bright time could be used 115 Department of Physics and Astronomy, Ohio University, for scientific observations with the arrival of two new Clippinger Labs, Athens, OH 45701, USA 116 Institute of Space Sciences (IEEC-CSIC), Carrer de Can instruments. MARVELS (Multi-object APO Radial Magrans, E-08193, Barcelona, Spain Velocity Exoplanet Large-area Survey; Paegert et al. 117 University of Texas at Austin, McDonald Observatory, 2015; Thomas et al. 2016) used a novel multiplexing in- McDonald Observatory, TX, 79734, USA 118 Brookhaven National Laboratory, Upton, NY 11973, USA terferometer to observe 60 stars simultaneously to search 119 National Optical Astronomy Observatories, Tucson, AZ, for radial velocity variations caused by hot Jupiters 85719, USA and close brown dwarf companions. APOGEE (Apache 120 School of Physics, The University of New South Wales, Point Observatory Galactic Evolution Experiment; Ma- Sydney NSW 2052, Australia 121 Department of Astrophysical Sciences, Princeton Univer- jewski et al. 2017) used a 300-fiber, R∼22,000 H-band sity, Princeton, NJ 08544, USA spectrograph (Majewski et al. 2017; Wilson et al., in 122 Instituto de F´ısica, Universidad Nacional Aut´onomade preparation) to measure stellar parameters, chemical M´exico,Apdo. Postal 20-364, Mexico abundances, and radial velocities, mainly for red giants 123 Department of Physical Sciences, The Open University, Milton Keynes MK7 6AA, UK (Zasowski et al. 2013). 124 Department of Physics, University of North Carolina Since routine operations started in 2000, there have Asheville, One University Heights, Asheville, NC 28804, USA been thirteen public data releases. All data releases are 125 Physics Department, Austin College, Sherman, TX 75092, cumulative, re-releasing the best reduction of all previ- USA 126 Department of Physics and Astronomy, University of ously taken data. The most recent of these was Data Kentucky, 505 Rose Street, Lexington, KY 40506, USA 127 Department of Physics, Bridgewater College, 402 E 1 https://airandspace.si.edu/collection-objects/ College St., Bridgewater, VA 22812, USA camera-imaging-digital-sloan-digital-sky-survey-ccd-array Data Release 13 5

Release 12 (Alam et al. 2015), which contained all of 1% at z = 0.7 (LRG), 2% at z = 0.85 (ELG), and 2% the SDSS-III data, as well as the re-reduced data from at z = 1.5 (quasar) while the enhanced Lyman-α for- SDSS-I and SDSS-II. SDSS-I to SDSS-III imaged 14,555 est sample will improve BOSS constraints by a factor degree2 in the five filters (Fukugita et al. 1996; Doi of 1.4. Furthermore, the clustering from eBOSS trac- et al. 2010). Most of the sky was surveyed once or ers will allow new measurements of redshift-space dis- twice, but regions in Stripe 82 were observed between tortions (RSD), non-Gaussianity in the primordial den- 70 and 90 times. By the time of their retirement, the sity field, and the summed mass of neutrino species. SDSS spectrographs had obtained R∼1800 optical spec- Extensively observing these redshift ranges for the first tra for 860,836 galaxies, 116,003 quasars, and 521,990 time in SDSS required re-evaluation of targeting strate- stars. With the BOSS spectrographs, the survey has gies. Preliminary targeting schemes for many of these added data with similar resolution for 1,372,737 galax- classes of objects were tested as part of SEQUELS (Sloan ies, 294,512 quasars, and 247,216 stars. APOGEE has Extended QUasar, Emission-line galaxy, Luminous red contributed high-resolution IR spectra of 156,593 stars. galaxy Survey), which used 126 plates observed across MARVELS had observed 3233 stars with at least 16 SDSS-III and -IV. DR13 includes all SEQUELS data, epochs of radial velocity measurements. giving the largest SDSS sample to date of spectra target- The success of the previous Sloan Digital Sky Surveys ing intermediate redshift ranges. SDSS-IV also allocated and the continuing importance of the wide-field, multi- a significant number of fibers on the eBOSS plates to plexing capability of the Sloan Foundation Telescope mo- two additional dark-time programs. TDSS(Time Domain tivated the organization of the fourth phase of the survey, Spectroscopic Survey; Morganson et al. 2015) seeks to SDSS-IV (Blanton et al. 2017). SDSS-IV extends SDSS understand the nature of celestial variables by delib- observations to many fibers covering the spatial extent of erately targeting objects that vary in combined SDSS nearby galaxies, to new redshift regimes, and to the parts DR9 and Pan-STARRS1 data (PS1; Kaiser et al. 2002) of the Milky Way and dwarf galaxies that are only visible . A large number of the likely quasar targets so se- from the Southern Hemisphere. The MaNGA (Mapping lected are also targeted by the main eBOSS algorithms Nearby Galaxies at APO) survey studies galaxy forma- and therefore meet the goals of both surveys. TDSS- tion and evolution across a wide range of masses and only targets fill ∼ 10 spectra per square degree. The morphological types by observing a substantial portion main goal of the SPIDERS (Spectroscopic Identification of the optical spatial extent of ∼104 galaxies (Bundy et of eROSITA Sources) survey is to characterize a sub- al. 2015). It accomplishes this goal by employing 17 bun- set of X-ray sources identified by eROSITA (extended dles ranging in size between 19 and 127 fibers to cover Roentgen Survey with an Imaging Telescope Array; targets selected from an extended version of the NASA- Predehl et al. 2014) . Until the first catalog of eROSITA Sloan Atlas2 and 12 bundles of 7 fibers for calibration sources is available, SPIDERS will target sources from stars. These integral field units (IFUs) feed the BOSS the RASS (Roetgen All Sky Survey; Voges et al. 1999) spectrographs, providing information on the properties and XMM (X-ray Multi-mirror Mission; Jansen et al. of gas and stars in galaxies out to 1.5–2.5 effective radii 2001). SPIDERS will also obtain on average ∼ 10 spec- (Re). tra per square degree over the course of SDSS-IV, but the Another survey, eBOSS (extended Baryon Oscillation number of fibers per square degree on a plate is weighted Spectroscopic Survey Dawson et al. 2016), shares the toward the later years to take advantage of the new data dark time equally with MaNGA. eBOSS will measure from eROSITA. with percent-level precision the distance-redshift relation In bright time at APO, APOGEE-2, the successor to with BAO in the clustering of matter over the relatively APOGEE (hereafter APOGEE-1) in SDSS-IV, will con- unconstrained redshift range 0.6 < z < 2.2. This red- tinue its survey of the Milky Way stellar populations. shift range probes the Universe during its transition from Critical areas of the Galaxy, however, cannot be observed matter-dominated to dark-energy-dominated. Multiple from APO, including the more distant half of the Galac- measurements of the angular diameter distance (dA(z)) tic bar, the fourth quadrant of the disk, and important and Hubble parameter (H(z)) from BAO over the red- dwarf satellites of the Milky Way, such as the Magellanic shifts covered by eBOSS are therefore crucial for under- Clouds and some dwarf spheroidals. SDSS-IV will for the standing the nature of dark energy. eBOSS will use first time include operations outside of APO as the result spectroscopic redshifts from more than 400,000 LRGs of Carnegie Observatories and the Chilean Participation and nearly 200,000 Emission-Line Galaxies (ELGs) to Group joining the collaboration. A second APOGEE extend the final BOSS galaxy clustering measurements spectrograph is being constructed for installation on the (Alam et al. 2017) by providing two new BAO dis- 2.5-meter du Pont Telescope (Bowen & Vaughan 1973) tance measurements over the redshift interval 0.6 < z < at Las Campanas Observatory (LCO) near La Serena, 1.1 Roughly 500,000 spectroscopically-confirmed quasars Chile. When APOGEE-2S begins survey operations in will be used as tracers of the underlying matter den- 2017, approximately 300 survey nights on the du Pont sity field at 0.9 < z < 2.2, providing the first mea- Telescope will be used to extend the APOGEE-2 survey surements of BAO in this redshift interval. Finally, the to the Southern Hemisphere. Lyman-α forest imprinted on approximately 120,000 new Data Release 13 (DR13) is the first data release for quasar spectra will give eBOSS an improved BAO mea- SDSS-IV, which will have regular public, documented surement over that achieved by BOSS (Delubac et al. data releases, in keeping with the philosophy of SDSS 2015; Bautista et al. 2017). The three new tracers will since its inception. In this paper, we describe the data provide BAO distance measurements with a precision of available in DR13, focusing on the data appearing here for the first time. We present overall descriptions of the 2 http://www.nsatlas.org sample sizes and targeting and provide a detailed bibli- 6 SDSS Collaboration ography of the technical papers available to understand tended Hα emission that fell into two broad classes based the data and the surveys in more detail. These technical on whether the intensity of emission is flat or increasing papers and the SDSS website (www.sdss.org) contain with radius. critical information about these data, which here is only The ages of stellar populations across CALIFA’s wide summarized. range of Hubble types show that galaxies with stellar 9 2. OVERVIEW OF THE SURVEY LANDSCAPE masses >∼ 5×10 M form from the inside out, with the more massive galaxies having older stellar populations in The release of DR13 coincides with the arrival of an each scaled radial bin (P´erezet al. 2013). Cano-D´ıazet astonishgly rich set of data from ongoing and recently al. (2016) showed the relationship between star forma- completed surveys outside of SDSS. Blanton et al. (2017), tion rate (SFR) and stellar mass that held for entire Bundy et al. (2015), and Majewski et al. (2017) describe star-forming galaxies was also true for spatially resolved how the SDSS-IV surveys compare in survey strategy, regions within a galaxy. In another intriguing clue to the size, and data within the wider arena of spectroscopic relation between small and large scales, S´anchez et al. surveys. We cite here some key science results for these (2014) found that the HII regions in 306 CALIFA galax- works to complement our brief history of SDSS. ies have a characteristic disk oxygen abundance gradient Spectroscopic redshift surveys of large scale structure when scaled to its effective radius. have resulted in BAO measurements over a range of red- Another large IFU survey, the SAMI survey (Bryant et shifts. The BAO signal has been detected at lower red- al. 2015), is currently underway at the Anglo-Australian shift (z∼0.1) from measurements of ∼75,000 galaxies in Telescope and will utimately observe 3400 galaxies. The the 6dF Galaxy Redshift Survey (Beutler et al. 2011). Early Data Release (Allen et al. 2015) provided data for WiggleZ measured the BAO signal at similar redshifts 107 galaxies, while the first major data release (Green to BOSS based on redshifts of ∼225,000 galaxies. The et al. 2017) includes 772 galaxies. Science results from final WiggleZ results for the 1-D BAO peak (Kazin et the SAMI survey so far include characterizing the galac- al. 2014) and the 2-D BAO peak (Hinton et al. 2017) at tic winds or extended diffuse ionizing gas in edge-on disk z=0.44, 0.6 and 0.73 agree with the BOSS results. Still galaxies Ho et al. (2016), mapping the quenching of star underway is the VIMOS Public Extragalactic Redshift formation proceeding from the outside-in in dense en- Survey, which focuses on higher redshifts than previous vironments Schaefer et al. (2017), and identifying stellar work (0.5

TABLE 1 Reduced spectroscopic data in DR13

Target Category # DR12 # DR12+13 MaNGA main galaxy sample: PRIMARY v1 2 0 600 SECONDARY v1 2 0 473 COLOR-ENHANCED v1 2 0 216 MaNGA ancillary galaxies1 0 31 MaNGA Other 0 62 SEQUELS LRG IZW 11781 21271 LRG RIW 11691 20967 QSO EBOSS CORE 19455 33367 QSO PTF 13227 22609 QSO REOBS 1357 2238 QSO EBOSS KDE 11836 20474 QSO EBOSS FIRST 293 519 QSO BAD BOSS 59 95 QSO BOSS TARGET 59 95 DR9 CALIB TARGET 28594 49765 SPIDERS RASS AGN 162 275 SPIDERS RASS CLUS 1533 2844 TDSS A 9412 17394 TDSS FES DE 40 70 TDSS FES DWARFC 19 29 TDSS FES NQHISN 73 148 TDSS FES MGII 1 2 TDSS FES VARBAL 55 103 SEQUELS PTF VARIABLE 700 1153 APOGEE-2 All Stars 164562 164562 NMSU 1-meter stars 894 894 Telluric stars 17293 17293 Commissioning stars 11917 11917

1 Many MaNGA ancillary targets were also targeted as part of the main galaxy sample, and are counted twice in this table.

TABLE 2 Value-Added Catalogs New in DR13

Catalog Description Reference http://data.sdss.org/sas/dr13/ SPIDERS Clusters demonstration sample catalog Clerc et al. (2016) eboss/spiders/analysis/ SPIDERS AGN target selection catalog Dwelly et al. (2017) eboss/spiders/target/ SPIDERS cluster target selection catalog Clerc et al. (2016) eboss/spiders/target/ WISE Forced Photometry Lang et al. (2016) eboss/photoObj/external/WISEForcedTarget/301/ Composite Spectra of Emission-line Galaxies Zhu et al. (2015) eboss/elg/composite/v1_0/ ELG Fisher selection catalog Delubac et al. (2017) eboss/target/elg/fisher-selection/ Redmonster redshift & spectral classification catalog Hutchinson et al. (2016) eboss/spectro/redux/redmonster/v5_9_0/v1_0_1/ QSO Variability Palanque-Delabrouille et al. (2016) eboss/qso/variability/ APOGEE DR13 red-clump catalog Bovy et al. (2014) apogee/vac/apogee-rc/cat/

New services are also available to users once they regis- torial examples for accessing SDSS data is provided at ter on the SciServer portal, and these services work seam- http://www.sdss.org/dr13/tutorials. All flat files lessly with the existing tools. Most notable among the are available for download at http://data.sdss.org/ new offerings are SciDrive, a distributed DropBox-like datamodel/ file store; SkyQuery, a federated cross-matching service 5. RECALIBRATION OF IMAGING DATA that compares and combines data from a collection (fed- eration) of multi-wavelength archives (SkyNodes); and DR13 includes a photometric recalibration of the SDSS SciServer Compute, a powerful new system for upload- imaging data. Finkbeiner et al. (2016) rederived the g, r, ing complex analysis scripts as Jupyter notebooks (using i and z band zero points using the PS1 photometric cali- Python, MatLab or R) running in Docker containers. brations of Schlafly et al. (2012), as well as rederiving the Links to all of these methods are provided at flat fields in all five bands (including u) . This effort im- the SDSS website (http://www.sdss.org/dr13/data_ proved the accuracy of the SDSS photometry previously access) The data processing software for APOGEE, hindered by a paucity of overlapping scans to perform BOSS, and SEGUE are publicly available at http:// ubercalibration¨ across the entire SDSS sky and by sub- www.sdss.org/dr13/software/products. A set of tu- optimal flatfields. The residual systematics, as measured Data Release 13 9

Fig. 1.— In grey are shown the locations in equatorial coordinates of all possible plates with MaNGA galaxies, each containing 17 galaxy targets. Because the MaNGA targets are selected to have SDSS photometry, this footprint corresponds to the Data Release 7 imaging data (Abazajian et al. 2009). Approximately 30% of these will be observed in the full planned MaNGA survey. The blue shows the plates observed in the first year of MaNGA for which data cubes are released in this paper.

Fig. 2.— Coverage of DR13 data from BOSS and SEQUELS in equatorial coordinates. The blue areas show the locations in equatorial coordinates of the five new BOSS plates (SGC) and the 126 SEQUELS plates (NGC) released in DR13. The green represents the area covered by BOSS in DR12. The SEQUELS plates released in DR12 lie in the same region as the new ones in DR13, providing complete coverage over roughly 400 square degrees. 10 SDSS Collaboration by comparison with PS1 photometry and spectral energy valuable first step in MaNGA’s goals to reveal the inter- distributions for stars with spectra, are reduced to 0.9, nal properties and dynamics of a statistically powerful 0.7, 0.7 and 0.8% in the griz bands, respectively; sev- sample of galaxies, that spans a broad range in stellar eral previously uncertain calibrations of specific runs are mass, local environment, morphology, and star formation also now much better constrained. The resulting recali- history. Individual observations across the sample are of brated imaging catalogs are the basis for the eBOSS and sufficient quality to characterize the spatially-dependent MaNGA targeting. composition of stars and gas as well as their internal For the MaNGA target selection, we are using the kinematics, thus providing important clues on growth NASA-Sloan Atlas (NSA; Blanton et al. 2011), a reanal- and star formation fueling, the build-up of spheroidal ysis of the SDSS photometric data using sky subtraction components, the cessation of star formation, and the in- and deblending better tuned for large galaxies. Rela- tertwined assembly history of galaxy subcomponents. tive to the originally distributed version of that cata- The survey was made possible through an instrumen- log, we have used the new calibrations mentioned above, tation initiative in SDSS-IV to develop a reliable and increased the redshift range up to z = 0.15, and have efficient way of bundling 1423 optical fibers into tightly- added an elliptical aperture Petrosian measurement of packed arrays that constitute MaNGA’s IFUs (Drory flux, which MaNGA targeting is based upon. DR13 in- et al. 2015). For each pointing, MaNGA observes 17 cludes the NSA catalog (version v1 0 1) associated with science targets with IFUs ranging from 19 to 127 fibers this reanalysis as the nsatlas CAS table and as a file on (with diameters of 12−3200). The IFU size distribu- the SAS. For the MaNGA galaxies released in DR13, we tion was optimized in concert with the sample design provide the actual images (referred to in MaNGA docu- (Wake et al. 2017) that targets SDSS-I/II main sam- mentation as “preimaging”) on the SAS as well. ple galaxies at a median redshift of 0.03 to obtain in Lang et al. (2016) reanalyzed data from the Wide-field 6 years a sample of 104 galaxies with uniform radial cov- Infrared Satellite Explorer (WISE; Wright et al. 2010) to erage and a roughly flat distribution in log M∗ limited 9 use for eBOSS targeting. They used positions and galaxy to M∗ > 10 M . Careful attention was paid to opti- profile measurements from SDSS photometry as input mizing hardware and an observing strategy that ensures structural models and constrained fluxes in the WISE high quality imaging spectroscopy (Law et al. 2015) and 3.4 µm and 4.6 µm bands. These results agree with to surface photometric flux calibration with a precision the standard WISE reductions to within 0.03 mag for better than 5% across most of the wavelength range, high signal-to-noise ratio, isolated point sources in WISE. 3,622-10,354A(Yan˚ et al. 2016). As described in these However, the new reductions provide flux measurements papers, salient aspects included protocols for constrain- for low signal-to-noise ratio (< 5σ) objects detected in ing hour-angles of observations to limit differential at- the SDSS but not in WISE (over 200 million objects). mospheric diffraction, dithering exposures to avoid gaps Despite the fact that the objects are undetected, their in the coverage of the targets because of space between flux measurements are nevertheless informative for tar- the fibers, and special calibration mini-bundles to ensure get selection, in particular for distinguishing stars from reliable absolute and relative flux calibration. An au- quasars. This photometry is provided as a value-added tomated pipeline delivers sky-subtracted, flux-calibrated catalog in the wiseForcedTarget CAS table and on the row-stacked-spectra (RSS) and datacubes for all sources SAS as described in Table 2. (Law et al. 2016). The Galactic extinction estimates published in the SDSS imaging tables (photoObj and related tables in the 6.1. MaNGA DR13 Main Galaxy Sample CAS) have been changed. The Galactic extinction still At the completion of SDSS-IV, the MaNGA survey’s uses the Schlegel et al. (1998) models of dust absorp- main galaxy sample will include ∼ 104 galaxies with 9 tion to estimate E(B − V ), but the Galactic extinction M∗ > 10 M and a roughly flat stellar mass distribu- coefficients for each band have been updated as recom- tion. DR13’s 1390 galaxy data cubes, corresponding to mended by Schlafly & Finkbeiner (2011). The extinction 1351 unique galaxies makes it the largest public sample of coefficients Ru, Rg, Rr, Ri, and Rz are changed from the galaxies with IFU spectroscopy. MaNGA’s main galaxy values used in BOSS (5.155, 3.793, 2.751, 2.086, 1.479) to sample consists of three major parts: Primary sample, updated values (4.239, 3.303, 2.285, 1.698, 1.263). The Secondary sample, and the Color-Enhanced supplement. corresponding numbers for the WISE bands are RW 1 = The Primary sample and the Secondary sample are 0.184 for the WISE 3.4µm band and RW 2 = 0.113 for selected from two luminosity-dependent redshift bands, the 4.6µm band (Fitzpatrick 1999). as shown in Figure 4. The selection for each sample is volume-limited at each absolute i-band magnitude. The 6. MANGA: INTEGRAL FIELD SPECTROSCOPIC DATA shape of the redshift bands is motivated by MaNGA’s MaNGA is investigating the internal kinematics and science requirements of having a uniform spatial cover- composition of gas and stars in low redshift (z ≤ 0.15) age in units of galaxy’s effective radius Re and having galaxies using fiber bundles to feed the BOSS spectro- a roughly flat stellar mass distribution (Yan et al. 2016; graphs. Bundy et al. (2015) describe the high-level sci- Wake et al. 2017). Figure 5 shows the distribution of the entific goals, scope, and context of the survey in inves- MaNGA DR13 galaxies in the stellar mass vs. dark mat- tigating galaxy formation while Yan et al. (2016) give a ter halo mass plane. Because more massive galaxies are detailed description of the survey design, execution, and on average larger, we observe them at a larger distance data quality relevant to DR13. With nearly 1390 data than low mass galaxies. We chose the redshift bands cubes released, MaNGA’s DR13 data products represent so that the great majority of the Primary (Secondary) the largest public sample to date of galaxies observed sample is covered by our fiber bundles to 1.5 Re (2.5 Re) with integral field spectroscopy. This data set signifies a along their major axes. This has the commensurate ef- Data Release 13 11

Fig. 3.— Coverage of APOGEE-2 DR13 data in Galactic coordinates; the raw data and its coverage is the same as in DR12, but it has been reprocessed through the latest reduction pipeline and ASPCAP versions. The color coding denotes the number of visits to each field, as indicated in the legend. 12 SDSS Collaboration fect of changing the physical resolution systematically as few got observed during the first year, while some a function of stellar mass, as illustrated in Figure 4. Po- have no observations during this period. More tar- tential deleterious effects of this change in sampling are gets for these programs will be observed in the fu- addressed by an ancillary program, described below. ture. To identify ancillary targets, one should use the We also designed a Color-Enhanced supplement, as a MANGA TARGET3 bitmask (or mngtarg3 in the dr- supplement to the Primary sample, to enhance the sam- pall file). All ancillary targets have MANGA TARGET3 pling of galaxies with rare color-magnitude combinations, greater than zero. Additional information on the scien- such as low-mass red galaxies, high mass blue galaxies, tific justification and targeting for each ancillary program and green valley galaxies. This is achieved by extending can be found at http://www.sdss.org/dr13/managa/ the redshift limits around the Primary sample redshift manga-target-selection/ancillary-targets. Here band for each underpopulated region in color-magnitude we provide some highlights and the corresponding bit- space. names. The combination of the Primary sample and the Color- Enhanced supplement is referred to as the Primary+ • Luminous AGN: This program increases the num- sample. We provide in our data release the redshift limits ber of host galaxies of the most luminous active over which each galaxy is selected. This permits a correc- galactic nuclei (AGN). The first source of targets is tion to the sample using 1/Vmax weight to reconstruct a the BAT 70-month Hard X-ray catalog (Baumgart- volume-limited representation of the galaxy population, ner et al. 2013). These have the bitname AGN BAT. provided that there is negligible galaxy evolution over To increase the sample size further we used the this limited redshift range. More details of how we ar- [OIII]-selected catalog of Mullaney et al. (2013) rived at this selection can be found in Yan et al. (2016) (bitname AGN OIII). To match the distribution of and Wake et al. (2017). Wake et al. (2017) provides the bolometric luminosities between the samples, we details of how to properly weight the sample to recon- selected 5 [OIII] objects at comparable Lbol to each struct a volume-limited representation. BAT object, within a redshift range of 0.01-0.08. Among the 1351 unique galaxies released in DR13, The bolometric corrections are from Shao et al. there are 600 Primary Sample galaxies, 473 Secondary (2013) and Vasudevan & Fabian (2009). Sample galaxies, and 216 Color-Enhanced supplement galaxies. There are 62 galaxies that do not belong to any • Edge-On Starbursts: We will use edge-on star- of the above. Some of these are ancillary program targets bursts to study the morphology and ionization (see below), some are filler objects on plates with spare state of large-scale outflows. To identify good tar- bundles, and others are galaxies selected using older, ob- gets, the specific star formation rate (sSFR) and solete versions of the selection and observed on early inclination of every object in the baseline MaNGA plates. For most statistical analyses, these 62 galaxies targeting catalog was calculated. The sSFR was should be excluded. determined using WISE photometry from Lang et Which sample a given target galaxy belongs to is given al. (2016) and the calibration between the W4 filter by the MANGA TARGET1 bitmask (or mngtarg1 in the and 22 µm emission in Jarrett et al. (2013) . We “drpall” file; see §5.2). Primary sample galaxies have bit then use a calibration from Cluver et al. (2014) to 10 set to 1, Secondary Sample galaxies have bit 11 set to derive the sSFR. The axis ratio SERSIC BA in the 1, and Color-Enhanced Supplement Galaxies have bit 12 targeting catalog was used to derive the inclination. set to 1. Bits 1-9 are for obsolete selections and should All targets in this program have log sSFR> −8.75 be ignored. and inclination > 75 deg.. The four galaxies in DR13 that have these properties, but were not in- 6.1.1. MaNGA Galaxy Ancillary Programs cluded in the main galaxy target sample, have bit- Roughly 5% of the IFUs are assigned to ancillary sci- name EDGE ON WINDS. ence programs defined by and allocated through internal competition and review. This assignment takes advan- • Close Pairs and Mergers: Interactions and merg- tage of sky regions with a low density of galaxies defined ers can play a key role in galaxy evolution, and by our main survey criteria (above) or where certain rare therefore an ancillary program was accepted that classes of objects, possibly outside our selection cuts, are either slightly adjusted the field centers for some of sufficiently high science value to re-allocate IFUs from targets already included in the main galaxy sam- the main program. Such high-value targets sometimes ple or placed new IFUs on galaxies. Close pairs come from the main sample, but the ancillary science are defined as galaxies in the NSA catalog or the goals require a different bundle size, a slightly different SDSS group catalog of Yang et al. (2007) and X. center positions, or higher prioritization over the random Yang (private communication) with projected sep- selection among the main sample galaxies. There are also aration < 50 kpc h−1 and line-of-sight velocity science cases where using observing strategies different (dV) < 500 km s−1, if both redshifts are available. from standard MaNGA observations are required. These If the bitname is PAIR ENLARGE, then to get the lead to dedicated plates. full pair required a larger IFU than the one origi- We solicited ancillary proposals in all these cate- nally scheduled by the targeting algorithm (Wake gories during the first year of survey observations, and et al. 2017). If the bitname is PAIR RECENTER this they started to be included in plate design half way means the original assigned MaNGA IFU is suffi- through this year. Consequently, the ancillary frac- ciently large, but requires re-centering. In addition tion in DR13 is smaller than 5%. Some ancillary to these already-planned galaxies, two sources of programs have tens of targets approved but only a new objects were used. The one galaxy in DR13 Data Release 13 13

−1 Fig. 4.— Principal selection cuts for the main MaNGA samples, where h = H0/100 km s . The colored bands show the selection cuts for the Primary (orange) and Secondary (green) samples illustrating the Mi dependence of the redshift limits. More luminous, and hence typically larger galaxies are selected at higher redshift than less luminous galaxies, ensuring that the angular size (1.5 Re or 2.5 Re) distribution is roughly independent of luminosity. The volume sampled also increases as the luminosity increases in such a way as to ensure a constant number density of galaxies at all luminosities. The points show the positions in this plane for the MaNGA galaxies in DR13, for the Primary (red triangles), Secondary (green squares) and Color-Enhanced samples (blue stars), although the Color-Enhanced selection also depends on NUV-icolor (see text for details). The numbers in the legend give the total number of observations of galaxies in each class, including repeat observations. The right hand y-axis gives an indication of the physical size of the mean spatial resolution element of the MaNGA data.

with PAIR SIM comes from the Galaxy Zoo Merg- cluding some of these analogs in the main galaxy ers Sample (Holincheck et al. 2016). A critical catalog, but they are slightly biased or deficient in sample comes from the ancillary program that re- certain regions of parameter space. Galaxies with quests that each galaxy be assigned a separate IFU the bitname MWA are drawn from the Licquia et al. PAIR 2IFU. Only four pairs of interacting galaxies (2015) catalog to provide galaxies in those under- are serendiptiously targeted in the main MaNGA represented parts of parameter space. galaxy sample with separate IFUs, and this sam- ple will compensate for the strong bias in the sin- • Dwarf Galaxies in MaNGA: The MaNGA main 9 gle IFU sample towards close separations or higher galaxy sample has galaxies with M∗ > 10 M , but redshifts. dwarf galaxies are the most numerous galaxies in the Universe. This ancillary program provides 2 • Massive Nearby Galaxies: Because the largest dwarf galaxies with MaNGA observations in DR13, MaNGA IFU covers 3200, more luminous, larger the first observations of a larger sample expected by galaxies observed to the same effective radius have the end of the survey covering a range of environ- poorer spatial resolution. The one MASSIVE ancil- ments. These galaxies are indicated by the bitname lary target in DR13 is part of a program to obtain DWARF and are drawn from the Geha et al. (2012) 9 a sample of nearby large galaxies with spatial reso- galaxy catalog with stellar masses < 10 M . lution better than 3 kpc and similar to the faintest galaxies in the MaNGA primary sample, at the cost • Brightest Cluster Galaxies: The brightest cluster of spatial extent. galaxies (BGCs) targeted here are brighter and in more massive halos than BCGs already in the • Milky Way Analogs: Licquia et al. (2015) defined MaNGA main sample and have the bitname BCG. a sample of Milky Way analogs based on M∗, SFR, We base our target selection on the updated Yang absolute magnitudes, and colors. MaNGA is in- et al. (2007) cluster catalog, created from the SDSS 14 SDSS Collaboration

Fig. 5.— The location of MaNGA galaxies (black triangles) in this data release area in the plane of stellar mass (M∗, x-axis) versus dark matter halo mass (Mh, y-axis). The two panels show red and blue galaxies separately,which are divided by a single color cut at g − r = 0.7. Stellar masses are taken from the NASA-Sloan Atlas, and halo masses from SDSS/DR7, using the method of Yang et al. (2007). Plotted as a colormap in the background are the number of MaNGA galaxies predicted for a 6-year full survey based on mock catalogs informed by the semi-analytic model of Guo et al. (2011) and constructed as in Li et al. (2006), which are the same as Figure 3 of Bundy et al. (2015). Normalized histograms show 1D marginalized M∗ distributions (top axes) and Mh distributions (right axes), with dashed lines for the full primary sample and solid lines for the red (left) and blue (right) populations.

DR7 NYU VAGC, an update of the DR4 version spectra for a given galaxy together onto a regularized of Blanton et al. (2005) 0.500 grid (see L16 for more detail). Both data cubes and RSS files are provided in a version with a log wavelength • MaNGA Resolved Stellar Populations: The ancil- scale, which is the standard extraction and is relatively lary program targets NGC 4163, a nearby dwarf smooth in velocity space, and in a version with a linear galaxy with existing HST imaging and high-quality wavelength scale, created directly from the native pixel color-magnitude diagrams selected from the ACS solution rather than by resampling the log-scaled spec- Nearby Galaxy Survey (Dalcanton et al. 2009). tra resampling. Each MaNGA data cube has associated This galaxy is flagged by the bitname ANGST. extensions corresponding to the estimated inverse vari- ance per pixel and a bad-pixel mask containing informa- tion about the quality of a given pixel within the cube 6.2. MaNGA Data Products: Individual Fiber Spectra (depth of coverage, bad values, presence of foreground and 3-D Data Cubes stars, etc). Additional extensions provide information In DR13, MaNGA is releasing both raw data (in the about the instrumental spectral resolution, individual ex- form of individual CCD frames) and reduced data pro- posures that went into the composite data cube, recon- duced by version 1 5 4 of the MaNGA Data Reduction structed griz broadband images created from the IFU Pipeline (DRP). Figure 6 illustrates the quality of the spectra, and estimates of the griz reconstructed point spectra from this pipeline. The MaNGA observing strat- spread function. egy is described by Law et al. (2015), and the flux cal- The objects observed by MaNGA for which data cubes ibration by Yan et al. (2016). Details of the MaNGA have been produced are summarized in the “drpall” file, DRP, data products, and data quality are given by Law a FITS binary table with one entry per object (including et al. (2016) (hereafter L16). All MaNGA data are in both galaxies and spectrophotometric standard stars ob- the form of multi-extension FITS files. served with 7-fiber IFUs to calibrate the MaNGA data). The DRP data products consist of intermediate re- This drpall file lists the name, coordinates, targeting in- duced data (sky-subtracted, flux-calibrated fiber spectra formation (e.g., redshift as given by the NASA Sloan with red and blue data combined for individual exposures Atlas), reduction quality, and other quantities of inter- of a plate) and final-stage data products (summary row- est to allow users to identify galaxy targets of interest. stacked spectra and data cubes) for each target galaxy. We note that MaNGA adopts two naming schemes. The The summary row stacked spectra (RSS files) are two- first, termed “mangaid” is an identifier unique to a given dimensional arrays provided for each galaxy in which astronomical object (e.g., 1-266039). The second, the each row corresponds to a single fiber spectrum. For “plate-ifu” uniquely identifies a given observation by con- a 127-fiber IFU with 9 exposures, there are thus 127 × catenating the plate id with the IFU number (e.g., 7443- 9 = 1143 rows in the RSS file. These RSS files have ad- 12701 identifies the first 127-fiber IFU on plate 7443). ditional extensions giving astrometric information about Since some galaxies are observed more than once on dif- the wavelength-dependent locations of each fiber on the ferent plates, a given mangaid can sometimes correspond sky. to more than one plate-ifu. The three-dimensional data cubes (axes R.A., Decl., The mangaid consists of 2 parts separated by a hyphen. wavelength) are created by combining the individual The first part is the id of the parent catalog from which Data Release 13 15

TABLE 3 Summary of MaNGA Ancillary Programs with Data in DR13

Ancillary Program # of Targets in DR13 1 # of Total Targets 1,2 BITNAME binary digit value Luminous AGN 1 267 AGN BAT 1 2 4 AGN OIII 2 4 Edge-On Star-forming Galaxies 4 166 EDGE ON WINDS 6 64 Close Pairs and Mergers 5 510 PAIR ENLARGE 7 128 10 PAIR RECENTER 8 256 1 PAIR SIM 9 512 1 PAIR 2IFU 10 1024 Massive Nearby Galaxies 1 310 MASSIVE 12 4096 Milky Way Analogs 2 250 MWA 13 8192 Dwarf Galaxies in MaNGA 2 247 DWARF 14 16384 Brightest Cluster Galaxies 2 378 BCG 17 131072 MaNGA Resolved Stellar Populations 1 4 ANGST 18 262144

1 An individual galaxy can be targeted by more than one ancillary program. 2 Number for each Ancillary Program refers to all targets in that program, regardless of bit name a target was selected. The second part is the position through the web browser and queried using any within that catalog. For most galaxy targets the catalog program able to read FITS binary tables. Once a id is 1 which refers to the NSA. For a small number of set of galaxies of interest has been identified, in- the early targets the catalog id is 12 and refers to an dividual data cubes, summary RSS files, interme- earlier version of the NSA (v1b 0 2). All galaxies from diate data products, etc. can be downloaded by this earlier version of the NSA are also in the final version browsing through the web directory tree. and so we release photometry etc for those targets from the final version of the NSA (v1 0 1), which is included • Large downloads of many DRP data prod- in the data release. Other catalogs referred to in the first ucts can be automated using rsync access part of the manga-id are for SDSS standard stars. to the SAS. For instance, to download all The full data model for all MaNGA DRP data products MaNGA data cubes with a logarithmic wave- can be found online at http://www.sdss.org/dr13/ length format: \rsync-aLrvz--include"*/ manga/manga-data/data-model/ and is also given in "--include"manga*LOGCUBE.fits. Appendix B of L16. gz"--exclude"*"rsync://data.sdss.org/ dr13/manga/spectro/redux/v1_5_4/./ 6.3. Retrieving MaNGA Data • The MaNGA drpall file can also be queried on- The raw data, reduced data, RSS, and 3-D data cubes line using the SDSS CASJobs system at http: for 1351 MaNGA galaxies are provided in DR13. From //skyserver.sdss.org/casjobs. While this can these data products, maps of line ratios, spectroscopic be useful for identifying MaNGA observations of indices, and kinematics can be made using standard soft- interest, CASJobs does not contain links to the ware. Because the first step in using the MaNGA data MaNGA data cubes and another method must be for science is to retrieve the spectra, we detail here and used to download the data themselves. on the SDSS website4 how to access the MaNGA spectra. • The SDSS SkyServer Explore tool at http:// 6.3.1. Reduced Data Products skyserver.sdss.org/dr13/en/tools/explore/ MaNGA DR13 reduced data products are stored on the will display basic information about MaNGA Science Archive Server at http://data.sdss.org/sas/ observations in DR13 that fall within a given cone dr13/manga/spectro/redux/v1_5_4/. The top level di- search on the sky. The relevant explore pages also rectory contains the summary drpall FITS table and provide direct links to the FITS data cubes on the subdirectories for each plate. Inside each plate direc- SAS. tory there are subdirectories for each MJD on which the plate was observed, containing intermediate (exposure 6.3.2. Raw Data level) data products. The ‘stack’ subdirectory within All MaNGA data taken in the first year of SDSS- each plate directory contains the final RSS and cube files IV observations are part of Data Release 13, in- for each MaNGA galaxy, sorted according to their plate- cluding data from special ancillary plates and co- ifu identifiers. Note that the ifu identifier in the filenames observing during APOGEE-2 time that are not part indicates the size of the IFU; everything in the 127 series of the DRP results. The raw data are stored on (e.g., 12701) is a 127-fiber bundle, etc. The 700 series the SAS in the directory http://data.sdss.org/sas/ ifus (e.g., 701) are the twelve spectrophotometric 7-fiber dr13/manga/spectro/data/ in subdirectories based on minibundles that target stars on each plate. the MJD when the data were taken. The mangacore These are the ways of getting at the data in DR13: directory5 contains the information needed to figure out the RA and Dec positions of fibers on plates, the dithered • Direct html browsing of the SAS at the above MJDs to be combined to make the final spectrum in link. The file drpall-v1 5 4.fits can be downloaded 5 http://svn.sdss.org/public/repo/manga/mangacore/tags/ 4 http:www.sdss.org/dr13/manga/manga-data/data-access/ v1_2_3/ 16 SDSS Collaboration apocomplete directory, and information on the calibra- it is important to account for when attempting to sub- tion files. L16 and the http://www.sdss.org/dr13/ tract the instrumental resolution from barely-resolved manga/manga-data/metadata/ website contain the rel- lines. There are two effects that combine to produce evant information about the file formats and the use of this overestimate. The first is that the impact of the the calibration files to get to the fully reduced spectra. wavelength rectification on the effective spectral resolu- Because these files are prepared for internal use, they re- tion was not accounted for when combining spectra from tain many old features that should be ignored, such as blue and red cameras. The second is that the Gaussian the names assigned to targeting bits, which still retain width of the LSF reported by the DRP is strictly the the names from SDSS-I. width of a pre-pixellization gaussian, while most end- user analysis routines adopt post-pixellization gaussians 6.4. Notes on using MaNGA data instead (i.e., the different between integrating a gaus- There are several important caveats to keep in mind sian profile over the pixel boundaries vs evaluating a when working with MaNGA data. In this discussion we gaussian at the pixel midpoint). This will be treated treat only the MaNGA data cubes. Summary RSS, inter- more completely in a future data release; in the mean- mediate, and even raw data products present some statis- time a 10% correction to the instrumental LSF quoted tical advantages to the data cubes, in particular reduced by the MaNGA data products appears to be a reason- covariance between adjacent data points and greater ease able correction factor if using post-pixellization analy- of forward modeling, but are substantially harder to use. sis routines. However, because this correction factor it- First and foremost, each MaNGA data cube has a FITS self is uncertain, derived astrophysical line-widths sub- header keyword DRP3QUAL indicating the quality of stantially below the instrumental resolution should be the reduction. 1-2% of the data cubes are flagged as sig- viewed to have unreliable accuracy at this time. A full nificantly problematic for various reasons, ranging from discussion of issues to consider is available at http: poor focus to flux calibration problems. Table B13 in //www.sdss.org/dr13/manga/manga-caveats/. L16 lists the bits that can be set with this flag that de- scribe why the end product was deemed unsatisfactory. 6.5. Highlights of MaNGA Science with DR13 Data Galaxies for which the CRITICAL quality bit (=30) is set The MaNGA survey has produced a number of scien- in DRP3QUAL should be treated with extreme caution. tific results based on early data, indicating the breadth of While there may be some spaxels in that data cube that research possible with the MaNGA data. Here we briefly are acceptable, in general the CRITICAL bit indicates summarize the results of papers completed within the widespread problems with the data reduction. Each data SDSS-IV collaboration using MaNGA data on spatially cube also has an associated MASK cube describing the resolved gas physics, stellar population properties, and quality of individual spaxels in the data cube. This stellar and gas kinematics. These MaNGA papers serve MASK extension can be used to identify areas of no cov- as a guide to prospective users of SDSS-IV data in these erage outside the fiber bundle footprint6, low coverage specific science areas. The larger context for this science near the edge of the dithered footprint, problematic ar- is provided to a small degree in the Introduction, and eas due to detector artifacts, regions known to contain the papers here also provide citations to the extensive bright Milky Way foreground stars, etc. A simple sum- literature on each topic, to which we refer the interested mary DONOTUSE bit is of particular importance indicating reader for additional context. Citations to the science elements that should be masked out for scientific analy- in each area should refer to the original papers, whether ses. MaNGA-based or not, rather than this brief executive Since the MaNGA data cubes adopt a 0.500 sampling summary. size (chosen based on Fourier analysis in optimal observ- ing conditions to not truncate the high-k modes of the 6.5.1. Gas physics observational PSF), while individual fibers have a 200 di- ameter footprint, there is significant covariance between The spatially resolved emission-line measurements adjacent MaNGA spaxels that must be taken into ac- have been used to understand the physical conditions of count whenever combining spectra. A simplified method the ionized gas in galaxies. Cheung et al. (2016a) iden- tified AGN winds as a surprisingly common occurrence for doing this is discussed in §9.3 of L16. The typical reconstructed point spread function of the MaNGA data in normal, quiescent galaxies, suggesting these winds as cubes has FWHM of 2.500. potentially critical in suppressing star formation. These As discussed by L16, the instrumental line spread func- winds may help address the question of how star forma- tion (LSF) in the wavelength direction reported by the tion remains suppressed in early-type galaxies. Cheung various extensions within the MaNGA DR13 data prod- et al. (2016a) report bisymmetric emission features co- ucts is underestimated by about 10 ± 2%. This cor- aligned with strong ionized-gas velocity gradients in a rection is comparable to the errors in the reported res- galaxy from which they infer the presence of centrally olution of the BOSS spectrographs seen in single-fiber driven winds in typical quiescent galaxies that host low- work (e.g., Shu et al. 2012; Palanque-Delabrouille et al. luminosity active nuclei. These galaxies account for as much as ten per cent of the quiescent population with 2013). Although this makes little difference when deter- 10 mining the astrophysical width of broad spectral lines, masses around 2 × 10 M . They calculate that the en- ergy input from the galaxies’ low-level active supermas- 6 The fiber footprint is a hexagon, but the standard FITS image sive black hole is capable of driving the observed wind, data structure is based around rectangular arrays. There must which contains sufficient mechanical energy to heat am- therefore be blank areas around the live IFU footprint in order to bient, cooler gas (also detected) and thereby suppress inscribe the hexagon inside a bounding rectangle. star formation. Data Release 13 17

Fig. 6.— Example spectra from a typical MaNGA data cube, adapted from Yan et al. (2016). The inset shows the SDSS color image with the hexagonal IFU footprint overlayed. The top spectrum is from the central spaxel; the bottom spectrum from a spaxel 1.2 Re from the center and is multiplied by a factor of 6 for easier comparison with the central spectrum. The differences in the stellar and gas components between the two regions can clearly be seen, as well as the large number of diagnostic features to understand the star formation history and the physical conditions of the gas.

The broader nature of ionized gas in early-types has tion of MaNGA, hence the interpretation of diagnostic also been the subject of papers by Belfiore et al. (2016) diagrams and the gas-phase metallicity measurements. and Belfiore et al. (2017) following up on the analysis They quantify for the MaNGA data how the contam- of a small sample observed with the MaNGA prototype ination by DIG moves HII regions towards composite instrument in Belfiore et al. (2015). By using spatially of LIER-like regions. DIG significantly biases measure- resolved maps of nebular diagnostics and stellar popula- ments of gas metallicity and metallicity gradients because tion ages, this work has added substantial support to the at different surface brightness, line ratios and line ratio notion that evolved stellar populations provide the ion- gradients can differ systematically. ization source for a galaxy class that arguably should be A careful treatment of gas-phase metallicities inferred renamed from LINER (Low Ionization Nuclear Emission from spectral maps of galaxies has suggested a key role Region) galaxies to LIER galaxies. LIERs, it turns out, for the dependence of metallicity on local stellar mass are ubiquitous in both quiescent galaxies and in the cen- surface density. Barrera-Ballesteros et al. (2016) present tral regions of galaxies where star formation takes place the stellar surface mass density vs gas metallicity relation at larger radii. The study of Belfiore et al. (2016) and for more than 500,000 spatially-resolved star-forming re- Belfiore et al. (2017) have put the occurrence of the LIER gions from a sample of 653 disk galaxies. These galaxies phenomenon into a physically relevant framework that span a larger range in mass than in previous samples can be directly tied to the diversity of the galaxy popu- where the correlations were first discovered. They con- lation as a whole. Specifically, they identify two classes firm a tight relation between these local properties, with of galaxies as extended LIER (eLIER) and central LIER higher metallicities as the surface density increases. They (cLIER), respectively, and study their kinematics and find that even over this larger mass range this local re- stellar population properties. cLIERs turn out to be lationship can simultaneously reproduce two well-known mostly late type galaxies located around the green val- properties of disk galaxies: their global mass-metallicity ley, while eLIERs are morphologically early types and relationship and their radial metallicity gradients. Their are indistinguishable from passive galaxies devoid of line results support the idea that in the present-day universe emission in terms of their stellar populations, morphol- local properties play a key role in determining the gas- ogy and central stellar velocity dispersion. phase metallicity in typical disk galaxies. The widespread ionization state of LIER gas might However, Cheung et al. (2016b) have found a galaxy in originally manifest as the Diffuse Intergalactic Gas (DIG) the middle of a gas accretion event, providing a detailed which is intermixed with star-forming regions. Zhang look at what appears to be a relatively rare occurrence in et al. (2016) studied galactic DIG emission and demon- the nearby Universe of this mode of galaxy growth. They strate how DIG in star-forming galaxies impacts the mea- present serendipitous observations of a large, asymmet- surements of emission line ratios at the spatial resolu- ric Hα complex that extends ∼ 800 (∼ 6.3 kpc) beyond 18 SDSS Collaboration the effective radius of a dusty, starbursting galaxy. This dients, Johnston et al. (2017) demonstrate a new tech- Hα extension is approximately three times the effective nique using MaNGA data that seeks to decompose the radius of the host galaxy and displays a tail-like morphol- underlying population into contributions from different ogy. This is suggestive of an inflow, which is consistent physical sub-components. They explore how the disk with its relatively lower gas-phase metallicities and its and bulge components in galaxies reached their current irregular gaseous kinematics. states with a new approach to fit the two-dimensional light profiles of galaxies as a function of wavelength and 6.5.2. Stellar populations to isolate the spectral properties of these galaxies’ disks Spatially resolved stellar population properties and and bulges. The MaNGA data have a spatial sampling stellar growth histories have been analyzed, following the of 0.5 00 per pixel, and successful decompositions were analysis of a small sample observed with the MaNGA carried out with galaxies observed with the 61- to 127- prototype instrument by Wilkinson et al. (2015) and Li fiber IFUs with fields of view of 2200 to 3200 in diameter et al. (2015). MaNGA has explored the role of environ- respectively. ment in shaping the radial distribution of stellar ages Rembold et al. (2017, submitted) have identified a and metallicities, with particular attention given to the “control sample” to the MaNGA luminous AGN host potential measurement systematics. Using different spec- galaxies, matched on mass, distance, morphology, and tral fitting techniques and complementary environmental inclination. Their conclusions based on SDSS-III spectra metrics, both Goddard et al. (2017b) and Zheng et al. of the central region can be tested with evaluation of the (2017) conclude that any environmental signal in the av- stellar populations throughout the galaxies with MaNGA erage shape of gradients is weak at best, with no obvious data. trends emerging in the initial MaNGA data. Goddard et al. (2017a) studied the internal gradients 6.5.3. Gas and stellar kinematics of the stellar population age and metallicity within 1.5 Several studies are investigating the kinematics of both Re obtained from full spectral fitting and confirm sev- stars and gas across the galaxy population. Li et al. eral key results of previous surveys. Age gradients tend (2017) perform dynamical modeling on a more exten- to be shallow for both early-type and late-type galax- sive and diverse sample of elliptical and spiral galaxies ies. As well known from previous studies, metallicity than had previously been done. By comparing the stel- gradients are often complex (and not well-modeled by lar mass-to-light ratios estimated from stellar population linear or log-linear functions of radius), varying in detail synthesis and from dynamics, they find a similar system- from galaxy to galaxy on small radial scales. On aver- atic variation of the initial mass function (IMF) to pre- age, however, over radial scales of order 1 Re, MaNGA vious investigations. Early-type galaxies (elliptical and data provide the strongest statistical constraints to date lenticular) with lower stellar mass-to-light ratios or veloc- that metallicity gradients are negative in both early and ity dispersions within one effective radius are consistent late-type galaxies, and are significantly steeper in disks. with a Chabrier-like IMF while galaxies with higher stel- These results continue to indicate that the radial depen- lar mass-to-light ratios or velocity dispersions are consis- dence of chemical enrichment processes are far more pro- tent with a more bottom heavy IMF such as the Salpeter nounced in disks than they are in spheroids, and indeed IMF. the relatively flat gradients in early-type galaxies are in- Two studies have taken advantage of the large MaNGA consistent with monolithic collapse. For both early and sample in DR13 to quantify the frequency and attributes late-type galaxies, more massive galaxies have steeper of galaxies with strong disparities between gas and stellar negative metallicity and age gradients. Since early-type kinematics. Chen et al. (2016) find that an appreciable galaxies tend to be more luminous, the overall steeper fraction of galaxies have counter-rotating gas and stars. age and metallicity gradients in late-type galaxies reflect Counter-rotation is found in about 2% of all blue galax- the fact that the trends in these gradients with galaxy ies. The central regions of blue counter-rotators show mass are more pronounced for late-type galaxies. God- younger stellar populations and more intense star for- dard et al. (2017a) take advantage of the unique MaNGA mation than in their outer parts. Jin et al. (2016) have sample size and mass range to characterize this correla- further studied the properties of 66 galaxies with kine- tion between metallicity and age gradients and galaxy matically misaligned gas and stars. They find that the mass, which explains the scatter in gradients values seen star-forming misaligned galaxies have positive gradients in previous studies. in Dn4000 and higher gas-phase metallicity, while the Ibarra-Medel et al. (2016) meanwhile infer spatially- green valley/quiescent ones have negative Dn4000 gra- resolved stellar mass assembly histories for the MaNGA dients and lower gas-phase metallicity on average. De- galaxies, extending previously known relations between spite this distinction, there is evidence that all types of galaxy type and assembly history to a larger mass kinematically-misaligned galaxies tend to be in more iso- range. Their findings are consistent with blue/star- lated environments. They propose that misaligned star forming/late-type galaxies assembling, on average, from forming galaxies originate from gas-rich progenitors ac- inside to out. Red/quiescent/early-type galaxies present creting external kinematic decoupled gas, while the mis- a more uniform spatial mass assembly, or at least aligned green valley/quiescent galaxies might be formed one that has been dynamically well mixed since star- from accreting kinematic decoupled gas from dwarf satel- formation ceased, consistent with the flatter gradients lites. seen, e.g., Goddard et al. (2017a). In general, low-mass Finally, Penny et al. (2016) examine the kinematics of galaxies show evidence of more irregular global and spa- a sample of 39 quenched low-mass galaxies. The major- tial assembly histories than massive galaxies. ity (87%) of these quenched low mass galaxies exhibit In a developing effort to model stellar population gra- coherent rotation in their stellar kinematics, and a num- Data Release 13 19 ber host distinct disks or spiral features. Just five (13%) al. 2017. The spectra from these tests were also used −1 are found to have rotation speeds vcirc < 15 km s at in one of the first science results from eBOSS, a study 1 Re, and are dominated by pressure support at all radii. of galactic-scale outflows traced by UV emission (Zhu et Two of the quenched low mass galaxies (5%) host kine- al. 2015). The selection algorithm used in these fields is matically distinct cores, with the stellar population at quite different from what will be used in eBOSS and these their centers counter-rotating with respect to the rest of targets will not contribute to the final clustering sample. the galaxy. The results support a picture in which the For this reason, we summarize the findings of the LRG majority of quenched low mass (dE) galaxies have a disk and quasar spectra below but reserve discussion of the origin. ELG spectra for future work. DR13 does contain value-added catalogs relevant to 7. SEQUELS: EBOSS, TDSS, AND SPIDERS DATA ELG sample. Raichoor et al. (2016) describe the Fischer DR13 includes the data from 126 plates observed un- discriminant used to select ELG targets using photome- der the SEQUELS program. This program was started try from SDSS, WISE, and SCUSS (Zou et al. 2016) for in SDSS-III as an ancillary program to take advantage the main ELG sample. Delubac et al. (2017) produce of some of the dark time released when BOSS was com- the catalogs used for ELG targeting with SDSS+WISE pleted early. The SEQUELS targets were quite different or SDSS+WISE+SCUSS data. Finally, Zhu et al. (2015) from BOSS targets because the program was designed to generated composite, continuum-normalized, spectra of finalize the eBOSS target selection algorithms. The pri- emission-line galaxies in the rest-UV from ELGs observed mary targets were defined by two different SDSS-WISE by BOSS in various programs. Table 2 lists the location selection algorithms to determine the eBOSS LRG pro- of these files on the SAS. gram (Prakash et al. 2016) and several optical-IR and variability selections to determine the eBOSS quasar pro- 7.1.1. Luminous Red Galaxies from WISE colors gram (Myers et al. 2015). Likewise, objects selected from The increase in redshift complicates selection both by a combination of X-ray and optical imaging were used to shifting the 4000 A˚ break into the i-band filter and by determine the final SPIDERS cluster (Clerc et al. 2016) requiring fainter targets than those observed in BOSS. and AGN (Dwelly et al. 2017) programs while variabil- WISE 3.4 µm photometry (W1 band) was introduced to ity in PanSTARRS imaging was used to determine the enable selection of this higher redshift sample. As part final TDSS program (Morganson et al. 2015; Ruan et al. of the SEQUELS program, two overlapping selections for 2016). LRGs at higher redshifts (0.6 < z < 1.0, vs 0.4 < z < 0.7 66 SEQUELS plates were completed in the final year of for CMASS) were employed, allowing tests of potential SDSS-III. The remaining 60 plates required to fill out the strategies for eBOSS. Color cuts that combine optical 400 square degree footprint were completed in the first and infrared photometry were employed, enabling the year of SDSS-IV. As mentioned above, these data served targeting of LRGs at these redshifts while maintaining a a crucial role for verification of the eBOSS, TDSS and high purity. This selection technique takes advantage of SPIDERS target samples. SEQUELS and eBOSS LRG the strong peak at a rest frame wavelength of 1.6 microns spectra were used to optimize the performance of a new that is present in the spectrum of most galaxies. This redshift classification scheme that now reliably classifies peak enters W1 as the redshift gets closer to 1, yielding more than 90% of eBOSS LRG spectra (Hutchinson et large differences between the optical/IR colors of z >0.6 al. 2016), thus meeting the ambitious goal set forth at galaxies and stars. the beginning of the program (Dawson et al. 2016). The SEQUELS selected a total of ∼ 70, 000 LRGs over an sample also seeds the eBOSS footprint to be used for clus- area of ∼ 700 deg2 with 120.0 0.9 according to the (Comparat et al. 2016; Raichoor et al. 2016; Delubac et relative density of stars and quasars as a function of color, 20 SDSS Collaboration magnitude, and photometric uncertainty. The selection information on the X-ray sources, including flux mea- alone results in a highly complete sample of quasars to be surements, and a quantitative measure of the reliability used for clustering studies, but with contamination from of the association to optical and AllWISE data. The stars that is too large to fit into the eBOSS fiber budget. SPIDERS Galaxy Cluster target list (Clerc et al. 2016) The XDQSO selection is supplemented by morphology contains 94,883 and 3,839 objects for RASS and XMM cuts to remove low redshift AGN and optical-IR colors respectively. to remove stellar objects with blackbody spectra. Vari- In SEQUELS, SPIDERS used similar targeting cata- ability data from the Palomar Transient Factory (PTF; logs, also available as value-added products, to test tar- Law et al. 2009) are used to supplement the selection of geting strategies and provide initial results. The selection Lyman-α forest quasars, producing tracers with a den- criteria are somewhat different than the final SPIDERS sity of 3.2 deg−2 where sufficient PTF imaging data are algorithms. Full details are available in Clerc et al. (2016) available. In addition to cosmological measurements, the and Dwelly et al. (2017). We summarize the SEQUELS quasar sample can be used to study quasar astrophysics SPIDERS data available in DR13 below. and galaxy evolution through studies of the quasar lumi- nosity function, composite spectra, and multi-wavelength 7.2.1. Optical spectra of Galaxies in X-ray Identified spectral energy distributions spanning the radio to the X- Clusters ray. Myers et al. (2015) found that ∼ 96% of all quasar The cluster pilot study (Clerc et al. 2016) takes ad- targets with r < 22 will be confidently classified. Section vantage of the CODEX (Constrain Dark Energy with 5 of Myers et al. (2015) provides information on the prop- X-ray clusters; Finoguenov et al., in preparation) can- erties of quasars observed with SEQUELS, including the didate cluster list, which is based on currently available numbers observed in the redshift ranges 0.9 < z < 2.1 RASS data. As part of DR13, we provide the catalog and z > 2.1. of X-ray detected galaxy clusters spectroscopically con- DR13 includes value-added catalogs with the variabil- firmed using SEQUELS-DR12 SPIDERS spectroscopic ity measurements (see Palanque-Delabrouille et al. 2011) data. Galaxy clusters were identified through the emis- using either PTF data (Myers et al. 2015) or Stripe 82 sion of their hot baryonic component as extended X-ray data (Palanque-Delabrouille et al. 2016) that are used sources in RASS. The optical counterparts were found for selecting quasars in eBOSS based on variability. The by optimally searching photometric data for the exis- locations of the FITS tables are given in Table 2. tence of a red-sequence formed by their member galax- ies. Spectroscopic redshifts obtained by SPIDERS pro- 7.1.3. Redmonster and Improved Redshifts for LRGs vide definitive confirmation of the clustered nature of The DR13 redshifts for all SEQUELS targets are de- these objects and their redshift (up to z ∼ 0.6). We as- termined in the usual fashion, with best-fitting combi- signed cluster membership combining an algorithm and nations of PCA eigenspectra. Hutchinson et al. (2016) visual validation of individual objects. The gas proper- describe a new pipeline, redmonster that uses a suite of ties derived from X-ray observations (luminosity, tem- 7 discrete theoretical galaxy spectra as a basis to determine perature, R200 ) are derived using precise cluster red- 2 the most likely redshift through a minimum χ , rather shifts (∆z ∼ 0.001). We compute galaxy cluster veloc- than linear combinations of templates used in Bolton ity dispersions using several methods adapted to the low et al. (2012). This pipeline achieves a 90.5% automated number of spectroscopic members per system (of the or- redshift and spectral classification success rate for the der 15-40) and we show that their values correlate with LRG target class, a significant improvement over the per- cluster X-ray luminosities, within expectations. Figure 7 formance of the previous pipeline. A value-added cata- shows the distribution of clusters with SEQUELS red- log using the new redmonster algorithms for the LRG shifts and membership in the redshift-X-ray luminosity sample is included in DR13. All spectra identified by plane. EBOSS TARGET0 bit 1 or 2 were classified by redmonster. The file is named redmonsterAll-v5 9 0.fits and is 7.2.2. Optical spectra of X-ray-identified AGN found on the SAS as described in Table 2. The addition of X-ray-identified AGN to the suite of AGN with well-sampled redshifts helps complete the in- 7.2. SPIDERS in SEQUELS ventory of AGN and trace black hole growth through- The main goal of the SPIDERS survey is to character- out cosmic history. SDSS has been observing optically- ize a subset of X-ray sources identified by eROSITA using identified AGN since its inception under the main large- optical spectra from the BOSS spectrograph. The ex- scale structure surveys, special targeting programs, and tended sources will mostly be galaxy clusters, which can as “mistakes” in other targeting classes. In addition, be used for cosmology. The point sources will mainly be there were several BOSS ancillary programs focused on AGN, which can be used to study the evolution of black X-ray follow-up, including a highly complete program on holes across cosmic time. For the first phase of SDSS-IV, the XMM-XXL north field (Menzel et al. 2016). The when eROSITA data are not yet available, SPIDERS will SPIDERS SEQUELS program has added spectra of 274 be targeting based on ROSAT and XMM data. The tar- ROSAT AGN to the SDSS sample, identified on the ba- get catalogs for galaxy clusters and AGN for SPIDERS sis of their SDSS colors only. The DR12 paper (Alam et from these two satellites have been included in DR13 as al. 2015) and Dwelly et al. (2017) provide details on the value-added catalogs. The SPIDERS AGN target cat- targeting of these AGN. alogs (Dwelly et al. 2017) contain 9,028 candidate tar- 7 gets from RASS and 819 from XMMSL (XMM-Newton R200 is defined as usual as the radius where the mean overden- Slew survey catalog; Saxton et al. 2008). They enclose sity is equal to 200× the critical density. Data Release 13 21

Fig. 7.— Distribution in the X-ray luminosity-redshift plane of galaxy clusters, adapted from Clerc et al. (2016). The gray dots mark the location of all candidate CODEX clusters using their photometric redshift. They form the main pool of cluster targets in SPIDERS. The CODEX candidate clusters that have already been spectrocopically confirmed and included in the MCXC meta-catalog (Piffaretti et al. 2011) are indicated by open symbols. The confirmed clusters lie mainly at higher X-ray luminosities. The 230 confirmed clusters in the SPIDERS DR13 value-added catalog (blue dots) extend the redshift range of known clusters in the SDSS area to z ∼ 0.6 in a systematic way. Among these 230, only 18 match clusters from the MCXC meta-catalog. The new clusters in general lie at lower X-ray luminosities and therefore probe lower X-ray masses than previous ROSAT-based analyses. The 6 clusters indicated by red diamonds are MCXC clusters in the SEQUELS-DR12 footprint not present in SPIDERS (note that ABELL 1361 is within a masked area of the CODEX survey, so it is not matched).

7.3. TDSS in SEQUELS SEQUELS sub-sample. The classification of spectra was Nearly 18,000 targets selected or co-selected by TDSS initially done using the BOSS pipeline (Bolton et al. have been observed among the 126 SEQUELS spectro- 2012), but the spectra were also visually inspected. Over- scopic plates. The targeting strategy for TDSS in SE- all, the pipeline performance is outstanding, with the QUELS was very similar to that ultimately chosen for highest-level spectral classification (e.g., star vs. galaxy the bulk of SDSS-IV (Morganson et al. 2015). Ruan vs. quasar) in agreement with our visual inspection for et al. (2016) present TDSS spectroscopic results from about 97% of the TDSS spectra and with only 2% of the the 66 initial DR12 SEQUELS plates, along with a de- pipeline redshifts for quasars requiring significant adjust- scription of the small differences in targeting within SE- ment. QUELS. Figure 8 depicts results for the initial TDSS About 90% of the TDSS spectroscopic fibers are aimed 22 SDSS Collaboration at initial classification spectra for variables chosen with- • Abundances are derived for several more species out primary bias based on color or specific lightcurve (C I, P, Ti II, Co, Cu, Ge, Rb) than in DR12. character. Their variability is determined from within PS1 multi-epoch imaging (Kaiser et al. 2010), or via • Results are now available for stars with Teff < longer-term photometric variability between SDSS and 3500 using a newly employed MARCS model at- PS1 imaging surveys (Morganson et al. 2015). The initial mosphere grid. SEQUELS results reveal comparable numbers of stars • Separate synthetic spectral grids are used for and quasars among these photometric variables. A sum- dwarfs and giants; results for dwarfs include rota- mary for TDSS quasars is that the PS1/SDSS variability tion and different isotope ratios are used for dwarfs criteria mitigate some of the (well-known) redshift biases and giants. of color-selection yielding a smooth and broad quasar redshift distribution, and that TDSS selects relatively • An initial, approximate, attempt has been made to redder quasars (e.g., than the eBOSS core quasar sam- account for variable line-spread functions (LSFs) as ple) as well as a higher fraction of some peculiar AGN a function of fiber number. (such as BALQSOs, and BL Lacs); and among variable stars, TDSS selects significant numbers of active late- • The correction for telluric absorption has been im- type stars, stellar pulsators such as RR Lyrae, and eclips- proved, primarily from a better LSF characteriza- ing and composite binaries, along with smaller subsets tion, leading to better recovery of the stellar spec- of white dwarfs, cataclysmic variables, and carbon stars tra. (see Ruan et al. 2016). The other ∼10% of the TDSS targets are objects al- • The relation adopted from microturbulence has ready having one or more earlier epochs of SDSS I-III been refined, and a relation for macroturbulence spectroscopy, and for which TDSS is taking a second has been incorporated. (or sometimes 3rd or 4th etc.) epoch to reveal antici- More details are provided in Holtzman et al (2017, pated spectroscopic variability. These “Few Epoch Spec- in preparation; H16). Overall, these changes have im- troscopy” (or FES) targets include various subsets of proved and enhanced the APOGEE stellar parameters quasars and stars (Morganson et al. 2015). Recent ex- and abundances, but users need to continue to be aware ample science papers representative of this FES category of potential issues and data flagging, as discussed on the of TDSS, include work on acceleration of broad absorp- SDSS web site and H16. In particular, parameters and tion lines in BALQSOs (e.g., Grier et al. 2016), and abundances for low metallicity stars and for cool stars recent discoveries adding to the rare class of “changing are more uncertain. look quasars” (e.g., Runnoe et al. 2016). Subsequent to the freezing of the DR13 release, sev- eral issues were discovered with the “calibrated” effec- 7.4. Retrieving SEQUELS data tive temperatures and surface gravities that are released All SDSS data releases are cumulative and therefore with DR13. Based on good agreement of the spectro- the SEQUELS data, whether taken in SDSS-III or SDSS- scopic effective temperature with photometric effective IV, have been reduced with the latest pipelines and in- temperatures for the bulk of the sample that is near cluded with all previous SDSS optical spectra data in this solar-metallicity, no calibration was applied to the DR13 data release. SEQUELS targets can be identified because spectroscopic effective temperatures. It now appears, the EBOSS TARGET0 bit will be set. The summary spAll- however, that these are systematically offset from photo- metric temperatures for stars of lower metallicity by as v5 9 0.fits datafile, which includes classification informa- 8 tion from the pipeline, is at https://data.sdss.org/ much as 200-300 K for stars at [Fe/H] ∼ −2. For surface sas/dr13/eboss/spectro/redux/v5_9_0/ on the SAS gravities, a calibration was derived based on asteroseis- or in the specObjAll table on the CAS. mic surface gravities. While this calibration is good for most of the sample, it now appears that it does not yield accurate surface gravities for metal-poor stars; the cali- 8. APOGEE-2: IMPROVED SPECTRAL EXTRACTION AND 9 SPECTROSCOPIC PARAMETERS FOR APOGEE-1 DATA brated surface gravities for stars with [M/H] < −1.5 are systematically too low, by as much as ∼ 0.5 dex. H16 The data released in DR13 are based on the same raw and the SDSS web site discuss both of these issues, and data as in DR12, and the pipelines for reduction and suggest post-release calibrations. analysis remain similar to those used in DR13. First, the data reduction pipeline (Nidever et al. 2015) re- 8.1. Improvements in Linelist and Data Analysis duced the 3-D raw data cubes into well-sampled, com- Several improvements were made with regard to the bined, sky-subtracted, telluric-corrected and wavelength- APOGEE LSF. The characterization of the LSF was im- calibrated 1-D spectra. The stellar parameters and abun- proved in one of the detectors, and an initial attempt was dances were derived using ASPCAP (APOGEE Stellar made at accommodating the LSF variability in the pa- Parameters and Chemical Abundances Pipeline; Garc´ıa rameter and abundance pipeline by constructing spectral P´erezet al. 2016) by finding the χ2 minimum when com- libraries for 4 different LSFs and using the most appro- paring the normalized observed spectra against a grid of priate one for the analysis of each star. synthetic spectra. However, the processing and analysis have been improved in several ways: 8 [X/H]=(log10(NX/NH))∗−(log10(NX/NH)) . [X/Fe]=[X/H]∗− [Fe/H]∗. • The linelist used for determining stellar parameters 9 [M/H] is closely related to [Fe/H]. See M´esz´aroset al. (2013) and abundances has been revised. for more details. Data Release 13 23

Fig. 8.— A representation of the TDSS spectroscopic characterization of imaging variables selected (from PanSTARRS 1 and SDSS imaging) without any primary bias based on color or specific lightcurve character. Nearly 16,000 TDSS photometric variables in the region of the sky covered by the 66 initial SEQUELS plates are plotted here in a traditional color-color diagram. Their classifications are based on spectra taken in SDSS-I-III, including SEQUELS data. Regions in color space traditionally considered as occupied by quasars, main-sequence stars, RR Lyrae, high-redshift quasars, and other miscellaneous subclasses are overlayed (see Sesar et al. 2007, but here using specific boundaries detailed in Morganson et al. 2015). The blue symbols depict TDSS photometric variables whose SDSS spectra verify them as quasars, while the red symbols depict TDSS variables whose SDSS spectra verify them as stars. A few hundred objects that were identified as galaxies or which could not be identified are not included in this plot. The wide distribution, as well as the overlap (in some regimes relatively densely), of both stars and quasars symbols verify that TDSS variability selection finds not only traditional quasars, but also those within color regimes commonly attributed to stars (and with analogous results for the spectroscopically confirmed variable stars).

The linelist adopted for DR13, linelist.20150714, Al and Co. We still rely heavily on “astrophysical” gf is an updated version of the one used for DR12 values for atomic lines, where the transition probabilities results (Shetrone et al. 2015) and is available at of individual lines are determined using the spectrum http://data.sdss.org/sas/dr13/apogee/spectro/ of a star with known parameters and abundances, to redux/speclib/linelists/. Shetrone et al. (2015) construct linelists in the H-band (Mel´endez& Barbuy noted a number of concerns with the DR12 linelist, 1999; Ryde & Schultheis 2015). For our linelist, we which have been corrected in the new version (H16). use the center-of-disk spectrum of the Sun (Livingston As in DR12, the molecular line list is a compilation of & Wallace 1991) and the full disk spectrum of the literature sources including transitions of CO, OH, CN, nearby, well-studied, red giant Arcturus (Hinkle et al. C2, H2, and SiH. The CN line list has been updated 1995). To calculate the astrophysical gf values, we used from the DR12 version using a compilation from C. Turbospectrum (Alvarez & Plez 1998; Plez 2012) to Sneden (private communication). All molecular data generate synthetic spectra of the center-of-disk for the are adopted without change, with the exception of a Sun and the integrated disk for Arcturus with varying few obvious typographical corrections. The atomic line oscillator strengths and damping values to fit the solar list was compiled using a number of literature sources and Arcturus spectra. In DR12, a different stellar for the lower excitation state of the transition and the atmosphere code was used for the gf determination and transition probability, which is usually characterized as the synthetic grid creation and synthetic integrated disk the “gf” or “log gf” value. Literature sources include spectra were used for the both the Sun and Arcturus. theoretical, astrophysical, and laboratory gf values. Although the change in gf values from these changes A few new lines were added from NIST510 and other is small (Shetrone et al. 2015), the DR13 linelist is literature publication since the DR12 line list was more self-consistently generated than previous versions. created, including hyperfine splitting components for In DR12, we derived final astrophysical gf values by weighting the solar astrophysical gfs at twice those of 10 http://www.nist.gov/pml/data/asd.cfm 24 SDSS Collaboration

Arcturus. The astrophysical gf solutions in DR13 are reported abundances for elements such as P that weighted according to line depth in Arcturus and in show ∼1 dex scatter in Figure 9 are based on fits the Sun, to give more weight to where the relatively to noise in the spectrum and do not reflect actual weak lines in the H-band produce the best signal, which abundances. Abundances for such elements should usually gives more weight to the Arcturus solution. For not be used unless they are confirmed by visual lines with laboratory oscillator strengths, the astrophys- inspection. They are reported because for a cer- ical log(gf) values were not allowed to vary beyond twice tain range in metallicity and temperature, the lines the error quoted by the source. can be detected and interesting chemistry exposed (e.g., Hawkins et al. 2016). 8.2. Additional Elements • The disagreement between the [Ti/Fe] value based In DR13, APOGEE provides abundances for elements on Ti I lines and Ti II lines is yet to be resolved. P, Cr, Co, Cu, Ge, and Rb for the first time. The abun- We are currently investigating the choice of lines dances of C, N, O, Na, Mg, Al, Si, S, Ca, Sc, Ti, V, (as noted by Hawkins et al. 2016) and the effect of Mn, Fe, and Ni were re-derived. We added two new 3-D and NLTE corrections. species: C I and Ti II, which provide valuable checks on the abundances derived from CH and Ti I lines (H16). • The mean and scatter shown in Figure 9 were calcu- The technique for calculating abundances of individual lated using stars with 4100K < Teff < 5000K. The elements is described in Holtzman et al. (2015). As for scatter becomes noticeably worse if the warmest or DR12 and as discussed in Holtzman et al. (2015), we coolest stars in the clusters are included. Warmer correct abundance ratios [X/Fe], except for [C/Fe] and stars have weaker lines in general in the H-band and [N/Fe], so there is no trend with temperature among the the coolest stars are affected by the issues with the members of a single star cluster. In Figure 9, we show cool grid (see below). box plots for the abundances for member stars in four clusters with a range of metallicity. Both the mean val- The elements included in DR13 are by no means the ues and the rms illustrate key points about the APOGEE only elements/species with lines present in the H-band DR13 abundances. amenable to abundance analysis. Hawkins et al. (2016), for example, independently derives Yb, along with many • The elements cover a wide range in nucleosynthe- other elements, for the APOKASC sample (Pinsonneault sis sites and atomic number. APOGEE is measur- et al. 2014). We expect to include additional elements in ing the abundances of α-elements, odd-Z elements, upcoming data releases. iron-peak elements, and neutron-capture elements, as well as the mass, mixing and AGB-contribution 8.3. Synthetic Spectral Grids at Teff < 3500K and with diagnostics C and N. Rotational Broadening • The high [X/Fe] for the α elements in the metal- In DR10-12, the synthetic spectral grid used by ASP- poor ([Fe/H]< −0.7) globular clusters is present as CAP (M´esz´aroset al. 2012; Zamora et al. 2015) was re- expected (Wallerstein 1962). The increased scat- stricted to temperatures hotter than 3500K because Ku- ter in O, Mg, and Al for these same clusters is in rucz model atmospheres are not available at cooler tem- part due to the well-known light element anomalies peratures. However, many important APOGEE targets in globular clusters (e.g., Kraft 1994; M´esz´aroset have cooler temperatures, including the luminous metal- al. 2015). The lower [Mn/Fe] in metal-poor stars rich M-giants in the bulge, cool asymptotic giant branch is consistent with previous work (e.g., Wallerstein stars, and M dwarf planet hosts. Therefore, we used 1962; Gratton 1989) using 1-dimensional models custom MARCS (Gustafsson et al. 2008) atmospheres that assume local thermodynamic equilibrium. provided by B. Edvardsson (private communication) to construct new synthetic spectra. For the stellar atmo- • The increased scatter at lower metallicities is not spheres for the giants, the atmospheres are spherical, entirely the result of actual inhomogeneities in clus- otherwise they are plane-parallel. The Kurucz model ter members, as there are increasing uncertainties grid and the new MARCS model grid overlap between associated with the weaker lines in more metal- 3500-4000K. In this region, there are some systematic poor stars. The error estimates reported in DR13 differences between the results from the two grids; to for these stars track the measured rms reasonably provide homogeneous results, DR13 adopts the parame- well, especially for errors in [X/Fe] > 0.1 dex, and ters and abundances from the analysis with the Kurucz therefore provide a useful guide for abundance ra- grid. tios that are not well-determined (H16). On top Figure 10 shows the current parameter space in Teff , of this, the ASPCAP methodology breaks down at log g and [M/H] covered by Data Release 13 stars. It is some level for stars in which abundances within a immediately apparent that the parameters derived from given “abundance group” (e.g., the alpha-elements) the MARCS grid do not match smoothly to the param- depart from solar abundance ratios, as is the case eters from the Kurucz grid. Possible explanations in- for second-generation stars in metal-poor globular clude the switch from plane-parallel Kurucz to spheri- clusters (M´esz´aroset al. 2015). cal MARCS for the giants and/or the large number of model atmospheres in the MARCS grid that failed to • ASPCAP reports the parameters of the best-fitting converge. Because ASPCAP requires a square grid of model in χ2 space. It does not currently report synthetic spectra, these “grid holes” were filled by ad- non-detections and upper limits. Therefore, most jacent model atmospheres. The linelist does not have Data Release 13 25

Fig. 9.— A box plot of the stars with 4100K < Teff < 5000K in 4 clusters spanning a range in metallicities (M13= −1.58, M71= −0.82, M67= 0.06, NGC 6791= 0.37 (Holtzman et al. 2015). C and N are not included because their values in a cluster change depending of the evolutionary state of the star as the result of first dredge-up (e.g., Masseron & Gilmore 2015). A boxed horizontal line indicates the interquartile range and median found for a particular element. The vertical tails extending from the boxes indicate the total range of abundances determined for each element, excluding outliers. Outliers are shown by × symbols.

FeH lines, which can be important features in the atmo- jewski et al. 2017). Fewer than 2% of red giants rotate spheres of cool M dwarfs. We are investigating the size at speeds detectable at the APOGEE spectrograph res- of the error caused by using incorrect model atmospheres olution (de Medeiros et al. 1996; Carlberg et al. 2011). in the grid, the possibilities of using alternative methods Therefore the first versions of ASPCAP did not include of interpolating and identifying the best-fit model, and rotational broadening as a dimension in the synthetic the addition of FeH lines. H16 provides more details on spectral grid, which substantially reduced the computing the construction of the cool grid and the resulting stel- time. However, dwarf stars, of which there are >20,000 lar parameters. Nonetheless, Figure 10 illustrates the in APOGEE-1, are frequently rapidly rotating, especially improved parameter space over DR12, which should aid if they are young. For DR13, we added a dimension to in classifying M stars correctly as dwarfs or giants and the synthetic spectral grid where spectra were broadened separate the early M from the late M stars. by a Gaussian kernel. To keep the computing time rea- The APOGEE-1 sample is dominated by giants used sonable, and in acknowledgment of the small effect that to probe the chemical cartography of the Galaxy (Ma- C and N abundances have on the atmospheres of warm 26 SDSS Collaboration

Fig. 10.— The Hertzsprung-Russell diagram for all DR13 APOGEE-1 stars. The calibrated gravities and temperatures are used. The cool stars now clearly fall into dwarf and giant categories, but the stellar parameters do not merge smoothly onto the warmer Kurucz-based grid. dwarfs, we fixed the [C/M] and [N/M] grids to solar val- DR14 is scheduled for July 31, 2017, and there will ues. The [C/Fe] and [N/Fe] values reported for dwarfs in be yearly data releases until the end of the survey in DR13 are calculated from windows after the best-fit pa- 2020. For MaNGA, future data releases will include more rameters are determined, in similar fashion to the other data cubes of galaxies that are currently being observed. individual elements. Figure 11 shows the improvement DR14 will include 2744 additional galaxies observed by in the stellar parameters for members of the Pleiades MaNGA. In DR15, a data interface system that includes star cluster. Some of the dwarfs in this young cluster web- and python-based tools to access a data analysis are rotating with v sin i> 50 km/s. In DR12, ASP- pipeline that carries out continuum and emission-line CAP found a best-fit solution for metal-poor stars. The fitting to provide estimates of stellar and gas kinemat- shallower the lines because of rotational broadening, the ics, emission line fluxes, and absorption line index mea- more metal-poor the star was reported to be. With the surements will be released. In addition, MaNGA has DR13 grid, there is no longer a prominent trend in [M/H] started a bright-time observing program, piggy-backing with v sin i, and the mean value of the cluster is now at on APOGEE-2, to build a stellar library. These reduced the expected metallicity. optical stellar spectra will be included in DR15. For eBOSS, future data releases will provide sufficient red- 8.4. Data Access shifts of LRG, ELG, and quasars to be of cosmological Data access via the CAS and SAS is similar to that for interest, on its way to the goals described in §1. For DR12 (Holtzman et al. 2015); “dr12” should be replaced example, DR14 contains the spectra for 2480 square de- with “dr13” in the pathname or DR13 as the context grees observed by eBOSS. TDSS and SPIDERS will also in CASJobs. Some of the tag/column names have been release many more spectra, redshifts, and classifications modified. While raw abundances are still given in the for variable and X-ray emitting objects, respectively. For FELEM arrays, “calibrated” abundances are now pre- APOGEE-2, the next data release (DR14) will contain sented in X H and X M arrays as well as in individual spectra for an additional 98,882 stars from the APOGEE- X FE tags/columns. For dwarfs, there is now a column 2N spectrograph at APO observed under SDSS-IV, and in the allStar fits file on the SAS or in the aspcapStar DR16 and subsequent DRs will provide spectra taken table on the CAS that reports the v sin i. DR13 also both with the APOGEE-2N and the APOGEE-2S spec- includes a new catalog of red clump stars based on the trograph at LCO. In addition to observing red giants as new ASPCAP parameters, following the method of Bovy the main target category, data will be obtained for RR et al. (2014), available at the location described in Table Lyrae stars in the bulge, dwarf spheroidal galaxies, the 2. Magellanic Clouds, Kepler Objects of Interest, and tar- gets in the K2 fields. 9. THE FUTURE SDSS-IV will continue to add to the SDSS legacy of data and data analysis tools in upcoming data releases. 10. ACKNOWLEDGEMENTS Data Release 13 27

Fig. 11.— A comparison of the [M/H] derived for Pleiades cluster dwarfs from DR12 (red open circles) and DR13 (black filled squares). The Pleiades has many rapidly spinning stars with rotationally broadened, shallow lines. When analyzed with the DR12 version of ASPCAP, which did not include rotational broadening, the best χ2 fit for the high v sin i stars was achieved for synthetic spectra with metallicity much lower than that known for the Pleiades. Therefore, the cluster average was sub-solar (red line). In DR13, where the dwarfs are run through a grid that includes a v sin i dimension, all Pleiades stars fall close to the correct value (black line), regardless of broadening.

We would like to thank the Center for Data Science pation Group, the Carnegie Institution for Science, (http://cds.nyu.edu) at New York University for their Carnegie Mellon University, the Chilean Participa- hospitality during DocuFeest 2016. tion Group, the French Participation Group, Harvard- Funding for the Sloan Digital Sky Survey IV has been Smithsonian Center for Astrophysics, Instituto de As- provided by the Alfred P. Sloan Foundation, the U.S. trof´ısica de Canarias, The Johns Hopkins University, Department of Energy Office of Science, and the Partici- Kavli Institute for the Physics and Mathematics of pating Institutions. SDSS-IV acknowledges support and the Universe (IPMU) / University of Tokyo, Lawrence resources from the Center for High-Performance Com- Berkeley National Laboratory, Leibniz Institut f¨urAs- puting at the University of Utah. The SDSS web site is trophysik Potsdam (AIP), Max-Planck-Institut f¨urAs- www.sdss.org. tronomie (MPIA Heidelberg), Max-Planck-Institut f¨ur SDSS-IV is managed by the Astrophysical Research Astrophysik (MPA Garching), Max-Planck-Institut f¨ur Consortium for the Participating Institutions of the Extraterrestrische Physik (MPE), National Astronomi- SDSS Collaboration including the Brazilian Partici- cal Observatory of China, New Mexico State University, 28 SDSS Collaboration

New York University, University of Notre Dame, Obser- Arizona, University of Colorado Boulder, University of vat´arioNacional / MCTI, The Ohio State University, Oxford, University of Portsmouth, University of Utah, Pennsylvania State University, Shanghai Astronomical University of Virginia, University of Washington, Univer- Observatory, United Kingdom Participation Group, Uni- sity of Wisconsin, Vanderbilt University, and Yale Uni- versidad Nacional Aut´onomade M´exico,University of versity.

REFERENCES Abazajian, K. N., Adelman-McCarthy, J. K., Ag¨ueros,M. A., De Silva, G. M., Freeman, K. C., Bland-Hawthorn, J., et al. 2015, et al. 2009, ApJS, 182, 543 MNRAS, 449, 2604 Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, Delubac, T., Bautista, J. E., Busca, N. G., et al. 2015, A&A, 574, 219, 12 A59 Alam, S., Ata, M., Bailey, S., et al. 2017, MNRAS, 470, 2617 Delubac, T., Raichoor, A., Comparat, J., et al. 2017, MNRAS, Allen, J. T., Croom, S. M., Konstantopoulos, I. S., et al. 2015, 465, 1831 MNRAS, 446, 1567 Doi, M., Tanaka, M., Fukugita, M., et al. 2010, AJ, 139, 1628 Alvarez, R., & Plez, B. 1998, A&A, 330, 1109 Dong, S., Xie, J.-W., Zhou, J.-L., Zheng, Z., & Luo, A. 2017, Anguiano, B., De Silva, G. M., Freeman, K., et al. 2016, MNRAS, arXiv:1706.07807 457, 2078 Dong, S., Zheng, Z., Zhu, Z., et al. 2014, ApJ, 789, L3 Alatalo, K., Blitz, L., Young, L. M., et al. 2011, ApJ, 735, 88 Drory, N., MacDonald, N., Bershady, M. A., et al. 2015, The Barrera-Ballesteros, J. K., Heckman, T. M., Zhu, G. B., et al. Astronomical Journal, 149, 77 2016, MNRAS, 463, 2513 Dwelly, T., Salvato, M., Merloni, A., et al. 2017, MNRAS, 469, Baumgartner, W. H., Tueller, J., Markwardt, C. B., et al. 2013, 1065 ApJS, 207, 19 Eisenstein, D. J., Zehavi, I., Hogg, D. W., et al. 2005, ApJ, 633, Bautista, J. E., Busca, N. G., Guy, J., et al. 2017, A&A, 603, A12 560 Belfiore, F., Maiolino, R., Bundy, K., et al. 2015, MNRAS, 449, Eisenstein, D. J., Weinberg, D. H., Agol, E., et al. 2011, AJ, 142, 867 72 Belfiore, F., Maiolino, R., Maraston, C., et al. 2016, MNRAS, Fang, X.-S., Zhao, G., Zhao, J.-K., Chen, Y.-Q., & Bharat 461, 3111 Kumar, Y. 2016, MNRAS, 463, 2494 Belfiore, F., Maiolino, R., Maraston, C., et al. 2017, MNRAS, Fern´andez-Trincado, J. G., Robin, A. C., Vieira, K., et al. 2015, 466, 2570 A&A, 583, A76 Beutler, F., Blake, C., Colless, M., et al. 2011, MNRAS, 416, 3017 Finkbeiner, D. P., Schlafly, E. F., Schlegel, D. J., et al. 2016, ApJ, Blanton, M. R., Schlegel, D. J., Strauss, M. A., et al. 2005, AJ, 822, 66 129, 2562 Fitzpatrick, E. L. 1999, PASP, 111, 63 Blanton, M. R., Kazin, E., Muna, D., Weaver, B. A., & Frasca, A., Molenda-Zakowicz,˙ J., De Cat, P., et al. 2016, A&A, Price-Whelan, A. 2011, AJ, 142, 31 594, A39 Blanton, M. R., Bershady, M. A., Abolfathi, B., et al. 2017, AJ, Frieman, J. A., Bassett, B., Becker, A., et al. 2008, AJ, 135, 338 154, 28 Fukugita, M., Ichikawa, T., Gunn, J. E., et al. 1996, AJ, 111, 1748 Bolton, A. S., Schlegel, D. J., Aubourg, E.,´ et al. 2012, AJ, 144, Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2016, 144 A&A, 595, A2 Bovy, J., Nidever, D. L., Rix, H.-W., et al. 2014, ApJ, 790, 127 Garc´ıaP´erez,A. E., Allende Prieto, C., Holtzman, J. A., et al. Bovy, Jo, Hogg, D. W., & Roweis, S. T. 2011, Annals of Applied 2016, AJ, 151, 144 Statistics, 5, Geha, M., Blanton, M. R., Yan, R., & Tinker, J. L. 2012, ApJ, Bovy, J., Hennawi, J. F., Hogg, D. W., et al. 2011, ApJ, 729, 141 757, 85 Bowen, I. S., & Vaughan, A. H., Jr. 1973, Appl. Opt., 12, 1430 Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Brough, S., van de Sande, J., Owers, M. S., et al. 2017, Messenger, 147, 25 arXiv:1704.01169 Goddard, D., Thomas, D., Maraston, C., et al. 2017, MNRAS, Bryant, J. J., Owers, M. S., Robotham, A. S. G., et al. 2015, 466, 4731 MNRAS, 447, 2857 Goddard, D., Thomas, D., Maraston, C., et al. 2017, MNRAS, Bundy, K., Bershady, M. A., Law, D. R., et al. 2015, ApJ, 798, 7 465, 688 Cano-D´ıaz,M., S´anchez, S. F., Zibetti, S., et al. 2016, ApJ, 821, Gomes, J. M., Papaderos, P., Kehrig, C., et al. 2016, A&A, 588, L26 A68 Cappellari, M., Emsellem, E., Krajnovi´c,D., et al. 2011, Gratton, R. G. 1989, A&A, 208, 171 MNRAS, 413, 813 Green, A. W., Croom, S. M., Scott, N., et al. 2017, Carlberg, J. K., Majewski, S. R., Patterson, R. J., et al. 2011, arXiv:1707.08402 ApJ, 732, 39 Gunn, J. E., Carr, M., Rockosi, C., et al. 1998, AJ, 116, 3040 Chang, H.-Y., Song, Y.-H., Luo, A.-L., et al. 2017, ApJ, 834, 92 Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, Chen, Y.-M., Shi, Y., Tremonti, C. A., et al. 2016, Nature 131, 2332 Communications, 7, 13269 Guo, Q., White, S., Boylan-Kolchin, M., et al. 2011, MNRAS, Cheung, E., Bundy, K., Cappellari, M., et al. 2016, Nature, 533, 413, 101 504 Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, Cheung, E., Stark, D. V., Huang, S., et al. 2016, ApJ, 832, 182 486, 951 Clerc, N., Merloni, A., Zhang, Y.-Y., et al. 2016, MNRAS, 463, Guzzo, L., Scodeggio, M., Garilli, B., et al. 2014, A&A, 566, A108 4490 Hawkins, K., Masseron, T., Jofr´e,P., et al. 2016, A&A, 594, A43 Cluver, M. E., Jarrett, T. H., Hopkins, A. M., et al. 2014, ApJ, Helmi, A., Veljanoski, J., Breddels, M. A., Tian, H., & Sales, 782, 90 L. V. 2017, A&A, 598, A58 Comparat, J., Delubac, T., Jouvel, S., et al. 2016, A&A, 592, Hinkle, K., Wallace, L., & Livingston, W. 1995, PASP, 107, 1042 A121 Hinton, S. R., Kazin, E., Davis, T. M., et al. 2017, MNRAS, 464, Dalcanton, J. J., Williams, B. F., Seth, A. C., et al. 2009, ApJS, 4807 183, 67 Ho, I.-T., Medling, A. M., Bland-Hawthorn, J., et al. 2016, Davis, T. A., Alatalo, K., Sarzi, M., et al. 2011, MNRAS, 417, 882 MNRAS, 457, 1257 Dawson, K. S., Schlegel, D. J., Ahn, C. P., et al. 2013, AJ, 145, 10 Holincheck, A. J., Wallin, J. F., Borne, K., et al. 2016, MNRAS, Dawson, K. S., Kneib, J.-P., Percival, W. J., et al. 2016, AJ, 151, 459, 720 44 Holtzman, J. A., Shetrone, M., Johnson, J. A., et al. 2015, AJ, de Medeiros, J. R., Da Rocha, C., & Mayor, M. 1996, A&A, 314, 150, 148 499 Hou, W., Luo, A., Yang, H., et al. 2015, MNRAS, 449, 1401 Data Release 13 29

Howes, L. M., Asplund, M., Casey, A. R., et al. 2014, MNRAS, Piffaretti, R., Arnaud, M., Pratt, G. W., Pointecouteau, E., & 445, 4241 Melin, J.- B. 2011, A&A, 534, A109 Hutchinson, T. A., Bolton, A. S., Dawson, K. S., et al. 2016, AJ, Piffl, T., Scannapieco, C., Binney, J., et al. 2014, A&A, 562, A91 152, 205 Pinsonneault, M. H., Elsworth, Y., Epstein, C., et al. 2014, ApJS, Ibarra-Medel, H. J., S´anchez, S. F., Avila-Reese, V., et al. 2016, 215, 19 MNRAS, 463, 2799 Plez, B. 2012, Astrophysics Source Code Library, ascl:1205.004 Jacobson, H. R., Friel, E. D., J´ılkov´a,L., et al. 2016, A&A, 591, Prakash, A., Licquia, T. C., Newman, J. A., et al. 2016, ApJS, A37 224, 34 Jansen, F., Lumb, D., Altieri, B., et al. 2001, A&A, 365, L1 Predehl, P., Andritschke, R., Becker, W., et al. 2014, in Jarrett, T. H., Masci, F., Tsai, C. W., et al. 2013, AJ, 145, 6 Proc. SPIE, Vol. 9144, Space Telescopes and Instrumentation Jin, Y., Chen, Y., Shi, Y., et al. 2016, MNRAS, 463, 913 2014: Ultraviolet to Gamma Ray, 91441T Johnston, E. J., H¨außler,B., Arag´on-Salamanca,A., et al. 2017, Raichoor, A., Comparat, J., Delubac, T., et al. 2016, A&A, 585, MNRAS, 465, 2317 A50 Kaiser, N., Burgett, W., Chambers, K., et al. 2010, Proc. SPIE, Recio-Blanco, A., Rojas-Arriagada, A., de Laverny, P., et al. 7733, 77330E 2017, A&A, 602, L14 Kaiser, N., Aussel, H., Burke, B. E., et al. 2002, Proc. SPIE, Ren, A., Fu, J., De Cat, P., et al. 2016, ApJS, 225, 28 4836, 154 Robin, A. C., Bienaym´e,O., Fern´andez-Trincado, J. G., & Reyl´e, Kazin, E. A., Koda, J., Blake, C., et al. 2014, MNRAS, 441, 3524 C. 2017, arXiv:1704.06274 Kehrig, C., Monreal-Ibero, A., Papaderos, P., et al. 2012, A&A, Rojas-Arriagada, A., Recio-Blanco, A., de Laverny, P., et al. 540, A11 2017, A&A, 601, A140 Kraft, R. P. 1994, PASP, 106, 553 Ruan, J. J., Anderson, S. F., Green, P. J., et al. 2016, ApJ, 825, Kunder, A., Kordopatis, G., Steinmetz, M., et al. 2017, AJ, 153, 137 75 Ruchti, G. R., Read, J. I., Feltzing, S., et al. 2015, MNRAS, 450, Kunder, A., Bono, G., Piffl, T., et al. 2014, A&A, 572, A30 2874 Lang, D., Hogg, D. W., & Schlegel, D. J. 2016, AJ, 151, 36 Runnoe, J. C., Cales, S., Ruan, J. J., et al. 2016, MNRAS, 455, Law, D. R., Cherinka, B., Yan, R., et al. 2016, AJ, 152, 83 1691 Law, D. R., Yan, R., Bershady, M. A., et al. 2015, The Ryde, N., & Schultheis, M. 2015, A&A, 573, A14 Astronomical Journal, 150, 19 S´anchez, S. F., Rosales-Ortega, F. F., Iglesias-P´aramo,J., et al. Law, N. M., Kulkarni, S. R., Dekany, R. G., et al. 2009, PASP, 2014, A&A, 563, A49 121, 1395 S´anchez, S. F., Garc´ıa-Benito,R., Zibetti, S., et al. 2016, A&A, Li, H., Ge, J., Mao, S., et al. 2017, ApJ, 838, 77 594, A36 Li, C., Wang, E., Lin, L., et al. 2015, ApJ, 804, 125 S´anchez, S. F., Kennicutt, R. C., Gil de Paz, A., et al. 2012, Li, C., Kauffmann, G., Jing, Y. P., et al. 2006, MNRAS, 368, 21 A&A, 538, A8 Li, H.-N., Zhao, G., Christlieb, N., et al. 2015, ApJ, 798, 110 Saxton, R. D., Read, A. M., Esquej, P., et al. 2008, A&A, 480, 611 Li, N., & Thakar, A. R. 2008, Computing in Science and Schaefer, A. L., Croom, S. M., Allen, J. T., et al. 2017, MNRAS, Engineering, 10, 18 464, 121 Licquia, T. C., Newman, J. A., & Brinchmann, J. 2015, ApJ, 809, Schlafly, E. F., Finkbeiner, D. P., Juri´c,M., et al. 2012, ApJ, 756, 96 158 Livingston, W., & Wallace, L. 1991, An atlas of the solar Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 spectrum in the infrared from 1850 to 9000 cm-1 (1.1 to 5.4 Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, micrometer) 525 Magrini, L., Randich, S., Kordopatis, G., et al. 2017, A&A, 603, Sesar, B., Ivezi´c, Z.,ˇ Lupton, R. H., et al. 2007, AJ, 134, 2236 A2 Shetrone, M., Bizyaev, D., Lawler, J. E., et al. 2015, ApJS, 221, Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2015, 24 AJ, in press, arXiv:1509.05420 Shao, L., Kauffmann, G., Li, C., Wang, J., & Heckman, T. M. Martell, S. L., Sharma, S., Buder, S., et al. 2017, MNRAS, 465, 2013, MNRAS, 436, 3451 3203 Shu, Y., Bolton, A. S., Schlegel, D. J., et al. 2012, AJ, 143, 90 Masseron, T., & Gilmore, G. 2015, MNRAS, 453, 1855 Siebert, A., Famaey, B., Minchev, I., et al. 2011, MNRAS, 412, Mel´endez,J., & Barbuy, B. 1999, ApJS, 124, 527 2026 Menzel, M.-L., Merloni, A., Georgakakis, A., et al. 2016, Smee, S. A., Gunn, J. E., Uomoto, A., et al. 2013, AJ, 146, 32 MNRAS, 457, 110 Steinmetz, M., Zwitter, T., Siebert, A., et al. 2006, AJ, 132, 1645 M´esz´aros, S., Martell, S. L., Shetrone, M., et al. 2015, AJ, 149, Smith, M. C., Ruchti, G. R., Helmi, A., et al. 2007, MNRAS, 379, 153 755 M´esz´aros, S., Holtzman, J., Garc´ıaP´erez,A. E., et al. 2013, AJ, Thakar, A. R. 2008, Computing in Science and Engineering, 10, 9 146, 133 Thakar, A. R., Szalay, A., Fekete, G., & Gray, J. 2008, M´esz´aros, S., Allende Prieto, C., Edvardsson, B., et al. 2012, AJ, Computing in Science and Engineering, 10, 30 144, 120 Thomas, N., Ge, J., Grieves, N., Li, R., & Sithajan, S. 2016, Morganson, E., Green, P. J., Anderson, S. F., et al. 2015, ApJ, PASP, 128, 045003 806, 244 Ting, Y.-S., Conroy, C., & Rix, H.-W. 2016, ApJ, 816, 10 Mullaney, J. R., Alexander, D. M., Fine, S., et al. 2013, MNRAS, Vasudevan, R. V., & Fabian, A. C. 2009, MNRAS, 392, 1124 433, 622 Voges, W., Aschenbach, B., Boller, T., et al. 1999, A&A, 349, 389 Myers, A. D., Palanque-Delabrouille, N., Prakash, A., et al. 2015, Wake, D. A., Bundy, K., Diamond-Stanic, A. M., et al. 2017, AJ, ApJS, 221, 27 in press, arXiv:1707.02989 Nidever, D. L., Holtzman, J. A., Allende Prieto, C., et al. 2015, Wallerstein, G. 1962, ApJS, 6, 407 AJ, 150, 173 Williams, M. E. K., Steinmetz, M., Binney, J., et al. 2013, Paegert, M., Stassun, K. G., De Lee, N., et al. 2015, AJ, 149, 186 MNRAS, 436, 101 Palanque-Delabrouille, N., Magneville, C., Y`eche, C., et al. 2016, Wilkinson, D. M., Maraston, C., Thomas, D., et al. 2015, A&A, 587, A41 MNRAS, 449, 328 Palanque-Delabrouille, N., Y`eche, C., Borde, A., et al. 2013, Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, A&A, 559, A85 AJ, 140, 1868-1881 Palanque-Delabrouille, N., Yeche, C., Myers, A. D., et al. 2011, Yan, R., Bundy, K., Law, D. R., et al. 2016, AJ, 152, 197 A&A, 530, A122 Yan, R., Tremonti, C., Bershady, M. A., et al. 2016, AJ, 151, 8 Penny, S. J., Masters, K. L., Weijmans, A.-M., et al. 2016, Yan, R., Tremonti, C., Bershady, M. A., et al. 2015, The MNRAS, 462, 3955 Astronomical Journal, 151, 8 P´erez,E., Cid Fernandes, R., Gonz´alezDelgado, R. M., et al. Yang, X., Mo, H. J., van den Bosch, F. C., et al. 2007, ApJ, 671, 2013, ApJ, 764, L1 153 Yanny, B., Rockosi, C., Newberg, H. J., et al. 2009, AJ, 137, 4377 30 SDSS Collaboration

York, D. G., Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, Zhao, G., Zhao, Y.-H., Chu, Y.-Q., Jing, Y.-P., & Deng, L.-C. 120, 1579 2012, Research in Astronomy and Astrophysics, 12, 723 Zamora, O., Garc´ıa-Hern´andez,D. A., Allende Prieto, C., et al. Zheng, Z., Wang, H., Ge, J., et al. 2017, MNRAS, 465, 4572 2015, AJ, 149, 181 Zhu, G. B., Comparat, J., Kneib, J.-P., et al. 2015, ApJ, 815, 48 Zasowski, G., Johnson, J. A., Frinchaboy, P. M., et al. 2013, AJ, Zou, H., Zhou, X., Jiang, Z., et al. 2016, AJ, 151, 37 146, 81 Zhai, Z., Tinker, J. L., Hahn, C., et al. 2016, arXiv:1607.05383