We investigate the environmental dependence of the stellar populations of galaxies in SDSS DR7. Echoing earlier works, we find that satellites are both more metal-rich (<0.1 dex) and older (<2 Gyr) than centrals of the same stellar mass. However, after separating star-forming, green valley and passive galaxies, we find that the true environmental dependence of both stellar metallicity (<0.03 dex) and age (<0.5 Gyr) is in fact much weaker. We show that the strong environmental effects found when galaxies are not differentiated result from a combination of selection effects brought about by the environmental dependence of the quenched fraction of galaxies, and thus we strongly advocate for the separation of star-forming, green valley and passive galaxies when the environmental dependence of galaxy properties are investigated. We also study further environmental trends separately for both central and satellite galaxies. We find that star-forming galaxies show no environmental effects, neither for centrals nor for satellites. In contrast, the stellar metallicities of passive and green valley satellites increase weakly with increasing halo mass, increasing local overdensity and decreasing projected distance from their central; this effect is interpreted in terms of moderate environmental starvation (`strangulation') contributing to the quenching of satellite galaxies. Finally, we find a unique feature in the stellar mass--stellar metallicity relation for passive centrals, where galaxies in more massive haloes have larger stellar mass at constant stellar metallicity; this effect is interpreted in terms of dry merging of passive central galaxies and/or progenitor bias.

Star-forming galaxies can in principle be transformed into passive systems by a multitude of processes that quench star formation, such as the halting of gas accretion (starvation) or the rapid removal of gas in AGN-driven outflows. However, it remains unclear which processes are the most significant, primary drivers of the SF-passive bimodality. We address this key issue in galaxy evolution by studying the chemical properties of 80,000 local galaxies in SDSS DR7. In order to distinguish between different quenching mechanisms, we analyse the stellar metallicities of star-forming, green valley and passive galaxies. We find that the significant difference in stellar metallicity between passive galaxies and their star-forming progenitors implies that for galaxies at all masses, quenching must have involved an extended phase of starvation. However, some form of gas ejection also has to be introduced into our models to best match the observed properties of local passive galaxies, indicating that, while starvation is likely to be the prerequisite for quenching, it is the combination of starvation and outflows that is responsible for quenching the majority of galaxies. Closed-box models indicate that the duration of the quenching phase is 2-3 Gyr, with an $e$-folding time of 2-4 Gyr, after which further star formation is prevented by an ejective/heating mode. Alternatively, leaky-box models find a longer duration for the quenching phase of 5-7 Gyr and an $e$-folding time of $\sim$1 Gyr, with outflows becoming increasingly important with decreasing stellar mass. Finally, our analysis of local green valley galaxies indicates that quenching is slower in the local Universe than at high-redshift.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is acquiring integral-field spectroscopy for the largest sample of galaxies to date. By 2020, the MaNGA Survey --- one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) --- will have observed a statistically representative sample of 10$^4$ galaxies in the local Universe ($z\lesssim0.15$). In addition to a robust data-reduction pipeline (DRP), MaNGA has developed a data-analysis pipeline (DAP) that provides higher-level data products. To accompany the first public release of its code base and data products, we provide an overview of the MaNGA DAP, including its software design, workflow, measurement procedures and algorithms, performance, and output data model. In conjunction with our companion paper Belfiore et al., we also assess the DAP output provided for 4718 observations of 4648 unique galaxies in the recent SDSS Data Release 15 (DR15). These analysis products focus on measurements that are close to the data and require minimal model-based assumptions. Namely, we provide stellar kinematics (velocity and velocity dispersion), emission-line properties (kinematics, fluxes, and equivalent widths), and spectral indices (e.g., D4000 and the Lick indices). We find that the DAP provides robust measurements and errors for the vast majority ($>$99%) of analyzed spectra. We summarize assessments of the precision and accuracy of our measurements as a function of signal-to-noise, and provide specific guidance to users regarding the limitations of the data. The MaNGA DAP software is publicly available and we encourage community involvement in its development.

We determine the stellar population properties - age, metallicity, dust reddening, stellar mass and the star formation history - for all spectra classified as galaxies that were published by the Sloan Digital Sky Survey (SDSS data release 14) and by the DEEP2 (data release 4) galaxy surveys. We perform full spectral fitting on individual spectra, making use of high spectral resolution stellar population models. Calculations are carried out for several choices of the model input, including three stellar initial mass functions and three input stellar libraries to the models. We study the accuracy of parameter derivation, in particular the stellar mass, as a function of the signal-to-noise of the galaxy spectra. We find that at low redshift, a signal to noise ratio per pixel around 20 (5) allows a statistical accuracy on $\log_{10}(M^{*}/M_{\odot})$ of 0.2 (0.4) dex, for the Chabrier IMF. For the first time, we study DEEP2 galaxies selected by their \OII luminosity in the redshift range $0.83<z<1.03$, finding that they are consistent with a flat number density in stellar mass in the range $10^9<M/M_{\odot}<10^{11.5}$. We find the resulting stellar mass function based on SDSS or eBOSS in agreement with previous studies (Maraston et al. 2013). We publish all catalogs of properties as well as model spectra of the continuum for these galaxies as a value added catalog of the fourteenth data release of the SDSS. This catalog is about twice as large as its predecessors (DR12) and will aid a variety of studies on galaxy evolution and cosmology.

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.

MaNGA provides the opportunity to make precise spatially resolved measurements of the IMF slope in galaxies owing to its unique combination of spatial resolution, wavelength coverage and sample size. We derive radial gradients in age, element abundances and IMF slope analysing optical and near-infrared absorption features from stacked spectra out to the half-light radius of 366 early-type galaxies with masses $9.9 - 10.8\;\log M/M_{\odot}$. We find flat gradients in age and [$\alpha$/Fe] ratio, as well as negative gradients in metallicity, consistent with the literature. We further derive significant negative gradients in the [Na/Fe] ratio with galaxy centres being well enhanced in Na abundance by up to 0.5 dex. Finally, we find a gradient in IMF slope with a bottom-heavy IMF in the centre (typical mass excess factor of 1.5) and a Milky Way-type IMF at the half-light radius. This pattern is mass-dependent with the lowest mass galaxies in our sample featuring only a shallow gradient around a Milky Way IMF. Our results imply the local IMF-$\sigma$ relation within galaxies to be even steeper than the global relation and hint towards the local metallicity being the dominating factor behind the IMF variations. We also employ different stellar population models in our analysis and show that a radial IMF gradient is found independently of the stellar population model used. A similar analysis of the Wing-Ford band provides inconsistent results and further evidence of the difficulty in measuring and modelling this particular feature.

In our previous work, we found that only two scenarios are capable of reproducing the observed integrated mass-metallicity relations for the gas and stellar components of local star-forming galaxies simultaneously. One scenario invokes a time-dependent metal outflow loading factor with stronger outflows at early times. The other scenario uses a time-dependent IMF slope with a steeper IMF at early times. In this work, we extend our study to investigate the radial profile of gas and stellar metallicity in local star-forming galaxies using spatially resolved spectroscopic data from the SDSS-IV MaNGA survey. We find that most galaxies show negative gradients in both gas and stellar metallicity with steeper gradients in stellar metallicity. The stellar metallicity gradient tend to be mass dependent with steeper gradients in more massive galaxies while no clear mass dependence is found for the gas metallicity gradient. Then we compare the observations with the predictions from a chemical evolution model of the radial profiles of gas and stellar metallicities. We confirm that the two scenarios proposed in our previous work are also required to explain the metallicity gradients. Based on these two scenarios we successfully reproduce the radial profiles of gas metallicity, stellar metallicity, stellar mass surface density, and star formation rate (SFR) surface density simultaneously. The origin of the negative gradient in stellar metallicity turns out to be driven by either radially dependent metal outflow or IMF slope. In contrast, the radial dependence of the gas metallicity is less constrained because of the degeneracy in model parameters.

We present a new spectral fitting code, Firefly, for deriving the stellar population properties of stellar systems. Firefly is a chi-squared minimisation fitting code that fits combinations of single-burst stellar population models to spectroscopic data, following an iterative best-fitting process controlled by the Bayesian Information Criterion. No priors are applied, rather all solutions within a statistical cut are retained with their weight. Moreover, no additive or multiplicative polynomia are employed to adjust the spectral shape. This fitting freedom is envisaged in order to map out the effect of intrinsic spectral energy distribution (SED) degeneracies, such as age, metallicity, dust reddening on galaxy properties, and to quantify the effect of varying input model components on such properties. Dust attenuation is included using a new procedure, which was tested on Integral Field Spectroscopic (IFS) data in a previous paper. The fitting method is extensively tested with a comprehensive suite of mock galaxies, real galaxies from the Sloan Digital Sky Survey and Milky Way globular clusters. We also assess the robustness of the derived properties as a function of signal-to-noise ratio and adopted wavelength range. We show that \FF\ is able to recover age, metallicity, stellar mass and even the star formation history remarkably well down to a $S/N\sim5$, for moderately dusty systems. Code and results are publicly available at www.icg.port.ac.uk/firefly.

We investigate the mass-metallicity relations for the gaseous (MZRgas) and stellar components (MZRstar) of local star-forming galaxies based on a representative sample from SDSS DR12. The mass-weighted average stellar metallicities are systematically lower than the gas metallicities. This difference in metallicity increases toward galaxies with lower masses and reaches 0.4-0.8 dex at 10^9 Msun (depending on the gas metallicity calibration). As a result, the MZRstar is much steeper than the MZRgas. The much lower metallicities in stars compared to the gas in low mass galaxies implies dramatic metallicity evolution with suppressed metal enrichment at early times. The aim of this paper is to explain the observed large difference in gas and stellar metallicity and to infer the origin of the mass-metallicity relations. To this end we develop a galactic chemical evolution model accounting for star formation, gas inflow and outflow. By combining the observed mass-metallicity relation for both gas and stellar components to constrain the models, we find that only two scenarios are able to reproduce the observations. Either strong metal outflow or a steep IMF slope at early epochs of galaxy evolution is needed. Based on these two scenarios, for the first time we successfully reproduce the observed MZRgas and MZRstar simultaneously, together with other independent observational constraints in the local universe. Our model also naturally reproduces the flattening of the MZRgas at the high mass end leaving the MZRstar intact, as seen in observational data.

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

The "kinematic" morphology-density relation for early-type galaxies posits that those galaxies with low angular momentum are preferentially found in the highest-density regions of the universe. We use a large sample of galaxy groups with halo masses 10^12.5 < M_halo < 10^14.5 M_sun/h observed with the Mapping Nearby Galaxies at APO (MaNGA) survey to examine whether there is a correlation between local environment and rotational support that is independent of stellar mass. We find no compelling evidence for a relationship between the angular momentum content of early-type galaxies and either local overdensity or radial position within the group at fixed stellar mass.

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 study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

We study the internal radial gradients of stellar population properties within $1.5\;R_{\rm e}$ and analyse the impact of galaxy environment. We use a representative sample of 721 galaxies with masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV survey MaNGA. We split this sample by morphology into early-type and late-type galaxies. Using the full spectral fitting code FIREFLY, we derive the light and mass-weighted stellar population properties age and metallicity, and calculate the gradients of these properties. We use three independent methods to quantify galaxy environment, namely the $N^{th}$ nearest neighbour, the tidal strength parameter $Q$ and distinguish between central and satellite galaxies. In our analysis, we find that early-type galaxies generally exhibit shallow light-weighted age gradients in agreement with the literature and mass-weighted median age gradients tend to be slightly positive. Late-type galaxies, instead, have negative light-weighted age gradients. We detect negative metallicity gradients in both early and late-type galaxies that correlate with galaxy mass, with the gradients being steeper and the correlation with mass being stronger in late-types. We find, however, that stellar population gradients, for both morphological classifications, have no significant correlation with galaxy environment for all three characterisations of environment. Our results suggest that galaxy mass is the main driver of stellar population gradients in both early and late-type galaxies, and any environmental dependence, if present at all, must be very subtle.

We study the internal gradients of stellar population properties within $1.5\;R_{\rm e}$ for a representative sample of 721 galaxies with stellar masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV MaNGA IFU survey. Through the use of our full spectral fitting code FIREFLY, we derive light and mass-weighted stellar population properties and their radial gradients, as well as full star formation and metal enrichment histories. We also quanfify the impact that different stellar population models and full spectral fitting routines have on the derived stellar population properties, and the radial gradient measurements. In our analysis, we find that age gradients tend to be shallow for both early-type and late-type galaxies. {\em Mass-weighted} age gradients of early-types are positive ($\sim 0.09\; {\rm dex}/R_{\rm e}$) pointing to "outside-in" progression of star formation, while late-type galaxies have negative {\em light-weighted} age gradients ($\sim -0.11\; {\rm dex}/R_{\rm e}$), suggesting an "inside-out" formation of discs. We detect negative metallicity gradients in both early and late-type galaxies, but these are significantly steeper in late-types, suggesting that radial dependence of chemical enrichment processes and the effect of gas inflow and metal transport are far more pronounced in discs. Metallicity gradients of both morphological classes correlate with galaxy mass, with negative metallicity gradients becoming steeper with increasing galaxy mass. The correlation with mass is stronger for late-type galaxies, with a slope of $d(\nabla [Z/H])/d(\log M)\sim -0.2\pm 0.05\;$, compared to $d(\nabla [Z/H])/d(\log M)\sim -0.05\pm 0.05\;$ for early-types. This result suggests that the merger history plays a relatively small role in shaping metallicity gradients of galaxies.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

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.