The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA).

Galaxy spectral energy distributions (SEDs) remain among the most challenging yet informative quantities to reproduce in simulations due to the large and complex mixture of physical processes that shape the radiation output of a galaxy. With the increasing number of surveys utilising broadband colours as part of their target selection criteria, the production of realistic SEDs in simulations is necessary for assisting in survey design and interpretation of observations. The recent success in reproducing the observed luminosity functions (LF) from far-UV to far-IR, using the state-of-the-art semi-analytic model \shark\ and the SED generator \prospect, represents a critical step towards better galaxy colour predictions. We show that with \shark\ and \prospect\ we can closely reproduce the optical colour distributions observed in the panchromatic GAMA survey. The treatment of feedback, star formation, central-satellite interactions and radiation re-processing by dust are critical for this achievement. The first three processes create a bimodal distribution, while dust attenuation defines the location and shape of the blue and red populations. While a naive comparison between observation and simulations displays the known issue of over-quenching of satellite galaxies, the introduction of empirically-motivated observational errors and classification from the same group finder used in GAMA greatly reduces this tension. The introduction of random re-assignment of $\sim 15\%$ of centrals/satellites as satellites/centrals on the simulation classification closely resembles the outcome of the group finder, providing a computationally less intensive method to compare simulations with observations.

Observers experience a series of limitations when measuring galaxy kinematics, such as variable seeing conditions and aperture size. These effects can be reduced using empirical corrections, but these equations are usually applicable within a restrictive set of boundary conditions (e.g. S\'ersic indices within a given range) which can lead to biases when trying to compare measurements made across a full kinematic survey. In this work, we present new corrections for two widely used kinematic parameters, $\lambda_R$ and $V/\sigma$, that are applicable across a broad range of galaxy shapes, measurement radii and ellipticities. We take a series of mock observations of N-body galaxy models and use these to quantify the relationship between the observed kinematic parameters, structural properties and different seeing conditions. Derived corrections are then tested using the full catalogue of galaxies, including hydro-dynamic models from the EAGLE simulation. Our correction is most effective for regularly-rotating systems, yet the kinematic parameters of all galaxies -- fast, slow and irregularly rotating systems -- are recovered successfully. We find that $\lambda_R$ is more easily corrected than $V/\sigma$, with relative deviations of 0.02 and 0.06 dex respectively. The relationship between $\lambda_R$ and $V/\sigma$, as described by the parameter $\kappa$, also has a minor dependence on seeing conditions. These corrections will be particularly useful for stellar kinematic measurements in current and future integral field spectroscopic (IFS) surveys of galaxies.

We use KiDS weak lensing data to measure variations in mean halo mass as a function of several key galaxy properties (namely: stellar colour, specific star formation rate, Sersic index, and effective radius) for a volume-limited sample of GAMA galaxies in a narrow stellar mass range ($M_* \sim 2$-$5 \times 10^{10}$ Msol). This mass range is particularly interesting, inasmuch as it is where bimodalities in galaxy properties are most pronounced, and near to the break in both the galaxy stellar mass function and the stellar-to-halo mass relation (SHMR). In this narrow mass range, we find that both size and Sersic index are better predictors of halo mass than colour or SSFR, with the data showing a slight preference for Sersic index. In other words, we find that mean halo mass is more strongly correlated with galaxy structure than either stellar populations or star formation rate. These results lead to an approximate lower bound on the dispersion in halo masses among $ \log M_* \approx {10.5}$ galaxies: we find that the dispersion is $\gtrsim 0.3$ dex. This would imply either that offsets from the mean SHMR are very tightly coupled to size/structure, or that the dispersion in the SHMR is larger than past results have suggested. Our results thus provide new empirical constraints on the relationship between stellar and halo mass assembly at this particularly interesting mass range.

We present the KMOS-CLASH (K-CLASH) survey, a K-band Multi-Object Spectrograph (KMOS) survey of the spatially-resolved gas properties and kinematics of 191 (predominantly blue) H$\alpha$-detected galaxies at $0.2 \lesssim z \lesssim 0.6$ in field and cluster environments. K-CLASH targets galaxies in four Cluster Lensing And Supernova survey with Hubble (CLASH) fields in the KMOS $IZ$-band, over $7'$ radius ($\approx2$-$3$ Mpc) fields-of-view. K-CLASH aims to study the transition of star-forming galaxies from turbulent, highly star-forming disc-like and peculiar systems at $z\approx1$-$3$, to the comparatively quiescent, ordered late-type galaxies at $z\approx0$, and to examine the role of clusters in the build-up of the red sequence since $z\approx1$. In this paper, we describe the K-CLASH survey, present the sample, and provide an overview of the K-CLASH galaxy properties. We demonstrate that our sample comprises star-forming galaxies typical of their stellar masses and epochs, residing both in field and cluster environments. We conclude K-CLASH provides an ideal sample to bridge the gap between existing large integral-field spectroscopy surveys at higher and lower redshifts. We find that star-forming K-CLASH cluster galaxies at intermediate redshifts have systematically lower stellar masses than their star-forming counterparts in the field, hinting at possible "downsizing" scenarios of galaxy growth in clusters at these epochs. We measure no difference between the star-formation rates of H$\alpha$-detected, star-forming galaxies in either environment after accounting for stellar mass, suggesting that cluster quenching occurs very rapidly during the epochs probed by K-CLASH, or that star-forming K-CLASH galaxies in clusters have only recently arrived there, with insufficient time elapsed for quenching to have occured.

We apply the spectral energy distribution-fitting code ProSpect to multiwavelength imaging for $\sim$7,000 galaxies from the GAMA survey at $z<0.06$, in order to extract their star-formation histories. We combine a parametric description of the star formation history with a closed-box evolution of metallicity where the present-day gas-phase metallicity of the galaxy is a free parameter. We show with this approach that we are able to recover the observationally-determined cosmic star formation history (CSFH), an indication that stars are being formed in the correct epoch of the Universe, on average, for the manner in which we are conducting SED fitting. We also show the contribution to the CSFH of galaxies of different present-day visual morphologies, and stellar masses. This analysis suggests that half of the mass in present-day elliptical galaxies was in place 11 Gyr ago, whereas in other morphological types the stellar mass formed later, up to 6 Gyr ago for present-day irregular galaxies. Similarly, the most massive galaxies in our sample were shown to have formed half their stellar mass by 10.5 Gyr ago, whereas the least massive galaxies formed half their stellar mass as late as 4 Gyr ago (the well-known effect of "galaxy downsizing"). Finally, our metallicity approach allows us to follow the average evolution in gas-phase metallicity for populations of galaxies, and extract the evolution of the cosmic metal mass density in stars and in gas, producing results in broad agreement with observations of metal densities in the Universe.

The Galaxy And Mass Assembly Survey (GAMA) covers five fields with highly complete spectroscopic coverage ($>95$ per cent) to intermediate depths ($r<19.8$ or $i < 19.0$ mag), and collectively spans 250 square degrees of Equatorial or Southern sky. Four of the GAMA fields (G09, G12, G15 and G23) reside in the ESO VST KiDS and ESO VISTA VIKING survey footprints, which combined with our GALEX, WISE and Herschel data provide deep uniform imaging in the $FUV\,NUV\,ugriZYJHK_s\,W1\,W2\,W3\,W4\,P100\,P160\,S250\,S350\,S500$ bands. Following the release of KiDS DR4, we describe the process by which we ingest the KiDS data into GAMA (replacing the SDSS data previously used for G09, G12 and G15), and redefine our core optical and near-IR catalogues to provide a complete and homogeneous dataset. The source extraction and analysis is based on the new ProFound image analysis package, providing matched-segment photometry across all bands. The data are classified into stars, galaxies, artefacts, and ambiguous objects, and objects are linked to the GAMA spectroscopic target catalogue. Additionally, a new technique is employed utilising ProFound to extract photometry in the unresolved MIR-FIR regime. The catalogues including the full FUV-FIR photometry are described and will be fully available as part of GAMA DR4. They are intended for both standalone science, selection for targeted follow-up with 4MOST, as well as an accompaniment to the upcoming and ongoing radio arrays now studying the GAMA $23^h$ field.

We present early science results from the First Large Absorption Survey in HI (FLASH), a spectroscopically blind survey for 21-cm absorption lines in cold hydrogen HI gas at cosmological distances using the Australian Square Kilometre Array Pathfinder (ASKAP). We have searched for HI absorption towards 1253 radio sources in the GAMA 23 field, covering redshifts between $z = 0.34$ and $0.79$ over a sky area of approximately 50 deg$^{2}$. In a purely blind search we did not obtain any detections of 21-cm absorbers above our reliability threshold. Assuming a fiducial value for the HI spin temperature of $T_{\rm spin}$ = 100 K and source covering fraction $c_{\rm f} = 1$, the total comoving absorption path length sensitive to all Damped Lyman $\alpha$ Absorbers (DLAs; $N_{\rm HI} \geq 2 \times 10^{20}$ cm$^{-2}$) is $\Delta{X} = 6.6 \pm 0.3$ ($\Delta{z} = 3.7 \pm 0.2$) and super-DLAs ($N_{\rm HI} \geq 2 \times 10^{21}$ cm$^{-2}$) is $\Delta{X} = 111 \pm 6$ ($\Delta{z} = 63 \pm 3$). We estimate upper limits on the HI column density frequency distribution function that are consistent with measurements from prior surveys for redshifted optical DLAs, and nearby 21-cm emission and absorption. By cross matching our sample of radio sources with optical spectroscopic identifications of galaxies in the GAMA 23 field, we were able to detect 21-cm absorption at $z = 0.3562$ towards NVSS J224500$-$343030, with a column density of $N_{\rm HI} = (1.2 \pm 0.1) \times 10^{20} (T_{\rm spin}/100~\mathrm{K})$ cm$^{-2}$. The absorber is associated with GAMA J22450.05$-$343031.7, a massive early-type galaxy at an impact parameter of 17 kpc with respect to the radio source and which may contain a massive ($M_{\rm HI} \gtrsim 3 \times 10^{9}$ M$_{\odot}$) gas disc. Such gas-rich early types are rare, but have been detected in the nearby Universe.

We introduce ProSpect, a generative galaxy spectral energy distribution (SED) package that encapsulates the best practices for SED methodologies in a number of astrophysical domains. ProSpect comes with two popular families of stellar population libraries (BC03 and EMILES), and a large variety of methods to construct star formation and metallicity histories. It models dust through the use of a Charlot & Fall attenuation model, with re-emission using Dale far-infrared templates. It also has the ability to model AGN through the inclusion of a simple AGN and hot torus model. Finally, it makes use of MAPPINGS-III photoionisation tables to produce line emission features. We test the generative and inversion utility of ProSpect through application to the Shark galaxy formation semi-analytic code, and informed by these results produce fits to the final ultraviolet to far-infrared photometric catalogues produces by the Galaxy and Mass Assembly Survey (GAMA). As part of the testing of ProSpect, we also produce a range of simple photometric stellar mass approximations covering a range of filters for both observed frame and rest frame photometry.

We present SimSpin, a new, public, software framework for generating integral field spectroscopy (IFS) data cubes from N-body/hydrodynamical simulations of galaxies, which can be compared directly with observational datasets. SimSpin provides a consistent method for studying a galaxy's stellar component. It can be used to explore how observationally inferred measurements of kinematics, such as the spin parameter $\lambda_R$, are impacted by the effects of, for example, inclination, seeing conditions, distance, etc. SimSpin is written in R and has been designed to be highly modular, flexible, and extensible. It is already being used by the astrophysics community to generate IFS-like cubes and FITS files for direct comparison of simulations to observations. In this paper, we explain the conceptual framework of SimSpin; how it is implemented in R; and we demonstrate SimSpin's current capabilities, providing as an example a brief investigation of how numerical resolution affects how reliably we can recover the intrinsic stellar kinematics of a simulated galaxy.

We use our catalogue of structural decomposition measurements for the extended GALEX Arecibo SDSS Survey (xGASS) to study the role of bulges both along and across the galaxy star-forming main sequence (SFMS). We show that the slope in the $sSFR$-$M_{\star}$ relation flattens by $\sim$0.1 dex per decade in $M_{\star}$ when re-normalising $sSFR$ by disc stellar mass instead of total stellar mass. However, recasting the $sSFR$-$M_{\star}$ relation into the framework of only disc-specific quantities shows that a residual trend remains against disc stellar mass with equivalent slope and comparable scatter to that of the total galaxy relation. This suggests that the residual declining slope of the SFMS is intrinsic to the disc components of galaxies. We further investigate the distribution of bulge-to-total ratios ($B/T$) as a function of distance from the SFMS ($\Delta SFR_{MS}$). At all stellar masses, the average $B/T$ of local galaxies decreases monotonically with increasing $\Delta SFR_{MS}$. Contrary to previous works, we find that the upper-envelope of the SFMS is not dominated by objects with a significant bulge component. This rules out a scenario in which, in the local Universe, objects with increased star formation activity are simultaneously experiencing a significant bulge growth. We suggest that much of the discrepancies between different works studying the role of bulges originates from differences in the methodology of structurally decomposing galaxies.

We present a detailed analysis of the radio galaxy PKS 2250-351, a giant of 1.2 Mpc projected size, its host galaxy, and its environment. We use radio data from the Murchison Widefield Array, the upgraded Giant Metre-wavelength Radio Telescope, the Australian Square Kilometre Array Pathfinder, and the Australia Telescope Compact Array to model the jet power and age. Optical and infra-red data come from the Galaxy And Mass Assembly (GAMA) survey and provide information on the host galaxy and environment. GAMA spectroscopy confirms that PKS 2250-351 lies at z=0.2115 in the irregular, and likely unrelaxed, cluster Abell 3936. We find its host is a massive, `red and dead' elliptical galaxy with negligible star formation but with a highly obscured active galactic nucleus dominating the mid-infrared emission. Assuming it lies on the local M-sigma relation it has an Eddington accretion rate of lambda_EDD~0.014. We find that the lobe-derived jet power (a time-averaged measure) is an order of magnitude greater than the hotspot-derived jet power (an instantaneous measure). We propose that over the lifetime of the observed radio emission (~300 Myr) the accretion has switched from an inefficient advection dominated mode to a thin-disc efficient mode, consistent with the decrease in jet power. We also suggest that the asymmetric radio morphology is due to its environment, with the host of PKS 2250-351 lying to the west of the densest concentration of galaxies in Abell 3936.

We use data from the GAMA and GALEX surveys to demonstrate that the UV upturn, an unexpected excess of ultraviolet flux from a hot stellar component, seen in the spectra of many early-type galaxies, arises from processes internal to individual galaxies with no measurable influence from the galaxies' larger environment. We first define a clean sample of passive galaxies without a significant contribution to their UV flux from low-level star formation. We confirm that galaxies with the optical colours of red sequence galaxies often have signs of residual star formation, which, without other information, would prevent a convincing demonstration of the presence of UV upturns. However, by including (NUV$-u$) and {\it WISE} (W2-W3) colours, and FUV data where it exists, we can convincingly constrain samples to be composed of non-star-forming objects. Using such a sample, we examine GALEX photometry of low redshift GAMA galaxies in a range of low-density environments, from groups to the general field, searching for UV upturns. We find a wide range of (NUV$-r$) colours, entirely consistent with the range seen -- and attributed to the UV upturn -- in low-redshift red sequence cluster galaxies. The range of colours is independent of group multiplicity or velocity dispersion, with isolated passive galaxies just as likely to have blue UV-to-optical colours, implying significant upturn components, as those in richer groups and in the previous data on clusters. This is supported by equivalent results for (FUV$-r$) colours which are clear indicators of upturn components.

Hierarchical models of structure formation predict that dark matter halo assembly histories are characterised by episodic mergers and interactions with other haloes. An accurate description of this process will provide insights into the dynamical evolution of haloes and the galaxies that reside in them. Using large cosmological N-body simulations, we characterise halo orbits to study the interactions between substructure haloes and their hosts, and how different evolutionary histories map to different classes of orbits. We use two new software tools - WhereWolf, which uses halo group catalogues and merger trees to ensure that haloes are tracked accurately in dense environments, and OrbWeaver, which quantifies each halo's orbital parameters. We demonstrate how WhereWolf improves the accuracy of halo merger trees, and we use OrbWeaver to quantify orbits of haloes. We assess how well analytical prescriptions for the merger timescale from the literature compare to measured merger timescales from our simulations and find that existing prescriptions perform well, provided the ratio of substructure-to-host mass is not too small. In the limit of small substructure-to-host mass ratio, we find that the prescriptions can overestimate the merger timescales substantially, such that haloes are predicted to survive well beyond the end of the simulation. This work highlights the need for a revised analytical prescription for the merger timescale that more accurately accounts for processes such as catastrophic tidal disruption.

We present a study of the dust, stars and atomic gas (HI) in an HI-selected sample of local galaxies (z<0.035) in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) fields. This HI-selected sample reveals a population of very high gas fraction (>80 per cent), low stellar mass sources that appear to be in the earliest stages of their evolution. We compare this sample with dust and stellar mass selected samples to study the dust and gas scaling relations over a wide range of gas fraction (proxy for evolutionary state of a galaxy). The most robust scaling relations for gas and dust are those linked to NUV-r (SSFR) and gas fraction, these do not depend on sample selection or environment. At the highest gas fractions, our additional sample shows the dust content is well below expectations from extrapolating scaling relations for more evolved sources, and dust is not a good tracer of the gas content. The specific dust mass for local galaxies peaks at a gas fraction of ~75 per cent. The atomic gas depletion time is also longer for high gas fraction galaxies, opposite to the trend found for molecular gas depletion timescale. We link this trend to the changing efficiency of conversion of HI to H2 as galaxies increase in stellar mass surface density as they evolve. Finally, we show that galaxies start out barely obscured and increase in obscuration as they evolve, yet there is no clear and simple link between obscuration and global galaxy properties.

We present a structural decomposition analysis of the galaxies in the extended GALEX Arecibo SDSS Survey (xGASS) using (gri) images from the Sloan Digital Sky Survey. Utilising the 2D Bayesian light profile fitting code ProFit, we fit single- and double-component models taking advantage of a robust Markov chain Monte Carlo optimisation algorithm in which we assume a Sersic profile for single-component models and a combination of a Sersic bulge and near-exponential disc (0.5 < n < 1.5) for double-component models. We investigate the effect of bulges on the atomic hydrogen (HI) content in galaxies by revisiting the HI-to-stellar mass scaling relations with the bulge-to-total ratio measured in the ProFit decompositions. We show that, at both fixed total and disc stellar mass, more bulge-dominated galaxies have systematically lower HI masses, implying that bulge-dominated galaxies with large HI reservoirs are rare in the local Universe. We see similar trends when separating galaxies by a bulge-to-total ratio based either on luminosity or stellar mass, however, the trends are more evident with luminosity. Importantly, when controlling for both stellar mass and star formation rate, the separation of atomic gas content reduces to within 0.3 dex between galaxies of different bulge-to-total ratios. Our findings suggest that the presence of a photometric bulge has little effect on the global HI gas reservoirs of local galaxies.

The Galaxy And Mass Assembly (GAMA) survey has morphologically identified a class of "Little Blue Spheroid" (LBS) galaxies whose relationship to other classes of galaxies we now examine in detail. Considering a sample of 868 LBSs, we find that such galaxies display similar but not identical colours, specific star formation rates, stellar population ages, mass-to-light ratios, and metallicities to Sd-Irr galaxies. We also find that LBSs typically occupy environments of even lower density than those of Sd-Irr galaxies, where ~65% of LBS galaxies live in isolation. Using deep, high-resolution imaging from VST KiDS and the new Bayesian, two-dimensional galaxy profile modeling code PROFIT, we further examine the detailed structure of LBSs and find that their S\'ersic indices, sizes, and axial ratios are compatible with those of low-mass elliptical galaxies. We then examine SAMI Galaxy survey integral field emission line kinematics for a subset of 62 LBSs and find that the majority (42) of these galaxies display ordered rotation with the remainder displaying disturbed/non-ordered dynamics. Finally, we consider potential evolutionary scenarios for a population with this unusual combination of properties, concluding that LBSs are likely formed by a mixture of merger and accretion processes still recently active in low-redshift dwarf populations. We also infer that if LBS-like galaxies were subjected to quenching in a rich environment, they would plausibly resemble cluster dwarf ellipticals.

We combine the Shark semi-analytic model of galaxy formation with the ProSpect software tool for spectral energy distribution (SED) generation to study the multi-wavelength emission of galaxies from the far-ultraviolet (FUV) to the far-infrared (FIR) at $0\le z\le 10$. We produce a physical model for the attenuation of galaxies across cosmic time by combining a local Universe empirical relation to compute the dust mass of galaxies from their gas metallicity and mass, attenuation curves derived from radiative transfer calculations of galaxies in the EAGLE hydrodynamic simulation suite, and the properties of Shark galaxies. We are able to produce a wide range of galaxies, from the $z=8$ star-forming galaxies with almost no extinction, $z=2$ submillimeter galaxies, down to the normal star-forming and red sequence galaxies at $z=0$. Quantitatively, we find that Shark reproduces the observed (i) the $z=0$ FUV-to-FIR, (ii) $0\le z\le 3$ rest-frame $K$-band, and (iii) $0\le z\le 10$ rest-frame FUV luminosity functions, (iv) $z\le 8$ UV slopes, (v) the FUV-to-FIR number counts (including the widely disputed 850$\mu$m), (vi) redshift distribution of bright $850\mu$m galaxies and (vii) the integrated cosmic SED from $z=0$ to $z=1$ to an unprecedented level. This is achieved without the need to invoke changes in the stellar initial mass function, dust-to-metal mass ratio, or metal enrichment timescales. Our model predicts star formation in galaxy disks to dominate in the FUV-to-optical, while bulges dominate at the NIR at all redshifts. The FIR sees a strong evolution in which disks dominate at $z\le 1$ and starbursts (triggered by both galaxy mergers and disk instabilities, in an even mix) dominate at higher redshifts, even out to $z=10$.

We present a simple model for galaxy attenuation by distilling SKIRT radiative transfer calculations for ~100,000 EAGLE galaxies at redshifts z=2-0. Our model adapts the two component screen model of Charlot & Fall (2000), parametrising the optical depth and slope of the ISM screen using the average dust surface density, $\Sigma_{\rm dust}$. We recover relatively tight relations between these parameters for the EAGLE sample, but also provide the scatter in these parameter owing to the morphological variation and orientation of galaxies. We also find that these relations are nearly independent of redshift in the EAGLE model. By pairing our model with an empirical prescription for birth clouds below the resolution scale of the simulation, we reproduce the observed relation between attenuation slope and optical depth for the first time in a cosmological simulation. We demonstrate that this result is remarkably independent of the attenuation properties assumed for birth cloud screen, merely requiring a boosted attenuation for infant stars. We present this model with a view to interpreting observations, as well as processing semi-analytic models and other hydrodynamic simulations.

In the current era of radio astronomy, continuum surveys observe a multitude of objects with complex morphologies and sizes, and are not limited to observing point sources. Typical radio source extraction software generates catalogues by using Gaussian components to form a model of the emission. This may not be well suited to complicated jet structures and extended emission, particularly in the era of interferometers with a high density of short baselines, that are sensitive to extended emission. In this paper, we investigate how the optically-motivated source detection package ProFound (Robotham et al. 2018) may be used to model radio emission of both complicated and point-like radio sources. We use a combination of observations and simulations to investigate how ProFound compares to other source extractor packages used for radio surveys. We find that ProFound can accurately recover both the flux densities of simulated Gaussian sources as well as extended radio galaxies. ProFound can create models that trace the complicated nature of these extended galaxies, which we show is not necessarily the case with other source extraction software. Our work suggests that our knowledge of the emission from extended radio objects may be both over or under-estimated using traditional software. We suggest that ProFound offers a useful alternative to the fitting of Gaussian components for generating catalogues from current and future radio surveys. Furthermore, ProFound's multi-wavelength capabilities will be useful in investigating radio sources in combination with multiwavelength data.

The structural properties of individual dark matter haloes, including shape, spin, concentration, and substructure, are linked to the halo's growth history, but the exact connection between the two is unclear. One open question, in particular, is the effect of major mergers on halo structure. We have performed a large set of simulations of binary equal-mass mergers between isolated haloes with various density profiles, to map out the relationship between the initial conditions and merger parameters and the structure of the final remnant. In this paper we describe our initial set-up and analysis methods, and report on the results for the size, shape, and spin of the merger remnant. The outcomes of mergers are most easily understood in terms of a scaled dimensionless energy parameter $\kappa$ and an angular momentum (or spin) parameter $\lambda$. We find that the axis ratio $c/a$ scales roughly linearly with energy $\kappa$ while the axis ratio $c/b$ scales linearly with spin $\lambda$. Qualitatively, mergers on radial orbits produce prolate remnants, while mergers on tangential orbits produce oblate remnants. The spin of the remnant can be predicted from angular momentum conservation, while the overall size changes as $\sim \kappa^{-5}$, as expected from self-similar scaling at constant mean density. We discuss potential cosmological applications for these simple patterns.

Several lines of evidence suggest that as dark matter haloes grow their scale radius increases, and that the density in their central region drops. Major mergers seem an obvious mechanism to explain both these phenomena, and the resulting patterns in the concentration--mass--redshift relation. To test this possibility, we have simulated equal-mass mergers between haloes with a variety of cosmological density profiles, placed on various different orbits. The remnants typically have higher densities than the initial conditions, but differ only slightly from self-similar scaling predictions. They are reasonably well fit by Einasto profiles, but have parameters distinct from those of the initial conditions. The net internal energy available to the merger remnant, relative to the internal energy of the initial conditions, $\kappa$, has the greatest influence on the properties of the final mass distribution. As expected, energetic encounters produce more extended remnants while mergers of strongly bound systems produce compact remnants. Surprisingly, however, the scale radius of the density profile shows the opposite trend, {\it increasing} in the remnants of low-energy encounters relative to energetic ones. Also even in the most energetic encounters, the density within the scale radius decreases only slightly (by 10--20\%), while for very low-energy systems it increases significantly after the merger. We conclude that while major mergers can produce remnants that are more diffuse at large radii, they are relatively ineffective at changing the central densities of haloes, and seem unlikely to explain the mean trends in the concentration--mass--redshift relation.

We present TreeFrog, a massively parallel halo merger tree builder that is capable comparing different halo catalogues and producing halo merger trees. The code is written in c++11, use the MPI and OpenMP API's for parallelisation, and includes python tools to read/manipulate the data products produced. The code correlates binding energy sorted particle ID lists between halo catalogues, determining optimal descendant/progenitor matches using multiple snapshots, a merit function that maximises the number of shared particles using pseudo-radial moments, and a scheme for correcting halo merger tree pathologies. Focusing on VELOCIraptor catalogues for this work, we demonstrate how searching multiple snapshots spanning a dynamical time significantly reduces the number of stranded halos, those lacking a descendant or a progenitor, critically correcting poorly resolved halos. We present a new merit function that improves the distinction between primary and secondary progenitors, reducing tree pathologies. We find FOF accretion rates and merger rates show similar mass ratio dependence. The model merger rates from Poole et al, (2017) agree with the measured net growth of halos through mergers.

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

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.

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 current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium.

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 VELOCIraptor, a massively parallel galaxy/(sub)halo finder that is also capable of robustly identifying tidally disrupted objects and separate stellar halos from galaxies. The code is written in c++11, use the MPI and OpenMP API's for parallelisation, and includes python tools to read/manipulate the data products produced. We demonstrate the power of the VELOCIraptor (sub)halo finder, showing how it can identify subhalos deep within the host that have negligible density contrasts to their parent halo. We find a subhalo mass-radial distance dependence: large subhalos with mass ratios of $\gtrsim10^{-2}$ are more common in the central regions that smaller subhalos, a result of dynamical friction and low tidal mass loss rates. This dependence is completely absent in (sub)halo finders in common use, which generally search for substructure in configuration space, yet is present in codes that track particles belonging to halos as they fall into other halos, such as HBT+. VELOCIraptor largely reproduces the dependence seen without tracking, finding a similar radial dependence to HBT+ in well resolved halos from our limited resolution fiducial simulation.

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as stripping, strangulation and starvation. Here we explore this further using the Galaxy And Mass Assembly (GAMA) survey which extends to ~2dex lower in stellar mass than SDSS, is more complete for closely-separated galaxies (>95% compared to >70%), and identifies lower-halo-mass groups outside of the very local Universe (M$_{\mathrm{halo}}\sim10^{12}$M$_{\odot}$ at $0.1<z<0.2$). As far as possible we aim to replicate the selections, completeness corrections and central/satellite division of one of the previous studies but find clear differences between passive fractions of centrals and satellites. We also find that our passive fractions increase with both halo-to-satellite mass ratio and central-to-second rank mass ratio. This suggests that quenching is more efficient in satellites that are low-mass for their halo ($i.e$ at high halo-to-satellite mass ratio in comparison to low halo-to-satellite mass ratio) and are more likely to be passive in older groups - forming a consistent picture of environmental quenching of satellites. We then discuss potential explanations for the previously observed similarity, such as dependence on the group-finding method.

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

Misalignment of gas and stellar rotation in galaxies can give clues to the origin and processing of accreted gas. Integral field spectroscopic observations of 1213 galaxies from the SAMI Galaxy Survey show that 11% of galaxies with fitted gas and stellar rotation are misaligned by more than 30 degrees in both field/group and cluster environments. Using SAMI morphological classifications and Sersic indices, the misalignment fraction is 45+/-6% in early-type galaxies, but only 5+/-1% in late-type galaxies. The distribution of position angle offsets is used to test the physical drivers of this difference. Slower dynamical settling time of the gas in elliptical stellar mass distributions accounts for a small increase in misalignment in early-type galaxies. However, gravitational dynamical settling time is insufficient to fully explain the observed differences between early- and late-type galaxies in the distributions of the gas/stellar position angle offsets. LTGs have primarily accreted gas close to aligned rather than settled from misaligned based on analysis of the skewed distribution of PA offsets compared to a dynamical settling model. Local environment density is less important in setting the misalignment fractions than morphology, suggesting that mergers are not the main source of accreted gas in these disks. Cluster environments are found to have gas misalignment driven primarily by cluster processes not by gas accretion.

A primary goal of integral field spectroscopic (IFS) surveys is to provide a statistical census of galaxies classified by their internal kinematics. As a result, the observational spin parameter, $\lambda_R$, has become one of the most popular methods of quantifying the relative importance of velocity dispersion and rotation in supporting a galaxy's inner structure. The goal of this paper is to examine the relationship between the observationally deduced $\lambda_R$ and one of the most commonly used theoretical spin parameters in the literature, the Bullock et al. (2001) $\lambda'$. Using a set of $N$-body realisations of galaxies from which we construct mock IFS observations, we measure $\lambda_R$ as an observer would, incorporating the effects of beam smearing and seeing conditions. Assuming parameters typical of current IFS surveys, we confirm that there are strong positive correlations between $\lambda_R$ and measurement radius, and strong negative correlations between $\lambda_R$ and size of the PSF, for late-type galaxies; these biases can be reduced using a recently proposed empirical correction. Once observational biases are corrected for, we find that $\lambda_R$ provides a good approximation to $\sim \sqrt{3}/2 \; \lambda'(\rm R_{\rm eff})$, where $\lambda'$ is evaluated for the galactic stellar component within 1 R$_{\rm eff}$.

Recently a number of studies have proposed that the dispersion along the star formation rate - stellar mass relation ($\sigma_{\mathrm{sSFR}}$-M$_{*}$) is indicative of variations in star-formation history (SFH) driven by feedback processes. They found a 'U'-shaped dispersion and attribute the increased scatter at low and high stellar masses to stellar and active galactic nuclei feed-back respectively. However, measuring $\sigma_{\mathrm{sSFR}}$ and the shape of the $\sigma_{\mathrm{sSFR}}$-M$_{*}$ relation is problematic and can vary dramatically depending on the sample selected, chosen separation of passive/star-forming systems, and method of deriving star-formation rates ($i.e.$ H$\alpha$ emission vs spectral energy distribution fitting). As such, any astrophysical conclusions drawn from measurements of $\sigma_{\mathrm{sSFR}}$ must consider these dependencies. Here we use the Galaxy And Mass Assembly survey to explore how $\sigma_{\mathrm{sSFR}}$ varies with SFR indicator for a variety of selections for disc-like `main sequence' star-forming galaxies including colour, star-formation rate, visual morphology, bulge-to-total mass ratio, S\'{e}rsic index and mixture modelling. We find that irrespective of sample selection and/or SFR indicator, the dispersion along the sSFR-M$_{*}$ relation does follow a 'U'-shaped distribution. This suggests that the shape is physical and not an artefact of sample selection or method. We then compare the $\sigma_{\mathrm{sSFR}}$-M$_{*}$ relation to state-of-the-art hydrodynamical and semi-analytic models and find good agreement with our observed results. Finally, we find that for group satellites this 'U'-shaped distribution is not observed due to additional high scatter populations at intermediate stellar masses.

We present Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY) early science results from the Australian Square Kilometre Array Pathfinder (ASKAP) observations of the NGC 7162 galaxy group. We use archival HIPASS and Australia Telescope Compact Array (ATCA) observations of this group to validate the new ASKAP data and the data reduction pipeline ASKAPsoft. We detect six galaxies in the neutral hydrogen (HI) 21-cm line, expanding the NGC 7162 group membership from four to seven galaxies. Two of the new detections are also the first HI detections of the dwarf galaxies, AM 2159-434 and GALEXASC J220338.65-431128.7, for which we have measured velocities of $cz=2558$ and $cz=2727$ km s$^{-1}$, respectively. We confirm that there is extended HI emission around NGC 7162 possibly due to past interactions in the group as indicated by the $40^{\circ}$ offset between the kinematic and morphological major axes for NGC 7162A, and its HI richness. Taking advantage of the increased resolution (factor of $\sim1.5$) of the ASKAP data over archival ATCA observations, we fit a tilted ring model and use envelope tracing to determine the galaxies' rotation curves. Using these we estimate the dynamical masses and find, as expected, high dark matter fractions of $f_{\mathrm{DM}}\sim0.81-0.95$ for all group members. The ASKAP data are publicly available.

