Context. Solar radio bursts originate mainly from high energy electrons accelerated in solar eruptions like solar flares, jets, and coronal mass ejections. A subcategory of solar radio bursts with short time duration may be used as a proxy to understand the wave generation and propagation within the corona. Aims. Complete case studies of the source size, position and kinematics of short term bursts are very rare due to instrumental limitations. A comprehensive multi-frequency spectroscopic and imaging study was carried out of a clear example of a solar type IIIb-III pair. Methods. In this work, the source of the radio burst was imaged with the interferometric mode, using the remote baselines of the LOw Frequency ARray (LOFAR). A detailed analysis of the fine structures in the spectrum and of the radio source motion with imaging was conducted. Results. The study shows how the fundamental and harmonic components have a significantly different source motion. The apparent source of the fundamental emission at 26.56MHz displaces away from the solar disk center at about 4 times the speed of light, while the apparent source of the harmonic emission at the same frequency shows a speed of < 0.02c. The source size of the harmonic emission, observed in this case, is smaller than that in previous studies, indicating the importance of the use of the remote baselines.

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported.

The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (<100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15 arcsec with a high dynamic range and good image fidelity. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna (LBA) system. We reduced the datasets and obtained an image for each A-team source. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further.

Solar S-bursts are short duration ($<1$ s at decameter wavelengths) radio bursts that have been observed during periods of moderate solar activity, where S stands for short. The frequency drift of S-bursts can reflect the density variation and the motion state of the electron beams. In this work, we investigate the frequency drift and the fine structure of the S-bursts with the LOw Frequency ARray (LOFAR). We find that the average frequency drift rate of the S-bursts within 110-180MHz could be described by $df/dt=-0.0077f^{1.59}$. With the high time and frequency resolution of LOFAR, we can resolve the fine structures of the observed solar S-bursts. A fine drift variation pattern was found in the structure of S-bursts (referred to as solar Sb-bursts in this paper) during the type-III storm on 2019 April 13, in the frequency band of 120-240 MHz. The Sb-bursts have a quasi-periodic segmented pattern, and the relative flux intensity tends to be large when the frequency drift rate is relatively large. This kind of structure exists in about 20\% of the solar S-burst events within the observed frequency range. We propose that the fine structure is due to the density fluctuations of the background coronal density. We performed a simulation based on this theory which can reproduce the shape and relative flux intensity of the Sb-bursts. This work shows that the fine structure of solar radio bursts can be used to diagnose the coronal plasma.

The Sun is an active star that can launch large eruptions of magnetised plasma into the heliosphere, called coronal mass ejections (CMEs). These ejections can drive shocks that accelerate particles to high energies, often resulting in radio emission at low frequencies (<200 MHz). To date, the relationship between the expansion of CMEs, shocks and particle acceleration is not well understood, partly due to the lack of radio imaging at low frequencies during the onset of shock-producing CMEs. Here, we report multi-instrument radio, white-light and ultraviolet imaging of the second largest flare in Solar Cycle 24 (2008-present) and its associated fast CME (3038+/-288 km/s). We identify the location of a multitude of radio shock signatures, called herringbones, and find evidence for shock accelerated electron beams at multiple locations along the expanding CME. These observations support theories of non-uniform, rippled shock fronts driven by an expanding CME in the solar corona.

We present an overview of the LOFAR Tied-Array All-Sky Survey (LOTAAS) for radio pulsars and fast transients. The survey uses the high-band antennas of the LOFAR Superterp, the dense inner part of the LOFAR core, to survey the northern sky (dec > 0 deg) at a central observing frequency of 135 MHz. A total of 219 tied-array beams (coherent summation of station signals, covering 12 square degrees), as well as three incoherent beams (covering 67 square degrees) are formed in each survey pointing. For each ofthe 222 beams, total intensity is recorded at 491.52 us time resolution. Each observation integrates for 1 hr and covers 2592 channels from 119 to 151 MHz. This instrumental setup allows LOTAAS to reach a detection threshold of 1 to 5 mJy for periodic emission. Thus far, the LOTAAS survey has resulted in the discovery of 73 radio pulsars. Among these are two mildly recycled binary millisecond pulsars (P = 13 and 33 ms), as well as the slowest-spinning radio pulsar currently known (P = 23.5 s). The survey has thus far detected 311 known pulsars, with spin periods ranging from 4 ms to 5.0 s and dispersion measures from 3.0 to 217 pc/cc. Known pulsars are detected at flux densities consistent with literature values. We find that the LOTAAS pulsar discoveries have, on average, longer spin periods than the known pulsar population. This may reflect different selection biases between LOTAAS and previous surveys, though it is also possible that slower-spinning pulsars preferentially have steeper radio spectra. LOTAAS is the deepest all-sky pulsar survey using a digital aperture array; we discuss some of the lessons learned that can inform the approach for similar surveys using future radio telescopes such as the Square Kilometre Array.

Dispersive delays due to the Solar wind introduce excess noise in high-precision pulsar timing experiments, and must be removed in order to achieve the accuracy needed to detect, e.g., low-frequency gravitational waves. In current pulsar timing experiments, this delay is usually removed by approximating the electron density distribution in the Solar wind either as spherically symmetric, or with a two-phase model that describes the contributions from both high- and low-speed phases of the Solar wind. However, no dataset has previously been available to test the performance and limitations of these models over extended timescales and with sufficient sensitivity. Here we present the results of such a test with an optimal dataset of observations of pulsar J0034-0534, taken with the German stations of LOFAR. We conclude that the spherical approximation performs systematically better than the two-phase model at almost all angular distances, with a residual root-mean-square (rms) given by the two-phase model being up to 28% larger than the result obtained with the spherical approximation. Nevertheless, the spherical approximation remains insufficiently accurate in modelling the Solar-wind delay (especially within 20 degrees of angular distance from the Sun), as it leaves timing residuals with rms values that reach the equivalent of 0.3 microseconds at 1400 MHz. This is because a spherical model ignores the large daily variations in electron density observed in the Solar wind. In the short term, broadband observations or simultaneous observations at low frequencies are the most promising way forward to correct for Solar-wind induced delay variations.

The source of high-energy protons (>500 MeV) responsible for the so-called ground level enhancements (GLEs) remains an open question in solar physics. One of the candidates is a shock wave driven by a coronal mass ejection, which is thought to accelerate particles via diffusive-shock acceleration. We perform physics-based simulations of proton acceleration using information on the shock and ambient plasma parameters derived from the observation of a real GLE event. We analyse the simulation results with the aim to find out which of the parameters are significant in controlling the acceleration efficiency and to get a better understanding of the conditions under which the shock can produce relativistic protons. We use results of the recently developed technique to determine the shock and ambient plasma parameters, applied to the 17 May 2012 GLE event, and carry out proton acceleration simulations with the Coronal Shock Acceleration model. We have performed proton acceleration simulations for nine individual magnetic field lines characterised by various plasma conditions. Analysis of the simulation results shows that the acceleration efficiency of the shock, i.e., its ability to accelerate particles to high energies, tends to be higher for those shock portions that are characterised by larger values of the scattering-centre compression ratio and/or the fast-mode Mach number. At the same time, the acceleration efficiency can be strengthened due to enhanced plasma density in the flux tube. Analysis of the delays between the flare onset and the production times of protons of 1 GV rigidity for different field lines in our simulations, and a subsequent comparison of those with the observed values indicate a possibility that quasi-perpendicular portions of the shock play the main role in producing relativistic protons.

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.

Coronal mass ejections (CMEs) and other solar eruptive phenomena can be physically linked by combining data from a multitude of ground-based and space-based instruments alongside models, however this can be challenging for automated operational systems. The EU Framework Package 7 HELCATS project provides catalogues of CME observations and properties from the Helio- spheric Imagers onboard the two NASA/STEREO spacecraft in order to track the evolution of CMEs in the inner heliosphere. From the main HICAT catalogue of over 2,000 CME detections, an automated algorithm has been developed to connect the CMEs observed by STEREO to any corresponding solar flares and active region (AR) sources on the solar surface. CME kinematic properties, such as speed and angular width, are compared with AR magnetic field properties, such as magnetic flux, area, and neutral line characteristics. The resulting LOWCAT catalogue is also compared to the extensive AR property database created by the EU Horizon 2020 FLARECAST project, which provides more complex magnetic field parameters derived from vector magnetograms. Initial statistical analysis has been undertaken on the new data to provide insight into the link between flare and CME events, and characteristics of eruptive ARs. Warning thresholds determined from analysis of the evolution of these parameters is shown to be a useful output for operational space weather purposes. Parameters of particular interest for further analysis include total unsigned flux, vertical current, and current helicity. The automated method developed to create the LOWCAT catalogue may also be useful for future efforts to develop operational CME forecasting.

We study the link between an expanding coronal shock and the energetic particles measured near Earth during the Ground Level Enhancement (GLE) of 17 May 2012. We developed a new technique based on multipoint imaging to triangulate the 3-D expansion of the shock forming in the corona. It uses images from three vantage points by mapping the outermost extent of the coronal region perturbed by the pressure front. We derive for the first time the 3-D velocity vector and the distribution of Mach numbers, $M_{FM}$, of the entire front as a function of time. Our approach uses magnetic field reconstructions of the coronal field, full magneto-hydrodynamic simulations and imaging inversion techniques. We find that the highest $M_{FM}$ values appear near the coronal neutral line within a few minutes of the Coronal Mass Ejection (CME) onset; this neutral line is usually associated with the source of the heliospheric current and plasma sheet. We illustrate the variability of the shock speed, shock geometry and Mach number along different modeled magnetic field lines. Despite the level of uncertainty in deriving the shock Mach numbers, all employed reconstruction techniques show that the release time of GeV particles occurs when the coronal shock becomes super-critical ($M_{FM}>3$). Combining in-situ measurements with heliospheric imagery, we also demonstrate that magnetic connectivity between the accelerator (the coronal shock of 17 May 2012) and the near-Earth environment is established via a magnetic cloud that erupted from the same active region roughly five days earlier.

Context. The Sun is an active source of radio emission ranging from long duration radio bursts associated with solar flares and coronal mass ejections to more complex, short duration radio bursts such as solar S bursts, radio spikes and fibre bursts. While plasma emission is thought to be the dominant emission mechanism for most radio bursts, the electron-cyclotron maser (ECM) mechanism may be responsible for more complex, short-duration bursts as well as fine structures associated with long-duration bursts. Aims. We investigate the conditions for ECM in the solar corona by considering the ratio of the electron plasma frequency {\omega}p to the electron-cyclotron frequency {\Omega}e. The ECM is theoretically possible when {\omega}p/{\Omega}e < 1. Methods. Two-dimensional electron density, magnetic field, plasma frequency, and electron cyclotron frequency maps of the off- limb corona were created using observations from SDO/AIA and SOHO/LASCO, together with potential field extrapolations of the magnetic field. These maps were then used to calculate {\omega}p/{\Omega}e and Alfven velocity maps of the off-limb corona. Results. We found that the condition for ECM emission ({\omega}p/{\Omega}e < 1) is possible at heights < 1.07 R_sun in an active region near the limb; that is, where magnetic field strengths are > 40 G and electron densities are greater than 3x10^8 cm-3. In addition, we found comparatively high Alfv\'en velocities (> 0.02 c or > 6000 km s-1) at heights < 1.07 R_sun within the active region. Conclusions. This demonstrates that the condition for ECM emission is satisfied within areas of the corona containing large magnetic fields, such as the core of a large active region. Therefore, ECM could be a possible emission mechanism for high-frequency radio and microwave bursts.

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.

As the most energetic eruptions in the solar system, coronal mass ejections (CMEs) can produce shock waves at both their front and flanks as they erupt from the Sun into the heliosphere. However, the amount of energy produced in these eruptions, and the proportion of their energy required to produce the waves, is not well characterised. Here we use observations of a solar eruption from 2014 February 25 to estimate the energy budget of an erupting CME and the globally-propagating "EIT wave" produced by the rapid expansion of the CME flanks in the low solar corona. The "EIT wave" is shown using a combination of radio spectra and extreme ultraviolet images to be a shock front with a Mach number greater than one. Its initial energy is then calculated using the Sedov-Taylor blast-wave approximation, which provides an approximation for a shock front propagating through a region of variable density. This approach provides an initial energy estimate of $\approx$2.8$\times$10$^{31}$ ergs to produce the "EIT wave", which is approximately 10% the kinetic energy of the associated CME (shown to be $\approx$2.5$\times$10$^{32}$ ergs). These results indicate that the energy of the "EIT wave" may be significant and must be considered when estimating the total energy budget of solar eruptions.

Using multi-wavelength imaging observations, in EUV, white light and radio, and radio spectral data over a large frequency range, we analyzed the triggering and development of a complex eruptive event. This one includes two components, an eruptive jet and a CME which interact during more than 30 min, and can be considered as physically linked. This was an unusual event. The jet is generated above a typical complex magnetic configuration which has been investigated in many former studies related to the build-up of eruptive jets; this configuration includes fan-field lines originating from a corona null point above a parasitic polarity, which is embedded in one polarity region of large Active Region (AR). The initiation and development of the CME, observed first in EUV, does not show usual signatures. In this case, the eruptive jet is the main actor of this event. The CME appears first as a simple loop system which becomes destabilized by magnetic reconnection between the outer part of the jet and the ambient medium. The progression of the CME is closely associated with the occurrence of two successive types II bursts from distinct origin. An important part of this study is the first radio type II burst for which the joint spectral and imaging observations allowed: i) to follow, step by step, the evolution of the spectrum and of the trajectory of the radio burst, in relationship with the CME evolution; ii) to obtain, without introducing an electronic density model, the B-field and the Alfven speed.

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

