We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

Submillimetre-selected galaxies (SMGs) at high redshift ($z$ $\sim$ 2) are potential host galaxies of active galactic nuclei (AGN). If the local Universe is a good guide, $\sim$ 50$\%$ of the obscured AGN amongst the SMG population could be missed even in the deepest X-ray surveys. Radio observations are insensitive to obscuration; therefore, very long baseline interferometry (VLBI) can be used as a tool to identify AGN in obscured systems. A well-established upper limit to the brightness temperature of 10$^5$ K exists in star-forming systems, thus VLBI observations can distinguish AGN from star-forming systems via brightness temperature measurements. We present 1.6 GHz European VLBI Network (EVN) observations of four SMGs (with measured redshifts) to search for evidence of compact radio components associated with AGN cores. For two of the sources, e-MERLIN images are also presented. Out of the four SMGs observed, we detect one source, J123555.14, that has an integrated EVN flux density of 201 $\pm$ 15.2 $\mu$Jy, corresponding to a brightness temperature of 5.2 $\pm$ 0.7 $\times$ 10$^5$ K. We therefore identify that the radio emission from J123555.14 is associated with an AGN. We do not detect compact radio emission from a possible AGN in the remaining sources (J123600.10, J131225.73, and J163650.43). In the case of J131225.73, this is particularly surprising, and the data suggest that this may be an extended, jet-dominated AGN that is resolved by VLBI. Since the morphology of the faint radio source population is still largely unknown at these scales, it is possible that with a $\sim$ 10 mas resolution, VLBI misses (or resolves) many radio AGN extended on kiloparsec scales.

The mJy variable extragalactic radio sky is known to be broadly non-changing with approximately $3\%$ of persistent radio sources exhibiting variability which is largely AGN-related. In the faint (<mJy) flux density regime, it is widely accepted that the radio source population begins to change from AGN dominated to star-formation dominated, together with an emergent radio-quiet AGN component. Very little is known about the variable source component in this sub-mJy regime. In this paper, we provide the first insight into the $\mu$Jy variable sky by performing a careful analysis using deep VLA data in the well studied GOODS-N field. Using five epochs spread across 22 years, we investigate approximately 480 radio sources finding 10 which show signs of variability. We attribute this variability to the presence of an AGN in these systems. We confirm and extend the results of previous surveys, finding that variability in the faint radio sky is rather modest with only $\leq$2\% of sources exhibiting significant variability between any two epochs. We find that 70\% of variable sources show variability on timescales of a few days whilst on longer decadal time-scales, the fraction of variable sources decreases to $<1\%$. This suggests that the radio variability peaks on shorter timescales as suggested by other studies. We find that 80\% of variable sources have VLBI counterparts, and we use multi-wavelength data to infer that these may well be core-dominated FR-I sources as postulated by wide-field VLBI surveys and semi-empirical simulations.

The $e$MERLIN Galaxy Evolution Survey ($e$MERGE) is an ambitious, multi-tiered extragalactic radio continuum survey being carried out with $e$MERLIN and the VLA at 1.4GHz and 6GHz. Exploiting the unique combination of high sensitivity and high angular resolution provided by radio interferometry, these observations will provide a powerful, obscuration-independent tool for tracing intense star-formation and AGN activity in galaxies out to $z\sim5$. In our first data release (DR1) we present $e$MERGE Tier 1, a 15-arcmin pointing centred on the GOODS-N field, imaged at 1.4GHz with the VLA and $e$MERLIN at $\sim 0.28''$ resolution down to an rms sensitivity of $\sim 1.2\,\mu$Jy beam$^{-1}$. This unique radio survey -- unrivaled at 1.4GHz in its combination of depth, areal coverage and angular resolution in the pre-SKA era -- allows us to localise and separate extended star-forming regions, nuclear starbursts and compact AGN core/jet systems in galaxies over the past two-thirds of cosmic history, a crucial step in tracing the apparently simultaneous growths of the stellar populations and central black holes in massive galaxies. In these proceedings we highlight some early science results from $e$MERGE DR1, including some examples of the sub-arcsecond morphologies and cold dust properties of 1.4GHz-selected galaxies. $e$MERGE Tier 1 will eventually reach sub-$\mu$Jy beam$^{-1}$ sensitivity at $0.28''$ resolution over a 30-arcmin field, providing crucial benchmarks for deep extragalactic surveys which will be undertaken with SKA in the next decade.

