
We investigate the use of estimators of weak lensing power spectra based on a
flatsky implementation of the PseudoCl (PCl) technique, where the masked
shear field is transformed without regard for masked regions of sky. This
masking mixes power, and Econvergence and Bmodes. To study the accuracy of
forwardmodelling and fullsky power spectrum recovery we consider both
largearea survey geometries, and smallscale masking due to stars and a
checkerboard model for fieldofview gaps. The power spectrum for the
largearea survey geometry is sparselysampled and highly oscillatory, which
makes modelling problematic. Instead, we derive an overall calibration for
largearea mask bias using simulated fields. The effects of smallarea star
masks can be accurately corrected for, while the checkerboard mask has
oscillatory and spiky behaviour which leads to percent biases. Apodisation of
the masked fields leads to increased biases and a loss of information. We find
that we can construct an unbiased forwardmodel of the raw PCls, and recover
the fullsky convergence power to within a few percent accuracy for both
Gaussian and lognormaldistributed shear fields. Propagating this through to
cosmological parameters using a FisherMatrix formalism, we find we can make
unbiased estimates of parameters for surveys up to 1,200 deg$^2$ with 30
galaxies per arcmin$^2$, beyond which the percent biases become larger than the
statistical accuracy. This implies a flatsky PCl analysis is accurate for
current surveys but a Euclidlike survey will require higher accuracy.

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.

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.

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

We measure the crosscorrelation signature between the Planck CMB lensing map
and the weak lensing observations from both the Redsequence Cluster Lensing
Survey (RCSLenS) and the CanadaFranceHawai Telescope Lensing Survey
(CFHTLenS). In addition to a Fourier analysis, we include the first
configurationspace detection, based on the estimators $\langle \kappa_{\rm
CMB} \kappa_{\rm gal} \rangle$ and $\langle \kappa_{\rm CMB} \gamma_{t}
\rangle$. Combining 747.2 deg$^2$ from both surveys, we find a detection
significance that exceeds $4.2\sigma$ in both Fourier and configurationspace
analyses. Scaling the predictions by a free parameter $A$, we obtain $A^{\rm
Planck}_{\rm CFHT}= 0.68\pm 0.31 $ and $A^{\rm Planck}_{\rm RCS}= 1.31\pm
0.33$. In preparation for the next generation of measurements similar to these,
we quantify the impact of different analysis choices on these results. First,
since none of these estimators probes the exact same dynamical range, we
improve our detection by combining them. Second, we carry out a detailed
investigation on the effect of apodization, zeropadding and mask
multiplication, validated on a suite of highresolution simulations, and find
that the latter produces the largest systematic bias in the cosmological
interpretation. Finally, we show that residual contamination from intrinsic
alignment and the effect of photometric redshift error are both largely
degenerate with the characteristic signal from massive neutrinos, however the
signature of baryon feedback might be easier to distinguish. The three lensing
datasets are now publicly available.

Future large scale cosmological surveys will provide huge data sets whose
analysis requires efficient data compression. Calculating accurate covariances
is extremely challenging with increasing number of statistics used. Here we
introduce a formalism for achieving efficient data compression, based on a
local expansion of statistical measures around a fiducial cosmological model.
We specifically apply and test this approach for the case of cosmic shear
statistics. We demonstrate the performance of our approach, using a Fisher
analysis on cosmic shear tomography described in terms of E/Bmode separating
statistics (COSEBIs). We show that our data compression is highly effective in
extracting essentially the full cosmological information from a strongly
reduced number of observables. Specifically, the number of statistics needed
decreases by at least one order of magnitude relative to the COSEBIs, which
already compress the data substantially compared to the shear twopoint
correlation functions. The efficiency appears to be affected only slightly if a
highly inaccurate covariance is used for defining the compressed statistics,
showing the robustness of the method. We conclude that an efficient data
compression is achievable and that the number of compressed statistics depends
on the number of model parameters.
In addition, we study how well band powers can be obtained from measuring
shear correlation functions over a finite interval of separations. We show the
strong limitations on the possibility to construct tophat filters in Fourier
space, for which the realspace analog has a finite support, yielding strong
bounds on the accuracy of band power estimates. The error on an estimated
bandpower is larger for a narrower filter and a smaller angular range which
for relevant cases can be as large as 10%.

Context. Gravitational lensing is one of the leading tools in understanding
the dark side of the Universe. The need for accurate, efficient and effective
methods which are able to extract this information along with other
cosmological parameters from cosmic shear data is ever growing. COSEBIs,
Complete Orthogonal Sets of E/BIntegrals, is a recently developed statistical
measure that encompasses the complete E/Bmode separable information contained
in the shear correlation functions measured on a finite angular range. Aims.
The aim of the present work is to test the properties of this newly developed
statistics for a higherdimensional parameter space and to generalize and test
it for shear tomography. Methods. We use Fisher analysis to study the
effectiveness of COSEBIs. We show our results in terms of figureofmerit
quantities, based on Fisher matrices. Results. We find that a relatively small
number of COSEBIs modes is always enough to saturate to the maximum information
level. This number is always smaller for 'logarithmic COSEBIs' than for 'linear
COSEBIs', and also depends on the number of redshift bins, the number and
choice of cosmological parameters, as well as the survey characteristics.
Conclusions. COSEBIs provide a very compact way of analyzing cosmic shear data,
i.e., all the E/Bmode separable secondorder statistical information in the
data is reduced to a small number of COSEBIs modes. Furthermore, with this
method the arbitrariness in data binning is no longer an issue since the
COSEBIs modes are discrete. Finally, the small number of modes also implies
that covariances, and their inverse, are much more conveniently obtainable,
e.g., from numerical simulations, than for the shear correlation functions
themselves.