In the currently debated context of using clusters of galaxies as cosmological probes, the need for well-defined cluster samples is critical. The XXL Survey has been specifically designed to provide a well characterised sample of some 500 X-ray detected clusters suitable for cosmological studies. The main goal of present article is to make public and describe the properties of the cluster catalogue in its present state, as well as of associated catalogues as super-clusters and fossil groups. We release a sample containing 365 clusters in total. We give the details of the follow-up observations and explain the procedure adopted to validate the cluster spectroscopic redshifts. Considering the whole XXL cluster sample, we have provided two types of selection, both complete in a particular sense: one based on flux-morphology criteria, and an alternative based on the [0.5-2] keV flux within one arcmin of the cluster centre. We have also provided X-ray temperature measurements for 80$\%$ of the clusters having a flux larger than 9$\times$10$^{-15}$$\rm \thinspace erg \, s^{-1} \, cm^{-2}$. Our cluster sample extends from z$\sim$0 to z$\sim$1.2, with one cluster at z$\sim$2. Clusters were identified through a mean number of six spectroscopically confirmed cluster members. Our updated luminosity function and luminosity-temperature relation are compatible with our previous determinations based on the 100 brightest clusters, but show smaller uncertainties. We also present an enlarged list of super-clusters and a sample of 18 possible fossil groups. This intermediate publication is the last before the final release of the complete XXL cluster catalogue when the ongoing C2 cluster spectroscopic follow-up is complete. It provides a unique inventory of medium-mass clusters over a 50~\dd\ area out to z$\sim$1.

