
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.

We study projected underdensities in the cosmic galaxy density field known as
'troughs', and their overdense counterparts, which we call 'ridges'. We
identify these regions using a bright sample of foreground galaxies from the
photometric KiloDegree Survey (KiDS), specifically selected to mimic the
spectroscopic Galaxy And Mass Assembly survey (GAMA). From an independent
sample of KiDS background galaxies, we measure the weak gravitational lensing
profiles of the troughs/ridges. We quantify their lensing strength $A$ as a
function of galaxy density percentile rank $P$ and overdensity $\delta$, and
find that the skewness in the galaxy density distribution is reflected in the
total mass distribution measured by weak lensing. We interpret our results
using the mock galaxy catalogue from the Marenostrum Institut de Ci\`encies de
l'Espai (MICE) Grand Challenge lightcone simulation, and find a good agreement
with our observations. Using signaltonoise weights derived from the Scinet
LIghtCone Simulations (SLICS) mock catalogue we optimally stack the lensing
signal of KiDS troughs with an angular radius $\theta_A$ = {5,10,15,20} arcmin,
resulting in {16.8,14.9,10.13,7.55} $\sigma$ detections. Finally, we select
troughs using a volumelimited sample of galaxies, split into two redshift bins
between 0.1 < z < 0.3. For troughs/ridges with transverse comoving radius $R_A$
= 1.9 Mpc/h, we find no significant difference between the comoving $A'(P)$ and
$A'(\delta)$ relation of the low and highredshift sample. Using the MICE and
SLICS mocks we predict that trough and ridge evolution could be detected with
gravitational lensing using deeper and wider lensing surveys, such as those
from the Large Synoptic Survey Telescope and Euclid.

We present cosmological parameter constraints from a joint analysis of three
cosmological probes: the tomographic cosmic shear signal in $\sim$450 deg$^2$
of data from the Kilo Degree Survey (KiDS), the galaxymatter crosscorrelation
signal of galaxies from the Galaxies And Mass Assembly (GAMA) survey determined
with KiDS weak lensing, and the angular correlation function of the same GAMA
galaxies. We use fast power spectrum estimators that are based on simple
integrals over the realspace correlation functions, and show that they are
practically unbiased over relevant angular frequency ranges. We test our full
pipeline on numerical simulations that are tailored to KiDS and retrieve the
input cosmology. By fitting different combinations of power spectra, we
demonstrate that the three probes are internally consistent. For all probes
combined, we obtain $S_8\equiv \sigma_8 \sqrt{\Omega_{\rm
m}/0.3}=0.800_{0.027}^{+0.029}$, consistent with Planck and the fiducial
KiDS450 cosmic shear correlation function results. Marginalising over wide
priors on the mean of the tomographic redshift distributions yields consistent
results for $S_8$ with an increase of $28\%$ in the error. The combination of
probes results in a $26\%$ reduction in uncertainties of $S_8$ over using the
cosmic shear power spectra alone. The main gain from these additional probes
comes through their constraining power on nuisance parameters, such as the
galaxy intrinsic alignment amplitude or potential shifts in the redshift
distributions, which are up to a factor of two better constrained compared to
using cosmic shear alone, demonstrating the value of largescale structure
probe combination.

We compute the sphericalsky weaklensing power spectrum of the shear and
convergence. We discuss various approximations, such as flatsky, and first
and second order Limber equations for the projection. We find that the impact
of adopting these approximations is negligible when constraining cosmological
parameters from current weak lensing surveys. This is demonstrated using data
from the CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS). We find that
the reported tension with Planck Cosmic Microwave Background (CMB) temperature
anisotropy results cannot be alleviated.
For future largescale surveys with unprecedented precision, we show that the
spherical secondorder Limber approximation will provide sufficient accuracy.
In this case, the cosmicshear power spectrum is shown to be in agreement with
the full projection at the subpercent level for l > 3, with the corresponding
errors an order of magnitude below cosmic variance for all l. When computing
the twopoint shear correlation function, we show that the flatsky fast Hankel
transformation results in errors below two percent compared to the full
spherical transformation. In the spirit of reproducible research, our numerical
implementation of all approximations and the full projection are publicly
available within the package nicaea at
http://www.cosmostat.org/software/nicaea.

We measure the projected galaxy clustering and galaxygalaxy lensing signals
using the Galaxy And Mass Assembly (GAMA) survey and KiloDegree Survey (KiDS)
to study galaxy bias. We use the concept of nonlinear and stochastic galaxy
biasing in the framework of halo occupation statistics to constrain the
parameters of the halo occupation statistics and to unveil the origin of galaxy
biasing. The bias function $\Gamma_{\text{gm}}(r_{\text{p}})$, where
$r_{\text{p}}$ is the projected comoving separation, is evaluated using the
analytical halo model from which the scale dependence of
$\Gamma_{\text{gm}}(r_{\text{p}})$, and the origin of the nonlinearity and
stochasticity in halo occupation models can be inferred. Our observations
unveil the physical reason for the nonlinearity and stochasticity, further
confirmed using hydrodynamical simulations, with the stochasticity mostly
originating from the nonPoissonian behaviour of satellite galaxies in the dark
matter haloes and their spatial distribution, which does not follow the spatial
distribution of dark matter in the halo. The observed nonlinearity is mostly
due to the presence of the central galaxies, as was noted from previous
theoretical work on the same topic. We also see that overall, more massive
galaxies reveal a stronger scale dependence, and out to a larger radius. Our
results show that a wealth of information about galaxy bias is hidden in halo
occupation models. These models should therefore be used to determine the
influence of galaxy bias in cosmological studies.

We derive the essentials of the skewed weak lensing likelihood via a simple
Hierarchical Model. Our likelihood passes four objective and
cosmologyindependent tests which a standard Gaussian likelihood fails. We
demonstrate that sound weak lensing analyses are naturally biased low, and this
does not indicate any new physics such as deviations from $\Lambda$CDM.
Mathematically, the biases arise because noisy twopoint functions follow
skewed distributions. This form of bias is already known from CMB analyses,
where the low multipoles have asymmetric error bars. Weak lensing is more
strongly affected by this asymmetry as galaxies form a discrete set of shear
tracer particles, in contrast to a smooth shear field. We demonstrate that the
biases can be up to 30 percent of the standard deviation per data point,
dependent on the properties of the weak lensing survey. Our likelihood provides
a versatile framework with which to address this bias in future weak lensing
analyses.

This paper is the first of a series of papers constraining cosmological
parameters with weak lensing peak statistics using $\sim 450~\rm deg^2$ of
imaging data from the Kilo Degree Survey (KiDS450). We measure high
signaltonoise ratio (SNR: $\nu$) weak lensing convergence peaks in the range
of $3<\nu<5$, and employ theoretical models to derive expected values. These
models are validated using a suite of simulations. We take into account two
major systematic effects, the boost factor and the effect of baryons on the
massconcentration relation of dark matter haloes. In addition, we investigate
the impacts of other potential astrophysical systematics including the
projection effects of large scale structures, intrinsic galaxy alignments, as
well as residual measurement uncertainties in the shear and redshift
calibration. Assuming a flat $\Lambda$CDM model, we find constraints for
$S_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.5}=0.746^{+0.046}_{0.107}$
according to the degeneracy direction of the cosmic shear analysis and
$\Sigma_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm
m}/0.3)^{0.38}=0.696^{+0.048}_{0.050}$ based on the derived degeneracy
direction of our highSNR peak statistics. The difference between the power
index of $S_{\rm 8}$ and in $\Sigma_{\rm 8}$ indicates that combining the two
probes has the potential to break the degeneracy in $\sigma_{\rm 8}$ and
$\Omega_{\rm m}$. Our results are consistent with the cosmic shear tomographic
correlation analysis of the same dataset and $\sim 2\sigma$ lower than the
Planck 2016 results.

We study the statistics of peaks in a weak lensing reconstructed mass map of
the first 450 square degrees of the Kilo Degree Survey. The map is computed
with aperture masses directly applied to the shear field with an NFWlike
compensated filter. We compare the peak statistics in the observations with
that of simulations for various cosmologies to constrain the cosmological
parameter $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$, which probes the
($\Omega_{\rm m}, \sigma_8$) plane perpendicularly to its main degeneracy. We
estimate $S_8=0.750\pm0.059$, using peaks in the signaltonoise range $0 \leq
{\rm S/N} \leq 4$, and accounting for various systematics, such as
multiplicative shear bias, mean redshift bias, baryon feedback, intrinsic
alignment, and shearposition coupling. These constraints are $\sim25\%$
tighter than the constraints from the high significance peaks alone ($3 \leq
{\rm S/N} \leq 4$) which typically trace singlemassive halos. This
demonstrates the gain of information from lowS/N peaks which correspond to the
projection of several smallmass halos along the lineofsight. Our results are
in good agreement with the tomographic shear twopoint correlation function
measurement in KiDS450. Combining shear peaks with nontomographic
measurements of the shear twopoint correlation functions yields an $\sim20\%$
improvement in the uncertainty on $S_8$ compared to the shear twopoint
correlation functions alone, highlighting the great potential of peaks as a
cosmological probe.

In recent years, multiple studies have reported substantial populations of
large, lowsurfacebrightness galaxies in local galaxy clusters. Various
theories that aim to explain the presence of such ultradiffuse galaxies (UDGs)
have since been proposed. A key question that will help to differentiate
between models is whether UDGs have counterparts in lowermass host haloes, and
what their abundance as a function of halo mass is. In this study we extend our
previous study of UDGs in galaxy clusters to galaxy groups. We measure the
abundance of UDGs in 325 spectroscopicallyselected groups from the Galaxy And
Mass Assembly (GAMA) survey. We make use of the overlapping imaging from the
ESO KiloDegree Survey (KiDS), from which we can identify galaxies with mean
surface brightnesses within their effective radii down to ~25.5 mag
arcsec$^{2}$ in the rband. We are able to measure a significant overdensity
of UDGs (with sizes r_eff > 1.5 kpc) in galaxy groups down to M200=10^12 Msun,
a regime where approximately only 1 in 10 groups contains a UDG that we can
detect. We combine measurements of the abundance of UDGs in haloes that cover
three orders of magnitude in halo mass, finding that their numbers scale quite
steeply with halo mass; N_UDG (R<R200) $\propto$ M200^(1.11+/0.07). To better
interpret this, we also measure the massrichness relation for brighter
galaxies down to $M^*_r$+2.5 in the same GAMA groups, and find a much shallower
relation of N_Bright (R<R200) $\propto$ M200^(0.78+/0.05). This shows that
UDGs are relatively more abundant, compared to bright galaxies, in massive
clusters than in groups. We discuss implications, but whether this difference
is related to a higher destruction rate of UDGs in groups, or whether massive
haloes have a positive effect on their formation, is not yet clear.

We perform a combined analysis of cosmic shear tomography, galaxygalaxy
lensing tomography, and redshiftspace multipole power spectra (monopole and
quadrupole) using 450 deg$^2$ of imaging data by the Kilo Degree Survey (KiDS)
overlapping with two spectroscopic surveys: the 2degree Field Lensing Survey
(2dFLenS) and the Baryon Oscillation Spectroscopic Survey (BOSS). We restrict
the galaxygalaxy lensing and multipole power spectrum measurements to the
overlapping regions with KiDS, and selfconsistently compute the full
covariance between the different observables using a large suite of $N$body
simulations. We methodically analyze different combinations of the observables,
finding that galaxygalaxy lensing measurements are particularly useful in
improving the constraint on the intrinsic alignment amplitude (by 30%, positive
at $3.5\sigma$ in the fiducial data analysis), while the multipole power
spectra are useful in tightening the constraints along the lensing degeneracy
direction (e.g. factor of two stronger matter density constraint in the
fiducial analysis). The fully combined constraint on $S_8 \equiv \sigma_8
\sqrt{\Omega_{\rm m}/0.3} = 0.742 \pm 0.035$, which is an improvement by 20%
compared to KiDS alone, corresponds to a $2.6\sigma$ discordance with Planck,
and is not significantly affected by fitting to a more conservative set of
scales. Given the tightening of the parameter space, we are unable to resolve
the discordance with an extended cosmology that is simultaneously favored in a
model selection sense, including the sum of neutrino masses, curvature,
evolving dark energy, and modified gravity. The complementarity of our
observables allow for constraints on modified gravity degrees of freedom that
are not simultaneously bounded with either probe alone, and up to a factor of
three improvement in the $S_8$ constraint in the extended cosmology compared to
KiDS alone.

We test extensions to the standard cosmological model with weak gravitational
lensing tomography using 450 deg$^2$ of imaging data from the Kilo Degree
Survey (KiDS). In these extended cosmologies, which include massive neutrinos,
nonzero curvature, evolving dark energy, modified gravity, and running of the
scalar spectral index, we also examine the discordance between KiDS and cosmic
microwave background measurements from Planck. The discordance between the two
datasets is largely unaffected by a more conservative treatment of the lensing
systematics and the removal of angular scales most sensitive to nonlinear
physics. The only extended cosmology that simultaneously alleviates the
discordance with Planck and is at least moderately favored by the data includes
evolving dark energy with a timedependent equation of state (in the form of
the $w_0w_a$ parameterization). In this model, the respective $S_8 = \sigma_8
\sqrt{\Omega_{\rm m}/0.3}$ constraints agree at the $1\sigma$ level, and there
is `substantial concordance' between the KiDS and Planck datasets when
accounting for the full parameter space. Moreover, the Planck constraint on the
Hubble constant is wider than in LCDM and in agreement with the Riess et al.
(2016) direct measurement of $H_0$. The dark energy model is moderately favored
as compared to LCDM when combining the KiDS and Planck measurements, and
remains moderately favored after including an informative prior on the Hubble
constant. In both of these scenarios, marginalized constraints in the $w_0w_a$
plane are discrepant with a cosmological constant at the $3\sigma$ level.
Moreover, KiDS constrains the sum of neutrino masses to 4.0 eV (95% CL), finds
no preference for time or scale dependent modifications to the metric
potentials, and is consistent with flatness and no running of the spectral
index. The analysis code is public at https://github.com/sjoudaki/kids450

