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).

We measured gas-phase metallicity, ionisation parameter and dust extinction for 1795 representative local star-forming galaxies using integral field spectroscopy from the SDSS-IV MaNGA survey. We self-consistently derive these quantities by comparing observed line fluxes with photoionisation models using a Bayesian framework. We also present the first comprehensive study of the [SIII]$\lambda\lambda$9069,9532 nebular lines, which have long been predicted to be ideal tracers of the ionisation parameter. Unfortunately, we find that current photoionisation models substantially over-predict [SIII] lines intensity, while broadly reproducing other observed optical line ratios. We discuss how to nonetheless make use of the information provided by [SIII] lines by setting a prior on the ionisation parameter. Following this approach, we derive spatially-resolved maps and radial profiles of metallicity and ionisation parameter. The metallicity radial profiles are comparable with previous works, with metallicity declining toward the outer parts and a flattening in the central regions, in agreement with infall models of galaxy formation, that predict that spiral discs build up through accretion of material, which leads to an inside-out growth. On the other hand, ionisation parameter radial profiles are flat for low-mass galaxies, while their slope becomes positive as galaxy mass increases. However, the ionisation parameter maps we obtain are clumpy, especially for low-mass galaxies. Ionisation parameter is tightly correlated with the H$\alpha$ equivalent width [EW(H$\alpha$)], following a nearly universal relation, which we attribute to the change of the spectral shape of ionising sources due to ageing of HII regions. We derive a positive correlation between ionisation parameter and metallicity at fixed EW(H$\alpha$), in disagreement with previous theoretical works expecting an anti-correlation.

In this paper we use strong line calibrations of N/O and O/H in MaNGA spaxel data to explore the systematics introduced by variations in N/O on various strong-line metallicity diagnostics. We find radial variations in N/O at fixed O/H which correlate with total galaxy stellar-mass; and which can induce $\sim 40 \%$ systematic uncertainties in oxygen abundance gradients when nitrogen-dependent abundance calibrations are used. Empirically, we find that these differences are associated with variation in the local star formation efficiency, as predicted by recent chemical evolution models for galaxies, but we cannot rule out other processes such as radial migration and the accretion of passive dwarf galaxies also playing a role.

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 investigate the environmental dependence of the local gas-phase metallicity in a sample of star-forming galaxies from the MaNGA survey. Satellite galaxies with stellar masses in the range $9<\log(M_{*}/M_{\odot})<10$ are found to be $\sim 0.05 \, \mathrm{dex}$ higher in metallicity than centrals of similar stellar mass. Within the low-mass satellite population, we find that the interstellar medium (ISM) metallicity depends most strongly on the stellar mass of the galaxy that is central to the halo, though there is no obvious difference in the metallicity gradients. At fixed total stellar mass, the satellites of high mass ($M_{*}>10^{10.5} \, \mathrm{M_{\odot}}$) centrals are $\sim 0.1 \, \mathrm{dex}$ more metal rich than satellites of low-mass ($M_{*} < 10^{10} \, \mathrm{M_{\odot}}$) centrals, controlling for local stellar mass surface density and gas fraction. Fitting a gas-regulator model to the spaxel data, we are able to account for variations in the local gas fraction, stellar mass surface density and local escape velocity-dependent outflows. We find that the best explanation for the metallicity differences is the variation in the average metallicity of accreted gas between different environments that depends on the stellar mass of the dominant galaxies in each halo. This is interpreted as evidence for the exchange of enriched gas between galaxies in dense environments that is predicted by recent simulations.

