Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and extragalactic surveys will be carried out simultaneously, so the combined target density will strongly vary. In this paper, we describe a new tiling algorithm that can naturally deal with the large target density variations on the sky and which automatically handles the different exposure times of targets. The tiling pattern is modelled as a marked point process, which is characterised by a probability density that integrates the requirements imposed by the 4MOST survey. The optimal tilling pattern with respect to the defined model is estimated by the tiles configuration that maximises the proposed probability density. In order to achieve this maximisation a simulated annealing algorithm is implemented. The algorithm automatically finds an optimal tiling pattern and assigns a tentative sky brightness condition and exposure time for each tile, while minimising the total execution time that is needed to observe the list of targets in the combined input catalogue of all surveys. Hence, the algorithm maximises the long-term observing efficiency and provides an optimal tiling solution for the survey. While designed for the 4MOST survey, the algorithm is flexible and can with simple modifications be applied to any other multi-object spectroscopic survey.

We present improved modelling of the redshift-space distortions of galaxy clustering that arise from peculiar velocities. We create mock galaxy catalogues in the framework of the halo model, using data from the Bolshoi project. These mock galaxy populations are inserted into the haloes with additional degrees of freedom that govern spatial and kinematical biases of the galaxy populations relative to the dark matter. We explore this generalised halo model with an MCMC algorithm, comparing the predictions to data from the Galaxy And Mass Assembly (GAMA) survey, and thus derive one of the first constraints on the detailed kinematic degrees of freedom for satellite galaxies within haloes. With this approach, the distortions of the redshift-space galaxy autocorrelations can be accounted for down to spatial separations close to 10 kpc, opening the prospect of improved RSD measurements of the perturbation growth rate by the inclusion of data from nonlinear scales.

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

We use multi-wavelength data from the Galaxy and Mass Assembly (GAMA) survey to explore the cause of red optical colours in nearby (0.002<z<0.06) spiral galaxies. We show that the colours of red spiral galaxies are a direct consequence of some environment-related mechanism(s) which has removed dust and gas, leading to a lower star formation rate. We conclude that this process acts on long timescales (several Gyr) due to a lack of morphological transformation associated with the transition in optical colour. The sSFR and dust-to-stellar mass ratio of red spiral galaxies is found to be statistically lower than blue spiral galaxies. On the other hand, red spirals are on average $0.9$ dex more massive, and reside in environments 2.6 times denser than their blue counterparts. We find no evidence of excessive nuclear activity, or higher inclination angles to support these as the major causes for the red optical colours seen in >= 47% of all spirals in our sample. Furthermore, for a small subsample of our spiral galaxies which are detected in HI, we find that the SFR of gas-rich red spiral galaxies is lower by ~1 dex than their blue counterparts.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

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 Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission.

