
It is common practice in cosmology to model largescale structure observables
as lognormal random fields, and this approach has been successfully applied in
the past to the matter density and weak lensing convergence fields separately.
We argue that this approach has fundamental limitations which prevent its use
for jointly modelling these two fields since the lognormal distribution's shape
can prevent certain correlations to be attainable. Given the need of ongoing
and future largescale structure surveys for fast joint simulations of
clustering and weak lensing, we propose two ways of overcoming these
limitations. The first approach slightly distorts the power spectra of the
fields using one of two algorithms that minimises either the absolute or the
fractional distortions. The second one is by obtaining more accurate
convergence marginal distributions, for which we provide a fitting function, by
integrating the lognormal density along the line of sight. The latter approach
also provides a way to determine directly from theory the skewness of the
convergence distribution and, therefore, the parameters for a lognormal fit. We
present the public code Fullsky Lognormal Astrofields Simulation Kit (FLASK)
which can make tomographic realisations on the sphere of an arbitrary number of
correlated lognormal or Gaussian random fields by applying either of the two
proposed solutions, and show that it can create joint simulations of clustering
and lensing with subpercent accuracy over relevant angular scales and
redshift ranges.

We present a method of calibrating the properties of photometric redshift
bins as part of a larger Markov Chain Monte Carlo (MCMC) analysis for the
inference of cosmological parameters. The redshift bins are characterised by
their mean and variance, which are varied as free parameters and marginalised
over when obtaining the cosmological parameters. We demonstrate that the
likelihood function for crosscorrelations in an angular power spectrum
framework tightly constrains the properties of bins such that they may be well
determined, reducing their influence on cosmological parameters and avoiding
the bias from poorly estimated redshift distributions. We demonstrate that even
with only three photometric and three spectroscopic bins, we can recover
accurate estimates of the mean redshift of a bin to within $\Delta\mu \approx
34 \times10^{3}$ and the width of the bin to $\Delta\sigma \approx
1\times10^{3}$ for galaxies near $z = 1$. This indicates that we may be able
to bring down the photometric redshift errors to a level which is in line with
the requirements for the next generation of cosmological experiments.

Host galaxy identification is a crucial step for modern supernova (SN)
surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey
Telescope (LSST), which will discover SNe by the thousands. Spectroscopic
resources are limited, so in the absence of realtime SN spectra these surveys
must rely on host galaxy spectra to obtain accurate redshifts for the Hubble
diagram and to improve photometric classification of SNe. In addition, SN
luminosities are known to correlate with hostgalaxy properties. Therefore,
reliable identification of host galaxies is essential for cosmology and SN
science. We simulate SN events and their locations within their host galaxies
to develop and test methods for matching SNe to their hosts. We use both real
and simulated galaxy catalog data from the Advanced Camera for Surveys General
Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe
residing in undetected faint hosts into our analysis, with an assumed hostless
rate of 5%. Our fully automated algorithm is run on catalog data and matches
SNe to their hosts with 91% accuracy. We find that including a machine learning
component, run after the initial matching algorithm, improves the accuracy
(purity) of the matching to 97% with a 2% cost in efficiency (true positive
rate). Although the exact results are dependent on the details of the survey
and the galaxy catalogs used, the method of identifying host galaxies we
outline here can be applied to any transient survey.

We present ANNz2, a new implementation of the public software for photometric
redshift (photoz) estimation of Collister and Lahav (2004), which now includes
generation of full probability distribution functions (PDFs). ANNz2 utilizes
multiple machine learning methods, such as artificial neural networks and
boosted decision/regression trees. The objective of the algorithm is to
optimize the performance of the photoz estimation, to properly derive the
associated uncertainties, and to produce both singlevalue solutions and PDFs.
In addition, estimators are made available, which mitigate possible problems of
nonrepresentative or incomplete spectroscopic training samples. ANNz2 has
already been used as part of the first weak lensing analysis of the Dark Energy
Survey, and is included in the experiment's first public data release. Here we
illustrate the functionality of the code using data from the tenth data release
of the Sloan Digital Sky Survey and the Baryon Oscillation Spectroscopic
Survey. The code is available for download at
https://github.com/IftachSadeh/ANNZ .