SDSS-IV MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is the largest integral-field spectroscopy survey to date, aiming to observe a statistically representative sample of 10,000 low-redshift galaxies. In this paper we study the reliability of the emission-line fluxes and kinematic properties derived by the MaNGA Data Analysis Pipeline (DAP). We describe the algorithmic choices made in the DAP with regards to measuring emission-line properties, and the effect of our adopted strategy of simultaneously fitting the continuum and line emission. The effect of random errors are quantified by studying various fit-quality metrics, idealized recovery simulations and repeat observations. This analysis demonstrates that the emission lines are well-fit in the vast majority of the MaNGA dataset and the derived fluxes and errors are statistically robust. The systematic uncertainty on emission-line properties introduced by the choice of continuum templates is also discussed. In particular, we test the effect of using different stellar libraries and simple stellar-population models on the derived emission-line fluxes and the effect of introducing different tying prescriptions for the emission-line kinematics. We show that these effects can generate large ($>$ 0.2 dex) discrepancies at low signal-to-noise and for lines with low equivalent width (EW); however, the combined effect is noticeable even for H$\alpha$ EW $>$ 6~\AA. We provide suggestions for optimal use of the data provided by SDSS data release 15 and propose refinements on the \DAP\ for future MaNGA data releases.

We present a method of fitting optical spectra of galaxies using a basis set of six vectors obtained from principal component analysis (PCA) of a library of synthetic spectra of 40000 star formation histories (SFHs). Using this library, we provide estimates of resolved effective stellar mass-to-light ratio ($\log \Upsilon^*$) for thousands of galaxies from the SDSS-IV/MaNGA integral-field spectroscopic survey. Using a testing framework built on additional synthetic SFHs, we show that the estimates of stellar mass-to-light ratio are reliable (as are their uncertainties) at a variety of signal-to-noise ratios, stellar metallicities, and dust attenuation conditions. Finally, we describe the future release of the resolved stellar mass-to-light ratios as a SDSS-IV/MaNGA Value-Added Catalog (VAC) and provide a link to the software used to conduct this analysis.

A galaxy's stellar mass is one of its most fundamental properties, but it remains challenging to measure reliably. With the advent of very large optical spectroscopic surveys, efficient methods that can make use of low signal-to-noise spectra are needed. With this in mind, we created a new software package for estimating effective stellar mass-to-light ratios $\log \Upsilon^*$ that uses principal component analysis(PCA) basis set to optimize the comparison between observed spectra and a large library of stellar population synthesis models. In Paper I, we showed that a with a set of six PCA basis vectors we could faithfully represent most optical spectra from the Mapping Nearby Galaxies at APO (MaNGA) survey;and we tested the accuracy of our M/L estimates using synthetic spectra. Here, we explore sources of systematic error in our mass measurements by comparing our new measurements to data from the literature. We compare our stellar mass surface density estimates to kinematics-derived dynamical mass surface density measurements from the DiskMass Survey and find some tension between the two which could be resolved if the disk scale-heights used in the kinematic analysis were overestimated by a factor of $\sim 1.5$. We formulate an aperture-corrected stellar mass catalog for the MaNGA survey, and compare to previous stellar mass estimates based on multi-band optical photometry, finding typical discrepancies of 0.1 dex. Using the spatially resolved MaNGA data, we evaluate the impact of estimating total stellar masses from spatially unresolved spectra, and we explore how the biases that result from unresolved spectra depend upon the galaxy's dust extinction and star formation rate. Finally, we describe a SDSS Value-Added Catalog which will include both spatially resolved and total (aperture-corrected) stellar masses for MaNGA galaxies.

Bars in galaxies are thought to stimulate both inflow of material and radial mixing along them. Observational evidence for this mixing has been inconclusive so far however, limiting the evaluation of the impact of bars on galaxy evolution. We now use results from the MaNGA integral field spectroscopic survey to characterise radial stellar age and metallicity gradients along the bar and outside the bar in 128 strongly barred galaxies. We find that age and metallicity gradients are flatter in the barred regions of almost all barred galaxies when compared to corresponding disk regions at the same radii. Our results re-emphasize the key fact that by azimuthally averaging integral field spectroscopic data one loses important information from non-axisymmetric galaxy components such as bars and spiral arms. We interpret our results as observational evidence that bars are radially mixing material in galaxies of all stellar masses, and for all bar morphologies and evolutionary stages.

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.