WAVES is designed to study the growth of structure, mass and energy on scales of ~1 kpc to ~10 Mpc over a 7 Gyr timeline. On the largest length scales (1-10 Mpc) WAVES will measure the structures defined by groups, filaments and voids, and their emergence over recent times. Comparisons with bespoke numerical simulations will be used to confirm, refine or refute the Cold Dark Matter paradigm. At intermediate length scales (10 kpc-1 Mpc) WAVES will probe the size and mass distribution of galaxy groups, as well as the galaxy merger rates, in order to directly measure the assembly of dark matter halos and stellar mass. On the smallest length scales (1-10 kpc) WAVES will provide accurate distance and environmental measurements to complement high-resolution space-based imaging to study the mass and size evolution of galaxy bulges, discs and bars. In total, WAVES will provide a panchromatic legacy dataset of ~1.6 million galaxies, firmly linking the very low ($z < 0.1$) and intermediate ($z \sim 0.8$) redshift Universe.

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 a new measurement of $E_{\rm G}$, which combines measurements of weak gravitational lensing, real-space galaxy clustering and redshift space distortions. This statistic was proposed as a consistency test of General Relativity (GR) that is insensitive to linear, deterministic galaxy bias and the matter clustering amplitude. We combine deep imaging data from KiDS with overlapping spectroscopy from 2dFLenS, BOSS DR12 and GAMA and find $E_{\rm G}(\overline{z}=0.267)=0.43 \pm 0.13$ (GAMA), $E_{\rm G}(\overline{z}=0.305)=0.27 \pm 0.08$ (LOWZ+2dFLOZ) and $E_{\rm G}(\overline{z}=0.554)=0.26 \pm 0.07$ (CMASS+2dFHIZ). We demonstrate that the existing tension in the value of the matter density parameter hinders the robustness of this statistic as solely a test of GR. We find that our $E_{\rm G}$ measurements, as well as existing ones in the literature, favour a lower matter density cosmology than the Cosmic Microwave Background. For a flat $\Lambda$CDM Universe and assuming GR, we find $\Omega_{\rm m}(z=0)=0.25\pm0.03$. With this paper we publicly release the 2dFLenS dataset at: \url{http://2dflens.swin.edu.au}.

The galaxies found in optical surveys fall in two distinct regions of a diagram of optical colour versus absolute magnitude: the red sequence and the blue cloud with the green valley in between. We show that the galaxies found in a submillimetre survey have almost the opposite distribution in this diagram, forming a `green mountain'. We show that these distinctive distributions follow naturally from a single, continuous, curved Galaxy Sequence in a diagram of specific star-formation rate versus stellar mass without there being the need for a separate star-forming galaxy Main Sequence and region of passive galaxies. The cause of the red sequence and the blue cloud is the geometric mapping between stellar mass/specific star-formation rate and absolute magnitude/colour, which distorts a continuous Galaxy Sequence in the diagram of intrinsic properties into a bimodal distribution in the diagram of observed properties. The cause of the green mountain is Malmquist bias in the submillimetre waveband, with submillimetre surveys tending to select galaxies on the curve of the Galaxy Sequence, which have the highest ratios of submillimetre-to-optical luminosity. This effect, working in reverse, causes galaxies on the curve of the Galaxy Sequence to be underrepresented in optical samples, deepening the green valley. The green valley is therefore not evidence (1) for there being two distinct populations of galaxies, (2) for galaxies in this region evolving more quickly than galaxies in the blue cloud and the red sequence, (c) for rapid quenching processes in the galaxy population.

Statistical studies of galaxy-galaxy interactions often utilise net change in physical properties of progenitors as a function of the separation between their nuclei to trace both the strength and the observable timescale of their interaction. In this study, we use two-point auto, cross and mark correlation functions to investigate the extent to which small-scale clustering properties of star forming galaxies can be used to gain physical insight into galaxy-galaxy interactions between galaxies of similar optical brightness and stellar mass. The Halpha star formers, drawn from the highly spatially complete Galaxy And Mass Assembly (GAMA) survey, show an increase in clustering on small separations. Moreover, the clustering strength shows a strong dependence on optical brightness and stellar mass, where (1) the clustering amplitude of optically brighter galaxies at a given separation is larger than that of optically fainter systems, (2) the small scale clustering properties (e.g. the strength, the scale at which the signal relative to the fiducial power law plateaus) of star forming galaxies appear to differ as a function of increasing optical brightness of galaxies. According to cross and mark correlation analyses, the former result is largely driven by the increased dust content in optically bright star forming galaxies. The latter could be interpreted as evidence of a correlation between interaction-scale and optical brightness of galaxies, where physical evidence of interactions between optically bright star formers, likely hosted within relatively massive halos, persist over larger separations than those between optically faint star formers.

We present the dust mass function (DMF) of 15,750 galaxies with redshift $z< 0.1$, drawn from the overlapping area of the GAMA and {\it H-}ATLAS surveys. The DMF is derived using the density corrected $V_{\rm max}$ method, where we estimate $V_{\rm max}$ using: (i) the normal photometric selection limit ($pV_{\rm max}$) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65\pm0.18)\times 10^{7}\,h^2_{70}\, M_{\odot}$, $\alpha=(1.22\pm 0.01)$, $\phi^{*}=(6.26\pm 0.28)\times 10^{-3}\,h^3_{70}\,\rm Mpc^{-3}\,dex^{-1}$. The resulting dust mass density parameter integrated down to $10^4\,M_{\odot}$ is $\Omega_{\rm d}=(1.11 \pm0.02)\times 10^{-6}$ which implies the mass fraction of baryons in dust is $f_{m_b}=(2.40\pm0.04)\times 10^{-5}$; cosmic variance adds an extra 7-17\,per\,cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long timescales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of $\Omega_{\rm d}=(0.88\pm0.03)\times 10^{-6}$ and $\Omega_{\rm d}=(0.060\pm0.005)\times 10^{-6}$ for late-types and early-types respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late-types is well matched by the GMSF scaled by $(8.07\pm0.35) \times 10^{-4}$.

We explore how the group environment may affect the evolution of star-forming galaxies. We select 1197 Galaxy And Mass Assembly (GAMA) groups at $0.05\leq z \leq 0.2$ and analyze the projected phase space (PPS) diagram, i.e. the galaxy velocity as a function of projected group-centric radius, as a local environmental metric in the low-mass halo regime $10^{12}\leq (M_{200}/M_{\odot})< 10^{14}$. We study the properties of star-forming group galaxies, exploring the correlation of star formation rate (SFR) with radial distance and stellar mass. We find that the fraction of star-forming group members is higher in the PPS regions dominated by recently accreted galaxies, whereas passive galaxies dominate the virialized regions. We observe a small decline in specific SFR of star-forming galaxies towards the group center by a factor $\sim 1.2$ with respect to field galaxies. Similar to cluster studies, we conclude for low-mass halos that star-forming group galaxies represent an infalling population from the field to the halo and show suppressed star formation.

Aims. We aim to investigate if the environment (characterised by the host dark matter halo mass) plays any role in shaping the galaxy star formation main sequence (MS). Methods. The Galaxy and Mass Assembly project (GAMA) combines a spectroscopic survey with photometric information in 21 bands from the far-ultraviolet (FUV) to the far-infrared (FIR). Stellar masses and dust-corrected star-formation rates (SFR) are derived from spectral energy distribution (SED) modelling using MAGPHYS. We use the GAMA galaxy group catalogue to examine the variation of the fraction of star-forming galaxies (SFG) and properties of the MS with respect to the environment. Results. We examine the environmental dependence for stellar mass selected samples without preselecting star-forming galaxies and study any dependence on the host halo mass separately for centrals and satellites out to z ~ 0.3. We find the SFR distribution at fixed stellar mass can be described by the combination of two Gaussians (referred to as the star-forming Gaussian and the quiescent Gaussian). Using the observed bimodality to define SFG, we investigate how the fraction of SFG F(SFG) and properties of the MS change with environment. For centrals, the position of the MS is similar to the field but with a larger scatter. No significant dependence on halo mass is observed. For satellites, the position of the MS is almost always lower (by ~0.2 dex) compared to the field and the width is almost always larger. F(SFG) is similar between centrals (in different halo mass bins) and field galaxies. However, for satellites F(SFG) decreases with increasing halo mass and this dependence is stronger towards lower redshift.

In galaxy clusters, efficiently accreting active galactic nuclei (AGN) are preferentially located in the infall regions of the cluster projected phase-space, and are rarely found in the cluster core. This has been attributed to both an increase in triggering opportunities for infalling galaxies, and a reduction of those mechanisms in the hot, virialised, cluster core. Exploiting the depth and completeness ($98\,$per cent at $r<19.8\,$mag) of the Galaxy And Mass Assembly survey (GAMA), we probe down the group halo mass function to assess whether AGN are found in the same regions in groups as they are in clusters. We select 451 optical AGN from 7498 galaxies with $\log_{10}(M_*/\text{M}_\odot) > 9.9$ in 695 groups with $11.53\leq \log_{10}(M_{200}/\text{M}_\odot) \leq 14.56$ at $z<0.15$. By analysing the projected phase-space positions of these galaxies we demonstrate that when split both radially, and into physically derived infalling and core populations, AGN position within group projected phase-space is dependent on halo mass. For groups with $\log_{10}(M_{200}/\text{M}_\odot)>13.5$, AGN are preferentially found in the infalling galaxy population with $3.6\sigma$ confidence. At lower halo masses we observe no difference in AGN fraction between core and infalling galaxies. These observations support a model where a reduced number of low-speed interactions, ram pressure stripping and intra-group/cluster medium temperature, the dominance of which increase with halo mass, work to inhibit AGN in the cores of groups and clusters with $\log_{10}(M_{200}/\text{M}_\odot)>13.5$, but do not significantly affect nuclear activity in cores of less massive structures.

We describe data release 3 (DR3) of the Galaxy And Mass Assembly (GAMA) survey. The GAMA survey is a spectroscopic redshift and multi-wavelength photometric survey in three equatorial regions each of 60.0 deg^2 (G09, G12, G15), and two southern regions of 55.7 deg^2 (G02) and 50.6 deg^2 (G23). DR3 consists of: the first release of data covering the G02 region and of data on H-ATLAS sources in the equatorial regions; and updates to data on sources released in DR2. DR3 includes 154809 sources with secure redshifts across four regions. A subset of the G02 region is 95.5% redshift complete to r<19.8 over an area of 19.5 deg^2, with 20086 galaxy redshifts, that overlaps substantially with the XXL survey (X-ray) and VIPERS (redshift survey). In the equatorial regions, the main survey has even higher completeness (98.5%), and spectra for about 75% of H-ATLAS filler targets were also obtained. This filler sample extends spectroscopic redshifts, for probable optical counterparts to H-ATLAS sub-mm sources, to 0.8 mag deeper (r<20.6) than the GAMA main survey. There are 25814 galaxy redshifts for H-ATLAS sources from the GAMA main or filler surveys. GAMA DR3 is available at the survey website (www.gama-survey.org/dr3/).

We apply four statistical learning methods to a sample of $7941$ galaxies ($z<0.06$) from the Galaxy and Mass Assembly (GAMA) survey to test the feasibility of using automated algorithms to classify galaxies. Using $10$ features measured for each galaxy (sizes, colours, shape parameters \& stellar mass) we apply the techniques of Support Vector Machines (SVM), Classification Trees (CT), Classification Trees with Random Forest (CTRF) and Neural Networks (NN), returning True Prediction Ratios (TPRs) of $75.8\%$, $69.0\%$, $76.2\%$ and $76.0\%$ respectively. Those occasions whereby all four algorithms agree with each other yet disagree with the visual classification (`unanimous disagreement') serves as a potential indicator of human error in classification, occurring in $\sim9\%$ of ellipticals, $\sim9\%$ of Little Blue Spheroids, $\sim14\%$ of early-type spirals, $\sim21\%$ of intermediate-type spirals and $\sim4\%$ of late-type spirals \& irregulars. We observe that the choice of parameters rather than that of algorithms is more crucial in determining classification accuracy. Due to its simplicity in formulation and implementation, we recommend the CTRF algorithm for classifying future galaxy datasets. Adopting the CTRF algorithm, the TPRs of the 5 galaxy types are : E, $70.1\%$; LBS, $75.6\%$; S0-Sa, $63.6\%$; Sab-Scd, $56.4\%$ and Sd-Irr, $88.9\%$. Further, we train a binary classifier using this CTRF algorithm that divides galaxies into spheroid-dominated (E, LBS \& S0-Sa) and disk-dominated (Sab-Scd \& Sd-Irr), achieving an overall accuracy of $89.8\%$. This translates into an accuracy of $84.9\%$ for spheroid-dominated systems and $92.5\%$ for disk-dominated systems.

The galaxy pairwise velocity dispersion (PVD) can provide important tests of non-standard gravity and galaxy formation models. We describe measurements of the PVD of galaxies in the Galaxy and Mass Assembly (GAMA) survey as a function of projected separation and galaxy luminosity. Due to the faint magnitude limit ($r < 19.8$) and highly-complete spectroscopic sampling of the GAMA survey, we are able to reliably measure the PVD to smaller scales ($r_\bot = 0.01$ Mpc/h) than previous work. The measured PVD at projected separations $r_\bot <~ 1$ Mpc/h increases near-monotonically with increasing luminosity from $\sigma \approx 200$ km/s at $M_r = -17$ mag to $\sigma \approx 600$ km/s at $M_r \approx -22$ mag. Analysis of the Gonzalez-Perez (2014) GALFORM semi-analytic model yields no such trend of PVD with luminosity: the model over-predicts the PVD for faint galaxies. This is most likely a result of the model placing too many low-luminosity galaxies in massive halos.

In this work we measure and study the cross-correlation signal between a foreground sample of GAMA galaxies with spectroscopic redshifts in the range $0.2<z<0.8$, and a background sample of H-ATLAS galaxies with photometric redshifts $\gtrsim1.2$. It constitutes a substantial improvement over the cross-correlation measurements made by Gonzalez-Nuevo et al. (2014) with updated catalogues and wider area (with $S/N\gtrsim 5$ below 10' and reaching $S/N\sim 20$ below 30"). The better statistics allow us to split the sample in different redshift bins and to perform a tomographic analysis (with $S/N\gtrsim 3$ below 10 arcmin and reaching $S/N\sim 15$ below 30"). Moreover, we implement a halo model to extract astrophysical information about the background galaxies and the deflectors that are producing the lensing link between the foreground (lenses) and background (sources) samples. In the case of the sources, we find typical mass values in agreement with previous studies: a minimum halo mass to host a central galaxy, $M_{min}\sim 10^{12.26} M_\odot$, and a pivot halo mass to have at least one sub-halo satellite, $M_1\sim 10^{12.84} M_\odot$. However, the lenses are massive galaxies or even galaxy groups/clusters, with minimum mass of $M_{min}^{lens}\sim 10^{13.06} M_\odot$. Above a mass of $M_1^{lens}\sim 10^{14.57} M_\odot$ they contain at least one additional satellite galaxy which contributes to the lensing effect. The tomographic analysis shows that, while $M_1^{lens}$ is almost redshift independent, there is a clear evolution of increase $M_{min}^{lens}$ with redshift in agreement with theoretical estimations. Finally, the halo modeling allows us to identify a strong lensing contribution to the cross-correlation for angular scales below 30". This interpretation is supported by the results of basic but effective simulations.

The Herschel Space Observatory has revealed a very different galaxyscape from that shown by optical surveys which presents a challenge for galaxy-evolution models. The Herschel surveys reveal (1) that there was rapid galaxy evolution in the very recent past and (2) that galaxies lie on a a single Galaxy Sequence (GS) rather than a star-forming `main sequence' and a separate region of `passive' or `red-and-dead' galaxies. The form of the GS is now clearer because far-infrared surveys such as the Herschel ATLAS pick up a population of optically-red star-forming galaxies that would have been classified as passive using most optical criteria. The space-density of this population is at least as high as the traditional star-forming population. By stacking spectra of H-ATLAS galaxies over the redshift range 0.001 < z < 0.4, we show that the galaxies responsible for the rapid low-redshift evolution have high stellar masses, high star-formation rates but, even several billion years in the past, old stellar populations - they are thus likely to be relatively recent ancestors of early-type galaxies in the Universe today. The form of the GS is inconsistent with rapid quenching models and neither the analytic bathtub model nor the hydrodynamical EAGLE simulation can reproduce the rapid cosmic evolution. We propose a new gentler model of galaxy evolution that can explain the new Herschel results and other key properties of the galaxy population.

We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star-formation rates for over 200,000 GAMA galaxies, 170,000 G10-COSMOS galaxies and 200,000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous dataset spanning a broad range in stellar mass (10^8---10^12 Msol), dust mass (10^6---10^9 Msol), and star-formation rates (0.01---100 Msol per yr), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star-formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous dataset using consistent mass and star-formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our cosmic star-formation history we precisely reproduce the stellar-mass density with an ISM replenishment factor of 0.50 +/- 0.07, consistent with our choice of Chabrier IMF plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065---0.004 units of dust mass is also formed; (2) Over the history of the Universe approximately 90 to 95 per cent of all dust formed has been destroyed and/or ejected.

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.

We present the evolution of the Cosmic Spectral Energy Distribution (CSED) from $z = 1 - 0$. Our CSEDs originate from stacking individual spectral energy distribution fits based on panchromatic photometry from the Galaxy and Mass Assembly (GAMA) and COSMOS datasets in ten redshift intervals with completeness corrections applied. Below $z = 0.45$, we have credible SED fits from 100 nm to 1 mm. Due to the relatively low sensitivity of the far-infrared data, our far-infrared CSEDs contain a mix of predicted and measured fluxes above $z = 0.45$. Our results include appropriate errors to highlight the impact of these corrections. We show that the bolometric energy output of the Universe has declined by a factor of roughly four -- from $5.1 \pm 1.0$ at $z \sim 1$ to $1.3 \pm 0.3 \times 10^{35}~h_{70}$~W~Mpc$^{-3}$ at the current epoch. We show that this decrease is robust to cosmic variance, SED modelling and other various types of error. Our CSEDs are also consistent with an increase in the mean age of stellar populations. We also show that dust attenuation has decreased over the same period, with the photon escape fraction at 150~nm increasing from $16 \pm 3$ at $z \sim 1$ to $24 \pm 5$ per cent at the current epoch, equivalent to a decrease in $A_\mathrm{FUV}$ of 0.4~mag. Our CSEDs account for $68 \pm 12$ and $61 \pm 13$ per cent of the cosmic optical and infrared backgrounds respectively as defined from integrated galaxy counts and are consistent with previous estimates of the cosmic infrared background with redshift.

We study the environments of low- and high- excitation radio galaxies (LERGs and HERGs respectively) in the redshift range $0.01 < z < 0.4$, using a sample of 399 radio galaxies and environmental measurements from the Galaxy And Mass Assembly (GAMA) survey. In our analysis we use the fifth nearest neighbour density ($\Sigma_{5}$) and the GAMA galaxy groups catalogue (G3Cv6) and construct control samples of galaxies matched in {\update stellar mass and colour} to the radio-detected sample. We find that LERGs and HERGs exist in different environments and that this difference is dependent on radio luminosity. High-luminosity LERGs ($L_{\rm NVSS} \gtrsim 10^{24}$ W Hz$^{-1}$) lie in much denser environments than a matched radio-quiet control sample (about three times as dense, as measured by $\Sigma_{5}$), and are more likely to be members of galaxy groups ($82^{+5}_{-7}$ percent of LERGs are in GAMA groups, compared to $58^{+3}_{-3}$ percent of the control sample). In contrast, the environments of the HERGs and lower luminosity LERGs are indistinguishable from that of a matched control sample. Our results imply that high-luminosity LERGs lie in more massive haloes than non-radio galaxies of similar stellar mass and colour, in agreement with earlier studies (Wake et al. 2008; Donoso et al. 2010). When we control for the preference of LERGs to be found in groups, both high- and low- luminosity LERGs are found in higher-mass haloes ($\sim 0.2$ dex; at least 97 percent significant) than the non-radio control sample.

We derive the low redshift galaxy stellar mass function (GSMF), inclusive of dust corrections, for the equatorial Galaxy And Mass Assembly (GAMA) dataset covering 180 deg$^2$. We construct the mass function using a density-corrected maximum volume method, using masses corrected for the impact of optically thick and thin dust. We explore the galactic bivariate brightness plane ($M_\star-\mu$), demonstrating that surface brightness effects do not systematically bias our mass function measurement above 10$^{7.5}$ M$_{\odot}$. The galaxy distribution in the $M-\mu$-plane appears well bounded, indicating that no substantial population of massive but diffuse or highly compact galaxies are systematically missed due to the GAMA selection criteria. The GSMF is {fit with} a double Schechter function, with $\mathcal M^\star=10^{10.78\pm0.01\pm0.20}M_\odot$, $\phi^\star_1=(2.93\pm0.40)\times10^{-3}h_{70}^3$Mpc$^{-3}$, $\alpha_1=-0.62\pm0.03\pm0.15$, $\phi^\star_2=(0.63\pm0.10)\times10^{-3}h_{70}^3$Mpc$^{-3}$, and $\alpha_2=-1.50\pm0.01\pm0.15$. We find the equivalent faint end slope as previously estimated using the GAMA-I sample, although we find a higher value of $\mathcal M^\star$. Using the full GAMA-II sample, we are able to fit the mass function to masses as low as $10^{7.5}$ $M_\odot$, and assess limits to $10^{6.5}$ $M_\odot$. Combining GAMA-II with data from G10-COSMOS we are able to comment qualitatively on the shape of the GSMF down to masses as low as $10^{6}$ $M_\odot$. Beyond the well known upturn seen in the GSMF at $10^{9.5}$ the distribution appears to maintain a single power-law slope from $10^9$ to $10^{6.5}$. We calculate the stellar mass density parameter given our best-estimate GSMF, finding $\Omega_\star= 1.66^{+0.24}_{-0.23}\pm0.97 h^{-1}_{70} \times 10^{-3}$, inclusive of random and systematic uncertainties.

We study the stellar-to-halo mass relation of central galaxies in the range 9.7<log_10(M_*/h^-2 M_sun)<11.7 and z<0.4, obtained from a combined analysis of the Kilo Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey. We use ~100 deg^2 of KiDS data to study the lensing signal around galaxies for which spectroscopic redshifts and stellar masses were determined by GAMA. We show that lensing alone results in poor constraints on the stellar-to-halo mass relation due to a degeneracy between the satellite fraction and the halo mass, which is lifted when we simultaneously fit the stellar mass function. At M_sun>5x10^10 h^-2 M_sun, the stellar mass increases with halo mass as ~M_h^0.25. The ratio of dark matter to stellar mass has a minimum at a halo mass of 8x10^11 h^-1 M_sun with a value of M_h/M_*=56_-10^+16 [h]. We also use the GAMA group catalogue to select centrals and satellites in groups with five or more members, which trace regions in space where the local matter density is higher than average, and determine for the first time the stellar-to-halo mass relation in these denser environments. We find no significant differences compared to the relation from the full sample, which suggests that the stellar-to-halo mass relation does not vary strongly with local density. Furthermore, we find that the stellar-to-halo mass relation of central galaxies can also be obtained by modelling the lensing signal and stellar mass function of satellite galaxies only, which shows that the assumptions to model the satellite contribution in the halo model do not significantly bias the stellar-to-halo mass relation. Finally, we show that the combination of weak lensing with the stellar mass function can be used to test the purity of group catalogues.

The merging history of galaxies can be traced with studies of dynamically close pairs. These consist of a massive primary galaxy and a less massive secondary (or satellite) galaxy. The study of the stellar populations of secondary (lower mass) galaxies in close pairs provides a way to understand galaxy growth by mergers. Here we focus on systems involving at least one massive galaxy - with stellar mass above $10^{11}M_\odot$ in the highly complete GAMA survey. Our working sample comprises 2,692 satellite galaxy spectra (0.1<z<0.3). These spectra are combined into high S/N stacks, and binned according to both an "internal" parameter, the stellar mass of the satellite galaxy (i.e. the secondary), and an "external" parameter, selecting either the mass of the primary in the pair, or the mass of the corresponding dark matter halo. We find significant variations in the age of the populations with respect to environment. At fixed mass, satellites around the most massive galaxies are older and possibly more metal rich, with age differences ~1-2Gyr within the subset of lower mass satellites ($\sim 10^{10}M_\odot$). These variations are similar when stacking with respect to the halo mass of the group where the pair is embedded. The population trends in the lower-mass satellites are consistent with the old stellar ages found in the outer regions of massive galaxies.

We present a robust calibration of the 1.4GHz radio continuum star formation rate (SFR) using a combination of the Galaxy And Mass Assembly (GAMA) survey and the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. We identify individually detected 1.4GHz GAMA-FIRST sources and use a late-type, non-AGN, volume-limited sample from GAMA to produce stellar mass-selected samples. The latter are then combined to produce FIRST-stacked images. This extends the robust parametrisation of the 1.4GHz-SFR relation to faint luminosities. For both the individually detected galaxies and our stacked samples, we compare 1.4GHz luminosity to SFRs derived from GAMA to determine a new 1.4GHz luminosity-to-SFR relation with well constrained slope and normalisation. For the first time, we produce the radio SFR-M* relation over 2 decades in stellar mass, and find that our new calibration is robust, and produces a SFR-M* relation which is consistent with all other GAMA SFR methods. Finally, using our new 1.4GHz luminosity-to-SFR calibration we make predictions for the number of star-forming GAMA sources which are likely to be detected in the upcoming ASKAP surveys, EMU and DINGO.

We present an analysis of the mid-infrared WISE sources seen within the equatorial GAMA G12 field, located in the North Galactic Cap. Our motivation is to study and characterize the behavior of WISE source populations in anticipation of the deep multi-wavelength surveys that will define the next decade, with the principal science goal of mapping the 3D large scale structures and determining the global physical attributes of the host galaxies. In combination with cosmological redshifts, we identify galaxies from their WISE W1 3.4um resolved emission, and by performing a star-galaxy separation using apparent magnitude, colors and statistical modeling of star-counts. The resultant galaxy catalog has ~590,000 sources in 60 deg^2, reaching a W1 5-sigma depth of 31 uJy. At the faint end, where redshifts are not available, we employ a luminosity function analysis to show that approximately 27% of all WISE extragalactic sources to a limit of 17.5 mag (31 uJy) are at high redshift, z > 1. The spatial distribution is investigated using two-point correlation functions and a 3D source density characterization at 5 Mpc and 20 Mpc scales. For angular distributions, we find brighter and more massive sources are strongly clustered relative to fainter and lower mass source; likewise, based on WISE colors, spheroidal galaxies have the strongest clustering, while late-type disk galaxies have the lowest clustering amplitudes. Along the radial direction, the strongest clustering is in the largest redshift shell, while the weakest is in the nearest redshift shell, consistent with the stellar mass and morphological type dependency results. In three dimensions, we find a number of distinct groupings, often bridged by filaments and super-structures. Using special visualization tools, we map these structures, exploring how clustering may play a role with stellar mass and galaxy type.

Abridged - We quantify the effect of the galaxy group environment (for 12.5 < log(M_group/Msun) < 14.0) on the star formation rates of the (morphologically-selected) population of disk-dominated local Universe spiral galaxies (z < 0.13) with stellar masses log(M*/Msun) > 9.5. Within this population, we find that, while a small minority of group satellites are strongly quenched, the group centrals, and the large majority of satellites exhibit levels of SFR indistinguishable from ungrouped "field" galaxies of the same M*, albeit with a higher scatter, and for all M*. Modelling these results, we deduce that disk-dominated satellites continue to be characterized by a rapid cycling of gas into and out of their ISM at rates similar to those operating prior to infall, with the on-going fuelling likely sourced from the group intrahalo medium (IHM) on Mpc scales, rather than from the circum-galactic medium on 100kpc scales. Consequently, the color-density relation of the galaxy population as a whole would appear to be primarily due to a change in the mix of disk- and spheroid-dominated morphologies in the denser group environment compared to the field, rather than to a reduced propensity of the IHM in higher mass structures to cool and accrete onto galaxies. We also suggest that the inferred substantial accretion of IHM gas by satellite disk-dominated galaxies will lead to a progressive reduction in their specific angular momentum, thereby representing an efficient secular mechanism to transform morphology from star-forming disk-dominated types to more passive spheroid-dominated types.

Verlinde (2016) proposed that the observed excess gravity in galaxies and clusters is the consequence of Emergent Gravity (EG). In this theory the standard gravitational laws are modified on galactic and larger scales due to the displacement of dark energy by baryonic matter. EG gives an estimate of the excess gravity (described as an apparent dark matter density) in terms of the baryonic mass distribution and the Hubble parameter. In this work we present the first test of EG using weak gravitational lensing, within the regime of validity of the current model. Although there is no direct description of lensing and cosmology in EG yet, we can make a reasonable estimate of the expected lensing signal of low redshift galaxies by assuming a background LambdaCDM cosmology. We measure the (apparent) average surface mass density profiles of 33,613 isolated central galaxies, and compare them to those predicted by EG based on the galaxies' baryonic masses. To this end we employ the ~180 square degrees overlap of the Kilo-Degree Survey (KiDS) with the spectroscopic Galaxy And Mass Assembly (GAMA) survey. We find that the prediction from EG, despite requiring no free parameters, is in good agreement with the observed galaxy-galaxy lensing profiles in four different stellar mass bins. Although this performance is remarkable, this study is only a first step. Further advancements on both the theoretical framework and observational tests of EG are needed before it can be considered a fully developed and solidly tested theory.

We present predictions for the galaxy-galaxy lensing profile from the EAGLE hydrodynamical cosmological simulation at redshift z=0.18, in the spatial range 0.02 < R/(Mpc/h) < 2, and for five logarithmically equi-spaced stellar mass bins in the range 10.3 < $\log_{10}$(Mstar/ $M_{\odot}$) < 11.8. We compare these excess surface density profiles to the observed signal from background galaxies imaged by the Kilo Degree Survey around spectroscopically confirmed foreground galaxies from the GAMA survey. Exploiting the GAMA galaxy group catalogue, the profiles of central and satellite galaxies are computed separately for groups with at least five members to minimise contamination. EAGLE predictions are in broad agreement with the observed profiles for both central and satellite galaxies, although the signal is underestimated at R$\approx$0.5-2 Mpc/h for the highest stellar mass bins. When central and satellite galaxies are considered simultaneously, agreement is found only when the selection function of lens galaxies is taken into account in detail. Specifically, in the case of GAMA galaxies, it is crucial to account for the variation of the fraction of satellite galaxies in bins of stellar mass induced by the flux-limited nature of the survey. We report the inferred stellar-to-halo mass relation and we find good agreement with recent published results. We note how the precision of the galaxy-galaxy lensing profiles in the simulation holds the potential to constrain fine-grained aspects of the galaxy-dark matter connection.

There exist conflicting observations on whether or not the environment of broad and narrow line AGN differ and this consequently questions the validity of the AGN unification model. The high spectroscopic completeness of the GAMA survey makes it ideal for a comprehensive analysis of the close environment of galaxies. To exploit this, and conduct a comparative analysis of the environment of broad and narrow line AGN within GAMA, we use a double-Gaussian emission line fitting method to model the more complex line profiles associated with broad line AGN. We select 209 type 1 (i.e., unobscured), 464 type 1.5-1.9 (partially obscured), and 281 type 2 (obscured) AGN within the GAMA II database. Comparing the fractions of these with neighbouring galaxies out to a pair separation of $350\,\text{kpc }h^{-1}$ and $\Delta z < 0.012$ shows no difference between AGN of different type, except at separations less than $20\,\text{kpc }h^{-1}$ where our observations suggest an excess of type 2 AGN in close pairs. We analyse the properties of the galaxies neighbouring our AGN and find no significant differences in colour or the star formation activity of these galaxies. Further to this we find that $\Sigma_5$ is also consistent between broad and narrow line AGN. We conclude that the observations presented here are consistent with AGN unification.

We investigate the contentious issue of the presence, or lack thereof, of satellites mass segregation in galaxy groups using the Galaxy And Mass Assembly (GAMA) survey, the GALFORM semi-analytic and the EAGLE cosmological hydrodynamical simulation catalogues of galaxy groups. We select groups with halo mass $12 \leqslant \log(M_{\text{halo}}/h^{-1}M_\odot) <14.5$ and redshift $z \leqslant 0.32$ and probe the radial distribution of stellar mass out to twice the group virial radius. All the samples are carefully constructed to be complete in stellar mass at each redshift range and efforts are made to regularise the analysis for all the data. Our study shows negligible mass segregation in galaxy group environments with absolute gradients of $\lesssim0.08$ dex and also shows a lack of any redshift evolution. Moreover, we find that our results at least for the GAMA data are robust to different halo mass and group centre estimates. Furthermore, the EAGLE data allows us to probe much fainter luminosities ($r$-band magnitude of 22) as well as investigate the three-dimensional spatial distribution with intrinsic halo properties, beyond what the current observational data can offer. In both cases we find that the fainter EAGLE data show a very mild spatial mass segregation at $z \leqslant 0.22$, which is again not apparent at higher redshift. Interestingly, our results are in contrast to some earlier findings using the Sloan Digital Sky Survey. We investigate the source of the disagreement and suggest that subtle differences between the group finding algorithms could be the root cause.

We compare common star-formation rate (SFR) indicators in the local Universe in the GAMA equatorial fields (around 160 sq. deg.), using ultraviolet (UV) photometry from GALEX, far-infrared (FIR) and sub-millimetre (sub-mm) photometry from H-ATLAS, and Halpha spectroscopy from the GAMA survey. With a high-quality sample of 745 galaxies (median redshift 0.08), we consider three SFR tracers: UV luminosity corrected for dust attenuation using the UV spectral slope beta (SFRUV,corr), Halpha line luminosity corrected for dust using the Balmer decrement (BD) (SFRHalpha,corr), and the combination of UV and IR emission (SFRUV+IR). We demonstrate that SFRUV,corr can be reconciled with the other two tracers after applying attenuation corrections by calibrating IRX (i.e. the IR to UV luminosity ratio) and attenuation in the Halpha (derived from BD) against beta. However, beta on its own is very unlikely to be a reliable attenuation indicator. We find that attenuation correction factors depend on parameters such as stellar mass, z and dust temperature (Tdust), but not on Halpha equivalent width (EW) or Sersic index. Due to the large scatter in the IRX vs beta correlation, when compared to SFRUV+IR, the beta-corrected SFRUV,corr exhibits systematic deviations as a function of IRX, BD and Tdust.

We cross-match the two currently largest all-sky photometric catalogs, mid-infrared WISE and SuperCOSMOS scans of UKST/POSS-II photographic plates, to obtain a new galaxy sample that covers 3pi steradians. In order to characterize and purify the extragalactic dataset, we use external GAMA and SDSS spectroscopic information to define quasar and star loci in multicolor space, aiding the removal of contamination from our extended-source catalog. After appropriate data cleaning we obtain a deep wide-angle galaxy sample that is approximately 95% pure and 90% complete at high Galactic latitudes. The catalog contains close to 20 million galaxies over almost 70% of the sky, outside the Zone of Avoidance and other confused regions, with a mean surface density of over 650 sources per square degree. Using multiwavelength information from two optical and two mid-IR photometric bands, we derive photometric redshifts for all the galaxies in the catalog, using the ANNz framework trained on the final GAMA-II spectroscopic data. Our sample has a median redshift of z_{med} = 0.2 but with a broad dN/dz reaching up to z>0.4. The photometric redshifts have a mean bias of |delta_z|~10^{-3}, normalized scatter of sigma_z = 0.033 and less than 3% outliers beyond 3sigma_z. Comparison with external datasets shows no significant variation of photo-z quality with sky position. Together with the overall statistics, we also provide a more detailed analysis of photometric redshift accuracy as a function of magnitudes and colors. The final catalog is appropriate for `all-sky' 3D cosmology to unprecedented depths, in particular through cross-correlations with other large-area surveys. It should also be useful for source pre-selection and identification in forthcoming surveys such as TAIPAN or WALLABY.

Aims. We aim to study the 250 micron luminosity function (LF) down to much fainter luminosities than achieved by previous efforts. Methods. We developed a modified stacking method to reconstruct the 250 micron LF using optically selected galaxies from the SDSS survey and Herschel maps of the GAMA equatorial fields and Stripe 82. Our stacking method not only recovers the mean 250 micron luminosities of galaxies that are too faint to be individually detected, but also their underlying distribution functions. Results. We find very good agreement with previous measurements in the overlapping luminosity range. More importantly, we are able to derive the LF down to much fainter luminosities (around 25 times fainter) than achieved by previous studies. We find strong positive luminosity evolution \propto (1 + z)^4.89\pm1.07 and moderate negative density evolution \propto (1 + z)^-1.02\pm0.54 over the redshift range z=[0.02, 0.5].

Using the L-Galaxies semi-analytic model we simultaneously fit the HI mass function, stellar mass function and galaxy colours. We find good fits to all three observations at z = 0 and to the stellar mass function and galaxy colours at z = 2. Using Markov Chain Monte Carlo (MCMC) techniques we adjust the L-Galaxies parameters to best fit the constraining data. In order to fit the HI mass function we must greatly reduce the gas surface density threshold for star formation, thus lowering the number of low HI mass galaxies. A simultaneous reduction in the star formation efficiency prevents the over production of stellar content. A simplified model in which the surface density threshold is eliminated altogether also provides a good fit to the data. Unfortunately, these changes weaken the fit to the Kennicutt-Schmidt relation and raise the star-formation rate density at recent times, suggesting that a change to the model is required to prevent accumulation of gas onto dwarf galaxies in the local universe.

We present a meta-analysis of star-formation rate (SFR) indicators in the GAMA survey, producing 12 different SFR metrics and determining the SFR-M* relation for each. We compare and contrast published methods to extract the SFR from each indicator, using a well-defined local sample of morphologically-selected spiral galaxies, which excludes sources which potentially have large recent changes to their SFR. The different methods are found to yield SFR-M* relations with inconsistent slopes and normalisations, suggesting differences between calibration methods. The recovered SFR-M* relations also have a large range in scatter which, as SFRs of the targets may be considered constant over the different timescales, suggests differences in the accuracy by which methods correct for attenuation in individual targets. We then recalibrate all SFR indicators to provide new, robust and consistent luminosity-to-SFR calibrations, finding that the most consistent slopes and normalisations of the SFR-M* relations are obtained when recalibrated using the radiation transfer method of Popescu et al. These new calibrations can be used to directly compare SFRs across different observations, epochs and galaxy populations. We then apply our calibrations to the GAMA II equatorial dataset and explore the evolution of star-formation in the local Universe. We determine the evolution of the normalisation to the SFR-M* relation from 0 < z < 0.35 - finding consistent trends with previous estimates at 0.3 < z < 1.2. We then provide the definitive z < 0.35 Cosmic Star Formation History, SFR-M* relation and its evolution over the last 3 billion years.

We present radio Active Galactic Nuclei (AGN) luminosity functions over the redshift range 0.005 < z < 0.75. The sample from which the luminosity functions are constructed is an optical spectroscopic survey of radio galaxies, identified from matched Faint Images of the Radio Sky at Twenty-cm survey (FIRST) sources and Sloan Digital Sky Survey (SDSS) images.The radio AGN are separated into Low Excitation Radio Galaxies (LERGs) and High Excitation Radio Galaxies (HERGs) using the optical spectra. We derive radio luminosity functions for LERGs and HERGs separately in the three redshift bins (0.005 < z < 0.3, 0.3 < z < 0.5 and 0.5 < z <0.75). The radio luminosity functions can be well described by a double power-law. Assuming this double power-law shape the LERG population displays little or no evolution over this redshift range evolving as ~$(1+z)^{0.06}$ assuming pure density evolution or ~ $(1+z)^{0.46}$ assuming pure luminosity evolution. In contrast, the HERG population evolves more rapidly, best fitted by ~$(1+z)^{2.93}$ assuming a double power-law shape and pure density evolution. If a pure luminosity model is assumed the best fitting HERG evolution is parameterised by ~$(1+z)^{7.41}$. The characteristic break in the radio luminosity function occurs at a significantly higher power (~1 dex) for the HERG population in comparison to the LERGs. This is consistent with the two populations representing fundamentally different accretion modes.

We present the Lambda Adaptive Multi-Band Deblending Algorithm in R (LAMBDAR), a novel code for calculating matched aperture photometry across images that are neither pixel- nor PSF-matched, using prior aperture definitions derived from high resolution optical imaging. The development of this program is motivated by the desire for consistent photometry and uncertainties across large ranges of photometric imaging, for use in calculating spectral energy distributions. We describe the program, specifically key features required for robust determination of panchromatic photometry: propagation of apertures to images with arbitrary resolution, local background estimation, aperture normalisation, uncertainty determination and propagation, and object deblending. Using simulated images, we demonstrate that the program is able to recover accurate photometric measurements in both high-resolution, low-confusion, and low-resolution, high-confusion, regimes. We apply the program to the 21-band photometric dataset from the Galaxy And Mass Assembly (GAMA) Panchromatic Data Release (PDR; Driver et al. 2016), which contains imaging spanning the far-UV to the far-IR. We compare photometry derived from LAMBDAR with that presented in Driver et al. (2016), finding broad agreement between the datasets. Nonetheless, we demonstrate that the photometry from LAMBDAR is superior to that from the GAMA PDR, as determined by a reduction in the outlier rate and intrinsic scatter of colours in the LAMBDAR dataset. We similarly find a decrease in the outlier rate of stellar masses and star formation rates using LAMBDAR photometry. Finally, we note an exceptional increase in the number of UV and mid-IR sources able to be constrained, which is accompanied by a significant increase in the mid-IR colour-colour parameter-space able to be explored.

We look for correlated changes in stellar mass and star formation rate along filaments in the cosmic web by examining the stellar masses and UV-derived star formation rates (SFR) of 1,799 ungrouped and unpaired spiral galaxies that reside in filaments. We devise multiple distance metrics to characterise the complex geometry of filaments, and find that galaxies closer to the cylindrical centre of a filament have higher stellar masses than their counterparts near the periphery of filaments, on the edges of voids. In addition, these peripheral spiral galaxies have higher specific star formation rates (SSFR) at a given mass. Complementing our sample of filament spiral galaxies with spiral galaxies in tendrils and voids, we find that the average SFR of these objects in different large scale environments are similar to each other with the primary discriminant in SFR being stellar mass, in line with previous works. However, the distributions of SFRs are found to vary with large-scale environment. Our results thus suggest a model in which in addition to stellar mass as the primary discriminant, the large-scale environment is imprinted in the SFR as a second order effect. Furthermore, our detailed results for filament galaxies suggest a model in which gas accretion from voids onto filaments is primarily in an orthogonal direction. Overall, we find our results to be in line with theoretical expectations of the thermodynamic properties of the intergalactic medium in different large-scale environments.

We present the first cosmological measurement derived from a galaxy density field subject to a `clipping' transformation. By enforcing an upper bound on the galaxy number density field in the Galaxy and Mass Assembly survey (GAMA), contributions from the nonlinear processes of virialisation and galaxy bias are greatly reduced. This leads to a galaxy power spectrum which is easier to model, without calibration from numerical simulations. We develop a theoretical model for the power spectrum of a clipped field in redshift space, which is exact for the case of anisotropic Gaussian fields. Clipping is found to extend the applicability of the conventional Kaiser prescription by more than a factor of three in wavenumber, or a factor of thirty in terms of the number of Fourier modes. By modelling the galaxy power spectrum on scales k < 0.3 h/Mpc and density fluctuations $\delta_g < 4$ we measure the normalised growth rate $f\sigma_8(z = 0.18) = 0.29 \pm 0.10$.

Measurement of the evolution of both active galactic nuclei (AGN) and star-formation in galaxies underpins our understanding of galaxy evolution over cosmic time. Radio continuum observations can provide key information on these two processes, in particular via the mechanical feedback produced by radio jets in AGN, and via an unbiased dust-independent measurement of star-formation rates. In this paper we determine radio luminosity functions at 325 MHz for a sample of AGN and star-forming galaxies by matching a 138 deg sq. radio survey conducted with the Giant Metrewave Radio Telescope (GMRT), with optical imaging and redshifts from the Galaxy And Mass Assembly (GAMA) survey. We find that the radio luminosity function at 325 MHz for star-forming galaxies closely follows that measured at 1.4 GHz. By fitting the AGN radio luminosity function out to $z = 0.5$ as a double power law, and parametrizing the evolution as ${\Phi} \propto (1 + z)^{k}$ , we find evolution parameters of $k = 0.92 \pm 0.95$ assuming pure density evolution and $k = 2.13 \pm 1.96$ assuming pure luminosity evolution. We find that the Low Excitation Radio Galaxies are the dominant population in space density at lower luminosities. Comparing our 325 MHz observations with radio continuum imaging at 1.4 GHz, we determine separate radio luminosity functions for steep and flat-spectrum AGN, and show that the beamed population of flat-spectrum sources in our sample can be shifted in number density and luminosity to coincide with the unbeamed population of steep-spectrum sources, as is expected in the orientation based unification of AGN.

We combine Herschel/SPIRE sub-millimeter (submm) observations with existing multi-wavelength data to investigate the characteristics of low redshift, optically red galaxies detected in submm bands. We select a sample of galaxies in the redshift range 0.01$\leq$z$\leq$0.2, having >5$\sigma$ detections in the SPIRE 250 micron submm waveband. Sources are then divided into two sub-samples of $red$ and $blue$ galaxies, based on their UV-optical colours. Galaxies in the $red$ sample account for $\approx$4.2 per cent of the total number of sources with stellar masses M$_{*}\gtrsim$10$^{10}$ Solar-mass. Following visual classification of the $red$ galaxies, we find that $\gtrsim$30 per cent of them are early-type galaxies and $\gtrsim$40 per cent are spirals. The colour of the $red$-spiral galaxies could be the result of their highly inclined orientation and/or a strong contribution of the old stellar population. It is found that irrespective of their morphological types, $red$ and $blue$ sources occupy environments with more or less similar densities (i.e., the $\Sigma_5$ parameter). From the analysis of the spectral energy distributions (SEDs) of galaxies in our samples based on MAGPHYS, we find that galaxies in the $red$ sample (of any morphological type) have dust masses similar to those in the $blue$ sample (i.e. normal spiral/star-forming systems). However, in comparison to the $red$-spirals and in particular $blue$ systems, $red$-ellipticals have lower mean dust-to-stellar mass ratios. Besides galaxies in the $red$-elliptical sample have much lower mean star-formation/specific-star-formation rates in contrast to their counterparts in the $blue$ sample. Our results support a scenario where dust in early-type systems is likely to be of an external origin.

We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to the far-IR. These data complement our spectroscopic campaign of over 300k galaxies, and are compiled from observations with a variety of facilities including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by VST and scheduled for observations by ASKAP. These data are processed to a common astrometric solution, from which photometry is derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to access and download data cutouts, or the full mosaics of the GAMA regions in each band. We focus, in particular, on the reduction and analysis of the VISTA VIKING data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the data and examining its integrity. Having derived the 21-band photometric catalogue we proceed to fit the data using the energy balance code MAGPHYS. These measurements are then used to obtain the first fully empirical measurement of the 0.1-500$\mu$m energy output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) across three time-intervals (0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More importantly, we identify significant and smooth evolution in the integrated photon escape fraction at all wavelengths, with the UV escape fraction increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06. The GAMA PDR will allow for detailed studies of the energy production and outputs of individual systems, sub-populations, and representative galaxy samples at $z<0.5$. The GAMA PDR can be found at: http://gama-psi.icrar.org/

We measure the projected 2-point correlation function of galaxies in the 180 deg$^2$ equatorial regions of the GAMA II survey, for four different redshift slices between z = 0.0 and z=0.5. To do this we further develop the Cole (2011) method of producing suitable random catalogues for the calculation of correlation functions. We find that more r-band luminous, more massive and redder galaxies are more clustered. We also find that red galaxies have stronger clustering on scales less than ~3 $h^{-1}$ Mpc. We compare to two different versions of the GALFORM galaxy formation model, Lacey et al (in prep.) and Gonzalez-Perez et al. (2014), and find that the models reproduce the trend of stronger clustering for more massive galaxies. However, the models under predict the clustering of blue galaxies, can incorrectly predict the correlation function on small scales and under predict the clustering in our sample of galaxies with ~3$L_r$ . We suggest possible avenues to explore to improve these cluster- ing predictions. The measurements presented in this paper can be used to test other galaxy formation models, and we make the measurements available online to facilitate this.

We examine the properties of galaxies in the Galaxies and Mass Assembly (GAMA) survey located in voids with radii $>10~h^{-1}$ Mpc. Utilising the GAMA equatorial survey, 592 void galaxies are identified out to z~0.1 brighter than $M_{r} = -18.4$, our magnitude completeness limit. Using the $W_{\rm{H\alpha}}$ vs. [NII]/H$\alpha$ (WHAN) line strength diagnostic diagram, we classify their spectra as star forming, AGN, or dominated by old stellar populations. For objects more massive than $5\times10^{9}$ M$_{\odot}$, we identify a sample of 26 void galaxies with old stellar populations classed as passive and retired galaxies in the WHAN diagnostic diagram, else they lack any emission lines in their spectra. When matched to WISE mid-IR photometry, these passive and retired galaxies exhibit a range of mid-IR colour, with a number of void galaxies exhibiting [4.6]-[12] colours inconsistent with completely quenched stellar populations, with a similar spread in colour seen for a randomly drawn non-void comparison sample. We hypothesise that a number of these galaxies host obscured star formation, else they are star forming outside of their central regions targeted for single fibre spectroscopy. When matched to a randomly drawn sample of non-void galaxies, the void and non-void galaxies exhibit similar properties in terms of optical and mid-IR colour, morphology, and star formation activity, suggesting comparable mass assembly and quenching histories. A trend in mid-IR [4.6]-[12] colour is seen, such that both void and non-void galaxies with quenched/passive colours <1.5 typically have masses higher than $10^{10}$ M$_{\odot}$, where internally driven processes play an increasingly important role in galaxy evolution.

The modification of star formation (SF) in galaxy interactions is a complex process, with SF observed to be both enhanced in major mergers and suppressed in minor pair interactions. Such changes likely to arise on short timescales and be directly related to the galaxy-galaxy interaction time. Here we investigate the link between dynamical phase and direct measures of SF on different timescales for pair galaxies, targeting numerous star-formation rate (SFR) indicators and comparing to pair separation, individual galaxy mass and pair mass ratio. We split our sample into the higher (primary) and lower (secondary) mass galaxies in each pair and find that SF is indeed enhanced in all primary galaxies but suppressed in secondaries of minor mergers. We find that changes in SF of primaries is consistent in both major and minor mergers, suggesting that SF in the more massive galaxy is agnostic to pair mass ratio. We also find that SF is enhanced/suppressed more strongly for short-time duration SFR indicators (e.g. H-alpha), highlighting recent changes to SF in these galaxies, which are likely to be induced by the interaction. We propose a scenario where the lower mass galaxy has its SF suppressed by gas heating or stripping, while the higher mass galaxy has its SF enhanced, potentially by tidal gas turbulence and shocks. This is consistent with the seemingly contradictory observations for both SF suppression and enhancement in close pairs.

The Kilo-Degree Survey (KiDS) is an optical wide-field survey designed to map the matter distribution in the Universe using weak gravitational lensing. In this paper, we use these data to measure the density profiles and masses of a sample of $\sim \mathrm{1400}$ spectroscopically identified galaxy groups and clusters from the Galaxy And Mass Assembly (GAMA) survey. We detect a highly significant signal (signal-to-noise-ratio $\sim$ 120), allowing us to study the properties of dark matter haloes over one and a half order of magnitude in mass, from $M \sim 10^{13}-10^{14.5} h^{-1}\mathrm{M_{\odot}}$. We interpret the results for various subsamples of groups using a halo model framework which accounts for the mis-centring of the Brightest Cluster Galaxy (used as the tracer of the group centre) with respect to the centre of the group's dark matter halo. We find that the density profiles of the haloes are well described by an NFW profile with concentrations that agree with predictions from numerical simulations. In addition, we constrain scaling relations between the mass and a number of observable group properties. We find that the mass scales with the total r-band luminosity as a power-law with slope $1.16 \pm 0.13$ (1-sigma) and with the group velocity dispersion as a power-law with slope $1.89 \pm 0.27$ (1-sigma). Finally, we demonstrate the potential of weak lensing studies of groups to discriminate between models of baryonic feedback at group scales by comparing our results with the predictions from the Cosmo-OverWhelmingly Large Simulations (Cosmo-OWLS) project, ruling out models without AGN feedback.

The Galaxy And Mass Assembly (GAMA) survey is one of the largest contemporary spectroscopic surveys of low-redshift galaxies. Covering an area of ~286 deg^2 (split among five survey regions) down to a limiting magnitude of r < 19.8 mag, we have collected spectra and reliable redshifts for 238,000 objects using the AAOmega spectrograph on the Anglo-Australian Telescope. In addition, we have assembled imaging data from a number of independent surveys in order to generate photometry spanning the wavelength range 1 nm - 1 m. Here we report on the recently completed spectroscopic survey and present a series of diagnostics to assess its final state and the quality of the redshift data. We also describe a number of survey aspects and procedures, or updates thereof, including changes to the input catalogue, redshifting and re-redshifting, and the derivation of ultraviolet, optical and near-infrared photometry. Finally, we present the second public release of GAMA data. In this release we provide input catalogue and targeting information, spectra, redshifts, ultraviolet, optical and near-infrared photometry, single-component S\'ersic fits, stellar masses, H$\alpha$-derived star formation rates, environment information, and group properties for all galaxies with r < 19.0 mag in two of our survey regions, and for all galaxies with r < 19.4 mag in a third region (72,225 objects in total). The database serving these data is available at http://www.gama-survey.org/.

Using results from the Herschel Astrophysical Terrahertz Large-Area Survey and the Galaxy and Mass Assembly project, we show that, for galaxy masses above approximately 1.0e8 solar masses, 51% of the stellar mass-density in the local Universe is in early-type galaxies (ETGs: Sersic n > 2.5) while 89% of the rate of production of stellar mass-density is occurring in late-type galaxies (LTGs: Sersic n < 2.5). From this zero-redshift benchmark, we have used a calorimetric technique to quantify the importance of the morphological transformation of galaxies over the history of the Universe. The extragalactic background radiation contains all the energy generated by nuclear fusion in stars since the Big Bang. By resolving this background radiation into individual galaxies using the deepest far-infrared survey with the Herschel Space Observatory and a deep near-infrared/optical survey with the Hubble Space Telescope (HST), and using measurements of the Sersic index of these galaxies derived from the HST images, we estimate that approximately 83% of the stellar mass-density formed over the history of the Universe occurred in LTGs. The difference between this and the fraction of the stellar mass-density that is in LTGs today implies there must have been a major transformation of LTGs into ETGs after the formation of most of the stars.

We present the catalogue of blended galaxy spectra from the Galaxy And Mass Assembly (GAMA) survey. These are cases where light from two galaxies are significantly detected in a single GAMA fibre. Galaxy pairs identified from their blended spectrum fall into two principal classes: they are either strong lenses, a passive galaxy lensing an emission-line galaxy; or occulting galaxies, serendipitous overlaps of two galaxies, of any type. Blended spectra can thus be used to reliably identify strong lenses for follow-up observations (high resolution imaging) and occulting pairs, especially those that are a late-type partly obscuring an early-type galaxy which are of interest for the study of dust content of spiral and irregular galaxies. The GAMA survey setup and its autoz automated redshift determination were used to identify candidate blended galaxy spectra from the cross-correlation peaks. We identify 280 blended spectra with a minimum velocity separation of 600 km/s, of which 104 are lens pair candidates, 71 emission-line-passive pairs, 78 are pairs of emission-line galaxies and and 27 are pairs of galaxies with passive spectra. We have visually inspected the candidates in the Sloan Digital Sky Survey (SDSS) and Kilo Degree Survey (KiDS) images. Many blended objects are ellipticals with blue fuzz (Ef in our classification). These latter "Ef" classifications are candidates for possible strong lenses, massive ellipticals with an emission-line galaxy in one or more lensed images. The GAMA lens and occulting galaxy candidate samples are similar in size to those identified in the entire SDSS. This blended spectrum sample stands as a testament of the power of this highly complete, second-largest spectroscopic survey in existence and offers the possibility to expand e.g., strong gravitational lens surveys.

We explore trends in galaxy properties with Mpc-scale structures using catalogues of environment and large scale structure from the Galaxy And Mass Assembly (GAMA) survey. Existing GAMA catalogues of large scale structure, group and pair membership allow us to construct galaxy stellar mass functions for different environmental types. To avoid simply extracting the known underlying correlations between galaxy properties and stellar mass, we create a mass matched sample of galaxies with stellar masses between $9.5 \leq \log{M_*/h^{-2} M_{\odot}} \leq 11$ for each environmental population. Using these samples, we show that mass normalised galaxies in different large scale environments have similar energy outputs, $u-r$ colours, luminosities, and morphologies. Extending our analysis to group and pair environments, we show galaxies that are not in groups or pairs exhibit similar characteristics to each other regardless of broader environment. For our mass controlled sample, we fail to see a strong dependence of S\'{e}rsic index or galaxy luminosity on halo mass, but do find that it correlates very strongly with colour. Repeating our analysis for galaxies that have not been mass controlled introduces and amplifies trends in the properties of galaxies in pairs, groups, and large scale structure, indicating that stellar mass is the most important predictor of the galaxy properties we examine, as opposed to environmental classifications.

We describe modifications to the joint stepwise maximum likelihood method of Cole (2011) in order to simultaneously fit the GAMA-II galaxy luminosity function (LF), corrected for radial density variations, and its evolution with redshift. The whole sample is reasonably well-fit with luminosity (Qe) and density (Pe) evolution parameters Qe, Pe = 1.0, 1.0 but with significant degeneracies characterized by Qe = 1.4 - 0.4Pe. Blue galaxies exhibit larger luminosity density evolution than red galaxies, as expected. We present the evolution-corrected r-band LF for the whole sample and for blue and red sub-samples, using both Petrosian and Sersic magnitudes. Petrosian magnitudes miss a substantial fraction of the flux of de Vaucouleurs profile galaxies: the Sersic LF is substantially higher than the Petrosian LF at the bright end.

We investigate the dependence of the galaxy luminosity function on geometric environment within the Galaxy And Mass Assembly (GAMA) survey. The tidal tensor prescription, based on the Hessian of the pseudo-gravitational potential, is used to classify the cosmic web and define the geometric environments: for a given smoothing scale, we classify every position of the surveyed region, $0.04<{z}<0.26$, as either a void, a sheet, a filament or a knot. We consider how to choose appropriate thresholds in the eigenvalues of the Hessian in order to partition the galaxies approximately evenly between environments. We find a significant variation in the luminosity function of galaxies between different geometric environments; the normalisation, characterised by $\phi^{*}$ in a Schechter function fit, increases by an order of magnitude from voids to knots. The turnover magnitude, characterised by $M^*$, brightens by approximately $0.5$ mag from voids to knots. However, we show that the observed modulation can be entirely attributed to the indirect local-density dependence. We therefore find no evidence of a direct influence of the cosmic web on the galaxy luminosity function.

We use data from the Galaxy And Mass Assembly (GAMA) survey in the redshift range 0.01$<$z$<$0.1 (8399 galaxies in $g$ to $K_s$ bands) to derive the stellar mass $-$ half-light radius relations for various divisions of 'early' and 'late'-type samples. We find the choice of division between early and late (i.e., colour, shape, morphology) is not particularly critical, however, the adopted mass limits and sample selections (i.e., the careful rejection of outliers and use of robust fitting methods) are important. In particular we note that for samples extending to low stellar mass limits ($<10^{10}\mathcal{M_{\odot}}$) the S\'ersic index bimodality, evident for high mass systems, becomes less distinct and no-longer acts as a reliable separator of early- and late-type systems. The final set of stellar mass $-$ half-light radius relations are reported for a variety of galaxy population subsets in 10 bands ($ugrizZYJHK_s$) and are intended to provide a comprehensive low-z benchmark for the many ongoing high-z studies. Exploring the variation of the stellar mass $-$ half-light radius relations with wavelength we confirm earlier findings that galaxies appear more compact at longer wavelengths albeit at a smaller level than previously noted: at $10^{10}\mathcal{M_{\odot}}$ both spiral systems and ellipticals show a decrease in size of 13% from $g$ to $K_s$ (which is near linear in log wavelength). Finally we note that the sizes used in this work are derived from 2D S\'ersic light profile fitting (using GALFIT3), i.e., elliptical semi-major half light radii, improving on earlier low-z benchmarks based on circular apertures.

We present bivariate luminosity and stellar mass functions of H$\alpha$ star forming galaxies drawn from the Galaxy And Mass Assembly (GAMA) survey. While optically deep spectroscopic observations of GAMA over a wide sky area enable the detection of a large number of $0.001<{SFR}_{H\alpha}$ (M$_{\odot}$ yr$^{-1}$)$<100$ galaxies, the requirement for an H$\alpha$ detection in targets selected from an $r$-band magnitude limited survey leads to an incompleteness due to missing optically faint star forming galaxies. Using $z<0.1$ bivariate distributions as a reference we model the higher-$z$ distributions, thereby approximating a correction for the missing optically faint star forming galaxies to the local SFR and stellar mass densities. Furthermore, we obtain the $r$-band LFs and stellar mass functions of H$\alpha$ star forming galaxies from the bivariate LFs. As our sample is selected on the basis of detected H$\alpha$ emission, a direct tracer of on-going star formation, this sample represents a true star forming galaxy sample, and is drawn from both photometrically classified blue and red sub-populations, though mostly from the blue population. On average 20-30% of red galaxies at all stellar masses are star forming, implying that these galaxies may be dusty star forming systems.

We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G$^3$Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, $M_h$, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass with group luminosity, virial volume and central galaxy stellar mass, $M_\star$, agree very well with predictions from mock group catalogues constructed from a GALFORM semi-analytical galaxy formation model implemented in the Millennium $\Lambda$CDM N-body simulation; 2) the measured relations of halo mass with velocity dispersion and projected half-abundance radius show weak tension with mock predictions, hinting at problems in the mock galaxy dynamics and their small scale distribution; 3) the median $M_h|M_\star$ measured from weak lensing depends more sensitively on the lognormal dispersion in $M_\star$ at fixed $M_h$ than it does on the median $M_\star|M_h$. Our measurements suggest an intrinsic dispersion of $\sigma_{\log(M_\star)}\sim 0.15$; 4) Comparing our mass estimates with those in the catalogue, we find that the G$^3$Cv5 mass can give biased results when used to select subsets of the group sample. Of the various new halo mass estimators that we calibrate using our weak lensing measurements, group luminosity is the best single-proxy estimator of group mass.

