
We use the results of previous work building a halo model formalism for the
distribution of neutral hydrogen, along with experimental parameters of future
radio facilities, to place forecasts on astrophysical and cosmological
parameters from next generation surveys. We consider 21 cm intensity mapping
surveys conducted using the BINGO, CHIME, FAST, TianLai, MeerKAT and SKA
experimental configurations. We work with the 5parameter cosmological dataset
of {$\Omega_m, \sigma_8, h, n_s, \Omega_b$} assuming a flat $\Lambda$CDM model,
and the astrophysical parameters {$v_{c,0}, \beta$} which represent the cutoff
and slope of the HI halo mass relation. We explore (i) quantifying the effects
of the astrophysics on the recovery of the cosmological parameters, (ii) the
dependence of the cosmological forecasts on the details of the astrophysical
parametrization, and (iii) the improvement of the constraints on probing
smaller scales in the HI power spectrum. For an SKA I MID intensity mapping
survey alone, probing scales up to $\ell_{\rm max} = 1000$, we find a factor of
$1.1  1.3$ broadening in the constraints on $\Omega_b$ and $\Omega_m$, and of
$2.4  2.6$ on $h$, $n_s$ and $\sigma_8$, if we marginalize over astrophysical
parameters without any priors. However, even the prior information coming from
the present knowledge of the astrophysics largely alleviates this broadening.
These findings do not change significantly on considering an extended HIHM
relation, illustrating the robustness of the results to the choice of the
astrophysical parametrization. Probing scales up to $\ell_{\rm max} = 2000$
improves the constraints by factors of 1.51.8. The forecasts improve on
increasing the number of tomographic redshift bins, saturating, in many cases,
with 4  5 redshift bins. We also forecast constraints for intensity mapping
with other experiments, and draw similar conclusions.

We present the gaugeinvariant formalism of cosmological weak lensing,
accounting for all the relativistic effects due to the scalar, vector, and
tensor perturbations at the linear order. While the light propagation is fully
described by the geodesic equation, the relation of the photon wavevector to
the physical quantities requires the specification of the frames, where they
are defined. By constructing the local tetrad bases at the observer and the
source positions, we clarify the relation of the weak lensing observables such
as the convergence, the shear, and the rotation to the physical size and shape
defined in the source restframe and the observed angle and redshift measured
in the observer restframe. Compared to the standard lensing formalism,
additional relativistic effects contribute to all the lensing observables. We
explicitly verify the gaugeinvariance of the lensing observables and compare
our results to previous work. In particular, we demonstrate that even in the
presence of the vector and tensor perturbations, the physical rotation of the
lensing observables vanishes at the linear order, while the tetrad basis
rotates along the light propagation compared to a FRW coordinate. Though the
latter is often used as a probe of primordial gravitational waves, the rotation
of the tetrad basis is indeed not a physical observable. We further clarify its
relation to the EB decomposition in weak lensing. Our formalism provides a
transparent and comprehensive perspective of cosmological weak lensing.

We study the distribution of dark matter in the nonlinear regime in a model
in which the primordial fluctuations include, in addition to the dominant
primordial Gaussian fluctuations generated by the standard $\Lambda CDM$
cosmological model, the effects of a cosmic string wake set up at the time of
equal matter and radiation, making use of cosmological $N$body simulations. At
early times the string wake leads to a planar overdensity of dark matter. We
study how this nonGaussian pattern of a cosmic string wake evolves in the
presence of the Gaussian perturbations, making use of wavelet and ridgeletlike
statistics specifically designed to extract string wake signals. At late times
the Gaussian fluctuations disrupt the string wake. We find that for a string
tension of $G \mu = 10^{7}$, a value just below the current observational
limit, the effects of a string wake can be identified in the dark matter
distribution, using the current level of the statistical analysis, down to a
redshift of $z = 10$.

