-
We investigate the use of estimators of weak lensing power spectra based on a
flat-sky implementation of the Pseudo-Cl (PCl) technique, where the masked
shear field is transformed without regard for masked regions of sky. This
masking mixes power, and E-convergence and B-modes. To study the accuracy of
forward-modelling and full-sky power spectrum recovery we consider both
large-area survey geometries, and small-scale masking due to stars and a
checkerboard model for field-of-view gaps. The power spectrum for the
large-area survey geometry is sparsely-sampled and highly oscillatory, which
makes modelling problematic. Instead, we derive an overall calibration for
large-area mask bias using simulated fields. The effects of small-area 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 forward-model of the raw PCls, and recover
the full-sky convergence power to within a few percent accuracy for both
Gaussian and lognormal-distributed shear fields. Propagating this through to
cosmological parameters using a Fisher-Matrix 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 flat-sky PCl analysis is accurate for
current surveys but a Euclid-like 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 galaxy-matter cross-correlation
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 real-space 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
KiDS-450 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 large-scale structure
probe combination.
-
We compute the spherical-sky weak-lensing power spectrum of the shear and
convergence. We discuss various approximations, such as flat-sky, 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 Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We find that
the reported tension with Planck Cosmic Microwave Background (CMB) temperature
anisotropy results cannot be alleviated.
For future large-scale surveys with unprecedented precision, we show that the
spherical second-order Limber approximation will provide sufficient accuracy.
In this case, the cosmic-shear power spectrum is shown to be in agreement with
the full projection at the sub-percent level for l > 3, with the corresponding
errors an order of magnitude below cosmic variance for all l. When computing
the two-point shear correlation function, we show that the flat-sky 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 (KiDS-450). We measure high
signal-to-noise 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
mass-concentration 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 high-SNR 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 NFW-like
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 signal-to-noise 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 shear-position 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 single-massive halos. This
demonstrates the gain of information from low-S/N peaks which correspond to the
projection of several small-mass halos along the line-of-sight. Our results are
in good agreement with the tomographic shear two-point correlation function
measurement in KiDS-450. Combining shear peaks with non-tomographic
measurements of the shear two-point correlation functions yields an $\sim20\%$
improvement in the uncertainty on $S_8$ compared to the shear two-point
correlation functions alone, highlighting the great potential of peaks as a
cosmological probe.
-
We perform a combined analysis of cosmic shear tomography, galaxy-galaxy
lensing tomography, and redshift-space 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 2-degree Field Lensing Survey
(2dFLenS) and the Baryon Oscillation Spectroscopic Survey (BOSS). We restrict
the galaxy-galaxy lensing and multipole power spectrum measurements to the
overlapping regions with KiDS, and self-consistently 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 galaxy-galaxy 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 cross-correlation between galaxy lensing measured
in the Kilo Degree Survey (KiDS-450) 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 best-fitting cosmological parameters from
the KiDS-450 cosmic shear and Planck CMB analyses. We find that our results are
consistent within $1\sigma$ with the KiDS-450 cosmology, with an amplitude
re-scaling 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 cross-correlation 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 stand-alone
competitive calibration tool.
-
We present a re-analysis of the CFHTLenS weak gravitational lensing survey
using Complete Orthogonal Sets of E/B-mode Integrals, known as COSEBIs. COSEBIs
provide a complete set of functions to efficiently separate E-modes from
B-modes and hence allow for robust and stringent tests for systematic errors in
the data. This analysis reveals significant B-modes on large angular scales
that were not previously seen using the standard E/B decomposition analyses. We
find that the significance of the B-modes 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
B-modes to be consistent with zero when the full range of angular scales are
considered.
-
We measure the cross-correlation signature between the Planck CMB lensing map
and the weak lensing observations from both the Red-sequence Cluster Lensing
Survey (RCSLenS) and the Canada-France-Hawai Telescope Lensing Survey
(CFHTLenS). In addition to a Fourier analysis, we include the first
configuration-space 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 configuration-space
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, zero-padding and mask
multiplication, validated on a suite of high-resolution 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-/B-mode 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 two-point
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 top-hat filters in Fourier
space, for which the real-space analog has a finite support, yielding strong
bounds on the accuracy of band power estimates. The error on an estimated
band-power 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-/B-Integrals, is a recently developed statistical
measure that encompasses the complete E-/B-mode 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 higher-dimensional 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 figure-of-merit
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-/B-mode separable second-order 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.