In this paper we aim to (i) test the number of statistically distinct classes required to classify the local galaxy population, and, (ii) identify the differences in the physical and star formation properties of visually-distinct galaxies. To accomplish this, we analyse the structural parameters (effective radius r_e, effective surface brightness within r_e (mu_e), central surface brightness (mu_0), and S'ersic index (n)), obtained by fitting the light profile of 432 galaxies (0.002<z<=0.02 Viking Z-band), and their spectral energy distribution using multi-band photometry in 18 broadbands to obtain the stellar mass (M*), the star formation rate (SFR), the specific SFR (sSFR) and the dust mass (M_{dust}), respectively. We show that visually distinct, star-forming dwarf galaxies (irregulars, blue spheroids and low surface brightness galaxies) form a unimodal population in a parameter space mapped by mu_e, mu_0, n, r_e, SFR, sSFR, M*, M_{dust} and (g-i). The SFR and sSFR distribution of passively evolving (dwarf) ellipticals on the other hand, statistically distinguish them from other galaxies with similar luminosity, while the giant galaxies clearly segregate into star-forming spirals and passive lenticulars. We therefore suggest that the morphology classification scheme(s) used in literature for dwarf galaxies only reflect the observational differences based on luminosity and surface brightness among the apparent distinct classes, rather than any physical differences between them.

We use 80922 galaxies in the Galaxy And Mass Assembly (GAMA) survey to measure the galaxy luminosity function (LF) in different environments over the redshift range 0.04<z<0.26. The depth and size of GAMA allows us to define samples split by colour and redshift to measure the dependence of the LF on environment, redshift and colour. We find that the LF varies smoothly with overdensity, consistent with previous results, with little environmental dependent evolution over the last 3 Gyrs. The modified GALFORM model predictions agree remarkably well with our LFs split by environment, particularly in the most overdense environments. The LFs predicted by the model for both blue and red galaxies are consistent with GAMA for the environments and luminosities at which such galaxies dominate. Discrepancies between the model and the data seen in the faint end of the LF suggest too many faint red galaxies are predicted, which is likely to be due to the over-quenching of satellite galaxies. The excess of bright blue galaxies predicted in underdense regions could be due to the implementation of AGN feedback not being sufficiently effective in the lower mass halos.

We use the Galaxy And Mass Assembly (GAMA) survey to measure the local Universe mass dependent merger fraction and merger rate using galaxy pairs and the CAS structural method, which identifies highly asymmetric merger candidate galaxies. Our goals are to determine which types of mergers produce highly asymmetrical galaxies, and to provide a new measurement of the local galaxy major merger rate. We examine galaxy pairs at stellar mass limits down to $M_{*} = 10^{8}M_{\odot}$ with mass ratios of $<$100:1 and line of sight velocity differences of $\Delta V<500$ km s$^{-1}$. We find a significant increase in mean asymmetries for projected separations less than the sum of the individual galaxy's Petrosian 90 radii. For systems in major merger pairs with mass ratios of $<$4:1 both galaxies in the pair show a strong increase in asymmetry, while in minor merger systems (with mass ratios of $>$4:1) the lower mass companion becomes highly asymmetric, while the larger galaxy is much less affected. The fraction of highly asymmetric paired galaxies which have a major merger companion is highest for the most massive galaxies and drops progressively with decreasing mass. We calculate that the mass dependent major merger fraction is fairly constant at $\sim1.3-2\%$ between $10^{9.5}<M_{*}<10^{11.5} M_{\odot}$, and increases to $\sim4\%$ at lower masses. When the observability time scales are taken into consideration, the major merger rate is found to approximately triple over the mass range we consider. The total co-moving volume major merger rate over the range $10^{8.0}<M_{*}<10^{11.5} M_{\odot}$ is $(1.2 \pm 0.5) \times 10^{-3}$ $h^{3}_{70}$ Mpc$^{-3}$ Gyr$^{-1}$.

We measure the mass functions for generically red and blue galaxies, using a z < 0.12 sample of log M* > 8.7 field galaxies from the Galaxy And Mass Assembly (GAMA) survey. Our motivation is that, as we show, the dominant uncertainty in existing measurements stems from how 'red' and 'blue' galaxies have been selected/defined. Accordingly, we model our data as two naturally overlapping populations, each with their own mass function and colour-mass relation, which enables us characterise the two populations without having to specify a priori which galaxies are 'red' and 'blue'. Our results then provide the means to derive objective operational definitions for the terms 'red' and 'blue', which are based on the phenomenology of the colour-mass diagrams. Informed by this descriptive modelling, we show that: 1.) after accounting for dust, the stellar colours of 'blue' galaxies do not depend strongly on mass; 2.) the tight, flat 'dead sequence' does not extend much below log M* ~ 10.5; instead, 3.) the stellar colours of 'red' galaxies vary rather strongly with mass, such that lower mass 'red' galaxies have bluer stellar populations; 4.) below log M* ~ 9.3, the 'red' population dissolves into obscurity, and it becomes problematic to talk about two distinct populations; as a consequence, 5.) it is hard to meaningfully constrain the shape, including the possibility of an upturn, of the 'red' galaxy mass function below log M* ~ 9. Points 1-4 provide meaningful targets for models of galaxy formation and evolution to aim for.

We use a highly complete subset of the GAMA-II redshift sample to fully describe the stellar mass dependence of close-pairs and mergers between 10^8 Msun and 10^12 Msun. Using the analytic form of this fit we investigate the total stellar mass accreting onto more massive galaxies across all mass ratios. Depending on how conservatively we select our robust merging systems, the fraction of mass merging onto more massive companions is 2.0%-5.6%. Using the GAMA-II data we see no significant evidence for a change in the close-pair fraction between redshift $z = 0.05-0.2$. However, we find a systematically higher fraction of galaxies in similar mass close-pairs compared to published results over a similar redshift baseline. Using a compendium of data and the function $\gamma_M =A(1+z)m$ to predict the major close-pair fraction, we find fitting parameters of $A = 0.021 \pm 0.001$ and $m = 1.53 \pm 0.08$, which represents a higher low-redshift normalisation and shallower power-law slope than recent literature values. We find that the relative importance of in-situ star-formation versus galaxy merging is inversely correlated, with star-formation dominating the addition of stellar material below Mstar and merger accretion events dominating beyond Mstar. We find mergers have a measurable impact on the whole extent of the GSMF, manifest as a deepening of the dip in the GSMF over the next Gyr and an increase in Mstar by as much as 0.01-0.05 dex.

We present an estimate of the galaxy stellar mass function and its division by morphological type in the local (0.025 < z < 0.06) Universe. Adopting robust morphological classifications as previously presented (Kelvin et al.) for a sample of 3,727 galaxies taken from the Galaxy And Mass Assembly survey, we define a local volume and stellar mass limited sub-sample of 2,711 galaxies to a lower stellar mass limit of M = 10^9.0 M_sun. We confirm that the galaxy stellar mass function is well described by a double Schechter function given by M* = 10^10.64 M_sun, {\alpha}1 = -0.43, {\phi}*1 = 4.18 dex^-1 Mpc^-3, {\alpha}2 = -1.50 and {\phi}*2 = 0.74 dex^-1 Mpc^-3. The constituent morphological-type stellar mass functions are well sampled above our lower stellar mass limit, excepting the faint little blue spheroid population of galaxies. We find approximately 71+3-4% of the stellar mass in the local Universe is found within spheroid dominated galaxies; ellipticals and S0-Sas. The remaining 29+4-3% falls predominantly within late type disk dominated systems, Sab-Scds and Sd-Irrs. Adopting reasonable bulge-to-total ratios implies that approximately half the stellar mass today resides in spheroidal structures, and half in disk structures. Within this local sample, we find approximate stellar mass proportions for E : S0-Sa : Sab-Scd : Sd-Irr of 34 : 37 : 24 : 5.

We derive the close pair fractions and volume merger rates as a function of luminosity and morphology for galaxies in the GAMA survey with -23 < M(r) < -17 at 0.01 < z < 0.22. The merger fraction is about 0.015 at all luminosities (assuming 1/2 of pairs merge) and the volume merger rate is about 0.00035 per cubic Mpc per Gyr. Dry mergers (between red or spheroidal galaxies) are uncommon and decrease with decreasing luminosity. Fainter mergers are wet, between blue or disky galaxies. Damp mergers (one of each type) follow the average of dry and wet mergers. In the brighter luminosity bin (-23 < M(r) < -20) the merger rate evolution is flat, irrespective of colour or morphology. The makeup of the merging population does not change since z = 0.2. Major mergers and dry mergers appear comparatively unimportant in the buildup of the red sequence over the past 2 Gyr. We compare the colour, morphology, environmental density and degree of activity of galaxies in pairs to those of more isolated objects in the same volume. Galaxies in close pairs tend to be both redder and slightly more spheroid-dominated. This may be due to "harassment" in multiple previous passes prior to the current interaction. Galaxy pairs do not appear to prefer significantly denser environments. There is no evidence of an enhancement in the AGN fraction in pairs, compared to other galaxies in the same volume.

We investigate the relationship between colour and structure within galaxies using a large, volume-limited sample of bright, low-redshift galaxies with optical to near-infrared imaging from the GAMA survey. We fit single-component, wavelength-dependent, elliptical Sersic models to all passbands simultaneously, using software developed by the MegaMorph project. Dividing our sample by Sersic index and colour, the recovered wavelength variations in effective radius (R_e) and Sersic index (n) reveal the internal structure, and hence formation history, of different types of galaxies. All these trends depend on n; some have an additional dependence on galaxy colour. Late-type galaxies (n_r < 2.5) show a dramatic increase in Sersic index with wavelength. This might be a result of their two-component (bulge-disk) nature, though stellar population gradients within each component and dust attenuation are likely to play a role. All galaxies show a substantial decrease in R_e with wavelength. This is strongest for early-types (n_r > 2.5), even though they maintain constant n with wavelength, revealing that ellipticals are a superimposition of different stellar populations associated with multiple collapse and merging events. Processes leading to structures with larger R_e must be associated with lower metallicity or younger stellar populations. This appears to rule out the formation of young cores through dissipative gas accretion as an important mechanism in the recent lives of luminous elliptical galaxies.