Weak lensing peak counts are a powerful statistical tool for constraining
cosmological parameters. So far, this method has been applied only to surveys
with relatively small areas, up to several hundred square degrees. As future
surveys will provide weak lensing datasets with size of thousands of square
degrees, the demand on the theoretical prediction of the peak statistics will
become heightened. In particular, large simulations of increased cosmological
volume are required. In this work, we investigate the possibility of using
simulations generated with the fast ComovingLagrangian acceleration (COLA)
method, coupled to the convergence map generator Ufalcon, for predicting the
peak counts. We examine the systematics introduced by the COLA method by
comparing it with a full TreePM code. We find that for a 2000 deg$^2$ survey,
the systematic error is much smaller than the statistical error. This suggests
that the COLA method is able to generate promising theoretical predictions for
weak lensing peaks. We also examine the constraining power of various
configurations of data vectors, exploring the influence of splitting the sample
into tomographic bins and combining different smoothing scales. We find the
combination of smoothing scales to have the most constraining power, improving
the constraints on the $S_8$ amplitude parameter by at least 40% compared to a
single smoothing scale, with tomography brining only limited increase in
measurement precision.

Galaxy spectra are essential to probe the spatial distribution of galaxies in
our Universe. To better interpret current and future spectroscopic galaxy
redshift surveys, it is important to be able to simulate these data sets. We
describe Uspec, a forward modeling tool to generate galaxy spectra taking into
account intrinsic galaxy properties as well as instrumental responses of a
given telescope. The model for the intrinsic properties of the galaxy
population was developed in an earlier work for broadband imaging surveys [1].
We apply Uspec to the SDSS/CMASS sample of Luminous Red Galaxies (LRGs). We
construct selection cuts that match those used to build this LRG sample, which
we then apply to data and simulations in the same way. The resulting real and
simulated average spectra show a very good agreement overall, with the
simulated one showing a slightly bluer galaxy population. For a quantitative
comparison, we perform Principal Component Analysis (PCA) of the sets of
spectra. By comparing the PCs constructed from simulations and data, we find
very good agreement for the first four components, and moderate for the fifth.
The distributions of the eigencoefficients also show an appreciable overlap. We
are therefore able to properly simulate the LRG sample taking into account the
SDSS/BOSS instrumental responses. The small residual differences between the
two samples can be ascribed to the intrinsic properties of the simulated galaxy
population, which can be reduced by adjusting the model parameters in the
future. This provides good prospects for the forward modeling of upcoming large
spectroscopic surveys.

Dark matter in the universe evolves through gravity to form a complex network
of halos, filaments, sheets and voids, that is known as the cosmic web.
Computational models of the underlying physical processes, such as classical
Nbody simulations, are extremely resource intensive, as they track the action
of gravity in an expanding universe using billions of particles as tracers of
the cosmic matter distribution. Therefore, upcoming cosmology experiments will
face a computational bottleneck that may limit the exploitation of their full
scientific potential. To address this challenge, we demonstrate the application
of a machine learning technique called Generative Adversarial Networks (GAN) to
learn models that can efficiently generate new, physically realistic
realizations of the cosmic web. Our training set is a small, representative
sample of 2D image snapshots from Nbody simulations of size 500 and 100 Mpc.
We show that the GANproduced results are qualitatively and quantitatively very
similar to the originals. Generation of a new cosmic web realization with a GAN
takes a fraction of a second, compared to the many hours needed by the Nbody
technique. We anticipate that GANs will therefore play an important role in
providing extremely fast and precise simulations of cosmic web in the era of
large cosmological surveys, such as Euclid and LSST.

Modeling the Point Spread Function (PSF) of widefield surveys is vital for
many astrophysical applications and cosmological probes including weak
gravitational lensing. The PSF smears the image of any recorded object and
therefore needs to be taken into account when inferring properties of galaxies
from astronomical images. In the case of cosmic shear, the PSF is one of the
dominant sources of systematic errors and must be treated carefully to avoid
biases in cosmological parameters. Recently, forward modeling approaches to
calibrate shear measurements within the MonteCarlo Control Loops ($MCCL$)
framework have been developed. These methods typically require simulating a
large amount of widefield images, thus, the simulations need to be very fast
yet have realistic properties in key features such as the PSF pattern. Hence,
such forward modeling approaches require a very flexible PSF model, which is
quick to evaluate and whose parameters can be estimated reliably from survey
data. We present a PSF model that meets these requirements based on a fast
deeplearning method to estimate its free parameters. We demonstrate our
approach on publicly available SDSS data. We extract the most important
features of the SDSS sample via principal component analysis. Next, we
construct our model based on perturbations of a fixed base profile, ensuring
that it captures these features. We then train a Convolutional Neural Network
to estimate the free parameters of the model from noisy images of the PSF. This
allows us to render a model image of each star, which we compare to the SDSS
stars to evaluate the performance of our method. We find that our approach is
able to accurately reproduce the SDSS PSF at the pixel level, which, due to the
speed of both the model evaluation and the parameter estimation, offers good
prospects for incorporating our method into the $MCCL$ framework.

