
As astronomical observations from space benefit from improved sensitivity,
the effectiveness of scientific programs is becoming limited by systematics
that often originate in poorly understood image sensor behavior. Traditional,
bottomup detector characterization methods provide one way to model underlying
detector physics, and generate ever more faithful numerical simulations, but
this approach is vulnerable to preconceptions and oversimplification. The
alternative topdown approach is laboratory emulation, which enables
observation, calibration, and analysis scenarios to be tested without relying
on a complete understanding of the underlying physics. This complements
detector characterization and simulation efforts by testing their validity. We
describe a laboratory facility and experimental testbed that supports the
emulation of a wide range of mission concepts such as gravitational weak
lensing measurements by WFIRST and high precision spectrophotometry of
transiting exoplanets by JWST. An Offner relay projects readily customizable
"scenes" (e.g. stars, galaxies, spectra) with very low optical aberration over
the full area of a typical optical or near infrared image sensor. f/8 and
slower focal ratios may be selected, spanning those of most proposed space
missions and approximating the point spread function (PSF) size of seeing
limited ground based surveys. Diffraction limited PSFs are projected over a
wide field of view and wavelength range to deliver highly predictable image
patterns down to subpixel scales with stable intensity and fine motion
control. The testbed enables realistic validation of detector performance on
sciencelike images, which aids mission design and survey strategy, as well as
targeted investigations of various detector effects.

Dark energy may be the first sign of new fundamental physics in the Universe,
taking either a physical form or revealing a correction to Einsteinian gravity.
Weak gravitational lensing and galaxy peculiar velocities provide complementary
probes of General Relativity, and in combination allow us to test modified
theories of gravity in a unique way. We perform such an analysis by combining
measurements of cosmic shear tomography from the CanadaFrance Hawaii Telescope
Lensing Survey (CFHTLenS) with the growth of structure from the WiggleZ Dark
Energy Survey and the SixdegreeField Galaxy Survey (6dFGS), producing the
strongest existing joint constraints on the metric potentials that describe
general theories of gravity. For scaleindependent modifications to the metric
potentials which evolve linearly with the effective dark energy density, we
find presentday cosmological deviations in the Newtonian potential and
curvature potential from the prediction of General Relativity to be (Delta
Psi)/Psi = 0.05 \pm 0.25 and (Delta Phi)/Phi = 0.05 \pm 0.3 respectively (68
per cent CL).

Interaction of charges in CCDs with the already accumulated charge
distribution causes both a flux dependence of the pointspread function (an
increase of observed size with flux, also known as the brighter/fatter effect)
and pixeltopixel correlations of the Poissonian noise in flat fields. We
describe these effects in the Dark Energy Camera (DECam) with charge dependent
shifts of effective pixel borders, i.e. the Antilogus et al. (2014) model,
which we fit to measurements of flatfield Poissonian noise correlations. The
latter fall off approximately as a powerlaw r^2.5 with pixel separation r,
are isotropic except for an asymmetry in the direct neighbors along rows and
columns, are stable in time, and are weakly dependent on wavelength. They show
variations from chip to chip at the 20% level that correlate with the silicon
resistivity. The charge shifts predicted by the model cause biased shape
measurements, primarily due to their effect on bright stars, at levels
exceeding weak lensing science requirements. We measure the flux dependence of
star images and show that the effect can be mitigated by applying the reverse
charge shifts at the pixel level during image processing. Differences in
stellar size, however, remain significant due to residuals at larger distance
from the centroid.

We present first results from the third GRavitational lEnsing Accuracy
Testing (GREAT3) challenge, the third in a sequence of challenges for testing
methods of inferring weak gravitational lensing shear distortions from
simulated galaxy images. GREAT3 was divided into experiments to test three
specific questions, and included simulated space and groundbased data with
constant or cosmologicallyvarying shear fields. The simplest (control)
experiment included parametric galaxies with a realistic distribution of
signaltonoise, size, and ellipticity, and a complex point spread function
(PSF). The other experiments tested the additional impact of realistic galaxy
morphology, multiple exposure imaging, and the uncertainty about a
spatiallyvarying PSF; the last two questions will be explored in Paper II. The
24 participating teams competed to estimate lensing shears to within systematic
error tolerances for upcoming StageIV dark energy surveys, making 1525
submissions overall. GREAT3 saw considerable variety and innovation in the
types of methods applied. Several teams now meet or exceed the targets in many
of the tests conducted (to within the statistical errors). We conclude that the
presence of realistic galaxy morphology in simulations changes shear
calibration biases by $\sim 1$ per cent for a wide range of methods. Other
effects such as truncation biases due to finite galaxy postage stamps, and the
impact of galaxy type as measured by the S\'{e}rsic index, are quantified for
the first time. Our results generalize previous studies regarding sensitivities
to galaxy size and signaltonoise, and to PSF properties such as seeing and
defocus. Almost all methods' results support the simple model in which additive
shear biases depend linearly on PSF ellipticity.