We constrain the average halo ellipticity of ~2 600 galaxy groups from the
Galaxy And Mass Assembly (GAMA) survey, using the weak gravitational lensing
signal measured from the overlapping Kilo Degree Survey (KiDS). To do so, we
quantify the azimuthal dependence of the stacked lensing signal around seven
different proxies for the orientation of the dark matter distribution, as it is
a priori unknown which one traces the orientation best. On small scales, the
major axis of the brightest group/cluster member (BCG) provides the best proxy,
leading to a clear detection of an anisotropic signal. In order to relate that
to a halo ellipticity, we have to adopt a model density profile. We derive new
expressions for the quadrupole moments of the shear field given an elliptical
model surface mass density profile. Modeling the signal with an elliptical
NavarroFrenkWhite (NFW) profile on scales < 250 kpc, which roughly
corresponds to half the virial radius, and assuming that the BCG is perfectly
aligned with the dark matter, we find an average halo ellipticity of e_h=0.38
+/ 0.12. This agrees well with results from colddarkmatteronly simulations,
which typically report values of e_h ~ 0.3. On larger scales, the lensing
signal around the BCGs does not trace the dark matter distribution well, and
the distribution of group satellites provides a better proxy for the halo's
orientation instead, leading to a 34 sigma detection of a nonzero halo
ellipticity at scales between 250 kpc and 750 kpc. Our results suggest that the
distribution of stars enclosed within a certain radius forms a good proxy for
the orientation of the dark matter within that radius, which has also been
observed in hydrodynamical simulations.

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