Upcoming weak lensing surveys will probe large fractions of the sky with
unprecedented accuracy. To infer cosmological constraints, a large ensemble of
survey simulations are required to accurately model cosmological observables
and their covariances. We develop a parallelized multilensplane pipeline
called UFalcon, designed to generate fullsky weak lensing maps from lightcones
within a minimal runtime. It makes use of LPICOLA, an approximate numerical
code, which provides a fast and accurate alternative to cosmological $N$Body
simulations. The UFalcon maps are constructed by nesting 2 simulations covering
a redshiftrange from $z=0.1$ to $1.5$ without replicating the simulation
volume. We compute the convergence and projected overdensity maps for LPICOLA
in the lightcone or snapshot mode. The generation of such a map, including the
LPICOLA simulation, takes about 3 hours walltime on 220 cores. We use the maps
to calculate the spherical harmonic power spectra, which we compare to
theoretical predictions and to UFalcon results generated using the full
$N$Body code GADGET2. We then compute the covariance matrix of the fullsky
spherical harmonic power spectra using 150 UFalcon maps based on LPICOLA in
lightcone mode. We consider the PDF, the higherorder moments and the variance
of the smoothed field variance to quantify the accuracy of the covariance
matrix, which we find to be a few percent for scales $\ell \sim 10^2$ to
$10^3$. We test the impact of this level of accuracy on cosmological
constraints using an optimistic survey configuration, and find that the final
results are robust to this level of uncertainty. The speed and accuracy of our
developed pipeline provides a basis to also include further important features
such as masking, varying noise and will allow us to compute covariance matrices
for models beyond $\Lambda$CDM. [abridged]

We explore a new technique to measure cosmic shear using Einstein rings. In
Birrer et al. (2017), we showed that the detailed modelling of Einstein rings
can be used to measure external shear to high precision. In this letter, we
explore how a collection of Einstein rings can be used as a statistical probe
of cosmic shear. We present a forecast of the cosmic shear information
available in Einstein rings for different strong lensing survey configurations.
We find that, assuming that the number density of Einstein rings in the COSMOS
survey is representative, future strong lensing surveys should have a
cosmological precision comparable to the current ground based weak lensing
surveys. We discuss how this technique is complementary to the standard cosmic
shear analyses since it is sensitive to different systematic and can be used
for crosscalibration.

Assessing the consistency of parameter constraints derived from different
cosmological probes is an important way to test the validity of the underlying
cosmological model. In an earlier work [Nicola et al., 2017], we computed
constraints on cosmological parameters for $\Lambda$CDM from an integrated
analysis of CMB temperature anisotropies and CMB lensing from Planck, galaxy
clustering and weak lensing from SDSS, weak lensing from DES SV as well as Type
Ia supernovae and Hubble parameter measurements. In this work, we extend this
analysis and quantify the concordance between the derived constraints and those
derived by the Planck Collaboration as well as WMAP9, SPT and ACT. As a measure
for consistency, we use the Surprise statistic [Seehars et al., 2014], which is
based on the relative entropy. In the framework of a flat $\Lambda$CDM
cosmological model, we find all data sets to be consistent with one another at
a level of less than 1$\sigma$. We highlight that the relative entropy is
sensitive to inconsistencies in the models that are used in different parts of
the analysis. In particular, inconsistent assumptions for the neutrino mass
break its invariance on the parameter choice. When consistent model assumptions
are used, the data sets considered in this work all agree with each other and
$\Lambda$CDM, without evidence for tensions.