We present the version of the point source catalogue of the XXL Survey that was used, in part, in the first series of XXL papers. In this paper we release, in our database in Milan and at CDS: (i) the X-ray source catalogue with 26056 objects in two areas of 25 deg2; (ii) the associated multiwavelength catalogues with candidate counterparts of the X-ray sources in the infrared, near-infrared, optical, and ultraviolet (plus spectroscopic redshift when available); and (iii) a catalogue of spectroscopic redshifts recently obtained in the southern XXL area. We also present the basic properties of the X-ray point sources and their counterparts. Other catalogues described in the second series of XXL papers will be released contextually, and will constitute the second XXL data release.

We present a new, open source, free semi-analytic model (SAM) of galaxy formation, Shark, designed to be highly flexible and modular, allowing easy exploration of different physical processes and ways of modelling them. We introduce the philosophy behind Shark and provide an overview of the physical processes included in the model. Shark is written in C++11 and has been parallelized with OpenMP. In the released version (v1.1), we implement several different models for gas cooling, active galactic nuclei, stellar and photo-ionisation feedback, and star formation (SF). We demonstrate the basic performance of Shark using the Planck15 cosmology SURFS simulations, by comparing against a large set of observations, including: the stellar mass function (SMF) and stellar-halo mass relation at z=0-4; the cosmic evolution of the star formation rate density (SFRD), stellar mass, atomic and molecular hydrogen; local gas scaling relations; and structural galaxy properties, finding excellent agreement. Significant improvements over previous SAMs are seen in the mass-size relation for disks/bulges, the gas-stellar mass and stellar mass-metallicity relations. To illustrate the power of Shark in exploring the systematic effects of the galaxy formation modelling, we quantify how the scatter of the SF main sequence and the gas scaling relations changes with the adopted SF law, and the effect of the starbursts H$_2$ depletion timescale on the SFRD and $\Omega_{\rm H_2}$. We compare Shark with other SAMs and the hydrodynamical simulation EAGLE, and find that SAMs have a much higher halo baryon fractions due to large amounts of intra-halo gas, which in the case of EAGLE is in the intergalactic medium.

Merger trees harvested from cosmological $N$-body simulations encode the assembly histories of dark matter halos over cosmic time, and are a fundamental component of semi-analytical models (SAMs) of galaxy formation. The ability to compare the tools used to construct merger trees, namely halo finders and tree building algorithms, in an unbiased and systematic manner is critical to assess the quality of merger trees. In this paper, we present the dendogram, a novel method to visualise merger trees, which provides a comprehensive characterisation of a halo's assembly history - tracking subhalo orbits, halo merger events, and the general evolution of halo properties. We show the usefulness of the dendogram as a diagnostic tool of merger trees by comparing halo assembly histories from a single $N$-Body simulation analysed with three different halo-finders -\textsc{VELOCIraptor}, \textsc{AHF} and \textsc{Rockstar} - and their associated tree-builders. Based on our analysis of the resulting dendograms, we highlight how they have been used to motivate improvements to \textsc{VELOCIraptor}. The dendogram software is publicly available online, at: https://github.com/rhyspoulton/MergerTree-Dendograms .

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.

Using a combined and consistently analysed GAMA, G10-COSMOS, and 3D-HST dataset we explore the evolution of the galaxy stellar-mass function over lookback times $t_{\rm L} \in \left[0.2,12.5\right] {\rm h^{-1}_{70} Gyr}$. We use a series of volume limited samples to fit Schechter functions in bins of $\sim\!$constant lookback time and explore the evolution of the best-fit parameters in both single and two-component cases. In all cases, we employ a fitting procedure that is robust to the effects of Eddington bias and sample variance. Surprisingly, when fitting a two-component Schechter function, we find essentially no evidence of temporal evolution in $M_\star$, the two $\alpha$ slope parameters, or the normalisation of the low-mass component. Instead, our fits suggest that the various shape parameters have been exceptionally stable over cosmic time, as has the normalisation of the low-mass component, and that the evolution of the stellar-mass function is well described by a simple build up of the high-mass component over time. When fitting a single component Schechter function, there is an observed evolution in both $M_\star$ and $\alpha$, however this is interpreted as being an artefact. Finally, we find that the evolution of the stellar-mass function, and the observed stellar mass density, can be well described by a simple model of constant growth in the high-mass source density over the last $11 {\rm h^{-1}_{70} Gyr}$.

The Deep Extragalactic VIsible Legacy Survey (DEVILS) is a large spectroscopic campaign at the Anglo-Australian Telescope (AAT) aimed at bridging the near and distant Universe by producing the highest completeness survey of galaxies and groups at intermediate redshifts ($0.3<z<1.0$). Our sample consists of $\sim$60,000 galaxies to Y$<$21.2mag, over $\sim$6deg$^{2}$ in three well-studied deep extragalactic fields (Cosmic Origins Survey field, COSMOS, Extended Chandra Deep Field South, ECDFS and the X-ray Multi-Mirror Mission Large-Scale Structure region, XMM-LSS - all Large Synoptic Survey Telescope deep-drill fields). This paper presents the broad experimental design of DEVILS. Our target sample has been selected from deep Visible and Infrared Survey Telescope for Astronomy (VISTA) Y-band imaging (VISTA Deep Extragalactic Observations, VIDEO and UltraVISTA), with photometry measured by ProFound. Photometric star/galaxy separation is done on the basis of NIR colours, and has been validated by visual inspection. To maximise our observing efficiency for faint targets we employ a redshift feedback strategy, which continually updates our target lists, feeding back the results from the previous night's observations. We also present an overview of the initial spectroscopic observations undertaken in late 2017 and early 2018.

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}$.

In this short note we publish the analytic quantile function for the Navarro, Frenk & White (NFW) profile. All known published and coded methods for sampling from the 3D NFW PDF use either accept-reject, or numeric interpolation (sometimes via a lookup table) for projecting random Uniform samples through the quantile distribution function to produce samples of the radius. This is a common requirement in N-body initial condition (IC), halo occupation distribution (HOD), and semi-analytic modelling (SAM) work for correctly assigning particles or galaxies to positions given an assumed concentration for the NFW profile. Using this analytic description allows for much faster and cleaner code to solve a common numeric problem in modern astronomy. We release R and Python versions of simple code that achieves this sampling, which we note is trivial to reproduce in any modern programming language.

Superclusters form from the largest enhancements in the primordial density perturbation field and extend for tens of Mpc, tracing the large-scale structure of the Universe. We characterise XLSSsCN01, a rich supercluster at z~0.3 detected in the XXL Survey, composed of X-ray clusters of different virial masses and luminosities. As one of the first studies on this topic, we investigate the stellar populations of galaxies in different environments in the supercluster region. We study a magnitude-limited (r<=20) and a mass-limited sample (log(M*/M_sun)>10.8) of galaxies in the virialised region and in the outskirts of 11 XLSSsCN01 clusters, in high- and low-density field. We compute the stellar population properties of galaxies using spectral energy distribution and spectral fitting techniques, and study the dependence of star formation rates (SFR), colours, and stellar ages on environment. For r<20, the fraction of SFing/blue galaxies, computed either from the specific-SFR (sSFR) or rest-frame(rf) colour, shows depletion within the cluster virial radii, where the number of galaxies with log(sSFR/yr^-1)>-12 and with (g-r)_rf<0.6 is lower than in the field. For log(M*/M_sun)>10.8, no trends with environment emerge, as massive galaxies are mostly already passive in all environments. No differences among low- and high-density field members and cluster members emerge in the sSFR-mass relation. The luminosity-weighted age-mass relation of the passive populations within cluster virial radii show signatures of recent environmental quenching. The study of luminous and massive galaxies in this supercluster shows that while environment has a prominent role in determining the fractions of SFing/blue galaxies, its effects on the star formation activity in SFing galaxies are negligible.

We present a detailed analysis of the specific star formation rate -- stellar mass ($\mathrm{sSFR}-M_*$) of $z\le 0.13$ disk central galaxies using a morphologically selected mass-complete sample ($M_* \ge 10^{9.5} M_{\odot}$). Considering samples of grouped and ungrouped galaxies, we find the $\mathrm{sSFR}-M_*$ relations of disk-dominated central galaxies to have no detectable dependence on host dark-matter halo (DMH) mass, even where weak-lensing measurements indicate a difference in halo mass of a factor $\gtrsim5$. We further detect a gradual evolution of the $\mathrm{sSFR}-M_*$ relation of non-grouped (field) central disk galaxies with redshift, even over a $\Delta z \approx 0.04$ ($\approx5\cdot10^{8}\mathrm{yr}$) interval, while the scatter remains constant. This evolution is consistent with extrapolation of the "main-sequence-of-star-forming-galaxies" from previous literature that uses larger redshift baselines and coarser sampling. Taken together, our results present new constraints on the paradigm under which the SFR of galaxies is determined by a self-regulated balance between gas inflows and outflows, and consumption of gas by star-formation in disks, with the inflow being determined by the product of the cosmological accretion rate and a fuelling-efficiency -- $\dot{M}_{\mathrm{b,halo}}\zeta$. In particular, maintaining the paradigm requires $\dot{M}_{\mathrm{b,halo}}\zeta$ to be independent of the mass $M_{\mathrm{halo}}$ of the host DMH. Furthermore, it requires the fuelling-efficiency $\zeta$ to have a strong redshift dependence ($\propto (1+z)^{2.7}$ for $M_*=10^{10.3} M_{\odot}$ over $z=0 - 0.13$), even though no morphological transformation to spheroids can be invoked to explain this in our disk-dominated sample. The physical mechanisms capable of giving rise to such dependencies of $\zeta$ on $M_{\mathrm{halo}}$ and $z$ for disks are unclear.

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.

Satellite galaxies are commonly used as tracers to measure the line-of-sight velocity dispersion ($\sigma_{\rm LOS}$) of the dark matter halo associated with their central galaxy, and thereby to estimate the halo's mass. Recent observational dispersion estimates of the Local Group, including the Milky Way and M31, suggest $\sigma\sim$50 km/s, which is surprisingly low when compared to the theoretical expectation of $\sigma\sim$100s km/s for systems of their mass. Does this pose a problem for $\Lambda$CDM? We explore this tension using the {\small{SURFS}} suite of $N$-body simulations, containing over 10000 (sub)haloes with well tracked orbits. We test how well a central galaxy's host halo velocity dispersion can be recovered by sampling $\sigma_{\rm LOS}$ of subhaloes and surrounding haloes. Our results demonstrate that $\sigma_{\rm LOS}$ is biased mass proxy. We define an optimal window in $v_{\rm LOS}$ and projected distance ($D_p$) -- $0.5\lesssim D_p/R_{\rm vir}\lesssim1.0$ and $v_{\rm LOS} \lesssim0.5V_{\rm esc}$, where $R_{\rm vir}$ is the virial radius and $V_{\rm esc}$ is the escape velocity -- such that the scatter in LOS to halo dispersion is minimised - $\sigma_{\rm LOS}=(0.5\pm0.1)\sigma_{v,{\rm H}}$. We argue that this window should be used to measure line-of-sight dispersions as a proxy for mass, as it minimises scatter in the $\sigma_{\rm LOS}-M_{\rm vir}$ relation. This bias also naturally explains the results from \cite{mcconnachie2012a}, who used similar cuts when estimating $\sigma_{\rm LOS,LG}$, producing a bias of $\sigma_{\rm LG}=(0.44\pm0.14)\sigma_{v,{\rm H}}$. We conclude that the Local Group's velocity dispersion does not pose a problem for $\Lambda$CDM and has a mass of $\log M_{\rm LG, vir}/M_\odot=12.0^{+0.8}_{-2.0}$.

