The Nancay Decameter Array (NDA) routinely observes low frequency (10-100 MHz) radio emissions of Jupiter and the Sun since 4 decades. The NDA observations, acquired with a variety of receivers with increasing performances, were the basis for numerous studies of jovian and solar radio emissions and now form a unique long-term database spanning >3 solar cycles and jovian revolutions. In addition, the NDA historically brought a fruitful support to space-based radio observatories of the heliosphere, to multi-wavelength analyses of solar activity and contributes to the development of space weather services. After having summarized the NDA characteristics, this article presents latest instrumental and database developments, some recent scientific results and perspectives for the next decade.

The high latitude radio emissions produced by the Cyclotron Maser Instability (CMI) in Jupiter's magnetosphere extend from a few kHz to 40 MHz. Part of the decametric emissions is of auroral origin, and part is driven by the moons Io, Europa and Ganymede. After summarizing the method used to identify Jupiter-satellite radio emissions, which consists in comparing space- and ground-based radio observations to ExPRES simulations of CMI-driven emissions in the time-frequency plane, we present a parametric study of the free parameters required by the ExPRES code (electron distribution function and resonant energy, magnetic field model, lead angle, and altitude of the ionospheric cut-off) in order to assess the accuracy of our simulations in the Io-Jupiter case. We find that Io-DAM arcs are fairly modeled by loss-cone driven CMI with electrons of 1-10 keV energy, using the ISaAC, VIPAL or VIP4 magnetic field model and a simple sinusoidal lead angle model. The altitude of the ionospheric cut-off has a marginal impact on the simulations. We discuss the impact of our results on the identification of Europa-DAM and Ganymede-DAM emissions.

With their wide fields of view and often relatively long coverage of any position in the sky in imaging survey mode, modern radio telescopes provide a data stream that is naturally suited to searching for rare transients. However, Radio Frequency Interference (RFI) can show up in the data stream in similar ways to such transients, and thus the normal pre-treatment of filtering RFI (flagging) may also remove astrophysical transients from the data stream before imaging. In this paper we investigate how standard flagging affects the detectability of such transients by examining the case of transient detection in an observing mode used for Low Frequency Array (LOFAR; \citep{LOFAR}) surveys. We quantify the fluence range of transients that would be detected, and the reduction of their SNR due to partial flagging. We find that transients with a duration close to the integration sampling time, as well as bright transients with durations on the order of tens of seconds, are completely flagged. For longer transients on the order of several tens of seconds to minutes, the flagging effects are not as severe, although part of the signal is lost. For these transients, we present a modified flagging strategy which mitigates the effect of flagging on transient signals. We also present a script which uses the differences between the two strategies, and known differences between transient RFI and astrophysical transients, to notify the observer when a potential transient is in the data stream.

Type II radio bursts are evidence of shocks in the solar atmosphere and inner heliosphere that emit radio waves ranging from sub-meter to kilometer lengths. These shocks may be associated with CMEs and reach speeds higher than the local magnetosonic speed. Radio imaging of decameter wavelengths (20-90 MHz) is now possible with LOFAR, opening a new radio window in which to study coronal shocks that leave the inner solar corona and enter the interplanetary medium and to understand their association with CMEs. To this end, we study a coronal shock associated with a CME and type II radio burst to determine the locations at which the radio emission is generated, and we investigate the origin of the band-splitting phenomenon.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

We present LOFAR high-band data over the frequency range 115-189 MHz for the X-ray binary SS 433, obtained in an observing campaign from 2013 February - 2014 May. Our results include a deep, wide-field map, allowing a detailed view of the surrounding supernova remnant W 50 at low radio frequencies, as well as a light curve for SS 433 determined from shorter monitoring runs. The complex morphology of W 50 is in excellent agreement with previously published higher-frequency maps; we find additional evidence for a spectral turnover in the eastern wing, potentially due to foreground free-free absorption. Furthermore, SS 433 is tentatively variable at 150 MHz, with both a debiased modulation index of 11 per cent and a $\chi^2$ probability of a flat light curve of $8.2 \times 10^{-3}$. By comparing the LOFAR flux densities with contemporaneous observations carried out at 4800 MHz with the RATAN-600 telescope, we suggest that an observed $\sim$0.5-1 Jy rise in the 150-MHz flux density may correspond to sustained flaring activity over a period of approximately six months at 4800 MHz. However, the increase is too large to be explained with a standard synchrotron bubble model. We also detect a wealth of structure along the nearby Galactic plane, including the most complete detection to date of the radio shell of the candidate supernova remnant G 38.7-1.4. This further demonstrates the potential of supernova remnant studies with the current generation of low-frequency radio telescopes.

