
Cosmic acceleration is the most surprising cosmological discovery in many
decades. Testing and distinguishing among possible explanations requires
cosmological measurements of extremely high precision probing the full history
of cosmic expansion and structure growth and, ideally, compare and contrast
matter and relativistic tracers of the gravity potential. This program is one
of the defining objectives of the WideField Infrared Survey Telescope
(WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010.
The WFIRST mission has the ability to improve these measurements by 12 orders
of magnitude compared to the current state of the art, while simultaneously
extending their redshift grasp, greatly improving control of systematic
effects, and taking a unified approach to multiple probes that provide
complementary physical information and crosschecks of cosmological results. We
describe in this annual report the activities of the Science Investigation Team
(SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017.
This team was selected by NASA in December 2015 in order to address the
stringent challenges of the WFIRST dark energy (DE) program through the
Project's formulation phase. This SIT has elected to jointly address Galaxy
Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the
fact that the imaging and spectroscopic elements of the HLS will be realized as
an integrated observing program, and they jointly impose requirements on
performance and operations. WFIRST is designed to be able to deliver a
definitive result on the origin of cosmic acceleration. It is not optimized for
Figure of Merit sensitivity but for control of systematic uncertainties and for
having multiple techniques each with multiple crosschecks. Our SIT work
focuses on understanding the potential systematics in the WFIRST DE
measurements.

We present cosmological constraints from the Dark Energy Survey (DES) using a
combined analysis of angular clustering of red galaxies and their
crosscorrelation with weak gravitational lensing of background galaxies. We
use a 139 square degree contiguous patch of DES data from the Science
Verification (SV) period of observations. Using large scale measurements, we
constrain the matter density of the Universe as Omega_m = 0.31 +/ 0.09 and the
clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/ 0.13
after marginalizing over seven nuisance parameters and three additional
cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16}
= 0.74 +/ 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 =
0.78 +/ 0.09 using two bins over the range 0.2 <z< 0.5. We study the
robustness of the results under changes in the data vectors, modelling and
systematics treatment, including photometric redshift and shear calibration
uncertainties, and find consistency in the derived cosmological parameters. We
show that our results are consistent with previous cosmological analyses from
DES and other data sets and conclude with a joint analysis of DES angular
clustering and galaxygalaxy lensing with Planck CMB data, Baryon Accoustic
Oscillations and Supernova type Ia measurements.

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.

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 .

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 explore the utility of Karhunen Loeve (KL) analysis in solving practical
problems in the analysis of gravitational shear surveys. Shear catalogs from
largefield weak lensing surveys will be subject to many systematic
limitations, notably incomplete coverage and pixellevel masking due to
foreground sources. We develop a method to use two dimensional KL eigenmodes of
shear to interpolate noisy shear measurements across masked regions. We explore
the results of this method with simulated shear catalogs, using statistics of
highconvergence regions in the resulting map. We find that the KL procedure
not only minimizes the bias due to masked regions in the field, it also reduces
spurious peak counts from shape noise by a factor of ~ 3 in the cosmologically
sensitive regime. This indicates that KL reconstructions of masked shear are
not only useful for creating robust convergence maps from masked shear
catalogs, but also offer promise of improved parameter constraints within
studies of shear peak statistics.

We present a new method for constructing threedimensional mass maps from
gravitational lensing shear data. We solve the lensing inversion problem using
truncation of singular values (within the context of generalized least squares
estimation) without a priori assumptions about the statistical nature of the
signal. This singular value framework allows a quantitative comparison between
different filtering methods: we evaluate our method beside the previously
explored Wiener filter approaches. Our method yields nearoptimal angular
resolution of the lensing reconstruction and allows cluster sized halos to be
deblended robustly. It allows for mass reconstructions which are 23
ordersofmagnitude faster than the Wiener filter approach; in particular, we
estimate that an allsky reconstruction with arcminute resolution could be
performed on a timescale of hours. We find however that linear, nonparametric
reconstructions have a fundamental limitation in the resolution achieved in the
redshift direction.

We present the results of the GREAT08 Challenge, a blind analysis challenge
to infer weak gravitational lensing shear distortions from images. The primary
goal was to stimulate new ideas by presenting the problem to researchers
outside the shear measurement community. Six GREAT08 Team methods were
presented at the launch of the Challenge and five additional groups submitted
results during the 6 month competition. Participants analyzed 30 million
simulated galaxies with a range in signal to noise ratio, pointspread function
ellipticity, galaxy size, and galaxy type. The large quantity of simulations
allowed shear measurement methods to be assessed at a level of accuracy
suitable for currently planned future cosmic shear observations for the first
time. Different methods perform well in different parts of simulation parameter
space and come close to the target level of accuracy in several of these. A
number of fresh ideas have emerged as a result of the Challenge including a
reexamination of the process of combining information from different galaxies,
which reduces the dependence on realistic galaxy modelling. The image
simulations will become increasingly sophisticated in future GREAT challenges,
meanwhile the GREAT08 simulations remain as a benchmark for additional
developments in shear measurement algorithms.

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.

In weak lensing investigations, galaxy shapes are deconvolved for the effects
of the point spread function (PSF) using stellar images. In this paper we use
physical models of the telescope optics to understand the spatial variation of
the PSF in the image plane. We introduce a set of parameters to model the key
aberrations, which include defocus, focal plane tilt, primary and offaxis
astigmatism. We also include the effects of guiding and seeing. We test our
model with data from the Blanco 4 meter telescope in Cerro Tololo, Chile. We
find that the physical model describes a substantial part of the PSF size and
anisotropy over the field of view (over 90 percent of it, based on a
chisquared metric). We identify the primary contributors to the PSF patterns
and study their covariances and principal components. We also identify
correlations with the effect of gravity on the telescope. Finally, we discuss
the improvements in PSF estimation that may be achieved by combining the
physical model in this study with the purely empirical approach of Jarvis and
Jain (2004).

We consider the possible gain in the measurement of lensing shear from
imaging data in multiple filters. Galaxy shapes may differ significantly across
filters, so that the same galaxy offers multiple samples of the shear. On the
other extreme, if galaxy shapes are identical in different filters, one can
combine them to improve the signaltonoise and thus increase the effective
number density of faint, high redshift galaxies. We use the GOODS dataset to
test these scenarios by calculating the covariance matrix of galaxy
ellipticities in four visual filters (B,V,i,z). We find that galaxy shapes are
highly correlated, and estimate the gain in galaxy number density by combining
their shapes.

The Shear TEsting Programme (STEP) is a collaborative project to improve the
accuracy and reliability of weak lensing measurement, in preparation for the
next generation of widefield surveys. We review sixteen current and emerging
shear measurement methods in a common language, and assess their performance by
running them (blindly) on simulated images that contain a known shear signal.
We determine the common features of algorithms that most successfully recover
the input parameters. We achieve previously unattained discriminatory precision
in our analysis, via a combination of more extensive simulations, and pairs of
galaxy images that have been rotated with respect to each other, thus removing
noise from their intrinsic ellipticities. The robustness of our simulation
approach is also confirmed by testing the relative calibration of methods on
real data.
Weak lensing measurement has improved since the first STEP paper. Several
methods now consistently achieve better than 2% precision, and are still being
developed. However, the simulations can now distinguish all methods from
perfect performance. Our main concern continues to be the potential for a
multiplicative shear calibration bias: not least because this can not be
internally calibrated with real data. We determine which galaxy populations are
responsible and, by adjusting the simulated observing conditions, we also
investigate the effects of instrumental and atmospheric parameters. We have
isolated several previously unrecognised aspects of galaxy shape measurement,
in which focussed development could provide further progress towards the
subpercent level of precision desired for future surveys.
[ABRIDGED]

We perform a cosmological parameter analysis of the 75 square degree CTIO
lensing survey in conjunction with CMB and Type Ia supernovae data. For Lambda
CDM cosmologies, we find that the amplitude of the power spectrum at low
redshift is given by sigma_8 = 0.81 (+0.15,0.10, 95% c.l.), where the error
bar includes both statistical and systematic errors. The total of all
systematic errors is smaller than the statistical errors, but they do make up a
significant fraction of the error budget. We find that weak lensing improves
the constraints on dark energy as well. The (constant) dark energy equation of
state parameter, w, is measured to be 0.89 (+0.16,0.21, 95% c.l.).
Marginalizing over a constant $w$ slightly changes the estimate of sigma_8 to
0.79 (+0.17, 0.14, 95% c.l.). We also investigate variable w cosmologies, but
find that the constraints weaken considerably; the next generation surveys are
needed to obtain meaningful constraints on the possible time evolution of dark
energy.

The Shear TEsting Programme, STEP, is a collaborative project to improve the
accuracy and reliability of all weak lensing measurements in preparation for
the next generation of widefield surveys. In this first STEP paper we present
the results of a blind analysis of simulated groundbased observations of
relatively simple galaxy morphologies. The most successful methods are shown to
achieve percent level accuracy. From the cosmic shear pipelines that have been
used to constrain cosmology, we find weak lensing shear measured to an accuracy
that is within the statistical errors of current weak lensing analyses, with
shear measurements accurate to better than 7%. The dominant source of
measurement error is shown to arise from calibration uncertainties where the
measured shear is over or underestimated by a constant multiplicative factor.
This is of concern as calibration errors cannot be detected through standard
diagnostic tests. The measured calibration errors appear to result from stellar
contamination, false object detection, the shear measurement method itself,
selection bias and/or the use of biased weights. Additive systematics (false
detections of shear) resulting from residual pointspread function anisotropy
are, in most cases, reduced to below an equivalent shear of 0.001, an order of
magnitude below cosmic shear distortions on the scales probed by current
surveys.
Our results provide a snapshot view of the accuracy of current groundbased
weak lensing methods and a benchmark upon which we can improve. To this end we
provide descriptions of each method tested and include details of the eight
different implementations of the commonly used Kaiser, Squires and Broadhurst
(1995) method (KSB+) to aid the improvement of future KSB+ analyses.

Given the basic parameters of a cosmic shear weak lensing survey, how well
can systematic errors due to anisotropy in the point spread function (PSF) be
corrected? The largest source of error in this correction to date has been the
interpolation of the PSF to the locations of the galaxies. To address this
error, we separate the PSF patterns into components that recur in multiple
exposures/pointings and those that vary randomly between different exposures
(such as those due to the atmosphere). In an earlier study we developed a
principal component approach to correct the recurring PSF patterns (Jarvis and
Jain 2004). In this paper we show how randomly varying PSF patterns can also be
circumvented in the measurement of shear correlations. For the twopoint
correlation function this is done by simply using pairs of galaxy shapes
measured in different exposures. Combining the two techniques allows us to
tackle generic combinations of PSF anisotropy patterns. The second goal of this
paper is to give a formalism for quantifying residual systematic errors due to
PSF patterns. We show how the main PSF corrections improve with increasing
survey area (and thus can stay below the reduced statistical errors), and we
identify the residual errors which do not scale with survey area. Our formalism
can be applied both to planned lensing surveys to optimize instrumental and
survey parameters and to actual lensing data to quantify residual errors.

We introduce a new algorithm for interpolating measurements of the
pointspread function (PSF) using stars from many exposures. The principal
components of the variation in the PSF pattern from multiple exposures are used
to solve for better fits for each individual exposure. These improved fits are
then used to correct the weak lensing shapes. Since we expect some degree of
correlation in the PSF anisotropy across exposures, using information from
stars in all exposures leads to a significant gain in the accuracy of the
estimated PSF. It means that in general, the accuracy of PSF reconstruction is
limited not by the number density of stars per exposure, but by the stacked
number density across all exposures in a given survey. This technique is
applied to the 75 square degree CTIO lensing survey, and we find that the PSF
variation is well described by a small number of components. There is a
significant improvement in the lensing measurements: the residual stellar PSF
correlations are reduced by several orders of magnitude, and the measured
Bmode in the twopoint correlation is consistent with zero down to 1
arcminute. We also discuss the applications of the PCA technique to future
surveys.