We report a highly significant ($>10\sigma$) spatial correlation between galaxies with $S_{350\mu\rm m}\ge 30\,$mJy detected in the equatorial fields of the \textsl{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) with estimated redshifts $\gtrsim 1.5$, and SDSS or GAMA galaxies at $0.2\le z\le 0.6$. The significance of the cross-correlation is much higher than those reported so far for samples with non-overlapping redshift distributions selected in other wavebands. Extensive, realistic simulations of clustered sub-mm galaxies amplified by foreground structures confirm that the cross-correlation is explained by weak gravitational lensing ($\mu<2$). The simulations also show that the measured amplitude and range of angular scales of the signal are larger than can be accounted for by galaxy-galaxy weak lensing. However, for scales $\lesssim 2\,$arcmin, the signal can be reproduced if SDSS/GAMA galaxies act as signposts of galaxy groups/clusters with halo masses in the range $10^{13.2}$--$10^{14.5} M_{\odot}$. The signal detected on larger scales appears to reflect the clustering of such halos.

The Galaxy And Mass Assembly (GAMA) survey has obtained spectra of over 230000 targets using the Anglo-Australian Telescope. To homogenise the redshift measurements and improve the reliability, a fully automatic redshift code was developed (autoz). The measurements were made using a cross-correlation method for both absorption-line and emission-line spectra. Large deviations in the high-pass filtered spectra are partially clipped in order to be robust against uncorrected artefacts and to reduce the weight given to single-line matches. A single figure of merit (FOM) was developed that puts all template matches onto a similar confidence scale. The redshift confidence as a function of the FOM was fitted with a tanh function using a maximum likelihood method applied to repeat observations of targets. The method could be adapted to provide robust automatic redshifts for other large galaxy redshift surveys. For the GAMA survey, there was a substantial improvement in the reliability of assigned redshifts and in the lowering of redshift uncertainties with a median velocity uncertainty of 33 km/s.

We present a large-scale clustering analysis of radio galaxies in the Very Large Array (VLA) Faint Images of the Radio Sky at Twenty-cm (FIRST) survey over the Galaxy And Mass Assembly (GAMA) survey area, limited to S1.4 GHz >1 mJy with spectroscopic and photometric redshift limits up to r < 19.8 and r < 22 mag, respectively. For the GAMA spectroscopic matches, we present the redshift-space and projected correlation functions, the latter of which yielding a correlation length r0 ~ 8.2 Mpc/h and linear bias of ~1.9 at z ~ 0.34. Furthermore, we use the angular two-point correlation function w({\theta}) to determine spatial clustering properties at higher redshifts. We find r0 to increase from ~6 to ~14 Mpc/h between z = 0.3 and z = 1.55, with the corresponding bias increasing from ~2 to ~10 over the same range. Our results are consistent with the bias prescription implemented in the SKADS simulations at low redshift, but exceed these predictions at z > 1. This is indicative of an increasing (rather than fixed) halo mass and/or AGN fraction at higher redshifts or a larger typical halo mass for the more abundant FRI sources.

We have analysed the growth of Brightest Group Galaxies and Brightest Cluster Galaxies (BGGs/BCGs) over the last 3 billion years using a large sample of 883 galaxies from the Galaxy And Mass Assembly Survey. By comparing the stellar mass of BGGs and BCGs in groups and clusters of similar dynamical masses, we find no significant growth between redshift $z=0.27$ and $z=0.09$. We also examine the number of BGGs/BCGs that have line emission, finding that approximately 65 per cent of BGGs/BCGs show H$\alpha$ in emission. From the galaxies where the necessary spectroscopic lines were accurately recovered (54 per cent of the sample), we find that half of this (i.e. 27 per cent of the sample) harbour on-going star formation with rates up to $10\,$M$_{\odot}$yr$^{-1}$, and the other half (i.e. 27 per cent of the sample) have an active nucleus (AGN) at the centre. BGGs are more likely to have ongoing star formation, while BCGs show a higher fraction of AGN activity. By examining the position of the BGGs/BCGs with respect to their host dark matter halo we find that around 13 per cent of them do not lie at the centre of the dark matter halo. This could be an indicator of recent cluster-cluster mergers. We conclude that BGGs and BCGs acquired their stellar mass rapidly at higher redshifts as predicted by semi-analytic models, mildly slowing down at low redshifts.

Based on data from the Galaxy and Mass Assembly (GAMA) survey, we report on the discovery of structures that we refer to as `tendrils' of galaxies: coherent, thin chains of galaxies that are rooted in filaments and terminate in neighbouring filaments or voids. On average, tendrils contain 6 galaxies and span 10 $h^{-1}$ Mpc. We use the so-called line correlation function to prove that tendrils represent real structures rather than accidental alignments. We show that voids found in the SDSS-DR7 survey that overlap with GAMA regions contain a large number of galaxies, primarily belonging to tendrils. This implies that void sizes are strongly dependent on the number density and sensitivity limits of a survey. We caution that galaxies in low density regions, that may be defined as `void galaxies,' will have local galaxy number densities that depend on such observational limits and are likely higher than can be directly measured.

We report the morphological classification of 3727 galaxies from the Galaxy and Mass Assembly survey with M_r < -17.4 mag and in the redshift range 0.025 < z < 0.06 (2.1 x 10^5 Mpc^3 ) into E, S0-Sa, SB0-SBa, Sab-Scd, SBab-SBcd, Sd-Irr and little blue spheroid classes. Approximately 70% of galaxies in our sample are disk dominated systems, with the remaining ~30% spheroid dominated. We establish the robustness of our classifications, and use them to derive morphological-type luminosity functions and luminosity densities in the ugrizYJHK passbands, improving on prior studies that split by global colour or light profile shape alone. We find that the total galaxy luminosity function is best described by a double-Schechter function while the constituent morphological-type luminosity functions are well described by a single-Schechter function. These data are also used to derive the star-formation rate densities for each Hubble class, and the attenuated and unattenuated (corrected for dust) cosmic spectral energy distributions, i.e., the instantaneous energy production budget. While the observed optical/near-IR energy budget is dominated 58:42 by galaxies with a significant spheroidal component, the actual energy production rate is reversed, i.e., the combined disk dominated populations generate ~1.3x as much energy as the spheroid dominated populations. On the grandest scale, this implies that chemical evolution in the local Universe is currently confined to mid-type spiral classes like our Milky Way.

From a volume limited sample of 45,542 galaxies and 6,000 groups with $z \leq 0.213$ we use an adapted minimal spanning tree algorithm to identify and classify large scale structures within the Galaxy and Mass Assembly (GAMA) survey. Using galaxy groups, we identify 643 filaments across the three equatorial GAMA fields that span up to 200 $h^{-1}$ Mpc in length, each with an average of 8 groups within them. By analysing galaxies not belonging to groups we identify a secondary population of smaller coherent structures composed entirely of galaxies, dubbed `tendrils' that appear to link filaments together, or penetrate into voids, generally measuring around 10 $h^{-1}$ Mpc in length and containing on average 6 galaxies. Finally we are also able to identify a population of isolated void galaxies. By running this algorithm on GAMA mock galaxy catalogues we compare the characteristics of large scale structure between observed and mock data; finding that mock filaments reproduce observed ones extremely well. This provides a probe of higher order distribution statistics not captured by the popularly used two-point correlation function.

We present the first application of a "multiple-tracer" redshift-space distortion (RSD) analysis to an observational galaxy sample, using data from the Galaxy and Mass Assembly survey (GAMA). Our dataset is an r < 19.8 magnitude-limited sample of 178,579 galaxies covering redshift interval z < 0.5 and area 180 deg^2. We obtain improvements of 10-20% in measurements of the gravitational growth rate compared to a single-tracer analysis, deriving from the correlated sample variance imprinted in the distributions of the overlapping galaxy populations. We present new expressions for the covariances between the auto-power and cross-power spectra of galaxy samples that are valid for a general survey selection function and weighting scheme. We find no evidence for a systematic dependence of the measured growth rate on the galaxy tracer used, justifying the RSD modelling assumptions, and validate our results using mock catalogues from N-body simulations. For multiple tracers selected by galaxy colour, we measure normalized growth rates in two independent redshift bins f*sigma_8(z=0.18) = 0.36 +/- 0.09 and f*sigma_8(z=0.38) = 0.44 +/- 0.06, in agreement with standard GR gravity and other galaxy surveys at similar redshifts.

We present observations of 18 galaxies from the Galaxy And Mass Assembly (GAMA) survey made with the SPIRAL optical integral field unit (IFU) on the Anglo-Australian Telescope. The galaxies are selected to have a narrow range in stellar mass (6x10^9Msolar < M* <2x10^10 Msolar) in order to focus on the effects of environment. Local galaxy environments are measured quantitatively using 5th nearest neighbour surface densities. We find that the total star formation rates (SFR) measured from the IFU data are consistent with total SFRs measured from aperture correcting either GAMA or Sloan Digital Sky Survey single-fibre observations. The mean differences are SFR_GAMA/SFR_IFU = 1.26+/-0.23, sigma=0.90 and for the Sloan Digital Sky Survey we similarly find SFR_Brinchmann/SFR_IFU = 1.34+/-0.17, sigma=0.67. Examining the relationships with environment, we find off-centre and clumpy Halpha emission is not significantly dependent on environment, being present in 2/7 (29^+20_-11 per cent) galaxies in high-density environments (>0.77 Mpc^-2), and 5/11 (45^+15_-13 per cent) galaxies in low-density environments (<0.77 Mpc^-2). We find a weak but not significant relationship of the total star formation rates of star-forming galaxies with environment. Due to the size of our sample and the scatter observed we do not draw a definitive conclusion about a possible SFR dependence on environment. Examining the spatial distribution of the Halpha emission, we find no evidence for a change in shape or amplitude of the radial profile of star-forming galaxies with environment. If these observations are borne out in larger samples this would infer that any environment-driven star-formation suppression must either act very rapidly (the `infall-and-quench' model) or that galaxies must evolve in a density-dependent manner (an `in-situ evolution' model).

We present a combined optical and X-ray analysis of the rich cluster Abell 1882 with the aim of identifying merging substructure and understanding the recent assembly history of this system. Our optical data consist of spectra drawn from the Galaxy and Mass Assembly (GAMA) survey, which lends itself to this kind of detailed study thanks to its depth and high spectroscopic completeness. We use 283 spectroscopically confirmed cluster members to detect and characterize substructure. We complement the optical data with X-ray data taken with both Chandra and XMM. Our analysis reveals that A1882 harbors two main components, A1882A and A1882B, which have a projected separation of 2Mpc and a line of sight velocity difference of v_{los}=-428km/s. The primary system, A1882A, has velocity dispersion sigma_v=500km/s and Chandra (XMM) temperature kT=3.57keV (3.31keV) while the secondary, A1882B, has sigma_v=457km/s and Chandra (XMM) temperature kT=2.39keV (2.12keV). The optical and X-ray estimates for the masses of the two systems are consistent within the uncertainties and indicate that there is twice as much mass in A1882A (M_{500}=1.5-1.9x10^{14}Msolar) when compared with A1882B (M_{500}=0.8-1.0x10^{14}Msolar). We interpret the A1882A/A1882B system as being observed prior to a core passage. Supporting this interpretation is the large projected separation of A1882A and A1882B and the dearth of evidence for a recent (<2Gyr) major interaction in the X-ray data. Two-body analyses indicate that A1882A and A1882B form a bound system with bound incoming solutions strongly favored. We compute blue fractions of f_b=0.28 and 0.18 for the spectroscopically confirmed member galaxies within r_{500} of the centers of A1882A and A1882B, respectively. These blue fractions do not differ significantly from the blue fraction measured from an ensemble of 20 clusters with similar mass and redshift.

We constrain the distributions of projected radial alignment angles of satellite galaxy shapes within the Galaxy And Mass Assembly survey group catalogue. We identify the galaxy groups using spectroscopic redshifts and measure galaxy projected ellipticities from Sloan Digital Sky Survey imaging. With a sample of 3,850 groups with 13,655 satellite galaxies with high quality shape measurements, we find a less than 2-sigma signal of radial alignments in the mean projected ellipticity components and the projected position angle when using galaxy shape estimates optimized for weak lensing measurements. Our radial alignment measurement increases to greater than 3-sigma significance relative to the expectation for no alignments if we use 2-D S\'ersic model fits to define galaxy orientations. Our weak measurement of radial alignments is in conflict with predictions from dark matter N-body simulations, which we interpret as evidence for large mis-alignments of baryons and dark matter in group and cluster satellites. Within our uncertainties, that are dominated by our small sample size, we find only weak and marginally significant trends of the radial alignment angle distributions on projected distance from the group centre, host halo mass, and redshift that could be consistent with a tidal torquing mechanism for radial alignments. Using our lensing optimized shape estimators, we estimate that intrinsic alignments of galaxy group members may contribute a systematic error to the mean differential projected surface mass density of groups inferred from weak lensing observations by -1 +/- 20% at scales around 300 kpc/h from the group centre assuming a photometric redshift r.m.s. error of 10%, and given our group sample with median redshift of 0.17 and median virial masses ~10^{13} h^{-1}M_{sun}.

We present evidence for stochastic star formation histories in low-mass (M* < 10^10 Msun) galaxies from observations within the Galaxy And Mass Assembly (GAMA) survey. For ~73,000 galaxies between 0.05<z<0.32, we calculate star formation rates (SFR) and specific star formation rates (SSFR = SFR/M*) from spectroscopic Halpha measurements and apply dust corrections derived from Balmer decrements. We find a dependence of SSFR on stellar mass, such that SSFRs decrease with increasing stellar mass for star-forming galaxies, and for the full sample, SSFRs decrease as a stronger function of stellar mass. We use simple parametrizations of exponentially declining star formation histories to investigate the dependence on stellar mass of the star formation timescale and the formation redshift. We find that parametrizations previously fit to samples of z~1 galaxies cannot recover the distributions of SSFRs and stellar masses observed in the GAMA sample between 0.05<z<0.32. In particular, a large number of low-mass (M* < 10^10 Msun) galaxies are observed to have much higher SSFRs than can be explained by these simple models over the redshift range of 0.05<z<0.32, even when invoking mass-dependent staged evolution. For such a large number of galaxies to maintain low stellar masses, yet harbour such high SSFRs, requires the late onset of a weak underlying exponentially declining SFH with stochastic bursts of star formation superimposed.

A full appreciation of the role played by gas metallicity (Z), star-formation rate (SFR), and stellar mass is fundamental to understanding how galaxies form and evolve. The connections between these three parameters at different redshifts significantly affect galaxy evolution, and thus provide important constraints for galaxy evolution models. Using data from the Sloan Digital Sky Survey-Data Release 7 (SDSS-DR7) and the Galaxy and Mass Assembly (GAMA) surveys we study the relationships and dependencies between SFR, Z, and stellar mass, as well as the Fundamental Plane for star-forming galaxies. We combine both surveys using volume-limited samples up to a redshift of z ~ 0.36. The GAMA and SDSS surveys complement each other when analyzing the relationships between SFR, Mass and Z. We present evidence for SFR and metallicity evolution to z ~ 0.2. We study the dependencies between SFR, Mass, Z, and specific star-formation rate (SSFR) on the M-Z, M-SFR, M-SSFR, Z-SFR, and Z-SSFR relations, finding strong correlations between all. Based on those dependencies, we propose a simple model that allows us to explain the different behaviour observed between low and high mass galaxies. Finally, our analysis allows us to confirm the existence of a Fundamental Plane, for which Mass=f(Z, SFR) in star-forming galaxies.

Measurements of the low-z Halpha luminosity function have a large dispersion in the local number density of sources, and correspondingly in the SFR density. The possible causes for these discrepancies include limited volume sampling, biases arising from survey sample selection, different methods of correcting for dust obscuration and AGN contamination. The Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey (SDSS) provide deep spectroscopic observations over a wide sky area enabling detection of a large sample of star-forming galaxies spanning 0.001<SFR(Halpha)<100 with which to robustly measure the evolution of the SFR density in the low-z universe. The large number of high SFR galaxies present in our sample allow an improved measurement of the bright end of the luminosity function, indicating that the decrease in number density of sources at bright luminosities is best described by a Saunders functional form rather than the traditional Schechter function. This result is consistent with other published luminosity functions in the FIR and radio. For GAMA and SDSS we find the r-band apparent magnitude limit, combined with the subsequent requirement for Halpha detection leads to an incompleteness due to missing bright Halpha sources with faint r-band magnitudes.

We present the design and performance of the multi-object fiber spectrographs for the Sloan Digital Sky Survey (SDSS) and their upgrade for the Baryon Oscillation Spectroscopic Survey (BOSS). Originally commissioned in Fall 1999 on the 2.5-m aperture Sloan Telescope at Apache Point Observatory, the spectrographs produced more than 1.5 million spectra for the SDSS and SDSS-II surveys, enabling a wide variety of Galactic and extra-galactic science including the first observation of baryon acoustic oscillations in 2005. The spectrographs were upgraded in 2009 and are currently in use for BOSS, the flagship survey of the third-generation SDSS-III project. BOSS will measure redshifts of 1.35 million massive galaxies to redshift 0.7 and Lyman-alpha absorption of 160,000 high redshift quasars over 10,000 square degrees of sky, making percent level measurements of the absolute cosmic distance scale of the Universe and placing tight constraints on the equation of state of dark energy. The twin multi-object fiber spectrographs utilize a simple optical layout with reflective collimators, gratings, all-refractive cameras, and state-of-the-art CCD detectors to produce hundreds of spectra simultaneously in two channels over a bandpass covering the near ultraviolet to the near infrared, with a resolving power R = \lambda/FWHM ~ 2000. Building on proven heritage, the spectrographs were upgraded for BOSS with volume-phase holographic gratings and modern CCD detectors, improving the peak throughput by nearly a factor of two, extending the bandpass to cover 360 < \lambda < 1000 nm, and increasing the number of fibers from 640 to 1000 per exposure. In this paper we describe the original SDSS spectrograph design and the upgrades implemented for BOSS, and document the predicted and measured performances.

We compare the environmental and star formation properties of far-infrared detected and non--far-infrared detected galaxies out to $z \sim0.5$. Using optical spectroscopy and photometry from the Galaxy And Mass Assembly (GAMA) and Sloan Digital Sky Survey (SDSS), with far-infrared observations from the {\em Herschel}-ATLAS Science Demonstration Phase (SDP), we apply the technique of Voronoi Tessellations to analyse the environmental densities of individual galaxies. Applying statistical analyses to colour, $r-$band magnitude and redshift-matched samples, we show there to be a significant difference at the 3.5$\sigma$ level between the normalized environmental densities of these two populations. This is such that infrared emission (a tracer of star formation activity) favours underdense regions compared to those inhabited by exclusively optically observed galaxies selected to be of the same $r-$band magnitude, colour and redshift. Thus more highly star-forming galaxies are found to reside in the most underdense environments, confirming previous studies that have proposed such a correlation. However, the degeneracy between redshift and far-infrared luminosity in our flux-density limited sample means that we are unable to make a stronger statement in this respect. We then apply our method to synthetic light cones generated from semi-analytic models, finding that over the whole redshift distribution the same correlations between star-formation rate and environmental density are found.

We study the interplay between gas phase metallicity (Z), specific star formation rate (SSFR) and neutral hydrogen gas (HI) for galaxies of different stellar masses. Our study uses spectroscopic data from GAMA and SDSS star forming galaxies, as well as HI-detection from the ALFALFA and GASS public catalogues. We present a model based on the Z-SSFR relation that shows that at a given stellar mass, depending on the amount of gas, galaxies will follow opposite behaviours. Low-mass galaxies with a large amount of gas will show high SSFR and low metallicities, while low-mass galaxies with small amounts of gas will show lower SSFR and high metallicities. In contrast, massive galaxies with a large amount of gas will show moderate SSFR and high metallicities, while massive galaxies with small amounts of gas will show low SSFR and low metallicities. Using ALFALFA and GASS counterparts, we find that the amount of gas is related to those drastic differences in Z and SSFR for galaxies of a similar stellar mass.

We present two large, nearby (0.013$\le$z$\le$0.06) samples of Early-Type Galaxies (ETGs): a visually classified sample of 220 ETGs, created using source-matched data from the Galaxy and Mass Assembly (GAMA) database with FIR/sub-mm detections from $Herschel$-ATLAS; and a visually classified sample of 551 ETGs which are undetected with $Herschel$-ATLAS. Active galactic nuclei (AGN) are removed from our samples using optical emission line diagnostics. These samples are scrutinised to determine characteristics of sub-mm detected versus undetected ETGs. We find similarities in the stellar mass distributions of the two ETG samples but testing other properties uncovers significant differences. The sub-mm detected sample is shown to have lower concentration and S\'ersic indices than those in the undetected sample - a result which may be linked to the presence of dust in the former. Optical and UV-optical colours are also shown to be much bluer, indicating that the dust is linked with recent star formation. The intrinsic effective radii are on average 1.5 times larger for the sub-mm detected ETGs. Surface densities and groups data from the GAMA database are examined for the two samples, leading to the conclusion that dusty ETGs inhabit sparser environments than non-dusty ETGs in the nearby universe, although environments of the brightest ETGs are shown to differ the least. Modified Planck functions are fit to the H-ATLAS detected PACS and SPIRE fluxes for ETGs with sub-mm flux densities of at least 3$\sigma$ in the 350$\mu$m SPIRE band, giving a resultant mean cold dust temperature of T$_{d}$=22.1K, with a range of 9-30K. The corresponding mean dust mass is 1.8$\times10^{7}$M$_{\odot}$, with a range of (0.08-35.0)$\times10^{7}$M$_{\odot}$. The dust masses calculated from these fits, normalised by stellar mass, are shown to increase with decreasing stellar mass and bluer colours.

We report the discovery of a well-defined correlation between B-band face-on central optical depth due to dust, \tau^f_B, and the stellar mass surface density, \mu_{*}, of nearby (z < 0.13) spiral galaxies: log(\tau^f_B) = 1.12(+-0.11)log(\mu_{*}/M_sol kpc^2)-8.6(+-0.8). This relation was derived from a sample of spiral galaxies taken from the Galaxy and Mass Assembly (GAMA) survey and detected in the FIR/submm in the Herschel-ATLAS survey. Using a quantitative analysis of the NUV attenuation-inclination relation for complete samples of GAMA spirals categorized according to \mu_{*} we demonstrate that this correlation can be used to statistically correct for dust attenuation purely on the basis of optical photometry and S'ersic-profile morphological fits. Considered together with previously established empirical relationships between stellar mass, metallicity and gas mass, the near linearity and high constant of proportionality of the \tau^f_B-\mu_{*} relation disfavors a stellar origin for the bulk of refractory grains in spiral galaxies, instead being consistent with the existence of a ubiquitous and very rapid mechanism for the growth of dust in the ISM. We use the \tau^f_B-\mu_{*} relation in conjunction with the radiation transfer model for spiral galaxies of Popescu & Tuffs (2011) to derive intrinsic scaling relations between specific star formation rate (sSFR), stellar mass, and \mu_{*}, in which the attenuation of the UV light used to measure the SFR is corrected on an object-to-object basis. A marked reduction in scatter in these relations is achieved which is demonstrably due to correction of both the inclination-dependent and face-on components of attenuation. Our results are consistent with a picture of spiral galaxies in which most of the submm emission originates from grains residing in translucent structures, exposed to UV in the diffuse interstellar radiation field.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

In this work we investigate in detail the effects local environment (groups and pairs) has on galaxies with stellar mass similar to the Milky-Way (L* galaxies). A volume limited sample of 6,150 galaxies is classified to determine emission features, morphological type and presence of a disk. This sample allows for characteristics of galaxies to be isolated (e.g. stellar mass and group halo mass), and their codependencies determined. We observe that galaxy-galaxy interactions play the most important role in shaping the evolution within a group halo, the main role of halo mass is in gathering the galaxies together to encourage such interactions. Dominant pair galaxies find their overall star formation enhanced when the pair's mass ratio is close to 1, otherwise we observe the same galaxies as we would in an unpaired system. The minor galaxy in a pair is greatly affected by its companion galaxy, and whilst the star forming fraction is always suppressed relative to equivalent stellar mass unpaired galaxies, it becomes lower still when the mass ratio of a pair system increases. We find that, in general, the close galaxy-galaxy interaction rate drops as a function of halo mass for a given amount of stellar mass. We find evidence of a local peak of interactions for Milky-Way stellar mass galaxies in Milky-Way halo mass groups. Low mass halos, and in particular Local Group mass halos, are an important environment for understanding the typical evolutionary path of a unit of stellar mass. We find compelling evidence for galaxy conformity in both groups and pairs, where morphological type conformity is dominant in groups, and emission class conformity is dominant in pairs. This suggests that group scale conformity is the result of many galaxy encounters over an extended period of time, whilst pair conformity is a fairly instantaneous response to a transitory interaction.

We have measured the clustering properties of low-redshift (z < 0.3) sub-mm galaxies detected at 250 micron in the Herschel-ATLAS Science Demonstration Phase (SDP) field. We selected a sample for which we have high-quality spectroscopic redshifts, obtained from reliably matching the 250-micron sources to a complete (for r < 19.4) sample of galaxies from the GAMA database. Both the angular and spatial clustering strength are measured for all z < 0.3 sources as well as for five redshift slices with thickness delta z=0.05 in the range 0.05 < z < 0.3. Our measured spatial clustering length r_0 is comparable to that of optically-selected, moderately star-forming (blue) galaxies: we find values around 5 Mpc. One of the redshift bins contains an interesting structure, at z = 0.164.