The observed 21cm signal from the epoch of reionization will be distorted
along the lineofsight by the peculiar velocities of matter particles. These
redshiftspace distortions will affect the contrast in the signal and will also
make it anisotropic. This anisotropy contains information about the
crosscorrelation between the matter density field and the neutral hydrogen
field, and could thus potentially be used to extract information about the
sources of reionization. In this paper, we study a collection of simulated
reionization scenarios assuming different models for the sources of
reionization. We show that the 21cm anisotropy is best measured by the
quadrupole moment of the power spectrum. We find that, unless the properties of
the reionization sources are extreme in some way, the quadrupole moment evolves
very predictably as a function of global neutral fraction. This predictability
implies that redshiftspace distortions are not a very sensitive tool for
distinguishing between reionization sources. However, the quadrupole moment can
be used as a modelindependent probe for constraining the reionization history.
We show that such measurements can be done to some extent by firstgeneration
instruments such as LOFAR, while the SKA should be able to measure the
reionization history using the quadrupole moment of the power spectrum to great
accuracy.

We present spectroscopic confirmation of two new lensed quasars via data
obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been
selected from the Dark Energy Survey (DES) and WISE based on their multiband
photometry and extended morphology in DES images. Images of DES J01155244 show
two blue point sources at either side of a red galaxy. Our longslit data
confirm that both point sources are images of the same quasar at $z_{s}=1.64.$
The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is
in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two
blue components are visible in the DES and SDSS images. The SDSS fiber spectrum
shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II
and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground
lens galaxy. The longslit Magellan spectra show that the blue components are
resolved images of the same quasar. The Einstein Radius is $0.68$"
corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other
candidates were observed and rejected, two being lowredshift pairs of
starburst galaxies, and one being a quasar behind a blue star. These first
confirmation results provide an important empirical validation of the
datamining and modelbased selection that is being applied to the entire DES
dataset.

We explore the potential of a future, ultrahigh energy cosmic ray (UHECR)
experiment, that is able to overcome the limitation of low statistics, to
detect anisotropy in the arrival directions of UHECRs. We concentrate on the
lower energy range of future instruments (E > 50 EeV), where, if the UHECR
source number density is not too low, the sources should be numerous enough to
imprint a clustering pattern in the sky, and thus possibly in the UHECR arrival
directions. Under these limits, the anisotropy signal should be dominated by
the clustering of astrophysical sources perse in the largescale structures,
and not the clustering of events around individual sources. We study the
potential for a statistical discrimination between different astrophysical
models which we parametrise by the number density of UHECR sources, the
possible bias of the UHECR accelerators with respect to the galaxy
distribution, and the unknown fraction of UHECRs that have been deflected by
large angles. We demonstrate that an orderofmagnitude increase in statistics
would allow to discriminate between a variety of astrophysical models, provided
that a subsample of light elements can be extracted, and that it represents at
least ~70% of the overall flux, sensitive to the UHECR source number density.
Even without knowledge of the composition, an anisotropy at the 99.7% level
should be detectable when the number of detected events exceeds 2000 beyond 50
EeV, as long as the composition is proton dominated, and the number density of
UHECR sources is relatively high. If the UHECR sources are strongly biased
relative to the galaxy distribution, as are for example galaxy clusters, an
anisotropy at the 99.7% should be detectable once the number of detected events
exceeds 1000, if the fraction of protons at the highest energies is ~60% or
higher.

