Issue |
J. Space Weather Space Clim.
Volume 10, 2020
Topical Issue - Scientific Advances from the European Commission H2020 projects on Space Weather
|
|
---|---|---|
Article Number | 10 | |
Number of page(s) | 16 | |
DOI | https://doi.org/10.1051/swsc/2020010 | |
Published online | 20 March 2020 |
Research Article
A LOFAR observation of ionospheric scintillation from two simultaneous travelling ionospheric disturbances
1
ASTRON – The Netherlands Institute for Radio Astronomy, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands
2
Department of Electronic and Electrical Engineering, University of Bath, Claverton Down, BA2 7AY Bath, UK
3
School of Science and Technology, Nottingham Trent University, Clifton Lane, NG11 8NS Nottingham, UK
4
Space Environment and Radio Engineering, School of Engineering, The University of Birmingham, Edgbaston, B15 2TT Birmingham, UK
5
Space Research Centre, Polish Academy of Sciences, Bartycka 18A, 00-716 Warsaw, Poland
6
Space Radio-Diagnostics Research Centre, University of Warmia and Mazury, ul. Romana Prawocheskiego 9, 10-719 Olsztyn, Poland
7
Technische Universität Berlin, Institut für Geodäsie und Geoinformationstechnik, Fakultät VI, Sekr. H 12, Hauptgebäude Raum H 5121, Straße des 17. Juni 135, 10623 Berlin, Germany
8
GFZ German Research Centre for Geosciences, Telegrafenberg, 14473 Potsdam, Germany
9
Shell Technology Center, 562149 Bangalore, India
10
Science and Technology B.V., 2616 LR Delft, The Netherlands
11
RAL Space, UKRI STFC, Rutherford Appleton Laboratory, Harwell Campus, OX11 0QX Oxfordshire, UK
12
Mt Stromlo Observatory, Research School of Astronomy and Astrophysics, Australian National University, Cotter Road, ACT 2611 Weston Creek, Australia
13
Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, 85748 Garching, Germany
14
Hamburger Sternwarte, Universität Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany
15
Thüringer Landessternwarte, Sternwarte 4, 07778 Tautenburg, Germany
16
Jodrell Bank Centre for Astrophysics (JBCA), Department of Physics & Astronomy, University of Manchester, Alan Turing Building, Oxford Road, M139PL Manchester, UK
17
Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
18
LPC2E – Université d’Orléans/CNRS, 45071 Orléans Cedex 2, France
19
Station de Radioastronomie de Nançay, Observatoire de Paris, PSL Research University, CNRS, Univ. Orléans, OSUC, 18330 Nançay, France
20
Department of Astrophysics/IMAPP, Radboud University, PO Box 9010, 6500 GL Nijmegen, The Netherlands
21
Nikhef, Science Park 105, 1098 XG Amsterdam, The Netherlands
22
Vrije Universiteit Brussel, Astronomy and Astrophysics Research Group, Pleinlaan 2, 1050 Brussel, Belgium
23
Kapteyn Astronomical Institute, University of Groningen, PO Box 800, 9700 AV Groningen, The Netherlands
24
Leibniz-Institut für Astrophysik Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany
25
ECAP, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erwin-Rommel-Str. 1, 91054 Erlangen, Germany
26
DESY, Platanenallee 6, 15738 Zeuthen, Germany
27
CIT, Rijksuniversiteit Groningen, Nettelbosje 1, 9747 AJ Groningen, The Netherlands
28
Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
29
Anton Pannekoek Institute, University of Amsterdam, Postbus 94249, 1090 GE Amsterdam, The Netherlands
30
Fakultät für Physik, Universität Bielefeld, Postfach 100131, 33501 Bielefeld, Germany
31
South African Radio Astronomy Observatory, 2 Fir Street, Black River Park, Observatory, 7925 Cape Town, South Africa
32
Department of Physics and Astronomy, University of the Western Cape, 7535 Cape Town, South Africa
33
Department of Physics and Electronics, Rhodes University, PO Box 94, 6140 Makhanda, South Africa
34
Jagiellonian University in Kraków, Astronomical Observatory, ul. Orla 171, 30-244 Kraków, Poland
35
Department of Physics, Khalifa University, PO Box 127788, Abu Dhabi, United Arab Emirates
36
Department of Space, Earth and Environment, Chalmers University of Technology, Onsala Space Observatory, 439 92 Onsala, Sweden
37
Joint Institute for VLBI-ERIC (JIVE), Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands
38
LESIA & USN, Observatoire de Paris, CNRS, PSL, SU/UP/UO, 92195 Meudon, France
* Corresponding author: fallows@astron.nl
Received:
2
December
2019
Accepted:
13
February
2020
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 Cassiopeia 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–80 MHz. 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 ~20–40 ms−1 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 ~110 ms−1 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.
Key words: Ionospheric scintillation / travelling ionospheric disturbances / instability mechanisms
© R.A Fallows et al., Published by EDP Sciences 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Radio waves from compact sources can be strongly affected by any ionised medium through which they pass. Refraction through large-scale density structures in the medium leads to strong lensing effects where the radio source appears, if imaged, to focus, de-focus and change shape as the density structures in the line of sight themselves move and change. Diffraction of the wavefront by small-scale density structures leads to variations building up in the intensity of the wavefront with distance from the scattering medium, due to interference between the scattered waves, an effect known as scintillation. Observations of all these effects thus contain a great deal of information on the medium through which the radio waves have passed, including the large-scale density, turbulence, and the movement of the medium across the line of sight. Since the second world war, a large number of studies have shown the effect of ionospheric density variations on radio signals, as reviewed by Aarons (1982), and this can lead to disruption for applications using Global Navigation Satellite Systems (GNSS, e.g., GPS), as thoroughly reviewed by Hapgood (2017). The Low-Frequency Array (LOFAR – van Haarlem et al., 2013) is Europe’s largest low-frequency radio telescope, operating across the frequency band 10–250 MHz, and with a dense array of stations in the Netherlands and, at the time of writing, 13 stations internationally from Ireland to Poland. It was conceived and designed for radio astronomy but, at these frequencies, the ionosphere can also have a strong effect on the radio astronomy measurement (de Gasperin et al., 2018). Ionospheric scintillation, which is rarely seen over the mid-latitudes on the high-frequency signals of GNSS, is seen almost continually in observations of strong natural radio sources by LOFAR.
The wide bandwidth available with LOFAR enables an easy and direct assessment of scattering conditions and how they change in a given observation, including whether scattering is weak or strong, or refractive effects dominate, and enables further information to be gleaned from delay-Doppler spectra (the 2-D FFT of a dynamic spectrum, termed variously as the “scattering function”, “generalised power spectrum”, or “secondary spectrum” – here we use the term “delay-Doppler” spectrum as this clearly describes what the spectrum shows). In observations of interstellar scintillation these spectra can exhibit discrete parabolic arcs which can be modelled to give information on the distance to the scattering “screen” giving rise to the scintillation and its velocity across the line of sight (Stinebring et al., 2001; Cordes et al., 2006). Broadband observations of ionospheric scintillation are not common, but such arcs have been observed using the Kilpisjärvi Atmospheric Imaging Receiver Array (KAIRA, McKay-Bukowski et al. (2014) – an independent station built using LOFAR hardware in arctic Finland) in a study by Fallows et al. (2014). Model spectra produced by Knepp & Nickisch (2009) have also illustrated parabolic arc structures, particularly in the case of scattering from a thin scattering screen.
The wide spatial distribution of LOFAR stations also enables scintillation conditions at these observing frequencies to be sampled over a large part of western Europe. A dense “core” of 24 stations, situated near Exloo in the north-east of the Netherlands, over an area with a diameter of ~3.5 km further provides a more detailed spatial view of the scintillation pattern in its field of view.
LOFAR thus enables detailed studies of ionospheric scintillation to be undertaken which can both reveal details which would be unavailable to discrete-frequency observations such as those taken using GNSS receivers, and act as a low-frequency complement to these observations to probe potentially different scattering scales.
A number of different phenomena can lead to scattering effects in radio wave propagation through the mid-latitude ionosphere: ionisation structures due to gradients in the spatial distribution of the plasma density can arise from a southward expansion of the auroral oval or from large- to small- scale travelling ionospheric disturbances (TIDs). Large-scale TIDs (LSTIDs) with wavelengths of about 200 km typically propagate southward after forming in the high-latitude ionosphere in response to magnetic disturbances (e.g., storms or sub-storms, Tsugawa et al., 2004). On the other hand, medium-scale TIDs (MSTIDs) seem to form in response to phenomena occurring in the neutral atmosphere triggering atmospheric gravity waves (AGWs), which then propagate upwards to generate TIDs at ionospheric heights (Kelley, 2009). The morphology of MSTIDs varies with local time, season, and magnetic longitude. Their propagation shows irregular patterns that vary on a case-by-case basis, although they commonly seem to propagate mainly equatorward during winter daytime and westward during summer night-time (Saito & Fukao, 1998; Hernández-Pajares et al., 2006, 2012; Tsugawa et al., 2007; Emardson et al., 2013). Smaller-scale ionisation gradients, likely associated with the Perkins instability (Kelley, 2009, 2011), can then form as a consequence of the presence of MSTIDs, potentially leading to scintillation at LOFAR frequencies.
In this paper, we perform an in-depth analysis of ionospheric scintillation seen in an observation of the strong natural radio source Cassiopeia A (Cas A) overnight on 18–19 August 2013. This observation was amongst the first of its kind taken with LOFAR and exhibited quite strong scattering effects across the 10–80 MHz band. The purpose of this paper is both technical and scientific: we first describe the observation itself, and then demonstrate several techniques to analyse LOFAR data and show how these can bring out the details of ionospheric structures. Finally, we use supporting data from GNSS and ionosondes to get a broader picture of conditions in the ionosphere at the time and how these give rise to the scintillation seen by LOFAR.
2 The LOFAR observation
LOFAR observed Cas A (Right Ascension 23h 23m 24s, Declination +58°48′54″) between 21:05 UT on 18 August 2013 and 04:05 UT on 19 August 2013, recording dynamic spectra from each individual station with a sampling time of 0.083 s over the band 2.24–97.55 MHz from each available station. The observing band was sampled with 7808 channels of 12.207 kHz each, but averaged over each successive 16-channel block to 488 subbands of 195.3125 kHz for the analyses described in this paper. At the time of observation the available stations were the 24 stations of the LOFAR “core”, 13 “remote” stations across the north-east of the Netherlands, and the international stations at Effelsburg, Unterweilenbach, Tautenburg, Potsdam, and Jülich (Germany), Nançay (France), Onsala (Sweden), and Chilbolton (UK). The reader is referred to van Haarlem et al. (2013) for full details of the LOFAR receiving system. The raw data for this observation can be obtained from the LOFAR long-term archive (https://lta.lofar.eu); observation ID L169059 under project “IPS”.
We first illustrate the data in a more traditional sense. Figure 1 shows time series’ at three discrete observing frequencies of the data taken by LOFAR station CS002, at the centre of the core, and their associated power spectra. The power spectra show a fairly typical shape for intensity scintillation: an initial flat section at the lowest spectral frequencies represents scattering from larger-scale density structures which are close enough to the observer that the scattered waves have not had the space to fully interfere to develop a full intensity scintillation pattern; the turnover (often termed the “Fresnel Knee”) indicates the largest density scales for which the intensity scintillation pattern has fully formed; this is followed by a power-law in the spectra illustrating the cascade from larger to smaller density scales, which is cut off in these spectra by white noise due to the receiving system (the flat section covering high spectral frequencies).
Fig. 1 (a) Time series of intensity received at three discrete frequencies by LOFAR station CS002 during the observation of Cas A on 18–19 August 2013, plus, (b) and (c) power spectra of two 10-min periods within these time series’. |
However, the advantage of observing a natural radio source with LOFAR is that full dynamic spectra can be produced covering the full observed band. Dynamic spectra of data taken by LOFAR station CS002 are presented in Figure 2, which includes a dynamic spectrum of the full observation, alongside more detailed views of three different single hours of the observation to illustrate the range of scattering conditions seen. The strength of the scattering can be seen much more clearly in this view, compared to time series’ from discrete observing frequencies. In general, scattering appears weak in this observation at the highest observing frequencies (where intensity remains highly correlated across the observing band) with a transition to strong scattering conditions as the observing frequency decreases. The frequency range displayed in these dynamic spectra is restricted to exclude the radio-frequency interference (RFI) which dominates below about 20 MHz and a fade in signal strength at the higher frequencies due to the imposition of a hard filter to exclude the FM waveband.
Fig. 2 Dynamic spectra of normalised intensity data taken by LOFAR station CS002 during the observation of Cas A on 18–19 August 2013. The dynamic spectrum of the entire observing period is given at the top, with zooms into three different hours of observation below to illustrate the range of conditions seen. White areas within the plots indicate where RFI was identified. |
RFI is still visible as white areas within the plots. These were identified by applying a median filter to the data using a window of (19.5 MHz × 4.2 s) to flatten out the scintillation pattern and then applying a threshold to identify the RFI. This method appears to be quite successful at identifying the RFI without also falsely identifying strong peaks in the scintillation as RFI. For subsequent analysis the RFI data points are replaced by an interpolation from nearby data, using the Python Astropy (Astropy Collaboration et al., 2013; Price-Whelan et al., 2018) library routine, “interpolate_replace_nans”. Normalisation of the data, to correct for long-period temporal variations in the system (e.g., gain variations resulting from the varying sensitivity of the receiving antenna array with source elevation), is carried out after RFI excision by dividing the intensities for each single frequency subband by a fitted 3rd-order polynomial.
When analysing the data, a variety of scattering conditions are observed during the course of the observation, as indicated in Figure 2. Different conditions also naturally occurred over the various international stations compared to those observed over the Dutch part of LOFAR. In this paper we therefore focus our analysis on only the first hour of observation and the measurements taken by the 24 core stations. This allows us to demonstrate the analysis techniques and to investigate the reason for the scintillation seen over this interval. Observations from later in this dataset undoubtedly show other effects and may be discussed in a subsequent publication.
3 LOFAR data analysis methods and results
3.1 Delay-Doppler spectra
The first stage of analysis was the calculation of delay-Doppler spectra: these were created from the dynamic spectra using 5-min time slices, advancing every minute through the observation, following the methods described in Fallows et al. (2014). To avoid regions more heavily contaminated by RFI, the frequency band used was restricted to 28.5–64.1 MHz. Example spectra from the first hour are presented in Figure 3.
Fig. 3 Example delay-Doppler spectra from the first hour of observation, taken using 5-min chunks of the dynamic spectrum from CS103 over the frequency band 28.5–64.1 MHz. |
The spectra show two clear arcs: the first is a steeper arc which varies in curvature throughout the first hour (henceforth labelled for convenience as the “primary arc”); the second is a very shallow arc (henceforth labelled as the “secondary arc”) which remains stable for the first 40 min of the observation before fading away. By the end of the first hour of observation the primary arc also becomes less distinctive for a short while before the delay-Doppler spectra again show distinctive structure, including a return of the secondary arc.
The variability of the curvature of the primary arc appears to follow a wave-like pattern during this part of the observation, as displayed in Figure 4. Here, simple parabolas involving only the square term (y = Cx 2 where C is the curvature) were plotted with various curvatures until a reasonable eyeball fit was achieved, and the resulting curvatures plotted for every minute of observation for the first hour. It proved impossible to achieve reasonable fits using least-squares methods due to confusion from non-arc structure in the spectra: fitting curvatures to these scintillation arcs is a well-known problem in the interstellar scintillation field and new methods of attempting this were presented at a recent workshop, but they are not easily described and have yet to be published. Hence, we do not attempt their application here.
Fig. 4 Curvatures of the steeper arc seen in delay-Doppler spectra calculated using data from CS103, from simple parabolas fitted by eye. The grey bounds represent an estimated error. |
The presence of two scintillation arcs likely indicates that scattering is dominated by two distinct layers in the ionosphere. A simple analysis, as described in Fallows et al. (2014), can be used to estimate the altitude of the scattering region with a basic formula relating arc curvature C to velocity V and distance L along the line of sight to the scattering region (Cordes et al., 2006):
The square term for the velocity illustrates the importance of gaining a good estimate of velocity to be able to accurately estimate the altitude of the scattering region via this method.
3.2 Scintillation pattern flow
The core area of LOFAR contains 24 stations within an area with a diameter of ~3.5 km. When viewing dynamic spectra from each of these stations it is clear that the scintillation pattern is mobile over the core (i.e., temporal shifts in the scintillation pattern are clear between stations) but does not necessarily evolve significantly. Therefore, the flow of the scintillation pattern over the core stations may be viewed directly by simply plotting the intensity received, for a single subband, by each station on a map of geographical station locations, for data from successive time steps. A movie (CasA_20130818_NL.mp4) of the scintillation pattern flow through the observation is published as an online supplement to this article (see Supplementary Material). The result, for 12 example time steps, is displayed in Figure 5, where a band of higher intensities can be seen to progress from north-west to south-east over the core. It should be noted that the data were integrated in time to 0.92 s for this purpose, to reduce both flicker due to noise and the duration of the movie. This does not average over any scintillation structure in this observation; structure with periodicities shorter than 1 s would be obvious in the delay-Doppler spectra as an extension of the arc(s) to greater than 0.5 Hz along the Doppler frequency axis.
Fig. 5 Normalised intensities received by all core stations at an observing frequency of 44.13 MHz, plotted on a geographical map of the stations. The intensities are colour-coded using a colour scale from yellow to purple with a range of 0.8–1.3 respectively. Times are at ~10 s intervals from 21:22:25 UT at top left to 21:24:15 UT at bottom right, and each plot uses data samples with an integration time of 0.92 s. Plot diameter is ~4.5 km. |
However, this is not the entire picture because the lines of sight from radio source to receivers are moving through the ionosphere as the Earth rotates, meaning that the scintillation pattern flow observed is a combination of flow due to the movement of density variations in the ionosphere and the movement of the lines of sight themselves through the ionosphere. Since the speed with which any single point on a line of sight passing through the ionosphere is dependent on the altitude of that point (the so-called ionosphere “pierce point”), this altitude needs to be either assumed or calculated to estimate a correction to the overall flow speed to obtain the natural ionospheric contribution. This introduces a natural uncertainty into estimates of velocity. Figure 6 shows the track of an ionospheric pierce-point at an assumed altitude of 200 km (an altitude chosen as representative of a typical F-region altitude where large-scale plasma structures are commonly observed) for the line of sight from core station CS002 to the radio source Cassiopeia A through the 7-h course of the observation to illustrate this movement. Although not the subject of this paper, it is worth noting that an east to west flow seen later in the observation appears to be solely due to the lines of sight moving across a mostly static ionospheric structure (see online movie of Supplementary Material), if the 200 km pierce point is assumed, further illustrating the necessity to take accurate care of the contribution from line of sight movement when assessing ionospheric speeds.
Fig. 6 Map showing the track of the 200 km pierce point of the line of sight from core station CS002 to Cassiopeia A from 2013-08-18T21:05:00 UT to 2013-08-19T04:05:00 UT. The thicker orange part of the track enhances the first hour of the observation. The black line winding a path across the centre of the image is the location of the border between the Netherlands and Germany. The location of CS002 is marked with a black star. |
The movie of the scintillation pattern flow, assuming a 200 km pierce point, shows a clear general north-west to south-east flow during the first hour of the observation, but also indicates some short (minutes) periods of confusion in which a north-east to south-west component might be just about discernable. Any second flow is likely to be associated with a second ionospheric layer and so warrants further investigation.
3.3 Estimating velocities
The representation of the scintillation pattern flow in movie form gives a direct and broad picture of the flow pattern and is very helpful in discovering short time-scale changes in speed and direction. However a cross-correlation analysis is still necessary to assess actual velocity(s). Correlation functions are calculated as follows:
-
Time series’ of intensity received by each station are calculated by averaging over the frequency band 55–65 MHz, with these frequencies chosen as the scintillation pattern remains highly correlated over this band.
-
For each 3-min data slice, advancing the start time of each successive slice by 10 s:
-
Calculate auto- and cross- power spectra using intensities from every station pair within the LOFAR core.
-
Apply low- and high-pass filters to exclude the DC-component and any slow system variation unlikely to be due to ionospheric effects, and white noise at the high spectral frequencies. The white noise is also subtracted using an average of spectral power over the high frequencies.
-
Inverse-FFT the power spectra back to the time domain to give auto- and cross-correlation functions.
-
In the analysis the high- and low-pass filter values were set to 0.01 Hz and 0.5 Hz respectively. This process results in a large set of cross-correlation functions for each time slice, each of which has an associated station–station baseline and a primary peak at, typically, a non-zero time delay from which a velocity can be calculated. However, the direction of the scintillation pattern flow still needs to be found for calculation of the actual velocity. For this, directions were assumed for each degree in the full 360° range of possible azimuth directions and the velocities re-calculated using the components of all baselines aligned with each assumed direction. This results, for each time slice, in 360 sets of velocities and from each set a median velocity and standard deviation about the median can be calculated (the median is used as this is less susceptible to rogue data points than the mean). The actual flow direction corresponds to the azimuth with the maximum median velocity and minimum standard deviation, as illustrated in Figure 7.
Fig. 7 Plots for single 3-min time slices of the median velocity and standard deviation of velocities about the median versus azimuth direction, calculated from the range of velocities found from all cross-correlation functions with the baselines within each station pair re-calculated for each assumed azimuth direction, in the usual form, counting clockwise from north. (a) Time slice commencing 21:05:00 UT using cross-correlations calculated after applying a high-pass filter at 0.01 Hz. (b) Time slice commencing 21:15:00 UT using cross-correlations calculated after applying a high-pass filter at 0.07 Hz. Note that the same y-axis is used for both velocity and standard deviation. |
From this analysis the primary velocity of ~20–40 ms−1 travelling from north-west to south-east is found, illustrated in Figure 7a, corresponding to the obvious scintillation pattern flow seen in the movie (see Supplementary Material). However, the presence of a second flow is still not obvious, although a hint of it can be seen in, for example, the second peak in the median velocity seen in Figure 7a.
A closer look at the auto-power spectra yielded the key to finding the second flow. Many spectra show a “bump” which can be viewed as being a second spectrum superposed on the main one. This is illustrated in Figure 8. To isolate this part of the spectrum, the spectra were re-filtered with a high-pass filter value of 0.07 Hz (the low-pass filter value remained the same), and correlation functions re-calculated. After following the same analysis as above to find median velocities and standard deviations, the second flow was found, as illustrated in Figure 7b.
Fig. 8 Example power spectrum calculated from 3 min of intensity data received by CS003. The black curve is the raw spectrum, the blue curve is the filtered and noise-subtracted spectrum. The locations of the low-pass filter and both high-pass filters used are illustrated. |
The analysis, using both high-pass filter values, has been carried out for the full data set. The velocities and associated directions in degrees azimuth for the first hour of the observation are given in Figure 9. Error bounds in the velocities are calculated as the standard deviation about the median of all velocity values available for the calculated azimuth direction.
Fig. 9 (a) Velocities calculated for the first hour of observation from cross-correlations created after filtering using the two different high-pass filter values. (b) Directions of these velocities, in degrees azimuth. |
The higher velocity (henceforth labelled as the “secondary velocity”) shows some scatter: periods where the secondary velocity drops to around the primary velocity values are due to the secondary velocity not being detected at these times; in these cases, it can still be detected in short-duration drops of velocity if correlation functions are re-calculated using an even higher high-pass filter value (the bump in these spectra appears shifted to slightly higher spectral frequencies). Values which decrease/increase towards/away from the primary velocity values likely represent a mix between the two velocities. The larger error bars seen in velocities may also be indicative in some instances of the standard deviation being broadened by some velocity values being more dominated by the other flow. The more extended period of scatter around 21:40–22:00 UT is a period where the secondary velocity is less apparent and the secondary scintillation arc fades from the delay-Doppler spectra. This indicates that the secondary structure is restricted in either space or time, either moving out of the field of view of the observation or ceasing for a period around 21:40 UT. It gives a first indication that the secondary velocity is associated with the secondary scintillation arc.
3.4 Estimating scattering altitudes
The velocities can now be used to estimate scattering altitudes, using the curvatures of the scintillation arcs and the simple formula given in equation (1). Initially the movement of the line of sight through the ionosphere is not accounted for, since this correction also requires an estimate of the pierce-point altitude to be reasonably calculated. Therefore an initial calculation of the scattering altitudes is made based on velocity values which are not corrected for this movement.
Using the primary velocities and combining these with the curvatures of the primary arc (Fig. 4) in equation (1), a range of distances, L, along the line of sight to the scattering region are found. These distances are converted to altitudes by accounting for source elevation (Cas A increased in elevation from 55°–64° during the first hour of observation). This process resulted in a range of altitudes to the scattering region of 200–900 km. Doing the same for the secondary velocities and applying an arc curvature of 3.2 ± 0.3 for the secondary scintillation arc gives estimated scattering altitudes of only ~70 km. If the primary/secondary velocities are combined vice-versa with the secondary/primary arc curvatures respectively, then the resulting scattering altitudes are clearly unreasonable (the secondary arc, primary velocity combination gives estimated altitudes of only ~10 km for example), lending further credence to the secondary velocity being associated with the secondary arc.
Velocity contributions from the line of sight movement are calculated as follows: for each time slice, t, the geographical locations beneath the pierce point of the line of sight through the ionosphere corresponding to the estimated scattering altitude at t are calculated, for both t and t + δt, where δt is taken as 3 min (the actual value is unimportant for this calculation). A velocity and its direction are found from the horizontal distance between these two locations and the direction of travel from one to the other. The general direction of the movement of the line of sight through the ionosphere is indicated by the orange line in Figure 6. Although the high scattering altitudes related to the primary scintillation arc and primary scintillation velocities lead to line-of-sight movements of up to ~35 ms−1, this movement is almost perpendicular to the direction of the primary scintillation velocity, limiting the actual contribution to ~5 ms−1. The line of sight movement is, however, in a very similar direction to the secondary velocities but the low corresponding scattering altitudes also limit the contribution in this case to ~5 ms−1.
An iterative procedure is then followed to correct the scintillation velocities for line-of-sight movement at the calculated scattering altitudes, re-calculate these altitudes, and re-calculate the line-of-sight movement. This procedure converges to a set of final scattering altitudes within five iterations. These are presented in Figure 10, with error bounds taken as the lowest and highest possible altitudes resulting from applying this procedure using the lower and upper limits of the arc curvature and scintillation velocity error bounds.
Fig. 10 Scattering altitudes estimated using equation (1), the primary velocities and primary scintillation arc curvatures (blue curve) and the secondary velocities and the curvature of the secondary scintillation arc (red dashed curves). |
The range of scattering altitudes encompassed by the error bounds is quite large in some instances, particularly where the calculated altitudes are higher. Although the square term for the velocity in equation (1) could lead to the natural conclusion that the error in the velocity dominates the error in scattering altitude, the errors in the velocity calculations are, for the most part, relatively small. Nevertheless, the error in the secondary velocity does appear to be the dominant error in the lower range of scattering altitudes (the red curves in Fig. 10). However, the dominant error for the higher range of scattering altitudes appears to be the scintillation arc curvatures, illustrating the importance of developing accurate fitting methods for these curvatures. Despite these concerns, it is clear that scattering is seen from two layers in the ionosphere; the primary scintillation arc arises from scattering in the F-region and the secondary scintillation arc arises from scattering much lower down in the D-region. Plasma decays by recombination with neutral species. In the F-region these densities are lower and so plasma lifetimes are longer than in the D-region. Typical plasma lifetimes in the F-region are of the order of hours, while they are of the order of minutes in the D-region. Hence the structures seen in each level may have a different source and time history.
4 Conditions in the ionosphere
We now investigate what the overall ionospheric conditions were at the time and hence the possible cause(s) of the scintillation seen by LOFAR at the time.
4.1 Geomagnetic conditions
The overall geomagnetic conditions at the time are given in Figure 11, which shows 24-h traces of the H-component of magnetic field for a representative set of magnetometers from the Norwegian magnetometer chain for 18 August 2013. Activity can be described as unsettled, with a minor substorm at high latitudes, peaking at the start of the LOFAR observation. However, geomagnetic activity remains quiet further south, and Kp took a value of 1 at 21 UT on 18th August 2013, indicating that this is unlikely to be a direct cause of the scintillation seen at LOFAR latitudes. We therefore investigate whether TIDs were present at the time and whether these could be consistent with the scintillation seen by LOFAR.
Fig. 11 Traces of the H-component of the geomagnetic field recorded on 18 August 2013 by a selection of magnetometer stations from the Norwegian chain. From top to bottom these are, along with their geomagnetic latitudes (2004, altitude 100 km): Longyearbyen (75.31° N), Bjørnøya (71.52° N), Nordkapp (67.87° N), Tromsø (66.69° N), Rørvik (62.28° N), and Karmøy (56.43° N). |
4.2 Ionosonde data
The presence of TIDs can be detected through the simultaneous appearance of wave-like structures on multiple sounding frequencies recorded by an ionosonde. This method is generally limited to a single point of observation and detection. The spatial extent of TIDs can be attempted by comparing multiple traces from different ionosondes, but this is limited by the low density of ionosondes in a given region. Measurements from the ionosonde in Chilton (UK) do indeed suggest the presence of wave-like patterns which, in principle, could be due to a large-scale TID propagating southward and/or MSTID triggered by a local Atmospheric Gravity Wave (Fig. 12).
Fig. 12 Multiple traces from the ionosonde in Chilton (UK) recorded between 20:00 18 August 2013 and 06:00 19 August 2013. |
4.3 GNSS data
However, measurements from ground-based GNSS receivers offer a more comprehensive view of the characteristics of any MSTIDs present (Kelley, 2009). In the present study, we focus on perturbations in the slant Total Electron Content (STEC) observed over the evening of 18 August 2013 from a network of GNSS stations around the LOFAR core stations (see Fig. 13). These stations are sufficient to infer the presence of TIDs and to infer the upper spatial scale-size limit of smaller-scale irregularities causing the intensity scintillation seen at LOFAR wavelengths.
Fig. 13 Map showing the locations of the GNSS stations used. |
The presence of TID-induced perturbations can be deduced from the presence of wave-like residuals on the STEC calculated for each satellite-receiver pair.
STEC was calculated and detrended following the methods of Hernández-Pajares et al. (2006), with the detrending carried out according to:
It is worth noting that the measured carrier phases L 1 and L 2 vary with time as a consequence of the motion of GNSS satellites relative to a given receiver on the Earth’s surface. As such, the spatial and temporal variabilities of ionisation gradients (such as those connected with TIDs and corresponding instabilities) become entangled. The various detrending methods (similar to Eq. (2)) lead to an estimate of ionisation gradients by considering temporal gradients only, with spatial and temporal variabilities intrinsically entangled in the GNSS observations.
Figure 14 shows examples of wave-like residuals on STEC for one pair of GNSS stations (Dentergem and Bruxelles in Belgium) observing the same GNSS satellite. The wave pattern is strongest over the first 2 h shown (18:00–20:00 UT) but then weakens considerably by the start of the LOFAR observation, although it remains evident. STEC from the observations of both stations appears well-correlated, with the Bruxelles dataset lagging behind that of Dentergem. Since Dentergem lies to the WNW of Bruxelles, this suggests a strong westerly component in the direction of travel, which could correspond with the secondary velocity seen by LOFAR.
Fig. 14 Example of a satellite-station pair. (a) PRN01 as observed on 18 August 2013 from Dentergem (DENT, blue line) and Bruxelles (BRUX, red line), both in Belgium, with baseline oriented from WNW to ESE. (b) Azimuth/elevation plot for PRN01 as observed from Dentergem. |
Figure 15 shows hourly plots of the overall geographical distribution of the STEC residuals calculated for all satellite passes seen within each hour by the GNSS stations used. The patterns shown in Figure 15 suggest a spatially and temporally varying propagation of MSTID wavefronts with components along the NE–SW as well as the NW–SE directions. Furthermore, the examples shown in Figure 15 also indicate the presence of smaller-scale ionisation structures in proximity to the wavefronts of the MSTIDs. This suggests that the scintillation seen by LOFAR is likely associated with the perpendicular propagation of two MSTIDs. However, the STEC variations here are also seen to fade by the start of the LOFAR observation.
Fig. 15 Hourly geographical distribution of all STEC perturbations in the evening of 18 August 2013: (a) 18:00–19:00 UT, (b) 19:00–20:00 UT, (c) 20:00–21:00 UT, and (d) 21:00–22:00 UT. |
A further illustration looks at the overall power spectral densities for the STEC residuals on all satellite-receiver pairs considered here over the hourly periods 20:00 UT–21:00 UT and 21:00 UT–22:00 UT (Fig. 16). The earlier hour is chosen alongside the hour covering the LOFAR data period as this better displays the components seen in the spectra. The temporal frequencies f can be converted into spatial scales L by assuming a given velocity V REL for the motion of the ionospheric structures across a GNSS raypath. That is:
(3)where V REL = V IONO−V SAT is the relative velocity between the velocity of the ionospheric structures and the scan velocity of a single raypath (at the same shell height). V SAT can be of the order of a few tens of ms−1 at 300 km.
Fig. 16 Power spectral densities of all the TEC residuals considered during the hours: (a) 20:00–21:00 UT and (b) 21:00–22:00 UT. The arrows indicate the two components considered in the text. |
There appear to be two main components in the energy cascade from larger to smaller ionisation scales: one with a period of 1666 s, and another component with a period of ~666 s. Taking V REL to be ~100 ms−1 (the secondary velocity seen by LOFAR as this is in a south-westerly direction and the example GNSS data in Figure 14 indicate a westerly component), these periodicities correspond to spatial scales of the order of 166 km and 66 km respectively. Beyond these scales the STEC analysis is limited by the sensitivity of the technique (Tsugawa et al., 2007), as the Power Spectral Densities reach the noise floor (Fig. 16). These orders of magnitudes suggest the presence of a larger-scale TID together with a smaller-scale TID (Kelley, 2009), while the energy cascade that can be observed through the Power Spectral Densities indicates that the large-scale structure breaks down into small-scale structures, likely owing to some instability mechanism.
4.4 Estimation of scale sizes of plasma structures
The scale sizes of the plasma structures causing the scintillation seen by LOFAR can also be calculated. The variations in the intensity of the received signal are caused by irregularities with a spatial scale size ranging from the Fresnel dimension to an order of magnitude below this value (Basu et al., 1998). The Fresnel length D F is related to the wavelength of the radio wave λ and the line of sight distance from the receiver to the scattering region L:
The Fresnel length was calculated for plasma structures at altitudes of 70 km, 200 km, 350 km, and 700 km, elevations of 55° and 64°, and at frequencies of 25.19 MHz, 35.15 MHz, and 60.15 MHz, and the results are shown in Table 1. The altitudes were chosen to cover the range of altitudes identified for the primary and secondary features in the LOFAR analysis, with the addition of 350 km as this altitude is commonly used within studies using GNSS satellites. The elevations of the radio source at the start and the end of the first hour of observation were used to establish the range of Fresnel scales for each altitude. The frequencies were chosen to match Figure 1.
The Fresnel length at altitudes of 70 km, 200 km, 350 km, and 700 km for three different frequencies received by LOFAR station CS002. The ranges represent calculation using the source elevation for the start and for the end of the first hour of observation. Values are in km.
Table 1 shows that the Fresnel length ranges between ~1 km and ~5 km and therefore the plasma structures causing the variations in signal intensity are likely to have a spatial scale size between ~100 m and ~5 km. The velocities calculated from the LOFAR data indicate that such structures would take tens of seconds to pass through the source-to-receiver line and the intensity variations in the observed signal occur on a similar timescale.
5 Further discussion
Geomagnetic activity was low in the mid-latitudes at the time, so enhanced activity was unlikely to be the direct cause of the scintillation observed. However, a weak sub-storm was seen at high latitudes and this reached its peak at the time of the start of the observation. An analysis of GNSS and ionosonde data reveals the presence of an MSTID travelling in the north-west to south-east direction. The larger-scale nature of this TID, and its direction of travel, are strongly consistent with the primary velocity and F-region scattering altitudes seen in the LOFAR observation. It is possible that this TID was caused by the geomagnetic activity at high latitude, but this is not confirmed. Simultaneously, an MSTID is also present travelling in a north-east to south-west direction which would most likely be associated with an atmospheric gravity wave propagating up from the neutral atmosphere. The smaller-scale nature of it, its direction of travel, and likely low-altitude source make it highly consistent with the secondary velocity and D-region scattering altitudes observed by LOFAR.
The amplitude of TID activity observed through GNSS STEC residuals decreased after 20:00 UT (as visible from Fig. 14 as well as from the comparison of hourly geographical maps in Fig. 15). However, the LOFAR observation did not start until 21:05 UT and the presence of scintillation on the radio frequencies observed by LOFAR remained significant for much of the first hour of observation. Whilst the presence of MSTIDs seems evident from the ionosonde multiple traces and GNSS STEC residuals in the region considered, their signatures do not appear simultaneously above the LOFAR core stations between 21:00 UT and 22:00 UT. This can be explained by the inability of GNSS to detect smaller amplitudes in STEC residuals, as the noise floor is encountered for observations with pierce points above the core LOFAR stations (Figs. 15 and 16). The scale sizes of plasma structures calculated for the LOFAR data indicate that these are an order of magnitude lower than those estimated from GNSS STEC. Smaller ionisation scales developing, for example, through the Perkins instability could induce scintillation on the VHF radio frequencies received by LOFAR but not on the L-band frequencies of GNSS. Hence, scintillation from these mid-latitude smaller-scale ionisation structures, formed through the Perkins instability in conjuction with the presence of TIDs, is likely to be what is detected through LOFAR.
6 Conclusions and outlook
This paper presents the results from one of the first observations of ionospheric scintillation taken using LOFAR, of the strong natural radio source Cassiopeia A taken overnight on 18–19 August 2013. The observation exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10–80 MHz. Delay-Doppler spectra from the first hour of observation showed two discrete parabolic arcs, one with a steep and variable curvature and the other with a shallow and static curvature, indicating that the scintillation was the result of scattering through two distinct layers in the ionosphere.
A cross-correlation analysis of the data received by stations in the LOFAR core reveals two different velocities in the scintillation pattern: a primary velocity of ~20–40 ms−1 is observed travelling in a north-west to south-east direction, which is associated with the primary parabolic arc and altitudes of the scattering layer varying in the range ~200–700 km. A secondary velocity of ~110 ms−1 is observed travelling in a north-east to south-west direction, which is associated with the secondary arc and a much lower scattering altitude of ~60–70 km. The latter velocity is associated with a secondary “bump” seen at higher spectral frequencies in power spectra calculated from time series’ of intensities, indicating that it is more strongly associated with smaller-scale structure in the ionosphere.
GNSS and ionosonde data from the time suggest the presence of two MSTIDs travelling in perpendicular directions. The F-region scattering altitudes calculated from the LOFAR primary scintillation arc and primary velocity, and the larger density scales associated with this, suggest that this is associated with a larger-scale TID seen in GNSS data potentially resulting from high-latitude geomagnetic activity. The D-region scattering altitudes of the secondary arc and secondary velocity suggest an atmospheric gravity wave source for a smaller-scale TID. These TIDs trigger an instability which leads to the breakdown of the large-scale density structure into smaller scales, giving rise to the scintillation observed. In the mid-latitude ionosphere the Perkins mechanism is the most likely instability and the features of the smaller-scale density variations observed seem consistent with this. To the best of our knowledge this is the first time that two TIDs have been directly observed simultaneously at different altitudes.
This observation demonstrates that LOFAR can be a highly valuable tool for observing ionospheric scintillation in the mid-latitudes over Europe and enables methods of analysis to be used which give greater insight into the likely sources of scattering and could be used to improve modelling of them. With a far greater range of frequencies (multi-octave if the LOFAR high-band is also used) and fine sampling both across the frequency band and in time, LOFAR observations offer a wider sensitivity than that available to GNSS measurements. The analysis techniques shown in this paper also demonstrate that LOFAR can observe ionospheric structures at different altitudes simultaneously; a capability not commonly available for GNSS observations. It also complements these measurements by probing potentially different scintillation regimes to those observed by GNSS.
Since this observation was taken, many more have been carried out under a number of projects, recording ionospheric scintillation data at times when the telescope would otherwise be idle. These demonstrate a wide range of scintillation conditions over LOFAR, some of which are seen only very occasionally and perhaps by only one or two of the international stations, illustrating the value to be had by monitoring the ionosphere at these frequencies. A design study, LOFAR4SpaceWeather (LOFAR4SW – funded from the European Community’s Horizon 2020 Programme H2020 INFRADEV-2017-1 under grant agreement 777442) currently underway will design a possible upgrade to LOFAR to enable, amongst other space weather observations, ionospheric monitoring in parallel with the regular radio astronomy observations. Such a design, if implemented, would enable a full statistical study of ionospheric scintillation at these frequencies, alongside the advances in scintillation modelling and our understanding of the ionospheric conditions causing it which can be gleaned in focussed studies such as that presented here.
Supplementary materials
CasA_20130818_NL.mp4: A movie depicting the scintillation pattern flow through the observation. Access here
Acknowledgments
This paper is based on data obtained with the International LOFAR Telescope (ILT) under project code “IPS”. LOFAR (van Haarlem et al., 2013) is the Low Frequency Array designed and constructed by ASTRON. It has observing, data processing, and data storage facilities in several countries, that are owned by various parties (each with their own funding sources), and that are collectively operated by the ILT foundation under a joint scientific policy. The ILT resources have benefitted from the following recent major funding sources: CNRS-INSU, Observatoire de Paris and Université d’Orléans, France; BMBF, MIWF-NRW, MPG, Germany; Science Foundation Ireland (SFI), Department of Business, Enterprise and Innovation (DBEI), Ireland; NWO, The Netherlands; The Science and Technology Facilities Council, UK; Ministry of Science and Higher Education, Poland. The work carried out at the University of Bath was supported by the Natural Environment Research Council (grant number NE/R009082/1) and by the European Space Agency/Thales Alenia Space Italy (H2020-MOM-TASI-016-00002). We thank Tromsø Geophysical Observatory, UiT the Arctic University of Norway, for providing the lyr, bjn, nor, tro, rvk, and kar magnetometer data. The Kp index and the Chilton ionosonde data were obtained from the U.K. Solar System Data Centre at the Rutherford Appleton Laboratory. Part of the research leading to these results has received funding from the European Community’s Horizon 2020 Programme H2020-INFRADEV-2017-1under grant agreement 777442. The editor thanks two anonymous reviewers for their assistance in evaluating this paper.
References
- Aarons J. 1982. Global morphology of ionospheric scintillations. Proc IEEE 70(4): 360–378. https://doi.org/10.1109/PROC.1982.12314. [CrossRef] [Google Scholar]
- Astropy Collaboration, Robitaille TP, Tollerud EJ, Greenfield P, Droettboom M, et al. 2013. Astropy: A community Python package for astronomy. A&A 558: A33. https://doi.org/10.1051/0004-6361/201322068. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Basu S, Weber E, Bullett T, Keskinen M, MacKenzie E, Doherty P, Sheehan R, Kuenzler H, Ning P, Bongiolatti J. 1998. Characteristics of plasma structuring in the cusp/cleft region at Svalbard. Radio Sci 33(6): 1885–1899. https://doi.org/10.1029/98RS01597. [CrossRef] [Google Scholar]
- Cordes JM, Rickett BJ, Stinebring DR, Coles WA. 2006. Theory of parabolic arcs in interstellar scintillation spectra. Astrophys J 637(1): 346. https://doi.org/10.1086/498332. [NASA ADS] [CrossRef] [Google Scholar]
- de Gasperin F, Mevius M, Rafferty D, Intema H, Fallows R. 2018. The effect of the ionosphere on ultra-low-frequency radio-interferometric observations. A&A 615: A179. https://doi.org/10.1051/0004-6361/201833012. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Emardson R, Jarlemark P, Johansson J, Scäfer S. 2013. Spatial variability in the ionosphere measured with GNSS networks. Radio Sci 48: 646–652. https://doi.org/10.1002/2013RS005152. [CrossRef] [Google Scholar]
- Fallows R, Coles W, McKay-Bukowski D, Vierinen J, Virtanen I, et al. 2014. Broadband meter-wavelength observations of ionospheric scintillation. J Geophys Res: Space Phys 119(12): 10,544–10,560. https://doi.org/10.1002/2014JA020406. [CrossRef] [Google Scholar]
- Hapgood M. 2017. Satellite navigation – amazing technology but insidious risk: Why everyone needs to understand space weather. Space Weather 15(4): 545–548. https://doi.org/10.1002/2017SW001638. [CrossRef] [Google Scholar]
- Hernández-Pajares M, Juan JM, Sanz J. 2006. Medium-scale travelling ionospheric disturbances affecting GPS measurements: Spatial and temporal analysis. J Geophys Res 111: A07S11. https://doi.org/10.1029/2005JA011474. [Google Scholar]
- Hernández-Pajares M, Juan JM, Sanz J, Aragón-Àngel A. 2012. Propagation of medium scale travelling ionospheric disturbances at different latitudes and solar cycle conditions. Radio Sci 47: RS0K05. https://doi.org/10.1029/2011RS004951. [CrossRef] [Google Scholar]
- Kelley MC. 2009. The Earth’s ionosphere: Plasma physics and electrodynamics, In: Vol. 96 of International Geophysics Series, 2nd edn., Elsevier, San Diego, CA. [Google Scholar]
- Kelley MC. 2011. On the origin of mesoscale TIDs at midlatitudes. Ann Geophys 29: 361–366. https://doi.org/10.5194/angeo-29-361-2011. [CrossRef] [Google Scholar]
- Knepp DL, Nickisch L. 2009. Multiple phase screen calculation of wide bandwidth propagation. Radio Sci 44(1): 1–11. https://doi.org/10.1029/2008RS004054. [CrossRef] [Google Scholar]
- McKay-Bukowski D, Vierinen J-P, Virtanen I, Fallows R, Postila M, et al. 2014. KAIRA: The Kilpisjärvi atmospheric imaging receiver array – system overview and first results. IEEE Trans Geosci Remote Sens 53(3): 1440–1451. https://doi.org/10.1109/TGRS.2014.2342252. [CrossRef] [Google Scholar]
- Price-Whelan AM, Sipöcz BM, Günther HM, Lim PL, Crawford SM, et al. 2018. The Astropy project: Building an open-science project and status of the v2.0 core package. Astron J 156: 123. https://doi.org/10.3847/1538-3881/aabc4f. [CrossRef] [Google Scholar]
- Saito A, Fukao S. 1998. High resolution mapping of TEC perturbations with the GSI GPS network over Japan. Geophys Res Lett 25(16): 3079–3082. [CrossRef] [Google Scholar]
- Stinebring D, McLaughlin M, Cordes J, Becker K, Goodman JE, Kramer M, Sheckard J, Smith C. 2001. Faint scattering around pulsars: probing the interstellar medium on solar system size scales. Astrophys J Lett 549(1): L97. https://doi.org/10.1086/319133. [NASA ADS] [CrossRef] [Google Scholar]
- Tsugawa T, Saito A, Otsuka Y. 2004. A statistical study of large-scale traveling ionospheric disturbances using the GPS network in Japan. J Geophys Res: Space Phys 109: A06302. [CrossRef] [Google Scholar]
- Tsugawa T, Otsuka Y, Coster AJ, Saito A. 2007. Medium-scale travelling ionospheric disturbances detected with dense and wide TEC maps over North America. Geophys Res Lett 34(L22): 101. https://doi.org/10.1029/2007GL031663. [CrossRef] [Google Scholar]
- van Haarlem MP, Wise MW, Gunst AW, Heald G, McKean JP, et al. 2013. LOFAR: The LOw-Frequency ARray. A&A 556: A2. https://doi.org/10.1051/0004-6361/201220873. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Cite this article as: Fallows RA, Forte B, Astin I, Allbrook T, Arnold A, et al. 2020. A LOFAR observation of ionospheric scintillation from two simultaneous travelling ionospheric disturbances. J. Space Weather Space Clim. 10, 10.
All Tables
The Fresnel length at altitudes of 70 km, 200 km, 350 km, and 700 km for three different frequencies received by LOFAR station CS002. The ranges represent calculation using the source elevation for the start and for the end of the first hour of observation. Values are in km.
All Figures
Fig. 1 (a) Time series of intensity received at three discrete frequencies by LOFAR station CS002 during the observation of Cas A on 18–19 August 2013, plus, (b) and (c) power spectra of two 10-min periods within these time series’. |
|
In the text |
Fig. 2 Dynamic spectra of normalised intensity data taken by LOFAR station CS002 during the observation of Cas A on 18–19 August 2013. The dynamic spectrum of the entire observing period is given at the top, with zooms into three different hours of observation below to illustrate the range of conditions seen. White areas within the plots indicate where RFI was identified. |
|
In the text |
Fig. 3 Example delay-Doppler spectra from the first hour of observation, taken using 5-min chunks of the dynamic spectrum from CS103 over the frequency band 28.5–64.1 MHz. |
|
In the text |
Fig. 4 Curvatures of the steeper arc seen in delay-Doppler spectra calculated using data from CS103, from simple parabolas fitted by eye. The grey bounds represent an estimated error. |
|
In the text |
Fig. 5 Normalised intensities received by all core stations at an observing frequency of 44.13 MHz, plotted on a geographical map of the stations. The intensities are colour-coded using a colour scale from yellow to purple with a range of 0.8–1.3 respectively. Times are at ~10 s intervals from 21:22:25 UT at top left to 21:24:15 UT at bottom right, and each plot uses data samples with an integration time of 0.92 s. Plot diameter is ~4.5 km. |
|
In the text |
Fig. 6 Map showing the track of the 200 km pierce point of the line of sight from core station CS002 to Cassiopeia A from 2013-08-18T21:05:00 UT to 2013-08-19T04:05:00 UT. The thicker orange part of the track enhances the first hour of the observation. The black line winding a path across the centre of the image is the location of the border between the Netherlands and Germany. The location of CS002 is marked with a black star. |
|
In the text |
Fig. 7 Plots for single 3-min time slices of the median velocity and standard deviation of velocities about the median versus azimuth direction, calculated from the range of velocities found from all cross-correlation functions with the baselines within each station pair re-calculated for each assumed azimuth direction, in the usual form, counting clockwise from north. (a) Time slice commencing 21:05:00 UT using cross-correlations calculated after applying a high-pass filter at 0.01 Hz. (b) Time slice commencing 21:15:00 UT using cross-correlations calculated after applying a high-pass filter at 0.07 Hz. Note that the same y-axis is used for both velocity and standard deviation. |
|
In the text |
Fig. 8 Example power spectrum calculated from 3 min of intensity data received by CS003. The black curve is the raw spectrum, the blue curve is the filtered and noise-subtracted spectrum. The locations of the low-pass filter and both high-pass filters used are illustrated. |
|
In the text |
Fig. 9 (a) Velocities calculated for the first hour of observation from cross-correlations created after filtering using the two different high-pass filter values. (b) Directions of these velocities, in degrees azimuth. |
|
In the text |
Fig. 10 Scattering altitudes estimated using equation (1), the primary velocities and primary scintillation arc curvatures (blue curve) and the secondary velocities and the curvature of the secondary scintillation arc (red dashed curves). |
|
In the text |
Fig. 11 Traces of the H-component of the geomagnetic field recorded on 18 August 2013 by a selection of magnetometer stations from the Norwegian chain. From top to bottom these are, along with their geomagnetic latitudes (2004, altitude 100 km): Longyearbyen (75.31° N), Bjørnøya (71.52° N), Nordkapp (67.87° N), Tromsø (66.69° N), Rørvik (62.28° N), and Karmøy (56.43° N). |
|
In the text |
Fig. 12 Multiple traces from the ionosonde in Chilton (UK) recorded between 20:00 18 August 2013 and 06:00 19 August 2013. |
|
In the text |
Fig. 13 Map showing the locations of the GNSS stations used. |
|
In the text |
Fig. 14 Example of a satellite-station pair. (a) PRN01 as observed on 18 August 2013 from Dentergem (DENT, blue line) and Bruxelles (BRUX, red line), both in Belgium, with baseline oriented from WNW to ESE. (b) Azimuth/elevation plot for PRN01 as observed from Dentergem. |
|
In the text |
Fig. 15 Hourly geographical distribution of all STEC perturbations in the evening of 18 August 2013: (a) 18:00–19:00 UT, (b) 19:00–20:00 UT, (c) 20:00–21:00 UT, and (d) 21:00–22:00 UT. |
|
In the text |
Fig. 16 Power spectral densities of all the TEC residuals considered during the hours: (a) 20:00–21:00 UT and (b) 21:00–22:00 UT. The arrows indicate the two components considered in the text. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.