\section{Introduction}
Meteors are well-known to reflect artificial radio emission in 
forward and backward scattering. Combined optical and radio meteor 
scatter observation of meteors go back to \citep{1947MNRAS.107..155P}, 
showing that meteor trains are overdense for most visible meteors. Forward
meteor scatter, where transmitter and receiver are not at the same location,
can be used to study meteors as is done in, e.g., the 
Belgian RAdio Meteor Stations (BRAMS) network \citep{2011MSSB}. 

Less well established is the detection of intrinsic natural radio emission 
from the meteor or meteor train itself. Early efforts to detect intrinsic
emisison relied on temporal coincidences (e.g., \citep{2000EMP82-83}).
The first more substantive reports based on radio imaging
came from all-sky imaging at low radio frequencies ($\sim$40\,MHz) with the Long
Wavelength Array, detecting intrinsic non-thermal radio emission
from fireballs 
\citep{2014ApJ...788L..26O,2016AGUFM.P33E..06O,2016JGRA..121.6808O}. The
emission persisted well after the meteor itself had faded. The
emission mechanism of these meteor radio afterglows is not fully
understood, and so far has not been independently confirmed with other
instruments. A survey with the Murchison Widefield Array
\citep{2018MNRAS.477.5167Z} did not show this intrinsic
emission at higher 72--103\,MHz frequencies.

Several mechanisms for broadband radio emission from meteors have been suggested in the
literature, namely reflected broadband terrestrial radio emission, 
intrinsic emission from chemically produced suprathermal electrons
\citep{2020JGRA..12528053O}, Langmuir waves
\citep{2015JGRA..120.9916O}, free-free emission
\citep{2018Ge&Ae..58..693F} and transition radiation
\citep{2020JGRA..12528053O}. Also, bright celestial radio emission shines on the meteor
trains and that radiation may be scattered, as suggested in \cite{2015JGRA..120.9916O}.

Here we report on low-frequency all-sky observations with the LOFAR
(``Low Frequency Array'') radio telescope \citep{2013A&A...556A...2V}
during the 2020 Perseid meteor shower \citep{2006MeteorShowers}. 
The radio data are complemented with
simultaneous, both temporally and spatially, optical video observations
from the CAMS BeNeLux network \citep{2011Icar..216...40J} in an effort to
study the altitude dependence of the proposed intrinsic radio emission.

\section{Observations}
The LOFAR radio telescope consists of thousands of dipole antennas
grouped into stations spread over Europe, with a dense core of
stations located in the North of the Netherlands. Each LOFAR station
has low-band antennas (LBAs) which can observe from 10 to 90\,MHz and
high-band antennas (HBAs) operating from 110 to 250\,MHz. During
regular LOFAR LBA observations radio signals are combined (beamformed)
hierarchically, by first combining radio signals from dipole antennas
in a station, and then combining or correlating the combined radio
signals from stations to increase sensitivity at the expense of
field-of-view.

The AARTFAAC piggyback instrument
\citep{2016JAI.....541008P,2019MNRAS.482.2502K,2021MNRAS.505.2966K}
provides a special observing mode whereby LOFAR can create all-sky
images of the radio sky. It does so by correlating the radio signals
from all 576 LBA dipoles from the inner 12 LOFAR stations located in
the dense core (with a 1.2\,km diameter) surrounding the ''super terp'' 
in Drenthe, the Netherlands, see Fig.\,\ref{fig:camsmap}.

LOFAR/AARTFAAC observations were obtained during the Perseid meteor
shower (2020 August 12-13) for 64 frequency channels of 48.8\,kHz
bandwidth (grouped into 16 subbands of 4 channels each), spread over
radio frequencies between 30 to 60\,MHz. These observations were used
to create all-sky images with $3\hbox{$^\prime$}$ spatial resolution and 2\,s
time resolution.

\begin{figure}
  \includegraphics[width=\linewidth]{figures/aartfaac.png}
  \includegraphics[width=\linewidth]{figures/meteormap}
  \label{fig:meteormap}
  \caption{Top: map of the LOFAR core with the location of all 512 dipole
    antennas used for AARTFAAC. Bottom: map with CAMS camera locations (red),
    the location of the LOFAR core (blue) and co-observed meteors.
  }
  \label{fig:camsmap}
\end{figure}