Context: The mass-metallicity relationship (MMR) of star-forming galaxies is well-established, however there is still some disagreement with respect to its exact shape and its possible dependence on other observables. Aims: We measure the MMR in the Galaxy And Mass Assembly (GAMA) survey. We compare our measured MMR to that measured in the Sloan Digital Sky Survey (SDSS) and study the dependence of the MMR on various selection criteria to identify potential causes for disparities seen in the literature. Methods: We use strong emission line ratio diagnostics to derive oxygen abundances. We then apply a range of selection criteria for the minimum signal-to-noise in various emission lines, as well as the apparent and absolute magnitude to study variations in the inferred MMR. Results: The shape and position of the MMR can differ significantly depending on the metallicity calibration and selection used. After selecting a robust metallicity calibration amongst those tested, we find that the mass-metallicity relation for redshifts 0.061< z<0.35 in GAMA is in reasonable agreement with that found in the SDSS despite the difference in the luminosity range probed. Conclusions: In view of the significant variations of the MMR brought about by reasonable changes in the sample selection criteria and method, we recommend that care be taken when comparing the MMR from different surveys and studies directly. We also conclude that there could be a modest level of evolution over 0.06<z<0.35 within the GAMA sample.

We use the GAMA I dataset combined with GALEX, SDSS and UKIDSS imaging to construct the low-redshift (z<0.1) galaxy luminosity functions in FUV, NUV, ugriz, and YJHK bands from within a single well constrained volume of 3.4 x 10^5 (Mpc/h)^{3}. The derived luminosity distributions are normalised to the SDSS DR7 main survey to reduce the estimated cosmic variance to the 5 per cent level. The data are used to construct the cosmic spectral energy distribution (CSED) from 0.1 to 2.1 \mum free from any wavelength dependent cosmic variance for both the elliptical and non-elliptical populations. The two populations exhibit dramatically different CSEDs as expected for a predominantly old and young population respectively. Using the Driver et al. (2008) prescription for the azimuthally averaged photon escape fraction, the non-ellipticals are corrected for the impact of dust attenuation and the combined CSED constructed. The final results show that the Universe is currently generating (1.8 +/- 0.3) x 10^{35} h W Mpc^{-3} of which (1.2 +/- 0.1) x 10^{35} h W Mpc^{-3} is directly released into the inter-galactic medium and (0.6 +/- 0.1) x 10^{35} h W Mpc^{-3} is reprocessed and reradiated by dust in the far-IR. Using the GAMA data and our dust model we predict the mid and far-IR emission which agrees remarkably well with available data. We therefore provide a robust description of the pre- and post dust attenuated energy output of the nearby Universe from 0.1micron to 0.6mm. The largest uncertainty in this measurement lies in the mid and far-IR bands stemming from the dust attenuation correction and its currently poorly constrained dependence on environment, stellar mass, and morphology.

Analysing all Galaxy and Mass Assembly (GAMA) galaxies within a factor two (+/- 0.3 dex) of the stellar mass of the Milky Way (MW), there is a 11.9% chance that one of these galaxies will have a close companion (within a projected separation of 70 kpc and radial separation of 400 km/s) that is at least as massive as the Large Magellanic Cloud (LMC). Two close companions at least as massive as the Small Magellanic Cloud (SMC) are rare at the 3.4% level. Two full analogues to the MW- LMC-SMC system were found in GAMA (all galaxies late-type and star forming), suggesting such a combination of close together, late-type, star-forming galaxies is rare: only 0.4% of MW mass galaxies (in the range where we could observe both the LMC and SMC) have such a system. In summary, the MW-LMC-SMC system is a 2.7? event (when recast into Gaussian statistics). Using cross-correlation comparisons we find that there is a preference for SMC- LMC binary pair analogues to be located within 2 Mpc of a range of different lumi- nosity groups. There is a particular preference is for such binaries to be located near LG luminosity systems. When these groups are subdivided into small magnitude gap and large magnitude gap subsets, the binaries prefer to be spatially associated with the small magnitude gap systems. These systems will be dynamically less evolved, but still offer the same amount of gravitational dark matter. This suggests that binaries such as the SMC-LMC might be transient systems, usually destroyed during vigorous merger events. Details of a particularly striking analogue to the MW-SMC-LMC and M31 complex are included.

We measure the two-point angular correlation function of a sample of 4,289,223 galaxies with r < 19.4 mag from the Sloan Digital Sky Survey as a function of photometric redshift, absolute magnitude and colour down to M_r - 5log h = -14 mag. Photometric redshifts are estimated from ugriz model magnitudes and two Petrosian radii using the artificial neural network package ANNz, taking advantage of the Galaxy and Mass Assembly (GAMA) spectroscopic sample as our training set. The photometric redshifts are then used to determine absolute magnitudes and colours. For all our samples, we estimate the underlying redshift and absolute magnitude distributions using Monte-Carlo resampling. These redshift distributions are used in Limber's equation to obtain spatial correlation function parameters from power law fits to the angular correlation function. We confirm an increase in clustering strength for sub-L* red galaxies compared with ~L* red galaxies at small scales in all redshift bins, whereas for the blue population the correlation length is almost independent of luminosity for ~L* galaxies and fainter. A linear relation between relative bias and log luminosity is found to hold down to luminosities L~0.03L*. We find that the redshift dependence of the bias of the L* population can be described by the passive evolution model of Tegmark & Peebles (1998). A visual inspection of a random sample of our r < 19.4 sample of SDSS galaxies reveals that about 10 per cent are spurious, with a higher contamination rate towards very faint absolute magnitudes due to over-deblended nearby galaxies. We correct for this contamination in our clustering analysis.

We present a detailed investigation into the effects of galaxy environment on their star formation rates (SFR) using galaxies observed in the Galaxy and Mass Assembly Survey (GAMA). We use three independent volume-limited samples of galaxies within z < 0.2 and Mr < -17.8. We investigate the known SFR-density relationship and explore in detail the dependence of SFR on stellar mass and density. We show that the SFR-density trend is only visible when we include the passive galaxy population along with the star-forming population. This SFR-density relation is absent when we consider only the star-forming population of galaxies, consistent with previous work. While there is a strong dependence of the EWH?a on density we find, as in previous studies, that these trends are largely due to the passive galaxy population and this relationship is absent when considering a "star-forming" sample of galaxies. We find that stellar mass has the strongest influence on SFR and EWH?a with the environment having no significant effect on the star-formation properties of the star forming population. We also show that the SFR-density relationship is absent for both early and late-type star-forming galaxies. We conclude that the stellar mass has the largest impact on the current SFR of a galaxy, and any environmental effect is not detectable. The observation that the trends with density are due to the changing morphology fraction with density implies that the timescales must be very short for any quenching of the SFR in infalling galaxies. Alternatively galaxies may in fact undergo predominantly in-situ evolution where the infall and quenching of galaxies from the field into dense environments is not the dominant evolutionary mode.

We have generated complementary halo mass estimates for all groups in the Galaxy And Mass Assembly Galaxy Group Catalogue (GAMA G3Cv1) using a modified caustic mass estimation algorithm, originally developed by Diaferio & Geller (1997). We calibrate the algorithm by applying it on a series of 9 GAMA mock galaxy light cones and investigate the effects of using different definitions for group centre and size. We select the set of parameters that provide median-unbiased mass estimates when tested on mocks, and generate mass estimates for the real group catalogue. We find that on average, the caustic mass estimates agree with dynamical mass estimates within a factor of 2 in 90.8 +/- 6.1% groups and compares equally well to velocity dispersion based mass estimates for both high and low multiplicity groups over the full range of masses probed by the G3Cv1.

We use the Herschel-ATLAS survey to conduct the first large-scale statistical study of the submm properties of optically selected galaxies. Using ~80,000 r-band selected galaxies from 126 deg^2 of the GAMA survey, we stack into submm imaging at 250, 350 and 500{\mu}m to gain unprecedented statistics on the dust emission from galaxies at z < 0.35. We find that low redshift galaxies account for 5% of the cosmic 250{\mu}m background (4% at 350{\mu}m; 3% at 500{\mu}m), of which approximately 60% comes from 'blue' and 20% from 'red' galaxies (rest-frame g - r). We compare the dust properties of different galaxy populations by dividing the sample into bins of optical luminosity, stellar mass, colour and redshift. In blue galaxies we find that dust temperature and luminosity correlate strongly with stellar mass at a fixed redshift, but red galaxies do not follow these correlations and overall have lower luminosities and temperatures. We make reasonable assumptions to account for the contaminating flux from lensing by red sequence galaxies and conclude that galaxies with different optical colours have fundamentally different dust emission properties. Results indicate that while blue galaxies are more luminous than red galaxies due to higher temperatures, the dust masses of the two samples are relatively similar. Dust mass is shown to correlate with stellar mass, although the dust/stellar mass ratio is much higher for low stellar mass galaxies, consistent with the lowest mass galaxies having the highest specific star formation rates. We stack the 250{\mu}m/NUV luminosity ratio, finding results consistent with greater obscuration of star formation at lower stellar mass and higher redshift. Submm luminosities and dust masses of all galaxies are shown to evolve strongly with redshift, indicating a fall in the amount of obscured star formation in ordinary galaxies over the last four billion years.

We determine the low-redshift field galaxy stellar mass function (GSMF) using an area of 143 deg^2 from the first three years of the Galaxy And Mass Assembly (GAMA) survey. The magnitude limits of this redshift survey are r < 19.4 mag over two thirds and 19.8 mag over one third of the area. The GSMF is determined from a sample of 5210 galaxies using a density-corrected maximum volume method. This efficiently overcomes the issue of fluctuations in the number density versus redshift. With H_0 = 70, the GSMF is well described between 10^8 and 10^11.5 Msun using a double Schechter function with mass^* = 10^10.66 Msun, phi_1^* = 3.96 x 10^-3 Mpc^-3, alpha_1 = -0.35, phi_2^* = 0.79 x 10^-3 Mpc^-3 and alpha_2 = -1.47. This result is more robust to uncertainties in the flow-model corrected redshifts than from the shallower Sloan Digital Sky Survey main sample (r < 17.8 mag). The upturn in the GSMF is also seen directly in the i-band and K-band galaxy luminosity functions. Accurately measuring the GSMF below 10^8 Msun is possible within the GAMA survey volume but as expected requires deeper imaging data to address the contribution from low surface-brightness galaxies.

We present single-S\'ersic two-dimensional model fits to 167,600 galaxies modelled independently in the ugrizYJHK bandpasses using reprocessed Sloan Digital Sky Survey Data Release Seven (SDSS DR7) and UKIRT Infrared Deep Sky Survey Large Area Survey (UKIDSS-LAS) imaging data available from the GAMA database. In order to facilitate this study we developed SIGMA, an R wrapper around several contemporary astronomy software packages including Source Extractor, PSF Extractor and GALFIT 3. SIGMA produces realistic 2D model fits to galaxies, employing automatic adaptive background subtraction and empirical PSF measurements on the fly for each galaxy in GAMA. Using these results, we define a common coverage area across the three GAMA regions containing 138,269 galaxies. We provide S\'ersic magnitudes truncated at 10 re which show good agreement with SDSS Petrosian and GAMA photometry for low S\'ersic index systems (n < 4), and much improved photometry for high S\'ersic index systems (n > 4), recovering as much as \Delta m = 0.5 magnitudes in the r band. We employ a K band S\'ersic index/u - r colour relation to delineate the massive (n > ~2) early-type galaxies (ETGs) from the late-type galaxies (LTGs). The mean S\'ersic index of these ETGs shows a smooth variation with wavelength, increasing by 30% from g through K. LTGs exhibit a more extreme change in S\'ersic index, increasing by 52% across the same range. In addition, ETGs and LTGs exhibit a 38% and 25% decrease respectively in half-light radius from g through K. These trends are shown to arise due to the effects of dust attenuation and stellar population/metallicity gradients within galaxy populations.

Galaxy and Mass Assembly (GAMA) is a project to study galaxy formation and evolution, combining imaging data from ultraviolet to radio with spectroscopic data from the AAOmega spectrograph on the Anglo-Australian Telescope. Using data from phase 1 of GAMA, taken over three observing seasons, and correcting for various minor sources of incompleteness, we calculate galaxy luminosity functions (LFs) and their evolution in the ugriz passbands. At low redshift, z < 0.1, we find that blue galaxies, defined according to a magnitude-dependent but non-evolving colour cut, are reasonably well fit over a range of more than ten magnitudes by simple Schechter functions in all bands. Red galaxies, and the combined blue-plus-red sample, require double power-law Schechter functions to fit a dip in their LF faintward of the characteristic magnitude M* before a steepening faint end. This upturn is at least partly due to dust-reddened disk galaxies. We measure evolution of the galaxy LF over the redshift range 0.002 < z < 0.5 using both a parametric fit and by measuring binned LFs in redshift slices. The characteristic luminosity L* is found to increase with redshift in all bands, with red galaxies showing stronger luminosity evolution than blue galaxies. The comoving number density of blue galaxies increases with redshift, while that of red galaxies decreases, consistent with prevailing movement from blue cloud to red sequence. As well as being more numerous at higher redshift, blue galaxies also dominate the overall luminosity density beyond redshifts z = 0.2. At lower redshifts, the luminosity density is dominated by red galaxies in the riz bands, by blue galaxies in u and g.

We present the dust properties and star-formation histories of local submillimetre-selected galaxies in Herschel-ATLAS, classified by optical morphology. The early-type galaxies (ETGs) that are detected contain as much dust as typical spirals, and form a unique sample that has been blindly selected at submillimetre wavelengths. Comparing H-ATLAS galaxies to a control sample of optically selected galaxies, we find 5.5% of luminous ETGs are detected in H-ATLAS. The H-ATLAS ETGs contain a significant mass of cold dust: the mean dust mass is 5.5x10^7 Msun, with individual galaxies ranging from 9x10^5-4x10^8 Msun. This is comparable to that of spirals in our sample, and is an order of magnitude more dust than that found for the control ETGs, which have a median dust mass inferred from stacking of (0.8-4.0)x10^6 Msun. The ETGs detected in H-ATLAS have bluer NUV-r colours, higher specific star-formation rates and younger stellar populations than ETGs which are optically selected, and may be transitioning from the blue cloud to the red sequence. We also find that H-ATLAS and control ETGs inhabit similar low-density environments. We conclude that the dust in H-ATLAS and control ETGs cannot be solely from stellar sources, and a large contribution from dust formed in the ISM or external sources is required. Alternatively, dust destruction may not be as efficient as predicted. We also explore the properties of the most passive spiral galaxies in our sample with SSFR<10^-11/yr. We find these passive spirals have lower dust-to-stellar mass ratios, higher stellar masses and older stellar population ages than normal spirals. The passive spirals inhabit low density environments similar to those of the normal spiral galaxies in our sample. This shows that the processes which turn spirals passive do not occur solely in the intermediate density environments of group and cluster outskirts. (Abridged)

This paper describes the first catalogue of photometrically-derived stellar mass estimates for intermediate-redshift (z < 0.65) galaxies in the Galaxy And Mass Assembly (GAMA) spectroscopic redshift survey. These masses, as well as the full set of ancillary stellar population parameters, will be made public as part of GAMA data release 2. Although the GAMA database does include NIR photometry, we show that the quality of our stellar population synthesis fits is significantly poorer when these NIR data are included. Further, for a large fraction of galaxies, the stellar population parameters inferred from the optical-plus-NIR photometry are formally inconsistent with those inferred from the optical data alone. This may indicate problems in our stellar population library, or NIR data issues, or both; these issues will be addressed for future versions of the catalogue. For now, we have chosen to base our stellar mass estimates on optical photometry only. In light of our decision to ignore the available NIR data, we examine how well stellar mass can be constrained based on optical data alone. We use generic properties of stellar population synthesis models to demonstrate that restframe colour alone is in principle a very good estimator of stellar mass-to-light ratio, M*/Li. Further, we use the observed relation between restframe (g-i) and M*/Li for real GAMA galaxies to argue that, modulo uncertainties in the stellar evolution models themselves, (g-i) colour can in practice be used to estimate M*/Li to an accuracy of < ~0.1 dex. This 'empirically calibrated' (g-i)-M*/Li relation offers a simple and transparent means for estimating galaxies' stellar masses based on minimal data, and so provides a solid basis for other surveys to compare their results to z < ~0.4 measurements from GAMA.

We present a technique to identify optical counterparts of 250 um-selected sources from the Herschel-ATLAS survey. Of the 6621 250 um > 32 mJy sources in our science demonstration catalogue we find that ~60 percent have counterparts brighter than r=22.4 mag in the Sloan Digital Sky Survey. Applying a likelihood ratio technique we are able to identify 2423 of the counterparts with a reliability R > 0.8. This is approximately 37 percent of the full 250 micron catalogue. We have estimated photometric redshifts for each of these 2423 reliable counterparts, while 1099 also have spectroscopic redshifts collated from several different sources, including the GAMA survey. We estimate the completeness of identifying counterparts as a function of redshift, and present evidence that 250 um-selected Herschel-ATLAS galaxies have a bimodal redshift distribution. Those with reliable optical identifications have a redshift distribution peaking at z ~ 0.25 +/- 0.05, while sub-mm colours suggest that a significant fraction with no counterpart above the r-band limit have z > 1. We also suggest a method for selecting populations of strongly-lensed high redshift galaxies. Our identifications are matched to UV--NIR photometry from the GAMA survey, and these data are available as part of the Herschel-ATLAS public data release.

We investigate the properties of satellite galaxies that surround isolated hosts within the redshift range 0.01 < z < 0.15, using data taken as part of the Galaxy And Mass Assembly survey. Making use of isolation and satellite criteria that take into account stellar mass estimates, we find 3514 isolated galaxies of which 1426 host a total of 2998 satellites. Separating the red and blue populations of satellites and hosts, using colour-mass diagrams, we investigate the radial distribution of satellite galaxies and determine how the red fraction of satellites varies as a function of satellite mass, host mass and the projected distance from their host. Comparing the red fraction of satellites to a control sample of small neighbours at greater projected radii, we show that the increase in red fraction is primarily a function of host mass. The satellite red fraction is about 0.2 higher than the control sample for hosts with 11.0 < log M_* < 11.5, while the red fractions show no difference for hosts with 10.0 < log M_* < 10.5. For the satellites of more massive hosts the red fraction also increases as a function of decreasing projected distance. Our results suggest that the likely main mechanism for the quenching of star formation in satellites hosted by isolated galaxies is strangulation.

We investigate the ultraviolet and optical properties and environment of low redshift galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration data. We use the Sloan Digital Sky Survey seventh release and the Galaxy And Mass Assembly database to select galaxies with r_Petro < 19.0 mag in the redshift range 0.02 < z < 0.2 and look for their submillimeter counterparts in H-ATLAS. Our results show that at low redshift, H-ATLAS detects mainly blue/star-forming galaxies with a minor contribution from red systems which are highly obscured by dust. In addition we find that the colour of a galaxy rather than the local density of its environment determines whether it is detectable by H-ATLAS. The average dust temperature of galaxies that are simultaneously detected by both PACS and SPIRE is 25K \pm 4K, independent of environment. This analysis provides a glimpse of the potential of the H-ATLAS data to investigate the submillimeter properties of galaxies in the local universe.

Using the complete GAMA-I survey covering ~142 sq. deg. to r=19.4, of which ~47 sq. deg. is to r=19.8, we create the GAMA-I galaxy group catalogue (G3Cv1), generated using a friends-of-friends (FoF) based grouping algorithm. Our algorithm has been tested extensively on one family of mock GAMA lightcones, constructed from Lambda-CDM N-body simulations populated with semi-analytic galaxies. Recovered group properties are robust to the effects of interlopers and are median unbiased in the most important respects. G3Cv1 contains 14,388 galaxy groups (with multiplicity >= 2$), including 44,186 galaxies out of a possible 110,192 galaxies, implying ~40% of all galaxies are assigned to a group. The similarities of the mock group catalogues and G3Cv1 are multiple: global characteristics are in general well recovered. However, we do find a noticeable deficit in the number of high multiplicity groups in GAMA compared to the mocks. Additionally, despite exceptionally good local spatial completeness, G3Cv1 contains significantly fewer compact groups with 5 or more members, this effect becoming most evident for high multiplicity systems. These two differences are most likely due to limitations in the physics included of the current GAMA lightcone mock. Further studies using a variety of galaxy formation models are required to confirm their exact origin.

The stellar initial mass function (IMF) describes the distribution in stellar masses produced from a burst of star formation. For more than fifty years, the implicit assumption underpinning most areas of research involving the IMF has been that it is universal, regardless of time and environment. We measure the high-mass IMF slope for a sample of low-to-moderate redshift galaxies from the Galaxy And Mass Assembly survey. The large range in luminosities and galaxy masses of the sample permits the exploration of underlying IMF dependencies. A strong IMF-star formation rate dependency is discovered, which shows that highly star forming galaxies form proportionally more massive stars (they have IMFs with flatter power-law slopes) than galaxies with low star formation rates. This has a significant impact on a wide variety of galaxy evolution studies, all of which rely on assumptions about the slope of the IMF. Our result is supported by, and provides an explanation for, the results of numerous recent explorations suggesting a variation of or evolution in the IMF.

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

We present an analysis of the properties of the lowest Halpha-luminosity galaxies (L_Halpha<4x10^32 W; SFR<0.02 Msun/yr) in the Galaxy And Mass Assembly (GAMA) survey. These galaxies make up the the rise above a Schechter function in the number density of systems seen at the faint end of the Halpha luminosity function. Above our flux limit we find that these galaxies are principally composed of intrinsically low stellar mass systems (median stellar mass =2.5x10^8 Msun) with only 5/90 having stellar masses M>10^10 Msun. The low SFR systems are found to exist predominantly in the lowest density environments (median density ~0.02 galaxy Mpc^-2 with none in environments more dense than ~1.5 galaxy Mpc^-2). Their current specific star formation rates (SSFR; -8.5 < log(SSFR[yr^-1])<-12.) are consistent with their having had a variety of star formation histories. The low density environments of these galaxies demonstrates that such low-mass, star-forming systems can only remain as low-mass and forming stars if they reside sufficiently far from other galaxies to avoid being accreted, dispersed through tidal effects or having their gas reservoirs rendered ineffective through external processes.

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

We use the Herschel-ATLAS science demonstration data to investigate the star-formation properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminosity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, we find that some individual objects are associated with high far-infrared luminosities. However, the far-infrared properties of the radio-loud population are statistically indistinguishable from those of a comparison population of radio-quiet galaxies matched in redshift and K-band absolute magnitude. There is thus no evidence that the host galaxies of these largely low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a population, have particularly unusual star-formation histories. Models in which the AGN activity in higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict a luminosity-dependent effect that is not seen in our data (which only span a limited range in radio luminosity) but which may well be detectable with the full Herschel-ATLAS dataset.

We use data from the Herschel-ATLAS to investigate the evolution of the far-infrared--radio correlation over the redshift range 0<z<0.5. Using the total far-infrared luminosity of all >5sigma sources in the Herschel-ATLAS Science Demonstration Field and cross-matching these data with radio data from the Faint Images of the Radio Sky at Twenty-Centimetres (FIRST) survey and the NRAO VLA Northern Sky Survey (NVSS), we obtain 104 radio counterparts to the Herschel sources. With these data we find no evidence for evolution in the far-infrared--radio correlation over the redshift range 0<z<0.5, where the median value for the ratio between far-infrared and radio luminosity, $q_{\rm IR}$, over this range is $q_{\rm IR} = 2.40\pm 0.12$ (and a mean of $q_{\rm IR}=2.52 \pm 0.03$ accounting for the lower limits), consistent with both the local value determined from {\em IRAS} and values derived from surveys targeting the high-redshift Universe. By comparing the radio fluxes of our sample measured from both FIRST and NVSS we show that previous results suggesting an increase in the value of $q_{\rm IR}$ from high to low redshift may be the result of resolving out extended emission of the low-redshift sources with relatively high-resolution interferometric data, although AGN contamination could still play a significant role. We also find tentative evidence that the longer wavelength, cooler dust is heated by an evolved stellar population which does not trace the star-formation rate as closely as the shorter wavelength $\ltsim 250~\mu$m emission or the radio emission, supporting suggestions based on detailed models of individual galaxies.

We present self-consistent star formation rates derived through pan-spectral analysis of galaxies drawn from the Galaxy and Mass Assembly (GAMA) survey. We determine the most appropriate form of dust obscuration correction via application of a range of extinction laws drawn from the literature as applied to Halpha, [O{II}] and UV luminosities. These corrections are applied to a sample of 31,508 galaxies from the GAMA survey at z < 0.35. We consider several different obscuration curves, including those of Milky Way, Calzetti (2001) and Fischera and Dopita (2005) curves and their effects on the observed luminosities. At the core of this technique is the observed Balmer decrement, and we provide a prescription to apply optimal obscuration corrections using the Balmer decrement. We carry out an analysis of the star formation history (SFH) using stellar population synthesis tools to investigate the evolutionary history of our sample of galaxies as well as to understand the effects of variation in the Initial Mass Function (IMF) and the effects this has on the evolutionary history of galaxies. We find that the Fischera and Dopita (2005) obscuration curve with an R_{v} value of 4.5 gives the best agreement between the different SFR indicators. The 2200A feature needed to be removed from this curve to obtain complete consistency between all SFR indicators suggesting that this feature may not be common in the average integrated attenuation of galaxy emission. We also find that the UV dust obscuration is strongly dependent on the SFR.

