Search

We present the Merian Survey, an optical imaging survey optimized for studying the physical properties of bright star-forming dwarf galaxies. Merian is carried out with two medium-band filters ($N708$ and $N540$, centered at $708$ and $540$ nm), custom-built for the Dark Energy Camera (DECam) on the Blanco telescope. Merian covers $\sim 750\,\mathrm{deg}^2$ of equatorial fields, overlapping with the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) wide, deep, and ultra-deep fields. When combined with the HSC-SSP imaging data ($grizy$), the new Merian DECam medium-band imaging allows for photometric redshift measurements via the detection of H$\rm\alpha$ and [OIII] line emission flux excess in the $N708$ and $N540$ filters, respectively, at $0.06

We present the phase one report of the Bright Star Subtraction (BSS) pipeline for the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). This pipeline is designed to create an extended PSF model by utilizing observed stars, followed by subtracting this model from the bright stars present in LSST data. Running the pipeline on Hyper Suprime-Cam (HSC) data shows a correlation between the shape of the extended PSF model and the position of the detector within the camera's focal plane. Specifically, detectors positioned closer to the focal plane's edge exhibit reduced circular symmetry in the extended PSF model. To mitigate this effect, we present an algorithm that enables users to account for the location dependency of the model. Our analysis also indicates that the choice of normalization annulus is crucial for modeling the extended PSF. Smaller annuli can exclude stars due to overlap with saturated regions, while larger annuli may compromise data quality because of lower signal-to-noise ratios. This makes finding the optimal annulus size a challenging but essential task for the BSS pipeline. Applying the BSS pipeline to HSC exposures allows for the subtraction of, on average, 100 to 700 stars brighter than 12th magnitude measured in g-band across a full exposure, with a full HSC exposure comprising ~100 detectors., Comment: 11 pages, 8 figures, published in Proc. SPIE 13101, Software and Cyberinfrastructure for Astronomy VIII, 131013N (25 July 2024)

We present the Bright Star Subtraction (BSS) pipeline for the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). This pipeline generates an extended PSF model using observed stars and subtracts the model from the bright stars in LSST data. When testing the pipeline on Hyper Suprime-Cam (HSC) data, we find that the shape of the extended PSF model depends on the location of the detector on the camera's focal plane. The closer a detector is to the edge of the focal plane, the less the extended PSF model is circularly symmetric. We introduce an algorithm that allows the user to consider the location dependency of the model., Comment: 4 pages, 1 figure; Astronomical Data Analysis Software & Systems XXXIII proceeding

The Merian survey is mapping $\sim$ 850 degrees$^2$ of the Hyper Suprime-Cam Strategic Survey Program (HSC-SSP) wide layer with two medium-band filters on the 4-meter Victor M. Blanco telescope at the Cerro Tololo Inter-American Observatory, with the goal of carrying the first high signal-to-noise (S/N) measurements of weak gravitational lensing around dwarf galaxies. This paper presents the design of the Merian filter set: N708 ($\lambda_c = 7080 \unicode{x212B}$, $\Delta\lambda = 275\unicode{x212B}$) and N540 ($\lambda_c = 5400\unicode{x212B}$, $\Delta\lambda = 210\unicode{x212B}$). The central wavelengths and filter widths of N708 and N540 were designed to detect the $\rm H\alpha$ and $\rm [OIII]$ emission lines of galaxies in the mass range $8<\rm \log M_*/M_\odot<9$ by comparing Merian fluxes with HSC broad-band fluxes. Our filter design takes into account the weak lensing S/N and photometric redshift performance. Our simulations predict that Merian will yield a sample of $\sim$ 85,000 star-forming dwarf galaxies with a photometric redshift accuracy of $\sigma_{\Delta z/(1+z)}\sim 0.01$ and an outlier fraction of $\eta=2.8\%$ over the redshift range $0.058