As widefield surveys yield ever more precise measurements, cosmology has
entered a phase of high precision requiring highly accurate and fast
theoretical predictions. At the heart of most cosmological model predictions is
a numerical solution of the EinsteinBoltzmann equations governing the
evolution of linear perturbations in the Universe. We present PyCosmo, a new
Pythonbased framework to solve this set of equations using a special pur pose
solver based on symbolic manipulations, automatic generation of C++ code and
sparsity optimisation. The code uses a consistency relation of the field
equations to adapt the time step and does not rely on physical approximations
for speedup. After reviewing the system of firstorder linear homogeneous
differential equations to be solved, we describe the numerical scheme
implemented in PyCosmo. We then compare the predictions and performance of the
code for the computation of the transfer functions of cosmological
perturbations and compare it to existing cosmological Boltzmann codes. We find
that we achieve comparable execution times for comparable accuracies. While
PyCosmo does not yet have all the features of other codes, our approach is
complementary to existing cosmological Boltzmann solvers and can be used as an
independent test of their numerical solutions. The symbolic representation of
the EinsteinBoltzmann equation system in PyCosmo provides a convenient
interface for implementing extended cosmological models. We also discuss how
the PyCosmo framework can also be used as a general framework to compute
cosmological quantities as well as observables for both interactive and
highperformance batch jobs applications. Information about the PyCosmo package
and future code releases are available at
http://www.cosmology.ethz.ch/research/softwarelab.html.

Determining the redshift distribution $n(z)$ of galaxy samples is essential
for several cosmological probes including weak lensing. For imaging surveys,
this is usually done using photometric redshifts estimated on an
objectbyobject basis. We present a new approach for directly measuring the
global $n(z)$ of cosmological galaxy samples, including uncertainties, using
forward modeling. Our method relies on image simulations produced using UFig
(Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation)
within the $MCCL$ (MonteCarlo Control Loops) framework. The galaxy population
is modeled using parametric forms for the luminosity functions, spectral energy
distributions, sizes and radial profiles of both blue and red galaxies. We
apply exactly the same analysis to the real data and to the simulated images,
which also include instrumental and observational effects. By adjusting the
parameters of the simulations, we derive a set of acceptable models that are
statistically consistent with the data. We then apply the same cuts to the
simulations that were used to construct the target galaxy sample in the real
data. The redshifts of the galaxies in the resulting simulated samples yield a
set of $n(z)$ distributions for the acceptable models. We demonstrate the
method by determining $n(z)$ for a cosmic shear like galaxy sample from the
4band Subaru SuprimeCam data in the COSMOS field. We also complement this
imaging data with a spectroscopic calibration sample from the VVDS survey. We
compare our resulting posterior $n(z)$ distributions to the one derived from
photometric redshifts estimated using 36 photometric bands in COSMOS and find
good agreement. This offers good prospects for applying our approach to current
and future large imaging surveys.

Approximate Bayesian Computation (ABC) is a method to obtain a posterior
distribution without a likelihood function, using simulations and a set of
distance metrics. For that reason, it has recently been gaining popularity as
an analysis tool in cosmology and astrophysics. Its drawback, however, is a
slow convergence rate. We propose a novel method, which we call qABC, to
accelerate ABC with Quantile Regression. In this method, we create a model of
quantiles of distance measure as a function of input parameters. This model is
trained on a small number of simulations and estimates which regions of the
prior space are likely to be accepted into the posterior. Other regions are
then immediately rejected. This procedure is then repeated as more simulations
are available. We apply it to the practical problem of estimation of redshift
distribution of cosmological samples, using forward modelling developed in
previous work. The qABC method converges to nearly same posterior as the basic
ABC. It uses, however, only 20\% of the number of simulations compared to basic
ABC, achieving a fivefold gain in execution time for our problem. For other
problems the acceleration rate may vary; it depends on how close the prior is
to the final posterior. We discuss possible improvements and extensions to this
method.

