We characterize the Gaia-Enceladus-Sausage kinematic structure recently discovered in the Galactic halo using photometric metallicities from the SkyMapper survey, and kinematics from Gaia radial velocities measurements. By examining the metallicity distribution functions (MDFs) of stars binned in kinematic/action spaces, we find that the $\sqrt{J_R}$ vs $L_z$ space allows for the cleanest selection of Gaia-Enceladus-Sausage stars with minimal contamination from disc or halo stars formed in situ or in other past mergers. Stars with $30 \leq \sqrt{J_R} \leq 50$ (kpc km s$^{-1})^{1/2}$ and $-500 \leq L_z \leq 500$ kpc km s$^{-1}$ have a narrow MDF centered at [Fe/H] $= -1.17$ dex with a dispersion of 0.34 dex. This [Fe/H] estimate is more metal-rich than literature estimates by $0.1-0.3$ dex. Based on the MDFs, we find that selection of Gaia-Enceladus-Sausage stars in other kinematic/action spaces without additional population information leads to contaminated samples. The clean Gaia-Enceladus-Sausage sample selected according to our criteria is slightly retrograde and lies along the blue sequence of the high $V_T$ halo CMD dual sequence. Using a galaxy mass-metallicity relation derived from cosmological simulations and assuming a mean stellar age of 10 Gyr we estimate the mass of the Gaia-Enceladus-Sausage progenitor satellite to be $10^{8.85-9.85}$ M$_{\odot}$, which is consistent with literature estimates based on disc dynamic and simulations. Additional information on detailed abundances and ages would be needed for a more sophisticated selection of purely Gaia-Enceladus-Sausage stars.

The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,M_\odot$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars.

Spectroscopic surveys of the Milky Way's stars have revealed spatial, chemical and kinematical structures that encode its history. In this work, we study their origins using a cosmological zoom simulation, VINTERGATAN, of a Milky Way-mass disc galaxy. We find that in connection to the last major merger at $z\sim 1.5$, cosmological accretion leads to the rapid formation of an outer, metal-poor, low-[$\alpha$/Fe] gas disc around the inner, metal-rich galaxy containing the old high-[$\alpha$/Fe] stars. This event leads to a bimodality in [$\alpha$/Fe] over a range of [Fe/H]. A detailed analysis of how the galaxy evolves since $z\sim 1$ is presented. We demonstrate the way in which inside-out growth shapes the radial surface density and metallicity profile and how radial migration preferentially relocates stars from the inner to the outer disc. Secular disc heating is found to give rise to increasing velocity dispersions and scaleheights with stellar age, which together with disc flaring explains several trends observed in the Milky Way, including shallower radial [Fe/H]-profiles above the midplane. We show how the galaxy formation scenario imprints non-trivial mappings between structural associations (i.e. thick and thin discs), velocity dispersions, $\alpha$-enhancements, and ages of stars, e.g. the most metal-poor stars in the low-[$\alpha$/Fe] sequence are found to have a scaleheight comparable to old high-[$\alpha$/Fe] stars. Finally, we illustrate how at low spatial resolution, comparable to the thickness of the galaxy, the proposed pathway to distinct sequences in [$\alpha$/Fe]-[Fe/H] cannot be captured.

Using the VINTERGATAN cosmological zoom simulation, we explore the contributions of the in situ and accreted material, and the effect of galaxy interactions and mergers in the assembly of a Milky Way-like galaxy. We find that the initial growth phase of galaxy evolution, dominated by repeated major mergers, provides the necessary physical conditions for the assembly of a thick, kinematically hot disk populated by high-[$\alpha$/Fe] stars, formed both in situ and in accreted satellite galaxies. We find that the diversity of evolutionary tracks followed by the simulated galaxy and its progenitors leads to very little overlap of the in situ and accreted populations for any given chemical composition. At a given age, the spread in [$\alpha$/Fe] abundance ratio results from the diversity of physical conditions in VINTERGATAN and its satellites, with an enhancement in [$\alpha$/Fe] found in stars formed during starburst episodes. Later, the cessation of the merger activity promotes the in situ formation of stars in the low-[$\alpha$/Fe] regime, in a radially extended, thin and overall kinematically colder disk, thus establishing chemically bimodal thin and thick disks, in line with observations. We draw links between notable features in the [Fe/H] - [$\alpha$/Fe] plane with their physical causes, and propose a comprehensive formation scenario explaining self-consistently, in the cosmological context, the main observed properties of the Milky Way.

