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 $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.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

We present 5 GHz e-MERLIN observations of the GOODS-N region at sub-arcsec resolution (0.2--0.5 arcsec). These data form part of the early commissioning observations for the e-MERLIN interferometer and a pilot for the e-MERLIN legacy program eMERGE. A total of 17 sources were detected with S/N>3. These observations provide unique information on the radio source morphology at sub-arcsec scales. For twelve of these sources, deeper 1.4 GHz MERLIN+VLA observations at the same spatial resolution are available, allowing radio spectral indices to be derived for ten sources on sub-arcsec angular scales. Via analysis of the spectral indices and radio morphologies, these sources have been identified as AGN cores in moderate-to-high redshift (1<z<4) galaxies. These results have provided AGN (or AGN candidate) classification for six previously unclassified sources and confirmed the AGN nature of the rest of the sample. Ultimately the eMERGE project will image the GOODS-N region at 1.4 and 5 GHz with higher resolution (about 50 mas at 5 GHz) and down to sub-microJy sensitivities. The unique combination of sensitivity and spatial resolution will be exploited to study star formation and AGN activity in distant galaxies.