We investigate possible signatures of halo assembly bias for
spectroscopically selected galaxy groups from the GAMA survey using weak
lensing measurements from the spatially overlapping regions of the deeper,
highimagingquality photometric KiDS survey. We use GAMA groups with an
apparent richness larger than 4 to identify samples with comparable mean host
halo masses but with a different radial distribution of satellite galaxies,
which is a proxy for the formation time of the haloes. We measure the weak
lensing signal for groups with a steeper than average and with a shallower than
average satellite distribution and find no sign of halo assembly bias, with the
bias ratio of $0.85^{+0.37}_{0.25}$, which is consistent with the $\Lambda$CDM
prediction. Our galaxy groups have typical masses of $10^{13} M_{\odot}/h$,
naturally complementing previous studies of halo assembly bias on galaxy
cluster scales.

We present the tomographic crosscorrelation between galaxy lensing measured
in the Kilo Degree Survey (KiDS450) with overlapping lensing measurements of
the cosmic microwave background (CMB), as detected by Planck 2015. We compare
our joint probe measurement to the theoretical expectation for a flat
$\Lambda$CDM cosmology, assuming the bestfitting cosmological parameters from
the KiDS450 cosmic shear and Planck CMB analyses. We find that our results are
consistent within $1\sigma$ with the KiDS450 cosmology, with an amplitude
rescaling parameter $A_{\rm KiDS} = 0.86 \pm 0.19$. Adopting a Planck
cosmology, we find our results are consistent within $2\sigma$, with $A_{\it
Planck} = 0.68 \pm 0.15$. We show that the agreement is improved in both cases
when the contamination to the signal by intrinsic galaxy alignments is
accounted for, increasing $A$ by $\sim 0.1$. This is the first tomographic
analysis of the galaxy lensing  CMB lensing crosscorrelation signal, and is
based on five photometric redshift bins. We use this measurement as an
independent validation of the multiplicative shear calibration and of the
calibrated source redshift distribution at high redshifts. We find that
constraints on these two quantities are strongly correlated when obtained from
this technique, which should therefore not be considered as a standalone
competitive calibration tool.

