Context. Observations of solar type II radio bursts provide a unique opportunity to analyze the non-thermal electrons accelerated by coronal shocks and also to diagnose the plasma density distribution in the corona. However, there are very rare high-frequency resolution interferometric observations for type II radio bursts that are capable of tracking these electrons. Aims. Recently, more spatially resolved high-resolution observations of type II radio bursts have been recorded with the Low-Frequency Array (LOFAR). Using these observations, we aim to track the location of a type II radio burst that experiences a sudden spectral bump. Methods. Here, we present the first radio imaging observations for a type II burst with a spectral bump. We measure the variation in source location and frequency drift of the type II burst, and deduct the density distribution along its propagation direction. Results. We identified a type II burst that experiences a sudden spectral bump in its frequency-time profile. The overall frequency drift rate is 0.06 MHz/s and it corresponds to an estimated speed of 295 km/s. The projected speed of the radio source obtained from imaging is 380 km/s towards the east direction. At the spectral bump, a deviation in the source locations of the type II split bands is observed. The band separation increases significantly in the north-south direction. Conclusions. The spectral bump shows an 8 MHz deviation at 60 MHz which corresponds to a 25% decrease in the plasma density. The estimated crossing distance during the spectrum bump in type II is 29 Mm suggesting that this density variation occurs in a confined area. This indicates that the shock most likely encounters the upper extent of a coronal hole.
M. Marongiu, A. Pellizzoni, S. Righini, S. Mulas, R. Nesti, A. Burtovoi, M. Romoli, G. Serra, G. Valente, E. Egron, G. Murtas, M. N. Iacolina, A. Melis, S. L. Guglielmino, S. Loru, P. Zucca, A. Zanichelli, M. Bachetti, A. Bemporad, F. Buffa, et al (26) One of the most important objectives of solar physics is the physical understanding of the solar atmosphere, the structure of which is also described in terms of the density (N) and temperature (T) distributions of the atmospheric matter. Several multi-frequency analyses show that the characteristics of these distributions are still debated, especially for the outer coronal emission. We aim to constrain the T and N distributions of the solar atmosphere through observations in the centimetric radio domain. We employ single-dish observations from two of the INAF radio telescopes at the K-band frequencies (18 - 26 GHz). We investigate the origin of the significant brightness temperature ($T_B$) level that we detected up to the upper corona ($\sim 800$ Mm of altitude with respect to the photospheric solar surface). To probe the physical origin of the atmospheric emission and to constrain instrumental biases, we reproduced the solar signal by convolving specific 2D antenna beam models. The analysis of the solar atmosphere is performed by adopting a physical model that assumes the thermal bremsstrahlung as the emission mechanism, with specific T and N distributions. The modelled $T_B$ profiles are compared with those observed by averaging solar maps obtained during the minimum of solar activity (2018 - 2020). The T and N distributions are compatible (within $25\%$ of uncertainty) with the model up to $\sim 60$ Mm and $\sim 100$ Mm of altitude, respectively. The analysis of the role of the antenna beam pattern on our solar maps proves the physical nature of the atmospheric emission in our images up to the coronal tails seen in our $T_B$ profiles. The challenging analysis of the coronal radio emission at higher altitudes, together with the data from satellite instruments will require further multi-frequency measurements.
Solar radio U-bursts are generated by electron beams traveling along closed magnetic loops in the solar corona. Low-frequency ($<$ 100 MHz) U-bursts serve as powerful diagnostic tools for studying large-sized coronal loops that extend into the middle corona. However, the positive frequency drift component (descending leg) of U-bursts has received less attention in previous studies, as the descending radio flux is weak. In this study, we utilized LOFAR interferometric solar imaging data from a U-burst that has a significant descending leg component, observed between 10 to 90 MHz on June 5th, 2020. By analyzing the radio source centroid positions, we determined the beam velocities and physical parameters of a large coronal magnetic loop that reached just about 1.3 $\rm{R_{\odot}}$ in altitude. At this altitude, we found the plasma temperature to be around 1.1 MK, the plasma pressure around 0.20 $\rm{mdyn,cm^{-2}}$, and the minimum magnetic field strength around 0.07 G. The similarity in physical properties determined from the image suggests a symmetric loop. The average electron beam velocity on the ascending leg was found to be 0.21 c, while it was 0.14 c on the descending leg. This apparent deceleration is attributed to a decrease in the range of electron energies that resonate with Langmuir waves, likely due to the positive background plasma density gradient along the downward loop leg.
M. Marongiu, A. Pellizzoni, S. Mulas, S. Righini, R. Nesti, G. Murtas, E. Egron, M. N. Iacolina, A. Melis, G. Valente, G. Serra, S. L. Guglielmino, A. Zanichelli, P. Romano, S. Loru, M. Bachetti, A. Bemporad, F. Buffa, R. Concu, G. L. Deiana, et al (15) The Sun is an extraordinary workbench, from which several fundamental astronomical parameters can be measured with high precision. Among these parameters, the solar radius $R_{\odot}$ plays an important role in several aspects, such as in evolutionary models. Despite the efforts in obtaining accurate measurements of $R_{\odot}$, the subject is still debated and measurements are puzzling and/or lacking in many frequency ranges. We aimed to determine the mean, equatorial, and polar radii of the Sun ($R_c$, $R_{eq}$, and $R_{pol}$) in the frequency range 18.1 - 26.1 GHz. We employed single-dish observations from the newly-appointed Medicina "Gavril Grueff" Radio Telescope and the Sardinia Radio Telescope (SRT) throughout 5 years, from 2018 to mid-2023, in the framework of the SunDish project for solar monitoring. Two methods to calculate the radius at radio frequencies are considered and compared. To assess the quality of our radius determinations, we also analysed the possible degrading effects of the antenna beam pattern on our solar maps, using two 2D-models. We carried out a correlation analysis with the evolution of the solar cycle through the calculation of Pearson's correlation coefficient $\rho$. We obtained several values for the solar radius - ranging between 959 and 994 arcsec - and $\rho$, with typical errors of a few arcsec. Our $R_{\odot}$ measurements, consistent with values reported in literature, suggest a weak prolatness of the solar limb ($R_{eq}$ > $R_{pol}$), although $R_{eq}$ and $R_{pol}$ are statistically compatible within 3$\sigma$ errors. The correlation analysis using the solar images from Grueff shows (1) a positive correlation between the solar activity and the temporal variation of $R_c$ (and $R_{eq}$) at all observing frequencies, and (2) a weak anti-correlation between the temporal variation of $R_{pol}$ and the solar activity at 25.8 GHz.
The Parker Solar Probe (PSP) provides us the unprecedentedly close approach observation to the Sun, and hence the possibility of directly understanding the "elementary process" which occurs in the kinetic scale of particles collective interactioin in solar coronal plasmas. We reported a kind of weak solar radio bursts (SRBs), which are detected by PSP when it passed a low-density magnetic channel during its second encounter phase. These weak SRBs have low starting frequecny $\sim 20$ MHz and narrow frequency range from a few tens MHz to a few hundres kHz. Their dynamic spectra display a strongly evolving feature of the intermediate relative drift rate decreasing rapidly from above 0.01/s to below 0.01/s. Analyses based on common empirical models of solar coronal plasmas indicate that these weak SRBs originate from the heliocentric distance $\sim 1.1-6.1~R_S$ (the solar radius), a typical solar wind acceleration region with a low-$\beta$ plasma, and indicate that their soruces have a typic motion velociy $\sim v_A$ (Alfvén velocity) obviously lower than that of fast electrons required by effectively exciting SRBs. We propose that solitary kinetic Alfvén waves with kinetic scales can be responsible for the generation of these small-scalevweak SRBs, called solitary wave radiation (SWR).
This study aims to investigate the ambiguous source and the underlying physical processes of the solar type III radio bursts that occurred on April 3, 2019, through the utilization of multiwavelength observations from the LOFAR radio telescope and the PSP space mission, as well as incorporating results from a PFSS and MHD models. The primary goal is to identify the spatial and temporal characteristics of the radio sources, as well as the plasma conditions along their trajectory. Data preprocessing techniques are applied to combine high- and low-frequency observations from LOFAR and PSP between 2.6 kHz and 80 MHz. We then extract information on the frequency drift and speed of the accelerated electron beams from the dynamic spectra. Additionally, we use LOFAR interferometric observations to image the sources of the radio emission at multiple frequencies and determine their locations and kinematics in the corona. Lastly, we analyze the plasma parameters and magnetic field along the trajectories of the radio sources using PFSS and MHD model results. We present several notable findings related to type III radio bursts. Firstly, through our automated implementation, we were able to effectively identify and characterize 9 type III radio bursts in the LOFAR-PSP combined dynamic spectrum and 16 type III bursts in the LOFAR dynamic spectrum. Secondly, our imaging observations show that the electrons responsible for these bursts originate from the same source and within a short time frame of fewer than 30 minutes. Finally, our analysis provides informative insights into the physical conditions along the path of the electron beams. For instance, we found that the plasma density obtained from the MAS model is significantly lower than the expected theoretical density.
In space weather studies and forecasting we employ magnetohydrodynamic (MHD) simulations which can provide rather accurate reconstruction of the solar wind dynamics and its evolution. However, all MHD simulations are restricted by the input data and the modelled solar wind characteristics need to be validated with different types of observations. That is very difficult, in particular for the solar wind characteristics close to the Sun, since the majority of in-situ observations are taken in the vicinity of the Earth. This is why all alternative methods for estimation of solar wind plasma characteristics are very important. In this study we utilise low radio frequency observations of pulsars to probe the total electron content along the line of sight. For the first time, we compare density estimates from pulsars with predictions from the 3D MHD modelling code; the EUropean Heliospheric FORecasting Information Asset (EUHFORIA). We find a very good correlation for the solar wind density along a given line of sight obtained by EUHFORIA and pulsar observations. We also demonstrate that the pulsar observations can be very useful not only for the model validation but also for understanding its limitations.
Solar-wind 3-D reconstruction tomography based on interplanetary scintillation (IPS) studies provides fundamental information for space-weather forecasting models, and gives the possibility to determine heliospheric column densities. Here we compare the time series of Solar-wind column densities derived from long-term observations of pulsars, and the Solar-wind reconstruction provided by the UCSD IPS tomography. This work represents a completely independent comparison and validation of these techniques to provide this measurement, and it strengthens confidence in the use of both in space-weather analyses applications.
Radio spectroscopy provides a unique inspection perspective for solar and space weather research, which can reveal the plasma and energetic electron information in the solar corona and inner heliosphere. However, Radio-Frequency Interference (RFI) from human activities affects sensitive radio telescopes, and significantly affects the quality of observation. Thus, RFI detection and mitigation for the observations is necessary to obtain high quality, science-ready data. The flagging of RFI is particularly challenging for the solar and space weather observations at low frequency, because the solar radio bursts can be brighter than the RFI, and may show similar temporal behavior. In this work, we investigate RFI flagging methods for solar and space weather observations, including a strategy for AOFlagger, and a novel method that makes use of a morphology convolution. These algorithms can effectively flag RFI while preserving solar radio bursts.
Large coronal loops around one solar radius in altitude are an important connection between the solar wind and the low solar corona. However, their plasma properties are ill-defined as standard X-ray and UV techniques are not suited to these low-density environments. Diagnostics from type J solar radio bursts at frequencies above 10 MHz are ideally suited to understand these coronal loops. Despite this, J-bursts are less frequently studied than their type III cousins, in part because the curvature of the coronal loop makes them unsuited for using standard coronal density models. We used LOw-Frequency-ARray (LOFAR) and Parker Solar Probe (PSP) solar radio dynamic spectrum to identify 27 type III bursts and 27 J-bursts during a solar radio noise storm observed on 10 April 2019. We found that their exciter velocities were similar, implying a common acceleration region that injects electrons along open and closed magnetic structures. We describe a novel technique to estimate the density model in coronal loops from J-burst dynamic spectra, finding typical loop apex altitudes around 1.3 solar radius. At this altitude, the average scale heights were 0.36 solar radius, the average temperature was around 1 MK, the average pressure was 0.7 mdyn cm$^{-2}$, and the average minimum magnetic field strength was 0.13 G. We discuss how these parameters compare with much smaller coronal loops.
Bartosz Dabrowski, Katarzyna Mikula, Pawel Flisek, Christian Vocks, PeiJin Zhang, Jasmina Magdalenić, Alexander Warmuth, Diana E. Morosan, Adam Froń, Richard A. Fallows, Mario M. Bisi, Andrzej Krankowski, Gottfried Mann, Leszek Blaszkiewicz, Eoin P. Carley, Peter T. Gallagher, Pietro Zucca, Pawel Rudawy, Marcin Hajduk, Kacper Kotulak, et al (1) The Sun is the source of different types of radio bursts that are associated with solar flares, for example. Among the most frequently observed phenomena are type III solar bursts. Their radio images at low frequencies (below 100 MHz) are relatively poorly studied due to the limitations of legacy radio telescopes. We study the general characteristics of types IIIb and U with stria structure solar radio bursts in the frequency range of 20 - 80 MHz, in particular the source size and evolution in different altitudes, as well as the velocity and energy of electron beams responsible for their generation. In this work types IIIb and U with stria structure radio bursts are analyzed using data from the LOFAR telescope including dynamic spectra and imaging observations, as well as data taken in the X-ray range (GOES and RHESSI satellites) and in the extreme ultraviolet (SDO satellite). In this study we determined the source size limited by the actual shape of the contour at particular frequencies of type IIIb and U solar bursts in a relatively wide frequency band from 20 to 80 MHz. Two of the bursts seem to appear at roughly the same place in the studied active region and their source sizes are similar. It is different in the case of another burst, which seems to be related to another part of the magnetic field structure in this active region. The velocities of the electron beams responsible for the generation of the three bursts studied here were also found to be different.
Observations of interplanetary scintillation (IPS - the scintillation of compact radio sources due to density variations in the solar wind) enable the velocity of the solar wind to be determined, and its bulk density to be estimated, throughout the inner heliosphere. A series of observations using the Low Frequency Array (LOFAR - a radio telescope centred on the Netherlands with stations across Europe) were undertaken using this technique to observe the passage of an ultra-fast CME which launched from the Sun following the X-class flare of 10 September 2017. LOFAR observed the strong radio source 3C147 at an elongation of 82 degrees from the Sun over a period of more than 30 hours and observed a strong increase in speed to 900km/s followed two hours later by a strong increase in the level of scintillation, interpreted as a strong increase in density. Both speed and density remained enhanced for a period of more than seven hours, to beyond the period of observation. Further analysis of these data demonstrates a view of magnetic-field rotation due to the passage of the CME, using advanced IPS techniques only available to a unique instrument such as LOFAR.
R.A. Fallows, B. Forte, M. Mevius, M. A. Brentjens, C. G. Bassa, M. M. Bisi, A. Offringa, G. Shaifullah, C. Tiburzi, H. Vedantham, P. Zucca Context. The occultation of a radio source by the plasma tail of a comet can be used to probe structure and dynamics in the tail. Such occultations are rare, and the occurrence of scintillation, due to small-scale density variations in the tail, remains somewhat controversial. Aims. A detailed observation taken with the Low-Frequency Array (LOFAR) of a serendipitous occultation of the compact radio source 3C196 by the plasma tail of comet C/2020 F3 (Neowise) is presented. 3C196 tracked almost perpendicularly behind the tail, providing a unique profile cut only a short distance downstream from the cometary nucleus itself. Methods. Interplanetary scintillation (IPS) is observed as the rapid variation of the intensity received of a compact radio source due to density variations in the solar wind. IPS in the signal received from 3C196 was observed for five hours, covering the full transit behind the plasma tail of comet C/2020 F3 (Neowise) on 16 July 2020, and allowing an assessment of the solar wind in which the comet and its tail are embedded. Results. The results reveal a sudden and strong enhancement in scintillation which is unequivocally attributable to the plasma tail. The strongest scintillation is associated with the tail boundaries, weaker scintillation is seen within the tail, and previously-unreported periodic variations in scintillation are noted, possibly associated with individual filaments of plasma. Furthermore, contributions from the solar wind and comet tail are separated to measure a sharp decrease in the velocity of material within the tail, suggesting a steep velocity shear resulting in strong turbulence along the tail boundary
Type IV radio burst has been studied for over 50 years. However, the specifics of the radio emission mechanisms is still an open question. In order to provide more information about the emission mechanisms, we studied a moving type IV radio burst with fine structures (spike group) by using the high resolution capability of Low-Frequency Array (LOFAR) on Aug 25, 2014\textbf (SOLA-D-21-00188). We present a comparison of Nançay RadioHeliograph (NRH) and the first LOFAR imaging data of type IV radio burst. The degree of circular polarization (DCP) is calculated at frequencies in the range 20$\sim$180 MHz using LOFAR data, and it was found that the value of DCP gradually increased during the event, with values of 10\%$\sim$20\%. LOFAR interferometric data were combined with white light observations in order to track the propagation of this type IV. The kinematics shows a westward motion of the radio sources, slower than the CME leading edge. The dynamic spectrum of LOFAR shows a large number of fine structures with duration of less than 1s and high brightness temperature ($T_\mathrm{B}$), i.e. $10^{12}$$\sim$$10^{13}$ K. The gradual increase of DCP supports gyrosynchrotron emission as the most plausible mechanism for the type IV. However, coherent emissions such as Electron Cyclotron Maser (ECM) instability can be responsible for small scale fine structures. Countless fine structures altogether were responsible for such high $T_\mathrm{B}$.
Radio emission of the quiet Sun is considered to be due to thermal bremsstrahlung emission of the hot solar atmosphere. The properties of the quiet Sun in the microwave band have been well studied, and they can be well described by the spectrum of bremsstrahlung emission. In the meter-wave and decameter-wave bands, properties of the quiet Sun have rarely been studied due to the instrumental limitations. In this work, we use the LOw Frequency ARray (LOFAR) telescope to perform high quality interferometric imaging spectroscopy observations of quiet Sun coronal emission at frequencies below 90~MHz. We present the brightness temperature spectrum, and size of the Sun in the frequency range of 20-80~MHz. We report on dark coronal regions with low brightness temperature that persist with frequency. The brightness temperature spectrum of the quiet Sun is discussed and compared with the bremsstrahlung emission of a coronal model and previous quiet Sun observations.
A. Pellizzoni, S. Righini, M. N. Iacolina, M. Marongiu, S. Mulas, G. Murtas, G. Valente, E. Egron, M. Bachetti, F. Buffa, R. Concu, G. L. Deiana, S. L. Guglielmino, A. Ladu, S. Loru, A. Maccaferri, P. Marongiu, A. Melis, A. Navarrini, A. Orfei, et al (11) We present a new solar radio imaging system implemented through the upgrade of the large single-dish telescopes of the Italian National Institute for Astrophysics (INAF), not originally conceived for solar observations. During the development and early science phase of the project (2018-2020), we obtained about 170 maps of the entire solar disk in the 18-26 GHz band, filling the observational gap in the field of solar imaging at these frequencies. These solar images have typical resolutions in the 0.7-2 arcmin range and a brightness temperature sensitivity <10 K. Accurate calibration adopting the Supernova Remnant Cas A as a flux reference, provided typical errors <3% for the estimation of the quiet-Sun level components and for active regions flux measurements. As a first early science result of the project, we present a catalog of radio continuum solar imaging observations with Medicina 32-m and SRT 64-m radio telescopes including the multi-wavelength identification of active regions, their brightness and spectral characterization. The interpretation of the observed emission as thermal bremsstrahlung components combined with gyro-magnetic variable emission pave the way to the use of our system for long-term monitoring of the Sun. We also discuss useful outcomes both for solar physics (e.g. study of the chromospheric network dynamics) and space weather applications (e.g. flare precursors studies).
Diana E. Morosan, Juska E. Räsänen, Anshu Kumari, Emilia K. J. Kilpua, Mario M. Bisi, Bartosz Dabrowski, Andrzej Krankowski, Jasmina Magdalenić, Gottfried Mann, Hanna Rothkaehl, Christian Vocks, Pietro Zucca The Sun is an active star that often produces numerous bursts of electromagnetic radiation at radio wavelengths. Low frequency radio bursts have recently been brought back to light with the advancement of novel radio interferometers. However, their polarisation properties have not yet been explored in detail, especially with the Low Frequency Array (LOFAR), due to difficulties in calibrating the data and accounting for instrumental leakage. Here, using a unique method to correct the polarisation observations, we explore the circular polarisation of different sub-types of solar type III radio bursts and a type I noise storm observed with LOFAR, which occurred during March-April 2019. We analysed six individual radio bursts from two different dates. We present the first Stokes V low frequency images of the Sun with LOFAR in tied-array mode observations. We find that the degree of circular polarisation for each of the selected bursts increases with frequency for fundamental emission, while this trend is either not clear or absent for harmonic emission. The type III bursts studied, that are part of a long--lasting type III storm, can have different senses of circular polarisation, occur at different locations and have different propagation directions. This indicates that the type III bursts forming a classical type III storm do not necessarily have a common origin but instead they indicate the existence of multiple, possibly unrelated, acceleration processes originating from solar minimum active regions.
The solar corona below 10 solar radii is an important region for early acceleration and transport of solar energetic particles (SEPs) by coronal mass ejection-driven shock waves. There, these waves propagate into a highly variable dynamic medium with steep gradients and rapidly expanding coronal magnetic fields, which modulates the particle acceleration near the shock or wave surfaces, and the way SEPs spread into the heliosphere. We present a study modelling the acceleration of SEPs in global coronal shock events in the corona, as well as their transport to 1 au, based on telescopic observations coupled with dynamic physical models.
L. K. Morabito, N. J. Jackson, S. Mooney, F. Sweijen, S. Badole, P. Kukreti, D. Venkattu, C. Groeneveld, A. Kappes, E. Bonnassieux, A. Drabent, M. Iacobelli, J. H. Croston, P. N. Best, M. Bondi, J. R. Callingham, J. E. Conway, A. T. Deller, M. J. Hardcastle, J. P. McKean, et al (54) [abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $\mu$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $\mu$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit.
Neal Jackson, Shruti Badole, John Morgan, Rajan Chhetri, Kaspars Prusis, Atvars Nikolajevs, Leah Morabito, Michiel Brentjens, Frits Sweijen, Marco Iacobelli, Emanuela Orrù, J. Sluman, R. Blaauw, H. Mulder, P. van Dijk, Sean Mooney, Adam Deller, Javier Moldon, J.R. Callingham, Jeremy Harwood, et al (52) The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS.
Eoin P. Carley, Baptiste Cecconi, Hamish A. Reid, Carine Briand, Sasikumar Raja, Sophie Masson, Vladimir V. Dorovskyy, Caterina Tiburzi, Nicole Vilmer, Pietro Zucca, Philippe Zarka, Michel Tagger, Jean-Mathias Griessmeier, Stéphane Corbel, Gilles Theureau, Alan Loh, Julien Girard Eruptive activity in the solar corona can often lead to the propagation of shock waves. In the radio domain the primary signature of such shocks are type II radio bursts, observed in dynamic spectra as bands of emission slowly drifting towards lower frequencies over time. These radio bursts can sometimes have inhomogeneous and fragmented fine structure, but the cause of this fine structure is currently unclear. Here we observe a type II radio burst on 2019-March-20th using the New Extension in Nançay Upgrading LOFAR (NenuFAR), a radio interferometer observing between 10-85 MHz. We show that the distribution of size-scales of density perturbations associated with the type II fine structure follows a power law with a spectral index in the range of $\alpha=-1.7$ to -2.0, which closely matches the value of $-5/3$ expected of fully developed turbulence. We determine this turbulence to be upstream of the shock, in background coronal plasma at a heliocentric distance of $\sim$2 R$_{\odot}$. The observed inertial size-scales of the turbulent density inhomogeneities range from $\sim$62 Mm to $\sim$209 km. This shows that type II fine structure and fragmentation can be due to shock propagation through an inhomogeneous and turbulent coronal plasma, and we discuss the implications of this on electron acceleration in the coronal shock.
A. Miraval Zanon, P. D'Avanzo, A. Ridolfi, F. Coti Zelati, S. Campana, C. Tiburzi, D. de Martino, T. Muñoz Darias, C. G. Bassa, L. Zampieri, A. Possenti, F. Ambrosino, A. Papitto, M. C. Baglio, M. Burgay, A. Burtovoi, D. Michilli, P. Ochner, P. Zucca We present simultaneous multiwavelength observations of the 4.66 ms redback pulsar PSR J1048+2339. We performed phase-resolved spectroscopy with the Very Large Telescope (VLT) searching for signatures of a residual accretion disk or intra-binary shock emission, constraining the companion radial velocity semi-amplitude ($K_2$), and estimating the neutron star mass ($M_{\rm NS}$). Using the FORS2-VLT intermediate-resolution spectra, we measured a companion velocity of $291 < K_2 < 348$ km s$^{-1}$ and a binary mass ratio of $0.209 < q < 0.250$. Combining our results for $K_2$ and $q$, we constrained the mass of the neutron star and the companion to $(1.0 < M_{\rm NS} < 1.6){\rm sin}^{-3}i\,M_{\odot}$ and $(0.24 < M_2 < 0.33){\rm sin}^{-3}i\,M_{\odot}$, respectively, where $i$ is the system inclination. The Doppler map of the H$\alpha$ emission line exhibits a spot feature at the expected position of the companion star and an extended bright spot close to the inner Lagrangian point. We interpret this extended emission as the effect of an intra-binary shock originating from the interaction between the pulsar relativistic wind and the matter leaving the companion star. The mass loss from the secondary star could be either due to Roche-lobe overflow or to the ablation of its outer layer by the energetic pulsar wind. Contrastingly, we find no evidence for an accretion disk. We report on the results of the SRT and the LOFAR simultaneous radio observations at three different frequencies (150 MHz, 336 MHz, and 1400 MHz). No pulsed radio signal is found in our search. This is probably due to both scintillation and the presence of material expelled from the system which can cause the absorption of the radio signal at low frequencies. Finally, we report on an attempt to search for optical pulsations using IFI+Iqueye mounted at the 1.2 m Galileo telescope at the Asiago Observatory.
A. M. Ryan, P. T. Gallagher, E. P. Carley, M. A. Brentjens, P. C. Murphy, C. Vocks, D. E. Morosan, H. Reid, J. Magdalenic, F. Breitling, P. Zucca, R. Fallows, G. Mann, A. Kerdraon, R. Halfwerk The solar corona is a highly-structured plasma which can reach temperatures of more than ~2 MK. At low frequencies (decimetric and metric wavelengths), scattering and refraction of electromagnetic waves are thought to considerably increase the imaged radio source sizes (up to a few arcminutes). However, exactly how source size relates to scattering due to turbulence is still subject to investigation. The theoretical predictions relating source broadening to propagation effects have not been fully confirmed by observations due to the rarity of high spatial resolution observations of the solar corona at low frequencies. Here, the LOw Frequency ARray (LOFAR) was used to observe the solar corona at 120-180 MHz using baselines of up to ~3.5 km (corresponding to a resolution of ~1-2') during the partial solar eclipse of 2015 March 20. A lunar de-occultation technique was used to achieve higher spatial resolution (~0.6') than that attainable via standard interferometric imaging (~2.4'). This provides a means of studying the contribution of scattering to apparent source size broadening. It was found that the de-occultation technique reveals a more structured quiet corona that is not resolved from standard imaging, implying scattering may be overestimated in this region when using standard imaging techniques. However, an active region source was measured to be ~4' using both de-occultation and standard imaging. This may be explained by the increased scattering of radio waves by turbulent density fluctuations in active regions, which is more severe than in the quiet Sun.
The Sun produces highly dynamic and eruptive events that can drive shocks through the corona. These shocks can accelerate electrons, which result in plasma emission in the form of a type II radio burst. Despite the large number of type II radio bursts observations, the precise origin of coronal shocks is still subject to investigation. Here we present a well observed solar eruptive event that occurred on 16 October 2015, focusing on a jet observed in the extreme ultraviolet (EUV) by the Atmospheric Imaging Assembly (SDO/AIA), a streamer observed in white-light by the Large Angle and Spectrometric Coronagraph (SOHO/LASCO), and a metric type II radio burst observed by the LOw Frequency Array (LOFAR). LOFAR interferometrically imaged the fundamental and harmonic sources of a type II radio burst and revealed that the sources did not appear to be co-spatial, as would be expected from the plasma emission mechanism. We correct for the separation between the fundamental and harmonic using a model which accounts for scattering of radio waves by electron density fluctuations in a turbulent plasma. This allows us to show the type II radio sources were located $\sim$0.5 R$_\odot$ above the jet and propagated at a speed of $\sim$1000 kms$^{-1}$, which was significantly faster than the jet speed of $\sim$200 kms$^{-1}$. This suggests that the type II burst was generated by a piston shock driven by the jet in the low corona.
C. Tiburzi, G. M. Shaifullah, C. G. Bassa, P. Zucca, J. P. W. Verbiest, N. K. Porayko, E. van der Wateren, R. A. Fallows, R. A. Main, G. H. Janssen, J. M. Anderson, A-.S. Bak Nielsen, J. Y. Donner, E. F. Keane, J. Künsemöller, S. Osłowski, J-.M. Grießmeier, M. Serylak, M. Brüggen, B. Ciardi, et al (5) High-precision pulsar timing requires accurate corrections for dispersive delays of radio waves, parametrized by the dispersion measure (DM), particularly if these delays are variable in time. In a previous paper we studied the Solar-wind (SW) models used in pulsar timing to mitigate the excess of DM annually induced by the SW, and found these to be insufficient for high-precision pulsar timing. Here we analyze additional pulsar datasets to further investigate which aspects of the SW models currently used in pulsar timing can be readily improved, and at what levels of timing precision SW mitigation is possible. Our goals are to verify: a) whether the data are better described by a spherical model of the SW with a time-variable amplitude rather than a time-invariant one as suggested in literature, b) whether a temporal trend of such a model's amplitudes can be detected. We use the pulsar-timing technique on low-frequency pulsar observations to estimate the DM and quantify how this value changes as the Earth moves around the Sun. Specifically, we monitor the DM in weekly to monthly observations of 14 pulsars taken with LOFAR across time spans of up to 6 years. We develop an informed algorithm to separate the interstellar variations in DM from those caused by the SW and demonstrate the functionality of this algorithm with extensive simulations. Assuming a spherically symmetric model for the SW density, we derive the amplitude of this model for each year of observations. We show that a spherical model with time-variable amplitude models the observations better than a spherical model with constant amplitude, but that both approaches leave significant SW induced delays uncorrected in a number of pulsars in the sample. The amplitude of the spherical model is found to be variable in time, as opposed to what has been previously suggested.
Low frequency radio wave scattering and refraction can have a dramatic effect on the observed size and position of radio sources in the solar corona. The scattering and refraction is thought to be due to fluctuations in electron density caused by turbulence. Hence, determining the true radio source size can provide information on the turbulence in coronal plasma. However, the lack of high spatial resolution radio interferometric observations at low frequencies, such as with the LOw Frequency ARray (LOFAR), has made it difficult to determine the true radio source size and level of radio wave scattering. Here we directly fit the visibilities of a LOFAR observation of a Type IIIb radio burst with an elliptical Gaussian to determine its source size and position. This circumvents the need to image the source and then de-convolve LOFAR's point spread function, which can introduce spurious effects to the source size and shape. For a burst at 34.76 MHz, we find full width at half maximum (FWHM) heights along the major and minor axes to be $18.8^\prime$ $\pm~0.1^\prime$ and $10.2^\prime$ $\pm~0.1^\prime$, respectively, at a plane of sky heliocentric distance of 1.75 R$_\odot$. Our results suggest that the level of density fluctuations in the solar corona is the main cause of the scattering of radio waves, resulting in large source sizes. However, the magnitude of $\varepsilon$ may be smaller than what has been previously derived in observations of radio wave scattering in tied-array images.
We present a python-based tool to detect the occultation of background sources by foreground Solar coronal mass ejections. The tool takes as input standard celestial coordinates of the source and translates those to the Helioprojective plane, and is thus well suited for use with a wide variety of background astronomical sources. This tool provides an easy means to search through a large archival dataset for such crossings and relies on the well-tested Astropy and Sunpy modules.
PeiJin Zhang, Pietro Zucca, Sarrvesh Seethapuram Sridhar, ChuanBing Wang, Diana E. Morosan, Bartosz Dabrowski, Andrzej Krankowski, Mario M. Bisi, Jasmina Magdalenic, Christian Vocks, Gottfried Mann 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.
J. W. Broderick, T. W. Shimwell, K. Gourdji, A. Rowlinson, S. Nissanke, K. Hotokezaka, P. G. Jonker, C. Tasse, M. J. Hardcastle, J. B. R. Oonk, R. P. Fender, R. A. M. J. Wijers, A. Shulevski, A. J. Stewart, S. ter Veen, V. A. Moss, M. H. D. van der Wiel, D. A. Nichols, A. Piette, M. E. Bell, et al (49) 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.
Richard A. Fallows, Biagio Forte, Ivan Astin, Tom Allbrook, Alex Arnold, Alan Wood, Gareth Dorrian, Maaijke Mevius, Hanna Rothkaehl, Barbara Matyjasiak, Andrzej Krankowski, James M. Anderson, Ashish Asgekar, I. Max Avruch, Mark Bentum, Mario M. Bisi, Harvey R. Butcher, Benedetta Ciardi, Bartosz Dabrowski, Sieds Damstra, et al (38) 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.
F. de Gasperin, J. Vink, J.P. McKean, A. Asgekar, M.J. Bentum, R. Blaauw, A. Bonafede, M. Bruggen, F. Breitling, W.N. Brouw, H.R. Butcher, B. Ciardi, V. Cuciti, M. de Vos, S. Duscha, J. Eisloffel, D. Engels, R.A. Fallows, T.M.O. Franzen, M.A. Garrett, et al (32) 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.
PeiJin Zhang, Pietro Zucca, ChuanBing Wang, Mario M. Bisi, Bartosz Dabrowski, Richard A. Fallows, Andrzej Krankowski, Jasmina Magdalenic, Gottfried Mann, Diana E. Morosan, Christian Vocks 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.
Diana E. Morosan, Eoin P. Carley, Laura A. Hayes, Sophie A. Murray, Pietro Zucca, Richard A. Fallows, Joe McCauley, Emilia K. J. Kilpua, Gottfried Mann, Christian Vocks, Peter T. Gallagher 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.
S. Sanidas, S. Cooper, C. G. Bassa, J. W. T. Hessels, V. I. Kondratiev, D. Michilli, B. W. Stappers, C. M. Tan, J. van Leeuwen, L. Cerrigone, R. A. Fallows, M. Iacobelli, E. Orru, R. F. Pizzo, A. Shulevski, M. C. Toribio, S. ter Veen, P. Zucca, L. Bondonneau, J.-M. Griessmeier, et al (3) 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.
C. Tiburzi, J. P. W. Verbiest, G. M. Shaifullah, G. H. Janssen, J. M. Anderson, A. Horneffer, J. Kuensemoeller, S. Oslowski, J. Y. Donner, M. Kramer, A. Kumari, N. K. Porayko, P. Zucca, B. Ciardi, R.-J. Dettmar, J.-M. Griessmeier, M. Hoeft, M. Serylak 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.
P. Zucca, D. E. Morosan, A. P. Rouillard, R. Fallows, P. T. Gallagher, J. Magdalenic, K-L. Klein, G. Mann, C. Vocks, E. P. Carley, M. M. Bisi, E. P. Kontar, H. Rothkaehl, B. Dabrowski, A. Krankowski, J. Anderson, A. Asgekar, M. E. Bell, M. J. Bentum, P. Best, et al (45) 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.
Alexis P. Rouillard, Illya Plotnikov, Rui F. Pinto, Margot Tirole, Michael Lavarra, Pietro Zucca, Rami Vainio, Allan J. Tylka, Angelos Vourlidas, Marc De Rosa, Jon Linker, Alexander Warmuth, Gottfried Mann, Christina M. Cohen, Robert A. Mewaldt 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 \omegap to the electron-cyclotron frequency \Omegae. The ECM is theoretically possible when \omegap/\Omegae < 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 \omegap/\Omegae and Alfven velocity maps of the off-limb corona. Results. We found that the condition for ECM emission (\omegap/\Omegae < 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én 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.
D. E. Morosan, P. T. Gallagher, P. Zucca, A. O'Flannagain, R. Fallows, H. Reid, J. Magdalenic, G. Mann, M. M. Bisi, A. Kerdraon, A. A. Konovalenko, A. L. MacKinnon, H. O. Rucker, B. Thide, C. Vocks, A. Alexov, J. Anderson, A. Asgekar, I. M. Avruch, M. J. Bentum, et al (33) 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.
D. E. Morosan, P. T. Gallagher, P. Zucca, R. Fallows, E. P. Carley, G. Mann, M. M. Bisi, A. Kerdraon, A. A. Konovalenko, A. L. MacKinnon, H. O. Rucker, B. Thidé, J. Magdalenić, C. Vocks, H. Reid, J. Anderson, A. Asgekar, I. M. Avruch, M. J. Bentum, G. Bernardi, et al (65) 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én 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énic 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én 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én 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én 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én 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én 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én 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.