Weak gravitational lensing measurements are traditionally made at optical wavelengths where many highly resolved galaxy images are readily available. However, the Square Kilometre Array (SKA) holds great promise for this type of measurement at radio
wavelengths owing to its greatly increased sensitivity and resolution over typical radio surveys. The key to successful weak lensing experiments is in measuring the shapes of detected sources to high accuracy. In this document we describe a simulation pipeline designed to simulate radio images of the quality required for weak lensing, and will be typical of SKA observations. We provide as input, images with realistic galaxy shapes which are then simulated to produce images as they would have been observed with a given radio interferometer. We exploit this pipeline to investigate various stages of a weak lensing experiment in order to better understand the effects that may impact shape measurement. We first show how the proposed SKA1-Mid array configurations perform when we compare the (known) input and output ellipticities. We then investigate how making small changes to these array configurations impact on this input-outut ellipticity comparison. We also demonstrate how alternative configurations for SKA1-Mid that are smaller in extent, and with a faster survey speeds produce similar performance to those originally proposed. We then show how a notional SKA configuration performs in the same shape measurement challenge. Finally, we describe ongoing efforts to utilise our simulation pipeline to address questions relating to how applicable current (mostly originating from optical data) shape measurement techniques are to future radio surveys. As an alternative to such image plane techniques, we lastly discuss a shape measurement technique based on the shapelets formalism that reconstructs the source shapes directly from the visibility data.
We present a study of weak lensing shear measurements for simulated galaxy images at radio wavelengths. We construct a simulation pipeline into which we can input galaxy images of known ellipticity, and with which we then simulate observations with e
MERLIN and the international LOFAR array. The simulations include the effects of the CLEAN algorithm, uv sampling, observing angle, and visibility noise, and produce realistic restored images of the galaxies. We apply a shapelet-based shear measurement method to these images and test our ability to recover the true source ellipticities. We model and deconvolve the effective PSF, and find suitable parameters for CLEAN and shapelet decomposition of galaxies. We demonstrate that ellipticities can be measured faithfully in these radio simulations, with no evidence of an additive bias and a modest (10%) multiplicative bias on the ellipticity measurements. Our simulation pipeline can be used to test shear measurement procedures and systematics for the next generation of radio telescopes.
We carry out an exploratory weak gravitational lensing analysis on a combined VLA and MERLIN radio data set: a deep (3.3 micro-Jy beam^-1 rms noise) 1.4 GHz image of the Hubble Deep Field North. We measure the shear estimator distribution at this rad
io sensitivity for the first time, finding a similar distribution to that of optical shear estimators for HST ACS data in this field. We examine the residual systematics in shear estimation for the radio data, and give cosmological constraints from radio-optical shear cross-correlation functions. We emphasize the utility of cross-correlating shear estimators from radio and optical data in order to reduce the impact of systematics. Unexpectedly we find no evidence of correlation between optical and radio intrinsic ellipticities of matched objects; this result improves the properties of optical-radio lensing cross-correlations. We explore the ellipticity distribution of the radio counterparts to optical sources statistically, confirming the lack of correlation; as a result we suggest a connected statistical approach to radio shear measurements.