The Vera C. Rubin Observatory Wide-Fast Deep (WFD) sky survey will reach unprecedented surface brightness depths over tens of thousands of square degrees. Surface brightness photometry has traditionally been a challenge. Current algorithms which combine object detection with sky estimation systematically over-subtract the sky, biasing surface brightness measurements at the faint end and destroying or severely compromising low surface brightness light. While it has recently been shown that properly accounting for undetected faint galaxies and the wings of brighter objects can in principle recover a more accurate sky estimate, this has not yet been demonstrated in practice. Obtaining a consistent spatially smooth underlying sky estimate is particularly challenging in the presence of representative distributions of bright and faint objects. In this paper we use simulations of crowded and uncrowded fields designed to mimic Hyper Suprime-Cam data to perform a series of tests on the accuracy of the recovered sky. Dependence on field density, galaxy type and limiting flux for detection are all considered. Several photometry packages are utilised: Source Extractor, Gnuastro, and the LSST Science Pipelines. Each is configured in various modes, and their performance at extreme low surface brightness analysed. We find that the combination of the Source Extractor software package with novel source model masking techniques consistently produce extremely faint output sky estimates, by up to an order of magnitude, as well as returning high fidelity output science catalogues., Comment: Accepted for publication in MNRAS. 41 pages: 9 figures and 1 table in the main text, 11 figures and 2 tables in the appendix

Galaxies broadly fall into two categories: star-forming (blue) galaxies and quiescent (red) galaxies. In between, one finds the less populated ``green valley". Some of these galaxies are suspected to be in the process of ceasing their star-formation through a gradual exhaustion of gas supply or already dead and are experiencing a rejuvenation of star-formation through fuel injection. We use the Galaxy And Mass Assembly database and the Galaxy Zoo citizen science morphological estimates to compare the morphology of galaxies in the green valley against those in the red sequence and blue cloud. Our goal is to examine the structural differences within galaxies that fall in the green valley, and what brings them there. Previous results found disc features such as rings and lenses are more prominently represented in the green valley population. We revisit this with a similar sized data set of galaxies with morphology labels provided by the Galaxy Zoo for the GAMA fields based on new KiDS images. Our aim is to compare qualitatively the results from expert classification to that of citizen science. We observe that ring structures are indeed found more commonly in green valley galaxies compared to their red and blue counterparts. We suggest that ring structures are a consequence of disc galaxies in the green valley actively exhibiting characteristics of fading discs and evolving disc morphology of galaxies. We note that the progression from blue to red correlates with loosening spiral arm structure., Comment: 11 pages, 21 figures, accepted to MNRAS

We derive single S\'ersic fits and bulge-disk decompositions for 13096 galaxies at redshifts z < 0.08 in the GAMA II equatorial survey regions in the Kilo-Degree Survey (KiDS) g, r and i bands. The surface brightness fitting is performed using the Bayesian two-dimensional profile fitting code ProFit. We fit three models to each galaxy in each band independently with a fully automated Markov-chain Monte Carlo analysis: a single S\'ersic model, a S\'ersic plus exponential and a point source plus exponential. After fitting the galaxies, we perform model selection and flag galaxies for which none of our models are appropriate (mainly mergers/Irregular galaxies). The fit quality is assessed by visual inspections, comparison to previous works, comparison of independent fits of galaxies in the overlap regions between KiDS tiles and bespoke simulations. The latter two are also used for a detailed investigation of systematic error sources. We find that our fit results are robust across various galaxy types and image qualities with minimal biases. Errors given by the MCMC underestimate the true errors typically by factors 2-3. Automated model selection criteria are accurate to > 90 % as calibrated by visual inspection of a subsample of galaxies. We also present g-r component colours and the corresponding colour-magnitude diagram, consistent with previous works despite our increased fit flexibility. Such reliable structural parameters for the components of a diverse sample of galaxies across multiple bands will be integral to various studies of galaxy properties and evolution. All results are integrated into the GAMA database., Comment: 36 pages, 33 figures. Accepted for publication in MNRAS