Analyzing a database of 26 years of observations of Jupiter from the Nan\c{c}ay Decameter Array, we study the occurrence of Io-independent emissions as a function of the orbital phase of the other Galilean satellites and Amalthea. We identify unambiguously the emissions induced by Ganymede and characterize their intervals of occurrence in CML and Ganymede phase and longitude. We also find hints of emissions induced by Europa and, surprisingly, by Amalthea. The signature of Callisto-induced emissions is more tenuous.

Context. The Sun is an active star that produces large-scale energetic events such as solar flares and coronal mass ejections and numerous smaller-scale events such as solar jets. These events are often associated with accelerated particles that can cause emission at radio wavelengths. The reconfiguration of the solar magnetic field in the corona is believed to be the cause of the majority of solar energetic events and accelerated particles. Aims. Here, we investigate a bright J-burst that was associated with a solar jet and the possible emission mechanism causing these two phenomena. Methods. We used data from the Solar Dynamics Observatory (SDO) to observe a solar jet, and radio data from the Low Frequency Array (LOFAR) and the Nan\c{c}ay Radioheliograph (NRH) to observe a J-burst over a broad frequency range (33-173 MHz) on 9 July 2013 at ~11:06 UT. Results. The J-burst showed fundamental and harmonic components and it was associated with a solar jet observed at extreme ultraviolet wavelengths with SDO. The solar jet occurred at a time and location coincident with the radio burst, in the northern hemisphere, and not inside a group of complex active regions in the southern hemisphere. The jet occurred in the negative polarity region of an area of bipolar plage. Newly emerged positive flux in this region appeared to be the trigger of the jet. Conclusions. Magnetic reconnection between the overlying coronal field lines and the newly emerged positive field lines is most likely the cause of the solar jet. Radio imaging provides a clear association between the jet and the J-burst which shows the path of the accelerated electrons.

HD189733 is an active K dwarf that is, with its transiting hot Jupiter, among the most studied exoplanetary systems. In this first paper of the Multiwavelength Observations of an eVaporating Exoplanet and its Star (MOVES) program, we present a 2-year monitoring of the large-scale magnetic field of HD189733. The magnetic maps are reconstructed for five epochs of observations, namely June-July 2013, August 2013, September 2013, September 2014, and July 2015, using Zeeman-Doppler Imaging. We show that the field evolves along the five epochs, with mean values of the total magnetic field of 36, 41, 42, 32 and 37 G, respectively. All epochs show a toroidally-dominated field. Using previously published data of Moutou et al. 2007 and Fares et al. 2010, we are able to study the evolution of the magnetic field over 9 years, one of the longest monitoring campaign for a given star. While the field evolved during the observed epochs, no polarity switch of the poles was observed. We calculate the stellar magnetic field value at the position of the planet using the Potential Field Source Surface extrapolation technique. We show that the planetary magnetic environment is not homogeneous over the orbit, and that it varies between observing epochs, due to the evolution of the stellar magnetic field. This result underlines the importance of contemporaneous multi-wavelength observations to characterise exoplanetary systems. Our reconstructed maps are a crucial input for the interpretation and modelling of our MOVES multi-wavelength observations.

This paper describes digital radio astronomical receivers used for decameter and meter wavelength observations. This paper describes digital radio astronomical receivers used for decameter and meter wavelength observations. Since 1998, digital receivers performing on-the-fly dynamic spectrum calculations or waveform data recording without data loss have been used at the UTR-2 radio telescope, the URAN VLBI system, and the GURT new generation radio telescope. Here we detail these receivers developed for operation in the strong interference environment that prevails in the decameter wavelength range. Data collected with these receivers allowed us to discover numerous radio astronomical objects and phenomena at low frequencies, a summary of which is also presented.

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

