The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured.

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.

We present the exact joint likelihood of pseudo-$C_\ell$ power spectrum estimates measured from an arbitrary number of Gaussian cosmological fields. Our method is applicable to both spin-0 fields and spin-2 fields, including a mixture of the two, and is relevant to Cosmic Microwave Background, weak lensing and galaxy clustering analyses. We show that Gaussian cosmological fields are mixed by a mask in such a way that retains their Gaussianity, without making any assumptions about the mask geometry. We then show that each auto- or cross-pseudo-$C_\ell$ estimator can be written as a quadratic form, and apply the known joint distribution of quadratic forms to obtain the exact joint likelihood of a set of pseudo-$C_\ell$ estimates in the presence of an arbitrary mask. Considering the polarisation of the Cosmic Microwave Background as an example, we show using simulations that our likelihood recovers the full, exact multivariate distribution of $EE$, $BB$ and $EB$ pseudo-$C_\ell$ power spectra. Our method provides a route to robust cosmological constraints from future Cosmic Microwave Background and large-scale structure surveys in an era of ever-increasing statistical precision.

We investigate cosmological parameter inference from realistic simulated weak lensing image data using only galaxy position angles, as opposed to full-ellipticity information. We demonstrate that input shear fields can be accurately reconstructed using only the statistics of source galaxy position angles and that, from these shear fields, we can successfully recover power spectra and infer the input cosmology. This paper builds on previous work on angle-only weak lensing estimation by extending the method to deal with variable and anisotropic PSF convolution and variable shear fields. Previous work employed a weighting scheme to downweight the contribution to shear estimates from sources aligned with the PSF. This work removes the need to downweight sources by convolving them with an image of the PSF rotated by $90^{\circ}$. We show thats this successfully undoes the rotation caused by PSF convolution, assuming we have reliable images of the PSF. We find that we can accurately recover the input shear signal from a simulated Stage III-like weak lensing data set using only the position angles to within an overall scale factor, and that the scale factor can be determined using a cosmology independent simulation with noise, galaxy, and PSF properties that match the observed data set. We then demonstrate that we can constrain cosmological parameters using angle-only shear estimates with a constraining power comparable to using current state-of-the-art shape measurement techniques that provide full-ellipticity information.

We present a weak lensing analysis of the 3 GHz VLA radio survey of the COSMOS field, which we correlate with overlapping HST-ACS optical observations using both intrinsic galaxy shape and cosmic shear correlation statistics. After cross-matching sources between the two catalogues, we measure the correlations of galaxy position angles and find a Pearson correlation coefficient of $0.14 \pm 0.03$. This is a marked improvement from previous studies which found very weak, or non-existent correlations, and gives insight into the emission processes of radio and optical galaxies. We also extract power spectra of averaged galaxy ellipticities (the primary observable for cosmic shear) from the two catalogues, and produce optical-optical, radio-optical and radio-radio spectra. The optical-optical auto-power spectrum was measured to a detection significance of 9.80$\sigma$ and is consistent with previous observations of the same field. For radio spectra (which we do not calibrate, given the unknown nature of their systematics), although we do not detect significant radio-optical (1.50$\sigma$) or radio-radio (1.45$\sigma$) $E$-mode power spectra, we do find the $E$-mode spectra to be more consistent with the shear signal expected from previous studies than with a null signal, and vice versa for $B$-mode and $EB$ cross-correlation spectra. Our results give promise that future radio weak lensing surveys with larger source number densities over larger areas will have the capability to measure significant weak lensing signals.

We develop a new method of extracting simultaneous measurements of weak lensing shear and a local rotation of the plane of polarization using observations of resolved radio sources. We show that the direction of polarization is statistically linked with that of the gradient of the total intensity field, and this provides the basis of our method. Using a number of sources spread over the sky, this method allows constraints to be placed on cosmic shear and birefringence, and it can be applied to any resolved radio sources for which such a correlation exists. Assuming that the rotation and shear are constant across the source, we use this relationship to construct a quadratic estimator and investigate its properties using simulated observations. We develop a calibration scheme using simulations based on the observed images to mitigate a bias which occurs in the presence of measurement errors and an astrophysical scatter on the polarization. The method is applied directly to archival data of radio galaxies where we measure a mean rotation signal of $\left<\omega\right>=-2.03^{\circ}\pm0.75^{\circ}$ and an average shear compatible with zero using 30 reliable sources. This level of constraint on an overall rotation is comparable with current leading constraints from CMB experiments and is expected to increase by at least an order of magnitude with future high precision radio surveys, such as those performed by the SKA. We also measure the shear and rotation two-point correlation functions and estimate the number of sources required to detect shear and rotation correlations in future surveys.