Using the cosmological zoom simulation VINTERGATAN, we present a new scenario for the onset of star formation at the metal-poor end of the low-[$\alpha$/Fe] sequence in a Milky Way-like galaxy. In this scenario, the galaxy is fueled by two distinct gas flows. One is enriched by outflows from massive galaxies, but not the other. While the former feeds the inner galactic region, the latter fuels an outer gas disk, inclined with respect to the main galactic plane, and with a significantly poorer chemical content. The first passage of the last major merger galaxy triggers tidal compression in the outer disk, which increases the gas density and eventually leads to star formation, at a metallicity 0.75 dex lower than the inner galaxy. This forms the first stars of the low-[$\alpha$/Fe] sequence. These in situ stars have halo-like kinematics, similarly to what is observed in the Milky Way, due to the inclination of the outer disk which eventually aligns with the inner one via gravitational torques. We show that this tilting disk scenario is likely to be common in Milky-Way like galaxies. This process implies that the low-[$\alpha$/Fe] sequence is populated in situ, simultaneously from two formation channels, in the inner and the outer galaxy, with distinct metallicities. This contrasts with purely sequential scenarios for the assembly of the Milky Way disk and could be tested observationally.

In the era of large spectroscopic surveys, massive databases of high-quality spectra provide tools to outline a new picture of our Galaxy. In this framework, an important piece of information is provided by our ability to infer stellar ages. We aim to provide empirical relations between stellar ages and abundance ratios for a sample of solar-like stars. We investigate the dependence on metallicity, and we apply our relations to Gaia-ESO samples of open clusters and field stars. We analyse high-resolution and high-S/N HARPS spectra of a sample of solar-like stars to obtain precise determinations of their atmospheric parameters and abundances through differential spectral analysis and age through isochrone fitting. We investigate the relations between ages and abundance ratios. For the abundance ratios with a steeper dependence on age, we perform multivariate linear regressions, including the dependence on metallicity. We apply our relations to a sample of open clusters located in 4<R$_{GC}$<16 kpc. Using them, we are able to recover the literature ages only for clusters located at R$_{GC}$>7 kpc. In these clusters, the content of s-elements is lower than expected from chemical evolution models, and consequently the [s/$\alpha$] are lower than in clusters of the same age located in the solar neighbourhood. With our chemical evolution model and a set of empirical yields, we suggest that a strong dependence on the star formation history and metallicity-dependent yields of s-elements can substantially modify the slope of the [s/$\alpha$]-[Fe/H]-age relation in different regions of the Galaxy. Our results point towards a non-universal relation [s/$\alpha$]-[Fe/H]-age, indicating the existence of relations at different R$_{GC}$ or for different star formation history. A better understanding of the s-process at high metallicity is necessary to fully understand the origin of these variations.

The HR 1614 overdensity in velocity space and has for a long time been known as an old (~2 Gyr) and metal-rich ([Fe/H]~0.2) nearby moving group that has a dissolving cluster origin. The existence of such old and metal-rich groups in the solar vicinity is quite unexpected since the vast majority of nearby moving groups are known to be young. In the light of new and significantly larger data sets we aim to re-investigate the properties and origin of the HR 1614 moving group. To identify and characterise the HR 1614 moving group we use astrometric data from Gaia DR2; distances, extinction, and reddening corrections from the StarHorse code; elemental abundances from the GALAH and APOGEE spectroscopic surveys; and photometric metallicities from the SkyMapper survey. Bayesian ages were estimated for the SkyMapper stars. Since the Hercules stream is the closest kinematical structure to the HR 1614 moving group in velocity space, we use it for comparison purposes. Stars that are likely to be members of the two groups were selected based on their space velocities. The HR 1614 moving group is located mainly at negative U velocities, does not form an arch of constant energy in the U-V space and is tilted in V. The overdensity is not chemically homogeneous but that its stars exist at a wide range of both metallicities, ages, and elemental abundance ratios. They are essentially similar to what is observed in the Galactic thin and thick disks, a younger population (~3 Gyr) that is metal-rich (-0.2<[Fe/H]<0.4) and alpha-poor. It should therefore not be considered as a dissolving open cluster, or an accreted population. We suggest that HR 1614 has a complex origin that could be explained by combining several different mechanisms such as resonances with the Galactic bar and spiral structure, phase-mixing of dissolving spiral structure, and phase-mixing due to an external perturbation.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

