Cosmic shear - the weak gravitational lensing effect generated by fluctuations of the gravitational tidal fields of the large-scale structure - is one of the most promising tools for current and future cosmological analyses. The spherical-Bessel decomposition of the cosmic shear field ("3D cosmic shear") is one way to maximise the amount of redshift information in a lensing analysis and therefore provides a powerful tool to investigate in particular the growth of cosmic structure that is crucial for dark energy studies. However, the computation of simulated 3D cosmic shear covariance matrices presents numerical difficulties, due to the required integrations over highly oscillatory functions. We present and compare two numerical methods and relative implementations to perform these integrations. We then show how to generate 3D Gaussian random fields on the sky in spherical coordinates, starting from the 3D cosmic shear covariances. To validate our field-generation procedure, we calculate the Minkowski functionals associated with our random fields, compare them with the known expectation values for the Gaussian case and demonstrate parameter inference from Minkowski functionals from a cosmic shear survey. This is a first step towards producing fully 3D Minkowski functionals for a lognormal field in 3D to extract Gaussian and non-Gaussian information from the cosmic shear field, as well as towards the use of Minkowski functionals as a probe of cosmology beyond the commonly used two-point statistics.

We present the first calculation of the cross-correlation between three-dimensional cosmic shear and the integrated Sachs-Wolfe (iSW) effect. Both signals are combined in a single formalism, which permits the computation of the full covariance matrix. In order to avoid the uncertainties presented by the non-linear evolution of the matter power spectrum and intrinsic alignments of galaxies, our analysis is restricted to large scales, i.e. multipoles below l=1000. We demonstrate in a Fisher analysis that this reduction compared to other studies of three-dimensional weak lensing extending to smaller scales is compensated by the information that is gained if the additional iSW signal and in particular its cross-correlation with lensing data are considered. Given the observational standards of upcoming weak lensing surveys like Euclid, marginal errors on cosmological parameters decrease by ten per cent compared to a cosmic shear experiment if both types of information are combined without a CMB prior. Once the constraining power of CMB data is added, the improvement becomes marginal.

We investigate how observations of strong lensing can be used to infer cosmological parameters, in particular the equation of state of dark energy. We focus on the growth of the critical lines of lensing clusters with the source redshift as this behaviour depends on the distance-redshift relation and is therefore cosmologically sensitive. Purely analytical approaches are generally insufficient because they rely on axisymmetric mass distributions and thus cannot take irregular critical curves into account. We devise a numerical method based on the Metropolis-Hastings algorithm: an elliptical generalization of the NFW density profile is used to fit a lens model to an observed configuration of giant luminous arcs while simultaneously optimizing the geometry. A semi-analytic method, which derives geometric parameters from critical points, is discussed as a faster alternative. We test the approaches on mock observations of gravitational lensing by a numerically simulated cluster. We find that no constraints can be derived from observations of individual clusters if no knowledge of the underlying mass distribution is assumed. Uncertainties are improved if a fixed lens model is used for a purely geometrical optimization, but the choice of a parametric model may produce strong biases.