We explore the radial distribution of star formation in galaxies in the SAMI Galaxy Survey as a function of their local group environment. Using a sample of galaxies in groups (with halo masses less than $ \simeq 10^{14} \, \mathrm{M_{\odot}}$) from the Galaxy And Mass Assembly Survey, we find signatures of environmental quenching in high-mass groups ($M_{G} > 10^{12.5} \, \mathrm{M_{\odot}}$). The mean integrated specific star formation rate of star-forming galaxies in high-mass groups is lower than for galaxies in low-mass groups or that are ungrouped, with $\Delta \log(sSFR/\mathrm{yr^{-1}}) = 0.45 \pm 0.07$. This difference is seen at all galaxy stellar masses. In high-mass groups, star-forming galaxies more massive than $M_{*} \sim 10^{10} \, \mathrm{M_{\odot}}$ have centrally-concentrated star formation. These galaxies also lie below the star-formation main sequence, suggesting they may be undergoing outside-in quenching. Lower mass galaxies in high-mass groups do not show evidence of concentrated star formation. In groups less massive than $M_{G} = 10^{12.5} \, \mathrm{M_{\odot}}$ we do not observe these trends. In this regime we find a modest correlation between centrally-concentrated star formation and an enhancement in total star formation rate, consistent with triggered star formation in these galaxies.

We present the second major release of data from the SAMI Galaxy Survey. Data Release Two includes data for 1559 galaxies, about 50% of the full survey. Galaxies included have a redshift range 0.004 < z < 0.113 and a large stellar mass range 7.5 < log (M_star/M_sun) < 11.6. The core data for each galaxy consist of two primary spectral cubes covering the blue and red optical wavelength ranges. For each primary cube we also provide three spatially binned spectral cubes and a set of standardised aperture spectra. For each core data product we provide a set of value-added data products. This includes all emission line value-added products from Data Release One, expanded to the larger sample. In addition we include stellar kinematic and stellar population value-added products derived from absorption line measurements. The data are provided online through Australian Astronomical Optics' Data Central. We illustrate the potential of this release by presenting the distribution of ~350,000 stellar velocity dispersion measurements from individual spaxels as a function of R/R_e, divided in four galaxy mass bins. In the highest stellar mass bin (log (M_star/M_sun)>11), the velocity dispersion strongly increases towards the centre, whereas below log (M_star/M_sun)<10 we find no evidence for a clear increase in the central velocity dispersion. This suggests a transition mass around log (M_star/M_sun) ~10 for galaxies with or without a dispersion-dominated bulge.

In order to determine the causes of kinematic asymmetry in the H$\alpha$ gas in the SAMI Galaxy Survey sample, we investigate the comparative influences of environment and intrinsic properties of galaxies on perturbation. We use spatially resolved H$\alpha$ velocity fields from the SAMI Galaxy Survey to quantify kinematic asymmetry ($\overline{v_{asym}}$) in nearby galaxies and environmental and stellar mass data from the GAMA survey. {We find that local environment, measured as distance to nearest neighbour, is inversely correlated with kinematic asymmetry for galaxies with $\mathrm{\log(M_*/M_\odot)}>10.0$, but there is no significant correlation for galaxies with $\mathrm{\log(M_*/M_\odot)}<10.0$. Moreover, low mass galaxies ($\mathrm{\log(M_*/M_\odot)}<9.0$) have greater kinematic asymmetry at all separations, suggesting a different physical source of asymmetry is important in low mass galaxies.} We propose that secular effects derived from gas fraction and gas mass may be the primary causes of asymmetry in low mass galaxies. High gas fraction is linked to high $\frac{\sigma_{m}}{V}$ (where $\sigma_m$ is H$\alpha$ velocity dispersion and $V$ the rotation velocity), which is strongly correlated with $\overline{v_{asym}}$, and galaxies with $\log(M_*/M_\odot)<9.0$ have offset $\overline{\frac{\sigma_{m}}{V}}$ from the rest of the sample. Further, asymmetry as a fraction of dispersion decreases for galaxies with $\log(M_*/M_\odot)<9.0$. Gas mass and asymmetry are also inversely correlated in our sample. We propose that low gas masses in dwarf galaxies may lead to asymmetric distribution of gas clouds, leading to increased relative turbulence.