The Vera C. Rubin Observatory will advance many areas of astronomy over the next decade with its unique wide-fast-deep multi-color imaging survey, the Legacy Survey of Space and Time (LSST). The LSST will produce approximately 20TB of raw data per night, which will be automatically processed by the LSST Science Pipelines to generate science-ready data products -- processed images, catalogs and alerts. To ensure that these data products enable transformative science with LSST, stringent requirements have been placed on their quality and scientific fidelity, for example on image quality and depth, astrometric and photometric performance, and object recovery completeness. In this paper we introduce faro, a framework for automatically and efficiently computing scientific performance metrics on the LSST data products for units of data of varying granularity, ranging from single-detector to full-survey summary statistics. By measuring and monitoring metrics, we are able to evaluate trends in algorithmic performance and conduct regression testing during development, compare the performance of one algorithm against another, and verify that the LSST data products will meet performance requirements by comparing to specifications. We present initial results using faro to characterize the performance of the data products produced on simulated and precursor data sets, and discuss plans to use faro to verify the performance of the LSST commissioning data products., Comment: 12 pages, 9 figues, Proc. SPIE 12189, Software and Cyberinfrastructure for Astronomy VII, Montreal, CA, July 2022

We present the capabilities of Galapagos--2 and Galfitm in the context of fitting 2-component profiles to galaxies, on the way to providing complete multi-band, multi-component fitting of large samples of galaxies in future surveys. We release both the code and the fit results to 234,239 objects from the DR3 of the Gama survey, a sample significantly deeper than previous works. We use stringent tests on both simulated and real data, as well as comparison to public catalogues to evaluate the advantages of using multi-band over single-band data. We show that multi-band fitting using Galfitm provides significant advantages when trying to decompose galaxies into their individual constituents, as more data are being used, by effectively being able to use the colour information buried in the individual exposures to its advantage. Using simulated data, we find that multi-band fitting significantly reduces the deviations from real parameter values, allows component sizes and S\'ersic indices to be recovered more accurately, and, by design, constrains the band-to-band variations of these parameters to more physical values. On both simulated and real data, we confirm that the SEDs of the 2 main components can be recovered to fainter magnitudes compared to using single-band fitting, which tends to recover disks and bulges to - on average - have identical SEDs when the galaxies become too faint, instead of the different SEDs they truly have. By comparing our results to those provided by other fitting codes, we confirm that they agree in general, but measurement errors can be significantly reduced by using the multi-band tools developed by the MegaMorph project. We conclude that the multi-band fitting employed by Galapagos-2 and Galfitm significantly improves the accuracy of structural galaxy parameters and enables much larger samples to be be used in a scientific analysis., Comment: 50 pages, 30 figures, accepted for publication in A&A

Galaxy populations show bimodality in a variety of properties: stellar mass, colour, specific star-formation rate, size, and S\'ersic index. These parameters are our feature space. We use an existing sample of 7556 galaxies from the Galaxy and Mass Assembly (GAMA) survey, represented using five features and the K-means clustering technique, showed that the bimodalities are the manifestation of a more complex population structure, represented by between 2 and 6 clusters. Here we use Self Organizing Maps (SOM), an unsupervised learning technique which can be used to visualize similarity in a higher dimensional space using a 2D representation, to map these five-dimensional clusters in the feature space onto two-dimensional projections. To further analyze these clusters, using the SOM information, we agree with previous results that the sub-populations found in the feature space can be reasonably mapped onto three or five clusters. We explore where the "green valley" galaxies are mapped onto the SOM, indicating multiple interstitial populations within the green valley population. Finally, we use the projection of the SOM to verify whether morphological information provided by GalaxyZoo users, for example, if features are visible, can be mapped onto the SOM-generated map. Voting on whether galaxies are smooth, likely ellipticals, or "featured" can reasonably be separated but smaller morphological features (bar, spiral arms) can not. SOMs promise to be a useful tool to map and identify instructive sub-populations in multidimensional galaxy survey feature space, provided they are large enough., Comment: 14 pages, 14 figures, accepted by MNRAS

In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population., Comment: Accepted for publication in MNRAS. GAMA Data Release 4 is available at: http://www.gama-survey.org/dr4/