Context: Active stars possess magnetized wind that has a direct impact on planets that can lead to radio emission. Mercury is a good test case to study the effect of the solar wind and interplanetary magnetic field on radio emission driven in the planet magnetosphere. Such studies could be used as proxies to characterize the magnetic field topology and intensity of exoplanets. Aims: The aim of this study is to quantify the radio emission in the Hermean magnetosphere. Methods: We use the MHD code PLUTO in spherical coordinates with an axisymmetric multipolar expansion for the Hermean magnetic field, to analyze the effect of the interplanetary magnetic field (IMF) orientation and intensity, as well as the hydrodynamic parameters of the solar wind (velocity, density and temperature), on the net power dissipated on the Hermean day and night side. We apply the formalism derived by Zarka [2001, 2007] to infer the radio emission level from the net dissipated power. We perform a set of simulations with different hydrodynamic parameters of the solar wind, IMF orientations and intensities, that allow us to calculate the dissipated power distribution and infer the existence of radio emission hot spots on the planet day side, and to calculate the integrated radio emission of the Hermean magnetosphere. Results: The obtained radio emission distribution of dissipated power is determined by the IMF orientation (associated with the reconnection regions in the magnetosphere), although the radio emission strength is dependent on the IMF intensity and solar wind hydro parameters. The calculated total radio emission level is in agreement with the one estimated in [Zarka, 2001], between $5 \times 10^{5}$ and $2 \times 10^{6}$ W.

The habitability of an exoplanet depends on many factors. One such factor is the impact of stellar eruptive events on nearby exoplanets. Currently this is poorly constrained due to heavy reliance on solar scaling relationships and a lack of experimental evidence. Potential impacts of Coronal Mass Ejections (CMEs), which are a large eruption of magnetic field and plasma from a star, are space weather and atmospheric stripping. A method for observing CMEs as they travel though the stellar atmosphere is the type II radio burst, and the new LOw Frequency ARray (LOFAR) provides a means for detection. We report on 15 hours of observation of YZ Canis Minoris (YZ CMi), a nearby M dwarf flare star, taken in LOFAR's beam-formed observation mode for the purposes of measuring transient frequency-dependent low frequency radio emission. The observations utilized Low-Band Antenna (10-90 MHz) or High-Band Antenna (110-190 MHz) for five three-hour observation periods. In this data set, there were no confirmed type II events in this frequency range. We explore the range of parameter space for type II bursts constrained by our observations Assuming the rate of shocks is a lower limit to the rate at which CMEs occur, no detections in a total of 15 hours of observation places a limit of $\nu_{type II} < 0.0667$ shocks/hr $ \leq \nu_{CME}$ for YZ CMi due to the stochastic nature of the events and limits of observational sensitivity. We propose a methodology to interpret jointly observed flares and CMEs which will provide greater constraints to CMEs and test the applicability of solar scaling relations.

We present first results from a LOFAR census of non-recycled pulsars. The census includes almost all such pulsars known (194 sources) at declinations ${\rm Dec}> 8^\circ$ and Galactic latitudes $|{\rm Gb}| > 3^\circ$, regardless of their expected flux densities and scattering times. Each pulsar was observed for $\geq 20$ minutes in the contiguous frequency range of 110--188 MHz. Full-Stokes data were recorded. We present the dispersion measures, flux densities, and calibrated total intensity profiles for the 158 pulsars detected in the sample. The median uncertainty in census dispersion measures ($1.5 \times 10^{-3}$ pc cm$^{-3}$) is ten times smaller, on average, than in the ATNF pulsar catalogue. We combined census flux densities with those in the literature and fitted the resulting broadband spectra with single or broken power-law functions. For 48 census pulsars such fits are being published for the first time. Typically, the choice between single and broken power-laws, as well as the location of the spectral break, were highly influenced by the spectral coverage of the available flux density measurements. In particular, the inclusion of measurements below 100 MHz appears essential for investigating the low-frequency turnover in the spectra for most of the census pulsars. For several pulsars, we compared the spectral indices from different works and found the typical spread of values to be within 0.5--1.5, suggesting a prevailing underestimation of spectral index errors in the literature. The census observations yielded some unexpected individual source results, as we describe in the paper. Lastly, we will provide this unique sample of wide-band, low-frequency pulse profiles via the European Pulsar Network Database.