GALSIM is a collaborative, opensource project aimed at providing an image
simulation tool of enduring benefit to the astronomical community. It provides
a software library for generating images of astronomical objects such as stars
and galaxies in a variety of ways, efficiently handling image transformations
and operations such as convolution and rendering at high precision. We describe
the GALSIM software and its capabilities, including necessary theoretical
background. We demonstrate that the performance of GALSIM meets the stringent
requirements of high precision image analysis applications such as weak
gravitational lensing, for current datasets and for the Stage IV dark energy
surveys of the Large Synoptic Survey Telescope, ESA's Euclid mission, and
NASA's WFIRSTAFTA mission. The GALSIM project repository is public and
includes the full code history, all open and closed issues, installation
instructions, documentation, and wiki pages (including a Frequently Asked
Questions section). The GALSIM repository can be found at
https://github.com/GalSimdevelopers/GalSim .

We present the cluster massrichness scaling relation calibrated by a weak
lensing analysis of >18000 galaxy cluster candidates in the
CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS). Detected using the
3DMatchedFilter clusterfinder of Milkeraitis et al., these cluster
candidates span a wide range of masses, from the small group scale up to
$\sim10^{15} M_{\odot}$, and redshifts 0.2 $\lesssim z\lesssim$ 0.9. The total
significance of the stacked shear measurement amounts to 54$\sigma$. We compare
cluster masses determined using weak lensing shear and magnification, finding
the measurements in individual richness bins to yield 1$\sigma$ compatibility,
but with magnification estimates biased low. This first direct mass comparison
yields important insights for improving the systematics handling of future
lensing magnification work. In addition, we confirm analyses that suggest
cluster miscentring has an important effect on the observed 3DMF halo
profiles, and we quantify this by fitting for projected cluster centroid
offsets, which are typically $\sim$ 0.4 arcmin. We bin the cluster candidates
as a function of redshift, finding similar cluster masses and richness across
the full range up to $z \sim$ 0.9. We measure the 3DMF massrichness scaling
relation $M_{200} = M_0 (N_{200} / 20)^\beta$. We find a normalization $M_0
\sim (2.7^{+0.5}_{0.4}) \times 10^{13} M_{\odot}$, and a logarithmic slope of
$\beta \sim 1.4 \pm 0.1$, both of which are in 1$\sigma$ agreement with results
from the magnification analysis. We find no evidence for a redshiftdependence
of the normalization. The CFHTLenS 3DMF cluster catalogue is now available at
cfhtlens.org.

Galaxygalaxy weak lensing is a direct probe of the mean matter distribution
around galaxies. The depth and sky coverage of the CFHT Legacy Survey yield
statistically significant galaxy halo mass measurements over a much wider range
of stellar masses ($10^{8.75}$ to $10^{11.3} M_{\odot}$) and redshifts ($0.2 <
z < 0.8$) than previous weak lensing studies. At redshift $z \sim 0.5$, the
stellartohalo mass ratio (SHMR) reaches a maximum of $4.0\pm0.2$ percent as a
function of halo mass at $\sim 10^{12.25} M_{\odot}$. We find, for the first
time from weak lensing alone, evidence for significant evolution in the SHMR:
the peak ratio falls as a function of cosmic time from $4.5 \pm 0.3$ percent at
$z \sim 0.7$ to $3.4 \pm 0.2$ percent at $z \sim 0.3$, and shifts to lower
stellar mass haloes. These evolutionary trends are dominated by red galaxies,
and are consistent with a model in which the stellar mass above which star
formation is quenched "downsizes" with cosmic time. In contrast, the SHMR of
blue, starforming galaxies is wellfit by a power law that does not evolve
with time. This suggests that blue galaxies form stars at a rate that is
balanced with their dark matter accretion in such a way that they evolve along
the SHMR locus. The redshift dependence of the SHMR can be used to constrain
the evolution of the galaxy population over cosmic time.