The Galaxy And Mass Assembly (GAMA) survey has been operating since February 2008 on the 3.9-m Anglo-Australian Telescope using the AAOmega fibre-fed spectrograph facility to acquire spectra with a resolution of R~1300 for 120,862 SDSS selected galaxies. The target catalogue constitutes three contiguous equatorial regions centred at 9h (G09), 12h (G12) and 14.5h (G15) each of 12 x 4 sq.deg to limiting fluxes of r < 19.4, r < 19.8, and r < 19.4 mag respectively (and additional limits at other wavelengths). Spectra and reliable redshifts have been acquired for over 98 per cent of the galaxies within these limits. Here we present the survey footprint, progression, data reduction, redshifting, re-redshifting, an assessment of data quality after 3 years, additional image analysis products (including ugrizYJHK photometry, Sersic profiles and photometric redshifts), observing mask, and construction of our core survey catalogue (GamaCore). From this we create three science ready catalogues: GamaCoreDR1 for public release, which includes data acquired during year 1 of operations within specified magnitude limits (February 2008 to April 2008); GamaCoreMainSurvey containing all data above our survey limits for use by the GAMA team and collaborators; and GamaCoreAtlasSv containing year 1, 2 and 3 data matched to Herschel-ATLAS Science Demonstration data. These catalogues along with the associated spectra, stamps and profiles can be accessed via the GAMA website: http://www.gama-survey.org/

In order to generate credible 0.1-2 {\mu}m SEDs, the GAMA project requires many Gigabytes of imaging data from a number of instruments to be re-processed into a standard format. In this paper we discuss the software infrastructure we use, and create self-consistent ugrizYJHK photometry for all sources within the GAMA sample. Using UKIDSS and SDSS archive data, we outline the pre-processing necessary to standardise all images to a common zeropoint, the steps taken to correct for seeing bias across the dataset, and the creation of Gigapixel-scale mosaics of the three 4x12 deg GAMA regions in each filter. From these mosaics, we extract source catalogues for the GAMA regions using elliptical Kron and Petrosian matched apertures. We also calculate S\'ersic magnitudes for all galaxies within the GAMA sample using SIGMA, a galaxy component modelling wrapper for GALFIT 3. We compare the resultant photometry directly, and also calculate the r band galaxy LF for all photometric datasets to highlight the uncertainty introduced by the photometric method. We find that (1) Changing the object detection threshold has a minor effect on the best-fitting Schechter parameters of the overall population (M* +/- 0.055mag, {\alpha} +/- 0.014, {\Phi}* +/- 0.0005 h^3 Mpc^{-3}). (2) An offset between datasets that use Kron or Petrosian photometry regardless of the filter. (3) The decision to use circular or elliptical apertures causes an offset in M* of 0.20mag. (4) The best-fitting Schechter parameters from total-magnitude photometric systems (such as SDSS modelmag or S\'ersic magnitudes) have a steeper faint-end slope than photometry dependent on Kron or Petrosian magnitudes. (5) Our Universe's total luminosity density, when calculated using Kron or Petrosian r-band photometry, is underestimated by at least 15%.

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

We describe the spectroscopic target selection for the Galaxy And Mass Assembly (GAMA) survey. The input catalogue is drawn from the Sloan Digital Sky Survey (SDSS) and UKIRT Infrared Deep Sky Survey (UKIDSS). The aim is to measure redshifts for galaxies in three 4x12 degree regions at 9h, 12h and 14.5h, on the celestial equator, with magnitude selections r < 19.4, z < 18.2 and K(AB) < 17.6 over all three regions, and r < 19.8 in the 12-h region. The target density is 1080 per sq. deg. in the 12-h region and 720 per sq. deg. in the other regions. The average GAMA target density and area are compared with completed and ongoing galaxy redshift surveys. The GAMA survey implements a highly complete star-galaxy separation that jointly uses an intensity-profile separator (delta(sg) = r-band psf mag - model mag) as per the SDSS) and a colour separator. The colour separator is defined as delta(sg,jk) = J - K - f(g-i), where f(g-i) is a quadratic fit to the J-K colour of the stellar locus over the range 0.3 < g-i < 2.3. All galaxy populations investigated are well separated with delta(sg,jk) > 0.2. From two years out of a three-year AAOmega program on the Anglo-Australian Telescope, we have obtained 79599 unique galaxy redshifts. Previously known redshifts in the GAMA region bring the total up to 98497. The median galaxy redshift is 0.2 with 99% at z < 0.5. We present some of the global statistical properties of the survey, including K-band galaxy counts, colour-redshift relations and preliminary n(z).

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

The Galaxy And Mass Assembly (GAMA) project is the latest in a tradition of large galaxy redshift surveys, and is now underway on the 3.9m Anglo-Australian Telescope at Siding Spring Observatory. GAMA is designed to map extragalactic structures on scales of 1kpc - 1Mpc in complete detail to a redshift of z~0.2, and to trace the distribution of luminous galaxies out to z~0.5. The principal science aim is to test the standard hierarchical structure formation paradigm of Cold Dark Matter (CDM) on scales of galaxy groups, pairs, discs, bulges and bars. We will measure (1) the Dark Matter Halo Mass Function (as inferred from galaxy group velocity dispersions); (2) baryonic processes, such as star formation and galaxy formation efficiency (as derived from Galaxy Stellar Mass Functions); and (3) the evolution of galaxy merger rates (via galaxy close pairs and galaxy asymmetries). Additionally, GAMA will form the central part of a new galaxy database, which aims to contain 275,000 galaxies with multi-wavelength coverage from coordinated observations with the latest international ground- and space-based facilities: GALEX, VST, VISTA, WISE, HERSCHEL, GMRT and ASKAP. Together, these data will provide increased depth (over 2 magnitudes), doubled spatial resolution (0.7"), and significantly extended wavelength coverage (UV through Far-IR to radio) over the main SDSS spectroscopic survey for five regions, each of around 50 deg^2. This database will permit detailed investigations of the structural, chemical, and dynamical properties of all galaxy types, across all environments, and over a 5 billion year timeline.

A heuristic greedy algorithm is developed for efficiently tiling spatially dense redshift surveys. In its first application to the Galaxy and Mass Assembly (GAMA) redshift survey we find it rapidly improves the spatial uniformity of our data, and naturally corrects for any spatial bias introduced by the 2dF multi object spectrograph. We make conservative predictions for the final state of the GAMA redshift survey after our final allocation of time, and can be confident that even if worse than typical weather affects our observations, all of our main survey requirements will be met.

We present luminosity and surface brightness distributions of 40,111 galaxies with K-band photometry from the United Kingdom Infrared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS), Data Release 3 and optical photometry from Data Release 5 of the Sloan Digital Sky Survey (SDSS). Various features and limitations of the new UKIDSS data are examined, such as a problem affecting Petrosian magnitudes of extended sources. Selection limits in K- and r-band magnitude, K-band surface brightness and K-band radius are included explicitly in the 1/Vmax estimate of the space density and luminosity function. The bivariate brightness distribution in K-band absolute magnitude and surface brightness is presented and found to display a clear luminosity--surface brightness correlation that flattens at high luminosity and broadens at low luminosity, consistent with similar analyses at optical wavelengths. Best fitting Schechter function parameters for the K-band luminosity function are found to be M*-5 log h=-23.19 +/- 0.04, alpha=-0.81 +/- 0.04 and phi*=(0.0166 +/- 0.0008)h^3 Mpc^{-3}, although the Schechter function provides a poor fit to the data at high and low luminosity, while the luminosity density in the K band is found to be j = (6.305 +/- 0.067) x 10^8 L_sun h Mpc^{-3}. However, we caution that there are various known sources of incompleteness and uncertainty in our results. Using mass-to-light ratios determined from the optical colours we estimate the stellar mass function, finding good agreement with previous results. Possible improvements are discussed that could be implemented when extending this analysis to the full LAS.

The nature of galaxy structures on large scales is a key observational prediction for current models of galaxy formation. The SDSS and 2dF galaxy surveys have revealed a number of structures on 40-150 h^-1 Mpc scales at low redshifts, and some even larger ones. To constrain galaxy number densities, luminosities, and stellar populations in large structures at higher redshift, we have investigated two sheet-like structures of galaxies at z=0.8 and 1.3 spanning 150 h^-1 comoving Mpc embedded in large quasar groups extending over at least 200 h^-1 Mpc. We present first results of an analysis of these sheet--like structures using two contiguous 1deg GALEX fields (FUV and NUV) cross-correlated with optical data from the Sloan Digital Sky Survey (SDSS). We derive a sample of 462 Lyman Break Galaxy (LBG) candidates coincident with the sheets. Using the GALEX and SDSS data, we show that the overall average spectral energy distribution of a LBG galaxy at z~1 is flat (in f_lambda) in the rest frame wavelength range from 1500A, to 4000A, implying evolved populations of stars in the LBGs. From the luminosity functions we get indications for overdensities in the two LQGs compared to their foreground regions. Similar conclusions come from the calculation of the 2-point correlation function, showing a 2sigma overdensity for the LBGs in the z~0.8 LQG on scales of 1.6 to 4.8 Mpc, indicating similar correlation scales for our LBG sample as their z~3 counterparts.

We present the final spectroscopic QSO catalogue from the 2dF-SDSS LRG and QSO (2SLAQ) Survey. This is a deep, 18<g<21.85 (extinction corrected), sample aimed at probing in detail the faint end of the broad line AGN luminosity distribution at z<2.6. The candidate QSOs were selected from SDSS photometry and observed spectroscopically with the 2dF spectrograph on the Anglo-Australian Telescope. This sample covers an area of 191.9 deg^2 and contains new spectra of 16326 objects, of which 8764 are QSOs, and 7623 are newly discovered (the remainder were previously identified by the 2QZ and SDSS surveys). The full QSO sample (including objects previously observed in the SDSS and 2QZ surveys) contains 12702 QSOs. The new 2SLAQ spectroscopic data set also contains 2343 Galactic stars, including 362 white dwarfs, and 2924 narrow emission line galaxies with a median redshift of z=0.22. We present detailed completeness estimates for the survey, based on modelling of QSO colours, including host galaxy contributions. This calculation shows that at g~21.85 QSO colours are significantly affected by the presence of a host galaxy up to redshift z~1 in the SDSS ugriz bands. In particular we see a significant reddening of the objects in g-i towards fainter g-band magnitudes. This reddening is consistent with the QSO host galaxies being dominated by a stellar population of age at least 2-3 Gyr. The full catalogue, including completeness estimates, is available on-line at http://www.2slaq.info/

We present an analysis of the small-to-intermediate scale clustering of samples of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey and the 2dF-SDSS LRG and QSO (2SLAQ) survey carefully matched to have the same rest-frame colours and luminosity. We study the spatial two-point auto-correlation function in both redshift-space and real-space of a combined sample of over 10,000 LRGs, which represent the most massive galaxies in the universe with stellar masses > 10^11 h^-1 M_sun and space densities 10^-4 h^-3 Mpc^-3. We find no significant evolution in the amplitude r_0 of the correlation function with redshift, but do see a slight decrease in the slope with increasing redshift over 0.19 < z < 0.55 and scales of 0.32 < r < 32 h^-1 Mpc. We compare our measurements with the predicted evolution of dark matter clustering and use the halo model to interpret our results. We find that our clustering measurements are inconsistent (>99.9% significance) with a passive model whereby the LRGs do not merge with one another; a model with a merger rate of 7.5 +/- 2.3% from z = 0.55 to z = 0.19 (i.e. an average rate of 2.4% Gyr^-1) provides a better fit to our observations. Our clustering and number density measurements are consistent with the hypothesis that the merged LRGs were originally central galaxies in different haloes which, following the merger of these haloes, merged to create a single Brightest Cluster Galaxy. In addition, we show that the small-scale clustering signal constrains the scatter in halo merger histories. When combined with measurements of the luminosity function, our results suggest that this scatter is sub-Poisson. While this is a generic prediction of hierarchical models, it has not been tested before.

We use the Fourth Data Release of the Sloan Digital Sky Survey to investigate the relation between galaxy rest frame u-r colour, morphology, as described by the concentration and Sersic indices, and environmental density, for a sample of 79,553 galaxies at z < ~0.1. We split the samples according to density and luminosity and recover the expected bimodal distribution in the colour-morphology plane, shown especially clearly by this subsampling. We quantify the bimodality by a sum of two Gaussians on the colour and morphology axes and show that, for the red/early-type population both colour and morphology do not change significantly as a function of density. For the blue/late-type population, with increasing density the colour becomes redder but the morphology again does not change significantly. Both populations become monotonically redder and of earlier type with increasing luminosity. There is no significant qualitative difference between the behaviour of the two morphological measures. We supplement the morphological sample with 13,655 galaxies assigned Hubble types by an artificial neural network. We find, however, that the resulting distribution is less well described by two Gaussians. Therefore, there are either more than two significant morphological populations, physical processes not seen in colour space, or the Hubble type, particularly the different subtypes of spirals Sa-Sd, has an irreducible fuzziness when related to environmental density. For each of the three measures of morphology, on removing the density relation due to it, we recover a strong residual relation in colour. However, on similarly removing the colour-density relation there is no evidence for a residual relation due to morphology. [Abridged]

We present a clustering analysis of Luminous Red Galaxies (LRGs) using nearly 9 000 objects from the final catalogue of the 2dF-SDSS LRG And QSO (2SLAQ) Survey. We measure the redshift-space two-point correlation function, xi(s), at the mean LRG redshift of z=0.55. A single power-law fits the deprojected correlation function, xi(r), with a correlation length of r_0=7.45+-0.35 Mpc and a power-law slope of gamma=1.72+-0.06 in the 0.4<r<50 Mpc range. But it is in the LRG angular correlation function that the strongest evidence for non-power-law features is found where a slope of gamma=-2.17+-0.07 is seen at 1<r<10 Mpc with a flatter gamma=-1.67+-0.03 slope apparent at r<~1 Mpc scales. We use the simple power-law fit to the galaxy xi(r) to model the redshift space distortions in the 2-D redshift-space correlation function, xi(sigma,pi). We fit for the LRG velocity dispersion, w_z, Omega_m and beta, where beta=Omega_m^0.6/b and b is the linear bias parameter. We find values of w_z=330kms^-1, Omega_m= 0.10+0.35-0.10 and beta=0.40+-0.05. These high redshift results, which incorporate the Alcock-Paczynski effect and the effects of dynamical infall, start to break the degeneracy between Omega_m and beta found in low-redshift galaxy surveys. This degeneracy is further broken by introducing an additional external constraint, the value of beta(z=0.1)=0.45 from 2dFGRS, and then considering the evolution of clustering from z~0 to z_LRG~0.55. With these combined methods we find Omega_m(z=0)=0.30+-0.15 and beta(z=0.55)=0.45+-0.05. Assuming these values, we find a value for b(z=0.55)=1.66+-0.35. We show that this is consistent with a simple ``high peaks'' bias prescription which assumes that LRGs have a constant co-moving density and their clustering evolves purely under gravity. [ABRIDGED]

We use empirical techniques to interpret the near-infrared colours of a sample of 5800 galaxies drawn from Sloan Digital Sky Survey (SDSS) main spectroscopic sample with YJHK photometry from the UK Infrared Deep Sky Survey (UKIDSS) data release one. We study correlations between near-IR colours measured within SDSS fibre and physical parameters derived from the spectra. These parameters include specific star formation rate, stellar age, metallicity and dust attenuation. All correlations are analyzed for samples of galaxies that are closely matched in redshift, in stellar mass and in concentration index. Whereas more strongly star-forming galaxies have bluer optical colours, the opposite is true at near-IR wavelengths -- galaxies with higher specific star formation rate have redder near-IR colours. This result agrees qualitatively with the predictions of models in which Thermally Pulsing Asymptotic Giant Branch (TP-AGB) stars dominate the H and K-band light of a galaxy following a burst of star formation. We also find a surprisingly strong correlation between the near-IR colours of star-forming galaxies and their dust attenuation as measured from the Balmer decrement. Unlike optical colours, however, near-IR colours exhibit very little dependence on galaxy inclination. This suggests that the correlation of near-IR colours with dust attenuation arises because TP-AGB stars are the main source of dust in the galaxy. Finally, we compare the near-IR colours of the galaxies in our sample to the predictions of three different stellar population models: the Bruzual & Charlot 2003 model, a preliminary version of a new model under development by Charlot & Bruzual, which includes a new prescription for AGB star evolution, and the Maraston 2005 model.

We describe the construction of MegaZ-LRG, a photometric redshift catalogue of over one million luminous red galaxies (LRGs) in the redshift range 0.4 < z < 0.7 with limiting magnitude i < 20. The catalogue is selected from the imaging data of the Sloan Digital Sky Survey Data Release 4. The 2dF-SDSS LRG and Quasar (2SLAQ) spectroscopic redshift catalogue of 13,000 intermediate-redshift LRGs provides a photometric redshift training set, allowing use of ANNz, a neural network-based photometric-redshift estimator. The rms photometric redshift accuracy obtained for an evaluation set selected from the 2SLAQ sample is sigma_z = 0.049 averaged over all galaxies, and sigma_z = 0.040 for a brighter subsample (i < 19.0). The catalogue is expected to contain ~5 per cent stellar contamination. The ANNz code is used to compute a refined star/galaxy probability based on a range of photometric parameters; this allows the contamination fraction to be reduced to 2 per cent with negligible loss of genuine galaxies. The MegaZ-LRG catalogue is publicly available on the World Wide Web from http://www.2slaq.info .

We measure the cosmological matter density by observing the positions of baryon acoustic oscillations in the clustering of galaxies in the Sloan Digital Sky Survey (SDSS). We jointly analyse the main galaxies and LRGs in the SDSS DR5 sample, using over half a million galaxies in total. The oscillations are detected with 99.74% confidence (3.0sigma assuming Gaussianity) compared to a smooth power spectrum. When combined with the observed scale of the peaks within the CMB, we find a best-fit value of Omega_m=0.256+0.029-0.024 (68% confidence interval), for a flat Lambda cosmology when marginalising over the Hubble parameter and the baryon density. This value of the matter density is derived from the locations of the baryon oscillations in the galaxy power spectrum and in the CMB, and does not include any information from the overall shape of the power spectra. This is an extremely clean cosmological measurement as the physics of the baryon acoustic oscillation production is well understood, and the positions of the oscillations are expected to be independent of systematics such as galaxy bias.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

We present a Fourier analysis of the clustering of galaxies in the combined Main galaxy and Luminous Red Galaxy (LRG) Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) sample. The aim of our analysis is to consider how well we can measure the cosmological matter density using the signature of the horizon at matter-radiation equality embedded in the large-scale power spectrum. The new data constrains the power spectrum on scales 100--600h^-1Mpc with significantly higher precision than previous analyses of just the SDSS Main galaxies, due to our larger sample and the inclusion of the LRGs. This improvement means that we can now reveal a discrepancy between the shape of the measured power and linear CDM models on scales 0.01<k<0.15hMpc^-1, with linear model fits favouring a lower matter density (Omega_m=0.22+/-0.04) on scales 0.01<k<0.06hMpc^-1 and a higher matter density (Omega_m=0.32+/-0.01) when smaller scales are included, assuming a flat LCDM model with h=0.73 and n_s=0.96. This discrepancy could be explained by scale-dependent bias and, by analysing subsamples of galaxies, we find that the ratio of small-scale to large-scale power increases with galaxy luminosity, so all of the SDSS galaxies cannot trace the same power spectrum shape over 0.01<k<0.2hMpc^-1. However, the data are insufficient to clearly show a luminosity-dependent change in the largest scale at which a significant increase in clustering is observed, although they do not rule out such an effect. Significant scale-dependent galaxy bias on large-scales, which changes with the r-band luminosity of the galaxies, could potentially explain differences in our Omega_m estimates and differences previously observed between 2dFGRS and SDSS power spectra and the resulting parameter constraints.

Bivariate luminosity functions (LFs) are computed for galaxies in the New York Value-Added Galaxy Catalogue, based on the Sloan Digital Sky Survey Data Release 4. The galaxy properties investigated are the morphological type, inverse concentration index, Sersic index, absolute effective surface brightness, reference frame colours, absolute radius, eClass spectral type, stellar mass and galaxy environment. The morphological sample is flux-limited to galaxies with r < 15.9 and consists of 37,047 classifications to an RMS accuracy of +/- half a class in the sequence E, S0, Sa, Sb, Sc, Sd, Im. These were assigned by an artificial neural network, based on a training set of 645 eyeball classifications. The other samples use r < 17.77 with a median redshift of z ~ 0.08, and a limiting redshift of z < 0.15 to minimize the effects of evolution. Other cuts, for example in axis ratio, are made to minimize biases. A wealth of detail is seen, with clear variations between the LFs according to absolute magnitude and the second parameter. They are consistent with an early type, bright, concentrated, red population and a late type, faint, less concentrated, blue, star forming population. This bimodality suggests two major underlying physical processes, which in agreement with previous authors we hypothesize to be merger and accretion, associated with the properties of bulges and discs respectively. The bivariate luminosity-surface brightness distribution is fit with the Choloniewski function (a Schechter function in absolute magnitude and Gaussian in surface brightness). The fit is found to be poor, as might be expected if there are two underlying processes.

We present a detailed investigation into the recent star formation histories of 5,697 Luminous Red Galaxies (LRGs) based on the Hdelta (4101A) and [OII] (3727A) lines. LRGs are luminous (L>3L*), galaxies which have been selected to have photometric properties consistent with an old, passively evolving stellar population. For this study we utilise LRGs from the recently completed 2dF-SDSS LRG and QSO survey (2SLAQ). Equivalent widths of the Hdelta and [OII] lines are measured and used to define three spectral types, those with only strong Hdelta absorption (k+a), those with strong [OII] in emission (em) and those with both (em+a). All other LRGs are considered to have passive star formation histories. The vast majority of LRGs are found to be passive (~80 per cent), however significant numbers of k+a (2.7 per cent), em+a (1.2 per cent) and em LRGs (8.6 per cent) are identified. An investigation into the redshift dependence of the fractions is also performed. A sample of SDSS MAIN galaxies with colours and luminosities consistent with the 2SLAQ LRGs is selected to provide a low redshift comparison. While the em and em+a fractions are consistent with the low redshift SDSS sample, the fraction of k+a LRGs is found to increase significantly with redshift. This result is interpreted as an indication of an increasing amount of recent star formation activity in LRGs with redshift. By considering the expected life time of the k+a phase, the number of LRGs which will undergo a k+a phase can be estimated. A crude comparison of this estimate with the predictions from semi-analytic models of galaxy formation shows that the predicted level of k+a and em+a activity is not sufficient to reconcile the predicted mass growth for massive early-types in a hierarchical merging scenario.

We present new measurements of the luminosity function (LF) of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) and the 2dF-SDSS LRG and Quasar (2SLAQ) survey. We have carefully quantified, and corrected for, uncertainties in the K and evolutionary corrections, differences in the colour selection methods, and the effects of photometric errors, thus ensuring we are studying the same galaxy population in both surveys. Using a limited subset of 6326 SDSS LRGs (with 0.17<z<0.24) and 1725 2SLAQ LRGs (with 0.5 <z<0.6), for which the matching colour selection is most reliable, we find no evidence for any additional evolution in the LRG LF, over this redshift range, beyond that expected from a simple passive evolution model. This lack of additional evolution is quantified using the comoving luminosity density of SDSS and 2SLAQ LRGs, brighter than M_r - 5logh = -22.5, which are 2.51+/-0.03 x 10^-7 L_sun Mpc^-3 and 2.44+/-0.15 x 10^-7 L_sun Mpc^-3 respectively (<10% uncertainty). We compare our LFs to the COMBO-17 data and find excellent agreement over the same redshift range. Together, these surveys show no evidence for additional evolution (beyond passive) in the LF of LRGs brighter than M_r - 5logh = -21 (or brighter than L*). We test our SDSS and 2SLAQ LFs against a simple ``dry merger'' model for the evolution of massive red galaxies and find that at least half of the LRGs at z=0.2 must already have been well-assembled (with more than half their stellar mass) by z=0.6. This limit is barely consistent with recent results from semi-analytical models of galaxy evolution.

We present a spectroscopic survey of almost 15,000 candidate intermediate-redshift Luminous Red Galaxies (LRGs) brighter than i=19.8, observed with 2dF on the Anglo-Australian Telescope. The targets were selected photometrically from the Sloan Digital Sky Survey (SDSS) and lie along two narrow equatorial strips covering 180 sq deg. Reliable redshifts were obtained for 92% of the targets and the selection is very efficient: over 90% have redshifts between 0.45 and 0.8. More than 80% of the ~11,000 red galaxies have pure absorption-line spectra consistent with a passively-evolving old stellar population. The redshift, photometric and spatial distributions of the LRGs are described. The 2SLAQ data will be released publicly from mid-2006, providing a powerful resource for observational cosmology and the study of galaxy evolution.

