
The CanadaFrance Imaging Survey (CFIS) will map the northern high Galactic
latitude sky in the $u$band ("CFISu", 10,000$\, {\rm deg^2}$) and in the
$r$band ("CFISr", 5,000$\, {\rm deg^2}$), enabling a host of standalone
science investigations, and providing some of the groundbased data necessary
for photometric redshift determination for the Euclid mission. In this first
contribution we present the $u$band component of the survey, describe the
observational strategy, and discuss some first highlight results, based on
approximately one third of the final area. We show that the Galactic anticenter
structure is distributed continuously along the line of sight, out to beyond 20
kpc, and possesses a metallicity distribution that is essentially identical to
that of the outer disk sampled by APOGEE. This suggests that it is probably a
buckled disk of old metalrich stars, rather than a stream or a flare. We also
discuss the future potential for CFISu in discovering starforming dwarf
galaxies around the Local Group, the characterization of the white dwarf and
blue straggler population of the Milky Way, as well as its sensitivity to
lowsurface brightness structures in external galaxies.

We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm
deg^2}$ of $u$band photometry taken as part of the CanadaFrance Imaging
Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the
Northern sky. By combing the CFHT $u$band photometry together with SDSS and
PanSTARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably
the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally
obtain good photometric distance estimates. This survey thus permits the
measurement of metallicities and distances of the dominant mainsequence
population out to approximately 30 kpc, and provides much higher number of
stars at large extraplanar distances than have been available from previous
surveys. We develop a nonparametric distancemetallicity decomposition
algorithm and apply it to the sky at $30\deg < b < 70\deg$ and to the North
Galactic Cap. We find that the metallicitydistance distribution is
wellrepresented by three populations whose metallicity distributions do not
vary significantly with vertical height above the disk. As traced in
mainsequence stars, the stellar halo component shows a vertical density
profile that is close to exponential, with a scale height of around 3 kpc. This
may indicate that the inner halo was formed partly from disk stars ejected in
an ancient minor merger.

The MICROSCOPE space mission, launched on April 25, 2016, aims to test the
weak equivalence principle (WEP) with a 10^15 precision. To reach this
performance requires an accurate and robust data analysis method, especially
since the possible WEP violation signal will be dominated by a strongly colored
noise. An important complication is brought by the fact that some values will
be missing therefore, the measured time series will not be strictly regularly
sampled. Those missing values induce a spectral leakage that significantly
increases the noise in Fourier space, where the WEP violation signal is looked
for, thereby complicating scientific returns. Recently, we developed an
inpainting algorithm to correct the MICROSCOPE data for missing values. This
code has been integrated in the official MICROSCOPE data processing pipeline
because it enables us to significantly measure an equivalence principle
violation (EPV) signal in a modelindependent way, in the inertial satellite
configuration. In this work, we present several improvements to the method that
may allow us now to reach the MICROSCOPE requirements for both inertial and
spin satellite configurations. The main improvement has been obtained using a
prior on the power spectrum of the colorednoise that can be directly derived
from the incomplete data. We show that after reconstructing missing values with
this new algorithm, a leastsquares fit may allow us to significantly measure
an EPV signal with a 0.96x10^15 precision in the inertial mode and 1.2x10^15
precision in the spin mode. Although, the inpainting method presented in this
paper has been optimized to the MICROSCOPE data, it remains sufficiently
general to be used in the general context of missing data in time series
dominated by an unknown colorednoise. The improved inpainting software, called
ICON, is freely available at http://www.cosmostat.org/software/icon.

This is the third in a series of papers that develop a new and flexible model
to predict weaklensing (WL) peak counts, which have been shown to be a very
valuable nonGaussian probe of cosmology. In this paper, we compare the
cosmological information extracted from WL peak counts using different
filtering techniques of the galaxy shear data, including linear filtering with
a Gaussian and two compensated filters (the starlet wavelet and the aperture
mass), and the nonlinear filtering method MRLens. We present improvements to
our model that account for realistic survey conditions, which are masks,
sheartoconvergence transformations, and nonconstant noise. We create
simulated peak counts from our stochastic model, from which we obtain
constraints on the matter density $\Omega_\mathrm{m}$, the power spectrum
normalisation $\sigma_8$, and the darkenergy parameter $w_0$. We use two
methods for parameter inference, a copula likelihood, and approximate Bayesian
computation (ABC). We measure the contour width in the
$\Omega_\mathrm{m}$$\sigma_8$ degeneracy direction and the figure of merit to
compare parameter constraints from different filtering techniques. We find that
starlet filtering outperforms the Gaussian kernel, and that including peak
counts from different smoothing scales helps to lift parameter degeneracies.
Peak counts from different smoothing scales with a compensated filter show very
little crosscorrelation, and adding information from different scales can
therefore strongly enhance the available information. Measuring peak counts
separately from different scales yields tighter constraints than using a
combined peak histogram from a single map that includes multiscale information.
Our results suggest that a compensated filter function with counts included
separately from different smoothing scales yields the tightest constraints on
cosmological parameters from WL peaks.