Higherorder, nonGaussian aspects of the largescale structure carry
valuable information on structure formation and cosmology, which is
complementary to secondorder statistics. In this work we measure second and
thirdorder weaklensing aperturemass moments from CFHTLenS and combine those
with CMB anisotropy probes. The third moment is measured with a significance of
$2\sigma$. The combined constraint on $\Sigma_8 = \sigma_8 (\Omega_{\rm
m}/0.27)^\alpha$ is improved by 10%, in comparison to the secondorder only,
and the allowed ranges for $\Omega_{\rm m}$ and $\sigma_8$ are substantially
reduced. Including general triangles of the lensing bispectrum yields tighter
constraints compared to probing mainly equilateral triangles. Second and
thirdorder CFHTLenS lensing measurements improve Planck CMB constraints on
$\Omega_{\rm m}$ and $\sigma_8$ by 26% for flat $\Lambda$CDM. For a model with
free curvature, the joint CFHTLenSPlanck result is $\Omega_{\rm m} = 0.28 \pm
0.02$ (68% confidence), which is an improvement of 43% compared to Planck
alone. We test how our results are potentially subject to three astrophysical
sources of contamination: sourcelens clustering, the intrinsic alignment of
galaxy shapes, and baryonic effects. We explore future limitations of the
cosmological use of thirdorder weak lensing, such as the nonlinear model and
the Gaussianity of the likelihood function.

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third
in a series of image analysis challenges, with a goal of testing and
facilitating the development of methods for analyzing astronomical images that
will be used to measure weak gravitational lensing. This measurement requires
extremely precise estimation of very small galaxy shape distortions, in the
presence of far larger intrinsic galaxy shapes and distortions due to the
blurring kernel caused by the atmosphere, telescope optics, and instrumental
effects. The GREAT3 challenge is posed to the astronomy, machine learning, and
statistics communities, and includes tests of three specific effects that are
of immediate relevance to upcoming weak lensing surveys, two of which have
never been tested in a community challenge before. These effects include
realistically complex galaxy models based on highresolution imaging from
space; spatially varying, physicallymotivated blurring kernel; and combination
of multiple different exposures. To facilitate entry by people new to the
field, and for use as a diagnostic tool, the simulation software for the
challenge is publicly available, though the exact parameters used for the
challenge are blinded. Sample scripts to analyze the challenge data using
existing methods will also be provided. See http://great3challenge.info and
http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.

We present a study of the relation between dark matter halo mass and the
baryonic content of host galaxies, quantified via luminosity and stellar mass.
Our investigation uses 154 deg2 of CanadaFranceHawaii Telescope Lensing
Survey (CFHTLenS) lensing and photometric data, obtained from the CFHT Legacy
Survey. We employ a galaxygalaxy lensing halo model which allows us to
constrain the halo mass and the satellite fraction. Our analysis is limited to
lenses at redshifts between 0.2 and 0.4. We express the relationship between
halo mass and baryonic observable as a power law. For the luminosityhalo mass
relation we find a slope of 1.32+/0.06 and a normalisation of
1.19+0.060.07x10^13 h70^1 Msun for red galaxies, while for blue galaxies the
bestfit slope is 1.09+0.200.13 and the normalisation is 0.18+0.040.05x10^13
h70^1 Msun. Similarly, we find a bestfit slope of 1.36+0.060.07 and a
normalisation of 1.43+0.110.08x10^13 h70^1 Msun for the stellar masshalo
mass relation of red galaxies, while for blue galaxies the corresponding values
are 0.98+0.080.07 and 0.84+0.200.16x10^13 h70^1 Msun. For red lenses, the
fraction which are satellites tends to decrease with luminosity and stellar
mass, with the sample being nearly all satellites for a stellar mass of 2x10^9
h70^2 Msun. The satellite fractions are generally close to zero for blue
lenses, irrespective of luminosity or stellar mass. This, together with the
shallower relation between halo mass and baryonic tracer, is a direct
confirmation from galaxygalaxy lensing that blue galaxies reside in less
clustered environments than red galaxies. We also find that the halo model,
while matching the lensing signal around red lenses well, is prone to
overpredicting the largescale signal for faint and less massive blue lenses.
This could be a further indication that these galaxies tend to be more isolated
than assumed. [abridged]