We present the ~800 star formation rate maps for the SAMI Galaxy Survey based on H{\alpha} emission maps, corrected for dust attenuation via the Balmer decrement, that are included in the SAMI Public Data Release 1. We mask out spaxels contaminated by non-stellar emission using the [O III]/H{\beta}, [N II]/H{\alpha}, [S II]/H{\alpha}, and [O I]/H{\alpha} line ratios. Using these maps, we examine the global and resolved star-forming main sequences of SAMI galaxies as a function of morphology, environmental density, and stellar mass. Galaxies further below the star-forming main sequence are more likely to have flatter star formation profiles. Early-type galaxies split into two populations with similar stellar masses and central stellar mass surface densities. The main sequence population has centrally-concentrated star formation similar to late-type galaxies, while galaxies >3{\sigma} below the main sequence show significantly reduced star formation most strikingly in the nuclear regions. The split populations support a two-step quenching mechanism, wherein halo mass first cuts off the gas supply and remaining gas continues to form stars until the local stellar mass surface density can stabilize the reduced remaining fuel against further star formation. Across all morphologies, galaxies in denser environments show a decreased specific star formation rate from the outside in, supporting an environmental cause for quenching, such as ram-pressure stripping or galaxy interactions.

We investigate the Tully-Fisher Relation (TFR) for a morphologically and kine- matically diverse sample of galaxies from the SAMI Galaxy Survey using 2 dimensional spatially resolved Halpha velocity maps and find a well defined relation across the stellar mass range of 8.0 < log(M*) < 11.5. We use an adaptation of kinemetry to parametrise the kinematic Halpha asymmetry of all galaxies in the sample, and find a correlation between scatter (i.e. residuals off the TFR) and asymmetry. This effect is pronounced at low stellar mass, corresponding to the inverse relationship between stellar mass and kinematic asymmetry found in previous work. For galaxies with log(M*) < 9.5, 25 +/- 3% are scattered below the root mean square (RMS) of the TFR, whereas for galaxies with log(M*) > 9.5 the fraction is 10 +/- 1% We use 'simulated slits' to directly compare our results with those from long slit spectroscopy and find that aligning slits with the photometric, rather than the kinematic, position angle, increases global scatter below the TFR. Further, kinematic asymmetry is correlated with misalignment between the photometric and kinematic position angles. This work demonstrates the value of 2D spatially resolved kinematics for accurate TFR studies; integral field spectroscopy reduces the underestimation of rotation velocity that can occur from slit positioning off the kinematic axis.

We present the first major release of data from the SAMI Galaxy Survey. This data release focuses on the emission-line physics of galaxies. Data Release One includes data for 772 galaxies, about 20% of the full survey. Galaxies included have the redshift range 0.004 < z < 0.092, a large mass range (7.6 < log(Mstellar/M$_\odot$) < 11.6), and star-formation rates of 10^-4 to 10^1\ M$_\odot$/yr. For each galaxy, we include two spectral cubes and a set of spatially resolved 2D maps: single- and multi-component emission-line fits (with dust extinction corrections for strong lines), local dust extinction and star-formation rate. Calibration of the fibre throughputs, fluxes and differential-atmospheric-refraction has been improved over the Early Data Release. The data have average spatial resolution of 2.16 arcsec (FWHM) over the 15~arcsec diameter field of view and spectral (kinematic) resolution R=4263 (sigma=30km/s) around Halpha. The relative flux calibration is better than 5\% and absolute flux calibration better than $\pm0.22$~mag, with the latter estimate limited by galaxy photometry. The data are presented online through the Australian Astronomical Observatory's Data Central.

