# Solar radio emission as a disturbance of aeronautical radionavigation

Christophe Marqué, Karl-Ludwig Klein, Christian Monstein, Hermann Opgenoorth, Antti Pulkkinen, Stephan Buchert, Säm Krucker, Rudiger van Hoof, Peter Thulesen

► To cite this version:

Christophe Marqué, Karl-Ludwig Klein, Christian Monstein, Hermann Opgenoorth, Antti Pulkkinen, et al.. Solar radio emission as a disturbance of aeronautical radionavigation. Journal of Space Weather and Space Climate, 2018, 8, pp.A42. 10.1051/swsc/2018029. hal-01902224

HAL Id: hal-01902224

https://hal.sorbonne-universite.fr/hal-01902224

Submitted on 23 Oct 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Solar radio emission as a disturbance of aeronautical radionavigation

Christophe Marqué<fnref n="1" />*, Karl-Ludwig Klein<fnref n="2" />, Christian Monstein<fnref n="3" />, Hermann Opgenoorth<fnref n="4" />, Antti Pulkkinen<fnref n="5" />, Stephan Buchert<fnref n="4" />, Säm Krucker<fnref n="6" />, Rudiger Van Hoof<fnref n="7" />, and Peter Thulesen<fnref n="8" />

Received 28 September 2017 / Accepted 6 August 2018

**Abstract** – On November 4th, 2015 secondary air traffic control radar was strongly disturbed in Sweden and some other European countries. The disturbances occurred when the radar antennas were pointing at the Sun. In this paper, we show that the disturbances coincided with the time of peaks of an exceptionally strong (~10⁵ Solar Flux Units) solar radio burst in a relatively narrow frequency range around 1 GHz. This indicates that this radio burst is the most probable space weather candidate for explaining the radar disturbances. The dynamic radio spectrum shows that the high flux densities are not due to synchrotron emission of energetic electrons, but to coherent emission processes, which produce a large variety of rapidly varying short bursts (such as pulsations, fiber bursts, and zebra patterns). The radio burst occurs outside the impulsive phase of the associated flare, about 30 min after the soft X-ray peak, and it is temporarily associated with fast evolving activity occurring in strong solar magnetic fields. While the relationship with strong magnetic fields and the coherent spectral nature of the radio burst provide hints towards the physical processes which generate such disturbances, we have so far no means to forecast them. Well-calibrated monitoring instruments of whole Sun radio fluxes covering the UHF band could at least provide a real-time identification of the origin of such disturbances, which reports in the literature show to also affect GPS signal reception.

**Keywords:** Sun: radio radiation / Sun: flares / Sun: solar-terrestrial relations

## 1 Introduction

Solar eruptive events such as flares or coronal mass ejections (CMEs) are often accompanied by radio emission produced by supra thermal electrons propagating in open or closed magnetic structures of the solar corona. Solar radio bursts have been observed basically in the whole radio spectrum that current technology makes accessible, from a few hundreds of GHz down to a few kHz. For technical reasons, the VHF-UHF domain (~30–3000 MHz) has been monitored

since the early days of solar radio astronomy, and it appeared particularly important for the understanding of fundamental processes occurring during eruptive events (Pick and Vilmer, 2008). From the meter to decimeter range, characteristic burst signatures correspond to well identified physical processes, such as shock waves (type II bursts), electron beams streaming along open or elongated magnetic field lines (type III bursts) or electron populations trapped in eruptive flux ropes and post-flare loops (type IV bursts). Type IV bursts are often accompanied by long duration events observed at EUV or soft X-ray wavelengths, and coronal mass ejections. Discovered in the 1950s in the meter range (Boischot, 1958), they have been

<footnote n="1">¹ Solar-Terrestrial Center of Excellence – SIDC, Royal Observatory of Belgium, Avenue Circulaire 3, 1180 Brussels, Belgium</footnote>

<footnote n="2">² LESIA-UMR 8109 and Station de radioastronomie de Nançay – Observatoire de Paris, CNRS (further affiliations: PSL Res. Univ., Univ. P & M Curie and Paris-Diderot, University of Orléans, OSUC), 92190 Meudon, France</footnote>

<footnote n="3">³ Institute for Particle Physics and Astrophysics, ETH Zürich, Wolfgang-Pauli-Strasse 27, 8093 Zürich, Switzerland</footnote>

<footnote n="4">⁴ Swedish Institute of Space Physics, Box 537, 75121 Uppsala, Sweden</footnote>

<footnote n="5">⁵ Goddard Space Flight Center Greenbelt, MD, USA</footnote>

<footnote n="6">⁶ University of Applied Sciences and Arts Northwestern Switzerland, Windisch, Switzerland</footnote>

<footnote n="7">⁷ Belgoccontrol, Tervuursesteenweg 303, 1820 Steenokkerzeel, Belgium</footnote>

<footnote n="8">⁸ Air Greenland – Engineering Department, P.O. Box 1012, Nuuk, Greenland</footnote>

<footnote>*Corresponding author: christophe.marque@oma.be</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** False echoes (in green) observed in the direction of the Sun at a Belgian A/C radar station, located to the top right. The border with France and Luxembourg appears in blue. The scale of the map, in Nautical Miles is indicated at the top left.</figcaption></figure>

since then observed from the decimeter to the decimeter range. They often consist of spectral fine structures (broadband pulsations, zebra bursts, fiber bursts) superimposed on a continuum emission (for recent reviews see Chernov, 2006 or Nindos et al., 2008). Finally, these bursts often show high brightness temperature and flux density, which represents a challenge in terms of emission mechanisms (Cliver et al., 2011).

The potential impact of such intense radio emission on technologies relying on the use of the radio spectrum has been an important topic of research in the field of space weather (Bala et al., 2002). Several technologies have been impacted by solar radio bursts occurring in the VHF–UHF bands: Global Navigation Satellite Systems (GNSS), where the solar emission affects the carrier-to-noise ratio of ground based receivers in scientific or industrial applications (Klobuchar et al., 1999; Cerruti et al., 2006; Demyanov et al., 2012), and radar systems used in military applications for the airspace surveillance (Knipp et al., 2016). It is also worth remembering that one of the first published reports of radio emission coming from the Sun concerns the jamming of military radars during World War II (Hey, 1946).

In this paper, we report on an intense solar radio burst that occurred on November 4th, 2015, when air traffic radars were disrupted in Sweden, and milder disturbances were observed