Extremely accurate shape measurements of galaxy images are needed to probe
dark energy properties with weak gravitational lensing surveys. To increase
survey area with a fixed observing time and pixel count, images from surveys
such as the Wide Field Infrared Survey Telescope (WFIRST) or Euclid will
necessarily be undersampled and therefore distorted by aliasing. Oversampled,
unaliased images can be obtained by combining multiple, dithered exposures of
the same source with a suitable reconstruction algorithm. Any such
reconstruction must minimally distort the reconstructed images for weak lensing
analyses to be unbiased. In this paper, we use the IMage COMbination (IMCOM)
algorithm of Rowe, Hirata, and Rhodes to investigate the effect of image
combination on shape measurements (size and ellipticity). We simulate dithered
images of sources with varying amounts of ellipticity and undersampling,
reconstruct oversampled output images from them using IMCOM, and measure shape
distortions in the output. Our simulations show that IMCOM creates no
significant distortions when the relative offsets between dithered images are
precisely known. Distortions increase with the uncertainty in those offsets but
become problematic only with relatively poor astrometric precision. E.g. for
images similar to those from the Astrophysics Focused Telescope Asset (AFTA)
implementation of WFIRST, combining eight undersampled images (sampling ratio
Q=1) with highly pessimistic uncertainty in astrometric registration
(\sigma_d~10^{3} pixels) yields an RMS shear error of O(10^{4}). Our analysis
pipeline is adapted from that of the Precision Projector Laboratory  a joint
project between NASA Jet Propulsion Laboratory and Caltech which characterizes
image sensors using laboratory emulations of astronomical data.

Weak gravitational lensing observations are a key science driver for the NASA
Wide Field Infrared Survey Telescope (WFIRST). To validate the performance of
the WFIRST infrared detectors, we have performed a laboratory emulation of weak
gravitational lensing measurements. Our experiments used a custom precision
projector system to image a target mask composed of a grid of pinholes,
emulating stellar point sources, onto a 1.7 micron cutoff Teledyne HgCdTe/H2RG
detector. We used a 880nm LED illumination source and f/22 pupil stop to
produce undersampled point spread functions similar to those expected from
WFIRST. We also emulated the WFIRST image reconstruction strategy, using the
IMage COMbination (IMCOM) algorithm to derive oversampled images from dithered,
undersampled input images. We created shear maps for this data and computed
shear correlation functions to mimic a real weak lensing analysis. After
removing only 2nd order polynomial fits to the shear maps, we found that the
correlation functions could be reduced to O(10^6). This places a conservative
upper limit on the detectorinduced bias to the correlation function (under our
test conditions). This bias is two orders of magnitude lower than the expected
weak lensing signal. Restricted to scales relevant to dark energy analyses (sky
separations > 0.5 arcmin), the bias is O(10^7): comparable to the requirement
for future weak lensing missions to avoid biasing cosmological parameter
estimates. Our experiment will need to be upgraded and repeated under different
configurations to fully characterize the shape measurement performance of
WFIRST IR detectors.

Cosmic shear is a powerful probe of cosmological parameters, but its
potential can be fully utilised only if galaxy shapes are measured with great
accuracy. Two major effects have been identified which are likely to account
for most of the bias for maximum likelihood methods in recent shear measurement
challenges. Model bias occurs when the true galaxy shape is not well
represented by the fitted model. Noise bias occurs due to the nonlinear
relationship between image pixels and galaxy shape. In this paper we
investigate the potential interplay between these two effects when an imperfect
model is used in the presence of high noise. We present analytical expressions
for this bias, which depends on the residual difference between the model and
real data. They can lead to biases not accounted for in previous calibration
schemes. By measuring the model bias, noise bias and their interaction, we
provide a complete statistical framework for measuring galaxy shapes with model
fitting methods from GRavitational lEnsing Accuracy Testing (GREAT) like
images. We demonstrate the noise and model interaction bias using a simple toy
model, which indicates that this effect can potentially be significant. Using
real galaxy images from the Cosmological Evolution Survey (COSMOS) we quantify
the strength of the model bias, noise bias and their interaction. We find that
the interaction term is often a similar size to the model bias term, and is
smaller than the requirements of the current and shortly upcoming galaxy
surveys.