We present Galaxy Zoo DECaLS: detailed visual morphological classifications for Dark Energy Camera Legacy Survey images of galaxies within the SDSS DR8 footprint. Deeper DECaLS images (r=23.6 vs. r=22.2 from SDSS) reveal spiral arms, weak bars, and tidal features not previously visible in SDSS imaging. To best exploit the greater depth of DECaLS images, volunteers select from a new set of answers designed to improve our sensitivity to mergers and bars. Galaxy Zoo volunteers provide 7.5 million individual classifications over 314,000 galaxies. 140,000 galaxies receive at least 30 classifications, sufficient to accurately measure detailed morphology like bars, and the remainder receive approximately 5. All classifications are used to train an ensemble of Bayesian convolutional neural networks (a state-of-the-art deep learning method) to predict posteriors for the detailed morphology of all 314,000 galaxies. When measured against confident volunteer classifications, the networks are approximately 99% accurate on every question. Morphology is a fundamental feature of every galaxy; our human and machine classifications are an accurate and detailed resource for understanding how galaxies evolve., Comment: Accepted by MNRAS July '21. Open access DOI below. Data at https://doi.org/10.5281/zenodo.4196266. Code at https://www.github.com/mwalmsley/zoobot. Docs at https://zoobot.readthedocs.io/. Interactive viewer at https://share.streamlit.io/mwalmsley/galaxy-poster/gz_decals_mike_walmsley.py

We estimate the Intracluster Light (ICL) component within a sample of 18 clusters detected in XMM Cluster Survey (XCS) data using deep ($\sim$ 26.8 mag) Hyper Suprime Cam Subaru Strategic Program DR1 (HSC-SSP DR1) $i$-band data. We apply a rest-frame ${\mu}_{B} = 25 \ \mathrm{mag/arcsec^{2}}$ isophotal threshold to our clusters, below which we define light as the ICL within an aperture of $R_{X,500}$ (X-ray estimate of $R_{500}$) centered on the Brightest Cluster Galaxy (BCG). After applying careful masking and corrections for flux losses from background subtraction, we recover $\sim$20% of the ICL flux, approximately four times our estimate of the typical background at the same isophotal level ($\sim$ 5%). We find that the ICL makes up about $\sim$ 24% of the total cluster stellar mass on average ($\sim$ 41% including the flux contained in the BCG within 50 kpc); this value is well-matched with other observational studies and semi-analytic/numerical simulations, but is significantly smaller than results from recent hydrodynamical simulations (even when measured in an observationally consistent way). We find no evidence for any links between the amount of ICL flux with cluster mass, but find a growth rate of $2-4$ for the ICL between $0.1 < z < 0.5$. We conclude that the ICL is the dominant evolutionary component of stellar mass in clusters from $z \sim 1$. Our work highlights the need for a consistent approach when measuring ICL alongside the need for deeper imaging, in order to unambiguously measure the ICL across as broad a redshift range as possible (e.g. 10-year stacked imaging from the Vera C. Rubin Observatory)., Comment: Accepted by MNRAS (05/01/2021), 20 pages, 17 figures

Strong gravitational lenses are a rare and instructive type of astronomical object. Identification has long relied on serendipity, but different strategies -- such as mixed spectroscopy of multiple galaxies along the line of sight, machine learning algorithms, and citizen science -- have been employed to identify these objects as new imaging surveys become available. We report on the comparison between spectroscopic, machine learning, and citizen science identification of galaxy-galaxy lens candidates from independently constructed lens catalogs in the common survey area of the equatorial fields of the GAMA survey. In these, we have the opportunity to compare high-completeness spectroscopic identifications against high-fidelity imaging from the Kilo Degree Survey (KiDS) used for both machine learning and citizen science lens searches. We find that the three methods -- spectroscopy, machine learning, and citizen science -- identify 47, 47, and 13 candidates respectively in the 180 square degrees surveyed. These identifications barely overlap, with only two identified by both citizen science and machine learning. We have traced this discrepancy to inherent differences in the selection functions of each of the three methods, either within their parent samples (i.e. citizen science focuses on low-redshift) or inherent to the method (i.e. machine learning is limited by its training sample and prefers well-separated features, while spectroscopy requires sufficient flux from lensed features to lie within the fiber). These differences manifest as separate samples in estimated Einstein radius, lens stellar mass, and lens redshift. The combined sample implies a lens candidate sky-density $\sim0.59$ deg$^{-2}$ and can inform the construction of a training set spanning a wider mass-redshift space., Comment: 22 figures, 6 tables, accepted for publication in The Astronomical Journal

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 either colour or SSFR, with the data showing a slight preference for Sersic index. In other words, we find that mean halo mass is more tightly correlated with galaxy structure than either past star formation history or current star formation rate. Our 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 closely 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., Comment: 16 pages, 7 figures; submitted to MNRAS; updated to reflect accepted version