We measure the crosscorrelation between FermiLAT gammaray photons and over
1000 deg$^2$ of weak lensing data from the CanadaFranceHawaii Telescope
Lensing Survey (CFHTLenS), the Red Cluster Sequence Lensing Survey (RCSLenS),
and the Kilo Degree Survey (KiDS). We present the first measurement of
tomographic weak lensing crosscorrelations and the first application of
spectral binning to crosscorrelations between gamma rays and weak lensing. The
measurements are performed using an angular power spectrum estimator while the
covariance is estimated using an analytical prescription. We verify the
accuracy of our covariance estimate by comparing it to two internal covariance
estimators. Based on the nondetection of a crosscorrelation signal, we derive
constraints on weakly interacting massive particle (WIMP) dark matter. We
compute exclusion limits on the dark matter annihilation crosssection
$\langle\sigma_\rm{ann} v \rangle$, decay rate $\Gamma_\rm{dec}$, and particle
mass $m_\rm{DM}$. We find that in the absence of a crosscorrelation signal,
tomography does not significantly improve the constraining power of the
analysis. Assuming a strong contribution to the gammaray flux due to
smallscale clustering of dark matter and accounting for known astrophysical
sources of gamma rays, we exclude the thermal relic crosssection for particle
masses of $m_\rm{DM}\lesssim 20$ GeV.

