
Observations of galaxies and galaxy clusters in the local universe can
account for only $\sim\,10\%$ of the total baryon content. Cosmological
simulations predict that the `missing baryons' are spread throughout
filamentary structures in the cosmic web, forming a lowdensity gas with
temperatures of $10^510^7\,\!$K. We search for this warmhot intergalactic
medium (WHIM) by stacking the Planck Compton $y$parameter map of the thermal
SunyaevZel'dovich (tSZ) effect for 1,002,334 pairs of CMASS galaxies from the
Sloan Digital Sky Survey. We model the contribution from the galaxy halo pairs
assuming spherical symmetry, finding a residual tSZ signal at the
$2.9\mbox{$\sigma$}$ level from a stacked filament of length $10.5\,h^{1}\,\rm
Mpc$ with a Compton parameter magnitude $y=(0.6\pm0.2)\times10^{8}$. We
consider possible sources of contamination and conclude that bound gas in
haloes may contribute only up to $20\%$ of the measured filamentary signal. To
estimate the filament gas properties we measure the gravitational lensing
signal for the same sample of galaxy pairs; in combination with the tSZ signal,
this yields an inferred gas density of $\rho_{\rm b}=(5.5\pm
2.9)\times\bar{\rho_{\rm b}}$ with a temperature $T=(2.7\pm 1.7) \times
10^6\,$K. This result is consistent with the predicted WHIM properties, and
overall the filamentary gas can account for $ 11\pm 7\%$ of the total baryon
content of the Universe. We also see evidence that the gas filament extends
beyond the galaxy pair. Averaging over this longer baseline boosts the
significance of the tSZ signal and increases the associated baryon content to
$28\pm 12\%$ of the global value.

When combining data sets to perform parameter inference, the results will be
unreliable if there are unknown systematics in data or models. Here we
introduce a flexible methodology, BACCUS: BAyesian Conservative Constraints and
Unknown Systematics, which deals in a conservative way with the problem of data
combination, for any degree of tension between experiments. We introduce
hyperparameters that describe a bias in each model parameter for each class of
experiments. A conservative posterior for the model parameters is then obtained
by marginalization both over these unknown shifts and over the width of their
prior. We contrast this approach with an existing hyperparameter method in
which each individual likelihood is scaled, comparing the performance of each
approach and their combination in application to some idealized models. Using
only these rescaling hyperparameters is not a suitable approach for the current
observational situation, in which internal null tests of the errors are passed,
and yet different experiments prefer models that are in poor agreement. The
possible existence of large shift systematics cannot be constrained with a
small number of data sets, leading to extended tails on the conservative
posterior distributions. We illustrate our method with the case of the $H_0$
tension between results from the cosmic distance ladder and physical
measurements that rely on the standard cosmological model.

Galaxies exhibit different clustering and quenching properties in clusters,
filaments, and the field, but it is still uncertain whether such differences
are imprints of the tidal environment on galaxy formation, or if they reflect
the variation of the underlying halo mass function across the cosmic web. We
measure the dependence of galaxy clustering and quenching on the cosmic web in
the Sloan Digital Sky Survey, characterized by the combination of spherical
overdensity $\delta_8$ and tidal anisotropy $\alpha_5$ centred on each galaxy.
We find that galaxy clustering is a strong function of either $\delta_8$ or
$\alpha_5$, and the largescale galaxy bias shows complex and rich behaviour on
the $\delta_8$ vs. $\alpha_5$ plane. Using the mean galaxy colour as a proxy
for the average quenched level of galaxies, we find that galaxy quenching is
primarily a function of $\delta_8$, with some subtle yet nontrivial dependence
on $\alpha_5$ at fixed $\delta_8$. The quenched galaxies generally show
stronger smallscale clustering than the starforming ones at fixed $\delta_8$
or $\alpha_5$, while the characteristic scale at which the amplitude of
clustering becomes comparable for both galaxy populations varies with
$\delta_8$ and $\alpha_5$. We compare these observed cosmic web dependences of
galaxy clustering and quenching with a mock galaxy catalogue constructed from
the iHOD model, which places quenched and starforming galaxies inside dark
matter haloes based on the stellartohalo mass relation and the halo quenching
model  the $\delta_8$ and $\alpha_5$ dependences of iHOD galaxies are thus
solely derived from the cosmic web modulation of the halo mass function. The
main observed trends are accounted for extremely well by the iHOD model. Thus
any additional direct effect of the largescale~(${>}5\,h^{1}{\rm Mpc}$) tidal
field on galaxy formation must be extremely weak. [Abridged]

In theories of modified gravity with the chameleon screening mechanism, the
strength of the fifth force depends on environment. This induces an environment
dependence of structure formation, which differs from $\Lambda$CDM. We show
that these differences can be captured by the marked correlation function. With
the galaxy correlation functions and number densities calibrated to match
between $f(R)$ and $\Lambda$CDM models in simulations, we show that the marked
correlation functions from using either the local density or halo mass as the
marks encode extra information, which can be used to test these theories. We
discuss possible applications of these statistics in observations.

We develop a new method, which is based on the optimal redshift weighting
scheme, to extract the maximal tomographic information of baryonic acoustic
oscillations (BAO) and redshift space distortions (RSD) from the extended
Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q)
survey. We validate our method using the EZ mocks, and apply our pipeline to
the eBOSS DR14Q sample in the redshift range of $0.8<z<2.2$. We report a joint
measurement of $f\sigma_8$ and twodimensional BAO parameters $D_{\rm A}$ and
$H$ at four effective redshifts of $z_{\rm eff}=0.98, 1.23, 1.52$ and $1.94$,
and provide the full data covariance matrix. Using our measurement combined
with BOSS DR12, MGS and 6dFGS BAO measurements, we find that the existence of
dark energy is supported by observations at a $7.4\sigma$ significance level.
Combining our measurement with BOSS DR12 and Planck observations, we constrain
the gravitational growth index to be $\gamma=0.580\pm0.082$, which is fully
consistent with the prediction of general relativity. This paper is part of a
set that analyses the eBOSS DR14 quasar sample.

The fourth generation of the Sloan Digital Sky Survey (SDSSIV) began
observations in July 2014. It pursues three core programs: APOGEE2, MaNGA, and
eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS.
This paper describes the first data release from SDSSIV, Data Release 13
(DR13), which contains new data, reanalysis of existing data sets and, like all
SDSS data releases, is inclusive of previously released data. DR13 makes
publicly available 1390 spatially resolved integral field unit observations of
nearby galaxies from MaNGA, the first data released from this survey. It
includes new observations from eBOSS, completing SEQUELS. In addition to
targeting galaxies and quasars, SEQUELS also targeted variabilityselected
objects from TDSS and Xray selected objects from SPIDERS. DR13 includes new
reductions of the SDSSIII BOSS data, improving the spectrophotometric
calibration and redshift classification. DR13 releases new reductions of the
APOGEE1 data from SDSSIII, with abundances of elements not previously
included and improved stellar parameters for dwarf stars and cooler stars. For
the SDSS imaging data, DR13 provides new, more robust and precise photometric
calibrations. Several valueadded catalogs are being released in tandem with
DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and
an updated redclump catalog for APOGEE. This paper describes the location and
format of the data now publicly available, as well as providing references to
the important technical papers that describe the targeting, observing, and data
reduction. The SDSS website, http://www.sdss.org, provides links to the data,
tutorials and examples of data access, and extensive documentation of the
reduction and analysis procedures. DR13 is the first of a scheduled set that
will contain new data and analyses from the planned ~6year operations of
SDSSIV.

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

We have searched for the signature of cosmic voids in the CMB, in both the
Planck temperature and lensingconvergence maps; voids should give decrements
in both. We use zobov voids from the DR12 SDSS CMASS galaxy sample. We base our
analysis on Nbody simulations, to avoid a posteriori bias. For the first time,
we detect the signature of voids in CMB lensing: the significance is
$3.2\sigma$, close to $\Lambda$CDM in both amplitude and projected
densityprofile shape. A temperature dip is also seen, at modest significance
($2.3\sigma$), with amplitude about 6 times the prediction. This temperature
signal is induced mostly by voids with radius between 100 and 150 Mpc/h, while
the lensing signal is mostly contributed by smaller voids  as expected;
lensing relates directly to density, while ISW depends on gravitational
potential. The void abundance in observations and simulations agree, as well.
We also repeated the analysis excluding lowersignificance voids: no lensing
signal is detected, with an upper limit of about twice the $\Lambda$CDM
prediction. But the mean temperature decrement now becomes nonzero at the
$3.7\sigma$ level (similar to that found by Granett et al.), with amplitude
about 20 times the prediction. However, the observed dependence of temperature
on void size is in poor agreement with simulations, whereas the lensing results
are consistent with $\Lambda$CDM theory. Thus, the overall tension between
theory and observations does not favour nonstandard theories of gravity,
despite the hints of an enhanced amplitude for the ISW effect from voids.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea galaxy
and quasar redshift survey. To trace the underlying dark matter distribution,
spectroscopic targets will be selected in four classes from imaging data. We
will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to
even higher redshift, DESI will target bright [O II] emission line galaxies up
to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying
dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the
Ly$\alpha$ forest absorption features in their spectra, which will be used to
trace the distribution of neutral hydrogen. When moonlight prevents efficient
observations of the faint targets of the baseline survey, DESI will conduct a
magnitudelimited Bright Galaxy Survey comprising approximately 10 million
galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy
and quasar redshifts will be obtained to measure the BAO feature and determine
the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
diameter field of view. The DESI collaboration will also deliver a
spectroscopic pipeline and data management system to reduce and archive all
data for eventual public use.

In this work we reformulate the forward modelling of the redshiftspace power
spectrum multipole moments for a masked density field, as encountered in galaxy
redshift surveys. Exploiting the symmetries of the redshiftspace correlation
function, we provide a maskedfield generalisation of the Hankel transform
relation between the multipole moments in real and Fourier space. Using this
result, we detail how a likelihood analysis requiring computation for a broad
range of desired $P(k)$ models may be executed $10^310^4$ times faster than
with other common approaches, together with significant gains in spectral
resolution. We present a concrete application to the complex angular geometry
of the VIPERS PDR1 release and discuss the validity of this technique for
finiteangle surveys.

We have derived estimators for the linear growth rate of density fluctuations
using the crosscorrelation function of voids and haloes in redshift space,
both directly and in Fourier form. In linear theory, this crosscorrelation
contains only monopole and quadrupole terms. At scales greater than the void
radius, linear theory is a good match to voids traced out by haloes in Nbody
simulations; smallscale random velocities are unimportant at these radii, only
tending to cause small and often negligible elongation of the redshiftspace
crosscorrelation function near its origin. By extracting the monopole and
quadrupole from the crosscorrelation function, we measure the linear growth
rate without prior knowledge of the void profile or velocity dispersion. We
recover the linear growth parameter $\beta$ to 9% precision from an effective
volume of 3(Gpc/h)^3 using voids with radius greater than 25Mpc/h. Smaller
voids are predominantly subvoids, which may be more sensitive to the random
velocity dispersion; they introduce noise and do not help to improve the
measurement. Adding velocity dispersion as a free parameter allows us to use
information at radii as small as half of the void radius. The precision on
$\beta$ is reduced to approximately 5%. Contrary to the simple redshiftspace
distortion pattern in overdensities, voids show diverse shapes in redshift
space, and can appear either elongated or flattened along the line of sight.
This can be explained by the competing amplitudes of the local density
contrast, plus the radial velocity profile and its gradient, with the latter
two factors being determined by the cumulative density profile of voids. The
distortion pattern is therefore determined solely by the void profile and is
different for voidincloud and voidinvoid. This diversity of redshiftspace
void morphology complicates measurements of the AlcockPaczynski effect using
voids.

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

We investigate the use of the crosscorrelation between galaxies and galaxy
groups to measure redshiftspace distortions (RSD) and thus probe the growth
rate of cosmological structure. This is compared to the classical approach
based on using galaxy autocorrelation. We make use of realistic simulated
galaxy catalogues that have been constructed by populating simulated dark
matter haloes with galaxies through halo occupation prescriptions. We adapt the
classical RSD dispersion model to the case of the groupgalaxy
crosscorrelation function and estimate the RSD parameter $\beta$ by fitting
both the full anisotropic correlation function $\xi(r_p,\pi)$ and its multipole
moments. In addition, we define a modified version of the latter statistics by
truncating the multipole moments to exclude strongly nonlinear distortions at
small transverse scales. We fit these three observable quantities in our set of
simulated galaxy catalogues and estimate statistical and systematic errors on
$\beta$ for the case of galaxygalaxy, groupgroup, and groupgalaxy
correlation functions. When ignoring offdiagonal elements of the covariance
matrix in the fitting, the truncated multipole moments of the groupgalaxy
crosscorrelation function provide the most accurate estimate, with systematic
errors below 3% when fitting transverse scales larger than $10Mpc/h$. Including
the full data covariance enlarges statistical errors but keep unchanged the
level of systematic error. Although statistical errors are generally larger for
groups, the use of groupgalaxy crosscorrelation can potentially allow the
reduction of systematics while using simple linear or dispersion models.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct
novel cosmological observations using the BOSS spectrograph at Apache Point
Observatory. Observations will be simultaneous with the Time Domain
Spectroscopic Survey (TDSS) designed for variability studies and the
Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for
studies of Xray sources. eBOSS will use four different tracers to measure the
distanceredshift relation with baryon acoustic oscillations (BAO). Using more
than 250,000 new, spectroscopically confirmed luminous red galaxies at a median
redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to
an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the
z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy
redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and
H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000
spectroscopicallyconfirmed quasars will provide the first BAO distance
measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8%
and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars
and reobservation of 60,000 quasars known from BOSS, we will obtain new
Lymanalpha forest measurements at redshifts z>2.1; these new data will enhance
the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS.
Furthermore, eBOSS will provide improved tests of General Relativity on
cosmological scales through redshiftspace distortion measurements, improved
tests for nonGaussianity in the primordial density field, and new constraints
on the summed mass of all neutrino species. Here, we provide an overview of the
cosmological goals, spectroscopic target sample, demonstration of spectral
quality from early data, and projected cosmological constraints from eBOSS.

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

We present the GAMA Panchromatic Data Release (PDR) constituting over
230deg$^2$ of imaging with photometry in 21 bands extending from the farUV to
the farIR. 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 21band 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.1500$\mu$m energy output of the Universe. Exploring the
Cosmic Spectral Energy Distribution (CSED) across three timeintervals
(0.31.1Gyr, 1.11.8~Gyr and 1.82.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, subpopulations, and representative galaxy
samples at $z<0.5$. The GAMA PDR can be found at: http://gamapsi.icrar.org/

A long standing question in cosmology is whether gravitational lensing
changes the distanceredshift relation $D(z)$ or the mean flux density of
sources. Interest in this has been rekindled by recent studies in nonlinear
relativistic perturbation theory that find biases in both the area of a surface
of constant redshift and in the mean distance to this surface, with a
fractional bias in both cases on the order of the mean squared convergence
$\langle \kappa^2 \rangle$. Any such area bias could alter CMB cosmology, and
the corresponding bias in mean flux density could affect supernova cosmology.
Here we show that, in an ensemble averaged sense, the perturbation to the area
of a surface of constant redshift is in reality much smaller, being on the
order of the cumulative bending angle squared, or roughly a partinamillion
effect. This validates the arguments of Weinberg (1976) that the mean
magnification $\mu$ of sources is unity and of Kibble \& Lieu (2005) that the
mean directionaveraged inverse magnification is unity. It also validates the
conventional treatment of lensing in analysis of CMB anisotropies. But the
existence of a scatter in magnification will cause any nonlinear function of
these conserved quantities to be statistically biased. The distance $D$, for
example, is proportional to $\mu^{1/2}$ so lensing will bias $\langle
D\rangle$ even if $\langle \mu \rangle=1$. The fractional bias in such
quantities is generally of order $\langle \kappa^2 \rangle$, which is orders of
magnitude larger than the area perturbation. Claims for large bias in area or
flux density of sources appear to have resulted from misinterpretation of such
effects: they do not represent a new nonNewtonian effect, nor do they
invalidate standard cosmological analyses.

We investigate the dependence of the mass function of darkmatter haloes on
their environment within the cosmic web of largescale structure. A dependence
of the halo mass function on largescale mean density is a standard element of
cosmological theory, allowing massdependent biasing to be understood via the
peakbackground split. On the assumption of a Gaussian density field, this
analysis can be extended to ask how the mass function depends on the
geometrical environment: clusters, filaments, sheets and voids, as classified
via the tidal tensor (the Hessian matrix of the gravitational potential). In
linear theory, the problem can be solved exactly, and the result is
attractively simple: the conditional mass function has no explicit dependence
on the local tidal field, and is a function only of the local density on the
filtering scale used to define the tidal tensor. There is nevertheless a strong
implicit predicted dependence on geometrical environment, because the local
density couples statistically to the derivatives of the potential. We compute
the predictions of this model and study the limits of their validity by
comparing them to results deduced empirically from $N$body simulations. We
have verified that, to a good approximation, the abundance of haloes in
different environments depends only on their densities, and not on their tidal
structure. In this sense we find relative differences between halo abundances
in different environments with the same density which are smaller than 13%.
Furthermore, for sufficiently large filtering scales, the agreement with the
theoretical prediction is good, although there are important deviations from
the Gaussian prediction at small, nonlinear scales. We discuss how to obtain
improved predictions in this regime, using the 'effectiveuniverse' approach.

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

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

We present a method for populating dark matter simulations with haloes of
mass below the resolution limit. It is based on stochastically sampling a field
derived from the density field of the halo catalogue, using constraints from
the conditional halo mass function n(m{\delta}). We test the accuracy of the
method and show its application in the context of building mock galaxy samples.
We find that this technique allows precise reproduction of the twopoint
statistics of galaxies in mock samples constructed with this method. Our
results demonstrate that the full information content of a simulation can be
communicated efficiently using only a catalogue of the more massive haloes.

By 1917, V.M. Slipher had singlehandedly established a tendency for 'spiral
nebulae' to be redshifted (21 out of 25 cases). From a modern perspective, it
could seem surprising that the expansion of the universe was not announced at
this point. Examination of Slipher's papers shows that he reached a more subtle
conclusion: the identification of cosmological peculiar velocities, including
the bulk motion of the Milky Way, leading to a beautiful argument in favour of
nebulae as distant stellar systems. Nevertheless, Slipher's data actually
contain evidence at >8sigma for a positive mean velocity, even after
subtracting the dipole owing to the motion of the observer. In 1929, Hubble
estimated distances for a sample of no greater depth, using redshifts due
almost entirely to Slipher. Hubble's distances were flawed in two distinct
ways: in addition to an incorrect absolute calibration, the largest distances
were systematically underestimated. Nevertheless, he claimed the detection of
a linear distanceredshift relation. Statistically, the evidence for such a
correlation is less strong than the simple evidence for a positive mean
velocity in Hubble's sample. Comparison with modern data shows that a sample of
more than twice Hubble's depth would generally be required in order to reveal
clearly the global linear expansion in the face of the 'noise' from peculiar
velocities. When the theoretical context of the time is examined, the role of
the de Sitter model and its prediction of a linear distanceredshift relation
looms large. A number of searches for this relation were performed prior to
Hubble over the period 19241928, with a similar degree of success. All were
based on the velocities measured by Slipher, whose work from a Century ago
stands out both for the precision of his measurements and for the subtle
clarity of the arguments he employed to draw correct conclusions from them.

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

In order to generate credible 0.12 {\mu}m SEDs, the GAMA project requires
many Gigabytes of imaging data from a number of instruments to be reprocessed
into a standard format. In this paper we discuss the software infrastructure we
use, and create selfconsistent ugrizYJHK photometry for all sources within the
GAMA sample. Using UKIDSS and SDSS archive data, we outline the preprocessing
necessary to standardise all images to a common zeropoint, the steps taken to
correct for seeing bias across the dataset, and the creation of Gigapixelscale
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 bestfitting 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
bestfitting Schechter parameters from totalmagnitude photometric systems
(such as SDSS modelmag or S\'ersic magnitudes) have a steeper faintend slope
than photometry dependent on Kron or Petrosian magnitudes. (5) Our Universe's
total luminosity density, when calculated using Kron or Petrosian rband
photometry, is underestimated by at least 15%.