In Sedgwick et al. (2019) we introduced and utilised a method to combat surface brightness and mass biases in galaxy sample selection, using core-collapse supernovae (CCSNe) as pointers towards their host galaxies, in order to: (i) search for low-surface brightness galaxies (LSBGs); (ii) assess the contributions of galaxies at a given mass to the star-formation-rate density (SFRD); and (iii) infer from this, using estimates of specific star-formation (SF) rate, the form of the SF-galaxy stellar mass function (GSMF). A CCSN-selection of SF-galaxies allows a probe of the form of the SFRD and GSMF deep into the dwarf galaxy mass regime. In the present work, we give improved constraints on our estimates of the SFRD and star-forming GSMF, in light of improved photometric redshift estimates required for estimates of galaxy stellar mass. The results are consistent with a power-law increase to SF-galaxy number density down to our low stellar mass limit of $\sim 10^{6.2}$ M$_{\odot}$. No deviation from the high-mass version of the surface brightness - mass relation is found in the dwarf mass regime. These findings imply no truncation to galaxy formation processes at least down to $\sim 10^{6.2}$ M$_{\odot}$., Comment: 6 pages, 1 figure. Submitted for publication in the proceedings of IAU Symposium No. 355, 2019, "The Realm of the Low-Surface-Brightness Universe"

We introduce a method for producing a galaxy sample unbiased by surface brightness and stellar mass, by selecting star-forming galaxies via the positions of core-collapse supernovae (CCSNe). Whilst matching $\sim$2400 supernovae from the SDSS-II Supernova Survey to their host galaxies using IAC Stripe 82 legacy coadded imaging, we find $\sim$150 previously unidentified low surface brightness galaxies (LSBGs). Using a sub-sample of $\sim$900 CCSNe, we infer CCSN-rate and star-formation rate densities as a function of galaxy stellar mass, and the star-forming galaxy stellar mass function. Resultant star-forming galaxy number densities are found to increase following a power-law down to our low mass limit of $\sim10^{6.4}$ M$_{\odot}$ by a single Schechter function with a faint-end slope of $\alpha = -1.41$. Number densities are consistent with those found by the EAGLE simulations invoking a $\Lambda$-CDM cosmology. Overcoming surface brightness and stellar mass biases is important for assessment of the sub-structure problem. In order to estimate galaxy stellar masses, a new code for the calculation of galaxy photometric redshifts, zMedIC, is also presented, and shown to be particularly useful for small samples of galaxies., Comment: 18 pages, 14 figures, accepted for publication in MNRAS