The eclipses of certain types of binary millisecond pulsars (i.e. `black widows' and `redbacks') are often studied using high-time-resolution, `beamformed' radio observations. However, they may also be detected in images generated from interferometric data. As part of a larger imaging project to characterize the variable and transient sky at radio frequencies <200 MHz, we have blindly detected the redback system PSR J2215+5135 as a variable source of interest with the Low-Frequency Array (LOFAR). Using observations with cadences of 2 weeks - 6 months, we find preliminary evidence that the eclipse duration is frequency dependent ($\propto \nu^{-0.4}$), such that the pulsar is eclipsed for longer at lower frequencies, in broad agreement with beamformed studies of other similar sources. Furthermore, the detection of the eclipses in imaging data suggests an eclipsing medium that absorbs the pulsed emission, rather than scattering it. Our study is also a demonstration of the prospects of finding pulsars in wide-field imaging surveys with the current generation of low-frequency radio telescopes.

The magnetospheric cusps are important sites of the coupling of a magnetosphere with the solar wind. The combination of both ground- and space-based observations at Earth have enabled considerable progress to be made in understanding the terrestrial cusp and its role in the coupling of the magnetosphere to the solar wind via the polar magnetosphere. Voyager 2 fully explored Neptune's cusp in 1989 but highly inclined orbits of the Cassini spacecraft at Saturn present the most recent opportunity to repeatedly studying the polar magnetosphere of a rapidly rotating planet. In this paper we discuss observations made by Cassini during two passes through Saturn's southern polar magnetosphere. Our main findings are that i) Cassini directly encounters the southern polar cusp with evidence for the entry of magnetosheath plasma into the cusp via magnetopause reconnection, ii) magnetopause reconnection and entry of plasma into the cusp can occur over a range of solar wind conditions, and iii) double cusp morphologies are consistent with the position of the cusp oscillating in phase with Saturn's global magnetospheric periodicities.

We report on the results of a search for radio transients between 115 and 190\,MHz with the LOw-Frequency ARray (LOFAR). Four fields have been monitored with cadences between 15 minutes and several months. A total of 151 images were obtained, giving a total survey area of 2275 deg$^2$. We analysed our data using standard LOFAR tools and searched for radio transients using the LOFAR Transients Pipeline (TraP). No credible radio transient candidate has been detected; however, we are able to set upper limits on the surface density of radio transient sources at low radio frequencies. We also show that low-frequency radio surveys are more sensitive to steep-spectrum coherent transient sources than GHz radio surveys. We used two new statistical methods to determine the upper limits on the transient surface density. One is free of assumptions on the flux distribution of the sources, while the other assumes a power-law distribution in flux and sets more stringent constraints on the transient surface density. Both of these methods provide better constraints than the approach used in previous works. The best value for the upper limit we can set for the transient surface density, using the method assuming a power-law flux distribution, is 1.3$\cdot$10$^{-3}$ deg$^{-2}$ for transients brighter than 0.3 Jy with a time-scale of 15 min, at a frequency of 150 MHz. We also calculated for the first time upper limits for the transient surface density for transients of different time-scales. We find that the results can differ by orders of magnitude from previously reported, simplified estimates.

Air showers induced by cosmic rays create nanosecond pulses detectable at radio frequencies. These pulses have been measured successfully in the past few years at the LOw Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete understanding of this phenomenon and the underlying physical processes, an absolute calibration of the detecting antenna system is needed. We present three approaches that were used to check and improve the antenna model of LOFAR and to provide an absolute calibration of the whole system for air shower measurements. Two methods are based on calibrated reference sources and one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short data-sets. An accuracy of 19% in amplitude is reached. The absolute calibration is also compared to predictions from air shower simulations. These results are used to set an absolute energy scale for air shower measurements and can be used as a basis for an absolute scale for the measurement of astronomical transients with LOFAR.