We synthesize the rest-frame Stroemgren colors using SDSS spectra for 99,088 galaxies selected from Data Release 1. This narrow-band ~200 AA photometric system (uz, vz, bz, yz), first designed for the determination of effective temperature, metallicity and gravity of stars, measures the continuum spectral slope of galaxies in the rest-frame 3200-5800 AA wavelength range. Galaxies form a remarkably narrow locus (~0.03 mag) in the resulting color-color diagram. The Bruzual & Charlot population synthesis models suggest that the position of a galaxy along this locus is controlled by a degenerate combination of metallicity and age of the dominant stellar population. Galaxy distribution along the locus is bimodal, with the local minimum corresponding to an ~1 Gyr old single stellar population. The position perpendicular to the locus is independent of metallicity and age, and reflects the galaxy's dust content, as implied by both the models and the statistics of IRAS detections. A comparison of this locus with the galaxy locus in the H_delta-D_n(4000) diagram, utilized by Kauffmann et al. (2003) to estimate stellar masses, reveals a tight correlation, although the two analyzed spectral ranges barely overlap. Overall, the galaxy spectral energy distribution in the entire UV to near-IR range can be described as a single-parameter family with an accuracy of 0.1 mag, or better. This nearly one-dimensional distribution of galaxies in the multi-dimensional space of measured parameters strongly supports the conclusion of Yip et al. (2004), based on a principal component analysis, that SDSS galaxy spectra can be described by a small number of eigenspectra. Apparently, the contributions of stellar populations that dominate the optical emission from galaxies are combined in a simple and well-defined way.

We present the 3D real space clustering power spectrum of a sample of \~600,000 luminous red galaxies (LRGs) measured by the Sloan Digital Sky Survey (SDSS), using photometric redshifts. This sample of galaxies ranges from redshift z=0.2 to 0.6 over 3,528 deg^2 of the sky, probing a volume of 1.5 (Gpc/h)^3, making it the largest volume ever used for galaxy clustering measurements. We measure the angular clustering power spectrum in eight redshift slices and combine these into a high precision 3D real space power spectrum from k=0.005 (h/Mpc) to k=1 (h/Mpc). We detect power on gigaparsec scales, beyond the turnover in the matter power spectrum, on scales significantly larger than those accessible to current spectroscopic redshift surveys. We also find evidence for baryonic oscillations, both in the power spectrum, as well as in fits to the baryon density, at a 2.5 sigma confidence level. The statistical power of these data to constrain cosmology is ~1.7 times better than previous clustering analyses. Varying the matter density and baryon fraction, we find \Omega_M = 0.30 \pm 0.03, and \Omega_b/\Omega_M = 0.18 \pm 0.04, The detection of baryonic oscillations also allows us to measure the comoving distance to z=0.5; we find a best fit distance of 1.73 \pm 0.12 Gpc, corresponding to a 6.5% error on the distance. These results demonstrate the ability to make precise clustering measurements with photometric surveys (abridged).

We determine the number counts and z=0-5 luminosity function for a well-defined, homogeneous sample of quasars from the Sloan Digital Sky Survey (SDSS). We conservatively define the most uniform statistical sample possible, consisting of 15,343 quasars within an effective area of 1622 deg^2 that was derived from a parent sample of 46,420 spectroscopically confirmed broad-line quasars in the 5282 deg^2 of imaging data from SDSS Data Release Three. The sample extends from i=15 to i=19.1 at z<3 and to i=20.2 for z>3. The number counts and luminosity function agree well with the results of the 2dF QSO Survey, but the SDSS data probe to much higher redshifts than does the 2dF sample. The number density of luminous quasars peaks between redshifts 2 and 3, although uncertainties in the selection function in this range do not allow us to determine the peak redshift more precisely. Our best fit model has a flatter bright end slope at high redshift than at low redshift. For z<2.4 the data are best fit by a redshift-independent slope of beta = -3.1 (Phi(L) propto L^beta). Above z=2.4 the slope flattens with redshift to beta=-2.37 at z=5. This slope change, which is significant at a >5-sigma level, must be accounted for in models of the evolution of accretion onto supermassive black holes.

We study the luminosity and color dependence of the galaxy 2-point correlation function in the Sloan Digital Sky Survey, starting from a sample of 200,000 galaxies over 2500 deg^2. We concentrate on the projected correlation function w(r_p), which is directly related to the real space \xi(r). The amplitude of w(r_p) grows continuously with luminosity, rising more steeply above the characteristic luminosity L_*. Redder galaxies exhibit a higher amplitude and steeper correlation function at all luminosities. The correlation amplitude of blue galaxies increases continuously with luminosity, but the luminosity dependence for red galaxies is less regular, with bright red galaxies more strongly clustered at large scales and faint red galaxies more strongly clustered at small scales. We interpret these results using halo occupation distribution (HOD) models assuming concordance cosmological parameters. For most samples, an HOD model with two adjustable parameters fits the w(r_p) data better than a power-law, explaining inflections at r_p ~ 1-3 Mpc/h as the transition between the 1-halo and 2-halo regimes of \xi(r). The implied minimum mass for a halo hosting a central galaxy above a luminosity threshold L grows as M_min ~ L at low luminosities and more steeply above L_*. The mass at which an average halo has one satellite galaxy brighter than L is M_1 ~ 23 M_min(L). These results imply a conditional luminosity function (at fixed halo mass) in which central galaxies lie far above a Schechter function extrapolation of the satellite population. HOD models nicely explain the joint luminosity-color dependence of w(r_p) in terms of the color fractions of central and satellite populations as a function of halo mass. The inferred HOD properties are in good qualitative agreement with theoretical predictions.

We discuss the construction of a photometric redshift catalogue of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS), emphasizing the principal steps necessary for constructing such a catalogue -- (i) photometrically selecting the sample, (ii) measuring photometric redshifts and their error distributions, (iii) and estimating the true redshift distribution. We compare two photometric redshift algorithms for these data and find that they give comparable results. Calibrating against the SDSS and SDSS-2dF spectroscopic surveys, we find that the photometric redshift accuracy is $\sigma \sim 0.03$ for redshifts less than 0.55 and worsens at higher redshift ($\sim 0.06$). These errors are caused by photometric scatter, as well as systematic errors in the templates, filter curves, and photometric zeropoints. We also parametrize the photometric redshift error distribution with a sum of Gaussians, and use this model to deconvolve the errors from the measured photometric redshift distribution to estimate the true redshift distribution. We pay special attention to the stability of this deconvolution, regularizing the method with a prior on the smoothness of the true redshift distribution. The methods we develop are applicable to general photometric redshift surveys.

We report the first result of the clustering analysis of Sloan Digital Sky Survey (SDSS) quasars. We compute the two-point correlation function (2PCF) of SDSS quasars in redshift space at $8h^{-1}{\rm Mpc} < s < 500h^{-1}{\rm Mpc}$, with particular attention to its baryonic signature. Our sample consists of 19986 quasars extracted from the SDSS Data Release 4 (DR4). The redshift range of the sample is $0.72 \le z \le 2.24$ (the mean redshift is $\bar z = 1.46$) and the reddening-corrected $i$-band apparent magnitude range is $15.0 \le m_{i,{\rm rc}} \le 19.1$. Due to the relatively low number density of the quasar sample, the bump in the power spectrum due to the baryon density, $\Omega_{\rm b}$, is not clearly visible. The effect of the baryon density is, however, to distort the overall shape of the 2PCF.The degree of distortion makes it an interesting alternate measure of the baryonic signature. Assuming a scale-independent linear bias and the spatially flat universe, i.e., $\Omega_{\rm b} + \Omega_{\rm d} + \Omega_\Lambda =1$, where $\Omega_{\rm d}$ and $\Omega_\Lambda$ denote the density parameters of dark matter and the cosmological constant, we combine the observed quasar 2PCF and the predicted matter 2PCF to put constraints on $\Omega_{\rm b}$ and $\Omega_\Lambda$. Our result is fitted as $0.80- 2.8\Omega_{\rm b} < \Omega_\Lambda < 0.90 - 1.4\Omega_{\rm b}$ at the 2$\sigma$ confidence level, which is consistent with results from other cosmological observations such as WMAP. (abridged)

We examine a sample of low redshift (10 < d < 150 Mpc) galaxies including galaxies with r-band absolute magnitudes as faint as -12.5 (for h=1), selected from the Sloan Digital Sky Survey Data Release 2 (SDSS). The sample is unique in containing galaxies of extremely low luminosities in a wide range of environments, selected with uniform and well-understood criteria. We present the luminosity function as well as the broad-band properties of low luminosity galaxies in this sample. A Schechter function is an insufficient parameterization of the r-band luminosity function; there is an upturn in the slope at low luminosities. The resulting slope at low luminosities in this sample is around -1.3. However, we almost certainly miss a large number of galaxies at very low luminosities due to low surface brightness selection effects, and we estimate that the true low luminosity slope may be as steep or steeper than -1.5. The results here are consistent with previous SDSS results and, in the g-band, roughly consistent with the results of the Two degree Field Galaxy Redshift Survey. Extremely low luminosity galaxies are predominantly blue, low surface brightness, exponential disks.

We present an 8 sigma detection of cosmic magnification measured by the variation of quasar density due to gravitational lensing by foreground large scale structure. To make this measurement we used 3800 square degrees of photometric observations from the Sloan Digital Sky Survey (SDSS) containing \~200,000 quasars and 13 million galaxies. Our measurement of the galaxy-quasar cross-correlation function exhibits the amplitude, angular dependence and change in sign as a function of the slope of the observed quasar number counts that is expected from magnification bias due to weak gravitational lensing. We show that observational uncertainties (stellar contamination, Galactic dust extinction, seeing variations and errors in the photometric redshifts) are well controlled and do not significantly affect the lensing signal. By weighting the quasars with the number count slope, we combine the cross-correlation of quasars for our full magnitude range and detect the lensing signal at >4 sigma in all five SDSS filters. Our measurements of cosmic magnification probe scales ranging from 60 kpc/h to 10 Mpc/h and are in good agreement with theoretical predictions based on the WMAP concordance cosmology. As with galaxy-galaxy lensing, future measurements of cosmic magnification will provide useful constraints on the galaxy-mass power spectrum.

We have used the 2dF instrument on the AAT to obtain redshifts of a sample of z<3, 18.0<g<21.85 quasars selected from SDSS imaging. These data are part of a larger joint programme: the 2dF-SDSS LRG and QSO Survey (2SLAQ). We describe the quasar selection algorithm and present the resulting luminosity function of 5645 quasars in 105.7 deg^2. The bright end number counts and luminosity function agree well with determinations from the 2dF QSO Redshift Survey (2QZ) data to g\sim20.2. However, at the faint end the 2SLAQ number counts and luminosity function are steeper than the final 2QZ results from Croom et al. (2004), but are consistent with the preliminary 2QZ results from Boyle et al. (2000). Using the functional form adopted for the 2QZ analysis, we find a faint end slope of beta=-1.78+/-0.03 if we allow all of the parameters to vary and beta=-1.45+/-0.03 if we allow only the faint end slope and normalization to vary. Our maximum likelihood fit to the data yields 32% more quasars than the final 2QZ parameterization, but is not inconsistent with other g>21 deep surveys. The 2SLAQ data exhibit no well defined ``break'' but do clearly flatten with increasing magnitude. The shape of the quasar luminosity function derived from 2SLAQ is in good agreement with that derived from type I quasars found in hard X-ray surveys. [Abridged]

We present the third edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog consists of the 46,420 objects in the SDSS Third Data Release that have luminosities larger than M_i = -22 (in a cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega_Lambda = 0.7), have at least one emission line with FWHM larger than 1000 km/s or are unambiguously broad absorption line quasars, are fainter than i = 15.0, and have highly reliable redshifts. The area covered by the catalog is 4188 sq. deg. The quasar redshifts range from 0.08 to 5.41, with a median value of 1.47; the high-redshift sample includes 520 quasars at redshifts greater than four, of which 17 are at redshifts greater than five. For each object the catalog presents positions accurate to better than 0.2 arcsec. rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800--9200A at a spectral resolution about 2000; the spectra can be retrieved from the public database using the information provided in the catalog. A total of 44,221 objects in the catalog were discovered by the SDSS; 28,400 of the SDSS discoveries are reported here for the first time.

We construct and analyze a u-band selected galaxy sample from the SDSS Southern Survey, which covers 275 sq. deg. The sample includes 43223 galaxies with spectroscopic redshifts in the range 0.005<z<0.3 and with 14.5<u<20.5. The S/N in the u-band Petrosian aperture is improved by coadding multiple epochs of imaging data and by including sky-subtraction corrections. Luminosity functions for the near-UV ^{0.1}u band (322+-26 nm) are determined in redshift slices of width 0.02, which show a highly significant evolution in M^{star} of -0.8+-0.1 mag between z=0 and 0.3; with M^{star} = -18.84+-0.05 (AB mag), log phi^{star} = -2.06+-0.03 (Mpc^{-3}) and log rho_L = 19.11+-0.02 (W Hz^{-1} Mpc^{-3}) at z=0.1. The faint-end slope determined for z<0.06 is given by alpha=-1.05+-0.08. This is in agreement with recent determinations from GALEX at shorter wavelengths. Comparing our z<0.3 luminosity density measurements with 0.2<z<1.2 from COMBO-17, we find that the 280-nm density evolves as rho_L proportional to (1+z)^{beta} with beta=2.1+-0.2; and find no evidence for any change in slope over this redshift range. By comparing with other measurements of cosmic star formation history, we estimate that the effective dust attenuation at 280 nm has increased by 0.8+-0.3 mag between z=0 and 1.

We report the intermediate-scale (0.3 to 40 Mpc/h) clustering of 35,000 luminous early-type galaxies at redshifts 0.16 to 0.44 from the Sloan Digital Sky Survey. We present the redshift-space two-point correlation function \xi(s), the projected correlation function w_p(r_p), and the deprojected real-space correlation function \xi(r), for approximately volume-limited samples. As expected, the galaxies are highly clustered, with the correlation length varying from 9.8 +/- 0.2 Mpc/h to 11.2 +/- 0.2 Mpc/h, dependent on the specific luminosity range. For the -23.2 < Mg < -21.2 sample, the inferred bias relative to that of L* galaxies is 1.84 +/- 0.11 for 1 Mpc/h < r_p < 10 Mpc/h, with yet stronger clustering on smaller scales. We detect luminosity-dependent bias within the sample but see no evidence for redshift evolution between z=0.2 and z=0.4. We find a clear indication for deviations from a power-law in the real-space correlation function, with a dip at ~ 2 Mpc/h scales and an upturn on smaller scales. The precision measurements of these clustering trends offer new avenues for the study of the formation and evolution of these massive galaxies.

We present the small-scale (0.2h^{-1} to 7h^{-1} Mpc) cross-correlations between 32,000 luminous early-type galaxies and a reference sample of 16 million normal galaxies from the Sloan Digital Sky Survey. Our method allows us to construct the spherically averaged, real-space cross-correlation function between the spectroscopic LRG sample and galaxies from the SDSS imaging. We report the cross-correlation as a function of scale, luminosity, and redshift. We find very strong luminosity dependences in the clustering amplitudes, up to a factor of 4 over a factor of 4 in luminosity, and measure this dependence with high signal-to-noise ratio. The luminosity dependence of bias is found to depend on scale, with more variation on smaller scales. The clustering as a function of scale is not a power law, but instead has a dip at 1h^{-1} Mpc and an excess on small scales. The fraction of red galaxies within the L* sample surrounding LRGs is a strong function of scale, as expected. However, the fraction of red galaxies evolves in redshift similarly on small and large scales, suggesting that cluster and field populations are changing in the same manner. The results highlight the advantage on small scales of using cross-correlation methods as a means of avoiding shot noise in samples of rare galaxies.

Compact groups (CGs) of galaxies -- relatively poor groups of galaxies in which the typical separations between members is of the order of a galaxy diameter -- offer an exceptional laboratory for the study of dense galaxian environments with short (<1Gyr) dynamical time-scales. In this paper, we present an objectively defined catalog of CGs in 153 sq deg of the Sloan Digital Sky Survey Early Data Release (SDSS EDR). To identify CGs, we applied a modified version of Hickson's (1982) criteria aimed at finding the highest density CGs and thus reducing the number of chance alignments. Our catalog contains 175 CGs down to a limiting galaxy magnitude of r* = 21. The resulting catalog has a median depth of approximately z = 0.13, substantially deeper than previous CG catalogs. Since the SDSS will eventually image up to one quarter of the celestial sphere, we expect our final catalog, based upon the completed SDSS, will contain on the order of 5,000 - 10,000 CGs. This catalog will be useful for conducting studies of the general characteristics of CGs, their environments, and their component galaxies.

We measure the redshift-dependent luminosity function and the comoving radial density of galaxies in the Sloan Digital Sky Survey Data Release 1 (SDSS DR1). Both measurements indicate that the apparent number density of bright galaxies increases by a factor ~3 as redshift increases from z=0 to z=0.3. This result is robust to the assumed cosmology, to the details of the K-correction and to direction on the sky. These observations are most naturally explained by significant evolution in the luminosity and/or number density of galaxies at redshifts z < 0.3. Such evolution is also consistent with the steep number-magnitude counts seen in the APM Galaxy Survey, without the need to invoke a local underdensity in the galaxy distribution distribution or magnitude scale errors.

We present CCD photometry in the Johnson-Kron-Cousins UBR bands for a sample of 320 galaxies selected from the Stromlo-APM Redshift Survey. We use this CCD data to estimate the galaxy luminosity function in the U, B, R and b_J bands, finding consistent results with earlier work. Fainter galaxies serendipitously observed on our CCD frames allow a check of the photometric calibration of the APM Galaxy Survey. We find no evidence of any significant scale error in the APM magnitudes.

In previous work on galaxy clusters, several authors reported a discovery of an unusual population of galaxies, which have spiral morphologies, but do not show any star formation activity. These galaxies are called ``passive spirals'', and have been interesting since it has been difficult to understand the existence of such galaxies. Using a volume limited sample (0.05<z<0.1 and Mr<-20.5; 25813 galaxies) of the Sloan Digital Sky Survey data, we have found 73 (0.28$\pm$0.03%) passive spiral galaxies and studied their environments. It is found that passive spiral galaxies live in local galaxy density 1-2 Mpc$^{-2}$ and 1-10 cluster-centric virial radius. Thus the origins of passive spiral galaxies are likely to be cluster related. These characteristic environments coincide with the previously reported environment where galaxy star formation rate suddenly declines and the so-called morphology-density relation turns. It is likely that the same physical mechanism is responsible for all of these observational results. The existence of passive spiral galaxies suggests that a physical mechanism that works calmly is preferred to dynamical origins such as major merger/interaction since such a mechanism can destroy spiral arm structures. Compared with observed cluster galaxy evolution such as the Butcher-Oemler effect and the morphological Butcher-Oemler effect, passive spiral galaxies are likely to be a galaxy population in transition between red, elliptical/S0 galaxies in low redshift clusters and blue, spiral galaxies more numerous in higher redshift clusters.

Supervised artificial neural networks are used to predict useful properties of galaxies in the Sloan Digital Sky Survey, in this instance morphological classifications, spectral types and redshifts. By giving the trained networks unseen data, it is found that correlations between predicted and actual properties are around 0.9 with rms errors of order ten per cent. Thus, given a representative training set, these properties may be reliably estimated for galaxies in the survey for which there are no spectra and without human intervention.

Understanding the clustering of galaxies has long been a goal of modern observational cosmology. Utilizing our photometric redshift technique a volume limited sample containing more than 2 million galaxies is constructed from the SDSS galaxy catalog. In the largest such analysis to date, we study the angular clustering as a function of luminosity and spectral type. Using Limber's equation we calculate the clustering length for the full data set as r0=5.77+/-0.10 Mpc/h. We find that r0 increases with luminosity by a factor of 1.6 over the sampled luminosity range, in agreement with previous redshift surveys. We also find that both the clustering length and the slope of the correlation function depend on the galaxy type. In particular, by splitting the galaxies in four groups by their rest-frame type we find a bimodal behavior in their clustering properties. Galaxies with spectral types similar to elliptical galaxies have a correlation length of 6.59 +/- 0.17 Mpc/h and a slope of the angular correlation function of 0.96 +/- 0.05 while blue galaxies have a clustering length of 4.51 +/- 0.19 Mpc/h and a slope of 0.68 +/- 0.09. The two intermediate color groups behave like their more extreme 'siblings', rather than showing a gradual transition in slope. We discuss these correlations in the context of current cosmological models for structure formation.

[Abridged] We present here a new and homogeneous sample of 3340 galaxies selected from the Sloan Digital Sky Survey (SDSS) based solely on the observed strength of their Hdelta absorption line. These galaxies are commonly known as ``post-starburst'' or ``E+A'' galaxies, and the study of these galaxies has been severely hampered by the lack of a large, statistical sample of such galaxies. In this paper, we rectify this problem by selecting a sample of galaxies which possess an absorption Hdelta equivalent width of EW(Hdelta_max) - Delta EW(Hdelta_max) > 4A from 106682 galaxies in the SDSS. We have performed extensive tests on our catalog including comparing different methodologies of measuring the Hdelta absorption and studying the effects of stellar absorption, dust extinction, emission-filling and measurement error. The measured abundance of our Hdelta-selected (HDS) galaxies is 2.6 +/- 0.1% of all galaxies within a volume-limited sample of 0.05<z<0.1 and M(r*)<-20.5, which is consistent with previous studies of such galaxies in the literature. We find that only 25 of our HDS galaxies in this volume-limited sample (3.5+/-0.7%) show no evidence for OII and Halpha emission, thus indicating that true E+A (or k+a) galaxies are extremely rare objects at low redshift, i.e., only 0.09+/-0.02% of all galaxies in this volume-limited sample are true E+A galaxies. In contrast, 89+/-5% of our HDS galaxies in the volume-limited sample have significant detections of the OII and Halpha emission lines. We find 27 galaxies in our volume-limited HDS sample that possess no detectable OII emission, but do however possess detectable Halpha emission. These galaxies may be dusty star-forming galaxies. We provide the community with this new catalog of Hdelta-selected galaxies to aid in the understanding of these galaxies.

We present a determination of the `Cosmic Optical Spectrum' of the Universe, i.e. the ensemble emission from galaxies, as determined from the red-selected Sloan Digital Sky Survey main galaxy sample and compare with previous results of the blue-selected 2dF Galaxy Redshift Survey. Broadly we find good agreement in both the spectrum and the derived star-formation histories. If we use a power-law star-formation history model where star-formation rate $\propto (1+z)^\beta$ out to z=1, then we find that $\beta$ of 2 to 3 is still the most likely model and there is no evidence for current surveys missing large amounts of star formation at high redshift. In particular `Fossil Cosmology' of the local universe gives measures of star-formation history which are consistent with direct observations at high redshift. Using the photometry of SDSS we are able to derive the cosmic spectrum in absolute units (i.e.$ W \AA$^{-1}$ Mpc$^{-3}$) at 2--5\AA resolution and find good agreement with published broad-band luminosity densities. For a Salpeter IMF the best fit stellar mass/light ratio is 3.7--7.5 $\Msun/\Lsun$ in the r-band (corresponding to $\omstars h = 0.0025$--0.0055) and from both the stellar emission history and the H$\alpha$ luminosity density independently we find a cosmological star-formation rate of 0.03--0.04 h $\Msun$ yr$^{-1}$ Mpc$^{-3}$ today.