We propose a new massmapping algorithm, specifically designed to recover
smallscale information from a combination of gravitational shear and flexion.
Including flexion allows us to supplement the shear on small scales in order to
increase the sensitivity to substructures and the overall resolution of the
convergence map without relying on strong lensing constraints. In order to
preserve all available small scale information, we avoid any binning of the
irregularly sampled input shear and flexion fields and treat the massmapping
problem as a general illposed inverse problem, regularised using a robust
multiscale wavelet sparsity prior. The resulting algorithm incorporates
redshift, reduced shear, and reduced flexion measurements for individual
galaxies and is made highly efficient by the use of fast Fourier estimators. We
test our reconstruction method on a set of realistic weak lensing simulations
corresponding to typical HST/ACS cluster observations and demonstrate our
ability to recover substructures with the inclusion of flexion which are lost
if only shear information is used. In particular, we can detect substructures
at the 15$^{\prime \prime}$ scale well outside of the critical region of the
clusters. In addition, flexion also helps to constrain the shape of the central
regions of the main dark matter halos. Our massmapping software, called
Glimpse2D, is made freely available at
http://www.cosmostat.org/software/glimpse .

Missing data are a common problem in experimental and observational physics.
They can be caused by various sources, either an instrument's saturation, or a
contamination from an external event, or a data loss. In particular, they can
have a disastrous effect when one is seeking to characterize a
colorednoisedominated signal in Fourier space, since they create a spectral
leakage that can artificially increase the noise. It is therefore important to
either take them into account or to correct for them prior to e.g. a
LeastSquare fit of the signal to be characterized. In this paper, we present
an application of the {\it inpainting} algorithm to mock MICROSCOPE data; {\it
inpainting} is based on a sparsity assumption, and has already been used in
various astrophysical contexts; MICROSCOPE is a French Space Agency mission,
whose launch is expected in 2016, that aims to test the Weak Equivalence
Principle down to the $10^{15}$ level. We then explore the {\it inpainting}
dependence on the number of gaps and the total fraction of missing values. We
show that, in a worstcase scenario, after reconstructing missing values with
{\it inpainting}, a LeastSquare fit may allow us to significantly measure a
$1.1\times10^{15}$ Equivalence Principle violation signal, which is
sufficiently close to the MICROSCOPE requirements to implement {\it inpainting}
in the official MICROSCOPE data processing and analysis pipeline. Together with
the previously published KARMA method, {\it inpainting} will then allow us to
independently characterize and crosscheck an Equivalence Principle violation
signal detection down to the $10^{15}$ level.

In asteroseismology, the observed time series often suffers from incomplete
time coverage due to gaps. The presence of periodic gaps may generate spurious
peaks in the power spectrum that limit the analysis of the data. Various
methods have been developed to deal with gaps in time series data. However, it
is still important to improve these methods to be able to extract all the
possible information contained in the data. In this paper, we propose a new
approach to handle the problem, the socalled inpainting method. This
technique, based on a sparsity prior, enables to judiciously fillin the gaps
in the data, preserving the asteroseismic signal, as far as possible. The
impact of the observational window function is reduced and the interpretation
of the power spectrum is simplified. This method is applied both on ground and
spacebased data. It appears that the inpainting technique improves the
oscillation modes detection and estimation. Additionally, it can be used to
study very long time series of many stars because its computation is very fast.
For a time series of 50 days of CoRoTlike data, it allows a speedup factor of
1000, if compared to methods of the same accuracy.

In this paper, we compare three methods to reconstruct galaxy cluster density
fields with weak lensing data. The first method called FLens integrates an
inpainting concept to invert the shear field with possible gaps, and a
multiscale entropy denoising procedure to remove the noise contained in the
final reconstruction, that arises mostly from the random intrinsic shape of the
galaxies. The second and third methods are based on a model of the density
field made of a multiscale grid of radial basis functions. In one case, the
model parameters are computed with a linear inversion involving a singular
value decomposition. In the other case, the model parameters are estimated
using a Bayesian MCMC optimization implemented in the lensing software
Lenstool. Methods are compared on simulated data with varying galaxy density
fields. We pay particular attention to the errors estimated with resampling. We
find the multiscale grid model optimized with MCMC to provide the best
results, but at high computational cost, especially when considering
resampling. The SVD method is much faster but yields noisy maps, although this
can be mitigated with resampling. The FLens method is a good compromise with
fast computation, high signal to noise reconstruction, but lower resolution
maps. All three methods are applied to the MACS J0717+3745 galaxy cluster
field, and reveal the filamentary structure discovered in Jauzac et al. 2012.
We conclude that sensitive priors can help to get high signal to noise, and
unbiased reconstructions.