elsewhere in Europe. A short account of the radar disturbances is presented in Section 2 to characterise the timing of the incidents. In Section 3 possible space weather causes are explored, including activity in the Earth's magnetosphere and at the Sun. The timing of the disturbances is shown to coincide with strong solar radio emission in the vicinity of 1000 MHz. The spectral characteristics of the radio burst are presented, together with context observations of the associated eruptive activity. The occurrence of such events and possible physical causes are discussed in Section 4, where recommendations for future observations are also formulated.

# 2 Disturbances of air traffic control radar on November 4th, 2015

## 2.1 Air traffic control

In the following, we outline the main aspects of Air Traffic Control (ATC) radar systems relevant to the present study. Most of the information comes from the Guidance Material on Comparison of Surveillance Technologies (GMST) document published by the International Civil Aviation Organization (ICAO, 2007).
Several radar systems are used to locate and gather information about aircraft in a given air space. Primary Surveillance Radars (PSR) provide basic information about azimuth and distance of a plane, by measuring the time difference between the emission of a pulse and the reception of the signal reflected on the frame of the aircraft. PSR are primarily set up at or near airports and generally operated worldwide either in L band (1215–1350 MHz) or S band (2700–3100 MHz).

Secondary Surveillance Radars (SSR) send coded queries to transponders aboard aircraft and get in return ancillary information on the plane: identification, barometric altitude of the aircraft and for some systems, selected technical parameters. They are “blind” to air planes having no functioning transponders. Two technologies are currently in use: mode A/C radars are an older generation of devices that send queries for name identification (mode A) and barometric altitudes (mode C) of the aircraft. These older radar systems present some limitation in terms of information being transmitted back by the aircraft and are essentially querying constantly all airplanes within their field of view. They are progressively superseded by mode-S radars, which can gather more information about the aircraft and interrogate specific planes only. Both generations of SSR operate in two bands: 1030 MHz for the uplink and 1090 MHz for the downlink; they have a range of about 200–250 Nautical Miles (370–463 km).

Other aeronautical services make use of the radio spectrum to help pilots operate the aircraft at critical times. For example, when a plane is about to land in difficult weather conditions (poor visibility), pilots may rely on an Instrument Landing System (ILS) to help them align correctly with the landing path and descent at the appropriate angle and speed towards the ground. ILS relies on VHF bands and provides via different subsystems the deviation in azimuth around the end point of the landing path (the localizer), and the position in elevation above the reference point (the gliding slope) in addition to other ancillary information (FAA, 2012).

## 2.2 Air traffic disturbances

On November 4th, 2015, several European air traffic authorities reported issues with operations of secondary air traffic radar systems at or near local sunset. In the following, we present problems affecting radars in Belgium and Sweden that were made available to us via a post event technical report analysis (Belgium), or that we compiled from several publicly available sources (Sweden). We also report on a possible ILS problem at the Thule Airbase in Greenland.

### 2.2.1 Disturbances in Belgium

Belgocontrol, the air traffic authority in Belgium, reported issues with a secondary A/C radar. False echoes, representing non-existing planes, were observed only in the direction of the Sun during two periods of time: 14:19:23–14:33:35 UT and 14:45:52–14:49:18 UT. Figure 1 shows the false plane positions aligned in the direction of the Sun observed at that time. No traffic interruption occurred due to this incident as the control software filtered out correctly the false information.

### 2.2.2 Disturbances in Sweden

Severe disruptions of the air traffic over the southern part of Sweden occurred on that day resulting in a de facto partial closure of the airspace and delayed arrivals and departures according to reports in the media (The Local, 2015). Based on publicly available information from the Swedish air traffic authority (Luftfartsverket; hereafter LFV), ATC secondary surveillance radar systems could not display proper information to the air traffic controllers, which lead the authorities to reduce aircraft movements for safety reasons (Luftfartsverket, 2015). More details were later made public by LFV during a workshop in November 2017, in Uppsala (A.J. Andersson, 2017; personal communication) indicating for this event that a simultaneous series of ATC radar disturbances occurred at several sites in Sweden starting around 14:19 UT. These disturbances were observed as ghost echoes in straight lines in the direction of the Sun and took place first between 14:19 and 14:34 UT, and then later on between 14:47 UT and 14:50 UT. The disturbances were more pronounced during the first time period. The magnitude of disturbances resulted in some stations being overloaded and triggering loss or delays of tracks.

Media reports and an early on press release pointed at a space weather event as the cause of the disruptions, evoking, confusingly a solar flare and a geomagnetic storm (The Local, 2015).

### 2.2.3 Other mentions of air traffic disturbances

Disturbances affecting SSR systems have been recently made public by air traffic authorities from Norway (Avinor Fly sikring) during the Birkeland Anniversary Conference 2017 conference, in Norway (A.D. Skjervold, 2017; personal com- munication<fnref n="1" />). Like in the Belgian and Swedish cases, ghost or false signals lining up in the direction of the Sun around 14:30 UT, were observed at several ATC radar stations in Norway without resulting in traffic perturbations. The same reference indicates that older generation radars (A/C types) were more “impacted”.

## 2.3 Disturbances in Greenland

On November 4th, 2015, an Air Greenland plane landing at Thule Airbase in Greenland at 14:49 UT experienced technical issues above 4000 feet altitude with a conflicting report between an ILS localizer (at 109.5 MHz) indicating a correct alignment with the runway and the autopilot being unable to hold on that same position information. The landing went without further complications. After the flight the ILS equipments were checked for malfunction but were cleared of any defect.

# 3 Solar and magnetospheric activity on November 4th, 2015

In the search of a candidate space weather event to explain the radar disturbances, we consider the timing of activity in the Earth magnetosphere and in the solar corona.

<footnote n="1">¹ https://www.mn.uio.no/english/about/news-and-events/events/The%20Birkeland%20Anniversary/andreas-d-skjervold_150617.pdf</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Overview of the November 4th 2015 event as observed in soft X-rays (GOES, bottom panel), at selected radio frequencies (Sagamore Hill station of RSTN, 2nd from bottom) and as dynamic spectrograms (HSRS spectrograph in Humain, Belgium; ORFEES spectrograph and Decameter Array, Nançay, France). Dark shading shows strong emission. Horizontal lines are from terrestrial emitters.</figcaption></figure>

