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

We show that correlations between the phases of the galaxy density field in redshift space provide additional information about the growth rate of large-scale structure that is complementary to the power spectrum multipoles. In particular, we consider the multipoles of the line correlation function (LCF), which correlates phases between three collinear points, and use the Fisher forecasting method to show that the LCF multipoles can break the degeneracy between the measurement of the growth rate of structure $f$ and the amplitude of perturbations $\sigma_8$ that is present in the power spectrum multipoles at large scales. This leads to an improvement in the measurement of $f$ and $\sigma_8$ by up to 220 per cent for $k_{\rm max} = 0.15 \, h\mathrm{Mpc}^{-1}$ and up to 50 per cent for $k_{\rm max} = 0.30 \, h\mathrm{Mpc}^{-1}$ at redshift $z=0.25$, with respect to power spectrum measurements alone for the upcoming generation of galaxy surveys like DESI and Euclid. The average improvements in the constraints on $f$ and $\sigma_8$ for $k_{\rm max} = 0.15 \, h\mathrm{Mpc}^{-1}$ are $\sim 90$ per cent for the DESI BGS sample with mean redshift $\overline{z}=0.25$, $\sim 40$ per cent for the DESI ELG sample with $\overline{z}=1.25$, and $\sim 40$ per cent for the Euclid H$\alpha$ galaxies with $\overline{z}=1.3$. For $k_{\rm max} = 0.30 \, h\mathrm{Mpc}^{-1}$, the average improvements are $\sim 40$ per cent for the DESI BGS sample and $\sim 20$ per cent for both the DESI ELG and Euclid H$\alpha$ galaxies.

We compare the performance of four state-of-the-art models for extracting isotropic measurements of the Baryon Acoustic Oscillation (BAO) scale. To do this, we created a new, public, modular code Barry, which contains datasets, model fitting tools, and model implementations incorporating different descriptions of non-linear physics and algorithms for isolating the BAO feature. These are then evaluated for bias, correlation, and fitting strength using mock power spectra and correlation functions developed for the Sloan Digital Sky Survey Data Release 12. Our main findings are as follows: 1) All of the models can recover unbiased constraints when fit to the pre- and post-reconstruction simulations. 2) Models that provide physical descriptions of the damping of the BAO feature (using e.g., standard perturbation or effective-field theory arguments) report smaller errors on average, although the distribution of mock $\chi^{2}$ values indicates these are underestimated. 3) Allowing the BAO damping scale to vary can provide tighter constraints for some mocks, but is an artificial improvement that only arises when noise randomly sharpens the BAO peak. 4) Unlike recent claims in the literature when utilising a BAO Extractor technique, we find no improvement in the accuracy of the recovered BAO scale. 5) We implement a procedure for combining all models into a single consensus result that improves over the standard method without obviously underestimating the uncertainties. Overall, Barry provides an excellent framework for performing the cosmological analyses for upcoming surveys, and for rapidly testing and validating new models.

We model the large-scale linear galaxy bias $b_g(x,z)$ as a function of redshift $z$ and observed absolute magnitude threshold $x$ for broadband continuum emission from the far infrared to ultra-violet, as well as for prominent emission lines, such as the H$\alpha$, H$\beta$, Lya and [OII] lines. The modelling relies on the semi-analytic galaxy formation model GALFORM, run on the state-of-the-art $N$-body simulation SURFS with the Planck 2015 cosmology. We find that both the differential bias at observed absolute magnitude $x$ and the cumulative bias for magnitudes brighter than $x$ can be fitted with a five-parameter model: $b_g(x,z)=a + b(1+z)^e(1 + \exp{[(x-c)d]})$. We also find that the bias for the continuum bands follows a very similar form regardless of wavelength due to the mixing of star-forming and quiescent galaxies in a magnitude limited survey. Differences in bias only become apparent when an additional colour separation is included, which suggest extensions to this work could look at different colours at fixed magnitude limits. We test our fitting formula against observations, finding reasonable agreement with some measurements within $1\sigma$ statistical uncertainties, and highlighting areas of improvement. We provide the fitting parameters for various continuum bands, emission lines and intrinsic galaxy properties, enabling a quick estimation of the linear bias in any typical survey of large-scale structure.