(Abridged) Conventional radio surveys of deep fields ordinarily have arc-second scale resolutions often insufficient to reliably separate radio emission in distant galaxies originating from star-formation and AGN-related activity. Very long baseline interferometry (VLBI) can offer a solution by identifying only the most compact radio emitting regions in galaxies at cosmological distances where the high brightness temperatures (in excess of $10^5$ K) can only be reliably attributed to AGN activity. We present the first in a series of papers exploring the faint compact radio population using a new wide-field VLBI survey of the GOODS-N field. The unparalleled sensitivity of the European VLBI Network (EVN) will probe a luminosity range rarely seen in deep wide-field VLBI observations, thus providing insights into the role of AGN to radio luminosities of the order $10^{22}~\mathrm{W\,Hz^{-1}}$ across cosmic time. The newest VLBI techniques are used to completely cover an entire 7'.5 radius area to milliarcsecond resolutions, while bright radio sources ($S > 0.1$ mJy) are targeted up to 25 arcmin from the pointing centre. Multi-source self-calibration, and a primary beam model for the EVN array are used to correct for residual phase errors and primary beam attenuation respectively. This paper presents the largest catalogue of VLBI detected sources in GOODS-N comprising of 31 compact radio sources across a redshift range of 0.11-3.44, almost three times more than previous VLBI surveys in this field. We provide a machine-readable catalogue and introduce the radio properties of the detected sources using complementary data from the e-MERLIN Galaxy Evolution survey (eMERGE).

We present new observations of the GOODS-N field obtained at 5.5 GHz with the Karl G. Jansky Very Large Array (VLA). The central region of the field was imaged to a median r.m.s. of 3 microJy/beam with a resolution of 0.5 arcsec. From a 14-arcmin diameter region we extracted a sample of 94 radio sources with signal-to-noise ratio greater than 5. Near-IR identifications are available for about 88 percent of the radio sources. We used different multi-band diagnostics to separate active galactic nuclei (AGN), both radiatively efficient and inefficient, from star-forming galaxies. From our analysis, we find that about 80 percent of our radio-selected sample is AGN-dominated, with the fraction raising to 92 percent when considering only the radio sources with redshift >1.5. This large fraction of AGN-dominated radio sources at very low flux densities (the median flux density at 5.5 GHz is 42 microJy), where star-forming galaxies are expected to dominate, is somewhat surprising and at odds with other results. Our interpretation is that both the frequency and angular resolution of our radio observations strongly select against radio sources whose brightness distribution is diffuse on scale of several kpc. Indeed, we find that the median angular sizes of the AGN-dominated sources is around 0.2-0.3 arcsec against 0.8 arcsec for star-forming galaxies. This highlights the key role that high frequency radio observations can play in pinpointing AGN-driven radio emission at microJy levels. This work is part of the eMERGE legacy project.

Context. Very Long Baseline Interferometry (VLBI) data are extremely sensitive to the phase stability of the VLBI array. This is especially important when we reach {\mu}Jy r.m.s. sensitivities. Calibration using standard phase referencing techniques is often used to improve the phase stability of VLBI data but the results are often not optimal. This is evident in blank fields that do not have in-beam calibrators. Aims. We present a calibration algorithm termed Multi-Source Self-Calibration (MSSC) which can be used after standard phase referencing on wide-field VLBI observations. This is tested on a 1.6 GHz wide-field VLBI data set of the Hubble Deep Field-North and the Hubble Flanking Fields. Methods. MSSC uses multiple target sources detected in the field via standard phase referencing techniques and modifies the visibili- ties so that each data set approximates to a point source. These are combined to increase the signal to noise and permit self-calibration. In principle, this should allow residual phase changes caused by the troposphere and ionosphere to be corrected. By means of faceting, the technique can also be used for direction dependent calibration. Results. Phase corrections, derived using MSSC, were applied to a wide-field VLBI data set of the HDF-N comprising of 699 phase centres. MSSC was found to perform considerably better than standard phase referencing and single source self-calibration. All detected sources exhibited dramatic improvements in dynamic range. Using MSSC, one source reached the detection threshold taking the total detected sources to twenty. 60% of these sources can now be imaged with uniform weighting compared to just 45% with standard phase referencing. The Parseltongue code which implements MSSC has been released and made publicly available to the astronomical community (https://github.com/jradcliffe5/multi_self_cal).