Recent and on-going large ground-based multi-object spectroscopic surveys allow to significantly increase the sample of spectroscopic binaries to get insight into their statistical properties. We investigate the repeated spectral observations of the Gaia-ESO Survey (GES) internal data release 5 to identify and characterize spectroscopic binaries with one visible component (SB1) in fields covering the discs, the bulge, the CoRot fields, and stellar clusters and associations. A statistical chi2-test is performed on spectra of the iDR5 sub-sample of approximately 43500 stars characterized by at least 2 observations and a S/N > 3. Our sample of RV variables is cleaned from contamination by pulsation/convection-induced variables using Gaia DR2 parallaxes and photometry. Monte-Carlo simulations using the SB9 catalogue of spectroscopic orbits allow to estimate our detection efficiency and to correct the SB1 rate to evaluate the GES SB1 binary fraction and its dependence with effective temperature and metallicity. We find 641 (resp., 803) FGK SB1 candidates at the 5 sigma (resp., 3 sigma) level. The orbital-period distribution is estimated from the RV standard-deviation distribution and reveals that the detected SB1 probe binaries with log(P[d]) < 4. We estimate the global GES SB1 fraction to be in the range 7-14% with a typical uncertainty of 4%. The GES SB1 frequency decreases with metallicity at a rate of -9+/-3%/dex in the metallicity range -2.7<FeH<+0.6. This anticorrelation is obtained with a confidence level higher than 93% on a homogeneous sample covering spectral types FGK and a large range of metallicities. When the present-day mass function is accounted for, this rate turns to 4+/-2%/dex with a confidence level higher than 88%. In addition we provide the variation of the SB1 fraction with metallicity separately for F, G, and K spectral types, as well as for dwarf and giant primaries.

Analyses of data from spectroscopic and astrometric surveys have led to conflicting results concerning the vertical characteristics of the Milky Way. Ages are often used to provide clarity, but typical uncertainties of $>$ 40\,\% restrict the validity of the inferences made. Using the \textit{Kepler} APOKASC sample for context, we explore the global population trends of two K2 campaign fields (3 and 6), which extend further vertically out of the Galactic plane than APOKASC. We analyse the properties of red giant stars utilising three asteroseismic data analysis methods to cross-check and validate detections. The Bayesian inference tool PARAM is used to determine the stellar masses, radii and ages. Evidence of a pronounced red giant branch bump and an [$\alpha$/Fe] dependence on the position of the red clump is observed from the radii distribution of the K2 fields. Two peaks in the age distribution centred at $\sim$5 and and $\sim$12 Gyr are found using a sample with $\sigma_{\rm{age}}$ $<$ 35\,\%. In a comparison with \textit{Kepler}, we find the older peak to be more prominent for K2. This age bimodality is also observed based on a chemical selection of low- ($\leq$ 0.1) and high- ($>$ 0.1) [$\alpha$/Fe] stars. As a function of vertical distance from the Galactic mid-plane ($|Z|$), the age distribution shows a transition from a young to old stellar population with increasing $|Z|$ for the K2 fields. Further coverage of campaign targets with high resolution spectroscopy is required to increase the yield of precise ages achievable with asteroseismology.

Stellar ages are a crucial component to studying the evolution of the Milky Way. Using Gaia DR2 distance estimates, it is now possible to estimate stellar ages for a larger volume of evolved stars through isochrone matching. This work presents [M/H]-age and [$\alpha$/M]-age relations derived for different spatial locations in the Milky Way disc. These relations are derived by hierarchically modelling the star formation history of stars within a given chemical abundance bin. For the first time, we directly observe that significant variation is apparent in the [M/H]-age relation as a function of both Galactocentric radius and distance from the disc mid-plane. The [M/H]-age relations support claims that radial migration has a significant effect in the plane of the disc. Using the [M/H] bin with the youngest mean age at each radial zone in the plane of the disc, the present-day metallicity gradient is measured to be $-0.059 \pm 0.010$ dex kpc$^{-1}$, in agreement with Cepheids and young field stars. We find a vertically flared distribution of young stars in the outer disc, confirming predictions of models and previous observations. The mean age of the [M/H]-[$\alpha$/M] distribution of the solar neighborhood suggests that the high-[M/H] stars are not an evolutionary extension of the low-$\alpha$ sequence. Our observational results are important constraints to Galactic simulations and models of chemical evolution.