\begin{figure}[!ht]
  \centering
  \parbox{.72\linewidth}{
  \includegraphics[width=\hsize]{77-large.png}}
  \parbox{.243\linewidth}{\includegraphics[width=\hsize]{144-poststamp.png}
  \includegraphics[width=\hsize]{235-poststamp.png}
  \includegraphics[width=\hsize]{249-poststamp.png}
  }
  \parbox{.72\linewidth}{
  \includegraphics[width=\hsize]{528-large.png}}
  \parbox{.243\linewidth}{\includegraphics[width=\hsize]{317-poststamp.png}
  \includegraphics[width=\hsize]{333-poststamp.png}
  \includegraphics[width=\hsize]{704-poststamp.png}
  }
  \caption{AARTFAAC image (integrated over all observing bands) of the
    meteors at 2020-08-12 21:25:48 UTC (top left) and 2020-08-13 01:41:24 UTC (bottom left). The large-scale diffuse emission is the Galactic plane, the bottom left image shows some residuals of the subtracted sources Cassiopeia~A and Cygnus~A. The right column shows zoomed in meteors, from top to bottom corresponding to numbers 144, 235, 249, 317, 333, and 704 in the CAMS data set. In red, the trajectory as computed from the CAMS optical observations is
    overlaid.}
    \label{fig:camsmatch77}
\end{figure}

All-sky radio images were created for each of the individual AARTFAAC
subbands (a set of 4 channels). This analysis used a newer version of the
pipeline described in \citet{2021arXiv210315160S}. As part of the
pipeline, the brightest radio sources Cassiopeia\,A and Cygnus\,A were
subtracted, to minimize imaging artefacts across the all-sky
images. Typically, in subbands containing strong narrow-band
terrestrial radio emission reflected by meteor ionization trails, such
as emission from the BRAMS meteor forward scatter network transmitter in Belgium near 50\,MHz
\citep{2016EGUGA..1811624L}, the imaging failed, and images from these
subbands were rejected. Images from the remaining subbands were
averaged in radio frequency across subbands. This resulted in 10713
all-sky radio images at a 2\,s cadence, corresponding to 5\,h of
data. Example all-sky images obtained with AARTFAAC are shown in
Fig.\,\ref{fig:camsmatch77}. A full time-lapse of these images is available as \cite{zenodo_mp4_perseids}.

The CAMS BeNeLux low-light video network is part of the global 
(''Cameras for Allsky Meteor Surveillance'') network and uses close to
100 low-light video cameras spread over the BeNeLux to triangulate optical 
meteors +4 and brighter using methods 
described in \cite{2011Icar..216...40J}. Weather was clear during 
the night of 2020 August 12-13 and the trajectory and orbits of 720 
meteors were measured \citep{2020eMetN...5..400R}.

Figure\,\ref{fig:camsmatch77} shows a subset of these meteor
trajectories overlaid on the all-sky radio images from the perspective
of the ''super terp''. All all-sky radio
images within a few seconds of CAMS detections were inspected manually
for radio emission coincident with the reconstructed meteor
trajectory. This resulted in 204 meteors where radio emission was
coincident with a CAMS detection. Of these meteors, 59 had a
discernible trail in the radio images, the others showed up as
point-like in the radio images. There were also dozens of radio
meteors without a counterpart in the CAMS data, mostly due to missing
sky coverage of CAMS in Germany at the time, see Fig.\,\ref{fig:camsmap}.

For each radio meteor coincident with a CAMS detection we have
determined the begin and end point of the radio trail in the all-sky
radio images and the duration for which the meteor was visible in the
radio images.


\section{Results}
The coincidence of optical and radio trails demonstrates that meteors
are a source of radio emission, be it scattering or intrinsic emission.
That emission can last up to minutes; the longest radio train we have observed lasted
6.5\,min, coincident with a CAMS meteor detection with optical magnitude of
$-9$. In general, we find that brighter optical meteors result in
radio trains that remain visible for longer, as shown in
Fig.\,\ref{fig:magnitude}a.

The detected radio emission is integrated over a 2-second time interval. 
Based on the decay of brightness in subsequent 2-s intervals, the short trails 
are mostly emission from the meteor's persistent train, rather than from the meteor itself.

The optical trails of CAMS detected meteors start and end at higher
altitudes for brighter meteors, and while the radio trains show a
similar dependency (Fig.\,\ref{fig:magnitude}b and c), the radio
emission is first detected at lower altitudes
($h=101\pm4$\,km) compared to optical emission ($h=107\pm6$\,km). The
altitude at which the optical and radio emission ends is comparable
for optical ($h=94\pm6$\,km) and radio ($h=92\pm6$\,km).

\begin{figure}[!ht]
  \centering
  \includegraphics[width=\linewidth]{figures/magnitude}
  \caption{The dependency of the radio visibility and the altitude of
    the beginning and end of the meteor trail with optical
    brightness and entry speed. Note that the CAMS-derived peak optical brightness is less
    reliable below -5 magnitude due to video blooming.}
  \label{fig:magnitude}
\end{figure}

Some of the radio trains that were persistent for several minutes were
distorted by high altitude winds, much like optical persistent trains. 
Figure\,\ref{fig:time-evolution} shows an example of the spatial 
evolution of the radio train with time. The peak intensity of the 
radio meteor train decreased over time.