We present a new scheme, $\it{galtag}$, for refining the photometric redshift measurements of faint galaxies by probabilistically tagging them to observed galaxy groups constructed from a brighter, magnitude-limited spectroscopy survey. First, this method is tested on the DESI light-cone data constructed on the GALFORM galaxy formation model to tests its validity. We then apply it to the photometric observations of galaxies in the Kilo-Degree Imaging Survey (KiDS) over a 1 deg$^2$ region centred at 15$^\mathrm{h}$. This region contains Galaxy and Mass Assembly (GAMA) deep spectroscopic observations (i-band<22) and an accompanying group catalogue to r-band<19.8. We demonstrate that even with some trade-off in sample size, an order of magnitude improvement on the accuracy of photometric redshifts is achievable when using $\it{galtag}$. This approach provides both refined photometric redshift measurements and group richness enhancement. In combination these products will hugely improve the scientific potential of both photometric and spectroscopic datasets. The $\it{galtag}$ software will be made publicly available at https://github.com/pkaf/galtag.git.

We introduce ProFound, a source finding and image analysis package. ProFound provides methods to detect sources in noisy images, generate segmentation maps identifying the pixels belonging to each source, and measure statistics like flux, size and ellipticity. These inputs are key requirements of ProFit, our recently released galaxy profiling package, where the design aim is that these two software packages will be used in unison to semi-automatically profile large samples of galaxies. The key novel feature introduced in ProFound is that all photometry is executed on dilated segmentation maps that fully contain the identifiable flux, rather than using more traditional circular or ellipse based photometry. Also, to be less sensitive to pathological segmentation issues, the de-blending is made across saddle points in flux. We apply ProFound in a number of simulated and real world cases, and demonstrate that it behaves reasonably given its stated design goals. In particular, it offers good initial parameter estimation for ProFit, and also segmentation maps that follow the sometimes complex geometry of resolved sources, whilst capturing nearly all of the flux. A number of bulge-disc decomposition projects are already making use of the ProFound and ProFit pipeline, and adoption is being encouraged by publicly releasing the software for the open source R data analysis platform under an LGPL-3 license on GitHub (github.com/asgr/ProFound).

The spherical Jeans equation is a widely used tool for dynamical study of gravitating systems in astronomy. Here we test its efficacy in robustly weighing the mass of Milky Way analogues, given they need not be in equilibrium or even spherical. Utilizing Milky Way stellar halos simulated in accordance with $\Lambda{\rm CDM}$ cosmology by Bullock and Johnston (2005) and analysing them under the Jeans formalism, we recover the underlying mass distribution of the parent galaxy, within distance $r/{\rm kpc}\in[10,100]$, with a bias of $\sim12\%$ and a dispersion of $\sim14\%$. Additionally, the mass profiles of triaxial dark matter halos taken from the SURFS simulation, within scaled radius $0.2<r/r_{\rm max}<3$, are measured with a bias of $\sim-2.4\%$ and a dispersion of $\sim10\%$. The obtained dispersion is not because of Poisson noise due to small particle numbers as it is twice the later. We interpret the dispersion to be due to the inherent nature of the $\Lambda{\rm CDM}$ halos, for example being aspherical and out-of-equilibrium. Hence the dispersion obtained for stellar halos sets a limit of about $12\%$ (after adjusting for random uncertainty) on the accuracy with which the mass profiles of the Milky Way-like galaxies can be reconstructed using the spherical Jeans equation. This limit is independent of the quantity and quality of the observational data. The reason for a non zero bias is not clear, hence its interpretation is not obvious at this stage.

Our nearest large cosmological neighbour, the Andromeda galaxy (M31), is a dynamical system, and an accurate measurement of its total mass is central to our understanding of its assembly history, the life-cycles of its satellite galaxies, and its role in shaping the Local Group environment. Here, we apply a novel approach to determine the dynamical mass of M31 using high velocity Planetary Nebulae (PNe), establishing a hierarchical Bayesian model united with a scheme to capture potential outliers and marginalize over tracers unknown distances. With this, we derive the escape velocity run of M31 as a function of galacto-centric distance, with both parametric and non-parametric approaches. We determine the escape velocity of M31 to be $470\pm{40}$ km s$^{-1}$ at a galacto-centric distance of 15 kpc, and also, derive the total potential of M31, estimating the virial mass and radius of the galaxy to be $0.8\pm{0.1}\times10^{12}\,M_\odot$ and $240\pm{10}$ kpc, respectively. Our M31 mass is on the low-side of the measured range, this supports the lower expected mass of the M31-Milky Way system from the timing and momentum arguments, satisfying the HI constraint on circular velocity between $10\lesssim R/\textrm{kpc}<35$, and agreeing with the stellar mass Tully-Fisher relation. To place these results in a broader context, we compare them to the key predictions of the $\Lambda{\rm CDM}$ cosmological paradigm, including the stellar-mass-halo-mass and the dark matter halo concentration-virial mass correlation, and finding it to be an outlier to this relation.

We present the Synthetic UniveRses For Surveys ({\sc surfs}) simulations, a set of N-body/Hydro simulations of the concordance $\Lambda$ Cold Dark Matter (\LCDM) cosmology. These simulations use Planck cosmology, contain up to 10 billion particles and sample scales & halo masses down to $1~$kpc & $10^8{\rm M}_\odot$. We identify and track haloes from $z=24$ to today using a state-of-the-art 6D halo finder and merger tree builder. We demonstrate that certain properties of halo merger trees are numerically converged for haloes composed of $\gtrsim100$ particles. Haloes smoothly grow in mass, $V_{\rm max}$, with the mass history characterised by $\log M(a)\propto\exp\left[-(a/\beta)^\alpha\right]$ where $a$ is the scale factor, $\alpha(M)\approx0.8$ \& $\beta(M)\approx0.024$, with these parameters decreasing with decreasing halo mass. Subhaloes follow power-law cumulative mass and velocity functions, i.e. $n(>f)\propto f^{-\alpha}$ with $\alpha_{M}=0.83\pm0.01$ and $\alpha_{V_{\rm max}}=2.13\pm0.03$ for mass \& velocity respectively, independent of redshift, as seen in previous studies. The halo-to-halo scatter in amplitude is $0.9$~dex. The number of subhaloes in a halo weakly correlates with a halo's concentration $c$ \& spin $\lambda$:haloes of high $c$ \& low $\lambda$ have $60\%$ more subhaloes than similar mass haloes of low $c$ \& high $\lambda$. High cadence tracking shows subhaloes are dynamic residents, with $25\%$ leaving their host halo momentarily, becoming a backsplash subhalo, and another $20\%$ changing hosts entirely, in agreement with previous studies. In general, subhaloes have elliptical orbits, $e\approx0.6$, with periods of $2.3^{+2.1}_{-1.7}$~Gyrs. Subhaloes lose most of their mass at pericentric passage with mass loss rates of $\sim40\%$~Gyr$^{-1}$. These catalogues will be made publicly available.

The halo mass function, encoding the comoving number density of dark matter halos of a given mass, plays a key role in understanding the formation and evolution of galaxies. As such, it is a key goal of current and future deep optical surveys to constrain the mass function down to mass scales which typically host $L_\star$ galaxies. Motivated by the proven accuracy of Press-Schechter-type mass functions, we introduce a related but purely empirical form consistent with standard formulae to better than 4\% in the medium-mass regime, $10^{10}-10^{13} h^{-1}M_\odot$. In particular, our form consists of 4 parameters, each of which has a simple interpretation, and can be directly related to parameters of the galaxy distribution, such as $L_\star$. Using this form within a hierarchical Bayesian likelihood model, we show how individual mass-measurement errors can be successfully included in a typical analysis whilst accounting for Eddington bias. We apply our form to a question of survey design in the context of a semi-realistic data model, illustrating how it can be used to obtain optimal balance between survey depth and angular coverage for constraints on mass function parameters. Open-source {\tt Python} and {\tt R} codes to apply our new form are provided at \url{http://mrpy.readthedocs.org} and \url{https://cran.r-project.org/web/packages/tggd/index.html} respectively.

We present the results of a VLT MUSE/FORS2 and Spitzer survey of a unique compact lensing cluster CLIO at z = 0.42, discovered through the GAMA survey using spectroscopic redshifts. Compact and massive clusters such as this are understudied, but provide a unique prospective on dark matter distributions and for finding background lensed high-z galaxies. The CLIO cluster was identified for follow up observations due to its almost unique combination of high mass and dark matter halo concentration, as well as having observed lensing arcs from ground based images. Using dual band optical and infra-red imaging from FORS2 and Spitzer, in combination with MUSE optical spectroscopy we identify 89 cluster members and find background sources out to z = 6.49. We describe the physical state of this cluster, finding a strong correlation between environment and galaxy spectral type. Under the assumption of a NFW profile, we measure the total mass of CLIO to be M$_{200} = (4.49 \pm 0.25) \times 10^{14}$ M$_\odot$. We build and present an initial strong-lensing model for this cluster, and measure a relatively low intracluster light (ICL) fraction of 7.21 $\pm$ 1.53% through galaxy profile fitting. Due to its strong potential for lensing background galaxies and its low ICL, the CLIO cluster will be a target for our 110 hour JWST 'Webb Medium-Deep Field' (WMDF) GTO program.

We present a general modified maximum likelihood (MML) method for inferring generative distribution functions from uncertain and biased data. The MML estimator is identical to, but easier and many orders of magnitude faster to compute than the solution of the exact Bayesian hierarchical modelling of all measurement errors. As a key application, this method can accurately recover the mass function (MF) of galaxies, while simultaneously dealing with observational uncertainties (Eddington bias), complex selection functions and unknown cosmic large-scale structure. The MML method is free of binning and natively accounts for small number statistics and non-detections. Its fast implementation in the R-package "dftools" is equally applicable to other objects, such as haloes, groups and clusters, as well as observables other than mass. The formalism readily extends to multi-dimensional distribution functions, e.g. a Choloniewski function for the galaxy mass-angular momentum distribution, also handled by dftools. The code provides uncertainties and covariances for the fitted model parameters and approximate Bayesian evidences. We use numerous mock surveys to illustrate and test the MML method, as well as to emphasise the necessity of accounting for observational uncertainties in MFs of modern galaxy surveys.

We study the scaling relations between the baryonic content and total mass of groups of galaxies, as these systems provide a unique way to examine the role of non-gravitational processes in structure formation. Using Planck and ROSAT data, we conduct detailed comparisons of the stacked thermal Sunyaev-Zel'dovich (tSZ) effect and X-ray scaling relations of galaxy groups found in the the Galaxy And Mass Assembly (GAMA) survey and the BAHAMAS hydrodynamical simulation. We use weak gravitational lensing data from the Kilo Degree Survey (KiDS) to determine the average halo mass of the studied systems. We analyse the simulation in the same way, using realistic weak lensing, X-ray, and tSZ synthetic observations. Furthermore, to keep selection biases under control, we employ exactly the same galaxy selection and group identification procedures to the observations and simulation. Applying this careful comparison, we find that the simulations are in agreement with the observations, particularly with regards to the scaling relations of the lensing and tSZ results. This finding demonstrates that hydrodynamical simulation have reached the level of realism that is required to interpret observational survey data and study the baryon physics within dark matter haloes, where analytical modelling is challenging. Finally, using simulated data, we demonstrate that our observational processing of the X-ray and tSZ signals is free of significant biases. We find that our optical group selection procedure has, however, some room for improvement.

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

An X-ray survey with the XMM-Newton telescope, XMM-XXL, has identified hundreds of galaxy groups and clusters in two 25 deg$^2$ fields. Combining spectroscopic and X-ray observations in one field, we determine how the kinetic energy of galaxies scales with hot gas temperature and also, by imposing prior constraints on the relative energies of galaxies and dark matter, infer a power-law scaling of total mass with temperature. Our goals are: i) to determine parameters of the scaling between galaxy velocity dispersion and X-ray temperature, $T_{\rm 300kpc}$, for the halos hosting XXL-selected clusters, and; ii) to infer the log-mean scaling of total halo mass with temperature, $\langle \ln M_{200} \, | \, T, z \rangle$. We apply an ensemble velocity likelihood to a sample of $> 1500$ spectroscopic redshifts within $132$ spectroscopically confirmed clusters with redshifts $z < 0.6$ to model, $\langle \ln \sigma_{\rm gal}\,|\,T,z\rangle$, where $\sigma_{\rm gal}$ is the velocity dispersion of XXL cluster member galaxies and $T$ is a 300 kpc aperture temperature. To infer total halo mass we use a precise virial relation for massive halos calibrated by N-body simulations along with a single degree of freedom summarizing galaxy velocity bias with respect to dark matter. For the XXL-N cluster sample, we find $\sigma_{\rm gal} \propto T^{0.63\pm0.05}$, a slope significantly steeper than the self-similar expectation of $0.5$. Assuming scale-independent galaxy velocity bias, we infer a mean logarithmic mass at a given X-ray temperature and redshift, $\langle\ln (E(z) M_{200}/10^{14}\,{\rm M}_{\odot})|T,z\rangle=\pi+\alpha \ln(T/T_p )+\beta\ln (E(z)/E(z_p) )$ using pivot values ${\rm k}T_{p}=2.2\,{\rm keV}$ and $z_p=0.25$, with normalization $\pi=0.45\pm0.24$ and slope $\alpha=1.89\pm0.15$. We obtain only weak constraints on redshift evolution, $\beta=-1.29\pm1.14$.

We present a phenomological model of the Cosmic Spectral Energy Distribution (CSED) and the integrated galactic light (IGL) over all cosmic time. This model, based on an earlier model by Driver et al. (2013), attributes the cosmic star formation history to two processes -- firstly, chaotic clump accretion and major mergers, resulting in the early-time formation of bulges and secondly, cold gas accretion, resulting in late-time disc formation. Under the assumption of a Universal Chabrier initial mass function, we combine the Bruzual & Charlot (2003) stellar libraries, the Charlot & Fall (2000) dust attenuation prescription and template spectra for emission by dust and active galactic nuclei to predict the CSED -- pre- and post-dust attenuation -- and the IGL throughout cosmic time. The phenomological model, as constructed, adopts a number of basic axioms and empirical results and has minimal free parameters. We compare the model output, as well as predictions from the semi-analytic model GALFORM to recent estimates of the CSED out to $z=1$. By construction, our empirical model reproduces the full energy output of the Universe from the ultraviolet to the far-infrared extremely well. We use the model to derive predictions of the stellar and dust mass densities, again finding good agreement. We find that GALFORM predicts the CSED for $z < 0.3$ in good agreement with the observations. This agreement becomes increasingly poor towards $z = 1$, when the model CSED is $\sim$50 per cent fainter. The latter is consistent with the model underpredicting the cosmic star formation history. As a consequence, GALFORM predicts a $\sim$30 per cent fainter IGL.

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 fraction of galaxies bound in groups in the nearby Universe is high (50% at z~0). Systematic studies of galaxy properties in groups are important in order to improve our understanding of the evolution of galaxies and of the physical phenomena occurring within this environment. We have built a complete spectrophotometric sample of galaxies within X-ray detected, optically spectroscopically confirmed groups and clusters (G&C), covering a wide range of halo masses at z<= 0.6. In the context of the XXL survey, we analyse a sample of 164 G&C in the XXL-North region (XXL-N), at z <= 0.6, with a wide range of virial masses (1.24 x 10^13 <=M_500 M_sun <= 6.63 x 10^14) and X-ray luminosities ( 2.27 x 10^41 <= L^XXL_500 (erg/s)<= 2.15 x10^44). The G&C are X-ray selected and spectroscopically confirmed. We describe the membership assignment and the spectroscopic completeness analysis, and compute stellar masses. As a first scientific exploitation of the sample, we study the dependence of the galaxy stellar mass function (GSMF) on global environment. We present a spectrophotometric characterisation of the G&C and their galaxies. The final sample contains 132 G&C, 22111 field galaxies and 2225 G&C galaxies with r-band magnitude <20. Of the G&C, 95% have at least three spectroscopic members, and 70% at least ten. The shape of the GSMF seems not to depend on environment (field versus G&C) or X-ray luminosity ( used as a proxy for the virial mass of the system). These results are confirmed by the study of the correlation between mean stellar mass of G&C members and L^XXL_500.We release the spectrophotometric catalogue of galaxies with all the quantities computed in this work. As a first homogeneous census of galaxies within X-ray spectroscopically confirmed G&C at these redshifts, this sample will allow environmental studies of the evolution of galaxy properties.