A fundamental bimodality of galaxies in the local Universe is apparent in many of the features used to describe them. Multiple sub-populations exist within this framework, each representing galaxies following distinct evolutionary pathways. Accurately identifying and characterising these sub-populations requires that a large number of galaxy features be analysed simultaneously. Future galaxy surveys such as LSST and Euclid will yield data volumes for which traditional approaches to galaxy classification will become unfeasible. To address this, we apply a robust $k$-means unsupervised clustering method to feature data derived from a sample of 7338 local-Universe galaxies selected from the Galaxy And Mass Assembly (GAMA) survey. This allows us to partition our sample into $k$ clusters without the need for training on pre-labelled data, facilitating a full census of our high dimensionality feature space and guarding against stochastic effects. We find that the local galaxy population natively splits into $2$, $3$, $5$ and a maximum of $6$ sub-populations, with each corresponding to a distinct ongoing evolutionary mechanism. Notably, the impact of the local environment appears strongly linked with the evolution of low-mass ($M_{*} < 10^{10}$ M$_{\odot}$) galaxies, with more massive systems appearing to evolve more passively from the blue cloud onto the red sequence. With a typical run time of $\sim3$ minutes per value of $k$ for our galaxy sample, we show how $k$-means unsupervised clustering is an ideal tool for future analysis of large extragalactic datasets, being scalable, adaptable, and providing crucial insight into the fundamental properties of the local galaxy population., Comment: 27 pages, 28 figures, 8 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Using a sample of 472 local Universe (z<0.06) galaxies in the stellar mass range 10.25 < log M*/M_sun < 10.75, we explore the variation in galaxy structure as a function of morphology and galaxy colour. Our sample of galaxies is sub-divided into red, green and blue colour groups and into elliptical and non-elliptical (disk-type) morphologies. Using KiDS and VIKING derived postage stamp images, a group of eight volunteers visually classified bars, rings, morphological lenses, tidal streams, shells and signs of merger activity for all systems. We find a significant surplus of rings ($2.3\sigma$) and lenses ($2.9\sigma$) in disk-type galaxies as they transition across the green valley. Combined, this implies a joint ring/lens green valley surplus significance of $3.3\sigma$ relative to equivalent disk-types within either the blue cloud or the red sequence. We recover a bar fraction of ~44% which remains flat with colour, however, we find that the presence of a bar acts to modulate the incidence of rings and (to a lesser extent) lenses, with rings in barred disk-type galaxies more common by ~20-30 percentage points relative to their unbarred counterparts, regardless of colour. Additionally, green valley disk-type galaxies with a bar exhibit a significant $3.0\sigma$ surplus of lenses relative to their blue/red analogues. The existence of such structures rules out violent transformative events as the primary end-of-life evolutionary mechanism, with a more passive scenario the favoured candidate for the majority of galaxies rapidly transitioning across the green valley., Comment: 17 pages, 9 figures, accepted for publication in MNRAS

We present a sample of 329 low to intermediate redshift ($0.05 < z < 0.3$) brightest cluster galaxies (BCGs) in X-ray selected clusters from the SPectroscopic IDentification of eRosita Sources (SPIDERS) survey, a spectroscopic survey within Sloan Digital Sky Survey-IV (SDSS-IV). We define our BCGs by simultaneous consideration of legacy X-ray data from ROSAT, maximum likelihood outputs from an optical cluster-finder algorithm and visual inspection. Using SDSS imaging data, we fit S\'ersic profiles to our BCGs in three bands (\textit{g}, \textit{r}, \textit{i}) with \textsc{SIGMA}, a \textsc{GALFIT}-based software wrapper. We examine the reliability of our fits by running our pipeline on ${\sim}10^{4}$ psf-convolved model profiles injected into 8 random cluster fields, we then use the results of this analysis to create a robust subsample of 198 BCGs. We outline three cluster properties of interest: overall cluster X-ray luminosity ($L_{X}$), cluster richness as estimated by \textsc{redMaPPer} ($ \lambda $) and cluster halo mass ($M_{200}$), which is estimated via velocity dispersion. In general, there are significant correlations with BCG stellar mass between all three environmental properties, but no significant trends arise with either S\'ersic index or effective radius. There is no major environmental dependence on the strength of the relation between effective radius and BCG stellar mass. Stellar mass therefore arises as the most important factor governing BCG morphology. Our results indicate that our sample consists of a large number of relaxed, mature clusters containing broadly homogeneous BCGs up to $z \sim 0.3$, suggesting that there is little evidence for much ongoing structural evolution for BCGs in these systems., Comment: 25 pages, 17 figures, 11 tables - accepted in MNRAS

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

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., Comment: Accepted for publication in MNRAS