We determine the accuracy of galaxy redshift distributions as estimated from
photometric redshift probability distributions $p(z)$. Our method utilises
measurements of the angular crosscorrelation between photometric galaxies and
an overlapping sample of galaxies with spectroscopic redshifts. We describe the
redshift leakage from a galaxy photometric redshift bin $j$ into a
spectroscopic redshift bin $i$ using the sum of the $p(z)$ for the galaxies
residing in bin $j$. We can then predict the angular crosscorrelation between
photometric and spectroscopic galaxies due to intrinsic galaxy clustering when
$i \neq j$ as a function of the measured angular crosscorrelation when $i=j$.
We also account for enhanced clustering arising from lensing magnification
using a halo model. The comparison of this prediction with the measured signal
provides a consistency check on the validity of using the summed $p(z)$ to
determine galaxy redshift distributions in cosmological analyses, as advocated
by the CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS). We present an
analysis of the photometric redshifts measured by CFHTLenS, which overlaps the
Baryon Oscillation Spectroscopic Survey (BOSS). We also analyse the
Redsequence Cluster Lensing Survey (RCSLenS), which overlaps both BOSS and the
WiggleZ Dark Energy Survey. We find that the summed $p(z)$ from both surveys
are generally biased with respect to the true underlying distributions. If
unaccounted for, this bias would lead to errors in cosmological parameter
estimation from CFHTLenS by less than $\sim 4\%$. For photometric redshift bins
which spatially overlap in 3D with our spectroscopic sample, we determine
redshift bias corrections which can be used in future cosmological analyses
that rely on accurate galaxy redshift distributions.