We introduce a method for modeling disk galaxies designed to take full advantage of data from integral field spectroscopy (IFS). The method fits equilibrium models to simultaneously reproduce the surface brightness, rotation and velocity dispersion profiles of a galaxy. The models are fully self-consistent 6D distribution functions for a galaxy with a Sersic-profile stellar bulge, exponential disk and parametric dark matter halo, generated by an updated version of GalactICS. By creating realistic flux-weighted maps of the kinematic moments (flux, mean velocity and dispersion), we simultaneously fit photometric and spectroscopic data using both maximum-likelihood and Bayesian (MCMC) techniques. We apply the method to a GAMA spiral galaxy (G79635) with kinematics from the SAMI Galaxy Survey and deep $g$- and $r$-band photometry from the VST-KiDS survey, comparing parameter constraints with those from traditional 2D bulge-disk decomposition. Our method returns broadly consistent results for shared parameters, while constraining the mass-to-light ratios of stellar components and reproducing the HI-inferred circular velocity well beyond the limits of the SAMI data. While the method is tailored for fitting integral field kinematic data, it can use other dynamical constraints like central fibre dispersions and \HI circular velocities, and is well-suited for modelling galaxies with a combination of deep imaging and HI and/or optical spectra (resolved or otherwise). Our implementation (MagRite) is computationally efficient and can generate well-resolved models and kinematic maps in under a minute on modern processors.

We investigate measuring the peculiar motions of galaxies up to $z=0.5$ using Type Ia supernovae (SNe Ia) from LSST, and predict the subsequent constraints on the growth rate of structure. We consider two cases. Our first is based on measurements of the volumetric SNe Ia rate and assumes we can obtain spectroscopic redshifts and light curves for varying fractions of objects that are detected pre-peak luminosity by LSST (some of which may be obtained by LSST itself and others which would require additional follow-up). We find that these measurements could produce growth rate constraints at $z<0.5$ that significantly outperform those using Redshift Space Distortions (RSD) with DESI or 4MOST, even though there are $\sim4\times$ fewer objects. For our second case, we use semi-analytic simulations and a prescription for the SNe Ia rate as a function of stellar mass and star formation rate to predict the number of LSST SNe IA whose host redshifts may already have been obtained with the Taipan+WALLABY surveys, or with a future multi-object spectroscopic survey. We find $\sim 18,000$ and $\sim 160,000$ SN Ia with host redshifts for these cases respectively. Whilst this is only a fraction of the total LSST-detected SNe Ia, they could be used to significantly augment and improve the growth rate constraints compared to only RSD. Ultimately, we find that combining LSST SNe Ia with large numbers of galaxy redshifts will provide the most powerful probe of large scale gravity in the $z<0.5$ regime over the coming decades.

We present an analysis of the global stellar populations of galaxies in the SAMI Galaxy Survey. Our sample consists of 1319 galaxies spanning four orders of magnitude in stellar mass and includes all morphologies and environments. We derive luminosity-weighted, single stellar population equivalent stellar ages, metallicities and alpha enhancements from spectra integrated within one effective radius apertures. Variations in galaxy size explain the majority of the scatter in the age--mass and metallicity--mass relations. Stellar populations vary systematically in the plane of galaxy size and stellar mass, such that galaxies with high stellar surface mass density are older, more metal-rich and alpha-enhanced than less dense galaxies. Galaxies with high surface mass densities have a very narrow range of metallicities, however, at fixed mass, the spread in metallicity increases substantially with increasing galaxy size (decreasing density). We identify residual correlations with morphology and environment. At fixed mass and size, galaxies with late-type morphologies, small bulges and low Sersic n are younger than early-type, high n, high bulge-to-total galaxies. Age and metallicity both show small residual correlations with environment; at fixed mass and size, galaxies in denser environments or more massive halos are older and somewhat more metal rich than those in less dense environments. We connect these trends to evolutionary tracks within the size--mass plane.

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

We present evidence for halo assembly bias as a function of geometric environment. By classifying GAMA galaxy groups as residing in voids, sheets, filaments or knots using a tidal tensor method, we find that low-mass haloes that reside in knots are older than haloes of the same mass that reside in voids. This result provides direct support to theories that link strong halo tidal interactions with halo assembly times. The trend with geometric environment is reversed at large halo mass, with haloes in knots being younger than haloes of the same mass in voids. We find a clear signal of halo downsizing - more massive haloes host galaxies that assembled their stars earlier. This overall trend holds independently of geometric environment. We support our analysis with an in-depth exploration of the L-Galaxies semi-analytic model, used here to correlate several galaxy properties with three different definitions of halo formation time. We find a complex relationship between halo formation time and galaxy properties, with significant scatter. We confirm that stellar mass to halo mass ratio, specific star-formation rate and mass-weighed age are reasonable proxies of halo formation time, especially at low halo masses. Instantaneous star-formation rate is a poor indicator at all halo masses. Using the same semi-analytic model, we create mock spectral observations using complex star-formation and chemical enrichment histories, that approximately mimic GAMA's typical signal-to-noise and wavelength range. We use these mocks to assert how well potential proxies of halo formation time may be recovered from GAMA-like spectroscopic data.

Observations with the Chandra X-ray Observatory are used to examine the hot gas properties within a sample of 10 galaxy groups selected from the Galaxy And Mass Assembly survey's optical Friends-of-Friends group catalogue. Our groups have been screened to eliminate spurious and unrelaxed systems, and the effectiveness of this procedure is demonstrated by the detection of intergalactic hot gas in 80 per cent of our sample. However, we find that 9 of the 10 are X-ray underluminous by a mean factor of $\sim$4 compared to typical X-ray-selected samples. Consistent with this, the majority of our groups have gas fractions that are lower and gas entropies somewhat higher than those seen in typical X-ray-selected samples. Two groups, which have high 2{\sigma} lower limits on their gas entropy, are candidates for the population of high-entropy groups predicted by some active galactic nucleus feedback models.

The first and second moments of stellar velocities encode important information about the formation history of the Galactic halo. However, due to the lack of tangential motion and inaccurate distances of the halo stars, the velocity moments in the Galactic halo have largely remained "known unknowns". Fortunately, our off-centric position within the Galaxy allows us to estimate these moments in the galacto-centric frame using the observed radial velocities of the stars alone. We use these velocities coupled with the Hierarchical Bayesian scheme, which allows easy marginalisation over the missing data (the proper-motion, and uncertainty-free distance and line-of-sight velocity), to measure the velocity dispersions, orbital anisotropy ($\beta$) and streaming motion ($v_{\rm rot}$) of the halo main-sequence turn-off (MSTO) and K-giant (KG) stars in the inner stellar halo (r $\lesssim 15$ kpc). We study the metallicity bias in kinematics of the halo stars and observe that the comparatively metal-rich ([Fe/H]$>-1.4$) and the metal-poor ([Fe/H]$\leq - 1.4$) MSTO samples show a clear systematic difference in $v_{\rm rot} \sim 20-40$ km s$^{-1}$, depending on how restrictive the spatial cuts to cull the disk contamination are. The bias is also detected in KG samples but with less certainty. Both MSTO and KG populations suggest that the inner stellar halo of the Galaxy is radially biased i.e. $\sigma_r>\sigma_\theta$ or $\sigma_\phi$ and $\beta \simeq 0.5$. The apparent metallicity contrariety in the rotation velocity among the halo sub-populations supports the co-existence of multiple populations in the galactic halo that may have formed through distinct formation scenarios, i.e. in-situ versus accretion.

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.

The SKA and its pathfinders will enable studies of HI emission at higher redshifts than ever before. In moving beyond the local Universe, this will require the use of cosmologically appropriate formulae that have traditionally been simplified to their low-redshift approximations. In this paper, we summarise some of the most important relations for tracing HI emission in the SKA era, and present an online calculator to assist in the planning and analysis of observations (hifi.icrar.org).

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

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.

We describe the selection of galaxies targeted in eight low redshift clusters (APMCC0917, A168, A4038, EDCC442, A3880, A2399, A119 and A85; $0.029 < z < 0.058$) as part of the Sydney-AAO Multi-Object integral field Spectrograph Galaxy Survey (SAMI-GS). We have conducted a redshift survey of these clusters using the AAOmega multi-object spectrograph on the 3.9m Anglo-Australian Telescope. The redshift survey is used to determine cluster membership and to characterise the dynamical properties of the clusters. In combination with existing data, the survey resulted in 21,257 reliable redshift measurements and 2899 confirmed cluster member galaxies. Our redshift catalogue has a high spectroscopic completeness ($\sim 94\%$) for $r_{\rm petro} \leq 19.4$ and clustercentric distances $R< 2\rm{R}_{200}$. We use the confirmed cluster member positions and redshifts to determine cluster velocity dispersion, $\rm{R}_{200}$, virial and caustic masses, as well as cluster structure. The clusters have virial masses $14.25 \leq {\rm log }({\rm M}_{200}/\rm{M}_{\odot}) \leq 15.19$. The cluster sample exhibits a range of dynamical states, from relatively relaxed-appearing systems, to clusters with strong indications of merger-related substructure. Aperture- and PSF-matched photometry are derived from SDSS and VST/ATLAS imaging and used to estimate stellar masses. These estimates, in combination with the redshifts, are used to define the input target catalogue for the cluster portion of the SAMI-GS. The primary SAMI-GS cluster targets have $R< \rm{R}_{200}$, velocities $|v_{\rm pec}| < 3.5\sigma_{200}$ and stellar masses $9.5 \leq {\rm log(M}^*_{approx}/\rm{M}_{\odot}) \leq 12$. Finally, we give an update on the SAMI-GS progress for the cluster regions.

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.

There are many proposed mechanisms driving the morphological transformation of disk galaxies to elliptical galaxies. In this paper, we determine if the observed transformation in low mass groups can be explained by the merger histories of galaxies. We measured the group mass-morphology relation for groups from the Galaxy and Mass Assembly group catalogue with masses from 10$^{11}$ - 10$^{15}$ M$_{\odot}$. Contrary to previous studies, the fraction of elliptical galaxies in our more complete group sample increases significantly with group mass across the full range of group mass. The elliptical fraction increases at a rate of 0.163$\pm$0.012 per dex of group mass for groups more massive than 10$^{12.5}$ M$_{\odot}$. If we allow for uncertainties in the observed group masses, our results are consistent with a continuous increase in elliptical fraction from group masses as low as 10$^{11}$M$_{\odot}$. We tested if this observed relation is consistent with merger activity using a GADGET-2 dark matter simulation of the galaxy groups. We specified that a simulated galaxy would be transformed to an elliptical morphology either if it experienced a major merger or if its cumulative mass gained from minor mergers exceeded 30 per cent of its final mass. We then calculated a group mass-morphology relation for the simulations. The position and slope of the simulated relation were consistent with the observational relation, with a gradient of 0.184$\pm$0.010 per dex of group mass. These results demonstrate a strong correlation between the frequency of merger events and disk-to-elliptical galaxy transformation in galaxy group environments.

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 ProFit, a new code for Bayesian two-dimensional photometric galaxy profile modelling. ProFit consists of a low-level C++ library (libprofit), accessible via a command-line interface and documented API, along with high-level R (ProFit) and Python (PyProFit) interfaces (available at github.com/ICRAR/ libprofit, github.com/ICRAR/ProFit, and github.com/ICRAR/pyprofit respectively). R ProFit is also available pre-built from CRAN, however this version will be slightly behind the latest GitHub version. libprofit offers fast and accurate two- dimensional integration for a useful number of profiles, including Sersic, Core-Sersic, broken-exponential, Ferrer, Moffat, empirical King, point-source and sky, with a simple mechanism for adding new profiles. We show detailed comparisons between libprofit and GALFIT. libprofit is both faster and more accurate than GALFIT at integrating the ubiquitous Serrsic profile for the most common values of the Serrsic index n (0.5 < n < 8). The high-level fitting code ProFit is tested on a sample of galaxies with both SDSS and deeper KiDS imaging. We find good agreement in the fit parameters, with larger scatter in best-fit parameters from fitting images from different sources (SDSS vs KiDS) than from using different codes (ProFit vs GALFIT). A large suite of Monte Carlo-simulated images are used to assess prospects for automated bulge-disc decomposition with ProFit on SDSS, KiDS and future LSST imaging. We find that the biggest increases in fit quality come from moving from SDSS- to KiDS-quality data, with less significant gains moving from KiDS to LSST.

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.

We present a quantitative measurement of the amount of clustering present in the inner $\sim30$ kpc of the stellar halo of the Andromeda galaxy (M31). For this we analyse the angular positions and radial velocities of the carefully selected Planetary Nebulae (PNe) in the M31 stellar halo. We study the cumulative distribution of pair-wise distances in angular position and line of sight velocity space, and find that the M31 stellar halo contains substantially more stars in the form of close pairs as compared to that of a featureless smooth halo. In comparison to a smoothed/scrambled distribution we estimate that the clustering excess in the M31 inner halo is roughly $40\%$ at maximum and on average $\sim 20\%$. Importantly, comparing against the 11 stellar halo models of \cite{2005ApJ...635..931B}, which were simulated within the context of the $\Lambda{\rm CDM}$ cosmological paradigm, we find that the amount of substructures in the M31 stellar halo closely resembles that of a typical $\Lambda{\rm CDM}$ halo.

We present a consistent total flux catalogue for a $\sim$1 deg$^2$ subset of the COSMOS region (R.A. $\in [149.55\degr, 150.65\degr]$, DEC $\in [1.80\degr, 2.73\degr]$) with near-complete coverage in 38 bands from the far-ultraviolet to the far-infrared. We produce aperture matched photometry for 128,304 objects with i < 24.5 in a manner that is equivalent to the Wright et al. (2016) catalogue from the low-redshift (z < 0.4) Galaxy and Mass Assembly (GAMA) survey. This catalogue is based on publicly available imaging from GALEX, CFHT, Subaru, VISTA, Spitzer and Herschel, contains a robust total flux measurement or upper limit for every object in every waveband and complements our re-reduction of publicly available spectra in the same region. We perform a number of consistency checks, demonstrating that our catalogue is comparable to existing data sets, including the recent COSMOS2015 catalogue (Laigle et al. 2016). We also release an updated Davies et al. (2015) spectroscopic catalogue that folds in new spectroscopic and photometric redshift data sets. The catalogues are available for download at http://cutout.icrar.org/G10/dataRelease.php. Our analysis is optimised for both panchromatic analysis over the full wavelength range and for direct comparison to GAMA, thus permitting measurements of galaxy evolution for 0 < z < 1 while minimising the systematic error resulting from disparate data reduction methods.

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 build on a recent photometric decomposition analysis of 7506 Galaxy and Mass Assembly (GAMA) survey galaxies to derive stellar mass function fits to individual spheroid and disk component populations down to a lower mass limit of log(M_*/M_sun)= 8. We find that the spheroid/disk mass distributions for individual galaxy morphological types are well described by single Schechter function forms. We derive estimates of the total stellar mass densities in spheroids (rho_spheroid = 1.24+/-0.49 * 10^8 M_sun Mpc^-3 h_0.7) and disks (rho_disk = 1.20+/-0.45 * 10^8 M_sun Mpc^-3 h_0.7), which translates to approximately 50% of the local stellar mass density in spheroids and 48% in disks. The remaining stellar mass is found in the dwarf "little blue spheroid" class, which is not obviously similar in structure to either classical spheroid or disk populations. We also examine the variation of component mass ratios across galaxy mass and group halo mass regimes, finding the transition from spheroid to disk mass dominance occurs near galaxy stellar mass ~10^11 M_sun and group halo mass ~10^12.5 M_sun/h. We further quantify the variation in spheroid-to-total mass ratio with group halo mass for central and satellite populations as well as the radial variation of this ratio within groups.