The kinematics of the Milky Way disc as a function of age are well measured at the solar radius, but have not been studied over a wider range of Galactocentric radii. Here, we measure the kinematics of mono-age, mono-$\mathrm{[Fe/H]}$ populations in the low and high $\mathrm{[\alpha/Fe]}$ discs between $4 \lesssim R \lesssim 13$ kpc and $|z| \lesssim 2$ kpc using 65,719 stars in common between APOGEE DR14 and $\it{Gaia}$ DR2 for which we estimate ages using a Bayesian neural network model trained on asteroseismic ages. We determine the vertical and radial velocity dispersions, finding that the low and high $\mathrm{[\alpha/Fe]}$ discs display markedly different age--velocity-dispersion relations (AVRs) and shapes $\sigma_z/\sigma_R$. The high $\mathrm{[\alpha/Fe]}$ disc has roughly flat AVRs and constant $\sigma_z/\sigma_R = 0.64\pm 0.04$, whereas the low $\mathrm{[\alpha/Fe]}$ disc has large variations in this ratio which positively correlate with the mean orbital radius of the population at fixed age. The high $\mathrm{[\alpha/Fe]}$ disc component's flat AVRs and constant $\sigma_z/\sigma_R$ clearly indicates an entirely different heating history. Outer disc populations also have flatter radial AVRs than those in the inner disc, likely due to the waning effect of spiral arms. Our detailed measurements of AVRs and $\sigma_z/\sigma_R$ across the disc indicate that low $\mathrm{[\alpha/Fe]}$, inner disc ($R \lesssim 10\,\mathrm{kpc}$) stellar populations are likely dynamically heated by both giant molecular clouds and spiral arms, while the observed trends for outer disc populations require a significant contribution from another heating mechanism such as satellite perturbations. We also find that outer disc populations have slightly positive mean vertical and radial velocities, likely because they are part of the warped disc.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

We present [C/N]-[Fe/H] abundance trends from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey, Data Release 14 (DR14), for red giant branch stars across the Milky Way Galaxy (MW, 3 kpc $<$ R $<$ 15 kpc). The carbon-to-nitrogen ratio (often expressed as [C/N]) can indicate the mass of a red giant star, from which an age can be inferred. Using masses and ages derived by Martig et al., we demonstrate that we are able to interpret the DR14 [C/N]-[Fe/H] abundance distributions as trends in age-[Fe/H] space. Our results show that an anti-correlation between age and metallicity, which is predicted by simple chemical evolution models, is not present at any Galactic zone. Stars far from the plane ($|$Z$|$ $>$ 1 kpc) exhibit a radial gradient in [C/N] ($\sim$ $-$0.04 dex/kpc). The [C/N] dispersion increases toward the plane ($\sigma_{[C/N]}$ = 0.13 at $|$Z$|$ $>$ 1 kpc to $\sigma_{[C/N]}$ = 0.18 dex at $|$Z$|$ $<$ 0.5 kpc). We measure a disk metallicity gradient for the youngest stars (age $<$ 2.5 Gyr) of $-$0.060 dex/kpc from 6 kpc to 12 kpc, which is in agreement with the gradient found using young CoRoGEE stars by Anders et al. Older stars exhibit a flatter gradient ($-$0.016 dex/kpc), which is predicted by simulations in which stars migrate from their birth radii. We also find that radial migration is a plausible explanation for the observed upturn of the [C/N]-[Fe/H] abundance trends in the outer Galaxy, where the metal-rich stars are relatively enhanced in [C/N].

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

Data from the SDSS-IV / Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) have been released as part of SDSS Data Releases 13 (DR13) and 14 (DR14). These include high resolution H-band spectra, radial velocities, and derived stellar parameters and abundances. DR13, released in August 2016, contained APOGEE data for roughly 150,000 stars, and DR14, released in August 2017, added about 110,000 more. Stellar parameters and abundances have been derived with an automated pipeline, the APOGEE Stellar Parameter and Chemical Abundance Pipeline (ASPCAP). We evaluate the performance of this pipeline by comparing the derived stellar parameters and abundances to those inferred from optical spectra and analysis for several hundred stars. For most elements -- C, Na, Mg, Al, Si, S, Ca, Cr, Mn, Ni -- the DR14 ASPCAP analysis have systematic differences with the comparisons samples of less than 0.05 dex (median), and random differences of less than 0.15 dex (standard deviation). These differences are a combination of the uncertainties in both the comparison samples as well as the ASPCAP-analysis. Compared to the references, magnesium is the most accurate alpha-element derived by ASPCAP, and shows a very clear thin/thick disk separation, while nickel is the most accurate iron-peak element (besides iron).