The fraction of galaxies supported by internal rotation compared to galaxies stabilized by internal pressure provides a strong constraint on galaxy formation models. In integral field spectroscopy surveys, this fraction is biased because survey instruments typically only trace the inner parts of the most massive galaxies. We present aperture corrections for the two most widely used stellar kinematic quantities $V/\sigma$ and $\lambda_{R}$. Our demonstration involves integral field data from the SAMI Galaxy Survey and the ATLAS$^{\rm{3D}}$ Survey. We find a tight relation for both $V/\sigma$ and $\lambda_{R}$ when measured in different apertures that can be used as a linear transformation as a function of radius, i.e., a first-order aperture correction. We find that $V/\sigma$ and $\lambda_{R}$ radial growth curves are well approximated by second order polynomials. By only fitting the inner profile (0.5$R_{\rm{e}}$), we successfully recover the profile out to one $R_{\rm{e}}$ if a constraint between the linear and quadratic parameter in the fit is applied. However, the aperture corrections for $V/\sigma$ and $\lambda_{R}$ derived by extrapolating the profiles perform as well as applying a first-order correction. With our aperture-corrected $\lambda_{R}$ measurements, we find that the fraction of slow rotating galaxies increases with stellar mass. For galaxies with $\log M_{*}/M_{\odot}>$ 11, the fraction of slow rotators is $35.9\pm4.3$ percent, but is underestimated if galaxies without coverage beyond one $R_{\rm{e}}$ are not included in the sample ($24.2\pm5.3$ percent). With measurements out to the largest aperture radius the slow rotator fraction is similar as compared to using aperture corrected values ($38.3\pm4.4$ percent). Thus, aperture effects can significantly bias stellar kinematic IFS studies, but this bias can now be removed with the method outlined here.

Recent cosmological hydrodynamical simulations suggest that integral field spectroscopy can connect the high-order stellar kinematic moments h3 (~skewness) and h4 (~kurtosis) in galaxies to their cosmological assembly history. Here, we assess these results by measuring the stellar kinematics on a sample of 315 galaxies, without a morphological selection, using 2D integral field data from the SAMI Galaxy Survey. A proxy for the spin parameter ($\lambda_{R_e}$) and ellipticity ($\epsilon_e$) are used to separate fast and slow rotators; there exists a good correspondence to regular and non-regular rotators, respectively, as also seen in earlier studies. We confirm that regular rotators show a strong h3 versus $V/\sigma$ anti-correlation, whereas quasi-regular and non-regular rotators show a more vertical relation in h3 and $V/\sigma$. Motivated by recent cosmological simulations, we develop an alternative approach to kinematically classify galaxies from their individual h3 versus $V/\sigma$ signatures. We identify five classes of high-order stellar kinematic signatures using Gaussian mixture models. Class 1 corresponds to slow rotators, whereas Classes 2-5 correspond to fast rotators. We find that galaxies with similar $\lambda_{R_e}-\epsilon_e$ values can show distinctly different h3-$V/\sigma$ signatures. Class 5 objects are previously unidentified fast rotators that show a weak h3 versus $V/\sigma$ anti-correlation. These objects are predicted to be disk-less galaxies formed by gas-poor mergers. From morphological examination, however, there is evidence for large stellar disks. Instead, Class 5 objects are more likely disturbed galaxies, have counter-rotating bulges, or bars in edge-on galaxies. Finally, we interpret the strong anti-correlation in h3 versus $V/\sigma$ as evidence for disks in most fast rotators, suggesting a dearth of gas-poor mergers among fast rotators.