We re-analyze data from the gravitational-wave event GW170817 and its host galaxy NGC4993 to demonstrate the importance of accurate total and peculiar velocities when measuring the Hubble constant using this nearby Standard Siren. We show that a number of reasonable choices can be made to estimate the velocities for this event, but that systematic differences remain between these measurements depending on the data used. This leads to significant changes in the Hubble constant inferred from GW170817. We present Bayesian Model Averaging as one way to account for these differences, and obtain $H_{0}=66.8^{+13.4}_{-9.2}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$. Adding additional information on the viewing angle from high resolution imaging of the radio counterpart refines this to $H_{0}=64.8^{+7.3}_{-7.2}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$. During this analysis we also present an alternative Bayesian model for the posterior on $H_{0}$ from Standard Sirens that works more closely with observed quantities from redshift and peculiar velocity surveys. Our results more accurately capture the true uncertainty on the total and peculiar velocities of NGC4993 and show that exploring how well different datasets characterize galaxy groups and the velocity field in the local Universe could improve this measurement further. These considerations impact any low-redshift distance measurement, and the improvements we suggest here can also be applied to standard candles like type Ia supernovae. GW170817 is particularly sensitive to peculiar velocity uncertainties because it is so close. For future standard siren measurements the importance of this error will decrease as (a) we will measure more distant standard sirens, and (b) the random direction of peculiar velocities will average out with more detections.

Redshifts have been so easy to measure for so long that we tend to neglect the fact that they too have uncertainties and are susceptible to systematic error. As we strive to measure cosmological parameters to better than 1% it is worth reviewing the accuracy of our redshift measurements. Surprisingly small systematic redshift errors, as low as 0.0001, can have a significant impact on the cosmological parameters we infer, such as $H_0$. Here we investigate an extensive (but not exhaustive) list of ways in which redshift estimation can go systematically astray. We review common theoretical errors, such as adding redshifts instead of multiplying by (1+z); using v=cz; and using only cosmological redshift in the estimates of luminosity and angular-diameter distances. We consider potential observational errors, such as rest wavelength precision, air to vacuum conversion at altitude, and spectrograph wavelength calibration. Finally, we explore physical effects, such as peculiar velocity corrections, galaxy internal velocities, gravitational redshifts, and overcorrecting within a bulk flow. We conclude that it would be quite easy for small systematic redshift errors to have infiltrated our data and be impacting our cosmological results. While it is unlikely that these errors are large enough to resolve the current $H_0$ tension, it remains possible, and redshift accuracy may become a limiting factor in near future experiments. With the enormous efforts going into calibrating the vertical axis of our plots (standard candles, rulers, clocks, and sirens) we argue that it is now worth paying a little more attention to the horizontal axis (redshifts).

Measurements of the growth rate of structure, $f\sigma_8$ in the low-redshift Universe allow stringent tests of the cosmological model. In this work, we provide new constraints on $f\sigma_{8}$ at an effective redshift of $z=0.03$ using the combined density and velocity fields measured by the 2MTF and 6dFGSv surveys. We do this by applying a new estimator of the redshift-space density and momentum (density-weighted velocity) power spectra, developed in the first paper of this series, to measured redshifts and peculiar velocities from these datasets. We combine this with models of the density and momentum power spectra in the presence of complex survey geometries and with an ensemble of simulated galaxy catalogues that match the survey selection functions and galaxy bias. We use these simulations to estimate the errors on our measurements and identify possible systematics. In particular, we are able to identify and remove biases caused by the non-Gaussianity of the power spectra by applying the Box-Cox transformation to the power spectra prior to fitting. After thorough validation of our methods we recover a constraint of $f\sigma_8(z_{\mathrm{eff}}=0.03)=0.404^{+0.082}_{-0.081}$ from the combined 2MTF and 6dFGSv data. This measurement is fully consistent with the expectations of General Relativity and the $\Lambda$ Cold Dark Matter cosmological model. It is also comparable and complementary to constraints using different techniques on similar data, affirming the usefulness of our method for extracting cosmology from velocity fields.

