The identification of strong gravitational lenses in large surveys has historically been a rather time consuming exercise. Early data from the Herschel Astrophysical Terahertz Large Area Survey (Herschel-ATLAS) demonstrate that lenses can be identified efficiently at submillimetre wavelengths using a simple flux criteria. Motivated by that development, this work considers the statistical properties of strong gravitational lens systems which have been, and will be, found by the Herschel-ATLAS. Analytical models of lens statistics are tested with the current best estimates for the various model ingredients. These include the cosmological parameters, the mass function and the lens density profile, for which we consider the singular isothermal sphere (SIS) and the Navarro, Frenk & White (NFW) approximations. The five lenses identified in the Herschel-ATLAS Science Demonstration Phase suggest a SIS density profile is preferred, but cannot yet constrain \Omega_\Lambda to an accuracy comparable with other methods. The complete Herschel-ATLAS data set should be sufficient to yield competitive constraints on \Omega_\Lambda. Whilst this huge number of lenses has great potential for constraining cosmological parameters, they will be most powerful in constraining uncertainty in astrophysical processes. Further investigation is needed to fully exploit this unprecedented data set.

We consider the issue of characterizing the coherent large-scale patterns from CMB temperature maps in globally anisotropic cosmologies. The methods we investigate are reasonably general; the particular models we test them on are the homogeneous but anisotropic relativistic cosmologies described by the Bianchi classification. Although the temperature variations produced in these models are not stochastic, they give rise to a "non-Gaussian" distribution of temperature fluctuations over the sky that is a partial diagnostic of the model. We explore two methods for quantifying non-Gaussian and/or non-stationary fluctuation fields in order to see how they respond to the Bianchi models.We first investigate the behavior of phase correlations between the spherical harmonic modes of the maps. Then we examine the behavior of the multipole vectors of the temperature distribution which, though defined in harmonic space, can indicate the presence of a preferred direction in real space, i.e. on the 2-sphere. These methods give extremely clear signals of the presence of anisotropy when applied to the models we discuss, suggesting that they have some promise as diagnostics of the presence of global asymmetry in the Universe.

We discuss and compare two alternative models for the two-point angular correlation function of galaxies detected through the sub-millimetre emission using the Herschel Space Observatory. The first, now-standard Halo Model, which represents the angular correlations as arising from one-halo and two-halo contributions, is flexible but complex and rather unwieldy. The second model is based on a much simpler approach: we incorporate a fitting function method to estimate the matter correlation function with approximate model of the bias inferred from the estimated redshift distribution to find the galaxy angular correlation function. We find that both models give a good account of the shape of the correlation functions obtained from published preliminary studies of the HerMES and H-ATLAS surveys performed using Herschel, and yield consistent estimates of the minimum halo mass within which the sub-millimetre galaxies must reside. We note also that both models predict an inflection in the correlation function at intermediate angular scales, so the presence of the feature in the measured correlation function does not unambiguously indicate the presence of intra-halo correlations. The primary barrier to more detailed interpretation of these clustering measurements lies in the substantial uncertainty surrounding the redshift distribution of the sources.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

All-sky maps of the cosmic microwave background temperature fluctuations are usually represented by a spherical harmonic decomposition involving modes labelled by their degree l and order m (where -l < m < +l). The zonal modes (i.e those with m = 0) are of particular interest because they vary only with galactic latitude; any anomalous behaviour in them might therefore be an indication of erroneous foreground substraction. We perform a simple statistical analysis of the modes with low l for sky maps derived via different cleaning procedures from the Wilkinson Microwave Anisotropy Probe (WMAP) and show that the zonal modes provide a useful diagnostic of possible systematics.