The morphological dependence of the luminosity function is studied using a sample containing approximately 1500 bright galaxies classified into Hubble types by visual inspections for a homogeneous sample obtained from the Sloan Digital Sky Survey (SDSS) northern equatorial stripes. Early-type galaxies are shown to have a characteristic magnitude by 0.45 mag brighter than spiral galaxies in the $r^{\ast}$ band, consistent with the `universal characteristic luminosity' in the $B$ band. The shape of the luminosity function differs rather little among different morphological types: we do not see any symptoms of the sharp decline in the faint end for the luminosity function for early-type galaxies at least 2 mag fainter than the characteristic magnitude, although the faint end behaviour shows a slight decline ($\alpha\lsim -1$) compared with the total sample. We also show that a rather flat faint end slope for early-type galaxies is not due to an increasing mixture of the dwarf galaxies which have softer cores. This means that there are numerous faint early-type galaxies with highly concentrated cores.

We combine Sloan Digital Sky Survey spectra of 22,000 luminous, red, bulge-dominated galaxies to get high S/N average spectra in the rest-frame optical and ultraviolet (2600A to 7000A). The average spectra of these massive, quiescent galaxies are early-type with weak emission lines and with absorption lines indicating an apparent excess of alpha elements over solar abundance ratios. We make average spectra of subsamples selected by luminosity, environment and redshift. The average spectra are remarkable in their similarity. What variations do exist in the average spectra as a function of luminosity and environment are found to form a nearly one-parameter family in spectrum space. We present a high signal-to-noise ratio spectrum of the variation. We measure the properties of the variation with a modified version of the Lick index system and compare to model spectra from stellar population syntheses. The variation may be a combination of age and chemical abundance differences, but the conservative conclusion is that the quality of the data considerably exceeds the current state of the models.

Using a catalog of 147,986 galaxy redshifts and fluxes from the Sloan Digital Sky Survey (SDSS) we measure the galaxy luminosity density at z=0.1 in five optical bandpasses corresponding to the SDSS bandpasses shifted to match their restframe shape at z=0.1. We denote the bands {0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}, with \lambda_{eff} = [3216, 4240, 5595, 6792, 8111] Angstroms respectively. We use a maximum likelihood method which allows for a general form for the shape of the luminosity function, simple luminosity and number evolution, incorporates flux uncertainties, and accounts for the flux limits of the survey. We find luminosity densities at z=0.1 in absolute AB magnitudes in a Mpc^3 of [-14.10 \pm 0.15, -15.18 \pm 0.03, -15.90 \pm 0.03, -16.24 \pm 0.03, -16.56 \pm 0.02] in [{0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}], respectively, using \Omega_0 =0.3, \Omega_\Lambda=0.7, and h=1, and using Petrosian magnitudes. Similar results are obtained using Sersic model magnitudes, suggesting that flux from outside the Petrosian apertures is not a major correction. In the {0.1}{r} band, the best fit Schechter function to our results has \phi_\ast = (1.49 \pm 0.04) \times 10^{-2} h^3 Mpc^{-3}, M_\ast - 5\log_{10} h = -20.44 \pm 0.01, and \alpha = -1.05\pm 0.01. In solar luminosities, the luminosity density in {0.1}{r} is (1.84 \pm 0.04) h 10^8 L_{{0.1}{r},\odot} Mpc^{-3}. Our results are consistent with other estimates of the luminosity density, from the Two-degree Field Galaxy Redshift Survey and the Millenium Galaxy Catalog. They represent a substantial change (\sim 0.5 mag) from earlier SDSS luminosity density results based on commissioning data, almost entirely because of the inclusion of evolution in the luminosity function model.

Using photometry and spectroscopy of 144,609 galaxies from the Sloan Digital Sky Survey, we present bivariate distributions of pairs of seven galaxy properties: four optical colors, surface brightness, radial profile shape as measured by the Sersic index, and absolute magnitude. In addition, we present the dependence of local galaxy density (smoothed on 8 h^{-1} Mpc scales) on all of these properties. Several classic, well-known relations among galaxy properties are evident at extremely high signal-to-noise ratio: the color-color relations of galaxies, the color-magnitude relations, the magnitude-surface brightness relation, and the dependence of density on color and absolute magnitude. We show that most of the i-band luminosity density in the universe is in the absolute magnitude and surface brightness ranges used. Some of the relationships between parameters, in particular the color--magnitude relations, show stronger correlations for exponential galaxies and concentrated galaxies taken separately than for all galaxies taken together. We provide a simple set of fits of the dependence of galaxy properties on luminosity for these two sets of galaxies.

The Sloan Digital Sky Survey (SDSS) is making a multi-colour, three dimensional map of the nearby Universe. The survey is in two parts. The first part is imaging one quarter of the sky in five colours from the near ultraviolet to the near infrared. In this imaging survey we expect to detect around 50 million galaxies to a magnitude limit g ~ 23. The second part of the survey, taking place concurrently with the imaging, is obtaining spectra for up to 1 million galaxies and 100,000 quasars. From these spectra we obtain redshifts and hence distances, in order to map out the three-dimensional distribution of galaxies and quasars in the Universe. These observations will be used to constrain models of cosmology and of galaxy formation and evolution. This article describes the goals and methods used by the SDSS, the current status of the survey, and highlights some exciting discoveries made from data obtained in the first two years of survey operations.

The angular distribution of galaxies encodes a wealth of information about large scale structure. Ultimately, the Sloan Digital Sky Survey (SDSS) will record the angular positions of order 10^8 galaxies in five bands, adding significantly to the cosmological constraints. This is the first in a series of papers analyzing a rectangular stripe 2.5x90 degrees from early SDSS data. We present the angular correlation function for galaxies in four separate magnitude bins on angular scales ranging from 0.003 degrees to 15 degrees. Much of the focus of this paper is on potential systematic effects. We show that the final galaxy catalog -- with the mask accounting for regions of poor seeing, reddening, bright stars, etc. -- is free from external and internal systematic effects for galaxies brighter than r* = 22. Our estimator of the angular correlation function includes the effects of the integral constraint and the mask. The full covariance matrix of errors in these estimates is derived using mock catalogs with further estimates using a number of other methods.

We describe the algorithm for selecting quasar candidates for optical spectroscopy in the Sloan Digital Sky Survey. Quasar candidates are selected via their non-stellar colors in "ugriz" broad-band photometry, and by matching unresolved sources to the FIRST radio catalogs. The automated algorithm is sensitive to quasars at all redshifts lower than z=5.8. Extended sources are also targeted as low-redshift quasar candidates in order to investigate the evolution of Active Galactic Nuclei (AGN) at the faint end of the luminosity function. Nearly 95% of previously known quasars are recovered (based on 1540 quasars in 446 square degrees). The overall completeness, estimated from simulated quasars, is expected to be over 90%, whereas the overall efficiency (quasars:quasar candidates) is better than 65%. The selection algorithm targets ultraviolet excess quasars to i^*=19.1 and higher-redshift (z>3) quasars to i^*=20.2, yielding approximately 18 candidates per square degree. In addition to selecting ``normal'' quasars, the design of the algorithm makes it sensitive to atypical AGN such as Broad Absorption Line quasars and heavily reddened quasars.

We compute the angular power spectrum C_l from 1.5 million galaxies in early SDSS data on large angular scales, l<600. The data set covers about 160 square degrees, with a characteristic depth of order 1 Gpc/h in the faintest (21<r<22) of our four magnitude bins. Cosmological interpretations of these results are presented in a companion paper by Dodelson et al (2001). The data in all four magnitude bins are consistent with a simple flat ``concordance'' model with nonlinear evolution and linear bias factors of order unity. Nonlinear evolution is particularly evident for the brightest galaxies. A series of tests suggest that systematic errors related to seeing, reddening, etc., are negligible, which bodes well for the sixtyfold larger sample that the SDSS is currently collecting. Uncorrelated error bars and well-behaved window functions make our measurements a convenient starting point for cosmological model fitting.

We present initial results for counts in cells statistics of the angular distribution of galaxies in early data from the Sloan Digital Sky Survey (SDSS). We analyze a rectangular stripe $2.5^\circ$ wide, covering approximately 160 sq. degrees, containing over $10^6$ galaxies in the apparent magnitude range $18 < r^\prime < 22$, with areas of bad seeing, contamination from bright stars, ghosts, and high galactic extinction masked out. This survey region, which forms part of the SDSS Early Data Release, is the same as that for which two-point angular clustering statistics have recently been computed. The third and fourth moments of the cell counts, $s_3$ (skewness) and $s_4$ (kurtosis), constitute the most accurate measurements to date of these quantities (for $r^\prime < 21$) over angular scales $0.015^\circ-0.3^\circ$. They display the approximate hierarchical scaling expected from non-linear structure formation models and are in reasonable agreement with the predictions of $\Lambda$-dominated cold dark matter models with galaxy biasing that suppresses higher order correlations at small scales. The results are in general consistent with previous measurements in the APM, EDSGC, and Deeprange surveys. These results suggest that the SDSS imaging data are free of systematics to a high degree and will therefore enable determination of the skewness and kurtosis to 1% and less then 10%, as predicted by Colombi, Szapudi, & Szalay (1998).

We positionally match sources observed by the Sloan Digital Sky Survey (SDSS), the Two Micron All Sky Survey (2MASS), and the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. Practically all 2MASS sources are matched to an SDSS source within 2 arcsec; ~11% of them are optically resolved galaxies and the rest are dominated by stars. About 1/3 of FIRST sources are matched to an SDSS source within 2 arcsec; ~80% of these are galaxies and the rest are dominated by quasars. Based on these results, we project that by the completion of these surveys the matched samples will include about 10^7 stars and 10^6 galaxies observed by both SDSS and 2MASS, and about 250,000 galaxies and 50,000 quasars observed by both SDSS and FIRST. Here we present a preliminary analysis of the optical, infrared and radio properties for the extragalactic sources from the matched samples. In particular, we find that the fraction of quasars with stellar colors missed by the SDSS spectroscopic survey is probably not larger than ~10%, and that the optical colors of radio-loud quasars are ~0.05 mag. redder (with 4-sigma significance) than the colors of radio-quiet quasars.

I summarize the scientific goals and current status of the Sloan Digital Sky Survey, briefly describe the Early Data Release, and discuss some recent scientific results obtained from commissioning data which are apposite to the distribution of luminous and dark matter in the Universe.

We present new measurements of scaling laws relating the luminosity of galaxies to the amplitude and shape of their dark matter halos. Early imaging and spectroscopic data from the Sloan Digital Sky Survey are used to make weak lensing measurements of the surface mass density contrast Delta Sigma_+ around classes of lens objects. This surface mass density contrast as a function of radius is a measure of the galaxy-mass correlation function (GMCF). Because spectroscopic redshifts are available for all lens objects, the mass and distance scales are well constrained. The GMCF measured around ~31,000 lenses is well fit by a power law of the form Delta Sigma_+ = (2.5+0.7-0.6) (R/1 Mpc)^{-0.8+-0.2} h M_sun pc^-2. We compare this GMCF to galaxy luminosity, type, and environment, and find that it varies strongly with all three. We quantify these variations by comparing the normalization of a fit to the inner 260 h^-1 kpc, M_260, to the galaxy luminosity. While M_260 is not strongly related to luminosity in bluest band u', there is a simple, linear relation between M_260 and luminosity in redder bands (g', r', i', and z'). We test the universality of these mass-to-light scalings by independently measuring them for spiral and elliptical galaxies,and for galaxies in a variety of environments. We find remarkable consistency in these determinations in the red bands, especially i' and z'. This consistency across a wide range of systems suggests that the measured scaling represents an excellent cosmic average, and that the integrated star formation history of galaxies is strongly related to the dark matter environments in which they form.

We present measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

The Sloan Digital Sky Survey is one of the first multicolor photometric and spectroscopic surveys designed to measure the statistical properties of galaxies within the local Universe. In this Letter we present some of the initial results on the angular 2-point correlation function measured from the early SDSS galaxy data. The form of the correlation function, over the magnitude interval 18<r*<22, is shown to be consistent with results from existing wide-field, photographic-based surveys and narrower CCD galaxy surveys. On scales between 1 arcminute and 1 degree the correlation function is well described by a power-law with an exponent of ~ -0.7. The amplitude of the correlation function, within this angular interval, decreases with fainter magnitudes in good agreement with analyses from existing galaxy surveys. There is a characteristic break in the correlation function on scales of approximately 1-2 degrees. On small scales, < 1', the SDSS correlation function does not appear to be consistent with the power-law form fitted to the 1'< theta <0.5 deg data. With a data set that is less than 2% of the full SDSS survey area, we have obtained high precision measurements of the power-law angular correlation function on angular scales 1' < theta < 1 deg, which are robust to systematic uncertainties. Because of the limited area and the highly correlated nature of the error covariance matrix, these initial results do not yet provide a definitive characterization of departures from the power-law form at smaller and larger angles. In the near future, however, the area of the SDSS imaging survey will be sufficient to allow detailed analysis of the small and large scale regimes, measurements of higher-order correlations, and studies of angular clustering as a function of redshift and galaxy type.

Early photometric data from the Sloan Digital Sky Survey (SDSS) contain angular positions for 1.5 million galaxies. In companion papers, the angular correlation function $w(\theta)$ and 2D power spectrum $C_l$ of these galaxies are presented. Here we invert Limber's equation to extract the 3D power spectrum from the angular results. We accomplish this using an estimate of $dn/dz$, the redshift distribution of galaxies in four different magnitude slices in the SDSS photometric catalog. The resulting 3D power spectrum estimates from $w(\theta)$ and $C_l$ agree with each other and with previous estimates over a range in wavenumbers $0.03 < k/{\rm h Mpc}^{-1} < 1$. The galaxies in the faintest magnitude bin ($21 < \rstar < 22$, which have median redshift $z_m=0.43$) are less clustered than the galaxies in the brightest magnitude bin ($18 < \rstar < 19$ with $z_m=0.17$), especially on scales where nonlinearities are important. The derived power spectrum agrees with that of Szalay et al. (2001) who go directly from the raw data to a parametric estimate of the power spectrum. The strongest constraints on the shape parameter $\Gamma$ come from the faintest galaxies (in the magnitude bin $21 < \rstar < 22$), from which we infer $\Gamma = 0.14^{+0.11}_{-0.06}$ (95% C.L.).

The Sloan Digital Sky Survey (SDSS) will observe around 10^6 spectra from targets distributed over an area of about 10,000 square degrees, using a multi-object fiber spectrograph which can simultaneously observe 640 objects in a circular field-of-view (referred to as a ``tile'') 1.49 degrees in radius. No two fibers can be placed closer than 55'' during the same observation; multiple targets closer than this distance are said to ``collide.'' We present here a method of allocating fibers to desired targets given a set of tile centers which includes the effects of collisions and which is nearly optimally efficient and uniform. Because of large-scale structure in the galaxy distribution (which form the bulk of the SDSS targets), a naive covering the sky with equally-spaced tiles does not yield uniform sampling. Thus, we present a heuristic for perturbing the centers of the tiles from the equally-spaced distribution which provides more uniform completeness. For the SDSS sample, we can attain a sampling rate greater than 92% for all targets, and greater than 99% for the set of targets which do not collide with each other, with an efficiency greater than 90% (defined as the fraction of available fibers assigned to targets).

We present bright galaxy number counts in five broad bands ($u', g', r', i', z'$) from imaging data taken during the commissioning phase of the Sloan Digital Sky Survey (SDSS). The counts are derived from two independent stripes of imaging scans along the Celestial Equator, one each toward the North and the South Galactic cap, covering about 230 and 210 square degrees, respectively. A careful study is made to verify the reliability of the photometric catalog. For galaxies brighter than $r^* = 16$, the catalog produced by automated software is examined against eye inspection of all objects. Statistically meaningful results on the galaxy counts are obtained in the magnitude range $12 \le r^* \le 21$, using a sample of 900,000 galaxies. The counts from the two stripes differ by about 30% at magnitudes brighter than $r^*= 15.5$, consistent with a local $2\sigma$ fluctuation due to large scale structure in the galaxy distribution. The shape of the number counts-magnitude relation brighter than $r^* = 16$ is well characterized by $N \propto 10^{0.6m}$, the relation expected for a homogeneous galaxy distribution in a ``Euclidean'' universe. In the magnitude range $ 16 < r^* < 21$, the galaxy counts from both stripes agree very well, and follow the prediction of the no-evolution model, although the data do not exclude a small amount of evolution. We use empirically determined color transformations to derive the galaxy number counts in the $B$ and $I_{814}$ bands. We compute the luminosity density of the universe at zero redshift in the five SDSS bands and in the $B$ band. We find ${\cal L}_{B} = 2.4 \pm 0.4 \times 10^8L_\odot h $Mpc$^{-3}$, for a reasonably wide range of parameters of the Schechter luminosity function in the $B$ band.

We present a measurement of the K-band luminosity function (LF) of field galaxies obtained from near-infrared imaging of a sample of 345 galaxies selected from the Stromlo-APM Redshift Survey. The LF is well-fitted over the ten magnitude range -26 < M_K < -16 by a Schechter function with parameters alpha = -1.16 +- 0.19, M* = -23.58 +- 0.42, phi* = 0.012 +- 0.008 Mpc^-3 assuming a Hubble constant of H_0 = 100 km/s/Mpc. We have also estimated the LF for two subsets of galaxies subdivided by the equivalent width of the Halpha emission line at EW(Halpha) = 10A. There is no significant difference in LF shape between the two samples, although there is a hint (~1 sigma significance) that emission line galaxies (ELGs) have M* roughly one magnitude fainter than non-ELGs. Contrary to the optical LF, there is no difference in faint-end slope alpha between the two samples.

We study the luminosity function and clustering properties of subsamples of local galaxies selected from the Stromlo-APM survey by the rest-frame equivalent width (EW) of the Halpha emission line. The b_J luminosity function of star-forming galaxies has a significantly steeper faint-end slope than that for quiescent galaxies: the majority of sub-L* galaxies are currently undergoing significant star formation. Emission line galaxies are less strongly clustered, both amongst themselves, and with the general galaxy population, than quiescent galaxies. Thus as well as being less luminous, star-forming galaxies also inhabit lower-density regions of the Universe than quiescent galaxies.

We analyze spectral properties of 1671 galaxies from the Stromlo-APM survey, selected to have 15 < b_J < 17.15 and having a mean redshift z = 0.05. This is a representative local sample of field galaxies, so the global properties of the galaxy population provide a comparative point for analysis of more distant surveys. We measure Halpha, Oii 3727, Sii 6716, 6731, Nii 6583 and Oi 6300 equivalent widths and the D_4000 break index. The 5A resolution spectra use an 8 arcsec slit, which typically covers 40-50% of the galaxy area. We find no evidence for systematic trends depending on the fraction of galaxy covered by the slit, and further analysis suggests that our spectra are representative of integrated galaxy spectra. We classify spectra according to their Halpha emission, which is closely related to massive star formation. Overall we find 61% of galaxies are Halpha emitters with rest-frame equivalent widths EW(Halpha) >= 2A. The emission-line galaxy (ELG) fraction is smaller than seen in the CFRS at z = 0.2 and is consistent with a rapid evolution of Halpha luminosity density. The ELG fraction, and EW(Halpha), increase at fainter absolute magnitudes, smaller projected area and smaller D_4000. In the local Universe, faint, small galaxies are dominated by star formation activity, while bright, large galaxies are more quiescent. This picture of the local Universe is quite different from the distant one, where bright galaxies appear to show rapidly-increasing activity back in time. (Abridged)

We study the luminosity function and clustering properties of subsamples of local galaxies selected from the Stromlo-APM survey by the rest-frame equivalent widths of their Halpha and Oii emission lines. The b_J luminosity function of star-forming galaxies has a significantly steeper faint-end slope than that for quiescent galaxies: the majority of sub-L* galaxies are currently undergoing significant star formation. Emission line galaxies are less strongly clustered, both amongst themselves, and with the general galaxy population, than quiescent galaxies. Thus as well as being less luminous, star-forming galaxies also inhabit lower-density regions of the Universe than quiescent galaxies.

We present preliminary results of a neutral hydrogen (HI) redshift survey to find low surface brightness (LSB) galaxies in the very nearby universe. Our sample consists of all galaxies in the APM catalog (Maddox et al 1990) with a mean surface brightness of mu > 24 mag/arcsec^2, down to a magnitude limit of b_j < 17. With the Parkes 64m Radiotelescope 35 objects were detected at v < 4300 km/s. The resulting luminosity function, HI mass function, and for the first time field mass function are presented. It is found that LSBs make a negligible contribution to the overall integrated luminosity, HI mass, and total mass contained in galaxies.

We present preliminary measurements of the b_J and K-band luminosity functions (LFs) of field galaxies obtained from optical and K-band imaging of a sample of galaxies selected from the Stromlo-APM Redshift Survey. The b_J LF is consistent with that previously published from photographic data. The K-band LF has been estimated over a range of 12 magnitudes and is reasonably well fit by a Schechter function with faint-end slope alpha = -1.2.

The Sloan Digital Sky Survey (SDSS) will carry out a digital photometric and spectroscopic survey over pi steradians in the northern Galactic cap. An array of CCD detectors used in drift-scan mode will image the sky in five passbands to a limiting magnitude of r' ~ 23. Selected from the imaging survey, 10^6 galaxies, 10^5 quasars and selected samples of stars will be observed spectroscopically. I describe the current status of the survey, which recently saw first light, and its prospects for constraining models for dark matter in the Universe.

The Sloan Digital Sky Survey (SDSS) will carry out a digital photometric and spectroscopic survey over pi steradians in the northern Galactic cap. An array of CCD detectors used in drift-scan mode will image the sky in five passbands to a limiting magnitude of r' ~ 23. Selected from the imaging survey, 10^6 galaxies, 10^5 quasars and selected samples of stars will be observed spectroscopically. In addition, a smaller (225 deg^2), deeper, southern survey will reach ~2.0 magnitudes fainter and will contain a wealth of information about variable sources, supernovae and proper motions. We describe the current status of the survey, which recently saw first light, and its prospects for constraining cosmological models.

I review the current observational status of the faint end of the optical luminosity function of field galaxies at low redshift. There is growing evidence for an excess number of dwarf galaxies that is not well fit by a single Schechter function. These dwarf galaxies tend to be of late morphological and spectral type, blue in colour, of low surface brightness and currently undergoing significant star formation.

We analyse emission-line properties of the bj-selected Stromlo-APM spectra (<z> = 0.05). Because this is a representative sample, we can study the global spectral properties of the local galaxy population. We classify spectra according to their H_alpha emission, which is closely related to massive star formation. This study gives a comparative local point for analysis of more distant surveys. We show that in the local universe, faint, small galaxies are dominated by star formation activity, while bright, large galaxies are more quiescent. Obviously this picture of the local universe is quite different from the distant one, where bright galaxies appear to show a rapidly-increasing activity back in time.

We estimate the luminosity function of field galaxies over a range of ten magnitudes (-22 < M_{B_J} < -12 for H_0 = 100 km/s/Mpc) by counting the number of faint APM galaxies around Stromlo-APM redshift survey galaxies at known distance. The faint end of the luminosity function rises steeply at M_{B_J} \approx -15, implying that the space density of dwarf galaxies is at least two times larger than predicted by a Schechter function with flat faint-end slope. Such a high abundance of dwarf galaxies at low redshift can help explain the observed number counts and redshift distributions of faint galaxies without invoking exotic models for galaxy evolution.

The Sloan Digital Sky Survey (SDSS) is a project to definitively map $\pi$ steradians of the local Universe. An array of CCD detectors used in drift-scan mode will digitally image the sky in five passbands to a limiting magnitude of $r' \sim 23$. Selected from the imaging survey, $10^6$ galaxies and $10^5$ quasars will be observed spectroscopically. I describe the current status of the survey, which is due to begin observations early in 1997, and its prospects for constraining models for dark matter in the Universe.

The APM Bright Galaxy Catalogue lists positions, magnitudes, shapes and morphological types for 14,681 galaxies brighter than $b_J$ magnitude 16.44 over a 4,180 square degree area of the southern sky. Galaxy and stellar images have been located from glass copy plates of the United Kingdom Schmidt Telescope (UKST) IIIaJ sky survey using the Automated Photographic Measuring (APM) facility in Cambridge, England. The majority of stellar images are rejected by the regularity of their image surface brightness profiles. Remaining images are inspected by eye on film copies of the survey material and classed as stellar, multiple stellar, galaxy, merger or noise. Galaxies are further classified as elliptical, lenticular, spiral, irregular or uncertain. The 180 survey fields are put onto a uniform photometric system by comparing the magnitudes of galaxies in the overlap regions between neighbouring plates. The magnitude zero-point, photometric uniformity and photographic saturation are checked with CCD photometry. Finally, the completeness and reliability of the catalogue is assessed using various internal tests and by comparing with several independently constructed galaxy catalogues.

Galaxy redshift surveys provide a distorted picture of the universe due to the non-Hubble component of galaxy motions. By measuring such distortions in the linear regime one can constrain the quantity $\beta = \Omega^{0.6}/b$ where $\Omega$ is the cosmological density parameter and $b$ is the (linear) bias factor for optically-selected galaxies. In this paper we apply two techniques for estimating $\beta$ from the Stromlo-APM redshift survey --- (1) measuring the anisotropy of the redshift space correlation function in spherical harmonics and (2) comparing the amplitude of the direction-averaged redshift space correlation function to the real space correlation function. We test the validity of these techniques, particularly whether the assumption of linear theory is justified, using two sets of large $N$-body simulations. We find that the first technique is affected by non-linearities on scales up to $\sim 30 \hMpc$. The second technique is less sensitive to non-linear effects and so is more useful for existing redshift surveys. The Stromlo-APM survey data favours a low value for $\beta$, with $\beta \simlt 0.6$. A bias parameter $b \approx 2$ is thus required if $\Omega \equiv 1$. However, higher-order correlations measured from the APM galaxy survey (Gazta\~{n}aga and Frieman 1994) indicate a low value for the bias parameter $b \approx 1$, requiring that $\Omega \simlt 0.5$. We also measure the relative bias for samples of galaxies of various luminosity and morphological type and find that low-luminosity galaxies are roughly three times less biased than $L^*$ galaxies. For the galaxy population as a whole, we measure a real space variance of galaxy counts in $8 \hMpc$ spheres of $(\sigma^2_8)_g = 0.89 \pm 0.05$.

We present clustering measurements for samples of galaxies selected by morphological type and luminosity from the recently completed Stromlo-APM Redshift Survey. Early type galaxies are clustered more strongly by a factor 3.5--5.5, than late type galaxies. Low-luminosity galaxies are clustered more weakly by a factor of $\sim 2$ than $L^*$ and brighter galaxies on scales $\simgt 1\hMpc$. Also the slope of the correlation function is steeper for low-luminosity galaxies, so that the amplitude is a factor 4 lower at $ 10 \hMpc$. No difference, however, is seen between the clustering of $L^*$ and more luminous galaxies, an observation which may be hard to reconcile with some theories of biased galaxy formation. Both redshift-space and real-space clustering estimates show a similar dependence on luminosity.