We study the properties of kinematically disturbed galaxies in the SAMI Galaxy Survey using a quantitative criterion, based on kinemetry (Krajnovic et al.). The approach, similar to the application of kinemetry by Shapiro et al. uses ionised gas kinematics, probed by H{\alpha} emission. By this method 23+/-7% of our 360-galaxy sub-sample of the SAMI Galaxy Survey are kinematically asymmetric. Visual classifications agree with our kinemetric results for 90% of asymmetric and 95% of normal galaxies. We find stellar mass and kinematic asymmetry are inversely correlated and that kinematic asymmetry is both more frequent and stronger in low-mass galaxies. This builds on previous studies that found high fractions of kinematic asymmetry in low mass galaxies using a variety of different methods. Concentration of star forma- tion and kinematic disturbance are found to be correlated, confirming results found in previous work. This effect is stronger for high mass galaxies (log(M*) > 10) and indicates that kinematic disturbance is linked to centrally concentrated star formation. Comparison of the inner (within 0.5Re) and outer H{\alpha} equivalent widths of asymmetric and normal galaxies shows a small but significant increase in inner equivalent width for asymmetric galaxies.

We use data from the Sydney-AAO Multi-Object Integral Field Spectrograph (SAMI) Galaxy Survey and the Galaxy And Mass Assembly (GAMA) survey to investigate the spatially-resolved signatures of the environmental quenching of star formation in galaxies. Using dust-corrected measurements of the distribution of H$\alpha$ emission we measure the radial profiles of star formation in a sample of 201 star-forming galaxies covering three orders of magnitude in stellar mass (M$_{*}$; $10^{8.1}$-$10^{10.95}\, $M$_{\odot}$) and in $5^{th}$ nearest neighbour local environment density ($\Sigma_{5}$; $10^{-1.3}$-$10^{2.1}\,$Mpc$^{-2}$). We show that star formation rate gradients in galaxies are steeper in dense ($\log_{10}(\Sigma_{5}/$Mpc$^{2})>0.5$) environments by $0.58\pm 0.29\, dex\, $r$_{e}^{-1}$ in galaxies with stellar masses in the range $10^{10}<$M$_{*}/$M$_{\odot}<10^{11}$ and that this steepening is accompanied by a reduction in the integrated star formation rate. However, for any given stellar mass or environment density the star-formation morphology of galaxies shows large scatter. We also measure the degree to which the star formation is centrally concentrated using the unitless scale-radius ratio ($r_{50,H\alpha}/r_{50,cont}$), which compares the extent of ongoing star formation to previous star formation. With this metric we find that the fraction of galaxies with centrally concentrated star formation increases with environment density, from $\sim 5\pm 4\%$ in low-density environments ($\log_{10}(\Sigma_{5}/$Mpc$^{2})<0.0$) to $30\pm 15\%$ in the highest density environments ($\log_{10}(\Sigma_{5}/$Mpc$^{2})>1.0$). These lines of evidence strongly suggest that with increasing local environment density the star formation in galaxies is suppressed, and that this starts in their outskirts such that quenching occurs in an outside-in fashion in dense environments and is not instantaneous.

Hector will be the new massively-multiplexed integral field spectroscopy (IFS) instrument for the Anglo-Australian Telescope (AAT) in Australia and the next main dark-time instrument for the observatory. Based on the success of the SAMI instrument, which is undertaking a 3400-galaxy survey, the integral field unit (IFU) imaging fibre bundle (hexabundle) technology under-pinning SAMI is being improved to a new innovative design for Hector. The distribution of hexabundle angular sizes is matched to the galaxy survey properties in order to image 90% of galaxies out to 2 effective radii. 50-100 of these IFU imaging bundles will be positioned by 'starbug' robots across a new 3-degree field corrector top end to be purpose-built for the AAT. Many thousand fibres will then be fed into new replicable spectrographs. Fundamentally new science will be achieved compared to existing instruments due to Hector's wider field of view (3 degrees), high positioning efficiency using starbugs, higher spectroscopic resolution (R~3000-5500 from 3727-7761A, with a possible redder extension later) and large IFUs (up to 30 arcsec diameter with 61-217 fibre cores). A 100,000 galaxy IFS survey with Hector will decrypt how the accretion and merger history and large-scale environment made every galaxy different in its morphology and star formation history. The high resolution, particularly in the blue, will make Hector the only instrument to be able to measure higher-order kinematics for galaxies down to much lower velocity dispersion than in current large IFS galaxy surveys, opening up a wealth of new nearby galaxy science.

