Auroral kilometric radiation (AKR) is a terrestrial radio emission excited by the same accelerated electrons which excite auroral emissions. Although it is well correlated with auroral and geomagnetic activity, the coupling timescales between AKR and different magnetospheric or ionospheric regions have yet to be determined. Estimation of these coupling timescales is non-trivial as a result of complex, non-linear processes which rarely occur in isolation. In this study, the mutual information between AKR intensity and different geomagnetic indices is used to assess the correlation between variables. Indices are shifted to different temporal lags relative to AKR intensity, and the lag at which the variables have the most shared information is found. This lag is interpreted as the coupling timescale. The AKR source region receives the effects of a shared driver before the auroral ionosphere. Conversely, the polar ionosphere reacts to a shared driver before the AKR source region. Bow shock interplanetary magnetic field BZ is excited about 1 h before AKR enhancements. This work provides quantitatively determined temporal context to the coupling timelines at Earth. The results suggest that there is a sequence of excitation following the onset of a shared driver: first, the polar ionosphere feels the effects, followed by the AKR source region and then the auroral ionosphere.
Science Foundation Ireland18/FRL/6199Irish Research CouncilGOIPD/2022/782Air Force Office of Scientific ResearchFA8655-22-7056Engineering and Physical Sciences Research CouncilEP/L015382/1Introduction
Auroral kilometric radiation (AKR) is the strongest natural terrestrial radio emission and is excited by the same electron precipitation which stimulates auroral emission in the ionosphere. AKR is generated by electron–cyclotron maser instability (ECMI), a seminal theory developed by . While describing the generation mechanism for AKR emission, this theory also stipulates that AKR emission frequency is inversely proportional to altitude and is beamed from the source region at near-right angles.
AKR was first observed in the 1960s by and has been studied in detail since by instruments on board spacecraft including IMP 6 and 8, Hawkeye, Wind, GEOTAIL, POLAR, IMAGE, the Cluster array and Cassini (e.g. ). With its main band appearing between 100 and 400 kHz, AKR is more broadly observed between 30 and 800 kHz , with powers up to 109 W (e.g. ). AKR observations have been shown to be well correlated with substorm activity (e.g. ) and indeed the Auroral Electrojet (AE) Index . Furthermore, showed enhancements in AKR power and frequency expansion at substorm onset over a decade of Wind/WAVES data. Finally, AKR is an excellent indicator of geomagnetic disturbance, particularly at nightside local times (LTs) .
AKR is generated in a region of low plasma density known as the auroral plasma cavity , within which precipitating electrons are observed . Energetic electrons within these auroral acceleration regions are propelled down magnetic field lines towards the ionosphere, often following a dipolarisation of the tail magnetic field. Some electrons precipitate into the ionosphere (exciting the aurora). However, depending on the pitch angle, other electrons may, by conservation of magnetic moments, undergo reflection at the magnetic mirror points.
These electrons travel upwards towards the auroral plasma cavity (e.g. ), where there is not enough plasma to contain the energy of the incoming electrons (e.g. ) (plasma density is sufficiently low that collisions are rare). Along with particles of other pitch angles, the electrons undergo wave–particle interactions and emit their energy via the ECMI as AKR. Simply put, the electron distribution inside the auroral plasma cavity has a variety of velocity–space gradients that can lead to the generation of waves.
AKR is anisotropically beamed in a hollow cone at angles nearly perpendicular to the source region , although some undergo strong refraction along the inner edges of the auroral cavity, which adds complexity to the beaming . This generates a statistical illumination region based on radio sources in both hemispheres, with generally right-handed circularly polarised emission from the northern magnetic (southern geographic) hemisphere and left-handed circularly polarised emission from the southern magnetic (northern geographic) hemisphere. The manner in which these two sources interact generates a statistical shadow zone at equatorial magnetic latitudes close to the planet (inside ≈5 Earth radii on the nightside, e.g. ), where AKR is not observed, and an observer at higher latitudes or further distance may see one or other source or a superposition of both. This anisotropic beaming creates challenges for observing AKR, since a spacecraft orbiting Earth will transit into and out of the illumination regions and observe changes in AKR emission as a result of solar wind–magnetosphere coupling.
Although AKR has been observed at all LTs , it is most reliably viewed on the nightside between 18:00 and 06:00 LT (e.g. ), relating to the statistical position of the source region. The observed power decreases as 1R2 (where R is the distance from the source region, ), resulting in higher powers being observed closer to the source region. Finally, a 24 h modulation of the AKR signal has been observed as a result of the diurnal precession of the tilted dipole magnetic field .
Previous AKR observations have highlighted the link between the emissions and geomagnetic activity. The intensity and frequency range of AKR have been shown to correlate with the AE Index , which is often explained by its strong link to substorm activity . Similarly, AKR enhancements are sometimes observed simultaneously with auroral brightenings (e.g. ). More broadly, longitudinal extensions of the source region are observed with increasing geomagnetic activity , enabling AKR observations at dayside LTs. Additionally, showed higher geomagnetic activity in AU, AL, PC(N) and SYM-H during AKR bursts and even greater enhancements in activity during the most intense AKR events. Finally, showed that prolonged southward interplanetary magnetic field (IMF) BZ enhanced AKR emission during the passing of a magnetic cloud event. Despite limitations with single spacecraft viewing, this study highlighted links between IMF BZ driving and the excitation of AKR emissions.
The terrestrial ionosphere is permeated by interwoven current systems, which are a barometer for energy transfer throughout the magnetosphere. As with any electrical circuit, when one current is excited, those that connect with it respond accordingly, ensuring current continuity. In this study, the horizontal currents in the polar and high-latitude ionosphere are characterised using the polar and auroral indices. For a detailed description of the role of ionospheric currents in the magnetosphere, the reader is directed to and (and references therein).
In this study, the geomagnetic indices AE, AU, AL, PC(N) and SYM-H (obtained via OMNIWeb; ) are utilised to examine coupling timescales within the magnetosphere. These geomagnetic indices are chosen for the study as they characterise auroral, polar and equatorial activity respectively and hence map to different regions of the magnetosphere. Comparing and contrasting their relationships with AKR will give a global view of AKR's role in the magnetosphere.
The PC index monitors activity in the polar ionosphere, and in the northern geographic (containing a southern magnetic pole) hemisphere it is derived from a magnetometer station at around 85° latitude in Greenland. Geomagnetic observatories such as these measure deflections in the magnetic field as a result of changes in overhead current systems, and as such the PC index can be used to probe the state of polar ionospheric currents. At such latitudes, a horizontal Hall current (pointing sunwards) characterises the speed of magnetic flux transport across the polar cap (e.g. ). Although the PC index is not an exact measure of the strength of these currents, an enhancement in the PC index demonstrates an excitation of these polar currents and hence an intensification of the tailward magnetic flux transport. For example, following the onset or enhancement of dayside reconnection , an increase in PC may be observed as flux builds up in the tail prior to substorm onset. Additionally, the PC index responds to changes in solar wind dynamic pressure (e.g. ) and magnitude of the interplanetary magnetic field and as such is often regarded as a monitor of solar wind energy input into the magnetosphere.
At auroral latitudes, the commonly used AE Index and its upper and lower envelopes (AU and AL respectively) characterise the strength of activity in the auroral zone. Similarly to the PC index, they are derived from magnetometer stations at auroral latitudes, which observe deflections in the geomagnetic field as a result of overhead currents, in this case the auroral electrojets. Notably an indicator of substorm activity, AL indicates the strength of westward electrojets, including the substorm electrojet (the ionospheric portion of the substorm current wedge). Similarly, AU indicates activity in eastward electrojets.
Finally, SYM-H is derived from magnetometer stations at equatorial latitudes and as such measures the strength of the ring current. The ring current is a westward-flowing equatorial magnetospheric current system formed from a combination of gradient and curvature drift of plasma in the dipole magnetic field. This system is excited by a large deposition of energy in the main phase of a geomagnetic storm. This results in well-observed characteristic signatures in the SYM-H index (e.g. ), which is often termed the ring current index. SYM-H also exhibits step-like increases as a result of rapid magnetospheric compressions known as sudden commencements (e.g. ).
In this paper, mutual information (MI) is used to characterise the strength of the relationship between two time series, without any knowledge of the order of their relationship. MI quantifies the correlation, linear or non-linear, between two sequences in a non-parametric (and hence model-independent) manner. Although this is an emerging technique in the field of solar system plasma physics, previous authors have demonstrated the use of MI as a powerful tool in understanding non-linear dynamics, and the technique is well established in the wider community . characterised the MI content between the product of IMF BZ with solar wind speed and the auroral electrojet index. They also investigated the effect of various different propagation techniques for the IMF and solar wind data. investigated the spatial correlation properties of the solar wind and used the normalised MI to confirm results without the assumption of linearity. used information theory to show that the northward IMF turnings observed around substorm onset are likely coincidental. Additionally, showed a relationship between solar flares and subsequent solar flares by calculating the MI of time-lagged waiting time distributions of flare events. Finally, used MI to characterise the coupling timescale between evidence of plasma injections in the Kronian magnetosphere and narrowband emissions.
In this study, the coupling timescales between the AKR intensity measured by Wind/WAVES and the geomagnetic indices described above (AE, AU, AL, PC(N), SYM-H) as well as IMF BZ are estimated. The shared information content between an AKR intensity time series paired with different geomagnetic parameters will be estimated using MI. The geomagnetic indices and IMF are time-lagged with respect to AKR intensity, with lags -60≤τ≤+60 min for indices and -120≤τ≤+30 min for IMF BZ. The lags for geomagnetic indices and IMF BZ are different due to expected differences in their relationship with AKR. Both the indices and AKR are measures of geomagnetic activity and hence will be more closely temporally aligned. IMF BZ however is a driver of the magnetosphere, so there may be a greater time lag between BZ changes and AKR excitation. Finally, in the initial stages of the work, the peak in the MI vs. lag data was found to be at lower (more negative) lags for IMF BZ. The lag at which the MI peaks gives an estimate of the time taken for information from a shared driver to propagate between different regions in the magnetosphere, i.e. the coupling timescale. The method and results of this analysis are described in Sect. , followed by concluding remarks.
Estimation of coupling timescalesData and method
AKR emission is selected from amongst a complex superposition of phenomena in Wind/WAVES/RAD1 data using the W21 technique (described in detail by ); a subset of these data is available online . Using these W21-selected data, the AKR intensity is integrated over 100–650 kHz (the same frequency band as , and ), utilising the technique described in detail by and . AKR intensity is calculated over a full frequency sweep in Wind/WAVES data, roughly 183 s (with some fraction of seconds that varies from integration to integration), normalising to 1 AU. For a detailed description of intensity calculation, the reader is directed to and . Since the geomagnetic indices which AKR will be compared with are provided at integer minute resolution, one or another of the datasets must be interpolated to allow point-to-point comparison. As the resolution of AKR power is uneven and not at integer minutes and interpolating the indices would be downsampling the data, the AKR power is interpolated to match the resolution of the geomagnetic indices. This results in AKR intensity measurements and the geomagnetic indices PC(N), AE, AU, AL, SYM-H, and IMF BZ all being on the same timescale for comparison.
AKR intensity is calculated over 10 years of data. AKR observations from 1995 to 2004 inclusive from Wind/WAVES/RAD1 are used and compared with available geomagnetic indices AE/AU/AL, PC(N), SYM-H and IMF BZ which are described in detail in the Introduction section. The years 1995 to 2004 are selected as this ranges from just after the launch of Wind (November 1994) to the year when Wind leaves the near-Earth environment and begins its journey to the L1 point. In the interim 10 years, Wind samples a broad parameter space of locations in the near-Earth environment, providing a range of AKR observations (it is important to note that AKR is strongly anisotropically beamed).
Where the relationship between two parameters is linear, the Pearson correlation coefficient is often used as a measure of the strength and direction of the relationship between two variables (e.g. ). However, initial investigations as well as previous studies (e.g. ) indicate a more complex relationship between excitations in auroral indices and AKR power. The MI (measured in bits) between two time series is a measure of the relationship between two variables and is sensitive to both linear and non-linear relationships. Simply put, the MI content between two time series describes the information learned about one time series by observing the other (e.g. ) and vice versa.
The standard definition of the MI content, Ia,b, between two signals a and b in terms of entropy, H, is Ia,b=Ha+Hb-Ha,b,
where entropies are defined in terms of the probability distributions of time series a and b. In this study, the mutual information is calculated using the function mutual_info_regression (standard options) within the scikit-learn Python package , which is based on the method by (following on from work by ). The estimate for MI is based on k-nearest neighbour statistics and is given by
Ia,b=ψk-〈ψna+1+ψnb+1〉+ψ(N),
where ψ is a digamma function, k is the number of nearest neighbours, na is the number of measurements in a, nb is the number of measurements in b, and N is the number of measurements overall. It is important to note at this point that the variables a and b used in this paper are discrete samples of continuous variables derived from empirical observations from both ground-based magnetometers and a space-based radio instrument. Rather than varying from -1 to 1 as for the Pearson correlation coefficient, I is positive, with higher values indicating more shared information content between two time series.
The method is used here as it has some advantages over the traditional method. In particular, it avoids bias which comes from the binning of the data into probability distribution functions. Additionally, the method is “data-efficient”, resolving structure to the smallest scales, and “adaptive”, meaning that the resolution improves as the amount of data increases. The reader is directed to for further details.
Random phase surrogates are a widely used method to test the results against the null hypothesis that there is no time structure in the time series. In order to assess the significance of calculated MI values, each calculation is repeated using a random phase surrogate of AKR intensity. The surrogate AKR intensity is generated using an amplitude-adjusted Fourier transform , which is calculated using the aaft Python package. The resulting surrogate AKR retains the amplitude information but not the temporal information in the AKR time series: both the surrogate and observed AKR intensities have very similar probability distribution functions (e.g. ). Since the surrogate AKR has no temporal relationship with magnetospheric activity, MI values between the surrogate and geomagnetic indices should be lower than for the real AKR intensity: the surrogate is used to quantify the MI content in correlations that occur “by chance”.
While using the MI estimate devised by (Eq. ), this investigation exploits a similar overall methodology to that employed by , who gave a detailed description of MI calculation and estimated the coupling timescales between the AE index and vxBZ (the product of the solar wind velocity x component and IMF BZ). Also comparing different propagation methods for vx and BZ measurements from the Wind spacecraft to Earth, they calculated the MI between the two time series for different values of “additional time lag”, which characterised the response time of the magnetospheric system.
In this study, a similar time lag approach is used to assess the strength of the relationship between time series of AKR and geomagnetic indices.
ResultsMI coupling analysis
Similarly to , the MI content between time series of AKR power and geomagnetic indices with different τ is calculated and presented as black crosses in Fig. . Note that the AKR power timestamps are fixed, and geomagnetic indices' time series are shifted by -60≤τ≤+60 min, while IMF BZ is shifted by -120≤τ≤+30 min. For each parameter, a piecewise linear (purple) and quadratic (gold) fit to the MI content between time series (black crosses) is performed, and the peak of each of those fits is indicated with a dashed vertical line in the corresponding colour. The piecewise linear curve is fitted to the data as one curve with free parameters including the position of the turnover. For all the parameters except PC(N) and SYM-H, the piecewise linear fit to the data has a lower root-mean-squared error (RMSE) than the quadratic fit, although the RMSE values are of the same order of magnitude for both fits in AU, PC(N) and SYM-H.
To determine whether the calculated MI values are significant, they are compared with the MI content between the surrogate AKR and the indices (see the description in Sect. ). The mean of the MI values between the surrogate AKR and each index across all the lag times is indicated with a blue arrow in each panel of Fig. . For all the indices, this “threshold” value falls at least 1 order of magnitude below the MI represented by the black crosses, confirming that the MI between the AKR and indices is statistically significant.
It is important to note that this analysis was also run on IMF total magnitude and BY as well as solar wind flow pressure, density and speed. No clear trend in MI as a function of lag was found, and indeed sometimes the MI was below the random phase surrogate threshold. Therefore, only results for parameters with clear MI vs. lag trends (i.e. AE, AU, AL, PC(N), SYM-H and BZ) are presented for brevity. Additionally, any time intervals with missing data (which are only found in PC(N) and BZ) were removed from both the parameter dataset and the AKR dataset before running the MI-lagging analysis.
Before discussing the MI results, it is important to note that this time-lagging technique was also applied to the data using the Pearson correlation coefficient r rather than MI, finding a line of best fit and its related r value. The Pearson correlation coefficient is a measure of the strength and direction of a linear relationship between two variables. r varies from -1 to 1, where the sign indicates whether the relationship is positive or negative. Values of r close to 0 indicate that y may increase or decrease as x increases . Conversely, values of r close to +1 (-1) suggest that y will increase (decrease) as x increases. For all the geomagnetic indices, the value of r across all the lag times was below 0.15. This suggests that the linear correlations between the data are low – since the values are much closer to 0 than 1. This emphasises the non-linear nature of the relationship between AKR intensity and geomagnetic indices (e.g. ) and the need for a non-linear measure of correlation, such as MI.
MI content between the labelled parameter and AKR-integrated intensity as a function of time lag between the two variables (black crosses). Lag time is the time shift applied to the geomagnetic index data. The parameters shown are (a) AE, (b) AU, (c) AL, (d) PC(N), (e) SYM-H and (f) IMF BZ. A piecewise linear (purple) and quadratic (gold) fit to the MI data is presented, with the position of the MI peak indicated with a dashed line in the corresponding colour and labelled with the x position. The corresponding shade represents the 80 % confidence interval. The mean of the MI content between each index and a random phase surrogate of AKR intensity is indicated with a blue arrow not in scale with the y axis where necessary.
For the auroral indices, the peak MI content between the indices and AKR power is within ±10 min. Taking only the peak of the better-fitting curve, each of the auroral indices has the most shared information with the AKR power with a small positive temporal lag. Peaks are found at +7 min for AE (Fig. a), +1 min for AU (Fig. b) and +8 min for AL (Fig. c). Errors in these values are found at the intersection of a horizontal line threading the peak and edges of the 80 % confidence interval, represented as a shade, and are presented in Table . Given previous work suggesting that AKR excitations are driven by the same phenomena which excite the auroral zone (which is characterised by AE/U/L), this suggests that the AKR source region (between 1.8 RE and 3 RE in altitude) feels these effects before the auroral region. It is interesting to note that the peak MI is much better defined for AE and AL, perhaps as a result of the strong relationship between AKR and substorm activity characterised by AL (unlike the study by , which showed the best correlation between AKR and AU). Additionally, the error on the calculated AU peak is larger (as detailed in Table ), and excitations in AKR and AU are more closely temporally aligned.
Fit parameters and errors on x2 and piecewise linear curve fitting to mutual information data.
The MI content between the northern geographic hemisphere polar cap index (PC(N)) and the AKR power is displayed in Fig. d. The PC(S) index is omitted to avoid interhemispheric asymmetries as a result of seasonal variations in ionospheric conductance. For both fits, the MI peaks at a negative lag ≥-5 min. For the slightly better-fitting quadratic fit, this suggests that the PC(N) index is excited 11 min before any corresponding enhancement in AKR power (the piecewise fit suggests τ=-5 min). This suggests that, 11 min before an excitation in the AKR power, a corresponding enhancement in the transfer of flux across the polar cap will be observed. This may relate to some portion of the building up of flux in the magnetotail preceding substorm onset, which is of the order of 30–90 min (e.g. ), and its related expansion onset has been shown to drive AKR enhancements. Alternatively, it could relate to the increase in AKR power in the 20 min prior to substorm onset shown by .
The ring current index SYM-H is compared with AKR intensity at different time lags and is presented in Fig. e. Although the best fit to these data is the quadratic curve (gold), the RMSE is similar for both fits. Compared with the other indices, there is greater variability across different time lags. The calculated peak varies from +3 (piecewise fit) to +9 (x2 fit) minutes, but both have large error values, which are detailed in Table . This may suggest that a shared driver between AKR and SYM-H may excite AKR emissions before the enhancement of the ring current (or indeed trigger a geomagnetic storm). It is important to note that the random phase surrogate threshold is only 1 order of magnitude below MI values for SYM-H.
Despite the large error values on the peak and the variability of the MI content across different lags, the MI content is above the threshold value set by the random phase surrogate and so may be interpreted as significant. Since AKR sources are generally centred at high latitudes (e.g. ) and the SYM-H index is measured by near-equatorial magnetometers, a strong relationship between the two was not expected. However, there may be more to investigate. For example, an assessment of AKR power with respect to storm phases detected in SYM-H (such as those presented by ) could allude to any relationship between AKR intensity and ring current activity.
Finally, bow shock IMF BZ is compared with AKR intensity, at lags ranging from -120 to +30 min, and is presented in Fig. f. Lags extending to more negative values are used because it is expected that changes in dayside IMF BZ may be considerably earlier than changes in AKR, a predominantly nightside phenomenon. The MI as a function of applied temporal lag is fitted best by the piecewise linear fit, with a peak at -53 min and a relatively tightly fitting confidence interval. This suggests that any change in IMF occurs some 53 min before any related AKR excitation or enhancement. This is in keeping with canonical timescales for propagation of dayside changes to nightside regions, i.e. substorm growth phases (e.g. ).
It is important to note at this point that, for all the parameters, the calculated MI values are below 0.05 for all the temporal lags. Although the significance of these values was tested using the random phase surrogate threshold, a discussion of potential reasons for these low values will be given here. Firstly, the geomagnetic indices and AKR intensity, although both measures of magnetospheric activity, sample different regions in the magnetosphere, particularly in terms of altitude. Additionally, due to the anisotropic beaming of AKR, its visibility and observed intensity also vary dramatically with latitude and local time (e.g. ). The AKR intensity data have not been corrected for this variation (as doing so is a non-trivial open question), which may contribute to low MI values. Finally, the measurement technique for AKR intensity and geomagnetic indices is very different: AKR intensity is observed by a single spacecraft moving through the magnetosphere, whereas the indices (excluding PC(N)) are averaged over multiple stations, so there may be some smoothing effects contributing to low MI values.
Application of a coupling timescale
The results presented in Sect. obtain a statistical estimate for the time lag between AKR intensity and geomagnetic indices. These were estimated over 10 years of AKR intensity data and as such represent a broad overall relationship between the two time series. To investigate this further, examples of substorms from November 2003 (where Wind was on the nightside of the planet and hence in the statistical AKR visibility region; e.g. ) were studied, and the relationship between AKR power and indices was examined. One such example is presented in Fig. . Substorm onsets were extracted from the Substorm Onsets and Phases from Indices of the Electrojet (SOPHIE, ) 75 % list, which determines substorm phases from percentiles of rates of change of SuperMAG AL (SML, equivalent to AL).
(a) Frequency–time–intensity spectrogram of W21-selected AKR data observed by the Wind satellite with time series of (b) AKR-integrated intensity between 100 and 650 kHz together with the (c) Auroral Electrojet Index (AE, black) and its upper (AU, purple) and lower (AL, green) envelopes. The green dotted line indicates the AL coupling timescale from the MI analysis (d) ring current index (SYM-H, black) and polar cap index (PC(N), pink). The pink dotted line indicates the PC(N) coupling timescale, and the grey shade indicates SYM-H coupling timescales from both fits. (e) IMF BZ in GSM coordinates; the orange dotted line indicates the coupling timescale for BZ. In all the panels, the black vertical dashed line denotes the substorm onset from the SOPHIE 75 % substorm list.
A frequency–time–intensity spectrogram of the W21-selected data is presented in Fig. a, with a black dashed line indicating substorm onset. AKR intensity begins to increase just before substorm onset (similarly to results by ), and this enhancement is also observed in a time series of AKR-integrated intensity presented in Fig. b. Additionally, the AKR band expands to lower frequencies shortly following substorm onset, although there is also some random noise uncorrelated in time and frequency around this low-frequency extension (LFE).
Focusing first on the auroral indices, AU, AL and AE are presented in purple, green and black respectively in Fig. c. AL shows a small but characteristic substorm signature starting with a rapid decrease shortly after the indicated substorm onset but simultaneously with the green dotted line. The green dotted line indicates the substorm onset plus the AL lag of the maximum MI (+8 min). This example shows that, although AKR is excited around substorm onset, the AL index begins its excitation around 8 min later. AE shows a similar enhancement dominated by the AL signature.
The polar cap index PC(N) and ring current index SYM-H are presented in pink and black respectively in Fig. d. The pink dotted line is drawn at substorm onset plus the PC(N) lag of maximum MI (-11 min). At the point at which the pink dotted line is drawn, PC(N) begins to reduce into trough. Recalling that MI assesses the relationship between variables independent of the order and direction of the relationship, this may be an indication of PC(N) decreasing about 11 min before an AKR intensity enhancement. The lag of the maximum MI from both fits (substorm onset +3 to +9 min) for SYM-H are presented as grey shade in Fig. d. This coupling indicator time period does not coincide with any obvious indicators of clear ring current activity and thus does not help to unravel the relationship between AKR and SYM-H.
Finally, IMF BZ is presented in orange in Fig. e. The orange dotted line indicates substorm onset (dashed black line) minus 53 min – the calculated lag of the maximum MI. This indicator coincides with a sharp southward IMF turning. This may relate to an enhancement in magnetic reconnection, leading to flux building up in the tail prior to substorm onset (and indeed AKR excitation).
Conclusions
In this study the MI content between AKR intensity and various relevant geomagnetic indices as well as IMF BZ at different temporal lags has been assessed. Both quadratic and piecewise linear functions were fitted to the MI as a function of lag data. The time lag at which maximum non-linear correlation occurs has been extracted and interpreted as the coupling timescale between AKR intensity and the phenomena represented by each geomagnetic index or BZ. It is important to note at this point that correlation does not necessarily imply causality. However, the wealth of literature on this subject suggests strong links between indices and AKR. The physical implication of this analysis and the key results are summarised below.
The AKR source region feels the effects of a shared driver before the auroral ionosphere.
This delay is more noticeable for westward electrojets (including the substorm electrojet) characterised by AL than eastward electrojets characterised by AU.
For AE/AU/AL the lag of the maximum MI is +7/+1/+8 min.
The PC(N) index is excited 5–11 min before any corresponding enhancement is observed in AKR.
This suggests an enhancement in anti-sunwards flux transport preceding AKR intensity: this confirms previous suggestions of flux building up in the magnetotail, leading to substorm onset, which excites AKR emission.
The relationship between temporal lag and mutual information between SYM-H and AKR intensity is less clear but may suggest that the AKR intensity feels the effect of a shared driver before the ring current.
IMF BZ changes 53 min before any related AKR excitation.
Further work to understand the coupling timescales between the AKR source region and the ionosphere could include a parameterisation of this analysis by different frequency regions within the AKR emission. Since emission frequency is inversely proportional to source altitude, this could provide more information on the propagation of driver effects along a magnetic field line. Additionally, further parametrisation by local time could help to explain viewing effects on these coupling timescales.
Code and data availability
Wind/WAVES data that have been empirically selected for AKR emissions using the technique by and a subset are available online (10.25935/wxv0-vr90, ). Geomagnetic indices AE, AU, AL, PC(N) and SYM-H as well as BZ were obtained via OMNIWeb (10.48322/45bb-8792, , https://omniweb.gsfc.nasa.gov/hw.html, last access: 28 January 2021, ). We gratefully acknowledge use of NASA/GSFC's Space Physics Data Facility's OMNIWeb service and OMNI data. The AU, AL and SYM-H indices used in this paper were provided by the WDC for Geomagnetism, Kyoto (https://doi.org/10.17593/15031-54800, ), via OMNIWeb. The PC(N) index was provided by the World Data Center for Geomagnetism, Copenhagen, via OMNIWeb. We acknowledge the use of Python packages scikit-learn , aaft (https://github.com/lneisenman/aaft, ) and scipy (https://docs.scipy.org/doc/scipy/index.html, last access: 16 July 2021). Code to estimate the coupling timescales using mutual information is available from (https://github.com/arfogg/generic_MI_lag_finder, last access: 30 June 2023; DOI: 10.5281/zenodo.10804655, ).
Author contributions
ARF, CMJ and SCC designed the method. ARF coded the method with testing from AB. ARF created all the plots with input from CMJ and SCC. JEW processed the data using the W21 technique. ARF prepared the manuscript with contributions from all the co-authors.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Acknowledgements
Alexandra Ruth Fogg's work was supported by Irish Research Council Government of Ireland Postdoctoral Fellowship GOIPD/2022/782 and Science Foundation Ireland Grant 18/FRL/6199. Caitríona M. Jackman's work was supported by Science Foundation Ireland Grant 18/FRL/6199. Sandra C. Chapman acknowledges AFOSR grant FA8655-22-7056 and an ISSI Johannes Geiss Fellowship. James E. Waters' work was supported by the EPSRC Centre for Doctoral Training in Next Generation Computational Modelling grant no. EP/L015382/1. The authors acknowledge the CNES (Centre National d'Études Spatiales) and CNRS (Centre National de la Recherche Scientifique)/INSU (Institut National des Sciences de l'Univers) programmes of planetology and heliophysics and the Observatoire de Paris for support to the Wind/WAVES team, together with the CDPP (Centre de Données de la Physique des Plasmas) for the provision of the Wind/WAVES RAD1 L2 data.
Financial support
This research has been supported by the Science Foundation Ireland (grant no. 18/FRL/6199), the Irish Research Council (grant no. GOIPD/2022/782), the Air Force Office of Scientific Research (grant no. FA8655-22-7056) and the Engineering and Physical Sciences Research Council (grant no. EP/L015382/1).
Review statement
This paper was edited by Bruce Tsurutani and reviewed by Rajkumar Hajra and one anonymous referee.
ReferencesAraki, T.: A physical model of the geomagnetic sudden commencement, Geoph. Monog. Series, 81, 183–200, 10.1029/GM081p0183, 1994.Benson, R. F. and Calvert, W.: ISIS 1 observations at the source of auroral kilometric radiation, Geophys. Res. Lett., 6, 479–482, 10.1029/GL006i006p00479, 1979.Benson, R. F., Calvert, W., and Klumpar, D. M.: Simultaneous wave and particle observations in the auroral kilometric radiation source region, Geophys. Res. Lett., 7, 959–962, 10.1029/GL007i011p00959, 1980.Bougeret, J. L., Kaiser, M. L., Kellogg, P. J., Manning, R., Goetz, K., Monson, S. J., Monge, N., Friel, L., Meetre, C. A., Perche, C., Sitruk, L., and Hoang, S.: WAVES: The radio and plasma wave investigation on the Wind spacecraft, Space Sci. Rev., 71, 231–263, 10.1007/BF00751331, 1995.Calvert, W.: The auroral plasma cavity, Geophys. Res. Lett., 8, 919–921, 10.1029/GL008i008p00919, 1981. Cowley, S. W. H.: Magnetosphere-ionosphere interactions: A tutorial review, Magnetospheric Current Systems, Geoph. Monog. Series, 118, 91–106, 2000.Davis, T. N. and Sugiura, M.: Auroral electrojet activity index AE and its universal time variations, J. Geophys. Res., 71, 785–801, 10.1029/JZ071i003p00785, 1966.Desch, M. D., Kaiser, M. L., and Farrell, W. M.: Control of terrestrial low frequency bursts by solar wind speed, Geophys. Res. Lett., 23, 1251–1254, 10.1029/96GL01352, 1996.Dunckel, N., Ficklin, B., Rorden, L., and Helliwell, R. A.: Low-Frequency Noise Observed in the Distant Magnetosphere with OGO 1, J. Geophys. Res.-Space, 75, 1854–1862, 10.1029/JA075i010p01854, 1970.Dungey, J. W.: Interplanetary magnetic field and the auroral zones, Phys. Rev. Lett., 6, 47, 10.1103/PhysRevLett.6.47, 1961.Ergun, R. E., Carlson, C. W., McFadden, J. P., Mozer, F. S., Delory, G. T., Peria, W., Chaston, C. C., Temerin, M., Elphic, R., Strangeway, R., Pfaff, R., Cattell, C. A., Klumpar, D., Shelly, E., Peterson, W., Moebius, E., and Kistley, L.: FAST satellite wave observations in the AKR source region, Geophys. Res. Lett., 25, 2061–2064, 10.1029/98GL00570, 1998.Fogg, A. R.: Generic MI lag finder, GitHub [code], https://github.com/arfogg/generic_MI_lag_finder, last access: 3 July 2023.Fogg, A. R.: arfogg/generic_MI_lag_finder: First release of the mutual information lag finding tool (v1.0.0), Zenodo [code], 10.5281/zenodo.10804655, 2024.Fogg, A. R., Lester, M., Yeoman, T. K., Burrell, A. G., Imber, S. M., Milan, S. E., Thomas, E. G., Sangha, H., and Anderson, B. J.: An Improved Estimation of SuperDARN Heppner-Maynard Boundaries using AMPERE data, J. Geophys. Res.-Space, 125, e2019JA027218, 10.1029/2019JA027218, 2020.Fogg, A. R., Jackman, C. M., Waters, J. E., Bonnin, X., Lamy, L., Cecconi, B., Issautier, K., and Louis, C. K.: Wind/WAVES Observations of Auroral Kilometric Radiation: Automated Burst Detection and Terrestrial Solar Wind – Magnetosphere Coupling Effects, J. Geophys. Res.-Space, 127, e2021JA030209, 10.1029/2021JA030209, 2022.Fogg, A. R., Lester, M., Yeoman, T. K., Carter, J. A., Milan, S. E., Sangha, H. K., Elsden, T., Wharton, S. J., James, M. K., Malone-Leigh, J., Paxton, L. J., Anderson, B. J., and Vines, S. K.: Multi-instrument observations of the effects of a solar wind pressure pulse on the high latitude ionosphere: a detailed case study of a geomagnetic sudden impulse, J. Geophys. Res.-Space, 128, e2022JA031136, 10.1029/2022JA031136, 2023.Fogg, A. R., Jackman, C. M., Coco, I., Douglas Rooney, L., Weigt, D. M., and Lester, M.: Why are some solar wind pressure pulses followed by geomagnetic storms?, J. Geophys. Res.-Space, 128, e2022JA031259, 10.1029/2022JA031259, 2024.Forsyth, C., Rae, I. J., Coxon, J. C., Freeman, M. P., Jackman, C. M., Gjerloev, J., and Fazakerley, A. N.: A new technique for determining Substorm Onsets and Phases from Indices of the Electrojet (SOPHIE), J. Geophys. Res.-Space, 120, 10592–10606, 10.1002/2015JA021343, 2015.Gallagher, D. L. and D'Angelo, N.: Correlations between solar wind parameters and auroral kilometric radiation intensity, Geophys. Res. Lett., 8, 1087–1089, 10.1029/GL008i010p01087, 1981.Green, J. L. and Gurnett, D. A.: A Correlation Between Auroral Kilometric Radiation and Inverted V Electron Precipitation, J. Geophys. Res., 84, 5216–5222, 10.1029/JA084iA09p05216, 1979.Green, J. L., Gurnett, D. A., and Shawhan, S. D.: The Angular Distribution of Auroral Kilometric Radiation, J. Geophys. Res., 82, 1825–1838, 10.1029/JA082i013p01825, 1977.Green, J. L., Boardsen, S., Garcia, L., Fung, S. F., and Reinisch, B. W.: Seasonal and solar cycle dynamics of the auroral kilometric radiation source region, J. Geophys. Res., 109, A05223, 10.1029/2003JA010311, 2003.Gurnett, D. A.: The Earth as a Radio Source: Terrestrial Kilometric Radiation, J. Geophys. Res., 79, 4227–4238, 10.1029/JA079i028p04227, 1974.Hashimoto, K., Matsumoto, H., Murata, T., Kaiser, M. L., and Bougeret, J.-L.: Comparison of AKR simultaneously observed by the GEOTAIL and WIND spacecraft, Geophys. Res. Lett., 25, 853–856, 10.1029/98GL00385, 1998.Hilgers, A.: The auroral radiating plasma cavities, Geophys. Res. Lett., 19, 237–240, 10.1029/91GL02938, 1992.Huff, R. L., Calvert, W., Craven, J. D., Frank, L. A., and Gurnett, D. A.: Mapping of Auroral Kilometric Radiation Sources to the Aurora, J. Geophys. Res., 93, 11445–11454, 10.1029/JA093iA10p11445, 1988.Iyemori, T.: Storm-Time Magnetospheric Currents Inferred from Mid-Latitude Geomagnetic Field Variations, J. Geomagn. Geoelectr., 42, 1249–1265, 10.5636/jgg.42.1249, 1990.Johnson, J. R. and Wing, S.: External versus internal triggering of substorms: An information-theoretical approach, Geophys. Res. Lett., 41, 5748–5754, 10.1002/2014GL060928, 2014.Johnson, M. T., Wygant, J. R., Cattell, C., Mozer, F. S., Temerin, M., and Scudder, J.: Observations of the seasonal dependence of the thermal plasma density in the Southern Hemisphere auroral zone and polar cap at 1 RE, J. Geophys. Res., 106, 19023–19033, 10.1029/2000JA900147, 2001.Kasaba, Y., Matsumoto, H., and Hashimoto, K.: The angular distribution of auroral kilometric radiation observed by the GEOTAIL spacecraft, Geophys. Res. Lett., 24, 2483–2486, 10.1029/97GL02599, 1997. Kozachenko, L. F. and Leonenko, N. N.: Sample estimate of the entropy of a random vector, Problemy Peredachi Informatsii, 23, 9–16, 1987.Kraskov, A., Stögbauer, H., and Grassberger, P.: Estimating mutual information, Phys. Rev. E, 69, 066138, 10.1103/PhysRevE.69.066138, 2004.Kurth, W. S., Murata, T., Lu, G., Gurnett, D. A., and Matsumoto, H.: Auroral kilometric radiation and the auroral electrojet index for the January 1997 magnetic cloud event, Geophys. Res. Lett., 25, 3027–3030, 10.1029/98GL00404, 1998.Lamy, L., Zarka, P., Cecconi, B., Prangé, R., Kurth, W. S., and Gurnett, D. A.: Saturn kilometric radiation: Average and statistical properties, J. Geophys. Res., 113, A07201, 10.1029/2007JA012900, 2008.Lamy, L., Zarka, P., Cecconi, B., and Prangé, R.: Auroral kilometric radiation diurnal, semidiurnal, and shorter term modulations disentangled by Cassini, J. Geophys. Res., 115, A09221, 10.1029/2010JA015434, 2010.Li, H., Wang, C., and Peng, Z.: Solar wind impacts on growth phase duration and substorm intensity: A statistical approach, J. Geophys. Res.-Space, 118, 4270–4278, 10.1002/jgra.50399, 2013.Lneiseman: aaft, GitHub [code], https://github.com/lneisenman/aaft, last access: 16 July 2021. MacKay, D. J. C.: Information Theory, Inference, and Learning Algorithms, Cambridge University Press, ISBN 9780521642989, 2003.March, T. K., Chapman, S. C., and Dendy, R. O.: Mutual information between geomagnetic indices and the solar wind as seen by WIND: Implications for propagation time estimates, Geophys. Res. Lett., 32, L04101, 10.1029/2004GL021677, 2005.Menietti, J. D., Mutel, R. L., Christopher, I. W., Hutchinson, K. A., and Sigwarth, J. B.: Simultaneous radio and optical observations of auroral structures: Implications for AKR beaming, J. Geophys. Res., 116, A12219, 10.1029/2011JA017168, 2011.Milan, S. E., Clausen, L. B. N., Coxon, J. C., Carter, J. A., Walach, M.-T., Laundal, K., Østgaard, N., Tenfjord, P., Reistad, J., Snekvik, K., Korth, H., and Anderson, B. J.: Overview of solar wind–magnetosphere–ionosphere–atmosphere coupling and the generation of magnetospheric currents, Space Sci. Rev., 206, 547–573, 10.1007/s11214-017-0333-0, 2017.Morioka, A., Miyoshi, Y., Tsuchiya, F., Misawa, H., Sakanoi, T., Yumoto, K., Anderson, R. R., Menietti, J. D., and Donovan, E. F.: Dual structure of auroral acceleration regions at substorm onsets as derived from auroral kilometric radiation spectra, J. Geophys. Res., 112, A06245, 10.1029/2006JA012186, 2007.Morioka, A., Miyoshi, Y., Tsuchiya, F., Misawa, H., Kasaba, Y., Asozu, T., Okano, S., Kadokura, A., Sato, N. Miyaoka, H., Yumoto, K., Parks, G. K., Honary, F., Trotignon, J. G., Décréau, P. M. E., and W., R. B.: On the simultaneity of substorm onset between two hemispheres, J. Geophys. Res., 116, A04211, 10.1029/2010JA016174, 2011.Morioka, A., Miyoshi, Y., Kurita, S., Kasaba, Y., Angelopoulos, V., Misawa, H., Kojima, H., and McFadden, J. P.: Universal time control of AKR: Earth as a spin-modulated variable radio source, J. Geophys. Res.-Space, 118, 1123–1131, 10.1002/jgra.50180, 2013.Morioka, A., Miyoshi, Y., Kasaba, Y., Sato, N., Kadokura, A., H., M., Miyashita, Y., and Mann, I.: Substorm onset process: Ignition of auroral acceleration and related substorm phases, J. Geophys. Res.-Space, 119, 1044–1059, 10.1002/2013JA019442, 2014.Mutel, R. L., Christopher, I. W., and Pickett, J. S.: Cluster multispacecraft determination of AKR angular beaming, Geophys. Res. Lett., 35, L07104, 10.1029/2008GL033377, 2008.Panchenko, M., Khodachenko, M. L., Kislyakov, A. G., Rucker, H. O., J., H., Kaiser, M. L., Bale, S. D., Lamy, L., Cecconi, B., Zarka, P., and Goetz, K.: Daily variations of auroral kilometric radiation observed by STEREO, Geophys. Res. Lett., 36, L06102, 10.1029/2008GL037042, 2009.Papitashvili, N.: OMNIWeb Plus, NASA [data set], https://omniweb.gsfc.nasa.gov/hw.html (last access: 28 January 2021), 2023.Papitashvili, N. E. and King, J. H.: OMNI 1-min Data, NASA Space Physics Data Facility [data set], 10.48322/45bb-8792, 2020. Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., Blondel, M., Prettenhofer, P., Weiss, R., Dubourg, V., Vanderplas, J., Passos, A., Cournapeau, D., Brucher, M., Perrot, M., and Duchesnay, E.: Scikit-learn: Machine Learning in Python, J. Mach. Learn. Res., 12, 2825–2830, 2011.Schreiber, T. and Schmitz, A.: Improved surrogate data for nonlinearity tests, Phys. Rev. Lett., 77, 635–638, 10.1103/PhysRevLett.77.635, 1996.Shannon, C. E.: A Mathematical Theory of Communication, Bell Syst. Tech. J., 27, 379–423, 10.1002/j.1538-7305.1948.tb01338.x, 1949.Snelling, J. M., Johnson, J. R., Willard, J., Nurhan, Y., Homan, J., and Wing, S.: Information Theoretical Approach to Understanding Flare Waiting Times, Astrophys. J., 899, 148, 10.3847/1538-4357/aba7b9, 2020.Stauning, P.: The Polar Cap index: A critical review of methods and a new approach, J. Geophys. Res.-Space, 118, 5021–5038, 10.1002/jgra.50462, 2013.Tindale, E., Chapman, S. C., Moloney, N. R., and Watkins, N. W.: The Dependence of Solar Wind Burst Size on Burst Duration and Its Invariance Across Solar Cycles 23 and 24, J. Geophys. Res.-Space, 123, 7196–7210, 10.1029/2018JA025740, 2018.Treumann, R. A. and Baumjohann, W.: Auroral Kilometric Radiation and Electron Pairing, Frontiers in Physics, 8, 386, 10.3389/fphy.2020.00386, 2020.Troshichev, O. A. and Andrezen, V. G.: The relationship between interplanetary quantities and magnetic activity in the southern polar cap, Planetary Space Science, 33, 415–419, 10.1016/0032-0633(85)90086-8, 1985.Troshichev, O. A., Dolgacheva, S., Stepanov, N. A., and Sormakov, D. A.: The PC index variations during 23/24 solar cycles: Relation to solar wind parameters and magnetospheric disturbances, J. Geophys. Res.-Space, 126, e2020JA028491, 10.1029/2020JA028491, 2021. Vaughan, S.: Scientific Inference: Learning from data, Cambridge University Press, ISBN 9781139176071, 2013.Voots, G. R., Gurnett, D. A., and Akasofu, S. I.: Auroral Kilometric Radiation as an Indicator of Auroral Magnetic Disturbances, J. Geophys. Res., 82, 2259–2266, 10.1029/JA082i016p02259, 1977.Walach, M.-T. and Grocott, A.: SuperDARN Observations During Geomagnetic Storms, Geomagnetically Active Times and Enhanced Solar Wind Driving, J. Geophys. Res.-Space, 124, 5828–5847, 10.1029/2019JA026816, 2019.Waters, J. E., Cecconi, B., Bonnin, X., and Lamy, L.: Wind/Waves flux density collection calibrated for Auroral Kilometric Radiation (Version 1.0), MASER [data set], 10.25935/wxv0-vr90, 2021a.Waters, J. E., Jackman, C. M., Lamy, L., Cecconi, B., Whiter, D., Bonnin, X., Issautier, K., and Fogg, A. R.: Empirical Selection of Auroral Kilometric Radiation During a Multipoint Remote Observation With Wind and Cassini, J. Geophys. Res.-Space, 126, e2021JA029425, 10.1029/2021JA029425, 2021b.Waters, J. E., Jackman, C. M., Whiter, D., Forsyth, C., Fogg, A. R., Lamy, L., Cecconi, B., Bonnin, X., and Issautier, K.: A perspective on substorm dynamics using 10 years of Auroral Kilometric Radiation observations from Wind, J. Geophys. Res.-Space, 127, e2022JA030449, 10.1029/2022JA030449, 2022.Wicks, R. T., Chapman, S. C., and Dendy, R. O.: Spatial correlation of solar wind fluctuations and their solar cycle dependence, Astrophys. J., 690, 734–742, 10.1088/0004-637X/690/1/734, 2009.Wilson III, L. B., Brosius, A. L., Gopalswamy, N., Nieves-Chinchilla, T., Szabo, A., Hurley, K., Phan, T., Kasper, J. C., Lugaz, N., Richardson, I. G., Chen, C. H. K., D., V., Wicks, R. T., and TenBarge, J. M.: A Quarter Century of Wind Spacecraft Discoveries, Rev. Geophys., 59, e2020RG000714, 10.1029/2020RG000714, 2021.Wing, S., Brandt, P. C., Mitchell, D. G., Johnson, J. R., Kurth, W. S., and Menietti, J. D.: Periodic Narrowband Radio Wave Emissions and Inward Plasma Transport at Saturn's Magnetosphere, Astron. J., 159, 249, 10.3847/1538-3881/ab818d, 2020. World Data Center for Geomagnetism Kyoto, Nose, M., Iyemori, T., Sugiura, M., and Kamei, T.: Geomagnetic AE index [data set], 10.17593/15031-54800, 2015.Wu, C. S.: Kinetic cyclotron and synchrotron maser instabilities: radio emission processes by direct amplification of radiation, Space Sci. Rev., 41, 215–298, 10.1007/BF00190653, 1985.Wu, C. S. and Lee, L. C.: A theory of the terrestrial kilometric radiation, Astrophys. J., 230, 621–626, 10.1086/157120, 1979.Zhao, W., Liu, S., Zhang, S., Zhou, Q., Yang, C., and He, Y.: Global Occurrences of Auroral Kilometric Radiation Related to Suprathermal Electrons in Radiation Belts, Geophys. Res. Lett., 46, 7230–7236, 10.1029/2019GL083944, 2019.