We present a simulation analysis of weak gravitational lensing flexion and
shear measurement using shapelet decomposition, and identify differences
between flexion and shear measurement noise in deep survey data. Taking models
of galaxies from the Hubble Space Telescope Ultra Deep Field (HUDF) and
applying a correction for the HUDF point spread function we generate lensed
simulations of deep, optical imaging data from Hubble's Advanced Camera for
Surveys (ACS), with realistic galaxy morphologies. We find that flexion and
shear estimates differ in our measurement pipeline: whereas intrinsic galaxy
shape is typically the dominant contribution to noise in shear estimates, pixel
noise due to finite photon counts and detector read noise is a major
contributor to uncertainty in flexion estimates, across a broad range of galaxy
signaltonoise. This pixel noise also increases more rapidly as galaxy
signaltonoise decreases than is found for shear estimates. We provide simple
power law fitting functions for this behaviour, for both flexion and shear,
allowing the effect to be properly accounted for in future forecasts for
flexion measurement. Using the simulations we also quantify the systematic
biases of our shapelet flexion and shear measurement pipeline for deep Hubble
data sets such as Galaxy Evolution from Morphology and SEDs, Space Telescope
A901/902 Galaxy Evolution Survey or the Cosmic Evolution Survey. Flexion
measurement biases are found to be significant but consistent with previous
studies.

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.

We present a finelybinned tomographic weak lensing analysis of the
CanadaFranceHawaii Telescope Lensing Survey, CFHTLenS, mitigating
contamination to the signal from the presence of intrinsic galaxy alignments
via the simultaneous fit of a cosmological model and an intrinsic alignment
model. CFHTLenS spans 154 square degrees in five optical bands, with accurate
shear and photometric redshifts for a galaxy sample with a median redshift of
zm =0.70. We estimate the 21 sets of cosmic shear correlation functions
associated with six redshift bins, each spanning the angular range of
1.5<theta<35 arcmin. We combine this CFHTLenS data with auxiliary cosmological
probes: the cosmic microwave background with data from WMAP7, baryon acoustic
oscillations with data from BOSS, and a prior on the Hubble constant from the
HST distance ladder. This leads to constraints on the normalisation of the
matter power spectrum sigma_8 = 0.799 +/ 0.015 and the matter density
parameter Omega_m = 0.271 +/ 0.010 for a flat Lambda CDM cosmology. For a flat
wCDM cosmology we constrain the dark energy equation of state parameter w =
1.02 +/ 0.09. We also provide constraints for curved Lambda CDM and wCDM
cosmologies. We find the intrinsic alignment contamination to be galaxytype
dependent with a significant intrinsic alignment signal found for earlytype
galaxies, in contrast to the latetype galaxy sample for which the intrinsic
alignment signal is found to be consistent with zero.

We present and describe im3shape, a new publicly available galaxy shape
measurement code for weak gravitational lensing shear. im3shape performs a
maximum likelihood fit of a bulgeplusdisc galaxy model to noisy images,
incorporating an applied point spread function. We detail challenges faced and
choices made in its design and implementation, and then discuss various
limitations that affect this and other maximum likelihood methods. We assess
the bias arising from fitting an incorrect galaxy model using simple noisefree
images and find that it should not be a concern for current cosmic shear
surveys.
We test im3shape on the GREAT08 Challenge image simulations, and meet the
requirements for upcoming cosmic shear surveys in the case that the simulations
are encompassed by the fitted model, using a simple correction for image noise
bias. For the fiducial branch of GREAT08 we obtain a negligible additive shear
bias and subtwo percent level multiplicative bias, which is suitable for
analysis of current surveys. We fall short of the subpercent level requirement
for upcoming surveys, which we attribute to a combination of noise bias and the
mismatch between our galaxy model and the model used in the GREAT08
simulations. We meet the requirements for current surveys across all branches
of GREAT08, except those with small or high noise galaxies, which we would cut
from our analysis. Using the GREAT08 metric we we obtain a score of Q=717 for
the usable branches, relative to the goal of Q=1000 for future experiments. The
code is freely available from https://bitbucket.org/joezuntz/im3shape