Low redshift surveys of galaxy peculiar velocities provide a wealth of cosmological information. We revisit the idea of extracting this information by directly measuring the redshift-space momentum power spectrum from such surveys. We provide a comprehensive theoretical and practical framework for estimating and fitting this from data, analogous to well understood techniques used to measure the galaxy density power spectrum from redshift surveys. We formally derive a new estimator, which includes the effects of shot noise and survey geometry; we evaluate the variance of the estimator in the Gaussian regime; we compute the optimal weights for the estimator; we demonstrate that the measurements are Gaussian distributed, allowing for easy extraction of cosmological parameters; and we explore the effects of peculiar velocity measurement errors. We finish with a proof-of-concept using realistic mock galaxy catalogues, which demonstrates that we can measure and fit both the redshift-space galaxy density and momentum power spectra from peculiar velocity surveys and that including the latter substantially improves our constraints on the growth rate of structure. We also provide theoretical descriptions for modelling the non-linear redshift-space density and momentum power spectrum multipoles, and forecasting the constraints on cosmological parameters using the Fisher information contained in these measurements for arbitrary weights. These may be useful for measurements of the galaxy density power spectrum even in the absence of peculiar velocities.

We present the final distance measurements for the 2MASS Tully-Fisher (2MTF) survey. The final 2MTF catalogue contains 2,062 nearby spiral galaxies in the CMB frame velocity range of 600 km s$^{-1}$ $< cz < 10,000$ km s$^{-1}$ with a mean velocity of 4,805 km s$^{-1}$. The main update in this release is the replacement of some archival HI data with newer ALFALFA data. Using the 2MTF template relation, we calculate the distances and peculiar velocities of all 2MTF galaxies. The mean uncertainties of the linear distance measurements are around 22\% in all three infrared bands. 2MTF measurements agree well with the distances from the Cosmicflows-3 compilation, which contains 1,117 common galaxies, including 28 with SNIa distance measurements. Using distances estimated from the `3-bands combined' 2MTF sample and a $\chi^2$ minimization method, we find best-fit bulk flow amplitudes of $308 \pm 26$ km s$^{-1}$, $318 \pm 29$ km s$^{-1}$, and $286 \pm 25$ km s$^{-1}$ at depths of $R_I = $ 20, 30 and 40 $h^{-1}{\rm Mpc}$, respectively, which is consistent with the $\Lambda$CDM model and with previous 2MTF results with different estimation techniques and a preliminary catalogue.

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

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

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

This document summarizes the science cases related to cosmology studies with the MaunaKea Spectroscopic Explorer (MSE), a highly-multiplexed (4332 fibers), wide FOV (1.5 sq deg), large aperture (11.25 m in diameter), optical/NIR (360nm to 1300nm) facility. The MSE High-z Cosmology Survey is designed to probe a large volume of the Universe with a galaxy density sufficient to measure the extremely-large-scale density fluctuations required to explore primordial non-Gaussianity and therefore inflation. We expect a measurement of the local parameter $f_{NL}$ to a precision $\sigma(f_{NL}) = 1.8$. Combining the MSE High-z Cosmology Survey data with data from a next generation CMB stage 4 experiment and existing DESI data will provide the first $5\sigma$ confirmation of the neutrino mass hierarchy from astronomical observations. In addition, the Baryonic Acoustic Oscillations (BAO) observed within the sample will provide measurements of the distance-redshift relationship in six different redshift bins between $z=1.6$ and 4.0, each with an accuracy of $\sim0.6\%$. The simultaneous measurements of Redshift Space Distortions (RSD) constrain the amplitude of the fluctuations, at a level ranging from $1.9\%$ to $3.6\%$. The proposed survey covers 10,000 ${\rm deg}^2$, measuring redshifts for three classes of target objects: Emission Line Galaxies (ELGs) with $1.6<z<2.4$, Lyman Break Galaxies (LBGs) with $2.4<z<4.0$, and quasars $2.1<z<3.5$. The ELGs and LBGs will be used as direct tracers of the underlying density field, while the Lyman-$\alpha$ forests in the quasar spectra will be utilized to probe structure. Exposures of duration 1,800sec will guarantee a redshift determination efficiency of $90\%$ for ELGS and at least $50\%$ for LBGs. The survey will represent 100 nights per year for a 5-year MSE program. Finally, three ideas for additional projects of cosmological interest are proposed.

In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity.

The low-order kinematic moments of galaxies, namely bulk flow and shear, enables us to test whether theoretical models can accurately describe the evolution of the mass density field in the nearby Universe. We use the so-called etaMLE maximum likelihood estimator in logdistance space to measure thesemoments from a combined sample of the 2MASS Tully-Fisher (2MTF) survey and the COSMICFLOWS-3 (CF3) compilation. Galaxies common between 2MTF and CF3 demonstrate a small zero-point difference of -0.016+-0.002 dex.We test the etaMLE on 16 mock 2MTF survey catalogues in order to explore how well the etaMLE recovers the true moments, and the effect of sample anisotropy. On the scale size of 37 Mpc/h, we find that the bulk flow of the local Universe is 259 +- 15 km/h in the direction is (l,b)=(300+-4, 23+-3) (Galactic coordinates). The average shear amplitude is 1.7+-0.4 h km/s/Mpc. We use a variable window function to explore the bulk and shear moments as a function of depth. In all cases, the measurements are consistent with the predictions of the L cold dark matter (LCDM) model.

Most statistical inference from cosmic large-scale structure relies on two-point statistics, i.e.\ on the galaxy-galaxy correlation function (2PCF) or the power spectrum. These statistics capture the full information encoded in the Fourier amplitudes of the galaxy density field but do not describe the Fourier phases of the field. Here, we quantify the information contained in the line correlation function (LCF), a three-point Fourier phase correlation function. Using cosmological simulations, we estimate the Fisher information (at redshift $z=0$) of the 2PCF, LCF and their combination, regarding the cosmological parameters of the standard $\Lambda$CDM model, as well as a Warm Dark Matter (WDM) model and the $f(R)$ and Symmetron modified gravity models. The galaxy bias is accounted for at the level of a linear bias. The relative information of the 2PCF and the LCF depends on the survey volume, sampling density (shot noise) and the bias uncertainty. For a volume of $1h^{-3}\rm Gpc^3$, sampled with points of mean density $\bar{n} = 2\times10^{-3} h^{3}\ \rm Mpc^{-3}$ and a bias uncertainty of 13\%, the LCF improves the parameter constraints by about 20\% in the $\Lambda$CDM cosmology and potentially even more in alternative models. Finally, since a linear bias only affects the Fourier amplitudes (2PCF), but not the phases (LCF), the combination of the 2PCF and the LCF can be used to break the degeneracy between the linear bias and $\sigma_8$, present in 2-point statistics.

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.

We create a combined sample of 10,904 late and early-type galaxies from the 2MTF and 6dFGSv surveys in order to accurately measure bulk flow in the local Universe. Galaxies and groups of galaxies common between the two surveys are used to verify that the difference in zero-points is $<0.02$ dex. We introduce a maximum likelihood estimator ($\eta$MLE) for bulk flow measurements which allows for more accurate measurement in the presence non-Gaussian measurement errors. To calibrate out residual biases due to the subtle interaction of selection effects, Malmquist bias and anisotropic sky distribution, the estimator is tested on mock catalogues generated from 16 independent large-scale GiggleZ and SURFS simulations. The bulk flow of the local Universe using the combined data set, corresponding to a scale size of 40 h$^{-1}$ Mpc, is $288\pm24$ km s$^{-1}$ in the direction $(l,b)=(296\pm6^{\circ}, 21\pm5^{\circ})$. This is the most accurate bulk flow measurement to date, and the amplitude of the flow is consistent with the $\Lambda$CDM expectation for similar size scales.

Taipan is a multi-object spectroscopic galaxy survey starting in 2017 that will cover 2pi steradians over the southern sky, and obtain optical spectra for about two million galaxies out to z<0.4. Taipan will use the newly-refurbished 1.2m UK Schmidt Telescope at Siding Spring Observatory with the new TAIPAN instrument, which includes an innovative 'Starbugs' positioning system capable of rapidly and simultaneously deploying up to 150 spectroscopic fibres (and up to 300 with a proposed upgrade) over the 6-deg diameter focal plane, and a purpose-built spectrograph operating from 370 to 870nm with resolving power R>2000. The main scientific goals of Taipan are: (i) to measure the distance scale of the Universe (primarily governed by the local expansion rate, H_0) to 1% precision, and the structure growth rate of structure to 5%; (ii) to make the most extensive map yet constructed of the mass distribution and motions in the local Universe, using peculiar velocities based on improved Fundamental Plane distances, which will enable sensitive tests of gravitational physics; and (iii) to deliver a legacy sample of low-redshift galaxies as a unique laboratory for studying galaxy evolution as a function of mass and environment. The final survey, which will be completed within 5 years, will consist of a complete magnitude-limited sample (i<17) of about 1.2x10^6 galaxies, supplemented by an extension to higher redshifts and fainter magnitudes (i<18.1) of a luminous red galaxy sample of about 0.8x10^6 galaxies. Observations and data processing will be carried out remotely and in a fully-automated way, using a purpose-built automated 'virtual observer' software and an automated data reduction pipeline. The Taipan survey is deliberately designed to maximise its legacy value, by complementing and enhancing current and planned surveys of the southern sky at wavelengths from the optical to the radio.

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 perform a detailed analysis of the covariance matrix of the spherically averaged galaxy power spectrum and present a new, practical method for estimating this within an arbitrary survey without the need for running mock galaxy simulations that cover the full survey volume. The method uses theoretical arguments to modify the covariance matrix measured from a set of small-volume cubic galaxy simulations, which are computationally cheap to produce compared to larger simulations and match the measured small-scale galaxy clustering more accurately than is possible using theoretical modelling. We include prescriptions to analytically account for the window function of the survey, which convolves the measured covariance matrix in a non-trivial way. We also present a new method to include the effects of supersample covariance and modes outside the small simulation volume which requires no additional simulations and still allows us to scale the covariance matrix. As validation, we compare the covariance matrix estimated using our new method to that from a brute force calculation using 500 simulations originally created for analysis of the Sloan Digital Sky Survey Main Galaxy Sample (SDSS-MGS). We find excellent agreement on all scales of interest for large scale structure analysis, including those dominated by the effects of the survey window, and on scales where theoretical models of the clustering normally break-down, but the new method produces a covariance matrix with significantly better signal-to-noise. Although only formally correct in real-space, we also discuss how our method can be extended to incorporate the effects of Redshift Space Distortions.

We present measurements of the velocity power spectrum and constraints on the growth rate of structure $f\sigma_{8}$, at redshift zero, using the peculiar motions of 2,062 galaxies in the completed 2MASS Tully-Fisher survey (2MTF). To accomplish this we introduce a model for fitting the velocity power spectrum including the effects of non-linear Redshift Space Distortions (RSD), allowing us to recover unbiased fits down to scales $k=0.2\,h\,{\rm Mpc}^{-1}$ without the need to smooth or grid the data. Our fitting methods are validated using a set of simulated 2MTF surveys. Using these simulations we also identify that the Gaussian distributed estimator for peculiar velocities of \cite{Watkins2015} is suitable for measuring the velocity power spectrum, but sub-optimal for the 2MTF data compared to using magnitude fluctuations $\delta m$, and that, whilst our fits are robust to a change in fiducial cosmology, future peculiar velocity surveys with more constraining power may have to marginalise over this. We obtain \textit{scale-dependent} constraints on the growth rate of structure in two bins, finding $f\sigma_{8} = [0.55^{+0.16}_{-0.13},0.40^{+0.16}_{-0.17}]$ in the ranges $k = [0.007-0.055, 0.55-0.150]\,h\,{\rm Mpc}^{-1}$. We also find consistent results using four bins. Assuming scale-\textit{independence} we find a value $f\sigma_{8} = 0.51^{+0.09}_{-0.08}$, a $\sim16\%$ measurement of the growth rate. Performing a consistency check of General Relativity (GR) and combining our results with CMB data only we find $\gamma = 0.45^{+0.10}_{-0.11}$, a remarkable constraint considering the small number of galaxies. All of our results are completely independent of the effects of galaxy bias, and fully consistent with the predictions of GR (scale-independent $f\sigma_{8}$ and $\gamma\approx0.55$).

Peculiar velocity surveys present a very promising route to measuring the growth rate of large-scale structure and its scale dependence. However, individual peculiar velocity surveys suffer from large statistical errors due to the intrinsic scatter in the relations used to infer a galaxy's true distance. In this context we use a Fisher Matrix formalism to investigate the statistical benefits of combining multiple peculiar velocity surveys. We find that for all cases we consider there is a marked improvement on constraints on the linear growth rate $f\sigma_{8}$. For example, the constraining power of only a few peculiar velocity measurements is such that the addition of the 2MASS Tully-Fisher survey (containing only $\sim2,000$ galaxies) to the full redshift and peculiar velocity samples of the 6-degree Field Galaxy Survey (containing $\sim 110,000$ redshifts and $\sim 9,000$ velocities) can improve growth rate constraints by $\sim20\%$. Furthermore, the combination of the future TAIPAN and WALLABY+WNSHS surveys has the potential to reach a $\sim3\%$ error on $f\sigma_{8}$, which will place tight limits on possible extensions to General Relativity. We then turn to look at potential systematics in growth rate measurements that can arise due to incorrect calibration of the peculiar velocity zero-point and from scale-dependent spatial and velocity bias. For next generation surveys, we find that neglecting velocity bias in particular has the potential to bias constraints on the growth rate by over $5\sigma$, but that an offset in the zero-point has negligible impact on the velocity power spectrum.

In this paper we investigate how observational effects could possibly bias cosmological inferences from peculiar velocity measurements. Specifically, we look at how bulk flow measurements are compared with theoretical predictions. Usually bulk flow calculations try to approximate the flow that would occur in a sphere around the observer. Using the Horizon Run 2 simulation we show that the traditional methods for bulk flow estimation can overestimate the magnitude of the bulk flow for two reasons: when the survey geometry is not spherical (the data do not cover the whole sky), and when the observations undersample the velocity distributions. Our results may explain why several bulk flow measurements found bulk flow velocities that seem larger than those expected in standard {\Lambda}CDM cosmologies. We recommend a different approach when comparing bulk flows to cosmological models, in which the theoretical prediction for each bulk flow measurement is calculated specifically for the geometry and sampling rate of that survey. This means that bulk flow values will not be comparable between surveys, but instead they are comparable with cosmological models, which is the more important measure.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

We present a fast iterative FFT-based reconstruction algorithm that allows for non- parallel redshift-space distortions (RSD). We test our algorithm on both N-body dark matter simulations and mock distributions of galaxies designed to replicate galaxy survey conditions. We compare solenoidal and irrotational components of the redshift distortion and show that an approximation of this distortion leads to a better estimate of the real-space potential (and therefore faster convergence) than ignoring the RSD when estimating the displacement field. Our iterative reconstruction scheme converges in two iterations for the mock samples corresponding to BOSS CMASS DR11 when we start with an approximation of the RSD. The scheme takes six iterations when the initial estimate, measured from the redshift-space overdensity, has no RSD correction. Slower convergence would be expected for surveys covering a larger angle on the sky. We show that this FFT based method provides a better estimate of the real space displacement field than a configuration space method that uses finite difference routines to compute the potential for the same grid resolution. Finally we show that a lognormal transform of the overdensity, used as a proxy for the linear overdensity, is beneficial in estimating the full displacement field from a dense sample of tracers. However the lognormal transform of the overdensity does not perform well when estimating the displacements from sparser simulations with a more realistic galaxy density.

Robust measurements based on current large-scale structure surveys require precise knowledge of statistical and systematic errors. This can be obtained from large numbers of realistic mock galaxy catalogues that mimic the observed distribution of galaxies within the survey volume. To this end we present a fast, distributed-memory, planar-parallel code, L-PICOLA, which can be used to generate and evolve a set of initial conditions into a dark matter field much faster than a full non-linear N-Body simulation. Additionally, L-PICOLA has the ability to include primordial non-Gaussianity in the simulation and simulate the past lightcone at run-time, with optional replication of the simulation volume. Through comparisons to fully non-linear N-Body simulations we find that our code can reproduce the $z=0$ power spectrum and reduced bispectrum of dark matter to within 2% and 5% respectively on all scales of interest to measurements of Baryon Acoustic Oscillations and Redshift Space Distortions, but 3 orders of magnitude faster. The accuracy, speed and scalability of this code, alongside the additional features we have implemented, make it extremely useful for both current and next generation large-scale structure surveys. L-PICOLA is publicly available at https://cullanhowlett.github.io/l-picola

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

We measure Redshift-Space Distortions (RSD) in the two-point correlation function of a sample of $63,163$ spectroscopically identified galaxies with $z < 0.2$, an epoch where there are currently only limited measurements, from the Sloan Digital Sky Survey (SDSS) Data Release 7 Main Galaxy Sample. Our sample, which we denote MGS, covers 6,813 deg$^2$ with an effective redshift $z_{eff}=0.15$ and is described in our companion paper (Paper I), which concentrates on BAO measurements. In order to validate the fitting methods used in both papers, and derive errors, we create and analyse 1000 mock catalogues using a new algorithm called PICOLA to generate accurate dark matter fields. Haloes are then selected using a friends-of-friends algorithm, and populated with galaxies using a Halo-Occupation Distribution fitted to the data. Using errors derived from these mocks, we fit a model to the monopole and quadrupole moments of the MGS correlation function. If we assume no Alcock-Paczynski (AP) effect (valid at $z=0.15$ for any smooth model of the expansion history), we measure the amplitude of the velocity field, $f\sigma_{8}$, at $z=0.15$ to be $0.49_{-0.14}^{+0.15}$. We also measure $f\sigma_{8}$ including the AP effect. This latter measurement can be freely combined with recent Cosmic Microwave Background results to constrain the growth index of fluctuations, $\gamma$. Assuming a background $\Lambda$CDM cosmology and combining with current Baryon Acoustic Oscillation data we find $\gamma = 0.64 \pm 0.09$, which is consistent with the prediction of General Relativity ($\gamma \approx 0.55$), though with a slight preference for higher $\gamma$ and hence models with weaker gravitational interactions.

We create a sample of spectroscopically identified galaxies with $z < 0.2$ from the Sloan Digital Sky Survey (SDSS) Data Release 7, covering 6813 deg$^2$. Galaxies are chosen to sample the highest mass haloes, with an effective bias of 1.5, allowing us to construct 1000 mock galaxy catalogs (described in Paper II), which we use to estimate statistical errors and test our methods. We use an estimate of the gravitational potential to "reconstruct" the linear density fluctuations, enhancing the Baryon Acoustic Oscillation (BAO) signal in the measured correlation function and power spectrum. Fitting to these measurements, we determine $D_{V}(z_{\rm eff}=0.15) = (664\pm25)(r_d/r_{d,{\rm fid}})$ Mpc; this is a better than 4 per cent distance measurement. This "fills the gap" in BAO distance ladder between previously measured local and higher redshift measurements, and affords significant improvement in constraining the properties of dark energy. Combining our measurement with other BAO measurements from BOSS and 6dFGS galaxy samples provides a 15 per cent improvement in the determination of the equation of state of dark energy and the value of the Hubble parameter at $z=0$ ($H_0$). Our measurement is fully consistent with the Planck results and the $\Lambda$CDM concordance cosmology, but increases the tension between Planck$+$BAO $H_0$ determinations and direct $H_0$ measurements.

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $\Lambda$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,\sigma$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

We present one thousand mock galaxy catalogues for the analysis of the Low Redshift Sample (LOWZ, effective redshift z ~ 10.32) of the Baryon Oscillation Spectroscopic Survey Data Releases 10 and 11. These mocks have been created following the PTHalos method of Manera13 et al. (2013) revised to include new developments. The main improvement is the introduction of a redshift dependence in the Halo Occupation Distribution in order to account for the change of the galaxy number density with redshift. These mock catalogues are used in the analyses of the LOWZ galaxy clustering by the BOSS collaboration.

Cosmological parameter constraints from the CMB power spectra alone suffer several well-known degeneracies. These degeneracies can be broken by numerical artefacts and also a variety of physical effects that become quantitatively important with high-accuracy data e.g. from the Planck satellite. We study degeneracies in models with flat and non-flat spatial sections, non-trivial dark energy and massive neutrinos, and investigate the importance of various physical degeneracy-breaking effects. We test the CAMB power spectrum code for numerical accuracy, and demonstrate that the numerical calculations are accurate enough for degeneracies to be broken mainly by true physical effects (the integrated Sachs-Wolfe effect, CMB lensing and geometrical and other effects through recombination) rather than numerical artefacts. We quantify the impact of CMB lensing on the power spectra, which inevitably provides degeneracy-breaking information even without using information in the non-Gaussianity. Finally we check the numerical accuracy of sample-based parameter constraints using CAMB and CosmoMC. In an appendix we document recent changes to CAMB's numerical treatment of massive neutrino perturbations, which are tested along with other recent improvements by our degeneracy exploration results.