LS I +61 303 is a gamma-ray binary that exhibits an outburst at GHz frequencies each orbital cycle of $\approx$ 26.5 d and a superorbital modulation with a period of $\approx$ 4.6 yr. We have performed a detailed study of the low-frequency radio emission of LS I +61 303 by analysing all the archival GMRT data at 150, 235 and 610 MHz, and conducting regular LOFAR observations within the Radio Sky Monitor (RSM) at 150 MHz. We have detected the source for the first time at 150 MHz, which is also the first detection of a gamma-ray binary at such a low frequency. We have obtained the light-curves of the source at 150, 235 and 610 MHz, all of them showing orbital modulation. The light-curves at 235 and 610 MHz also show the existence of superorbital variability. A comparison with contemporaneous 15-GHz data shows remarkable differences with these light-curves. At 15 GHz we see clear outbursts, whereas at low frequencies we see variability with wide maxima. The light-curve at 235 MHz seems to be anticorrelated with the one at 610 MHz, implying a shift of $\sim$ 0.5 orbital phases in the maxima. We model the shifts between the maxima at different frequencies as due to changes in the physical parameters of the emitting region assuming either free-free absorption or synchrotron self-absorption, obtaining expansion velocities for this region close to the stellar wind velocity with both mechanisms.

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

We report the detection of 48 millisecond pulsars (MSPs) out of 75 observed thus far using the LOFAR in the frequency range 110-188 MHz. We have also detected three MSPs out of nine observed in the frequency range 38-77 MHz. This is the largest sample of MSPs ever observed at these low frequencies, and half of the detected MSPs were observed for the first time at frequencies below 200 MHz. We present the average pulse profiles of the detected MSPs, their effective pulse widths, and flux densities and compare these with higher observing frequencies. The flux-calibrated, multifrequency LOFAR pulse profiles are publicly available via the EPN Database of Pulsar Profiles. We also present average values of dispersion measures (DM) and discuss DM and profile variations. About 35% of the MSPs show strong narrow profiles, another 25% exhibit scattered profiles, and the rest are only weakly detected. A qualitative comparison of LOFAR profiles with those at higher radio frequencies shows constant separation between profile components. Similarly, the profile widths are consistent with those observed at higher frequencies, unless scattering dominates at the lowest frequencies. This is very different from what is observed for normal pulsars and suggests a compact emission region in the MSP magnetosphere. The amplitude ratio of the profile components, on the other hand, can dramatically change towards low frequencies, often with the trailing component becoming dominant. As previously demonstrated this can be caused by aberration and retardation. This data set enables high-precision studies of pulse profile evolution with frequency, dispersion, Faraday rotation, and scattering in the interstellar medium. Characterising and correcting these systematic effects may improve pulsar-timing precision at higher observing frequencies, where pulsar timing array projects aim to directly detect gravitational waves.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

Context. The Sun is an active source of radio emission that is often associated with energetic phenomena ranging from nanoflares to coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), numerous millisecond duration radio bursts have been reported, such as radio spikes or solar S bursts (where S stands for short). To date, these have neither been studied extensively nor imaged because of the instrumental limitations of previous radio telescopes. Aims. Here, Low Frequency Array (LOFAR) observations were used to study the spectral and spatial characteristics of a multitude of S bursts, as well as their origin and possible emission mechanisms. Methods. We used 170 simultaneous tied-array beams for spectroscopy and imaging of S bursts. Since S bursts have short timescales and fine frequency structures, high cadence (~50 ms) tied-array images were used instead of standard interferometric imaging, that is currently limited to one image per second. Results. On 9 July 2013, over 3000 S bursts were observed over a time period of ~8 hours. S bursts were found to appear as groups of short-lived (<1 s) and narrow-bandwidth (~2.5 MHz) features, the majority drifting at ~3.5 MHz/s and a wide range of circular polarisation degrees (2-8 times more polarised than the accompanying Type III bursts). Extrapolation of the photospheric magnetic field using the potential field source surface (PFSS) model suggests that S bursts are associated with a trans-equatorial loop system that connects an active region in the southern hemisphere to a bipolar region of plage in the northern hemisphere. Conclusions. We have identified polarised, short-lived solar radio bursts that have never been imaged before. They are observed at a height and frequency range where plasma emission is the dominant emission mechanism, however they possess some of the characteristics of electron-cyclotron maser emission.