Weak Gravitational Lensing is a powerful probe of the dark sector of the
Universe. One of the main challenges for this technique is the treatment of
systematics in the measurement of cosmic shear from galaxy shapes. In an
earlier work, Refregier & Amara (2014) have proposed the Monte Carlo Control
Loops (MCCL) to overcome these effects using a forward modeling approach. We
focus here on one of the control loops in this method, the task of which is the
calibration of the shear measurement. For this purpose, we first consider the
requirements on the shear systematics for a given survey and propagate them to
different systematics terms. We use two onepoint statistics to calibrate the
shear measurement and six further onepoint statistics as diagnostics. We also
propagate the systematics levels that we estimate from the onepoint functions
to the twopoint functions for the different systematic error sources. This
allows us to assess the consistency between the systematics levels measured in
different ways. To test the method, we construct synthetic sky surveys with an
area of 1,700 deg$^2$. With some simplifying assumptions, we are able to meet
the requirements on the shear calibration for this survey configuration.
Furthermore, we account for the total residual shear systematics in terms of
the contributing sources. We discuss how this MCCL framework can be applied to
current and future weak lensing surveys.

We demonstrate the potential of Deep Learning methods for measurements of
cosmological parameters from density fields, focusing on the extraction of
nonGaussian information. We consider weak lensing mass maps as our dataset. We
aim for our method to be able to distinguish between five models, which were
chosen to lie along the $\sigma_8$  $\Omega_m$ degeneracy, and have nearly the
same twopoint statistics. We design and implement a Deep Convolutional Neural
Network (DCNN) which learns the relation between five cosmological models and
the mass maps they generate. We develop a new training strategy which ensures
the good performance of the network for high levels of noise. We compare the
performance of this approach to commonly used nonGaussian statistics, namely
the skewness and kurtosis of the convergence maps. We find that our
implementation of DCNN outperforms the skewness and kurtosis statistics,
especially for high noise levels. The network maintains the mean discrimination
efficiency greater than $85\%$ even for noise levels corresponding to ground
based lensing observations, while the other statistics perform worse in this
setting, achieving efficiency less than $70\%$. This demonstrates the ability
of CNNbased methods to efficiently break the $\sigma_8$  $\Omega_m$
degeneracy with weak lensing mass maps alone. We discuss the potential of this
method to be applied to the analysis of real weak lensing data and other
datasets.

We study the substructure content of the strong gravitational lens
RXJ11311231 through a forward modelling approach that relies on generating an
extensive suite of realistic simulations. We use a semianalytic merger tree
prescription that allows us to stochastically generate substructure populations
whose properties depend on the dark matter particle mass. These synthetic halos
are then used as lenses to produce realistic mock images that have the same
features, e.g. luminous arcs, quasar positions, instrumental noise and PSF, as
the data. We then analyse the data and the simulations in the same way with
summary statistics that are sensitive to the signal being targeted and are able
to constrain models of dark matter statistically using Approximate Bayesian
Computing (ABC) techniques. In this work, we focus on the thermal relic mass
estimate and fix the semianalytic descriptions of the substructure evolution
based on recent literature. We are able, based on the HST data for
RXJ11311231, to rule out a warm dark matter thermal relic mass below 2 keV at
the 2$\sigma$ confidence level.

We extend the results of previous analyses towards constraining the abundance
and clustering of postreionization ($z \sim 05$) neutral hydrogen (HI)
systems using a halo model framework. We work with a comprehensive HI dataset
including the smallscale clustering, column density and mass function of HI
galaxies at low redshifts, intensity mapping measurements at intermediate
redshifts and the UV/optical observations of Damped Lyman Alpha (DLA) systems
at higher redshifts. We use a Markov Chain Monte Carlo (MCMC) approach to
constrain the parameters of the bestfitting models, both for the HIhalo mass
relation and the HI radial density profile. We find that a radial exponential
profile results in a good fit to the lowredshift HI observations, including
the clustering and the column density distribution. The form of the profile is
also found to match the highredshift DLA observations, when used in
combination with a threeparameter HIhalo mass relation and a redshift
evolution in the HI concentration. The halo model predictions are in good
agreement with the observed HI surface density profiles of lowredshift
galaxies, and the general trends in the the impact parameter and covering
fraction observations of highredshift DLAs. We provide convenient tables
summarizing the bestfit halo model predictions.