We present a reanalysis of the CFHTLenS weak gravitational lensing survey
using Complete Orthogonal Sets of E/Bmode Integrals, known as COSEBIs. COSEBIs
provide a complete set of functions to efficiently separate Emodes from
Bmodes and hence allow for robust and stringent tests for systematic errors in
the data. This analysis reveals significant Bmodes on large angular scales
that were not previously seen using the standard E/B decomposition analyses. We
find that the significance of the Bmodes is enhanced when the data is split by
galaxy type and analysed in tomographic redshift bins. Adding tomographic bins
to the analysis increases the number of COSEBIs modes, which results in a less
accurate estimation of the covariance matrix from a set of simulations. We
therefore also present the first compressed COSEBIs analysis of survey data,
where the COSEBIs modes are optimally combined based on their sensitivity to
cosmological parameters. In this tomographic CCOSEBIs analysis we find the
Bmodes to be consistent with zero when the full range of angular scales are
considered.

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

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

We present a new training set for estimating empirical photometric redshifts
of galaxies, which was created as part of the 2dFLenS project. This training
set is located in a 700 sq deg area of the KiDS South field and is randomly
selected and nearly complete at r<19.5. We investigate the photometric redshift
performance obtained with ugriz photometry from VSTATLAS and W1/W2 from WISE,
based on several empirical and template methods. The best redshift errors are
obtained with kerneldensity estimation, as are the lowest biases, which are
consistent with zero within statistical noise. The 68th percentiles of the
redshift scatter for magnitudelimited samples at r<(15.5, 17.5, 19.5) are
(0.014, 0.017, 0.028). In this magnitude range, there are no known ambiguities
in the colourredshift map, consistent with a small rate of redshift outliers.
In the fainter regime, the KDE method produces p(z) estimates per galaxy that
represent unbiased and accurate redshift frequency expectations. The p(z) sum
over any subsample is consistent with the true redshift frequency plus Poisson
noise. Further improvements in redshift precision at r<20 would mostly be
expected from filter sets with narrower passbands to increase the sensitivity
of colours to small changes in redshift.

We present high signaltonoise galaxygalaxy lensing measurements of the
BOSS CMASS sample using 250 square degrees of weak lensing data from CFHTLenS
and CS82. We compare this signal with predictions from mock catalogs trained to
match observables including the stellar mass function and the projected and two
dimensional clustering of CMASS. We show that the clustering of CMASS, together
with standard models of the galaxyhalo connection, robustly predicts a lensing
signal that is 2040% larger than observed. Detailed tests show that our
results are robust to a variety of systematic effects. Lowering the value of
$S_{\rm 8}=\sigma_{\rm 8} \sqrt{\Omega_{\rm m}/0.3}$ compared to Planck2015
reconciles the lensing with clustering. However, given the scale of our
measurement ($r<10$ $h^{1}$ Mpc), other effects may also be at play and need
to be taken into consideration. We explore the impact of baryon physics,
assembly bias, massive neutrinos, and modifications to general relativity on
$\Delta\Sigma$ and show that several of these effects may be nonnegligible
given the precision of our measurement. Disentangling cosmological effects from
the details of the galaxyhalo connection, the effects of baryons, and massive
neutrinos, is the next challenge facing joint lensing and clustering analyses.
This is especially true in the context of large galaxy samples from Baryon
Acoustic Oscillation surveys with precise measurements but complex selection
functions.