[Abridged] Typical disc galaxies forming in a LambdaCDM cosmology encounter a violent environment, where they often experience mergers with massive satellites. The fact that disc galaxies are ubiquitous in the local Universe suggests that a quiescent history is not necessary for their formation. Modern cosmological simulations can now obtain relatively realistic populations of disc galaxies, but it still remains to be clarified how discs manage to survive massive mergers. Here we use a suite of high-resolution hydrodynamical simulations set in a LambdaCDM cosmology to elucidate the fate of discs encountering massive mergers. We extract a sample of approximately 100 disc galaxies and follow the changes in their post-merger morphologies, as tracked by their disc-to-total ratios (D/T). We also examine the relations between their present-day morphology, assembly history and gas fractions. We find that approximately half of present-day disc galaxies underwent at least one merger with a satellite of total mass exceeding the host system's stellar mass, a third had mergers with satellites of mass exceeding 3 times the host's stellar mass, and approximately one-sixth had mergers with satellites of mass exceeding 10 times of the host's stellar mass. These mergers lead to a sharp, but often temporary, decrease in the D/T of the hosts, implying that discs are usually disrupted but then quickly re-grow. To do so, high cold gas fractions are required post-merger, as well as a relatively quiescent recent history (over a few Gyrs before z=0). Our results show that discs can form via diverse merger pathways and that quiescent histories are not the dominant mode of disc formation., Comment: 21 pages, PASA, in press

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., Comment: 19 pages, 14 figures, accepted by MNRAS

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., Comment: 14 pages; 10 figures; MNRAS accepted

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., Comment: 33 pages, 14 + 7 Figures, MNRAS accepted

We study the influence of environment on the structure of disc galaxies, using \texttt{IMFIT} to measure the g- and r-band structural parameters of the surface-brightness profiles for $\sim$700 low-redshift (z$<$0.063) cluster and field disc galaxies with intermediate stellar mass (0.8 $\times$ 10$^{10}$ $M_{\odot}$ $<$ $M_{\star}$ $<$ 4 $\times$ 10$^{10}$ $M_{\odot}$) from the Sloan Digital Sky Survey, DR7. Based on this measurement, we assign each galaxy to a surface-brightness profile type (Type I $\equiv$ single-exponential, Type II $\equiv$ truncated, Type III $\equiv$ anti-truncated). In addition, we measure (g-r) restframe colour for disc regions separated by the break radius. Cluster disc galaxies (at the same stellar mass) have redder (g-r) colour by $\sim$0.2 mag than field galaxies. This reddening is slightly more pronounced outside the break radius. Cluster disc galaxies also show larger global S\'ersic-indices and are more compact than field discs, both by $\sim$15\%. This change is connected to a flattening of the (outer) surface-brightness profile of Type I and - more significantly - of Type III galaxies by $\sim$8\% and $\sim$16\%, respectively, in the cluster environment compared to the field. We find fractions of Type I, II and III of (6$\pm$2)\%, (66$\pm$4)\% and (29$\pm$4)\% in the field and (15$_{-4}^{+7}$)\%, (56$\pm$7)\% and (29$\pm$7)\% in the cluster environment, respectively. We suggest that the larger abundance of Type I galaxies in clusters (matched by a corresponding decrease in the Type II fraction) could be the signature of a transition between Type II and Type I galaxies produced/enhanced by environment-driven mechanisms., Comment: accepted, MNRAS, January 19, 2017

With a large sample of bright, low-redshift galaxies with optical$-$near-IR imaging from the GAMA survey we use bulge-disc decompositions to understand the wavelength-dependent behavior of single-S\'ersic structural measurements. We denote the variation in single-S\'ersic index with wavelength as $\mathcal{N}$, likewise for effective radius we use $\mathcal{R}$. We find that most galaxies with a substantial disc, even those with no discernable bulge, display a high value of $\mathcal{N}$. The increase in S\'ersic index to longer wavelengths is therefore intrinsic to discs, apparently resulting from radial variations in stellar population and/or dust reddening. Similarly, low values of $\mathcal{R}$ ($<$ 1) are found to be ubiquitous, implying an element of universality in galaxy colour gradients. We also study how bulge and disc colour distributions vary with galaxy type. We find that, rather than all bulges being red and all discs being blue in absolute terms, both components become redder for galaxies with redder total colours. We even observe that bulges in bluer galaxies are typically bluer than discs in red galaxies, and that bulges and discs are closer in colour for fainter galaxies. Trends in total colour are therefore not solely due to the colour or flux dominance of the bulge or disc., Comment: 16 pages, 13 figures, accepted for publication in MNRAS

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., Comment: 16 pages, 15 figures (including one 3D figure that can be interacted with when reading the PDF in Adobe Reader). Accepted for publication in MNRAS on January 13th, 2016

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., Comment: 13 pages, 9 figures (high-resolution figures available in journal version), MNRAS accepted