Recent progress in cosmology has relied on combining different cosmological
probes. In earlier work, we implemented an integrated approach to cosmology
where the probes are combined into a common framework at the map level. This
has the advantage of taking full account of the correlations between the
different probes, to provide a stringent test of systematics and of the
validity of the cosmological model. We extend this analysis to include not only
CMB temperature, galaxy clustering, weak lensing from SDSS but also CMB
lensing, weak lensing from the DES SV survey, Type Ia SNe and $H_{0}$
measurements. This yields 12 auto and crosspower spectra as well as
background probes. Furthermore, we extend the treatment of systematic
uncertainties. For $\Lambda$CDM, we find results that are consistent with our
earlier work. Given our enlarged data set and systematics treatment, this
confirms the robustness of our analysis and results. Furthermore, we find that
our bestfit cosmological model gives a good fit to the data we consider with
no signs of tensions within our analysis. We also find our constraints to be
consistent with those found by WMAP9, SPT and ACT and the KiDS weak lensing
survey. Comparing with the Planck Collaboration results, we see a broad
agreement, but there are indications of a tension from the marginalized
constraints in most pairs of cosmological parameters. Since our analysis
includes CMB temperature Planck data at $10 < \ell < 610$, the tension appears
to arise between the Planck high$\ell$ and the other measurements.
Furthermore, we find the constraints on the probe calibration parameters to be
in agreement with expectations, showing that the data sets are mutually
consistent. In particular, this yields a confirmation of the amplitude
calibration of the weak lensing measurements from SDSS, DES SV and Planck CMB
lensing from our integrated analysis. [abridged]

The Xray emission from active galactic nucleus (AGN) is a major component of
extragalactic Xray sky. In this paper, we use the Xray luminosity function
(XLF) and halo occupation distribution (HOD) formalism to construct a halo
model for the Xray emission from AGNs. Verifying that the two inputs (XLF and
HOD) are in agreement with each other, we compute the autocorrelation power
spectrum in the soft Xray band (0.52 keV) due to the AGNs potentially
resolved by eROSITA (extended ROentgen Survey with an Imaging Telescope Array)
mission and explore the redshift and mass dependence of the power spectrum.
Studying the relative contribution of the Poisson and the clustering terms to
the total power, we find that at multipoles $l\lesssim 1000$ (i.e. large
scales), the clustering term is larger than the Poisson term. We also forecast
the potential of Xray autocorrelation power spectrum and Xraylensing
crosscorrelation power spectrum using eROSITA and eROSITALSST (Large Synoptic
Survey Telescope) surveys, respectively, to constrain the HOD parameters and
their redshift evolution. In addition, we compute the power spectrum of the
AGNs lying below the flux resolution limit of eROSITA, which is essential to
understand in order to extract the Xray signal from the hot diffuse gas
present in galaxies and clusters.

We propose a novel approach for mitigating radio frequency interference (RFI)
signals in radio data using the latest advances in deep learning. We employ a
special type of Convolutional Neural Network, the UNet, that enables the
classification of clean signal and RFI signatures in 2D timeordered data
acquired from a radio telescope. We train and assess the performance of this
network using the HIDE & SEEK radio data simulation and processing packages, as
well as early Science Verification data acquired with the 7m singledish
telescope at the Bleien Observatory. We find that our UNet implementation is
showing competitive accuracy to classical RFI mitigation algorithms such as
SEEK's SumThreshold implementation. We publish our UNet software package on
GitHub under GPLv3 license.