We report 20 new lithium-rich giants discovered within the Gaia-ESO Survey, including the first Li-rich giant with evolutionary stage confirmed by CoRoT data. Atmospheric parameters and abundances were derived in model atmosphere analyses using medium-resolution GIRAFFE or high-resolution UVES spectra. These results are part of the fifth internal data release of Gaia-ESO. The Li abundances were corrected for non-LTE effects. We used Gaia DR2 parallaxes to estimate distances and luminosities. The giants have A(Li) > 2.2 dex. The majority of them (14 out of 20 stars) are in the CoRoT fields. Four giants are located in the field of three open clusters but are not members. Two giants were observed in fields towards the Galactic bulge but are likely in the inner disk. One of the bulge field giants is super Li-rich with A(Li) = 4.0 dex. We identified one giant with infrared excess at 22 microns. Two other giants, with large vsin i, might be Li-rich because of planet engulfment. Another giant is found to be barium enhanced and thus could have accreted material from a former AGB companion. Otherwise, besides the Li enrichment, the evolutionary stages are the only other connection between these new Li-rich giants. The CoRoT data confirm that one Li-rich giant is at the core-He burning stage. The other giants are concentrated in close proximity to the RGB luminosity bump, the core-He burning stages, or the early-AGB. This is very clear when looking at the Gaia-based luminosities of the Li-rich giants. This is also seen when the CoRoT Li-rich giants are compared to a larger sample of 2252 giants observed in the CoRoT fields by the Gaia-ESO Survey, which are distributed all over the RGB in the Teff-logg diagram. These observations show that evolutionary stage is a major factor behind the Li enrichment in giants. Other processes, like planet accretion, contribute to a smaller scale. [abridged]

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

In the age of high-resolution spectroscopic stellar surveys of the Milky Way, the number of stars with detailed abundances of multiple elements is rapidly increasing. These elemental abundances are directly influenced by the evolutionary history of the Galaxy, but this can be difficult to interpret without an absolute timeline of the abundance enrichment. We present age-abundance trends for [M/H], [{\alpha}/M], and 17 individual elements using a sample of 721 solar neighbourhood Hipparcos red giant stars observed by APOGEE. These age trends are determined through a Bayesian hierarchical modelling method presented by Feuillet et al. (2016). We confirm that the [{\alpha}/M]- age relation in the solar neighbourhood is steep and relatively narrow (0.20 dex age dispersion), as are the [O/M]- and [Mg/M]-age relations. The age trend of [C/N] is steep and smooth, consistent with stellar evolution. The [M/H]-age relation has a mean age dispersion of 0.28 dex and a complex overall structure. The oldest stars in our sample are those with the lowest and highest metallicities, while the youngest stars are those with solar metallicity. These results provide strong constraints on theoretical models of Galactic chemical evolution (GCE). We compare them to the predictions of one-zone GCE mod- els and multi-zone mixtures, both analytic and numerical. These comparisons support the hypothesis that the solar neighbourhood is composed of stars born at a range of Galactocentric radii, and that the most metal-rich stars likely migrated from a region with earlier and more rapid star formation such as the inner Galaxy.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

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