We present cosmological constraints from 2D weak gravitational lensing by the
largescale structure in the CanadaFrance Hawaii Telescope Lensing Survey
(CFHTLenS) which spans 154 square degrees in five optical bands. Using accurate
photometric redshifts and measured shapes for 4.2 million galaxies between
redshifts of 0.2 and 1.3, we compute the 2D cosmic shear correlation function
over angular scales ranging between 0.8 and 350 arcmin. Using nonlinear models
of the darkmatter power spectrum, we constrain cosmological parameters by
exploring the parameter space with Population Monte Carlo sampling. The best
constraints from lensing alone are obtained for the smallscale
densityfluctuations amplitude sigma_8 scaled with the total matter density
Omega_m. For a flat LambdaCDM model we obtain sigma_8(Omega_m/0.27)^0.6 =
0.79+0.03.
We combine the CFHTLenS data with WMAP7, BOSS and an HST distanceladder
prior on the Hubble constant to get joint constraints. For a flat LambdaCDM
model, we find Omega_m = 0.283+0.010 and sigma_8 = 0.813+0.014. In the case
of a curved wCDM universe, we obtain Omega_m = 0.27+0.03, sigma_8 =
0.83+0.04, w_0 = 1.10+0.15 and Omega_K = 0.006+0.0060.004.
We calculate the Bayesian evidence to compare flat and curved LambdaCDM and
darkenergy CDM models. From the combination of all four probes, we find models
with curvature to be at moderately disfavoured with respect to the flat case. A
simple darkenergy model is indistinguishable from LambdaCDM. Our results
therefore do not necessitate any deviations from the standard cosmological
model.

The CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS) comprises deep
multicolour (u*g'r'i'z') photometry spanning 154 square degrees, with accurate
photometric redshifts and shape measurements. We demonstrate that the redshift
probability distribution function summed over galaxies provides an accurate
representation of the galaxy redshift distribution accounting for random and
catastrophic errors for galaxies with best fitting photometric redshifts z_p <
1.3.
We present cosmological constraints using tomographic weak gravitational
lensing by largescale structure. We use two broad redshift bins 0.5 < z_p <=
0.85 and 0.85 < z_p <= 1.3 free of intrinsic alignment contamination, and
measure the shear correlation function on angular scales in the range ~140
arcmin. We show that the problematic redshift scaling of the shear signal,
found in previous CFHTLS data analyses, does not afflict the CFHTLenS data. For
a flat LambdaCDM model and a fixed matter density Omega_m=0.27, we find the
normalisation of the matter power spectrum sigma_8=0.771 \pm 0.041. When
combined with cosmic microwave background data (WMAP7), baryon acoustic
oscillation data (BOSS), and a prior on the Hubble constant from the HST
distance ladder, we find that CFHTLenS improves the precision of the fully
marginalised parameter estimates by an average factor of 1.52. Combining our
results with the above cosmological probes, we find Omega_m=0.2762 \pm 0.0074
and sigma_8=0.802 \pm 0.013.

We present the CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS) that
accurately determines a weak gravitational lensing signal from the full 154
square degrees of deep multicolour data obtained by the CFHT Legacy Survey.
Weak gravitational lensing by largescale structure is widely recognised as one
of the most powerful but technically challenging probes of cosmology. We
outline the CFHTLenS analysis pipeline, describing how and why every step of
the chain from the raw pixel data to the lensing shear and photometric redshift
measurement has been revised and improved compared to previous analyses of a
subset of the same data. We present a novel method to identify data which
contributes a nonnegligible contamination to our sample and quantify the
required level of calibration for the survey. Through a series of
cosmologyinsensitive tests we demonstrate the robustness of the resulting
cosmic shear signal, presenting a scienceready shear and photometric redshift
catalogue for future exploitation.

Weak gravitational lensing has the potential to constrain cosmological
parameters to high precision. However, as shown by the Shear TEsting Programmes
(STEP) and GRavitational lEnsing Accuracy Testing (GREAT) Challenges, measuring
galaxy shears is a nontrivial task: various methods introduce different
systematic biases which have to be accounted for. We investigate how pixel
noise on the image affects the bias on shear estimates from a
MaximumLikelihood forward modelfitting approach using a sum of coelliptical
S\'{e}rsic profiles, in complement to the theoretical approach of an an
associated paper. We evaluate the bias using a simple but realistic galaxy
model and find that the effects of noise alone can cause biases of order 110%
on measured shears, which is significant for current and future lensing
surveys. We evaluate a simulationbased calibration method to create a bias
model as a function of galaxy properties and observing conditions. This model
is then used to correct the simulated measurements. We demonstrate that this
method can effectively reduce noise bias so that shear measurement reaches the
level of accuracy required for estimating cosmic shear in upcoming lensing
surveys.