Cosmic rays and solar energetic particles may be accelerated to relativistic energies by shock waves in astrophysical plasmas. On the Sun, shocks and particle acceleration are often associated with the eruption of magnetized plasmoids, called coronal mass ejections (CMEs). However, the physical relationship between CMEs and shock particle acceleration is not well understood. Here, we use extreme ultraviolet, radio and white-light imaging of a solar eruptive event on 22 September 2011 to show that a CME-induced shock (Alfv\'en Mach number 2.4$^{+0.7}_{-0.8}$) was coincident with a coronal wave and an intense metric radio burst generated by intermittent acceleration of electrons to kinetic energies of 2-46 keV (0.1-0.4 c). Our observations show that plasmoid-driven quasi-perpendicular shocks are capable of producing quasi-periodic acceleration of electrons, an effect consistent with a turbulent or rippled plasma shock surface.

Super-Alfv\'enic shock waves associated with coronal mass ejections (CMEs) can produce radio emission known as Type II bursts. In the absence of direct imaging, accurate estimates of coronal electron densities, magnetic field strengths and Alfv\'en speeds are required in order to calculate the kinematics of shocks. To date, 1D radial models have been used, but these are not appropriate for shocks propagating in non-radial directions. Here, we study a coronal shock wave associated with a CME and Type II radio burst using 2D electron density and Alfv\'en speed maps to determine the locations that shocks are excited as the CME expands through the corona. Coronal density maps were obtained from emission measures derived from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) and polarized brightness measurements from the Large Angle and Spectrometric Coronagraph (LASCO) on board the Solar and Heliospheric Observatory (SOHO). Alfv\'en speed maps were calculated using these density maps and magnetic field extrapolations from the Helioseismic and Magnetic Imager (SDO/HMI). The computed density and Alfv\'en speed maps were then used to calculate the shock kinematics in non-radial directions. Using the kinematics of the Type II burst and associated shock, we find our observations to be consistent with the formation of a shock located at the CME flanks where the Alfv\'en speed has a local minimum. 1D density models are not appropriate for shocks that propagate non-radially along the flanks of a CME. Rather, the 2D density, magnetic field and Alfv\'en speed maps described here give a more accurate method for determining the fundamental properties of shocks and their relation to CMEs.

The Rosse Solar-Terrestrial Observatory (RSTO; www.rosseobservatory.ie) was established at Birr Castle, Co. Offaly, Ireland (53 05'38.9", 7 55'12.7") in 2010 to study solar radio bursts and the response of the Earth's ionosphere and geomagnetic field. To date, three Compound Astronomical Low-cost Low-frequency Instrument for Spectroscopy and Transportable Observatory (CALLISTO) spectrometers have been installed, with the capability of observing in the frequency range 10-870 MHz. The receivers are fed simultaneously by biconical and log-periodic antennas. Nominally, frequency spectra in the range 10-400 MHz are obtained with 4 sweeps per second over 600 channels. Here, we describe the RSTO solar radio spectrometer set-up, and present dynamic spectra of a sample of Type II, III and IV radio bursts. In particular, we describe fine-scale structure observed in Type II bursts, including band splitting and rapidly varying herringbone features.