This paper describes the definition of a typical nextgeneration spacebased
weak gravitational lensing experiment. We first adopt a set of toplevel
science requirements from the literature, based on the scale and depth of the
galaxy sample, and the avoidance of systematic effects in the measurements
which would bias the derived shear values. We then identify and categorise the
contributing factors to the systematic effects, combining them with the correct
weighting, in such a way as to fit within the toplevel requirements. We
present techniques which permit the performance to be evaluated and explore the
limits at which the contributing factors can be managed. Besides the modelling
biases resulting from the use of weighted moments, the main contributing
factors are the reconstruction of the instrument point spread function (PSF),
which is derived from the stellar images on the image, and the correction of
the charge transfer inefficiency (CTI) in the CCD detectors caused by radiation
damage.

The first half of this paper explores the origin of systematic biases in the
measurement of weak gravitational lensing. Compared to previous work, we expand
the investigation of PSF instability and fold in for the first time the effects
of nonidealities in electronic imaging detectors and imperfect galaxy shape
measurement algorithms. Together, these now explain the additive A(l) and
multiplicative M(l) systematics typically reported in current lensing
measurements. We find that overall performance is driven by a product of a
telescope/camera's *absolute performance*, and our *knowledge about its
performance*.
The second half of this paper propagates any residual shear measurement
biases through to their effect on cosmological parameter constraints. Fully
exploiting the statistical power of Stage IV weak lensing surveys will require
additive biases A<1.8e12 and multiplicative biases M<4.0e3. These can be
allocated between individual budgets in hardware, calibration data and
software, using results from the first half of the paper.
If instrumentation is stable and wellcalibrated, we find extant shear
measurement software from GREAT10 already meet requirements on galaxies
detected at S/N=40. Averaging over a population of galaxies with a realistic
distribution of sizes, it also meets requirements for a 2D cosmic shear
analysis from space. If used on fainter galaxies or for 3D cosmic shear
tomography, existing algorithms would need calibration on simulations to avoid
introducing bias at a level similar to the statistical error. Requirements on
hardware and calibration data are discussed in more detail in a companion
paper. Our analysis is intentionally general, but is specifically being used to
drive the hardware and ground segment performance budget for the design of the
European Space Agency's recentlyselected Euclid mission.

The aperture mass statistic is a common tool used in weak lensing studies. By
convolving lensing maps with a filter function of a specific scale, chosen to
be larger than the scale on which the noise is dominant, the lensing signal may
be boosted with respect to the noise. This allows for detection of structures
at increased fidelity. Furthermore, higherorder statistics of the aperture
mass (such as its skewness or kurtosis), or counting of the peaks seen in the
resulting aperture mass maps, provide a convenient and effective method to
constrain the cosmological parameters. In this paper, we more fully explore the
formalism underlying the aperture mass statistic. We demonstrate that the
aperture mass statistic is formally identical to a wavelet transform at a
specific scale. Further, we show that the filter functions most frequently used
in aperture mass studies are not ideal, being nonlocal in both real and
Fourier space. In contrast, the wavelet formalism offers a number of wavelet
functions that are localized both in real and Fourier space, yet similar to the
'optimal' aperture mass filters commonly adopted. Additionally, for a number of
wavelet functions, such as the starlet wavelet, very fast algorithms exist to
compute the wavelet transform. This offers significant advantages over the
usual aperture mass algorithm when it comes to image processing time,
demonstrating speedup factors of ~ 5  1200 for aperture radii in the range 2
to 64 pixels on an image of 1024 x 1024 pixels.

