
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 demonstrate that deep goodseeing VLT/HAWKI $K_\mathrm{s}$ images
complemented with $g$+$z$band photometry can yield a sensitivity for weak
lensing studies of massive galaxy clusters at redshifts \mbox{$0.7\lesssim z
\lesssim 1.1$}, which is almost identical to the sensitivity of HST/ACS mosaics
of singleorbit depth. Key reasons for this good performance are the excellent
image quality frequently achievable for $K_\mathrm{s}$ imaging from the ground,
a highly effective photometric selection of background galaxies, and a galaxy
ellipticity dispersion that is noticeably lower than for optically observed
highredshift galaxy samples. Incorporating results from the 3DHST and
UltraVISTA surveys we also obtained a more accurate calibration of the source
redshift distribution than previously achieved for similar optical weak lensing
data sets. Here we studied the extremely massive galaxy cluster
RCS2$J$232727.7$$020437 (\mbox{$z=0.699$}), combining deep VLT/\mbox{HAWKI}
$K_\mathrm{s}$ images (point spread function with a 0\farcs35 full width at
half maximum) with LBT/LBC photometry. The resulting weak lensing mass
reconstruction suggests that the cluster consists of a single overdensity,
which is detected with a peak significance of $10.1\sigma$. We constrained the
cluster mass to \mbox{$M_\mathrm{200c}/(10^{15} \mathrm{M}_\odot)
=2.06^{+0.28}_{0.26}(\mathrm{stat.})\pm 0.12 (\mathrm{sys.})$} assuming a
spherical Navarro, Frenk \& White model and simulationbased priors on the
concentration, making it one of the most massive galaxy clusters known in the
\mbox{$z\gtrsim 0.7$} Universe. We also crosschecked the HAWKI measurements
through an analysis of overlapping HST/ACS images, yielding fully consistent
estimates of the lensing signal.

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.

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 build a background cluster candidate catalog from the Next Generation
Virgo Cluster Survey, using our detection algorithm RedGOLD. The NGVS covers
104$deg^2$ of the Virgo cluster in the $u^*,g,r,i,z$bandpasses to a depth of $
g \sim 25.7$~mag (5$\sigma$). Part of the survey was not covered or has shallow
observations in the $r$band. We build two cluster catalogs: one using all
bandpasses, for the fields with deep $r$band observations ($\sim 20 \
deg^2$), and the other using four bandpasses ($u^*,g,i,z$) for the entire NGVS
area. Based on our previous CFHTLS W1 studies, we estimate that both of our
catalogs are $\sim100\%$($\sim70\%$) complete and $\sim80\%$ pure, at $z\le
0.6$($z\lesssim1$), for galaxy clusters with masses of $M\gtrsim10^{14}\
M_{\odot}$. We show that when using four bandpasses, though the photometric
redshift accuracy is lower, RedGOLD detects massive galaxy clusters up to
$z\sim 1$ with completeness and purity similar to the fiveband case. This is
achieved when taking into account the bias in the richness estimation, which is
$\sim40\%$ lower at $0.5\le z<0.6$ and $\sim20\%$ higher at $0.6<z< 0.8$, with
respect to the fiveband case. RedGOLD recovers all the Xray clusters in the
area with mass $M_{500} > 1.4 \times 10^{14} \rm M_{\odot}$ and $0.08<z<0.5$.
Because of our different cluster richness limits and the NGVS depth, our
catalogs reach to lower masses than the published redMaPPer cluster catalog
over the area, and we recover $\sim 90100\%$ of its detections.

We set out to quantify the number density of quiescent massive compact
galaxies at intermediate redshifts. We determine structural parameters based on
iband imaging using the CFHT equatorial SDSS Stripe 82 (CS82) survey (~170 sq.
degrees) taking advantage of an exquisite median seeing of ~0.6''. We select
compact massive (M > 5x10^10 M_sun) galaxies within the redshift range of
0.2<z<0.6. The large volume sampled allows to decrease the effect of cosmic
variance that has hampered the calculation of the number density for this
enigmatic population in many previous studies. We undertake an exhaustive
analysis in an effort to untangle the various findings inherent to the diverse
definition of compactness present in the literature. We find that the absolute
number of compact galaxies is very dependent on the adopted definition and can
change up to a factor of >10. We systematically measure a factor of ~5 more
compacts at the same redshift than what was previously reported on smaller
fields with HST imaging, which are more affected by cosmic variance. This means
that the decrease in number density from z ~ 1.5 to z ~ 0.2 might be only of a
factor of ~25, significantly smaller than what previously reported. This
supports progenitor bias as the main contributor to the size evolution. This
milder decrease is roughly compatible with the predictions from recent
numerical simulations. Only the most extreme compact galaxies, with Reff <
1.5x( M/10^11 M_sun)^0.75 and M > 10^10.7 M_sun, appear to drop in number by a
factor of ~20 and hence likely experience a noticeable size evolution.