\begin{figure*}
  \centering
  \scalebox{1.04}{
  \includegraphics[width=0.45\linewidth]{figures/time_evolution}\hspace{5mm}
  \includegraphics[width=0.47\linewidth]{freq_evolution}}
  \caption{(Left:) The temporal and spatial evolution of a persistent
    radio meteor train, showing distortion by high altitude winds. The
    time is referenced to the beginning of the optical meteor
    trail. (Right:) Frequency evolution of a meteor train from a
    bright optical meteor. In channels where the emission from the
    meteor was very bright, the calibration has failed. That is the
    case in the bands where narrow-band terrestrial radio emission is
    reflected from the ionization train. Known transmitters operate at
    32.54\,MHz (TV transmitter), 49.97 and 49.99\,MHz (reflection of
    the BRAMS meteor radars; \citealt{2016EGUGA..1811624L}) and
    50.00\,MHz (amateur radio transmissions).}
  \label{fig:time-evolution}
\end{figure*}

Some frequency channels show reflected narrow-band artificial radio sources, but all
channels contain a broadband radio component that is not likely from
artificial radio sources. For bright radio meteors, the 
signal-to-noise ratio was sufficient to
detect this broadband emission in almost all individual channels for 
frequencies below $\sim50$\,MHz. An example is shown
in Fig.\,\ref{fig:time-evolution}.

Further insight was obtained when repeating the AARTFAAC observations 
during the 2020 Geminid meteor and the 2021 Quadrantid meteor shower. 
Unfortunately, cloudy weather prevented simultaneous CAMS
video observations. During that campaign, the all-sky radio images 
showed significantly fewer radio meteors compared to the observations of the
Perseids. 


\section{Discussion}
We can exclude an important contribution from forward scattered
narrow band terrestrial radio emission, which is frequency resolved in
our observations. Broadband emission would require a source of very broad band 
(30--50\,MHz) terrestrial radio emission, which we consider unlikely.

Of the proposed intrinsic emission scenarios, chemically produced 
suprathermal electrons \citep{2020JGRA..12528053O} could cause lingering
radiation over minutes timescale if those electrons are captured by
atoms and molecules that slowly diffuse into the train. Other proposed mechanisms such as
Langmuir waves \citep{2015JGRA..120.9916O}, free-free emission
\citep{2018Ge&Ae..58..693F} and transition radiation
\citep{2020JGRA..12528053O}, would be expected strong in the meteor
itself, but we see the radiation increase in intensity following the
meteor head before fading.

There are several possible reasons why fewer meteor trails were detected during
the Geminids and Quadrantids. Geminids and Quadrantids are known to reach
peak brightness at lower altitudes than the Perseids and show generally weaker
optical persistent trains  \citep{2006MeteorShowers}. In our case, the lower 
entry velocities of the
Geminids and Quadrantids in comparison to the Perseids is perhaps not to blame.
We note that several slower sporadic meteors of similar brightness
were observed in radio during the Perseids meteor shower, see Fig.\,\ref{fig:magnitude}a. 

The atmospheric conditions may be important. During the Perseids, the
summer weather provided hot and dry conditions, while the
Geminids and Quadrantids were observed during the rainy winter season.

Finally, it is possible that the local sidereal time was important. 
The scattering of bright celestial radio sources, as suggested in 
\cite{2015JGRA..120.9916O} was rejected in that paper with an argument involving the
fact the earlier LWA detections are very bright. Since we have detected
much weaker radio meteors, that argument does not hold here.

During the Perseids, the Galactic center and Galactic plane, which is the
source of most low-frequency radio emission, was above the LOFAR
horizon, while this was not the case during the Geminids and
Quadrantids. 


\section{Conclusions}
The AARTFAAC observations presented here show that persistent radio emission 
from meteors coincides with their optical
trajectories, and can be detected for meteors of magnitude $\sim0$ and
brighter. Aside from a narrow-band scattering of artificial radio sources, there is also a 
generally broadband emission detected for frequencies below
50\,MHz. The continuum radio emission is first detected at lower altitudes 
compared to the optical emission, while both radio and optical emission end 
at similar altitudes. For the brightest optical meteors,
we find that the radio emission can persist for several minutes. 
Persistent radio trains are affected by high altitude winds like
the optical persistent trains.

Analysis of the AARTFAAC observations is ongoing, so it is currently
not possible to distinguish between the suggested origins of the broadband 
radio emission. 

\vspace{5mm}
\noindent
\textbf{Acknowledgements:} We thank Albert van Duin and Edwin van Dijk
for the use of facilities at the Edgar Getreuer Observatory in
Burlage. We acknowledge interesting discussions with the researchers
in the BRAMS project. We thank all camera operators and other
contributors to CAMS-BeNeLux. PJ acknowledges support from NASA grant
80NSSC19K0563 (Solar System Workings). 

\vspace{5mm}
\noindent
\textbf{Supplemental materials:} Properties of the radio meteors observed
simultaneously with CAMS are available as \citep{zenodo_aartfaac_cams_csv}.