We present a sample of local red giant stars observed using the New Mexico State University 1 m telescope with the APOGEE spectrograph, for which we estimate stellar ages and the age distribution from the high-resolution spectroscopic stellar parameters and accurate distance measurements from Hipparcos. The high-resolution (R ~ 23,000), near infrared (H-band, 1.5-1.7 micron) APOGEE spectra provide measurements of the stellar atmospheric parameters (temperature, surface gravity, [M/H], and [alpha/M]). Due to the smaller uncertainties in surface gravity possible with high-resolution spectra and accurate Hipparcos distance measurements, we are able to calculate the stellar masses to within 40%. For red giants, the relatively rapid evolution of stars up the red giant branch allows the age to be constrained based on the mass. We examine methods of estimating age using both the mass-age relation directly and a Bayesian isochrone matching of measured parameters, assuming a constant star formation history (SFH). To improve the prior on the SFH, we use a hierarchical modeling approach to constrain the parameters of a model SFH from the age probability distribution functions of the data. The results of an alpha dependent Gaussian SFH model shows a clear relation between age and [alpha/M] at all ages. Using this SFH model as the prior for an empirical Bayesian analysis, we construct a full age probability distribution function and determine ages for individual stars. The age-metallicity relation is flat, with a slight decrease in [M/H] at the oldest ages and a ~ 0.5 dex spread in metallicity. For stars with ages < 1 Gyr we find a smaller spread, consistent with radial migration having a smaller effect on these young stars than on the older stars.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) has built the largest moderately high-resolution (R=22, 500) spectroscopic map of the stars across the Milky Way, and including dust-obscured areas. The APOGEE Stellar Parameter and Chemical Abundances Pipeline (ASPCAP) is the software developed for the automated analysis of these spectra. ASPCAP determines atmospheric parameters and chemical abundances from observed spectra by comparing observed spectra to libraries of theoretical spectra, using chi-2 minimization in a multidimensional parameter space. The package consists of a fortran90 code that does the actual minimization, and a wrapper IDL code for book-keeping and data handling. This paper explains in detail the ASPCAP components and functionality, and presents results from a number of tests designed to check its performance. ASPCAP provides stellar effective temperatures, surface gravities, and metallicities precise to 2%, 0.1 dex, and 0.05 dex, respectively, for most APOGEE stars, which are predominantly giants. It also provides abundances for up to 15 chemical elements with various levels of precision, typically under 0.1 dex. The final data release (DR12) of the Sloan Digital Sky Survey III contains an APOGEE database of more than 150,000 stars. ASPCAP development continues in the SDSS-IV APOGEE-2 survey.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE), one of the programs in the Sloan Digital Sky Survey III (SDSS-III), has now completed its systematic, homogeneous spectroscopic survey sampling all major populations of the Milky Way. After a three year observing campaign on the Sloan 2.5-m Telescope, APOGEE has collected a half million high resolution (R~22,500), high S/N (>100), infrared (1.51-1.70 microns) spectra for 146,000 stars, with time series information via repeat visits to most of these stars. This paper describes the motivations for the survey and its overall design---hardware, field placement, target selection, operations---and gives an overview of these aspects as well as the data reduction, analysis and products. An index is also given to the complement of technical papers that describe various critical survey components in detail. Finally, we discuss the achieved survey performance and illustrate the variety of potential uses of the data products by way of a number of science demonstrations, which span from time series analysis of stellar spectral variations and radial velocity variations from stellar companions, to spatial maps of kinematics, metallicity and abundance patterns across the Galaxy and as a function of age, to new views of the interstellar medium, the chemistry of star clusters, and the discovery of rare stellar species. As part of SDSS-III Data Release 12, all of the APOGEE data products are now publicly available.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

The SDSS-III/APOGEE survey operated from 2011-2014 using the APOGEE spectrograph, which collects high-resolution (R~22,500), near-IR (1.51-1.70 microns) spectra with a multiplexing (300 fiber-fed objects) capability. We describe the survey data products that are publicly available, which include catalogs with radial velocity, stellar parameters, and 15 elemental abundances for over 150,000 stars, as well as the more than 500,000 spectra from which these quantities are derived. Calibration relations for the stellar parameters (Teff, log g, [M/H], [alpha/M]) and abundances (C, N, O, Na, Mg, Al, Si, S, K, Ca, Ti, V, Mn, Fe, Ni) are presented and discussed. The internal scatter of the abundances within clusters indicates that abundance precision is generally between 0.05 and 0.09 dex across a broad temperature range; within more limited ranges and at high S/N, it is smaller for some elemental abundances. We assess the accuracy of the abundances using comparison of mean cluster metallicities with literature values, APOGEE observations of the solar spectrum and of Arcturus, comparison of individual star abundances with other measurements, and consideration of the locus of derived parameters and abundances of the entire sample, and find that it is challenging to determine the absolute abundance scale; external accuracy may be good to 0.1-0.2 dex. Uncertainties may be larger at cooler temperatures (Teff<4000K). Access to the public data release and data products is described, and some guidance for using the data products is provided.

