The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120–168 MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45°00′00″ to 57°00′00″) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325 694 sources are detected with a signal of at least five times the noise, and the source density is a factor of ∼10 higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S144 MHz = 71 μJy beam−1 and the point-source completeness is 90% at an integrated flux density of 0.45 mJy. The resolution of the images is 6″ and the positional accuracy is within 0.2″. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further ∼20 articles that highlight the scientific potential of LoTSS.

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire northern sky. Each of the 3170 pointings will be observed for 8 h, which, at most declinations, is sufficient to produce ∼5 resolution images with a sensitivity of ∼100 µJy/beam and accomplish the main scientific aims of the survey, which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Owing to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44 000 sources are detected in the images that have a resolution of 25 , typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45• 00 00 to 57• 00 00 ).

Context. Radio relics are diffuse extended synchrotron sources that originate from shock fronts induced by galaxy cluster mergers. The actual particle acceleration mechanism at the shock fronts is still under debate. The galaxy cluster 1RXS J0603.3+4214 hosts one of the most intriguing examples of radio relics, known as the Toothbrush. Aims. To understand the mechanism(s) that accelerate relativistic particles in the intracluster medium (ICM), we investigate the spectral properties of large scale diffuse extended sources in the merging galaxy cluster 1RXS J0603.3+4214. Methods. We present new wideband radio continuum observations made with uGMRT and VLA. Our new observations in combination with previously published data, allowed us to carry out a detailed high spatial resolution spectral and curvature analysis of the known diffuse radio emission sources, over a broad range of frequencies.Results. The integrated spectrum of the Toothbrush follows closely a power law over close to 2 decades in frequency, with a spectral index of −1.16 ± 0.02. We do not find any evidence of spectral steepening below 8 GHz. The subregions of the main Toothbrush also exhibit near-perfect power laws, implying a very regular combination of shock properties across the shock front. Recent numerical simulations show an intriguing similar spectral index, suggesting that the radio spectrum is dominated by the average over the inhomogeneities within the shock, with most of the emission coming from the tail of the Mach number distribution. In contrast to the Toothbrush, the spectrum of the fainter relics show a high frequency steepening. Moreover, also the integrated spectrum of the halo follows a power law from 150 MHz to 3 GHz with a spectral index of −1.16 ± 0.04. We do not find any evidence for spectral curvature, not even in subareas of the halo. This suggest a homogeneous acceleration throughout the cluster volume. Between the "brush" region of the Toothbrush and the halo, the color-color analysis revealed emission that was consistent with an overlap between the two different spectral regions. Conclusions. None of the relic structures, the Toothbrush as a whole or its subregions or the other two fainter relics, show spectral shapes consistent with a single injection of relativistic electrons, such as at a shock, followed by synchrotron aging in a relatively homogeneous environment. Inhomogeneities in some combination of Mach number, magnetic field strengths and projection effects dominate the observed spectral shapes.

Context. New generation low-frequency telescopes are exploring a new parameter space in terms of depth and resolution. The data taken with these interferometers, for example with the LOw Frequency ARray (LOFAR), are often calibrated in a low signal-to-noise ratio regime and the removal of critical systematic effects is challenging. The process requires an understanding of their origin and properties. Aims. In this paper we describe the major systematic effects inherent to next generation low-frequency telescopes, such as LOFAR. With this knowledge, we introduce a data processing pipeline that is able to isolate and correct these systematic effects. The pipeline will be used to calibrate calibrator observations as the first step of a full data reduction process. Methods. We processed two LOFAR observations of the calibrator 3C 196: the first using the Low Band Antenna (LBA) system at 42-66 MHz and the second using the High Band Antenna (HBA) system at 115-189 MHz. Results. We were able to isolate and correct for the effects of clock drift, polarisation misalignment, ionospheric delay, Faraday rotation, ionospheric scintillation, beam shape, and bandpass. The designed calibration strategy produced the deepest image to date at 54 MHz. The image has been used to confirm that the spectral energy distribution of the average radio source population tends to flatten at low frequencies.Conclusions. We prove that LOFAR systematic effects can be described by a relatively small number of parameters. Furthermore, the identification of these parameters is fundamental to reducing the degrees of freedom when the calibration is carried out on fields that are not dominated by a strong calibrator.

In this data release from the ongoing LOw-Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) we present 120-168 MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44 • 30 and 1h00m +28 • 00 and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451 hrs (7.6 PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6 resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144 MHz have: a median rms sensitivity of 83 µJy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2 ; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8 mJy/beam. By creating three 16 MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of > ±0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20 resolution 120-168 MHz continuum images have a median rms sensitivity of 95 µJy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 × 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8 mJy/beam at 4 and 2.2 mJy/beam at 20 ; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data to facilitate the thorough scientific exploitation of this unique dataset.