Weak lensing has become an increasingly important tool in cosmology and the use of galaxy shapes to measure cosmic shear has become routine. The weak-lensing distortion tensor contains two other effects in addition to the two components of shear: the convergence and rotation. The rotation mode is not measurable using the standard cosmic shear estimators based on galaxy shapes, as there is no information on the original shapes of the images before they were lensed. Due to this, no estimator has been proposed for the rotation mode in cosmological weak-lensing surveys, and the rotation mode has never been constrained. Here, we derive an estimator for this quantity, which is based on the use of radio polarisation measurements of the intrinsic position angles of galaxies. The rotation mode can be sourced by physics beyond $\Lambda$CDM, and also offers the chance to perform consistency checks of $\Lambda$CDM and of weak-lensing surveys themselves. We present simulations of this estimator and show that, for the pedagogical example of cosmic string spectra, this estimator could detect a signal that is consistent with the constraints from Planck. We examine the connection between the rotation mode and the shear $B$-modes and thus how this estimator could help control systematics in future radio weak-lensing surveys.

We develop and apply the position angle-only shear estimator of Whittaker et al. (2014) to realistic galaxy images. This is done by demonstrating the method on the simulations of the GREAT3 challenge (Mandelbaum et al. 2014b), which include contributions from anisotropic PSFs. We measure the position angles of the galaxies using three distinct methods - the integrated light method, quadrupole moments of surface brightness, and using model-based ellipticity measurements provided by IM3SHAPE. A weighting scheme is adopted to address biases in the position angle measurements which arise in the presence of an anisotropic PSF. Biases on the shear estimates, due to measurement errors on the position angles and correlations between the measurement errors and the true position angles, are corrected for using simulated galaxy images and an iterative procedure. The properties of the simulations are estimated using the deep field images provided as part of the challenge. A method is developed to match the distributions of galaxy fluxes and half-light radii from the deep fields to the corresponding distributions in the field of interest. We recover angle-only shear estimates with a performance close to current well-established model and moments-based methods for all three angle measurement techniques. The Q-values for all three methods are found to be Q~400. The code is freely available online at http://www.jb.man.ac.uk/~mbrown/angle_only_shear/.

We discuss methods for performing weak lensing using radio observations to recover information about the intrinsic structural properties of the source galaxies. Radio surveys provide unique information that can benefit weak lensing studies, such as HI emission, which may be used to construct galaxy velocity maps, and polarized synchrotron radiation; both of which provide information about the unlensed galaxy and can be used to reduce galaxy shape noise and the contribution of intrinsic alignments. Using a proxy for the intrinsic position angle of an observed galaxy, we develop techniques for cleanly separating weak gravitational lensing signals from intrinsic alignment contamination in forthcoming radio surveys. Random errors on the intrinsic orientation estimates introduce biases into the shear and intrinsic alignment estimates. However, we show that these biases can be corrected for if the error distribution is accurately known. We demonstrate our methods using simulations, where we reconstruct the shear and intrinsic alignment auto and cross-power spectra in three overlapping redshift bins. We find that the intrinsic position angle information can be used to successfully reconstruct both the lensing and intrinsic alignment power spectra with negligible residual bias.

We develop a method for performing a weak lensing analysis using only measurements of galaxy position angles. By analysing the statistical properties of the galaxy orientations given a known intrinsic ellipticity distribution, we show that it is possible to obtain estimates of the shear by minimizing a $\chi^2$ statistic. The method is demonstrated using simulations where the components of the intrinsic ellipticity are taken to be Gaussian distributed. Uncertainties in the position angle measurements introduce a bias into the shear estimates which can be reduced to negligible levels by introducing a correction term into the formalism. We generalize our approach by developing an algorithm to obtain direct shear estimators given any azimuthally symmetric intrinsic ellipticity distribution. We introduce a method of measuring the position angles of the galaxies from noisy pixelized images, and propose a method to correct for biases which arise due to pixelization and correlations between measurement errors and galaxy ellipticities. We also develop a method to constrain the sample of galaxies used to obtain an estimate of the intrinsic ellipticity distribution such that fractional biases in the resulting shear estimates are below a given threshold value. We demonstrate the angle only method by applying it to simulations where the ellipticities are taken to follow a log-normal distribution. We compare the performance of the position angle only method with the standard method based on full ellipticity measurements by reconstructing lensing convergence maps from both numerical simulations and from the CFHTLenS data. We find that the difference between the convergence maps reconstructed using the two methods is consistent with noise.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.