We employ the first two years of data from the near-infrared, high-resolution SDSS-III/APOGEE spectroscopic survey to investigate the distribution of metallicity and alpha-element abundances of stars over a large part of the Milky Way disk. Using a sample of ~10,000 kinematically-unbiased red-clump stars with ~5% distance accuracy as tracers, the [alpha/Fe] vs. [Fe/H] distribution of this sample exhibits a bimodality in [alpha/Fe] at intermediate metallicities, -0.9<[Fe/H]<-0.2, but at higher metallicities ([Fe/H]=+0.2) the two sequences smoothly merge. We investigate the effects of the APOGEE selection function and volume filling fraction and find that these have little qualitative impact on the alpha-element abundance patterns. The described abundance pattern is found throughout the range 5<R<11 kpc and 0<|Z|<2 kpc across the Galaxy. The [alpha/Fe] trend of the high-alpha sequence is surprisingly constant throughout the Galaxy, with little variation from region to region (~10%). Using simple galactic chemical evolution models we derive an average star formation efficiency (SFE) in the high-alpha sequence of ~4.5E-10 1/yr, which is quite close to the nearly-constant value found in molecular-gas-dominated regions of nearby spirals. This result suggests that the early evolution of the Milky Way disk was characterized by stars that shared a similar star formation history and were formed in a well-mixed, turbulent, and molecular-dominated ISM with a gas consumption timescale (1/SFE) of ~2 Gyr. Finally, while the two alpha-element sequences in the inner Galaxy can be explained by a single chemical evolutionary track this cannot hold in the outer Galaxy, requiring instead a mix of two or more populations with distinct enrichment histories.

We present new $BVI$ photometry containing 40,000 stars of the Galactic globular cluster M15 (NGC 7078) covering a $25 \times 25$ arcminute area centered on the cluster with a magnitude range from the tip of the Red Giant Branch to three magnitudes below the main sequence turn-off. Using $\alpha$-enhanced Dartmouth Stellar Evolution Program models, we find an age of $13.0 \ \pm \ 1.0$ Gyr and distance modulus of $(m-M)_V = 15.4 \pm 0.1$ through isochrone fitting. Unlike previous works, we find good agreement between the observed completeness-corrected stellar luminosity function and models.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) - one of the Sloan Digital Sky Survey III programs -- is using near-infrared spectra of $\sim 100,000$ red giant branch star candidates to study the structure of the Milky Way. In the course of the survey, APOGEE also acquires spectra of hot field stars to serve as telluric calibrators for the primary science targets. We report the serendipitous discovery of two rare, fast-rotating B stars of the $\sigma$ Ori E type among those blue field stars observed during the first year of APOGEE operations. Both of the discovered stars display the spectroscopic signatures of the rigidly rotating magnetospheres (RRM) common to this class of highly-magnetized ($B \sim 10$ kiloGauss) stars, increasing the number of known RRM stars by $\sim 10 \%$. One (HD 345439) is a main-sequence B star with unusually strong He absorption (similar to $\sigma$ Ori E), while the other (HD 23478) fits a "He-normal" B3IV classification. We combine the APOGEE discovery spectra with other optical and near-infrared spectra of these two stars, and of $\sigma$ Ori E itself, to show how near-infrared spectroscopy can be a uniquely powerful tool for discovering more of these rare objects, which may show little/no RRM signatures in their optical spectra. We discuss the potential for further discovery of $\sigma$ Ori E type stars, as well as the implications of our discoveries for the population of these objects and insights into their origin and evolution.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We present Galactic mean metallicity maps derived from the first year of the SDSS-III APOGEE experiment. Mean abundances in different zones of Galactocentric radius (0 < R < 15 kpc) at a range of heights above the plane (0 < |z| < 3 kpc), are derived from a sample of nearly 20,000 stars with unprecedented coverage, including stars in the Galactic mid-plane at large distances. We also split the sample into subsamples of stars with low and high-[{\alpha}/M] abundance ratios. We assess possible biases in deriving the mean abundances, and find they are likely to be small except in the inner regions of the Galaxy. A negative radial gradient exists over much of the Galaxy; however, the gradient appears to flatten for R < 6 kpc, in particular near the Galactic mid-plane and for low-[{\alpha}/M] stars. At R > 6 kpc, the gradient flattens as one moves off of the plane, and is flatter at all heights for high-[{\alpha}/M] stars than for low-[{\alpha}/M] stars. Alternatively, these gradients can be described as vertical gradients that flatten at larger Galactocentric radius; these vertical gradients are similar for both low and high-[{\alpha}/M] populations. Stars with higher [{\alpha}/M] appear to have a flatter radial gradient than stars with lower [{\alpha}/M]. This could suggest that the metallicity gradient has grown steeper with time or, alternatively, that gradients are washed out over time by migration of stars.