The exceptional sensitivity of the SKA will allow observations of the Cosmic
Dawn and Epoch of Reionization (CD/EoR) in unprecedented detail, both
spectrally and spatially. This wealth of information is buried under Galactic
and extragalactic foregrounds, which must be removed accurately and precisely
in order to reveal the cosmological signal. This problem has been addressed
already for the previous generation of radio telescopes, but the application to
SKA is different in many aspects.
In this chapter we summarise the contributions to the field of foreground
removal in the context of high redshift and high sensitivity 21cm
measurements. We use a stateoftheart simulation of the SKA Phase 1
observations complete with cosmological signal, foregrounds and
frequencydependent instrumental effects to test both parametric and
nonparametric foreground removal methods. We compare the recovered
cosmological signal using several different statistics and explore one of the
most exciting possibilities with the SKA  imaging of the ionized bubbles.
We find that with current methods it is possible to remove the foregrounds
with great accuracy and to get impressive power spectra and images of the
cosmological signal. The frequencydependent PSF of the instrument complicates
this recovery, so we resort to splitting the observation bandwidth into smaller
segments, each of a common resolution.
If the foregrounds are allowed a random variation from the smooth power law
along the line of sight, methods exploiting the smoothness of foregrounds or a
parametrization of their behaviour are challenged much more than nonparametric
ones. However, we show that correction techniques can be implemented to restore
the performances of parametric approaches, as long as the firstorder
approximation of a power law stands.

We provide an overview of the science benefits of combining information from
the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope
(LSST). We first summarise the capabilities and timeline of the LSST and
overview its science goals. We then discuss the science questions in common
between the two projects, and how they can be best addressed by combining the
data from both telescopes. We describe how weak gravitational lensing and
galaxy clustering studies with LSST and SKA can provide improved constraints on
the causes of the cosmological acceleration. We summarise the benefits to
galaxy evolution studies of combining deep optical multiband imaging with
radio observations. Finally, we discuss the excellent match between one of the
most unique features of the LSST, its temporal cadence in the optical waveband,
and the time resolution of the SKA.

HI intensity mapping (IM) is a novel technique capable of mapping the
largescale structure of the Universe in three dimensions and delivering
exquisite constraints on cosmology, by using HI as a biased tracer of the dark
matter density field. This is achieved by measuring the intensity of the
redshifted 21cm line over the sky in a range of redshifts without the
requirement to resolve individual galaxies. In this chapter, we investigate the
potential of SKA1 to deliver HI intensity maps over a broad range of
frequencies and a substantial fraction of the sky. By pinning down the baryon
acoustic oscillation and redshift space distortion features in the matter power
spectrum  thus determining the expansion and growth history of the Universe
 these surveys can provide powerful tests of dark energy models and
modifications to General Relativity. They can also be used to probe physics on
extremely large scales, where precise measurements of spatial curvature and
primordial nonGaussianity can be used to test inflation; on small scales, by
measuring the sum of neutrino masses; and at high redshifts where nonstandard
evolution models can be probed. We discuss the impact of foregrounds as well as
various instrumental and survey design parameters on the achievable
constraints. In particular we analyse the feasibility of using the SKA1
autocorrelations to probe the largescale signal.

The Square Kilometer Array (SKA) has the potential to produce galaxy redshift
surveys which will be competitive with other state of the art cosmological
experiments in the next decade. In this chapter we summarise what capabilities
the first and the second phases of the SKA will be able to achieve in its
current state of design. We summarise the different cosmological experiments
which are outlined in further detail in other chapters of this Science Book.
The SKA will be able to produce competitive Baryonic Oscillation (BAOs)
measurements in both its phases. The first phase of the SKA will provide
similar measurements as optical and IR experiments with completely different
systematic effects whereas the second phase being transformational in terms of
its statistical power. The SKA will produce very accurate Redshift Space
Distortions (RSD) measurements, being superior to other experiments at lower
redshifts, due to the large number of galaxies. Cross correlations of the
galaxy redshift data from the SKA with radio continuum surveys and optical
surveys will provide extremely good calibration of photometric redshifts as
well as extremely good bounds on modifications of gravity. Basing on a
Principle Component Analysis (PCA) approach, we find that the SKA will be able
to provide competitive constraints on dark energy and modified gravity models.
Due to the large area covered the SKA it will be a transformational experiment
in measuring physics from the largest scales such as nonGaussian signals from
$\textrm{f}_{\textrm{nl}}$. Finally, the SKA might produce the first real time
measurement of the redshift drift. The SKA will be a transformational machine
for cosmology as it grows from an early Phase 1 to its full power.