This paper is the second in a pair of articles presenting data release 1 (DR1) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), the largest single open-time key project carried out with the Herschel Space Observatory. The H-ATLAS is a wide-area imaging survey carried out in five photometric bands at 100, 160, 250, 350 and 500$\mu$m covering a total area of 600deg$^2$. In this paper we describe the identification of optical counterparts to submillimetre sources in DR1, comprising an area of 161 deg$^2$ over three equatorial fields of roughly 12$^\circ$x4.5$^\circ$ centred at 9$^h$, 12$^h$ and 14.5$^h$ respectively. Of all the H-ATLAS fields, the equatorial regions benefit from the greatest overlap with current multi-wavelength surveys spanning ultraviolet (UV) to mid-infrared regimes, as well as extensive spectroscopic coverage. We use a likelihood-ratio technique to identify SDSS counterparts at r<22.4 for 250-$\mu$m-selected sources detected at $\geq$ 4$\sigma$ ($\approx$28mJy). We find `reliable' counterparts (reliability R$\geq$0.8) for 44,835 sources (39 per cent), with an estimated completeness of 73.0 per cent and contamination rate of 4.7 per cent. Using redshifts and multi-wavelength photometry from GAMA and other public catalogues, we show that H-ATLAS-selected galaxies at $z<0.5$ span a wide range of optical colours, total infrared (IR) luminosities, and IR/UV ratios, with no strong disposition towards mid-IR-classified AGN in comparison with optical selection. The data described herein, together with all maps and catalogues described in the companion paper (Valiante et al. 2016), are available from the H-ATLAS website at www.h-atlas.org.

We perform automated bulge + disc decomposition on a sample of $\sim$7500 galaxies from the Galaxy And Mass Assembly (GAMA) survey in the redshift range of 0.002$<$z$<$0.06 using SIGMA, a wrapper around GALFIT3. To achieve robust profile measurements we use a novel approach of repeatedly fitting the galaxies, varying the input parameters to sample a large fraction of the input parameter space. Using this method we reduce the catastrophic failure rate significantly and verify the confidence in the fit independently of $\chi^2$. Additionally, using the median of the final fitting values and the 16$^{th}$ and 84$^{th}$ percentile produces more realistic error estimates than those provided by GALFIT, which are known to be underestimated. We use the results of our decompositions to analyse the stellar mass - half-light radius relations of bulges, discs and spheroids. We further investigate the association of components with a parent disc or elliptical relation to provide definite $z=0$ disc and spheroid $\mathcal{M_\star}-R_{\rm e}$ relations. We conclude by comparing our local disc and spheroid $\mathcal{M_\star}-R_{\rm e}$ to simulated data from EAGLE and high redshift data from CANDELS-UDS. We show the potential of using the mass-size relation to study galaxy evolution in both cases but caution that for a fair comparison all data sets need to be processed and analysed in the same manner.

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.

There is strong evidence that cosmological N-body simulations dominated by Warm Dark Matter (WDM) contain spurious or unphysical haloes, most readily apparent as regularly spaced low-mass haloes strung along filaments. We show that spurious haloes are a feature of traditional N-body simulations of cosmological structure formation models, including WDM and Cold Dark Matter (CDM) models, in which gravitational collapse proceeds in an initially anisotropic fashion, and arises naturally as a consequence of discreteness-driven relaxation. We demonstrate this using controlled N-body simulations of plane-symmetric collapse and show that spurious haloes are seeded at shell crossing by localised velocity perturbations induced by the discrete nature of the density field, and that their characteristic separation should be approximately the mean inter-particle separation of the N-body simulation, which is fixed by the mass resolution within the volume. Using cosmological N-body simulations in which particles are split into two collisionless components of fixed mass ratio, we find that the spatial distribution of the two components show signatures of discreteness-driven relaxation in their spatial distribution on both large and small scales. Adopting a spline kernel gravitational softening that is of order the comoving mean inter-particle separation helps to suppress the effect of discreteness-driven relaxation, but cannot eliminate it completely. These results provide further motivation for recent developments of new algorithms, which include, for example, revisions of the traditional N-body approach by means of spatially adaptive anistropric gravitational softenings or explicit solutions for the evolution of dark matter in phase space.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

This short document is intended as a companion and introduction to the Detailed Science Case (DSC) for the Maunakea Spectroscopic Explorer. It provides a concise summary of the essential characteristics of MSE from the perspective of the international astronomical community. MSE is a wide field telescope (1.5 square degree field of view) with an aperture of 11.25m. It is dedicated to multi-object spectroscopy at several different spectral resolutions in the range R ~ 2500 - 40000 over a broad wavelength range (0.36 - 1.8{\mu}m). MSE will enable transformational science in areas as diverse as exoplanetary host characterization; stellar monitoring campaigns; tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of the distant Galaxy; connecting galaxies to the large scale structure of the Universe; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars. MSE is the largest ground based optical and near infrared telescope in its class, and it will occupy a unique and critical role in the emerging network of astronomical facilities active in the 2020s. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, eROSITA, SKA, and WFIRST, and is an ideal feeder facility for E-ELT, TMT and GMT.

Deep observations of galaxies reveal faint extended stellar components (hereafter ESCs) of streams, shells, and halos. These are a natural prediction of hierarchical galaxy formation, as accreted satellite galaxies are tidally disrupted by their host. We investigate whether or not global properties of the ESC could be used to test of dark matter, reasoning that they should be sensitive to the abundance of low-mass satellites, and therefore the underlying dark matter model. Using cosmological simulations of galaxy formation in the favoured Cold Dark Matter (CDM) and Warm Dark Matter (WDM) models ($m_{\rm WDM}$=0.5,1,2 keV/$c^2$), which suppress the abundance of low-mass satellites, we find that the kinematics and orbital structure of the ESC is consistent across models. However, we find striking differences in its spatial structure, as anticipated -- a factor of $\sim$10 drop in spherically averaged mass density between $\sim$10% and $\sim$75% of the virial radius in the more extreme WDM runs ($m_{\rm WDM}$=0.5, 1 keV/$c^2$) relative to the CDM run. These differences are consistent with the mass assembly histories of the different components, and are present across redshifts. However, even the least discrepant of the WDM models is incompatible with current observational limits on $m_{\rm WDM}$. Importantly, the differences we observe when varying the underlying dark matter are comparable to the galaxy-to-galaxy variation we expect within a fixed dark matter model. This suggests that it will be challenging to place limits on dark matter using only the unresolved spatial structure of the the ESC.

We combine wide and deep galaxy number-count data from GAMA, COSMOS/G10, HST ERS, HST UVUDF and various near-, mid- and far- IR datasets from ESO, Spitzer and Herschel. The combined data range from the far-UV (0.15microns) to far-IR (500microns), and in all cases the contribution to the integrated galaxy light (IGL) of successively fainter galaxies converges. Using a simple spline fit, we derive the IGL and the extrapolated-IGL in all bands. We argue undetected low surface brightness galaxies and intra-cluster/group light is modest, and that our extrapolated-IGL measurements are an accurate representation of the extra-galactic background light. Our data agree with most earlier IGL estimates and with direct measurements in the far-IR, but disagree strongly with direct estimates in the optical. Close agreement between our results and recent very high-energy experiments (H.E.S.S. and MAGIC), suggest that there may be an additional foreground affecting the direct estimates. The most likely culprit could be the adopted Zodiacal light model. Finally we use a modified version of the two-component model to integrate the EBL and obtain measurements of the Cosmic Optical Background (COB) and Cosmic Infrared Background (CIB) of: $24^{+4}_{-4}$nWm$^{-2}$sr$^{-1}$ and $26^{+5}_{-5}$nWm$^{-2}$sr$^{-1}$ respectively (48:52%). Over the next decade, upcoming space missions such as Euclid and WFIRST, have the capacity to reduce the COB error to $<1\%$, at which point comparisons to the very high energy data could have the potential to provide a direct detection and measurement of the reionisation field.

Galaxies and their dark matter haloes are part of a complex network of mass structures, collectively called the cosmic web. Using the tidal tensor prescription these structures can be classified into four cosmic environments: voids, sheets, filaments and knots. As the cosmic web may influence the formation and evolution of dark matter haloes and the galaxies they host, we aim to study the effect of these cosmic environments on the average mass of galactic haloes. To this end we measure the galaxy-galaxy lensing profile of 91,195 galaxies, within 0.039 < z < 0.263, from the spectroscopic Galaxy And Mass Assembly (GAMA) survey, using ~100 square degrees of overlapping data from the Kilo-Degree Survey (KiDS). In each of the four cosmic environments we model the contributions from group centrals, satellites and neighbouring groups to the stacked galaxy-galaxy lensing profiles. After correcting the lens samples for differences in the stellar mass distribution, we find no dependence of the average halo mass of central galaxies on their cosmic environment. We do find a significant increase in the average contribution of neighbouring groups to the lensing profile in increasingly dense cosmic environments. We show, however, that the observed effect can be entirely attributed to the galaxy density at much smaller scales (within 4 Mpc/h), which is correlated with the density of the cosmic environments. Within our current uncertainties we find no direct dependence of galaxy halo mass on their cosmic environment.

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 report an expanded sample of visual morphological classifications from the Galaxy and Mass Assembly (GAMA) survey phase two, which now includes 7,556 objects (previously 3,727 in phase one). We define a local (z <0.06) sample and classify galaxies into E, S0-Sa, SB0-SBa, Sab-Scd, SBab-SBcd, Sd-Irr, and "little blue spheroid" types. Using these updated classifications, we derive stellar mass function fits to individual galaxy populations divided both by morphological class and more general spheroid- or disk-dominated categories with a lower mass limit of log(Mstar/Msun) = 8 (one dex below earlier morphological mass function determinations). We find that all individual morphological classes and the combined spheroid-/bulge-dominated classes are well described by single Schechter stellar mass function forms. We find that the total stellar mass densities for individual galaxy populations and for the entire galaxy population are bounded within our stellar mass limits and derive an estimated total stellar mass density of rho_star = 2.5 x 10^8 Msun Mpc^-3 h_0.7, which corresponds to an approximately 4% fraction of baryons found in stars. The mass contributions to this total stellar mass density by galaxies that are dominated by spheroidal components (E and S0-Sa classes) and by disk components (Sab-Scd and Sd-Irr classes) are approximately 70% and 30%, respectively.

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/

From two very simple axioms: (1) that AGN activity traces spheroid formation, and (2) that the cosmic star-formation history is dominated by spheroid formation at high redshift, we derive simple expressions for the star-formation histories of spheroids and discs, and their implied metal enrichment histories. Adopting a Baldry-Glazebrook initial mass function we use these relations and apply PEGASE.2 to predict the z=0 cosmic spectral energy distributions (CSEDs) of spheroids and discs. The model predictions compare favourably to the dust-corrected CSED recently reported by the Galaxy And Mass Assembly (GAMA) team from the FUV through to the K band. The model also provides a reasonable fit to the total stellar mass contained within spheroid and disc structures as recently reported by the Millennium Galaxy Catalogue team. Three interesting inferences can be made following our axioms: (1) there is a transition redshift at z ~ 1.7 at which point the Universe switches from what we refer to as "hot mode evolution" (i.e., spheroid formation/growth via mergers and/or collapse) to what we term "cold mode evolution" (i.e., disc formation/growth via gas infall and minor mergers); (2) there is little or no need for any pre-enrichment prior to the main phase of star-formation; (3) in the present Universe mass-loss is fairly evenly balanced with star-formation holding the integrated stellar mass density close to a constant value. The model provides a simple prediction of the energy output from spheroid and disc projenitors, the build-up of spheroid and disc mass, and the mean metallicity enrichment of the Universe.

Both theoretical predictions and observations of the very nearby Universe suggest that low-mass galaxies (log$_{10}$[M$_{*}$/M$_{\odot}$]<9.5) are likely to remain star-forming unless they are affected by their local environment. To test this premise, we compare and contrast the local environment of both passive and star-forming galaxies as a function of stellar mass, using the Galaxy and Mass Assembly survey. We find that passive fractions are higher in both interacting pair and group galaxies than the field at all stellar masses, and that this effect is most apparent in the lowest mass galaxies. We also find that essentially all passive log$_{10}$[M$_{*}$/M$_{\odot}$]<8.5 galaxies are found in pair/group environments, suggesting that local interactions with a more massive neighbour cause them to cease forming new stars. We find that the effects of immediate environment (local galaxy-galaxy interactions) in forming passive systems increases with decreasing stellar mass, and highlight that this is potentially due to increasing interaction timescales giving sufficient time for the galaxy to become passive via starvation. We then present a simplistic model to test this premise, and show that given our speculative assumptions, it is consistent with our observed results.

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 use the first 100 sq. deg. of overlap between the Kilo-Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey to determine the galaxy halo mass of ~10,000 spectroscopically-confirmed satellite galaxies in massive ($M > 10^{13}h^{-1}{\rm M}_\odot$) galaxy groups. Separating the sample as a function of projected distance to the group centre, we jointly model the satellites and their host groups with Navarro-Frenk-White (NFW) density profiles, fully accounting for the data covariance. The probed satellite galaxies in these groups have total masses $\log M_{\rm sub} /(h^{-1}{\rm M}_\odot) \approx 11.7 - 12.2$ consistent across group-centric distance within the errorbars. Given their typical stellar masses, $\log M_{\rm \star,sat}/(h^{-2}{\rm M}_\odot) \sim 10.5$, such total masses imply stellar mass fractions of $M_{\rm \star,sat} /M_{\rm sub} \approx 0.04 h^{-1}$ . The average subhalo hosting these satellite galaxies has a mass $M_{\rm sub} \sim 0.015M_{\rm host}$ independent of host halo mass, in broad agreement with the expectations of structure formation in a $\Lambda$CDM universe.

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.

Astronomical data is often uncertain with errors that are heteroscedastic (different for each data point) and covariant between different dimensions. Assuming that a set of D-dimensional data points can be described by a (D - 1)-dimensional plane with intrinsic scatter, we derive the general likelihood function to be maximised to recover the best fitting model. Alongside the mathematical description, we also release the hyper-fit package for the R statistical language (github.com/asgr/hyper.fit) and a user-friendly web interface for online fitting (hyperfit.icrar.org). The hyper-fit package offers access to a large number of fitting routines, includes visualisation tools, and is fully documented in an extensive user manual. Most of the hyper-fit functionality is accessible via the web interface. In this paper we include applications to toy examples and to real astronomical data from the literature: the mass-size, Tully-Fisher, Fundamental Plane, and mass-spin-morphology relations. In most cases the hyper-fit solutions are in good agreement with published values, but uncover more information regarding the fitted model.

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 "Wide Area VISTA Extra-galactic Survey" (WAVES) is a 4MOST Consortium Design Reference Survey which will use the VISTA/4MOST facility to spectroscopically survey ~2million galaxies to $r_{\rm AB} < 22$ mag. WAVES consists of two interlocking galaxy surveys ("WAVES-Deep" and "WAVES-Wide"), providing the next two steps beyond the highly successful 1M galaxy Sloan Digital Sky Survey and the 250k Galaxy And Mass Assembly survey. WAVES will enable an unprecedented study of the distribution and evolution of mass, energy, and structures extending from 1-kpc dwarf galaxies in the local void to the morphologies of 200-Mpc filaments at $z\sim1$. A key aim of both surveys will be to compare comprehensive empirical observations of the spatial properties of galaxies, groups, and filaments, against state-of-the-art numerical simulations to distinguish between various Dark Matter models.

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 properties of the first 250 $\mu$m blind sample of nearby galaxies (15 < D < 46 Mpc) containing 42 objects from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Herschel's sensitivity probes the faint end of the dust luminosity function for the first time, spanning a range of stellar mass (7.4 < log$_{10}$ M$_{\star}$ < 11.3 M$_{\odot}$), star formation activity (-11.8 < log$_{10}$ SSFR < -8.9 yr$^{-1}$), gas fraction (3-96 per cent), and colour (0.6 < FUV-Ks < 7.0 mag). The median cold dust temperature is 14.6 K, colder than in the Herschel Reference Survey (18.5 K) and Planck Early Release Compact Source Catalogue 17.7 K. The mean dust-to-stellar mass ratio in our sample is higher than these surveys by factors of 3.7 and 1.8, with a dust mass volume density of (3.7 $\pm$ 0.7) x 10$^{5}$ M$_{\odot}$ Mpc$^{-3}$. Counter-intuitively, we find that the more dust rich a galaxy, the lower its UV attenuation. Over half of our dust-selected sample are very blue in FUV-Ks colour, with irregular and/or highly flocculent morphology, these galaxies account for only 6 per cent of the sample's stellar mass but contain over 35 per cent of the dust mass. They are the most actively star forming galaxies in the sample, with the highest gas fractions and lowest UV attenuation. They also appear to be in an early stage of converting their gas into stars, providing valuable insights into the chemical evolution of young galaxies.

The Herschel Space Observatory has had a tremendous impact on the study of extragalactic dust. Specifically, early-type galaxies (ETG) have been the focus of several studies. In this paper we combine results from two Herschel studies - a Virgo cluster study HeViCS and a broader, low-redshift H-ATLAS/GAMA study - and contrast the dust and associated properties for similar mass galaxies. This comparison is motivated by differences in results exhibited between multiple Herschel studies of early-type galaxies. A comparison between consistent modified blackbody derived dust mass is carried out, revealing strong differences between the two samples in both dust mass and dust-to-stellar mass ratio. In particular, the HeViCS sample lacks massive ETG with as high a specific dust content as found in H-ATLAS. This is most likely connected with the difference in environment for the two samples. We calculate nearest neighbour environment densities in a consistent way, showing that H-ATLAS ETG occupy sparser regions of the local Universe, whereas HeViCS ETG occupy dense regions. This is also true for ETG that are not Herschel-detected but are in the Virgo and GAMA parent samples. Spectral energy distributions are fit to the panchromatic data. From these we find that in H-ATLAS the specific star formation rate anticorrelates with stellar mass and reaches values as high as in our Galaxy. On the other hand HeViCS ETG appear to have little star formation. Based on the trends found here, H-ATLAS ETG are thought to have more extended star formation histories and a younger stellar population than HeViCS ETG.

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 present a simple method of generating approximate synthetic halo catalogues: HALOGEN. This method uses a combination of $2^{nd}$-order Lagrangian Perturbation Theory (2LPT) in order to generate the large-scale matter distribution, analytical mass functions to generate halo masses, and a single-parameter stochastic model for halo bias to position haloes. HALOGEN represents a simplification of similar recently published methods. Our method is constrained to recover the 2-point function at intermediate ($10Mpc/h<r<50Mpc/h$) scales, which we show is successful to within 2 per cent. Larger scales ($\sim100Mpc/h$) are reproduced to within 15 per cent. We compare several other statistics (e.g. power spectrum, point distribution function, redshift space distortions) with results from N-Body simulations to determine the validity of our method for different purposes. One of the benefits of HALOGEN is its flexibility, and we demonstrate this by showing how it can be adapted to varying cosmologies and simulation specifications. A driving motivation for the development of such approximate schemes is the need to compute covariance matrices and study the systematic errors for large galaxy surveys, which requires thousands of simulated realisations. We discuss the applicability of our method in this context, and conclude that it is well suited to mass production of appropriate halo catalogues. The code is publicly available at https://github.com/savila/halogen

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

The science achievable with SKA HI surveys will be greatly increased through the combination of HI data with that at other wavelengths. These multiwavelength datasets will enable studies to move beyond an understanding of HI gas in isolation to instead understand HI as an integral part of the highly complex baryonic processes that drive galaxy evolution. As they evolve, galaxies experience a host of environmental and feedback influences, many of which can radically impact their gas content. Important processes include: accretion (hot and cold mode, mergers), depletion (star formation, galactic winds, AGN), phase changes (ionised/atomic/molecular), and environmental effects (ram pressure stripping, tidal effects, strangulation). Governing all of these to various extents is the underlying dark matter distribution. In turn, the result of these processes can significantly alter the baryonic states in which material is finally observed (stellar populations, dust, chemistry) and its morphology (galaxy type, bulge/disk ratio, bars, warps, radial profile). To fully understand the evolution of HI and the role it plays in galactic evolution requires the ability to quantify each of these separate processes, and hence to coordinate SKA HI surveys with extensive multi-band photometric and spectroscopic campaigns. In addition, multiwavelength data is essential for statistical methods of HI analysis such as HI stacking and intensity mapping cross-correlations. In this chapter, we examine some of the principal science motivations for acquiring multiwavelength data to match that from the extragalactic SKA HI surveys, and review the currently planned capacity to achieve this (eg. LSST, Euclid, W-FIRST, SPICA, ALMA, and 4MOST).

We discuss the construction of the Galaxy And Mass Assembly (GAMA) 10h region (G10) using publicly available data in the Cosmic Evolution Survey region (COSMOS) in order to extend the GAMA survey to z~1 in a single ~1 deg^2 field. In order to obtain the maximum number of high precision spectroscopic redshifts we re-reduce all archival zCOSMOS-bright data and use the GAMA automatic cross-correlation redshift fitting code autoz. We use all available redshift information (autoz, zCOSMOS-bright 10k, PRIMUS, VVDS, SDSS and photometric redshifts) to calculate robust best-fit redshifts for all galaxies and visually inspect all 1D and 2D spectra to obtain 16,583 robust redshifts in the full COSMOS region. We then define the G10 region to be the central ~1deg^2 of COSMOS, which has relatively high spectroscopic completeness, and encompasses the CHILES VLA region. We define a combined r < 23.0mag & i < 22.0mag G10 sample (selected to have the highest bijective overlap) with which to perform future analysis, containing 9,861 sources with reliable high precision VLT-VIMOS spectra. All tables, spectra and imaging are available at: http://ict.icrar.org/cutout/G10 .

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.

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

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 use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.

Comparison of observed satellite galaxies of the Milky Way (hereafter MW) with dark matter subhaloes in cosmological $N$-body simulations of MW-mass haloes suggest that such subhaloes, if they exist, are occupied by satellites in a stochastic fashion. We examine how inefficient massive star formation and associated supernova feedback in high-redshift progenitors of present-day low-mass subhaloes might contribute to this stochasticity. Using a Monte Carlo approach to follow the assembly histories of present-day low-mass haloes with $10^7 \lesssim M \leq 10^{10}$ ${\rm M}_{\odot}$, we identify when cooling and star formation is likely to proceed, and observe that haloes with present-day masses $\lesssim 10^9 {\rm M}_{\odot}$ never grow sufficiently massive to support atomic hydrogen line cooling. Noting that the star formation timescale decreases sharply with stellar mass as $t_{\rm PMS} \propto m_{\ast}^{-2.5}$, we argue that, should the conditions for high mass star formation arise in low-mass haloes, the ensuing supernovae are likely to disrupt ongoing lower-mass star formation and unbind gas within the halo. This potentially star-forming gas is unlikely to be replenished in lower mass haloes because of, e.g. cosmological reionization, and so we expect galaxy formation to be stymied in a manner that depends on host halo assembly history and the efficiency and timing of star formation in proto-galaxies, which we illustrate using a Monte Carlo model. Based on these simple physical arguments, we assert that stochasticity of star formation and feedback is an essential but overlooked ingredient in modelling galaxy formation on the smallest scales.

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.

We present a comparison of the physical properties of a rest-frame $250\mu$m selected sample of massive, dusty galaxies from $0<z<5.3$. Our sample comprises 29 high-redshift submillimetre galaxies (SMGs) from the literature, and 843 dusty galaxies at $z<0.5$ from the Herschel-ATLAS, selected to have a similar stellar mass to the SMGs. The $z>1$ SMGs have an average SFR of $390^{+80}_{-70}\,$M$_\odot$yr$^{-1}$ which is 120 times that of the low-redshift sample matched in stellar mass to the SMGs (SFR$=3.3\pm{0.2}$ M$_\odot$yr$^{-1}$). The SMGs harbour a substantial mass of dust ($1.2^{+0.3}_{-0.2}\times{10}^9\,$M$_\odot$), compared to $(1.6\pm0.1)\times{10}^8\,$M$_\odot$ for low-redshift dusty galaxies. At low redshifts the dust luminosity is dominated by the diffuse ISM, whereas a large fraction of the dust luminosity in SMGs originates from star-forming regions. At the same dust mass SMGs are offset towards a higher SFR compared to the low-redshift H-ATLAS galaxies. This is not only due to the higher gas fraction in SMGs but also because they are undergoing a more efficient mode of star formation, which is consistent with their bursty star-formation histories. The offset in SFR between SMGs and low-redshift galaxies is similar to that found in CO studies, suggesting that dust mass is as good a tracer of molecular gas as CO.

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.

The Galaxy And Mass Assembly (GAMA) survey furnishes a deep redshift catalog that, when combined with the Wide-field Infrared Explorer ($WISE$), allows us to explore for the first time the mid-infrared properties of $> 110, 000$ galaxies over 120 deg$^2$ to $z\simeq 0.5$. In this paper we detail the procedure for producing the matched GAMA-$WISE$ catalog for the G12 and G15 fields, in particular characterising and measuring resolved sources; the complete catalogs for all three GAMA equatorial fields will be made available through the GAMA public releases. The wealth of multiwavelength photometry and optical spectroscopy allows us to explore empirical relations between optically determined stellar mass (derived from synthetic stellar population models) and 3.4micron and 4.6micron WISE measurements. Similarly dust-corrected Halpha-derived star formation rates can be compared to 12micron and 22micron luminosities to quantify correlations that can be applied to large samples to $z<0.5$. To illustrate the applications of these relations, we use the 12micron star formation prescription to investigate the behavior of specific star formation within the GAMA-WISE sample and underscore the ability of WISE to detect star-forming systems at $z\sim0.5$. Within galaxy groups (determined by a sophisticated friends-of-friends scheme), results suggest that galaxies with a neighbor within 100$\,h^{-1} $kpc have, on average, lower specific star formation rates than typical GAMA galaxies with the same stellar mass.

(Abridged) We present a non-parametric cell-based method of selecting highly pure and largely complete samples of spiral galaxies using photometric and structural parameters as provided by standard photometric pipelines and simple shape fitting algorithms, demonstrably superior to commonly used proxies. Furthermore, we find structural parameters derived using passbands longwards of the $g$ band and linked to older stellar populations, especially the stellar mass surface density $\mu_*$ and the $r$ band effective radius $r_e$, to perform at least equally well as parameters more traditionally linked to the identification of spirals by means of their young stellar populations. In particular the distinct bimodality in the parameter $\mu_*$, consistent with expectations of different evolutionary paths for spirals and ellipticals, represents an often overlooked yet powerful parameter in differentiating between spiral and non-spiral/elliptical galaxies. We investigate the intrinsic specific star-formation rate - stellar mass relation ($\psi_* - M_*$) for a morphologically defined volume limited sample of local universe spiral galaxies, defined using the cell-based method with an appropriate parameter combination. The relation is found to be well described by $\psi_* \propto M_*^{-0.5}$ over the range of $10^{9.5} M_{\odot} \le M_* \le 10^{11} M_{\odot}$ with a mean interquartile range of $0.4\,$dex. This is somewhat steeper than previous determinations based on colour-selected samples of star-forming galaxies, primarily due to the inclusion in the sample of red quiescent disks.

The 2MASS, UKIDSS-LAS and VISTA VIKING surveys have all now observed the GAMA 9hr region in the $K_s$ band. Here we compare the detection rates, photometry, basic size measurements, and single-component GALFIT structural measurements for a sample of 37,591 galaxies. We explore the sensitivity limits where the data agree for a variety of issues including: detection, star-galaxy separation, photometric measurements, size and ellipticity measurements, and S\'{e}rsic measurements. We find that 2MASS fails to detect at least 20% of the galaxy population within all magnitude bins, however for those that are detected we find photometry is robust ($\pm$ 0.2 mag) to 14.7 AB mag and star-galaxy separation to 14.8 AB mag. For UKIDSS-LAS we find incompleteness starts to enter at a flux limit of 18.9 AB mag, star-galaxy separation is robust to 16.3 AB mag and structural measurements are robust to 17.7 AB mag. VISTA VIKING data is complete to approximately 20.0 AB mag and structural measurements appear robust to 18.8 AB mag.

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

The dark matter halo mass function (HMF) is a characteristic property of cosmological structure formation models, quantifying the number density of dark matter haloes per unit mass in the Universe. A key goal of current and planned large galaxy surveys is to measure the HMF and to use it to test theories of dark matter and dark energy. We present a new web application for calculating the HMF -- the frontend HMFcalc and the engine hmf. HMFcalc has been designed to be flexible, efficient and easy to use, providing observational and theoretical astronomers alike with the means to explore standard functional forms of the HMF or to tailor their own. We outline the theoretical background needed to compute the HMF, we show how it has been implemented in hmf, and finally we provide worked examples that illustrate HMFcalc's versatility as an analysis tool.

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.

The parameters governing the standard \Lambda Cold Dark Matter cosmological model have been constrained with unprecedented accuracy by precise measurements of the cosmic microwave background by the Wilkinson Microwave Anisotropy Probe (WMAP) and Planck satellites. Each new data release has refined further our knowledge of quantities -- such as the matter density parameter \Omega_M -- that are imprinted on the dark matter halo mass function (HMF), a powerful probe of dark matter and dark energy models. In this letter we trace how changes in the cosmological parameters over the last decade have influenced uncertainty in our knowledge of the HMF. We show that this uncertainty has reduced significantly since the 3rd WMAP data release, but the rate of this reduction is slowing. This is limited by uncertainty in the normalisation \sigma_8, whose influence is most pronounced at the high mass end of the mass function. Interestingly, we find that the accuracy with which we can constrain the HMF in terms of the cosmological parameters has now reached the point at which it is comparable to the scatter in HMF fitting functions. This suggests that the power of the HMF as a precision probe of dark matter and dark energy hinges on more accurate determination of the theoretical HMF. Finally, we assess prospects of using the HMF to differentiate between Cold and Warm Dark Matter models based on ongoing improvements in measurements of \Omega_M, and we comment briefly on optimal survey strategies for constraining dark matter and dark energy models using the HMF.

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

In this paper, we demonstrate a new method for fitting galaxy profiles which makes use of the full multi-wavelength data provided by modern large optical-near-infrared imaging surveys. We present a new version of GALAPAGOS, which utilises a recently-developed multi-wavelength version of GALFIT, and enables the automated measurement of wavelength dependent S\'ersic profile parameters for very large samples of galaxies. Our new technique is extensively tested to assess the reliability of both pieces of software, GALFIT and GALAPAGOS on both real ugrizY JHK imaging data from the GAMA survey and simulated data made to the same specifications. We find that fitting galaxy light profiles with multi-wavelength data increases the stability and accuracy of the measured parameters, and hence produces more complete and meaningful multi-wavelength photometry than has been available previously. The improvement is particularly significant for magnitudes in low S/N bands and for structural parameters like half-light radius re and S\'ersic index n for which a prior is used by constraining these parameters to a polynomial as a function of wavelength. This allows the fitting routines to push the magnitude of galaxies for which sensible values can be derived to fainter limits. The technique utilises a smooth transition of galaxy parameters with wavelength, creating more physically meaningful transitions than single-band fitting and allows accurate interpolation between passbands, perfect for derivation of rest-frame values.

We have constructed a sample of radio-loud objects with optical spectroscopy from the Galaxy and Mass Assembly (GAMA) project over the Herschel-ATLAS Phase 1 fields. Classifying the radio sources in terms of their optical spectra, we find that strong-emission-line sources (`high-excitation radio galaxies') have, on average, a factor ~4 higher 250-micron Herschel luminosity than weak-line (`low-excitation') radio galaxies and are also more luminous than magnitude-matched radio-quiet galaxies at the same redshift. Using all five H-ATLAS bands, we show that this difference in luminosity between the emission-line classes arises mostly from a difference in the average dust temperature; strong-emission-line sources tend to have comparable dust masses to, but higher dust temperatures than, radio galaxies with weak emission lines. We interpret this as showing that radio galaxies with strong nuclear emission lines are much more likely to be associated with star formation in their host galaxy, although there is certainly not a one-to-one relationship between star formation and strong-line AGN activity. The strong-line sources are estimated to have star-formation rates at least a factor 3-4 higher than those in the weak-line objects. Our conclusion is consistent with earlier work, generally carried out using much smaller samples, and reinforces the general picture of high-excitation radio galaxies as being located in lower-mass, less evolved host galaxies than their low-excitation counterparts.