Sustained geomagnetic activity is revealed by the Kp index on and around November 4th. Kp exceeds values of 5, the lower limit for the NOAA qualification² as a “minor” geomagnetic storm, on 3rd, 4th, and 7th of November. But Kp stays below 7, the threshold of the “strong” category. Following the NOAA classification 1700 minor geomagnetic storms occur during a solar activity cycle, and 200 strong ones. There is hence no exceptional geomagnetic disturbance that could be

a candidate space weather event explaining the air traffic radar problems on 4th November.

## 3.1 Solar eruptive activity

A solar flare (SOL2015-11-04T13:52) took place in NOAA AR 12443 while close to the central meridian. This region began its transit on the solar disk on October 29th as a complex sunspot group, with a McIntosh classification evolving from Ehc on October 30th to Fkc on November 4th. On the day of the flare, its magnetic type, according to the Mount-Wilson

<footnote n="1">¹ https://www.mn.uio.no/english/about/news-and-events/events/The%20Birkeland%20Anniversary/andreas-d-skjervold_150617.pdf</footnote>
classification was βδ. From October 29th till November 4th noon, an intense flaring activity took place with two strong B class flares, 55 C class flares, and one M class flare. None of these flares was accompanied by intense broadband radio emissions. Similarly, no significant noise storm activity was recorded at metric wavelength during that time frame.

SOL2015-11-04T13:52 was a long duration soft X-ray M class flare, also observed in Hz by the Kanzelhöhe Observatory as an 1B optical flare. It was accompanied by a halo CME first observed at 15:12 UT in the field of view of the SoHO/LASCO-C2 coronagraph (Brueckner et al., 1995). A streamer deflection suggests a CME-driven shock wave, which is indeed confirmed by metric and decametric type II burst emissions.

In the following we examine in more detail the soft X-ray and radio observations of the event.

### 3.1.1 Instrumentation

We make use of soft X-ray observations provided by the XRS instrument (Garcia, 1994) aboard the Geostationary Operational Environmental Satellites (GOES-13, and GOES-15), and EUV contextual images from the AIA instrument aboard the Solar Dynamics Observatory spacecraft (Lemen et al., 2012).

The radio spectra discussed in this paper were collected at three facilities in Belgium (Humain Radio Astronomy Station), France (Nançay Radio Astronomy Station), and Switzerland (Bleien and Zürich Observatories).

In Humain (Belgium), a 6-m dish radio telescope is used to monitor the solar activity on a daily basis. A dual-polarised log-periodic feed antenna covers the band 300–3000 MHz, but currently the setup is such that flux density measurements are performed only between 275 and 1495 MHz (in particular, no polarisation information is available). The Humain Solar Radio Spectrometer (HSRS) attached to the antenna is a Software Defined Radio receiver which is scanning the whole band in steps of 20 MHz. For each step, a Fast Fourier Transform spectrum is calculated and stacked with the others to produce a wide band spectrum with a constant frequency resolution of ~98 kHz and a time resolution of ~0.250 s.