We measured stacked weak lensing cluster masses for a sample of 1325 galaxy
clusters detected by the RedGOLD algorithm in the CanadaFranceHawaii
Telescope Legacy Survey W1 and the Next Generation Virgo Cluster Survey at
$0.2<z<0.5$, in the optical richness range $10<\lambda<70$. After a selection
of our best richness subsample ($20<\lambda<50$), this is the most
comprehensive lensing study of a $\sim 100\%$ complete and $\sim 90\%$ pure
optical cluster catalogue in this redshift range, with a total of 346 clusters
in $\sim164~deg^2$. We test three different mass models, and our best model
includes a basic halo model, with a Navarro Frenk and White profile, and
correction terms that take into account cluster miscentering, nonweak shear,
the twohalo term, the contribution of the Brightest Cluster Galaxy, and an a
posteriori correction for the intrinsic scatter in the massrichness relation.
With this model, we obtain a massrichness relation of $\log{M_{\rm
200}/M_{\odot}}=(14.48\pm0.04)+(1.14\pm0.23)\log{(\lambda/40)}$ (statistical
uncertainties). This result is consistent with other published lensing
massrichness relations. When compared to Xray masses and mass proxies, we
find that on average weak lensing masses are $\sim 10\%$ higher than those
derived in the Xray in the range $2\times10^{13}M_{\rm \odot}<E(z) M^{X}_{\rm
200}<2\times10^{14}M_{\rm \odot}$, in agreement with most previous results and
simulations. We also give the coefficients of the scaling relations between the
lensing mass and Xray mass proxies, $L_X$ and $T_X$, and compare them with
previous results.

RedGOLD searches for redsequence galaxy overdensities while minimizing
contamination from dusty starforming galaxies. It imposes an NFW profile and
calculates cluster detection significance and richness. We optimize these
latter two parameters using both simulations and Xray detected cluster
catalogs, and obtain a catalog $\sim 80\%$ pure up to $z \sim 1$, and $\sim
100\%$ ($\sim 70\%$) complete at $z\le 0.6$ ( $z\lesssim1$) for galaxy clusters
with $M \gtrsim 10^{14}\ {\rm M_{\odot}}$ at the CFHTLS Wide depth. In the
CFHTLS W1, we detect 11 cluster candidates per $\rm deg^2$ out to $z\sim1.1$.
When we optimize both completeness and purity, RedGOLD obtains a cluster
catalog with higher completeness and purity than other public catalogs,
obtained using CFHTLS W1 observations, for $M \gtrsim 10^{14}\ {\rm
M_{\odot}}$. We use Xray detected cluster samples to extend the study of the
Xray temperatureoptical richness relation to a lower mass threshold, and find
a mass scatter at fixed richness of $\sigma_{lnM\lambda}=0.39\pm0.07$ and
$\sigma_{lnM\lambda}=0.30\pm0.13$ for the Gozaliasl et al. (2014) and Mehrtens
et al. (2012) samples. When considering similar mass ranges as previous work,
we recover a smaller scatter in mass at fixed richness. We recover $93\%$ of
the redMaPPer detections, and find that its richness estimates is on average
$\sim 4050\%$ larger than ours at $z>0.3$. RedGOLD recovers Xray cluster
spectroscopic redshifts at better than $5\%$ up to $z\sim1$, and the centers
within a few tens of arcseconds.

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 present ThewiZZ, an open source and userfriendly software for estimating
the redshift distributions of photometric galaxies with unknown redshifts by
spatially crosscorrelating them against a reference sample with known
redshifts. The main benefit of ThewiZZ is in separating the angular pair
finding and correlation estimation from the computation of the output
clustering redshifts allowing anyone to create a clustering redshift for their
sample without the intervention of an "expert". It allows the end user of a
given survey to select any subsample of photometric galaxies with unknown
redshifts, match this sample's catalog indices into a valueadded data file,
and produce a clustering redshift estimation for this sample in a fraction of
the time it would take to run all the angular correlations needed to produce a
clustering redshift. We show results with this software using photometric data
from the KiloDegree Survey (KiDS) and spectroscopic redshifts from the Galaxy
and Mass Assembly (GAMA) survey and the Sloan Digital Sky Survey (SDSS). The
results we present for KiDS are consistent with the redshift distributions used
in a recent cosmic shear analysis from the survey. We also present results
using a hybrid machine learningclustering redshift analysis that enables the
estimation of clustering redshifts for individual galaxies. ThewiZZ can be
downloaded at http://github.com/morriscb/ThewiZZ/.

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.

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.