The stellar mass and metallicity are among the fundamental parameters of galaxies. An understanding of the interplay between those properties as well as their environmental dependence will give us a general picture of the physics and feedback processes ongoing in groups of galaxies. We study the relationships and environmental dependencies between the stellar mass, and gas metallicity for more than 1900 galaxies in groups up to redshift 0.35 using the Galaxy And Mass Assembly (GAMA) survey. Using a control sample of more than 28,000 star-forming field galaxies, we find evidence for a decrement of the gas metallicity for galaxies in groups.

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 present a pan-chromatic analysis of an unprecedented sample of 1402 250 micron-selected galaxies at z < 0.5 (mean z = 0.24) from the Herschel-ATLAS survey. We complement our Herschel 100-500 micron data with UV-K-band photometry from the Galaxy And Mass Assembly (GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic SEDs for a representative sample of 250 micron selected galaxies spanning the most recent 5 Gyr of cosmic history. We derive estimates of physical parameters, including star formation rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at z < 0.5 has a far-infrared luminosity in the range 10^10 - 10^12 Lsolar (SFR: 1-50 Msolar/yr) thus is broadly representative of normal star forming galaxies over this redshift range. We show that 250 micron-selected galaxies contain a larger mass of dust at a given infra-red luminosity or star formation rate than previous samples selected at 60 micron from IRAS. We derive typical SEDs for H-ATLAS galaxies, and show that the emergent SED shape is most sensitive to specific star formation rate. The optical-UV SEDs also become more reddened due to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red templates. They may therefore be most appropriate for inferring total IR luminosities from moderate redshift submillimetre selected samples and for inclusion in models of the lower redshift submillimetre galaxy populations.

The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work we apply twenty published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the inter-halo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe super-halo scales accurately identify high-density regions corresponding to high mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed.

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 identify near-infrared Ks band counterparts to Herschel-ATLAS sub-mm sources, using a preliminary object catalogue from the VISTA VIKING survey. The sub-mm sources are selected from the H-ATLAS Phase 1 catalogue of the GAMA 9h field, which includes all objects detected at 250, 350 or 500 um with the SPIRE instrument. We apply and discuss a likelihood ratio (LR) method for VIKING candidates within a search radius of 10" of the 22,000 SPIRE sources with a 5 sigma detection at 250 um. We find that 11,294(51%) of the SPIRE sources have a best VIKING counterpart with a reliability $R\ge 0.8$, and the false identification rate of these is estimated to be 4.2%. We expect to miss ~5% of true VIKING counterparts. There is evidence from Z-J and J-Ks colours that the reliable counterparts to SPIRE galaxies are marginally redder than the field population. We obtain photometric redshifts for ~68% of all (non-stellar) VIKING candidates with a median redshift of 0.405. Comparing to the results of the optical identifications supplied with the Phase I catalogue, we find that the use of medium-deep near-infrared data improves the identification rate of reliable counterparts from 36% to 51%.

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 describe the infrared properties of sources detected over ~36 deg^2 of sky in the GAMA 15-hr equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H-ATLAS) and Wide-field Infrared Survey (WISE). With 5-sigma point-source depths of 34 and 0.048 mJy at 250 micron and 3.4 micron, respectively, we are able to identify 50.6% of the H-ATLAS sources in the WISE survey, corresponding to a surface density of ~630 deg^{-2}. Approximately two-thirds of these sources have measured spectroscopic or optical/near-IR photometric redshifts of z<1. For sources with spectroscopic redshifts at z<0.3, we find a linear correlation between the infrared luminosity at 3.4 micron and that at 250 micron, with +-50% scatter over ~1.5 orders of magnitude in luminosity, ~10^9 - 10^{10.5} L_sun. By contrast, the matched sources without previously measured redshifts (r>~20.5) have 250-350 micron flux density ratios that suggest either high-redshift galaxies (z>~1.5) or optically faint low-redshift galaxies with unusually low temperatures (T<~20). Their small 3.4-250 micron flux ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large AGN fraction (~30%) in a 12 micron flux-limited subsample of H-ATLAS sources, also consistent with there being a significant population of high-redshift sources in the no-redshift sample.

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 near-IR surface photometry (2D-profiling) for a sample of 29 nearby galaxies for which super-massive black hole (SMBH) masses are constrained. The data is derived from the UKIDSS-LASS survey representing a significant improvement in image quality and depth over previous studies based on 2MASS data. We derive the spheroid luminosity and spheroid S\'ersic index for each galaxy with GALFIT3 and use these data to construct SMBH mass -bulge luminosity ($M_{\rm bh}$--$L$) and SMBH - S\'ersic index ($M_{\rm bh}$--$n$) relations. The best fit K-band relation for elliptical and disk galaxies is $\log(M_{\rm bh}/M_{\odot})= -0.36(\pm 0.03) (M_{\rm K} + 18) + 6.17(\pm 0.16)$ with an intrinsic scatter of 0.4$^{+0.09}_{-0.06}$dex whilst for elliptical galaxies we find $\log(M_{\rm bh}/M_{\odot})= -0.42(\pm 0.06) (M_{\rm K} + 22) + 7.5(\pm 0.15)$ with an intrinsic scatter of 0.31$^{+0.087}_{-0.047}$dex. Our revised $M_{\rm bh}$--$L$ relation agrees closely with the previous near-IR constraint by \citet{tex:G07}. The lack of improvement in the intrinsic scatter in moving to higher quality near-IR data suggests that the SMBH relations are not currently limited by the quality of the imaging data but is either intrinsic or a result of uncertainty in the precise number of required components required in the profiling process. Contrary to expectation (see \citealt{tex:GD07a}) a relation between SMBH mass and the S\'ersic index was not found at near-IR wavelengths. This latter outcome is believed to be explained by the generic inconsistencies between 1D and 2D galaxy profiling which are currently under further investigation.

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 the first direct and unbiased measurement of the evolution of the dust mass function of galaxies over the past 5 billion years of cosmic history using data from the Science Demonstration Phase of the Herschel-ATLAS. The sample consists of galaxies selected at 250{\mu}m which have reliable counterparts from SDSS at z < 0.5, and contains 1867 sources. Dust masses are calculated using both a single temperature grey-body model for the spectral energy distribution and also using a model with multiple temperature components. The dust temperature for either model shows no trend with redshift. Splitting the sample into bins of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At z = 0.4 - 0.5, massive galaxies had dust masses about five times larger than in the local Universe. At the same time, the dust-to-stellar mass ratio was about 3-4 times larger, and the optical depth derived from fitting the UV-sub-mm data with an energy balance model was also higher. This increase in the dust content of massive galaxies at high redshift is difficult to explain using standard dust evolution models and requires a rapid gas consumption timescale together with either a more top-heavy IMF, efficient mantle growth, less dust destruction or combinations of all three. This evolution in dust mass is likely to be associated with a change in overall ISM mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.

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 describe the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) Early Release Science (ERS) observations in the Great Observatories Origins Deep Survey (GOODS) South field. The new WFC3 ERS data provide calibrated, drizzled mosaics in the UV filters F225W, F275W, and F336W, as well as in the near-IR filters F098M (Ys), F125W (J), and F160W (H) with 1-2 HST orbits per filter. Together with the existing HST Advanced Camera for Surveys (ACS) GOODS-South mosaics in the BViz filters, these panchromatic 10-band ERS data cover 40-50 square arcmin at 0.2-1.7 {\mu}m in wavelength at 0.07-0.15" FWHM resolution and 0.090" Multidrizzled pixels to depths of AB\simeq 26.0-27.0 mag (5-{\sigma}) for point sources, and AB\simeq 25.5-26.5 mag for compact galaxies. In this paper, we describe: a) the scientific rationale, and the data taking plus reduction procedures of the panchromatic 10-band ERS mosaics; b) the procedure of generating object catalogs across the 10 different ERS filters, and the specific star-galaxy separation techniques used; and c) the reliability and completeness of the object catalogs from the WFC3 ERS mosaics. The excellent 0.07-0.15" FWHM resolution of HST/WFC3 and ACS makes star- galaxy separation straightforward over a factor of 10 in wavelength to AB\simeq 25-26 mag from the UV to the near-IR, respectively.