Weak lensing experiments are a powerful probe of cosmology through their
measurement of the mass distribution of the universe. A challenge for this
technique is to control systematic errors that occur when measuring the shapes
of distant galaxies. In this paper we investigate noise bias, a systematic
error that arises from second order noise terms in the shape measurement
process. We first derive analytical expressions for the bias of general Maximum
Likelihood Estimators (MLEs) in the presence of additive noise. We then find
analytical expressions for a simplified toy model in which galaxies are modeled
and fitted with a Gaussian with its size as a single free parameter. Even for
this very simple case we find a significant effect. We also extend our analysis
to a more realistic 6parameter elliptical Gaussian model. We find that the
noise bias is generically of the order of the inversesquared signaltonoise
ratio (SNR) of the galaxies and is thus of the order of a percent for galaxies
of SNR of 10, i.e. comparable to the weak lensing shear signal. This is nearly
two orders of magnitude greater than the systematics requirements for future
allsky weak lensing surveys. We discuss possible ways to circumvent this
effect, including a calibration method using simulations discussed in an
associated paper.

High precision cosmology with weak gravitational lensing requires a precise
measure of the Point Spread Function across the imaging data where the accuracy
to which high spatial frequency variation can be modelled is limited by the
stellar number density across the field. We analyse dense stellar fields imaged
at the CanadaFranceHawaii Telescope to quantify the degree of high spatial
frequency variation in groundbased imaging Point Spread Functions and compare
our results to models of atmospheric turbulence. The data shows an anisotropic
turbulence pattern with an orientation independent of the wind direction and
wind speed. We find the amplitude of the high spatial frequencies to decrease
with increasing exposure time as $t^{1/2}$, and find a negligibly small
atmospheric contribution to the Point Spread Function ellipticity variation for
exposure times $t>180$ seconds. For future surveys analysing shorter exposure
data, this anisotropic turbulence will need to be taken into account as the
amplitude of the correlated atmospheric distortions becomes comparable to a
cosmological lensing signal on scales less than $\sim 10$ arcminutes. This
effect could be mitigated, however, by correlating galaxy shear measured on
exposures imaged with a time separation greater than 50 seconds, for which we
find the spatial turbulence patterns to be uncorrelated.

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image
analysis challenge aimed at the development of algorithms to analyze
astronomical images. Specifically, the challenge is to measure varying image
distortions in the presence of a variable convolution kernel, pixelization and
noise. This is the second in a series of challenges set to the astronomy,
computer science and statistics communities, providing a structured environment
in which methods can be improved and tested in preparation for planned
astronomical surveys. GREAT10 extends upon previous work by introducing
variable fields into the challenge. The "Galaxy Challenge" involves the precise
measurement of galaxy shape distortions, quantified locally by two parameters
called shear, in the presence of a known convolution kernel. Crucially, the
convolution kernel and the simulated gravitational lensing shape distortion
both now vary as a function of position within the images, as is the case for
real data. In addition, we introduce the "Star Challenge" that concerns the
reconstruction of a variable convolution kernel, similar to that in a typical
astronomical observation. This document details the GREAT10 Challenge for
potential participants. Continually updated information is also available from
http://www.greatchallenges.info.

A simple, yet general, formalism for the optimized linear combination of
astrophysical images is constructed and demonstrated. The formalism allows the
user to combine multiple undersampled images to provide oversampled output at
high precision. The proposed method is general and may be used for any
configuration of input pixels and point spread function; it also provides the
noise covariance in the output image along with a powerful metric for
describing undesired distortion to the image convolution kernel. The method
explicitly provides knowledge and control of the inevitable compromise between
noise and fidelity in the output image. We present a first prototype
implementation of the method, outlining steps taken to generate an efficient
algorithm. This implementation is then put to practical use in reconstructing
fullysampled output images using simulated, undersampled input exposures that
are designed to mimic the proposed \emph{Wide Field InfraRed Survey Telescope}
(\emph{WFIRST}). We examine results using randomly rotated and dithered input
images, while also assessing betterknown "ideal" dither patterns: comparing
results we illustrate the use of the method as a survey design tool. Finally,
we use the method to test the robustness of linear image combination when
subject to practical realities such as missing input pixels and focal plane
plate scale variations.