The new frontier of cosmology will be led by threedimensional surveys of the
largescale structure of the Universe. Based on its allsky surveys and
redshift depth, the SKA is destined to revolutionize cosmology, in combination
with future optical/ infrared surveys such as Euclid and LSST. Furthermore, we
will not have to wait for the full deployment of the SKA in order to see
transformational science. In the first phase of deployment (SKA1), allsky HI
intensity mapping surveys and allsky continuum surveys are forecast to be at
the forefront on the major questions of cosmology. We give a broad overview of
the major contributions predicted for the SKA. The SKA will not only deliver
precision cosmology  it will also probe the foundations of the standard model
and open the door to new discoveries on largescale features of the Universe.

21cm intensity mapping experiments aim to observe the diffuse neutral
hydrogen (HI) distribution on large scales which traces the Cosmic structure.
The Square Kilometre Array (SKA) will have the capacity to measure the 21cm
signal over a large fraction of the sky. However, the redshifted 21cm signal in
the respective frequencies is faint compared to the Galactic foregrounds
produced by synchrotron and freefree electron emission. In this article, we
review selected foreground subtraction methods suggested to effectively
separate the 21cm signal from the foregrounds with intensity mapping
simulations or data. We simulate an intensity mapping experiment feasible with
SKA phase 1 including extragalactic and Galactic foregrounds. We give an
example of the residuals of the foreground subtraction with a independent
component analysis and show that the angular power spectrum is recovered within
the statistical errors on most scales. Additionally, the scale of the Baryon
Acoustic Oscillations is shown to be unaffected by foreground subtraction.

By the time that the first phase of the Square Kilometre Array is deployed it
will be able to perform state of the art Large Scale Structure (LSS) as well as
Weak Gravitational Lensing (WGL) measurements of the distribution of matter in
the Universe. In this chapter we concentrate on the synergies that result from
crosscorrelating these different SKA data products as well as external
correlation with the weak lensing measurements available from CMB missions. We
show that the Dark Energy figures of merit obtained individually from WGL/LSS
measurements and their independent combination is significantly increased when
their full crosscorrelations are taken into account. This is due to the
increased knowledge of galaxy bias as a function of redshift as well as the
extra information from the different cosmological dependences of the
crosscorrelations. We show that the crosscorrelation between a spectroscopic
LSS sample and a weak lensing sample with photometric redshifts can calibrate
these same photometric redshifts, and their scatter, to high accuracy by
modelling them as nuisance parameters and fitting them simultaneously
cosmology. Finally we show that Modified Gravity parameters are greatly
constrained by this crosscorrelations because weak lensing and redshift space
distortions (from the LSS survey) break strong degeneracies in common
parameterisations of modified gravity.

Several experiments are underway to detect the cosmic redshifted 21cm signal
from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very
low signaltonoise ratio, these observations aim for a statistical detection
of the signal by measuring its power spectrum. We investigate the extraction of
the variance of the signal as a first step towards detecting and constraining
the global history of the EoR. Signal variance is the integral of the signal's
power spectrum, and it is expected to be measured with a high significance. We
demonstrate this through results from a simulation and parameter estimation
pipeline developed for the Low Frequency Array (LOFAR)EoR experiment. We show
that LOFAR should be able to detect the EoR in 600 hours of integration using
the variance statistic. Additionally, the redshift ($z_r$) and duration
($\Delta z$) of reionization can be constrained assuming a parametrization. We
use an EoR simulation of $z_r = 7.68$ and $\Delta z = 0.43$ to test the
pipeline. We are able to detect the simulated signal with a significance of 4
standard deviations and extract the EoR parameters as $z_r =
7.72^{+0.37}_{0.18}$ and $\Delta z = 0.53^{+0.12}_{0.23}$ in 600 hours,
assuming that systematic errors can be adequately controlled. We further show
that the significance of detection and constraints on EoR parameters can be
improved by measuring the crossvariance of the signal by crosscorrelating
consecutive redshift bins.