We calculate the local UV flux density in the GALEX MIS FUV and NUV bands using redshifts provided by SDSS DR7. Luminosity functions are calculated for the overlapping MIS and SDSS sample, allowing flux densities to be measured and the local star formation rate (SFR) to be calculated using volumes much larger than previous FUV based estimates. We calculate flux densities for a dust corrected low redshift (0.013 < z < 0.1) sample of f[FUV] = 22.24 \pm 3.13 \times 10^25 h ergs s^-1 Hz^-1 Mpc^-3, f[NUV] = 38.54 \pm 5.30 \times 10^25 h ergs s^-1 Hz^-1 Mpc^-3. The star formation rate density found is 0.0312\pm0.0045 h M\odot yr^-1 Mpc^-3. This is larger than published rates recently found using the UV implied SFR, though the major discrepancy is the correction made for dust attenuation and once this is dealt with consistently the results agree well. These values are also consistent with recent Halpha derived SFRs. Once cosmic variance is taken into account most of the recent SFRs at low redshift (z < 0.3) found in the literature can be brought into agreement, however the lowest redshift values (z < 0.045) do appear to be significantly lower.

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.

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

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.

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 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/

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 determine an expression for the cosmic variance of any "normal" galaxy survey based on examination of M* +/- 1 mag galaxies in the SDSS DR7 data cube. We find that cosmic variance will depend on a number of factors principally: total survey volume, survey aspect ratio, and whether the area surveyed is contiguous or comprised of independent sight-lines. As a rule of thumb cosmic variance falls below 10% once a volume of 10^7h_0.7^-3Mpc^3 is surveyed for a single contiguous region with a 1:1 aspect ratio. Cosmic variance will be lower for higher aspect ratios and/or non-contiguous surveys. Extrapolating outside our test region we infer that cosmic variance in the entire SDSS DR7 main survey region is ~7% to z < 0.1. The equation obtained from the SDSS DR7 region can be generalised to estimate the cosmic variance for any density measurement determined from normal galaxies (e.g., luminosity densities, stellar mass densities and cosmic star-formation rates) within the volume range 10^3 to 10^7 h^-3_0.7Mpc^3. We apply our equation to show that 2 sightlines are required to ensure cosmic variance is <10% in any ASKAP galaxy survey (divided into dz ~0.1 intervals, i.e., ~1 Gyr intervals for z <0.5). Likewise 10 MeerKAT sightlines will be required to meet the same conditions. GAMA, VVDS, and zCOSMOS all suffer less than 10% cosmic variance (~3%-8%) in dz intervals of 0.1, 0.25, and 0.5 respectively. Finally we show that cosmic variance is potentially at the 50-70% level, or greater, in the HST Ultra Deep Field depending on assumptions as to the evolution of clustering. 100 or 10 independent sightlines will be required to reduce cosmic variance to a manageable level (<10%) for HST ACS or HST WFC3 surveys respectively (in dz ~ 1 intervals). Cosmic variance is therefore a significant factor in the z>6 HST studies currently underway.

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 explore the shape of the galaxy luminosity function (LF) in groups of different mass by creating composite LFs over large numbers of groups. Following previous work using total group luminosity as the mass indicator, here we split our groups by multiplicity and by estimated virial (group halo) mass, and consider red (passive) and blue (star forming) galaxies separately. In addition we utilise two different group catalogues (2PIGG and Yang et al.) in order to ascertain the impact of the specific grouping algorithm and further investigate the environmental effects via variations in the LF with position in groups. Our main results are that LFs show a steepening faint end for early type galaxies as a function of group mass/ multiplicity, with a much suppressed trend (evident only in high mass groups) for late type galaxies. Variations between LFs as a function of group mass are robust irrespective of which grouping catalogue is used, and broadly speaking what method for determining group `mass' is used. We find in particular that there is a significant deficit of low-mass passive galaxies in low multiplicity groups, as seen in high redshift clusters. Further to this, the variation in the LF appears to only occur in the central regions of systems, and in fact seems to be most strongly dependent on the position in the group relative to the virial radius. Finally, distance-rank magnitude relations were considered. Only the Yang groups demonstrated any evidence of a correlation between a galaxy's position relative to the brightest group member and its luminosity. 2PIGG possessed no such gradient, the conclusion being the FOF algorithm suppresses the signal for weak luminosity--position trends and the Yang grouping algorithm naturally enhances it.

We combine data from the MGC, SDSS and UKIDSS LAS surveys to produce ugrizYJHK luminosity functions and densities from within a common, low redshift volume (z<0.1, ~71,000 h_1^-3 Mpc^3 for L* systems) with 100 per cent spectroscopic completeness. In the optical the fitted Schechter functions are comparable in shape to those previously reported values but with higher normalisations (typically 0, 30, 20, 15, 5 per cent higher phi*-values in u, g, r, i, z respectively over those reported by the SDSS team). We attribute these to differences in the redshift ranges probed, incompleteness, and adopted normalisation methods. In the NIR we find significantly different Schechter function parameters (mainly in the M* values) to those previously reported and attribute this to the improvement in the quality of the imaging data over previous studies. This is the first homogeneous measurement of the extragalactic luminosity density which fully samples both the optical and near-IR regimes. Unlike previous compilations that have noted a discontinuity between the optical and near-IR regimes our homogeneous dataset shows a smooth cosmic spectral energy distribution (CSED). After correcting for dust attenuation we compare our CSED to the expected values based on recent constraints on the cosmic star-formation history and the initial mass function.

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

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.

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.

We derive probability density functions for the projected axial ratios of the real and mock 2PIGG galaxy groups, and use this data to investigate the intrinsic three dimensional shape of the dark matter ellipsoids that they trace. As well as analysing the raw data for groups of varying multiplicities, a convolution corrected form of the data is also considered which weights the probability density function according to the results of multiple Monte-Carlo realizations of discrete samples from the input spatial distributions. The important effect observed is that the best fit distribution for all the raw data is a prolate ellipsoid with a Gaussian distribution of axial ratios with $\bar{\beta}=0.36$ and $\sigma=0.14$, whilst for the convolved data the best fit solution is that of an oblate ellipsoid $\bar{\beta}=0.22$ and $\sigma=0.1$. Previously only prolate distributions were thought compatible with the data, this being interprated as evidence of filamentary collapse at nodes. We also find that even after allowing for the sampling effects, the corrected data is better fit using separate multiplicity bins, which display a trend towards more spherical halos in higher multiplicity groups. Finally, we find that all results in the real data are in good agreement with the mock data from $\Lambda$CDM simulations, KS tests showing that all comparative data have been drawn from the same distributions within the $1\sigma$ confidence limits.

We derive composite luminosity functions (LF) for galaxies in groups and examine the behaviour of the LF as a function of group luminosity (used as an indicator of group or halo mass). We consider both the entire galaxy population and split galaxies into red and blue (quiescent and star forming) samples, in order to examine possible mechanisms behind observed variations of galaxy properties with environment. We find evidence that $M^*$ brightens and $\alpha$ steepens with group luminosity, until a threshold value where the LF parameters stabilize at those found in rich clusters. The effect is seen in the total LF and for the blue and red galaxies separately. The behaviour of the quiescent and star-forming samples is qualitatively consistent with variations resulting from interactions and mergers, where mergers build the bright end of the luminosity function at the same time as dwarf irregulars have their star formation quenched and evolve into dwarf ellipticals. These processes appear to take place preferentially in low luminosity groups and to be complete at a group luminosity of -22.5 in $B$, corresponding to a halo mass of order $10^{13.5} M_{\odot}$.