We study how the sizes and radial profiles of galaxies vary with wavelength, by fitting Sersic functions simultaneously to imaging in nine optical and near-infrared bands. To quantify the wavelength dependence of effective radius we use the ratio, $\mathcal{R}$, of measurements in two restframe bands. The dependence of Sersic index on wavelength, $\mathcal{N}$, is computed correspondingly. Vulcani et al. (2014) have demonstrated that different galaxy populations present sharply contrasting behaviour in terms of $\mathcal{R}$ and $\mathcal{N}$. Here we study the luminosity dependence of this result. We find that at higher luminosities, early-type galaxies display a more substantial decrease in effective radius with wavelength, whereas late-types present a more pronounced increase in Sersic index. The structural contrast between types thus increases with luminosity. By considering samples at different redshifts, we demonstrate that lower data quality reduces the apparent difference between the main galaxy populations. However, our conclusions remain robust to this effect. We show that accounting for different redshift and luminosity selections partly reconciles the size variation measured by Vulcani et al. with the weaker trends found by other recent studies. Dividing galaxies by visual morphology confirms the behaviour inferred using morphological proxies, although the sample size is greatly reduced. Finally, we demonstrate that varying dust opacity and disc inclination can account for features of the joint distribution of $\mathcal{R}$ and $\mathcal{N}$ for late-type galaxies. However, dust does not appear to explain the highest values of $\mathcal{R}$ and $\mathcal{N}$. The bulge-disc nature of galaxies must also contribute to the wavelength-dependence of their structure., Comment: 14 pages, 11 figures, accepted by MNRAS; v2 corrected to match journal version

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/, Comment: 31 pages and 30 figures. Accepted in MNRAS. High-resolution copy available from our release site: http://gama-psi.icrar.org/ or directly via http://www.simondriver.org/mwavev05.pdf

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., Comment: Accepted for publication in MNRAS

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., Comment: 22 pages, 18 figures, accepted for publication in MNRAS on May 20, 2015

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.

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

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}

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., Comment: Accepted by MNRAS - 16 figures, 42 pages including 1 appendix - all results will be made available online via MNRAS, source code can be provided on request

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., Comment: 16 pages, 6 figures, accepted for publication in MNRAS (see MNRAS version for high-quality figures)

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., Comment: 24 pages, 21 figures, published in MNRAS (v2 fixes author names; v3 matches published version)

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., Comment: 6 pages, 3 figures, accepted for publication in MNRAS Letters on Tue 28 Jan 2014

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., Comment: 27 pages, 22 figures, accepted for publication in MNRAS (see MNRAS version for high-quality figures)

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., Comment: 16 pages, 15 figures, 7 tables. Accepted for publication in Publications of the Astronomical Society of Australia

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., Comment: Accepted for MNRAS; 41 pages, 32 figures, 6 tables. This is Paper III in a logical series of papers. Paper I (Bamford et al.) and Paper II (Vika et al.) are in prep and will be available soon. Several Figures have been degraded for file size reasons. Please follow http://www.nottingham.ac.uk/~ppzbh2/MM_III_Haeussler.pdf for a higher quality version of this paper

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., Comment: 16 pages, 12 figures

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., Comment: 36 pages, 25 figures, accepted for publication in MNRAS (7 December 2011). High resolution version may be found at http://www.gama-survey.org/pubs/papers/Kelvin_2011.pdf

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., Comment: Accepted for publication in MNRAS. 35 pages; 21 figures including 6 in 2 appendices -- v2. corrected three extremely embarrassing typos in coauthors' names as they appeared in the arXiv metadata; manuscript text is unchanged

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%., Comment: 38 pages, 10 Tables, 26 figures. Submitted to MNRAS (revised once). Image resolution has been lowered. For higher resolution, see http://www.gama-survey.org/