In this paper we introduce PkANN, a freely available software package for
interpolating the nonlinear matter power spectrum, constructed using
Artificial Neural Networks (ANNs). Previously, using Halofit to calculate
matter power spectrum, we demonstrated that ANNs can make extremely quick and
accurate predictions of the power spectrum. Now, using a suite of 6380 Nbody
simulations spanning 580 cosmologies, we train ANNs to predict the power
spectrum over the cosmological parameter space spanning $3\sigma$ confidence
level (CL) around the concordance cosmology. When presented with a set of
cosmological parameters ($\Omega_{\rm m} h^2, \Omega_{\rm b} h^2, n_s, w,
\sigma_8, \sum m_\nu$ and redshift $z$), the trained ANN interpolates the power
spectrum for $z\leq2$ at subper cent accuracy for modes up to
$k\leq0.9\,h\textrm{Mpc}^{1}$. PkANN is faster than computationally expensive
Nbody simulations, yet provides a worstcase error $<1$ per cent fit to the
nonlinear matter power spectrum deduced through Nbody simulations. The
overall precision of PkANN is set by the accuracy of our Nbody simulations, at
5 per cent level for cosmological models with $\sum m_\nu<0.5$ eV for all
redshifts $z\leq2$. For models with $\sum m_\nu>0.5$ eV, predictions are
expected to be at 5 (10) per cent level for redshifts $z>1$ ($z\leq1$). The
PkANN interpolator may be freely downloaded from
http://zuserver2.star.ucl.ac.uk/~fba/PkANN

The combination of multiple cosmological probes can produce measurements of
cosmological parameters much more stringent than those possible with any
individual probe. We examine the combination of two highly correlated probes of
latetime structure growth: (i) weak gravitational lensing from a survey with
photometric redshifts and (ii) galaxy clustering and redshift space distortions
from a survey with spectroscopic redshifts. We choose generic survey designs so
that our results are applicable to a range of current and future photometric
redshift (e.g. KiDS, DES, HSC, Euclid) and spectroscopic redshift (e.g. DESI,
4MOST, Sumire) surveys. Combining the surveys greatly improves their power to
measure both dark energy and modified gravity. An independent, nonoverlapping
combination sees a dark energy figure of merit more than 4 times larger than
that produced by either survey alone. The powerful synergies between the
surveys are strongest for modified gravity, where their constraints are
orthogonal, producing a nonoverlapping joint figure of merit nearly 2 orders
of magnitude larger than either alone. Our projected angular power spectrum
formalism makes it easy to model the crosscorrelation observable when the
surveys overlap on the sky, producing a joint data vector and full covariance
matrix. We calculate a samesky improvement factor, from the inclusion of these
crosscorrelations, relative to nonoverlapping surveys. We find nearly a
factor of 4 for dark energy and more than a factor of 2 for modified gravity.
The exact forecast figures of merit and samesky benefits can be radically
affected by a range of forecasts assumption, which we explore methodically in a
sensitivity analysis. We show that that our fiducial assumptions produce robust
results which give a good average picture of the science return from combining
photometric and spectroscopic surveys.

One of the most promising ways to study the epoch of reionization (EoR) is
through radio observations of the redshifted 21cm line emission from neutral
hydrogen. These observations are complicated by the fact that the mapping of
redshifts to lineofsight positions is distorted by the peculiar velocities of
the gas. Such distortions can be a source of error if they are not properly
understood, but they also encode information about cosmology and astrophysics.
We study the effects of redshift space distortions on the power spectrum of
21cm radiation from the EoR using large scale $N$body and radiative transfer
simulations. We quantify the anisotropy introduced in the 21cm power spectrum
by redshift space distortions and show how it evolves as reionization
progresses and how it relates to the underlying physics. We go on to study the
effects of redshift space distortions on LOFAR observations, taking instrument
noise and foreground subtraction into account. We find that LOFAR should be
able to directly observe the power spectrum anisotropy due to redshift space
distortions at spatial scales around $k \sim 0.1$ Mpc$^{1}$ after $\gtrsim$
1000 hours of integration time. At larger scales, sample errors become a
limiting factor, while at smaller scales detector noise and foregrounds make
the extraction of the signal problematic. Finally, we show how the
astrophysical information contained in the evolution of the anisotropy of the
21cm power spectrum can be extracted from LOFAR observations, and how it can
be used to distinguish between different reionization scenarios.