We develop a statistical estimator to infer the redshift probability
distribution of a photometric sample of galaxies from its angular
crosscorrelation in redshift bins with an overlapping spectroscopic sample.
This estimator is a minimum variance weighted quadratic function of the data: a
quadratic estimator. This extends and modifies the methodology presented by
McQuinn & White (2013). The derived source redshift distribution is degenerate
with the source galaxy bias, which must be constrained via additional
assumptions. We apply this estimator to constrain source galaxy redshift
distributions in the KiloDegree imaging survey through crosscorrelation with
the spectroscopic 2degree Field Lensing Survey, presenting results first as a
binned stepwise distribution in the range z < 0.8, and then building a
continuous distribution using a Gaussian process model. We demonstrate the
robustness of our methodology using mock catalogues constructed from Nbody
simulations, and comparisons with other techniques for inferring the redshift
distribution.

We investigate the impact of astrophysical systematics on cosmic shear
cosmological parameter constraints from the CanadaFranceHawaii Telescope
Lensing Survey (CFHTLenS), and the concordance with cosmic microwave background
measurements by Planck. We present updated CFHTLenS cosmic shear tomography
measurements extended to degree scales using a covariance calibrated by a new
suite of Nbody simulations. We analyze these measurements with a new model
fitting pipeline, accounting for key systematic uncertainties arising from
intrinsic galaxy alignments, baryonic effects in the nonlinear matter power
spectrum, and photometric redshift uncertainties. We examine the impact of the
systematic degrees of freedom on the cosmological parameter constraints, both
independently and jointly. When the systematic uncertainties are considered
independently, the intrinsic alignment amplitude is the only degree of freedom
that is substantially preferred by the data. When the systematic uncertainties
are considered jointly, there is no consistently strong preference in favor of
the more complex models. We quantify the level of concordance between the
CFHTLenS and Planck datasets by employing two distinct data concordance tests,
grounded in Bayesian evidence and information theory. We find that the two data
concordance tests largely agree with one another, and that the level of
concordance between the CFHTLenS and Planck datasets is sensitive to the exact
details of the systematic uncertainties included in our analysis, ranging from
decisive discordance to substantial concordance as the treatment of the
systematic uncertainties becomes more conservative. The least conservative
scenario is the one most favored by the cosmic shear data, but it is also the
one that shows the greatest degree of discordance with Planck. The data and
analysis code are public at https://github.com/sjoudaki/cfhtlens_revisited

We present measurements of the spatial mapping between (hot) baryons and the
total matter in the Universe, via the crosscorrelation between the thermal
SunyaevZeldovich (tSZ) map from Planck and the weak gravitational lensing maps
from the Red Sequence Cluster Survey (RCSLenS). The crosscorrelations are
performed on the map level where all the sources (including diffuse
intergalactic gas) contribute to the signal. We consider two
configurationspace correlation function estimators, $\xi^{ y\kappa}$ and
$\xi^ {y\gamma_{t}}$, and a Fourier space estimator, $C_{\ell}^{y\kappa}$, in
our analysis. We detect a significant correlation out to three degrees of
angular separation on the sky. Based on statistical noise only, we can report
13$\sigma$ and 17$\sigma$ detections of the crosscorrelation using the
configurationspace $y\kappa$ and $y\gamma_{t}$ estimators, respectively.
Including a heuristic estimate of the sampling variance yields a detection
significance of 6$\sigma$ and 8$\sigma$, respectively. A similar level of
detection is obtained from the Fourierspace estimator, $C_{\ell}^{y\kappa}$.
As each estimator probes different dynamical ranges, their combination improves
the significance of the detection. We compare our measurements with predictions
from the cosmoOWLS suite of cosmological hydrodynamical simulations, where
different galactic feedback models are implemented. We find that a model with
considerable AGN feedback that removes large quantities of hot gas from galaxy
groups and WMAP7yr bestfit cosmological parameters provides the best match to
the measurements. All baryonic models in the context of a Planck cosmology
overpredict the observed signal. Similar cosmological conclusions are drawn
when we employ a halo model with the observed `universal' pressure profile.