We present LZIFU (LaZy-IFU), an IDL toolkit for fitting multiple emission lines simultaneously in integral field spectroscopy (IFS) data. LZIFU is useful for the investigation of the dynamical, physical and chemical properties of gas in galaxies. LZIFU has already been applied to many world-class IFS instruments and large IFS surveys, including the Wide Field Spectrograph, the new Multi Unit Spectroscopic Explorer (MUSE), the Calar Alto Legacy Integral Field Area (CALIFA) survey, the Sydney-Australian-astronomical-observatory Multi-object Integral-field spectrograph (SAMI) Galaxy Survey. Here we describe in detail the structure of the toolkit, and how the line fluxes and flux uncertainties are determined, including the possibility of having multiple distinct kinematic components. We quantify the performance of LZIFU, demonstrating its accuracy and robustness. We also show examples of applying LZIFU to CALIFA and SAMI data to construct emission line and kinematic maps, and investigate complex, skewed line profiles presented in IFS data. The code is made available to the astronomy community through github. LZIFU will be further developed over time to other IFS instruments, and to provide even more accurate line and uncertainty estimates.

In the low redshift Universe (z<0.3), our view of galaxy evolution is primarily based on fibre optic spectroscopy surveys. Elaborate methods have been developed to address aperture effects when fixed aperture sizes only probe the inner regions for galaxies of ever decreasing redshift or increasing physical size. These aperture corrections rely on assumptions about the physical properties of galaxies. The adequacy of these aperture corrections can be tested with integral-field spectroscopic data. We use integral-field spectra drawn from 1212 galaxies observed as part of the SAMI Galaxy Survey to investigate the validity of two aperture correction methods that attempt to estimate a galaxy's total instantaneous star formation rate. We show that biases arise when assuming that instantaneous star formation is traced by broadband imaging, and when the aperture correction is built only from spectra of the nuclear region of galaxies. These biases may be significant depending on the selection criteria of a survey sample. Understanding the sensitivities of these aperture corrections is essential for correct handling of systematic errors in galaxy evolution studies.

We have observed two kinematically offset active galactic nuclei (AGN), whose ionised gas is at a different line-of-sight velocity to their host galaxies, with the SAMI integral field spectrograph (IFS). One of the galaxies shows gas kinematics very different to the stellar kinematics, indicating a recent merger or accretion event. We demonstrate that the star formation associated with this event was triggered within the last 100 Myr. The other galaxy shows simple disc rotation in both gas and stellar kinematics, aligned with each other, but in the central region has signatures of an outflow driven by the AGN. Other than the outflow, neither galaxy shows any discontinuity in the ionised gas kinematics at the galaxy's centre. We conclude that in these two cases there is no direct evidence of the AGN being in a supermassive black hole binary system. Our study demonstrates that selecting kinematically offset AGN from single-fibre spectroscopy provides, by definition, samples of kinematically peculiar objects, but IFS or other data are required to determine their true nature.

The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey includes four volume limited galaxy samples based on cuts in a proxy for stellar mass, along with low-stellar mass dwarf galaxies all selected from the Galaxy And Mass Assembly (GAMA) survey. The GAMA regions were selected because of the vast array of ancillary data available, including ultraviolet through to radio bands. These fields are on the celestial equator at 9, 12, and 14.5 hours, and cover a total of 144 square degrees (in GAMA-I). Higher density environments are also included with the addition of eight clusters. The clusters have spectroscopy from 2dFGRS and SDSS and photometry in regions covered by the Sloan Digital Sky Survey (SDSS) and/or VLT Survey Telescope/ATLAS. The aim is to cover a broad range in stellar mass and environment, and therefore the primary survey targets cover redshifts 0.004 < z < 0.095, magnitudes r$_{pet}$ < 19.4, stellar masses $10^{7} - 10^{12}$ M$_{sol}$, and environments from isolated field galaxies through groups to clusters of $10^{15}$ M$_{sol}$.