We study the arrival directions of 69 ultrahigh energy cosmic rays (UHECRs)
observed at the Pierre Auger Observatory (PAO) with energies exceeding 55 EeV.
We investigate whether the UHECRs exhibit the anisotropy signal expected if the
primary particles are protons that originate in galaxies in the local universe,
or in sources correlated with these galaxies. We crosscorrelate the UHECR
arrival directions with the positions of IRASPSCz and 2MASS6dF galaxies
taking into account particle energy losses during propagation. This is the
first time that the 6dF survey is used in a search for the sources of UHECRs
and the first time that the PSCz survey is used with the full 69 PAO events.
The observed crosscorrelation signal is larger for the PAO UHECRs than for 94%
(98%) of realisations from an isotropic distribution when crosscorrelated with
the PSCz (6dF). On the other hand the observed crosscorrelation signal is
lower than that expected from 85% of realisations, had the UHECRs originated in
galaxies in either survey. The observed crosscorrelation signal does exceed
that expected by 50% of the realisations if the UHECRs are randomly deflected
by intervening magnetic fields by 5 degrees or more. We propose a new method of
analysing the expected anisotropy signal, by dividing the predicted UHECR
source distribution into equal predicted flux radial shells, which can help
localise and constrain the properties of UHECR sources. We find that the 69 PAO
events are consistent with isotropy in the nearest of three shells we define,
whereas there is weak evidence for correlation with the predicted source
distribution in the two more distant shells in which the galaxy distribution is
less anisotropic.

We study the impact of the spread spectrum effect in radio interferometry on
the quality of image reconstruction. This spread spectrum effect will be
induced by the wide fieldofview of forthcoming radio interferometric
telescopes. The resulting chirp modulation improves the quality of
reconstructed interferometric images by increasing the incoherence of the
measurement and sparsity dictionaries. We extend previous studies of this
effect to consider the more realistic setting where the chirp modulation varies
for each visibility measurement made by the telescope. In these first
preliminary results, we show that for this setting the quality of
reconstruction improves significantly over the case without chirp modulation
and achieves almost the reconstruction quality of the case of maximal, constant
chirp modulation.

The accurate and precise removal of 21cm foregrounds from Epoch of
Reionization redshifted 21cm emission data is essential if we are to gain
insight into an unexplored cosmological era. We apply a nonparametric
technique, Generalized Morphological Component Analysis or GMCA, to simulated
LOFAREoR data and show that it has the ability to clean the foregrounds with
high accuracy. We recover the 21cm 1D, 2D and 3D power spectra with high
accuracy across an impressive range of frequencies and scales. We show that
GMCA preserves the 21cm phase information, especially when the smallest
spatial scale data is discarded. While it has been shown that LOFAREoR image
recovery is theoretically possible using image smoothing, we add that wavelet
decomposition is an efficient way of recovering 21cm signal maps to the same
or greater order of accuracy with more flexibility. By comparing the GMCA
output residual maps (equal to the noise, 21cm signal and any foreground
fitting errors) with the 21cm maps at one frequency and discarding the smaller
wavelet scale information, we find a correlation coefficient of 0.689, compared
to 0.588 for the equivalently smoothed image. Considering only the central 50%
of the maps, these coefficients improve to 0.905 and 0.605 respectively and we
conclude that wavelet decomposition is a significantly more powerful method to
denoise reconstructed 21cm maps than smoothing.

