We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $\mu$Jy/beam rms noise. Its J2000 center position $\alpha=04^h 13^m 26.4^s$, $\delta=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $\mu$Jy with counts of individual DEEP2 sources between $10$ $\mu$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $\mu$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $\mu$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.

Radio interferometric gain calibration can be biased by incomplete sky models and radio frequency interference, resulting in calibration artefacts that can restrict the dynamic range of the resulting images. It has been suggested that calibration algorithms employing heavy-tailed likelihood functions are less susceptible to this due to their robustness against outliers in the data. We present an algorithm based on a Student's t-distribution which leverages the framework of complex optimisation and Wirtinger calculus for efficient and robust interferometric gain calibration. We integrate this algorithm as an option in the newly released calibration software package, CubiCal. We demonstrate that the algorithm can mitigate some of the biases introduced by incomplete sky models and radio frequency interference by applying it to both simulated and real data. Our results show significant improvements compared to a conventional least-squares solver which assumes a Gaussian likelihood function. Furthermore, we provide some insight into why the algorithm outperforms the conventional solver, and discuss specific scenarios (for both direction-independent and direction-dependent self-calibration) where this is expected to be the case.

We use a direct observational approach to investigate the possibility of testing the Copernican principle with data from upcoming radio surveys. In particular we illustrate the importance of measuring derivatives transverse to the past light-cone when prior knowledge of the value of the cosmological constant is not available.

The construction of the KAT-7 array in the Karoo region of the Northern Cape in South Africa was intended primarily as an engineering prototype for technologies and techniques applicable to the MeerKAT telescope. This paper looks at the main engineering and scien- tific highlights from this effort, and discusses their applicability to both MeerKAT and other next-generation radio telescopes. In particular we found that the composite dish surface works well, but it becomes complicated to fabricate for a dish lacking circular symmetry; the Stir- ling cycle cryogenic system with ion pump to achieve vacuum works but demands much higher maintenance than an equivalent Gifford-McMahon cycle system; the ROACH (Recon- figurable Open Architecture Computing Hardware)-based correlator with SPEAD (Stream- ing Protocol for Exchanging Astronomical Data) protocol data transfer works very well and KATCP (Karoo Array Telescope Control Protocol) control protocol has proven very flexible and convenient. KAT-7 has also been used for scientific observations where it has a niche in mapping low surface-brightness continuum sources, some extended HI halos and OH masers in star-forming regions. It can also be used to monitor continuum source variability, observe pulsars, and make VLBI observations

We give an outline of an algorithm designed to reconstruct the background cosmological metric within the class of spherically symmetric dust universes that may include a cosmological constant. Luminosity and age data are used to derive constraints on the geometry of the universe up to a redshift of $z = 1.75$. It is shown that simple radially inhomogeneous void models that are sometimes used as alternative explanations for the apparent acceleration of the late time Universe cannot be ruled out by these data alone.

We present a new algorithm that can reconstruct the full distributions of metric components within the class of spherically symmetric dust universes that may include a cosmological constant. The algorithm is capable of confronting this class of solutions with arbitrary data and opens a new observational window to determine the value of the cosmological constant. In this work we use luminosity and age data to constrain the geometry of the universe up to a redshift of $z = 1.75$. We show that, although current data are perfectly compatible with homogeneous models of the universe, simple radially inhomogeneous void models that are sometimes used as alternative explanations for the apparent acceleration of the late time universe cannot yet be ruled out. In doing so we reconstruct the density of cold dark matter out to $z = 1.75$ and derive constraints on the metric components when the universe was 10.5 Gyr old within a comoving volume of approximately 1 Gpc$^{3}$.

We investigate the possibility of using Gaussian process regression to smooth data on the current past null-cone for use as the input to a relativistic integration scheme. The algorithm we present is designed to reconstruct the metric of spacetime within the class of spherically symmetric dust universes, with or without a cosmological constant. Assuming that gravity is well described by General Relativity, we demonstrate how the algorithm can be employed to test the Copernican principle based on currently available observations. It is shown that currently available data is not sufficient for a conclusive result. The intrinsic noise present in realistic data presents a challenge for our smoothing algorithm and we discuss some of its limitations as well as possible extensions to it. We conclude by demonstrating how a direct determination of the cosmological constant is possible using redshift drift data.