The ORFEES (Observations Radio pour FEDOME et l'Étude des Éruptions Solaires) spectrograph setup in Nançay observes the flux density of the entire Sun between 100 and 1000 MHz, using a parabolic dish of diameter 5 m and a log-periodic dipole array. The spectral range is covered by five bands of about 200 MHz bandwidth. Each band is scanned within 0.1 s. The observations were complemented towards lower frequencies (20–80 MHz) by the Nançay Decameter Array (NDA; Lecacheux, 2000), a sensitive spectrograph connected to a field of 144 helicoidal antennas.

While routine solar radio observations were interrupted recently in Switzerland, the event under study has fortunately been observed by two instruments: one in Bleien, a Callisto spectrometer (Benz et al., 2005) connected to a cylindrical horn-antenna covering the band 1000–1300 MHz, and mounted on a 5-m dish telescope, and the other one in Zürich: a Callisto receiver connected to a cylindrical horn antenna at 1427 MHz, mounted on a 5-m dish telescope. In both cases, heterodyne down-converters were used to accommodate for

the original input frequency range of the Callisto receivers (45–870 MHz).

In addition to dynamic spectra, we also used fixed-frequency time histories observed by the Sagamore Hill station of the Radio Solar Telescope Network (RSTN; Kennewell, 2008) operated by the US Air Force.

### 3.1.2 Soft X-ray and radio emission

Figure 2 shows the time histories of the full Sun soft X-ray (bottom panel) and radio emission. The radio spectra cover a large band including the radar frequencies. The event looks unremarkable in GOES and in radio bands monitored by the RSTN Sagamore Hill station. The soft X-ray burst started at 13:35 UT and peaked at the M3.7 level at 13:52 UT. This burst is significant, but moderate: judging from the time since the beginning of the GOES soft X-ray observations, 31 such bursts or stronger are expected to occur per year. Moreover the burst had its most intense phase well before the radar incidents. An impulsive microwave burst was observed during the soft X-ray rise (second panel from bottom, 2695 to 8800 MHz), with faint counterparts at frequencies below 1000 MHz.

This is a typical impulsive burst, due to gyrosynchrotron emission of mildly relativistic electrons. Hard X-ray emission was observed up to 100 keV by the Fermi mission.<fnref n="3" />

While the soft X-rays and the microwaves decayed, a new rise to much higher flux density than during the impulsive phase occurred between about 14:05 and 15:10 UT at frequencies below 1500 MHz. The emission was broadband, with a low-frequency limit at or below 30 MHz. This limit is affected by the decrease of sensitivity of the Nançay Decameter Array antennas towards low frequencies, the emission probably extends to still lower frequencies. The large bandwidth and long duration show that this is a type IV radio burst.

The flux densities measured by RSTN are moderate, with peak values not exceeding 1400 Solar Flux Units (SFU)<fnref n="4" /> in any band monitored by this network but at 1415 MHz. The type IV burst is by no means a strong radio event, with the marked exception of the frequency range 800–1500 MHz seen in the HSRS, ORFEES and Callisto dynamic spectra. The peak flux density at 1415 MHz measured by RSTN is about 5800 SFU. Between 1000 and ~1200 MHz, two periods of particularly strong emission occur at 14:15–14:30 and 14:45–15:00 UT. This timing is in remarkable agreement with the strongest radar disturbances (~14:19–14:34 and 14:46–14:49 UT in Belgium), and points to the solar radio burst as the cause of these disturbances.

## 3.2 Calibration of the radio flux densities

As shown in Figure 2, the peak time of the light curves corresponds to the development of the type IV burst emission as observed by ORFEES and HSRS, and in particular to the strengthening of the high frequency part of the Type IV burst (above ~900 MHz). Even uncalibrated measurements provided by Nançay, Humain and Bleien observations hinted at a much stronger emission peak in that high frequency range. We performed a cross-calibration procedure based on relative

<footnote n="3">³ http://sprg.ssl.berkeley.edu/~tohban/browser/</footnote>

<footnote n="4">⁴ 1 SFU = $10^{-22}$ W m$^{-2}$ Hz$^{-1}$.</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Top: light curves in SFU of the November 4th 2015 event from the HSRS and ORFEES spectrographs. Bottom: light curves in relative dB from the HSRS spectrograph at several frequencies between 600 and 1500 MHz.</figcaption></figure>

measurements with respect to the quiet Sun level observed prior to that event. A spectrum of the quiet Sun level has first been calculated by fitting the Quiet Sun level observations provided by the Nobeyama Radio Polarimeters (Torii et al., 1979), the US-Air Force RSTN network and the F10.7 measurements at Penticton (Tapping, 2013).

The fitting was performed using a Gaussian Processes regression (Pedregosa et al., 2011) to cope with the large spread of RSTN values. For each frequency, we calculated a daily average and its associated standard deviation. For Nobeyama, we assumed a 1% standard deviation. The Quiet Sun flux

density in the mid of the radar band (1060 MHz) is found to be 75 SFU. This is the value used in our cross-calibration. As no polarisation measurements are available in the radar band, we assumed an unpolarised emission.

### 3.2.1 ORFEES calibration

The ORFEES spectrograph upper limit (1000 MHz) falls short of the radar band but is close enough to provide pertinent information on the flux density. The antenna was off-pointed by about half a beam at 1000 MHz since about 10:20 UT. The flux
<table><caption>**Table 1.** Summary of the magnitude of the radio burst (peak flux) at several frequencies.</caption><thead><tr><th>Frequency (MHz)</th><th>ORFEES (Nançay)</th><th>HSRS (Humain)</th><th>Callisto (Bleien)</th><th>Callisto Zürich</th><th>NOAA report</th></tr></thead><tbody><tr><td>610</td><td>820</td><td>1034</td><td></td><td></td><td>1000<sup>a</sup></td></tr><tr><td>1000</td><td>10<sup>5b</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>1060</td><td></td><td>1.57 × 10<sup>5c</sup></td><td></td><td></td><td></td></tr><tr><td>1000–1250</td><td></td><td></td><td>1.23 × 10<sup>5</sup></td><td></td><td></td></tr><tr><td>1415, 1427</td><td></td><td>5245<sup>d</sup></td><td></td><td>6310</td><td>5800<sup>e</sup></td></tr></tbody></table>

density was recalibrated at 1000 MHz using the following procedure: the flux densities measured in the morning, before the antenna was pointed at the Sun, are considered to represent the sum of the cold sky and contributions from the hardware (antenna, receiver, cables). Their average value was subtracted from the flux density time profile. It was then assumed that, after this subtraction, the flux density between 12:10 and 12:20 UT represented the quiet-Sun level. The Quiet Sun flux density inferred from the fit of the quiet Sun spectrum just mentioned was used to evaluate the calibrated flux density at 1000.11 MHz. As a result, we find that the 1000 MHz emission during the burst on November 4th, 2015, between 13:30 and 15:30 UT, reached a peak value of $10^5$ SFU (0.98 × $10^5$ SFU at 14:28:12 UT).

### 3.2.2 HSRS calibration

In Humain, the two radar bands 1030 and 1090 MHz are effectively observed by the instrument, but due to the proximity of an ATC radar (Saint Hubert Airport) about 21 km to the South East of the Humain station, the light curves in these two bands are heavily disturbed by the radar emission. For this event a light curve at 1060 MHz (central frequency) was calculated as the median, in spectral dimension, of observations between 1030 and 1090 MHz.

Before the raw data were scaled to SFU, estimates of the sky and ground contributions were made. Figure 3, top panel, shows the light curves recorded at Humain between 13:30 and 15:30 UT at 1060 MHz and two other bands which are registered in Belgium for Radioastronomy, 608–614 MHz and 1400–1427 MHz. Observations stopped shortly before 15:30 UT, about 45 min before local sunset, when the telescope went back automatically to its night time “parking position”. For all these bands, quiet sun flux densities were derived from the modelled quiet sun spectrum mentioned earlier. For comparison, light curves observed by ORFEES in the band 608–614 MHz and at 1000 MHz are also shown. For HSRS, the peak flux density at 1060 MHz (band 1030–1090 MHz) is approximately 157 kSFU at 14:28:18 UT, with two time intervals where the flux is above 100 kSFU: 14:26:12–14:29:33 and 14:47:26–14:47:51 UT. The peak flux values in the two other bands are 1034 and 5245 SFU at 608–614 and 1400–1427 MHz, which can be compared to the peak fluxes reported

by NOAA in the same bands (1000 and 5800 SFU at 610 and 1415 MHz respectively).

### 3.2.3 Swiss observations calibrations

The calibration of the Callisto observations in Bleien and Zürich followed the same procedure as for the two other radio data sets, and made use of the quiet sun spectrum aforementioned to derive a calibrated event spectrum for each instrument. The Bleien spectrum shows some signs of saturation as the upper limit of the receiver dynamic range is reached. Taking this into account, plus the horn frequency characteristics, the best estimate for the peak flux density value is $1.23 \times 10^5$ SFU over the frequency range 1000–1250 MHz. Peak values in the 1427 MHz band are rather similar in Humain, Zürich and RSTN observations (5245, 6300, 5800 SFU respectively).

### 3.2.4 Summary of radio calibrations

Table 1 gives an overview of the peak flux values estimated from the three radio observations used in the present study. Within the approximations and known sources of uncertainty, the peak of emission in the radar bands is higher than $10^5$ SFU and possibly as high as $\sim 1.6 \times 10^5$ SFU.

## 3.3 Spectrum and fine structure of the radio emission in the 600–1500 MHz range

The time history of the strong radio emission around 1000 MHz is shown in more detail in Figure 3, bottom panel, at selected frequencies across the spectrum. The time histories at the lower (600–620 MHz) and upper (1400–1420 and 1480–1500 MHz) limits of the plotted range are similar, with a rather smooth rise starting 14:15 UT and a maximum near 14:27 UT, which is followed by a decay with superposed similar fluctuations. The emission is clearly broadband. In the frequency range between these limits, in the range 700–1220 MHz, the emission is stronger. In addition, the initial rise is much steeper than at the lower and upper limit of the frequency range. After the steep initial rise, however, the time profiles of the strongest emission again resemble those at the lower and higher frequencies. They have a similar slow further rise, and peak at similar times. So the intense emission in the range 700–1220 MHz is superposed on a broadband continuum.

The steep rise is not an instrumental artefact, but a solar phenomenon, as seen by a comparison between the ORFEES and HSRS time profiles at respectively 1000 and 1060 MHz shown in the top panel of Figure 3. It is a second characteristic property of the radio burst spectrum around 1000 MHz, besides its high flux density.

More detailed spectra of the strong burst in the range 600–1500 MHz are shown in Figures 4–6. The spectra have been processed to increase the visibility of the fine structures. During this process the relative intensity from one frequency band to the next is lost. The spectrum in Figure 4 comprises broadband pulsations at the beginning of the first intense event, 14:15:30 UT, in the range 600–1150 MHz. Comparison of time profiles in fixed-frequency ranges plotted below the dynamic spectrum shows that it is indeed the pulsations which create the steep rise of the flux densities noted above. The 250 ms time evolution of

<footnote>a at 14:27 UT;</footnote>

<footnote>b at 14:28:12 UT;</footnote>

<footnote>c at 14:28:18 UT;</footnote>

<footnote>d at 14:27:25–30 UT;</footnote>

<footnote>e at 14:27 UT.</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Fine structure of the radio spectrum in the first minutes of the strong emission between 14:15 and 15:00 UT.</figcaption></figure>

HSRS probably does not resolve individual pulsations when they are too densely packed. We cannot rule out that the intensity observed during that flare arises from numerous spike bursts as observed during the 2006 December 06 event (Cliver et al., 2011).

The visual impression from the spectrum in Figure 4 is that the bright emission starts with individual pulses near 14:15:15, which thereafter become brighter and more densely packed. In the centre of the band the spectrum shows a more or less continuous emission after some time, but the succession of pulsations is clearly visible at the high-frequency limit throughout the spectral plot. So one may suspect that the intense emission is a superposition of many broadband pulsations. Other types of spectral fine structure become visible: fiber bursts (for instance, 14:15:40–14:16:00 between 960 and 1200 MHz), zebra pattern (e.g., 14:15:30–14:15:40, same frequency range). These types of structures are observed over the whole event as shown in Figure 5. These fine structures are discussed in Isliker and Benz (1994), Jiřička et al. (2001), Aurass et al. (2003), Chernov (2006).

The dynamic spectrum during the maximum of the radio burst near 1000 MHz is shown in Figure 6. Again broadband pulsations are clearly discernible. The spectral fine structures identified above are all well known from the observations of

type IV bursts at frequencies of several hundreds of MHz. They are usually ascribed to micro instabilities of non thermal electron populations trapped in coronal magnetic fields. A sustained feature in the range of the radar frequencies are the broadband pulsations.

**3.4 The eruptive event: EUV and coronographic observations**

The impulsive phase of the flare between ~13:25 and 13:50 UT is associated with several signatures in the EUV range: expanding loops, a coronal wave and the appearance of a dimming region to the North of the active region. A partial halo CME going essentially southward appears in the LASCO C2 field of view at 14:24 UT. The gradual phase of the flare, which overlaps with the radio event, shows classical post flare loops in the core of the active region. The time of the radio event corresponds to a late phase in this post-flare loop period, with new loops forming to the North-West of the active region.

The penumbra and umbra of the trailing sunspot (see Fig. 7) are swept over by a late expanding flare ribbon during at least the duration of the strongest and sudden radio flux increase near 1000 MHz. The radio sudden brightening starts
<figure><img src="image_6.png" /><figcaption>Fig. 5. Fine structures consisting of fiber bursts and zebra pattern at different stages of the event as observed by the Callisto spectrometer in Bleien.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 6. Dynamic radio spectrum observed with HSRS around the peak of the strong emission. Note that the flux density scale was set independently for the individual frequencies. The apparent maximum near 1400 MHz is an artifact of this procedure. The horizontal line at 1090 MHz, is actually the downlink part of a secondary ATC radar transmission close the station of Humain.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** AIA observations at 1600 Å at different times during the flaring event. The two panels to the bottom show a bright ribbon sweeping through the umbra of the sunspot circled in red, close to the time of the sharp increase of the radio flux density.</figcaption></figure>

<table><caption>**Table 2.** Parameters for the calculation of radar sensitivity to interferencesª.</caption><thead><tr><td>Parameter</td><td>Description</td><td>Value</td></tr></thead><tbody><tr><td>G</td><td>Maximum antenna gain</td><td>27 dBi</td></tr><tr><td>PT</td><td>Pointing in elevation</td><td>~6–8°</td></tr><tr><td>BW</td><td>Beam width</td><td>10° at 3 dB</td></tr><tr><td>S<sub>min</sub></td><td>Minimum detectable signal at antenna output</td><td>-103 dBm</td></tr><tr><td>CL</td><td>Cable loss</td><td>2 dB</td></tr><tr><td>RBW</td><td>Typical receiver bandwidth</td><td>5.5 MHz</td></tr><tr><td>S/I</td><td>Signal to interference ratio</td><td>12 dB</td></tr><tr><td>MG</td><td>Aeronautical safety margin and multiple technology allowance</td><td>12 dB</td></tr></tbody></table>

slightly earlier (14:15:30 UT) than the beginning of this flare ribbon expansion (14:17 UT). The penumbra appears profoundly modified by this intrusion of the flare ribbon, at least up to the end of the radio event. Earlier during the event, the first and weaker radio flux increase (from 13:40 till ~13:49 UT) corresponds also to a flare ribbon expanding into the same sunspot umbra.

## 4 Discussion and conclusion

The radar disturbances, which in Belgium, Sweden and Norway consisted of false echoes received when the antennas looked sunwards, had no coincident geomagnetic activity

enhancement. But the striking temporal association between the two periods of strongest disturbances and the two peaks of extreme solar radio emission in a band covering the radar frequencies points at the Sun as the origin of the perturbations.

### 4.1 Impact of the solar burst on radar systems and radionavigation

Detailed technical information concerning secondary radars is not publicly available. However, compatibility studies for the coexistence of different services in L band mention generic information about radar susceptibilities to interferences. Table 2 summarises basic facts extracted from Eurocontrol compatibility studies (Frequentis, 2007; Micallef, 2009) and private communications (B. Collard, DSNA & L. Lapène, DGAC).

<footnote>ª We follow here the methodology by Frequentis (2007), where a conservative interference threshold at receiver input is calculated as: $S_{min} + G - CL - (S/I) - MG$, which amounts to -102 dBm.</footnote>
<table><caption>**Table 3.** Peak flux densities in SFU for the strongest radio bursts since 2000 (peak flux greater than 50 000 SFU at 1415 MHz) tabulated by NOAA<sup>i</sup> (RSTN network at 1415 MHz) and by the Nobeyama Observatory<sup>ii</sup> (1000 and 2000 MHz).</caption><thead><tr><th>Date</th><th>Flux at 1000 MHz</th><th>Flux at 1415 MHz</th><th>Flux at 2000 MHz</th></tr></thead><tbody><tr><td>2001 April 15</td><td>N/A</td><td>54 000</td><td>N/A</td></tr><tr><td>2002 April 21</td><td>150 000</td><td>110 000</td><td>9000</td></tr><tr><td>2006 December 06</td><td>N/A</td><td>139 000<sup>a,b</sup></td><td>N/A</td></tr><tr><td>2006 December 13</td><td>440 000</td><td>130 000<sup>a</sup></td><td>302 000</td></tr><tr><td>2006 December 14</td><td>N/A</td><td>55 600</td><td>N/A</td></tr><tr><td>2011 February 15</td><td>46 000</td><td>54 000</td><td>1500</td></tr><tr><td>2011 September 24</td><td>N/A</td><td>110 000</td><td>N/A</td></tr><tr><td>2012 March 05</td><td>502 000</td><td>20 000<sup>c</sup></td><td>19 000</td></tr></tbody></table>

From these elements it appears that a conservative interfer-
ence level threshold at a radar receiver input is about
-102 dBm (see Table 2 for details). This can be compared
to the quiet sun flux density and the peak flux density during
the event. At 1090 MHz, we considered a quiet sun flux density
of 75 SFU. The power at receiver level is therefore:

$$
\begin{align*}
P_{\text{rec}} &= 10 \log (75 \times 10^{-22} \text{ W m}^{-2} \text{ Hz}^{-1} \times A_{\text{eff}} \times \text{RBW}) - \text{CL} \\
&\approx -101 \text{ dBm}
\end{align*}
$$

with the effective area:

$$
A_{\text{eff}} = \frac{G\lambda^2}{4\pi} \approx 3 \, \text{m}^2.
$$

This power is close to the threshold of interference, but is
not detrimental in normal operation. During the peak of the
radio event, for which flux densities above 100 kSFU were
observed, the power at receiver level amounted to at least
$P_{\text{rec}} \approx -68$ dBm, when radars pointed at the Sun. This is
about 34 dB above the interference threshold. We can therefore
conclude that radars, for which the Sun was between 2 and 12°
elevation, could easily pick up the radio burst, at power levels
several orders of magnitude above the interference threshold.

We also reported an issue (see Sect. 2.3) with ILS equip-
ments operating at 109 MHz for a plane landing at Thule Air-
base at the end of the radio event (14:49 UT). Post-event
discussions with technicians from the company at Thule Air
Base indicate that this was unusual and considered as “first
time observations”. The aircraft suffering the incident had its
ILS equipment checked afterwards for any defect. At the time
of the event, and at the altitude of the plane when the problem
was noticed (above 4000 feet), the Sun was 2.9 degrees below
the horizon. In the VHF band, dynamic spectra of the radio
emission (Fig. 2, top panel) show a low frequency extension
of the type IV burst, with occasional type III bursts. One was
in progress at 14:49 UT and was approximately 10–30 times
stronger than the background as measured by the spectra from
the Sagamore Hill and San Vito RSTN stations.

Refraction is needed in order to make the burst visible to
the airplane. At these low frequencies, strong refraction in

the ionosphere due to atmospheric gravity waves is frequently
observed, and it goes along with focussing effects that are read-
ily visible in solar radio spectra (Mercier et al., 1989). Raytrac-
ing calculations performed for Air Greenland for that event
seem to corroborate this conclusion (Per Høeg, 2016, DTU
Space, private communication). We conclude that the distur-
bance at 109 MHz can in principle be related to solar radio
emission. However, Figure 3 shows that the burst around
1 GHz, which we identified as the source of the radar distur-
bances, clearly stands out above the broad band radio event
at frequencies above 600 MHz.

## 4.2 Mechanism of the radio emission

Unlike the VHF burst, the radio emission around
1000 MHz is about 2000 times stronger than the Quiet Sun.
Explaining such a magnitude remains challenging. The bright-
est radio emission in the vicinity of the radar frequencies on
November 4th, 2015 is a collection of spectral fine structures,
especially broadband pulsations. The type IV burst has an over-
all evolution that is similar over a much more extended spectral
range than the bright emission. The fine structure with its sud-
den onset is a distinctive feature at those frequencies where the
burst is bright. In addition, fiber bursts and zebra patterns are
observed during some time intervals. They are ascribed to
the interaction between upper hybrid waves and whistlers,
which are generated in microinstabilities of trapped non ther-
mal electron populations (Zlotnik, 2013). In the corona, where
the plasma frequency is usually believed to exceed the cyclo-
tron frequency, such instabilities create plasma waves, which
then couple to radio waves that can escape.

An alternative process that creates directly escaping radio
waves, and is known as the mechanism of very intense radio
bursts in planetary magnetospheres, is the electron cyclotron
maser. This mechanism is usually expected to not operate in
the solar corona, because of the high electron plasma fre-
quency. The exceptional occurrence of the strong radio burst
near 1000 MHz could of course be explained by an exceptional
situation of plasma parameters in the corona, so that electron
cyclotron maser emission might arise. The argument was put
forward by Régnier (2015) and Cliver et al. (2011). However,
the identification of well-known fine structures, like fiber bursts

<footnote>i ftp://ftp.swpc.noaa.gov/pub/warehouse/</footnote>

<footnote>ii http://solar.nro.nao.ac.jp/norp/index.html</footnote>

<footnote>a Saturation limit;</footnote>

<footnote>b Cliver et al. (2011) report for that event a peak flux density of ~106 SFU from OVSA observations between 1 and 1.6 GHz;</footnote>

<footnote>c End of observations at peak flux; probably underestimated.</footnote>
and zebra patterns, that are observed in type IV bursts, but usu-
ally at lower frequencies where it is still more unlikely that the
cyclotron frequency exceed the plasma frequency, makes such
an exceptional situation improbable. Instead, it rather argues for
coherent plasma emission.

The sudden rise in flux density observed first around 13:40 UT and especially later around 14:15 UT corresponds approximately in time to the expansion of one of the flare ribbons through the umbra of one of the sunspots of the active region (13:40 and 14:17 UT). This behaviour is observed in a surprisingly small number of M and X class flares. In a survey of 588 M and X class flares observed by the TRACE and Hinode/SOT telescopes, Leping and Zhang (2009) have noticed only 20 events displaying a flare ribbon sweeping across the umbra of one of the two main sunspots of the active region. However, only two of their 20 events are associated with intense solar radio bursts (flux density above 100 000 SFU) and both occurred in the same region (NOAA AR 10930, on December 6th and 13th, 2006). We checked a more recent event, which took place in September 24th, 2011, with a peak flux of 110 000 SFU at 1415 MHz. It also shows an abrupt intensification of the flux density approximately as a flare ribbon sweeps across the umbra of one of the sunspots in the flaring active region. One physical implication could be that particularly strong magnetic fields are involved during the bursts of interest, which in turn could be a new support for the interpretation by a cyclotron maser mechanism. However, at this stage, one cannot draw any firm conclusion linking these two phenomena.

## 4.3 The occurrence frequency of strong solar radio bursts near 1000 MHz

A question that arises from this analysis is how frequently
such intense radio emission can trigger noticeable or damaging
consequences on human technologies. Nita et al. (2002) anal-
ysed the occurrence of strong radio bursts over nearly 40 years,
deriving the distribution function of burst occurrence depend-
ing on frequency, flux and phase of the solar cycle. Following
their analysis, events of magnitude greater than 10⁵ SFU in the
band 1–1.7 GHz occur every ~8 months during solar maxi-
mum, every ~38 months during solar minimum and every
~23 months over the 40 years of the survey. That does not nec-
essarily mean that a noticeable effect will be observed at these
frequencies. At a given location on Earth, air traffic radar sys-
tems will be most affected at or near Sun rise or Sun set. If we
suppose that radar systems are vulnerable within a 2 h window
at dawn or dusk, the occurrence rate will be multiplied by ~6:
every ~48 months, every ~228 months (19 years) and every
~138 months (11.5 years) for solar maximum, solar minimum
and the 40 years period. We note that the 2015 annual report of
the Swedish air traffic authority does mention that ATC distur-
bances linked to solar radio bursts have already occurred in
Sweden in 1999 and 2003, but with milder effects on the oper-
ations (Luftfartsverket, 2015).

A survey of strong solar radio bursts observed since the
beginning of the millennium (see Table 3), made from lists
of events from NOAA and Nobeyama Observatory, reveals at
least 8 events with peak flux density greater than 50 000
SFU, and two events peaking between 400 000 and 500 000

SFU at 1000 MHz. It is worth noticing that one of these two
events occurred close to the minimum of cycle 23–24.

The electron plasma frequency in flaring active regions has a range of values, and bursts like the one on November 4th, 2015 can occur at neighboring frequencies, especially in a range comprising GNSS frequencies. GNSS disturbances were indeed observed on November 4th<fnref n="5" />, but to our knowledge, without consequences on related services. If air traffic radars are affected when the Sun is low above the horizon, ground based GNSS systems will be more impacted when the Sun is close to the local zenith, which corresponds to the maximum of the gain of a typical GNSS antenna (Carrano et al., 2009). Services like Air Traffic management, which rely both on radars and GNSS systems, like the WAAS (Cerruti et al., 2008) in the USA, or EGNOS in Europe could therefore be affected at any time during the day.

Table 3 shows also, like for the present event, that peak flux
density can differ by several orders of magnitude in frequency
bands that are relatively close (1000 and 1415 MHz). As a mat-
ter of facts, the event investigated in this study was rather unre-
markable in the radio bands typically used for space weather
monitoring (RSTN) and its magnitude and potential impact
went unnoticed. A monitoring at 1415 MHz solely is likely,
in the future, to miss potentially problematic radio bursts for
air traffic radars and this alone would justify to establish
world-wide, real time solar monitoring facilities covering the
whole UHF band. At least, bands allocated to important ser-
vices used in air traffic activities should be closely monitored,
but comprehensive studies about the impact of solar radio
bursts on air traffic radio services will only be possible if more
technical and operational details are made available.

*Acknowledgements.* We thank Kristoffer Leer DTU Space, Lyngby, Denmark for making all personal contacts between several parties involved in the Greenland event. KLK acknowledges helpful discussions with B. Collard, B. de Courville, L. Lapène, J.-Y. Prado, R. Rosso, and B. Roturier. Research at Paris Observatory was supported by the Agence Nationale pour la Recherche (ANR/ASTRID, DGA) project Outils radioastronomiques pour la météorologie de l'espace (ORME, contract No. ANR-14-ASTR-0027). The construction of the CALLISTO spectrometer was financed by the Swiss National Science Foundation SNSF (Grant Nr. 20-67995.02). The Belgian solar radio instrumentation is funded by the Solar Terrestrial Center of Excellence. The editor thanks Bryn Jones, Dale E. Gary and an anonymous referee for their assistance in evaluating this paper.

References

Aurass H, Klein K-L, Zlotnik EY, Zaitsev VV. 2003. Solar type IV burst spectral fine structures. I. Observations. *A&A* **410**: 1001–1010.

Bala B, Lanzerotti LJ, Gary DE, Thomson DJ. 2002. Noise in wireless systems produced by solar radio bursts. *Radio Sci* **37(2)**: 1018.

Benz AO, Monstein C, Meyer H. 2005. Callisto – a new concept for solar radio spectrometers. *Sol Phys* **226**: 143–151.

<footnote n="5">⁵ see http://gnss.be/Atmospheric_Maps/srb_event.php?date=2015-11-04</footnote>
Boischot A. 1958. Étude du rayonnement radioélectrique solaire sur 169 MHz à l'aide d'un grand interféromètre à réseau. *Ann Astrophys* **21**: 273.

Brueckner GE, Howard RA, Koomen MJ, Korendyke CM, Michels DJ, et al. 1995. The large angle spectroscopic coronagraph (LASCO). *Sol Phys* **162**: 357–402.

Carrano CS, Bridgwood CT, Groves KM. 2009. Impacts of the December 2006 solar radio bursts on the performance of GPS. *Radio Sci* **44**.

Cerruti AP, Kintner PM, Gary DE, Lanzerotti LJ, de Paula ER, Vo HB. 2006. Observed solar radio burst effects on GPS/Wide Area Augmentation System carrier-to-noise ratio. *Space Weather* **4**: 10006.

Cerruti AP, Kintner PM, Gary DE, Mannucci AJ, Meyer RF, Doherty P, Coster AJ. 2008. Effect of intense December 2006 solar radio bursts on GPS receivers. *Space Weather* **6**: S10D07.

Chernov GP. 2006. Solar radio bursts with drifting stripes in emission and absorption. *Space Sci Rev* **127**: 195–326.

Cliver EW, White SM, Balasubramaniam KS. 2011. The solar decimetric spike burst of 2006 December 6: Possible evidence for field-aligned potential drops in post-eruption loops. *Astrophys J* **743**: 145.

Demyanov VV, Afraimovich EL, Jin S. 2012. An evaluation of potential solar radio emission power threat on GPS and GLONASS performance. *GPS Solut* **16**: 411.

FAA. 2012. *Instrument flying handbook*. Tech. Rep., US Department of Transportation – Federal Aviation Administration, https://www.faa.gov/regulations_policies/handbooks_manuals/aviation/media/faa-h-8083-15b.pdf.

Frequentis AG. 2007. *B-AMC interference analysis and spectrum requirements*. Tech. Rep., Eurocontrol, https://www.eurocontrol.int/sites/default/files/article/content/documents/communications/22102007-b-amc-project-deliverable-d4-v11.pdf

Garcia HA. 1994. Temperature and emission measure from GOES soft X-ray measurements. *Sol Phys* **154**: 275–308.

Hey JS. 1946. Solar radiations in the 4–6 metre radio wave-length band. *Nature* **157**: 47–48.

ICAO. 2007. *Guidance material on comparison of surveillance technologies (GMST)*. Tech. Rep., International Civil Aviation Organization, http://www.icao.int/RO_APAC/Documents/edocs/cns/gmsttechnology.pdf

Isliker H, Benz AO. 1994. Catalogue of 1–3 GHz solar flare radio emission. *A&AS* **104**: 145–160.

Jiřička K, Karlický M, Mészárosová H, Snížek V. 2001. Global statistics of 0.8–2.0 GHz radio bursts and fine structures observed during 1992–2000 by the Ondřejov radiospectrograph. *A&A* **375**: 243–250.

Kennewell JA. 2008. *RSTN Solar Radio Telescopes (Discrete Frequency) and Data*. URL: http://www.deepsouthern skies.com/LSO/RSTN.pdf.

Klobuchar JA, Kunches JM, VanDierendonck AJ. 1999. Eye on the ionosphere: potential solar radio burst effects on GPS signal to noise. *GPS Solut*, 3(2), 69–71.

Knipp DJ, Ramsay AC, Beard ED, Boright AL, Cade WB, et al. 2016. The May 1967 great storm and radio disruption event: Extreme space weather and extraordinary responses. *Space Weather* **14**: 614–633.

Lecacheux A. 2000. The Nançay Decameter Array: a useful step towards giant, new generation radio telescopes for long wavelength radio astronomy. Geophysical Monograph Series, Vol. 119, American Geophysical Union, Washington DC, p. 321.

Lemen JR, Title AM, Akin DJ, Boerner PF, Chou C, et al. 2012. The Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO). *Sol Phys* **275**: 17–40.