We present a methodology for the regularisation and combination of sparse sampled and irregularly gridded observations from fibre-optic multi-object integral-field spectroscopy. The approach minimises interpolation and retains image resolution on combining sub-pixel dithered data. We discuss the methodology in the context of the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey underway at the Anglo-Australian Telescope. The SAMI instrument uses 13 fibre bundles to perform high-multiplex integral-field spectroscopy across a one degree diameter field of view. The SAMI Galaxy Survey is targeting 3000 galaxies drawn from the full range of galaxy environments. We demonstrate the subcritical sampling of the seeing and incomplete fill factor for the integral-field bundles results in only a 10% degradation in the final image resolution recovered. We also implement a new methodology for tracking covariance between elements of the resulting datacubes which retains 90% of the covariance information while incurring only a modest increase in the survey data volume.

We present the Early Data Release of the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The SAMI Galaxy Survey is an ongoing integral field spectroscopic survey of ~3400 low-redshift (z<0.12) galaxies, covering galaxies in the field and in groups within the Galaxy And Mass Assembly (GAMA) survey regions, and a sample of galaxies in clusters. In the Early Data Release, we publicly release the fully calibrated datacubes for a representative selection of 107 galaxies drawn from the GAMA regions, along with information about these galaxies from the GAMA catalogues. All datacubes for the Early Data Release galaxies can be downloaded individually or as a set from the SAMI Galaxy Survey website. In this paper we also assess the quality of the pipeline used to reduce the SAMI data, giving metrics that quantify its performance at all stages in processing the raw data into calibrated datacubes. The pipeline gives excellent results throughout, with typical sky subtraction residuals in the continuum of 0.9-1.2 per cent, a relative flux calibration uncertainty of 4.1 per cent (systematic) plus 4.3 per cent (statistical), and atmospheric dispersion removed with an accuracy of 0."09, less than a fifth of a spaxel.

We present the discovery of a luminous unresolved H II complex on the edge of dwarf galaxy GAMA J141103.98-003242.3 using data from the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. This dwarf galaxy is situated at a distance of ~100 Mpc and contains an unresolved region of H II emission that contributes ~70 per cent of the galaxy's H_alpha luminosity, located at the top end of established H II region luminosity functions. For the H II complex, we measure a star-formation rate of 0.147\pm0.041 M_solar yr^-1 and a metallicity of 12+log(O/H) = 8.01\pm0.05 that is lower than the rest of the galaxy by ~0.2 dex. Data from the H I Parkes All-Sky Survey (HIPASS) indicate the likely presence of neutral hydrogen in the galaxy to potentially fuel ongoing and future star-forming events. We discuss various triggering mechanisms for the intense star-formation activity of this H II complex, where the kinematics of the ionised gas are well described by a rotating disc and do not show any features indicative of interactions. We show that SAMI is an ideal instrument to identify similar systems to GAMA J141103.98-003242.3, and the SAMI Galaxy Survey is likely to find many more of these systems to aid in the understanding of their formation and evolution.

Optical positions from the SuperCOSMOS Sky Survey have been compared in detail with accurate radio positions that define the second realisation of the International Celestial Reference Frame (ICRF2). The comparison was limited to the IIIaJ plates from the UK/AAO and Oschin (Palomar) Schmidt telescopes. A total of 1373 ICRF2 sources was used, with the sample restricted to stellar objects brighter than $B_J=20$ and Galactic latitudes $|b|>10^{\circ}$. Position differences showed an rms scatter of $0.16''$ in right ascension and declination. While overall systematic offsets were $<0.1''$ in each hemisphere, both the systematics and scatter were greater in the north.