Weak gravitational lensing has become a common tool to constrain the
cosmological model. The majority of the methods to derive constraints on
cosmological parameters use secondorder statistics of the cosmic shear.
Despite their success, secondorder statistics are not optimal and degeneracies
between some parameters remain. Tighter constraints can be obtained if
secondorder statistics are combined with a statistic that is efficient to
capture nonGaussianity. In this paper, we search for such a statistical tool
and we show that there is additional information to be extracted from
statistical analysis of the convergence maps beyond what can be obtained from
statistical analysis of the shear field. For this purpose, we have carried out
a large number of cosmological simulations along the {\sigma}8{\Omega}m
degeneracy, and we have considered three different statistics commonly used for
nonGaussian features characterization: skewness, kurtosis and peak count. To
be able to investigate nonGaussianity directly in the shear field we have used
the aperture mass definition of these three statistics for different scales.
Then, the results have been compared with the results obtained with the same
statistics estimated in the convergence maps at the same scales. First, we show
that shear statistics give similar constraints to those given by convergence
statistics, if the same scale is considered. In addition, we find that the peak
count statistic is the best to capture nonGaussianities in the weak lensing
field and to break the {\sigma}8{\Omega}m degeneracy. We show that this
statistical analysis should be conducted in the convergence maps: first,
because there exist fast algorithms to compute the convergence map for
different scales, and secondly because it offers the opportunity to denoise the
reconstructed convergence map, which improves nonGaussian features extraction.

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses
on a problem that is of crucial importance for future observations in
cosmology. The shapes of distant galaxies can be used to determine the
properties of dark energy and the nature of gravity, because light from those
galaxies is bent by gravity from the intervening dark matter. The observed
galaxy images appear distorted, although only slightly, and their shapes must
be precisely disentangled from the effects of pixelisation, convolution and
noise. The worldwide gravitational lensing community has made significant
progress in techniques to measure these distortions via the Shear TEsting
Program (STEP). Via STEP, we have run challenges within our own community, and
come to recognise that this particular image analysis problem is ideally
matched to experts in statistical inference, inverse problems and computational
learning. Thus, in order to continue the progress seen in recent years, we are
seeking an infusion of new ideas from these communities. This document details
the GREAT08 Challenge for potential participants. Please visit
http://www.great08challenge.info for the latest information.

We have performed a 70 billion darkmatter particles Nbody simulation in a 2
$h^{1}$ Gpc periodic box, using the concordance, cosmological model as favored
by the latest WMAP3 results. We have computed a fullsky convergence map with a
resolution of $\Delta \theta \simeq 0.74$ arcmin$^{2}$, spanning 4 orders of
magnitude in angular dynamical range. Using various highorder statistics on a
realistic cut sky, we have characterized the transition from the linear to the
nonlinear regime at $\ell \simeq 1000$ and shown that realistic galactic
masking affects highorder moments only below $\ell < 200$. Each domain
(Gaussian and nonGaussian) spans 2 decades in angular scale. This map is
therefore an ideal tool for testing mapmaking algorithms on the sphere. As a
first step in addressing the full map reconstruction problem, we have
benchmarked in this paper two denoising methods: 1) Wiener filtering applied to
the Spherical Harmonics decomposition of the map and 2) a new method, called
MRLens, based on the modification of the Maximum Entropy Method on a Wavelet
decomposition. While the latter is optimal on large spatial scales, where the
signal is Gaussian, MRLens outperforms the Wiener method on small spatial
scales, where the signal is highly nonGaussian. The simulated fullsky
convergence map is freely available to the community to help the development of
new mapmaking algorithms dedicated to the next generation of weaklensing
surveys.

Ordinary baryonic particles (such as protons and neutrons) account for only
onesixth of the total matter in the Universe. The remainder is a mysterious
"dark matter" component, which does not interact via electromagnetism and thus
neither emits nor reflects light. As dark matter cannot be seen directly using
traditional observations, very little is currently known about its properties.
It does interact via gravity, and is most effectively probed through
gravitational lensing: the deflection of light from distant galaxies by the
gravitational attraction of foreground mass concentrations. This is a purely
geometrical effect that is free of astrophysical assumptions and sensitive to
all matter  whether baryonic or dark. Here we show high fidelity maps of the
largescale distribution of dark matter, resolved in both angle and depth. We
find a loose network of filaments, growing over time, which intersect in
massive structures at the locations of clusters of galaxies. Our results are
consistent with predictions of gravitationally induced structure formation, in
which the initial, smooth distribution of dark matter collapses into filaments
then into clusters, forming a gravitational scaffold into which gas can
accumulate, and stars can be built.

This paper presents a new method for the reconstruction of weak lensing mass
maps. It uses the multiscale entropy concept, which is based on wavelets, and
the False Discovery Rate which allows us to derive robust detection levels in
wavelet space. We show that this new restoration approach outperforms several
standard techniques currently used for weak shear mass reconstruction. This
method can also be used to separate E and B modes in the shear field, and thus
test for the presence of residual systematic effects. We concentrate on large
blind cosmic shear surveys, and illustrate our results using simulated shear
maps derived from NBody LambdaCDM simulations with added noise corresponding
to both groundbased and spacebased observations.