We investigate the interpolation of power spectra of matter fluctuations
using Artificial Neural Network (PkANN). We present a new approach to confront
smallscale nonlinearities in the power spectrum of matter fluctuations. This
everpresent and pernicious uncertainty is often the Achilles' heel in
cosmological studies and must be reduced if we are to see the advent of
precision cosmology in the latetime Universe. We show that an optimally
trained artificial neural network (ANN), when presented with a set of
cosmological parameters (Omega_m h^2, Omega_b h^2, n_s, w_0, sigma_8, m_nu and
redshift z), can provide a worstcase error <=1 per cent (for z<=2) fit to the
nonlinear matter power spectrum deduced through Nbody simulations, for modes
up to k<=0.7 h/Mpc. Our power spectrum interpolator is accurate over the entire
parameter space. This is a significant improvement over some of the current
matter power spectrum calculators. In this paper, we detail how an accurate
interpolation of the matter power spectrum is achievable with only a sparsely
sampled grid of cosmological parameters. Unlike largescale Nbody simulations
which are computationally expensive and/or infeasible, a welltrained ANN can
be an extremely quick and reliable tool in interpreting cosmological
observations and parameter estimation. This paper is the first in a series. In
this method paper, we generate the nonlinear matter power spectra using
HaloFit and use them as mock observations to train the ANN. This work sets the
foundation for Paper II, where a suite of Nbody simulations will be used to
compute the nonlinear matter power spectra at subper cent accuracy, in the
quasinonlinear regime 0.1 h/Mpc <= k <= 0.9 h/Mpc. A trained ANN based on
this Nbody suite will be released for the scientific community.

We introduce a new implementation of the FastICA algorithm on simulated LOFAR
EoR data with the aim of accurately removing the foregrounds and extracting the
21cm reionization signal. We find that the method successfully removes the
foregrounds with an average fitting error of 0.5 per cent and that the 2D and
3D power spectra are recovered across the frequency range. We find that for
scales above several PSF scales the 21cm variance is successfully recovered
though there is evidence of noise leakage into the reconstructed foreground
components. We find that this blind independent component analysis technique
provides encouraging results without the danger of prior foreground
assumptions.

We present forecasts for constraints on cosmological models which can be
obtained by forthcoming radio continuum surveys: the wide surveys with the LOw
Frequency ARray (LOFAR), Australian Square Kilometre Array Pathfinder (ASKAP)
and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). We
use simulated catalogues appropriate to the planned surveys to predict
measurements obtained with the source autocorrelation, the crosscorrelation
between radio sources and CMB maps (the Integrated SachsWolfe effect), the
crosscorrelation of radio sources with foreground objects due to cosmic
magnification, and a joint analysis together with the CMB power spectrum and
supernovae. We show that near future radio surveys will bring complementary
measurements to other experiments, probing different cosmological volumes, and
having different systematics. Our results show that the unprecedented sky
coverage of these surveys combined should provide the most significant
measurement yet of the Integrated SachsWolfe effect. In addition, we show that
using the ISW effect will significantly tighten constraints on modified gravity
parameters, while the best measurements of dark energy models will come from
galaxy autocorrelation function analyses. Using the combination of EMU and
WODAN to provide a full sky survey, it will be possible to measure the dark
energy parameters with an uncertainty of \{$\sigma (w_0) = 0.05$, $\sigma (w_a)
= 0.12$\} and the modified gravity parameters \{$\sigma (\eta_0) = 0.10$,
$\sigma (\mu_0) = 0.05$\}, assuming Planck CMB+SN(current data) priors.
Finally, we show that radio surveys would detect a primordial nonGaussianity
of $f_{\rm NL}$ = 8 at 1$\sigma$ and we briefly discuss other promising
probes.

We observe a large excess of power in the statistical clustering of Luminous
Red Galaxies in the photometric SDSS galaxy sample called MegaZ DR7. This is
seen over the lowest multipoles in the angular power spectra C_{\ell} in four
equally spaced redshift bins between 0.45 < z < 0.65. However, it is most
prominent in the highest redshift band at ~ 4 sigma and it emerges at an
effective scale k ~ 0.01 h Mpc^{1}. Given that MegaZ DR7 is the largest cosmic
volume galaxy survey to date (3.3 (Gpc h^{1})^3) this implies an anomaly on
the largest physical scales probed by galaxies. Alternatively, this signature
could be a consequence of it appearing at the most systematically susceptible
redshift. There are several explanations for this excess power that range from
systematics to new physics. This could have important consequences for the next
generation of galaxy surveys or the LCDM model. We test the survey, data and
excess power, as well as possible origins.