As several large singledish radio surveys begin operation within the coming
decade, a wealth of radio data will become available and provide a new window
to the Universe. In order to fully exploit the potential of these data sets, it
is important to understand the systematic effects associated with the
instrument and the analysis pipeline. A common approach to tackle this is to
forwardmodel the entire system  from the hardware to the analysis of the data
products. For this purpose, we introduce two newly developed, opensource
Python packages: the HI Data Emulator (HIDE) and the Signal Extraction and
Emission Kartographer (SEEK) for simulating and processing singledish radio
survey data. HIDE forwardmodels the process of collecting astronomical radio
signals in a singledish radio telescope instrument and outputs pixellevel
timeordereddata. SEEK processes the timeordereddata, removes artifacts from
Radio Frequency Interference (RFI), automatically applies flux calibration, and
aims to recover the astronomical radio signal. The two packages can be used
separately or together depending on the application. Their modular and flexible
nature allows easy adaptation to other instruments and data sets. We describe
the basic architecture of the two packages and examine in detail the noise and
RFI modeling in HIDE, as well as the implementation of gain calibration and RFI
mitigation in SEEK. We then apply HIDE & SEEK to forwardmodel a Galactic
survey in the frequency range 990  1260 MHz based on data taken at the Bleien
Observatory. For this survey, we expect to cover 70% of the full sky and
achieve a median signaltonoise ratio of approximately 5  6 in the cleanest
channels including systematic uncertainties. However, we also point out the
potential challenges of high RFI contamination and baseline removal when
examining the early data from the Bleien Observatory.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea galaxy
and quasar redshift survey. To trace the underlying dark matter distribution,
spectroscopic targets will be selected in four classes from imaging data. We
will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to
even higher redshift, DESI will target bright [O II] emission line galaxies up
to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying
dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the
Ly$\alpha$ forest absorption features in their spectra, which will be used to
trace the distribution of neutral hydrogen. When moonlight prevents efficient
observations of the faint targets of the baseline survey, DESI will conduct a
magnitudelimited Bright Galaxy Survey comprising approximately 10 million
galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy
and quasar redshifts will be obtained to measure the BAO feature and determine
the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
diameter field of view. The DESI collaboration will also deliver a
spectroscopic pipeline and data management system to reduce and archive all
data for eventual public use.

Recent observational progress has led to the establishment of the standard
$\Lambda$CDM model for cosmology. This development is based on different
cosmological probes that are usually combined through their likelihoods at the
latest stage in the analysis. We implement here an integrated scheme for
cosmological probes, which are combined in a common framework starting at the
map level. This treatment is necessary as the probes are generally derived from
overlapping maps and are thus not independent. It also allows for a thorough
test of the cosmological model and of systematics through the consistency of
different physical tracers. As a first application, we combine current
measurements of the Cosmic Microwave Background (CMB) from the Planck
satellite, and galaxy clustering and weak lensing from SDSS. We consider the
spherical harmonic power spectra of these probes including all six auto and
crosscorrelations along with the associated full Gaussian covariance matrix.
This provides an integrated treatment of different analyses usually performed
separately including CMB anisotropies, cosmic shear, galaxy clustering,
galaxygalaxy lensing and the Integrated SachsWolfe (ISW) effect with galaxy
and shear tracers. We derive constraints on $\Lambda$CDM parameters that are
compatible with existing constraints and highlight tensions between data sets,
which become apparent in this integrated treatment. We discuss how this
approach provides a complete and powerful integrated framework for probe
combination and how it can be extended to include other tracers in the context
of current and future wide field cosmological surveys.

We present extended modeling of the strong lens system RXJ11311231 with
archival data in two HST bands in combination with existing lineofsight
contribution and velocity dispersion estimates. Our focus is on source size and
its influence on timedelay cosmography. We therefore examine the impact of
masssheet degeneracy and especially the degeneracy pointed out by Schneider &
Sluse (2013) using the source reconstruction scale. We also extend on previous
work by further exploring the effects of priors on the kinematics of the lens
and the external convergence in the environment of the lensing system. Our
results coming from RXJ11311231 are given in a simple analytic form so that
they can be easily combined with constraints coming from other cosmological
probes. We find that the choice of priors on lens model parameters and source
size are subdominant for the statistical errors for $H_0$ measurements of this
systems. The choice of prior for the source is subdominant at present (2%
uncertainty on $H_0$) but may be relevant for future studies. More importantly,
we find that the priors on the kinematic anisotropy of the lens galaxy have a
significant impact on our cosmological inference. When incorporating all the
above modeling uncertainties, we find $H_0 = 86.6^{+6.8}_{6.9}$ km
s$^{1}$Mpc$^{1}$, when using kinematic priors similar to other studies. When
we use a different kinematic prior motivated by Barnab\`e et al. (2012) but
covering the same anisotropic range, we find $H_0 = 74.5^{+8.0}_{7.8}$ km
s$^{1}$Mpc$^{1}$. This means that the choice of kinematic modeling and priors
have a significant impact on cosmographic inferences. The way forward is either
to get better velocity dispersion measures which would down weight the impact
of the priors or to construct physically motivated priors for the velocity
dispersion model.