Leping L, Zhang J. 2009. Statistics of flares sweeping across sunspots. *Astrophys J Lett* **706**: L17–L21.

Luftfartsverket. 2015. Annual Report, 2015. https://www.lfv.se/glob-alassets/nyheter/nyheter-2016/eng_lfv2015_lores.pdf.

Mercier C, Genova F, Aubier MG. 1989. Radio observations of atmospheric gravity waves. *Ann Geophys* **7**: 195–202.

Micallef J. 2009. Compatibility criteria and interference scenarios for SSR systems. Tech. rep., Eurocontrol, https://www.eurocontrol.int/sites/default/files/article/content/documents/communications/24082009-lcis-c3-criteria-and-tests-v10.pdf.

Nindos A, Aurass H, Klein K-L, Trottet G. 2008. Radio emission of flares and coronal mass ejections. Invited review. *Sol Phys* **253**: 3–41.

Nita GM, Gary DE, Lanzerotti LJ, Thomson DJ. 2002. The peak flux distribution of solar radio bursts. *Astrophys J* **570**: 423–438.

Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, et al. 2011. Scikit-learn: machine learning in Python. *J Mach Learn Res* **12**: 2825–2830.

Pick M, Vilmer N. 2008. Sixty-five years of solar radioastronomy: flares, coronal mass ejections and Sun Earth connection. *A&ARv* **16**: 1–153.

Régnier S. 2015. A new approach to the maser emission in the solar corona. *A&A* **581**: A9.

Tapping KF. 2013. The 10.7 cm solar radio flux (F10.7). *Space Weather* **11**: 394–406.

The Local. 2015. *Solar storm grounds Swedish air traffic*. URL: https://www.thelocal.se/20151104/solar-storm-grounds-swedish-air-traffic.

Torii C, Tsukiji Y, Kobayashi S, Yoshimi N, Tanaka H, Enome S. 1979. Full-automatic radiopolarimeters for solar patrol at microwave frequencies. *Proceedings of the Research Institute of Atmospherics, Nagoya University* **26**: 129–132.

Zlotnik EY. 2013. Instability of electrons trapped by the coronal magnetic field and its evidence in the fine structure (zebra pattern) of solar radio spectra. *Sol Phys* **284**: 579–588.

**Cite this article as:** Marqué C, Klein K, Monstein C, Opgenoorth H, Pulkkinen A, et al. 2018. Solar radio emission as a disturbance of aeronautical radionavigation. *J. Space Weather Space Clim.* **8**, A42.