Articles | Volume 40, issue 4
Ann. Geophys., 40, 433–461, 2022
Ann. Geophys., 40, 433–461, 2022
Regular paper
04 Jul 2022
Regular paper | 04 Jul 2022

Diagnostic study of geomagnetic storm-induced ionospheric changes over very low-frequency signal propagation paths in the mid-latitude D region

Diagnostic study of geomagnetic storm-induced ionospheric changes over very low-frequency signal propagation paths in the mid-latitude D region
Victor U. J. Nwankwo1, William Denig2, Sandip K. Chakrabarti3, Olugbenga Ogunmodimu4,5, Muyiwa P. Ajakaiye1, Johnson O. Fatokun1, Paul I. Anekwe1, Omodara E. Obisesan1, Olufemi E. Oyanameh1, and Oluwaseun V. Fatoye1 Victor U. J. Nwankwo et al.
  • 1Space, Atmospheric Physics and Radio wave Propagation Laboratory, Anchor University, Lagos, Nigeria
  • 2Department of Sciences, St. Joseph's College of Maine, Standish, ME 04084, USA
  • 3Indian Centre for Space Physics, Kolkata 700084, India
  • 4Department of Electrical Engineering, Manchester Metropolitan University, Manchester, UK
  • 5Department of Electrical Engineering, University of Bolton, Bolton, UK

Correspondence: Victor U. J. Nwankwo (


We performed a diagnostic study of geomagnetic storm-induced disturbances that are coupled to the mid-latitude D region by quantifying the propagation characteristics of very low-frequency (VLF) radio signals from transmitters located in Cumbria, UK (call sign GQD), and Rhauderfehn, Germany (DHO), and received in southern France (A118). We characterised the diurnal VLF amplitudes from two propagation paths into five metrics, namely the mean amplitude before sunrise (MBSR), the midday amplitude peak (MDP), the mean amplitude after sunset (MASS), the sunrise terminator (SRT) and the sunset terminator (SST). We analysed and monitored trends in the variation of signal metrics for up to 20 storms to relate the deviations in the signal amplitudes that were attributable to the storms. Five storms and their effects on the signals were examined in further detail. Our results indicate that relative to pre-storm levels the storm day MDP exhibited characteristic decreases in about 80 % (67 %) of the events for the DHO-A118 (GQD-A118) propagation path. The MBSR showed decreases of about 60 % (77 %), whereas the MASS decreased by 67 % (58 %). Conversely, the SRT and SST showed amplitude decreases of 33 % (25 %) and 47 % (42 %), respectively. Of the two propagation paths, the amplitude decreases for the DHO-A118 propagation path signal were greater, as previously noted by Nwankwo et al. (2016). To better understand the state of the ionosphere over the signal propagation paths and how it might have affected the VLF amplitudes, we further analysed the virtual heights (h'E, h'F1 and h'F2) and critical frequencies (foE, foF1 and foF2) from ionosondes located near the transmitters. The results of this analysis showed significant increases and fluctuations in both the F-region critical frequencies and virtual heights during the geomagnetic storms. The largest increases in the virtual heights occurred near the DHO transmitter in Rhauderfehn (Germany), suggesting a strong storm response over the region which might account for the larger MDP decrease along the DHO-A118 propagation path.

1 Introduction

The terrestrial magnetosphere is formed by the interaction between the solar wind and the Earth's magnetic field (McPherron et al.2008). In contrast, the ionosphere is largely the result of solar photoionisation of the neutral atmosphere balanced against chemical recombination and particle transport (Prolss2004; Kelley2009). The magnetosphere and ionosphere are coupled via the geomagnetic field, effectively tying these seemingly disparate regions into a global magnetosphere–ionosphere (M–I) system (Blanc1988; Nwankwo et al.2016). Within the concept of an open magnetosphere (Dungey1961), energy is transferred from the solar wind to the M–I system via magnetic reconnection (Cassak2016). Variations within the interplanetary environment, driven by solar disturbances, affect the M–I system, particularly when the interplanetary magnetic field (IMF), embedded within the solar wind, is oriented southward relative to the outer northward-directed geomagnetic field (Gonzalez et al.1999; Liu and Li2002). A geomagnetic storm is triggered when the solar–terrestrial interaction is sufficiently intense to energise the ring current (Jordanova et al.2020) and to solicit a negative response in ground-based magnetometers in terms of the disturbance storm time (Dst) index (Russel et al.1974; Mayaud1980; Borovsky and Shprits2017). The strength of a geomagnetic storm is typically classified as small, moderate or intense for Dst values less than (<) −30, <-50 and <-100 nT, respectively (Gonzalez et al.1994). The two leading drivers of geomagnetic storms are coronal mass ejections (CMEs) and corotating interactive regions (CIRs) (Baker2000). A CME is the impulsive release of solar material into interplanetary space from a solar active region that may, or may not, be associated with a solar flare (Youssef2012). Conversely, a CIR is the result of a high-speed stream (HSS) emitted from a solar coronal hole overtaking the background solar wind (Choi et al.2009). In both cases the shock front at the leading edge of the interplanetary disturbance increases the ram pressure imposed on the dayside magnetopause, causing a reconfiguration of the dayside Chapman–Ferraro current system (Chapman and Ferraro1930) and, in turn, the magnetosphere as a whole (Ganushkina et al.2018). A typical geomagnetic storm has three phases consisting of (1) a sudden storm commencement (SSC) at the time of the increased ram pressure, (2) a main phase as the ring current is energised and (3) a recovery phase as the magnetosphere returns to a more quiescent ground state (Akasofu2018; Gonzalez et al.1994). There are distinct differences between how and when CME-driven versus CIR-driven storms affect the Earth (Borovsky and Denton2006). Typically, CME-driven storms are stronger with regards to Dst (Tsurutani et al.2006) and occur more frequently near the maximum in the Sun's 11-year activity cycle (Gonzalez et al.1999), whereas CIR-driven storms are weaker and occur predominantly on the decreasing phase of solar activity (Tsurutani et al.1995). While the impacts of both CME- and CIR-driven geomagnetic storms on the middle to upper atmosphere have been extensively studied and are well known (Fuller-Rowell et al.1994; Burch2016; Heelis and Maute2020), less certain are the geomagnetic storm effects in the lower ionosphere (Lastovicka1996; Kumar and Kumar2014). For the purposes of this report our attention is limited to CME-driven, moderate geomagnetic storms and the resulting impacts on D-region very low-frequency (VLF) radio-wave propagation within the Earth–ionosphere wave guide (EIWG). The VLF band spans the range from 3 kHz to 30 kHz.

The discovery of the ionosphere is generally attributed to Carl Gauss, who, in 1838, speculated on the existence of an ionised atmospheric layer to explain variations in the measured geomagnetic field (Gauss1838; Glassmeier and Tsurutani2014). Heaviside (1902) and Kennelly (1902) both independently speculated on the existence of an electrically reflecting layer following Guglielmo Marconi's demonstration (Marconi1901) of long-range radio-wave transmission. Appleton and Barnett (1925a, b) were the first to prove the existence of this “electrical layer”, or simply E layer, which they determined was located approximately 80 to 90 km above the Earth's surface. The E layer corresponded to a local maximum in the vertical electron density profile which was effective at reflecting radio waves of up to several megahertz (Beynon1969). Measurements of reflecting layers both above (Appleton1927) and below (Colwell and Friend1936) the E layer, or more generally E region, soon followed, which were non-coincidently named the F region and D region, respectively. Further studies soon found a separation in the F region that was then categorised as the F1 and F2 regions (Appleton and Naismith1935). These layers are the primary regions which comprise the “ionosphere”, a term that was first coined by Robert Watson-Watt in 1926 and broadly adopted into the scientific lexicon after 1929 (Gardiner1969; Appleton1933). It has been long recognised that the ionosphere has both deleterious and enabling effects on terrestrial radio-wave communications, dependent on the transmission frequency band in use (Poole1999; Bennington1944). The D region is uniquely relevant to VLF communication as a means of communication with submarines and related underwater vehicles (Moore1967; Waheed and Yousufzai2011; Lanzagorta2012; Sun et al.2021).

The nominal altitudes for the D- (daytime), E-, F1- and F2-region peaks are 60, 110, 170 and 300 km, respectively, although there is considerable variability within each layer depending on the solar-geophysical conditions (Mangla and Yadov2011). The various ionospheric layers are the result of ionisation production and loss. Lodge (1902) was the first to suggest that photoionisation of the background neutral atmosphere by solar ultraviolet (UV) radiation was responsible for establishing Marconi's electrically reflecting layer. Chapman (1931) codified this primary production source for the E and F layers. Solar UV radiation cannot effectively penetrate to altitudes below about 100 km (Marr1965). An exception is Lyman-α radiation at a wavelength of 121.5 nm within the atmospheric window of low absorption (Machol et al.2019) and is largely responsible for maintaining the quiescent dayside D region (Nicolet and Aikin1960). The complex chemical processes maintaining the ionospheric layers were initially discussed by Bates and Massey (1946, 1947) and more recently by numerous authors (Rishbeth1973; Schunk and Nagy2009; Pavlov2012). Ionospheric loss processes are mostly recombination, charge exchange and diffusion (Banks and Kockarts1973). Above about 150 km the dominant ion species is atomic oxygen, whereas molecular ions are more abundant at the lower altitudes (Johnson1966). Consequently, dissociative recombination is the dominant loss process for the D and E regions, whereas radiative recombination and charge exchange (followed by dissociative recombination) are more prevalent closer to the F2 peak (Mangla and Yadov2011). Within the topside ionosphere, diffusion becomes the dominant loss process (Rishbeth1973). Solar flares emit photons across the electromagnetic spectrum, producing a transient ionisation source affecting the delicate balance of ionospheric production and loss (Mitra1974). X-rays, as well as gamma rays, emitted by solar flares are sufficiently energetic to penetrate the lower ionosphere and briefly dominate the ionisation production rate of the dayside D region (Hayes et al.2021). At night the D-region electron density is greatly diminished due to the loss of solar radiant photons, whilst diffusive recombination continues unabated such that the D region seamlessly coalesces into the lower E region at about 90 km in the absence of other significant ionisation sources (Thomas et al.2007; Thomson and McRae2009). Prior to the space age, the detection of sudden ionospheric disturbances (SIDs) in the amplitude and phase of VLF radio-wave transmissions within the D region was used as an established proxy technique for monitoring the occurrence of solar flares (Lincoln1964; Moral et al.2013; Hegde et al.2018), which were known to have a deleterious impact on radio-wave communications (Sauer and Wilkinson2008; Dellinger1937). The widely accepted standard for specifying the ionospheric electron density profile (EDP) is the empirically based International Reference Ionosphere (IRI) model which has evolved over time (Rawer et al.1978; Rawer1981; Bilitza1990, 2001; Bilitza and Reinisch2008; Bilitza2018), with special consideration given to the D region of the lower ionosphere (Bilitza1981; Friedrich and Torkar1992; Danilov and Smirnova1995; Bilitza1998). An alternative approach for specifying the D region is the use of specialised atmospheric models, such as the Whole Atmosphere Community Climate Model with D-region ion chemistry (WACCM-D), which are focused on the neutral atmosphere (Verronen et al.2016; Andersson et al.2016; Siskind et al.2017), which, of course, is coupled to the ionised atmosphere via chemistry (Turunen et al.1996; Schunk1996, 1999; Verronen et al.2005; Turunen et al.2009; Kovacs et al.2016; Verronen et al.2016; Turunen et al.2016; Miyoshi et al.2021) within the overall M–I system.

Techniques used to probe the ionosphere include both ground-based and space-based approaches. The earliest methods used by Appleton and Barnett (1925b, 1926) in short-range transmitter-to-receiver trials differentiated the reflected sky wave from the direct ground wave to determine the height of the E layer. About the same time, Merve Tuve and Gregory Breit (Tuve and Breit1925; Breit and Tuve1925) proposed a methodology of using pulsed radio-wave transmissions for measuring the heights of overhead reflecting layers. The Tuve–Breit methodology was the basis for ionosondes and for many years was the preeminent scientific technique for “sounding” the ionosphere (Bibl1998). Ionosondes typically operate in the high-frequency (HF) domain from 3 to 30 MHz to derive the vertical ionisation profiles of the E and F regions. Ionosondes are quite affordable, leading to their widespread use in ionospheric characterisation (Stamper et al.2005). However, a key limitation is that ground-based ionosondes can effectively measure only the “bottom-side” ionosphere at and below the F-region peaks (Reinisch and Xueqin1983). A variant is the space-based “topside sounder” approach from which measurements of the topside ionosphere, above the F2 peak, can be obtained (Chapman and Warren1968; Benson2010). Ionospheric profiles of the E and F regions can also be obtained using incoherent scatter radars (ISRs) (Robinson et al.2009), which are typically operated at several hundred megahertz within the ultra-high-frequency (UHF) range (Belehaki et al.2016) and rely on the principle of Thomson electron scattering (Farley et al.1961; Dougherty and Farley1961, 1963). Conversely, coherent scatter HF radars (Greenwald et al.1995) rely on Bragg scattering (Takefu1989) from ionospheric structures to probe the ionosphere (Belehaki et al.2016). Examples of these advanced ionospheric measuring technologies include the Alouette topside sounder (Jackson1986), the Millstone Hill ISR (Evans1969a, b) and the Super Dual Auroral (HF) Radar Network (SuperDARN) (Greenwald2021). More recently, researchers have leveraged the outstanding capabilities of the Global Navigation Satellite System (GNSS), initially the Global Positioning System (GPS), for ionospheric characterisation (Davies and Hartmann1997) using receivers on the ground (Mannucci et al.1998; Prol et al.2021) and in space (Mannucci et al.2020). The Continuously Operating Reference Stations (CORS) (Snay and Soler2008) is a good example of a ground-based GNSS network, whereas the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) (Yue et al2014) is an example of the related space-based approach. With regard to the lower E and D regions, none of the aforementioned technologies is particularly effective at monitoring the bottommost ionospheric region (Cummer et al.1998; Kumar et al.2015), which is the focus of the present work. However, the perspective provided by considering the state of the local ionosphere allows us to assess our findings within the context of a coupled M–I system.

Reliance on VLF waves has proven to be an effective tool for monitoring and characterising the lower ionosphere (Sechrist1974; Inan et al.2010; Gross and Cohen2020). Early research within the VLF band was focused on lightning-induced “sferics” (Pierce1969), a technique that formed the basis for modern lightning detection networks (Betz et al.2009). VLF sferics can travel significant distances within the EIWG (Wait and Spies1964) and be detected as sound “tweeks” due to the dispersive nature of the ionosphere (Singh et al.2016). D-region characteristics can be derived from these events, although some care is required to account for their sporadic nature, both temporally and spatially (McCormick and Morris2018). Conversely, controlled experiments using known VLF frequencies and transmitter–receiver great circle paths (TRGCPs) can be used to characterise the lower ionosphere nominally to and from fixed locations (Kumar and Kumar2020). Well-known features in VLF TRGCP propagation are diurnal variations in amplitude and phase (Yokoyama and Tanimura1933; Pierce1955; Taylor1960; Chilton et al.1964; Lynn1978; Barr et al.2000; McRae and Thomson2000; Sharma and More2017) and characteristic signatures of sunrise and sunset (Walker1965; Crombie1964, 1966; Ries1967; Samanes et al.2015; Sharma and More2017; Gu and Xu2020). These features can be explained within the context of the wave-mode theory of VLF propagation attributed to Budden (1951, 1953, 1957) and promulgated by Wait (1960, 1961, 1963, 1964, 1968, 1970) and collaborators (Spies and Wait1961; Wait and Spies1964). While the transient response of VLF TRGCP propagation to solar flares has been well documented (Mitra1974; Thomson and Clilverd2001; McRae and Thomson2004; Abd Rashid et al.2013; Palit et al.2013; Rozhnoi et al.2019), the associated characterisation to geomagnetic storms has been less quantified due, perhaps, to the mixed ionospheric responses in the lower atmosphere from geoeffective CMEs and CIRs (Turner et al. 2006; Kim et al.2008; Laughlin et al.2008; Verbanac et al.2011; Soni et al.2020). The focus of this effort is then to augment the limited body of research related to the impact of geomagnetic storms on VLF propagation along TRGCPs within the EIWG (Tatsuta et al.2015).

The ionospheric response to geomagnetic storms is varied, and interpreting the response in terms of a regional or global specification of electron density often requires the use of sophisticated environmental models (Schunk et al.2004; Immel and Mannucci2013; Greer et al.2017). However, these models are mostly focused on the upper ionosphere, E region and above, and inclusion of the D region involves separate modules that, in turn, are mostly focused on the prompt D-region response to solar transient events in the forms of flares and related energetic particle events (Bilitza1998; Eccles et al.2005; Sauer and Wilkinson2008; Rogers and Honary2014; Kulyamin and Dymnikov2016). The coupling between the dayside upper- and lower-ionospheric regions is tenuous as the electron density of the upper ionosphere is largely driven by solar extreme ultraviolet (EUV) radiation, whereas chemistry, mostly involving nitric oxide (NO), controls the quiescent D region (Siskind et al.2017). To facilitate the development of improved D-region models, the impacts of geomagnetic storms must be considered (Spjeldvik and Thorne1975; Dickinson and Bennett1978). Particularly germane to this present discussion is the impact of post-storm energetic particle precipitation (EPP) on the chemistry of the D region (Seppala et al.2015; Rodger et al.2015). Indices used to specify the D region include the reflection height, H (km), and bottom-side profile sharpness, β (km−1), parameters originally developed by James Wait (Wait and Spies1964) and commonly used in research applications (Thomson1993; Thomas et al.2007; Nina et al.2021) as well as operations involving VLF propagation within the EIWG (Nunn et al.2004). Another pair of indices refer to the sunset D-layer disappearance time (DLDT) and the sunrise D-layer preparation time (DLPT) associated with anomalous ionospheric behaviour during geomagnetic storms (Choudhury et al.2015) and earthquakes (Sasmal and Chakrabarti2009; Chakrabarti et al.2010). Within this paper we choose to adopt a modified set of indices originally introduced by Nwankwo et al. (2016), which corresponds to the midday signal amplitude peak (MDP), the mean signal amplitude before sunrise (MBSR) and the mean signal amplitude after sunset (MASS). The use of these indices within the study by Nwankwo et al. (2016) revealed trends, albeit inconsistent, in the signal strength of VLF radio waves propagating within the EIWG in response to several geomagnetic storms of moderate intensity. The intent of the current effort is to expand on the findings of Nwankwo et al. (2016) to further elucidate the effects of geomagnetic storms on the physics of the mid-latitude D region. The following paragraphs review the status of geomagnetic storm-related impacts on the high-, middle- and low-latitude D-region ionosphere (Lastovicka1996).

It was previously noted that solar Lyman-α radiation is primarily responsible for maintaining the quiescent dayside D region (Nicolet and Aikin1960), whereas flare-associated X-rays are the dominant D-region ionisation source during solar flares (Thomson and Clilverd2001; Quan et al.2021). Solar flares are transient events which actively emit ionising X-rays lasting for up to several tens of minutes (Veronig et el.2002). A related class of solar transient is a solar particle event (SPE) resulting from an interplanetary shock (Tsurutani et al.2003; Mittal et al.2011; Chandra et al.2013; Dierckxsens et al.2015; Gopalswamy2018), wherein charged particles, mostly protons, are accelerated to high energies and impact the lower atmosphere at the higher latitudes within the open magnetosphere (Zawedde et al.2018). The precipitation of greater than 10 MeV solar energetic protons (SEPs) both affects the chemistry of the lower mesosphere (Ahrens and Henson 2021), between 50 and 85 km (Turunen et al.2009), and acts as a source of ionisation that can temporarily increase the D-region electron density (Hunsucker1992; Sauer and Wilkinson2008; Neal et al.2013). A sufficiently intense and energetic flux of SEPs can impact radio-wave communications in the form of a polar-cap absorption (PCA) event with a delayed onset following a solar flare, assuming the flare has a related CME, and a duration lasting for up to several days (Rose and Ziauddin1962; Potemra et al.1970; Mitra1974; Rogers and Honary2014; Rogers et al.2016). The Antarctic-Arctic Radiation-belt (Dynamic) Deposition-VLF Atmospheric Research Konsortium (AARDDVARK) network was established to probe the D region with extreme sensitivity (Clilverd et al.2009, 2014; Neal et al.2015). An example of an early use of the AARDDVARK network was monitoring changes in the polar D region from an SPE (Clilverd et al.2007). Strictly speaking, an SPE is quite distinct from a geomagnetic storm, although they both have a common originating source. However, as it relates to our objective of better quantifying the D-region response to geomagnetic storms, we are unaware of any efforts to separate out the storm-specific responses from other, albeit related, sources.

During and following geomagnetic storms, enhanced fluxes of energetic electrons from the outer van Allen radiation belts (van Allen et al.1958) precipitate into the sub-auroral atmosphere (Peter et al.2006), contributing to the formation of the storm time mid-latitude D region (Pedersen1962; Grafe et al.1980; Horne et al.2009; Zawedde et al.2018; George et al.2020). The nominal L-shell location (McIlwain1961) for the outer radiation belt is from L3 to L10 (George et al.2020), which, for an idealised Earth dipole field, maps to invariant latitudes of 55 and 72, respectively (Kilfoyle and Jacka1968). The processes which regulate the electron populations within the outer belts during storm conditions are complicated and interdependent (Reeves and Daglis2016; Baker2019). An aspect not yet discussed is the fundamental role that aurora/magnetospheric substorms play in the dynamics of the magnetosphere (Akasofu1964; McPherron1979; Rostoker et al.1980; Spence1996; Akasofu2020) and, germane to the topic at hand, the impacts on the mid-latitude D region (Guerrero et al.2017). A substorm is described as a “transient process initiated on the night side of the Earth in which a significant amount of energy derived from the solar wind-magnetosphere interaction is deposited in the auroral ionosphere and magnetosphere” (McPherron1979). This description is consistent with the concept of an open magnetosphere wherein the IMF and geomagnetic field lines merge at the dayside magnetopause and are then swept tailward with the solar wind into the nightside, or geotail, where the geomagnetic field lines and IMF, respectively, reconnect (Russell1991). When the IMF has a southward component, the magnetopause standoff distance at noon is nominally located at about 10 Earth radii (RE) (Aubry1970; Fairfield1971; Shue et al.1997; Suvorova and Dmitriev2015; Bonde et al.2018; Samsonov et al.2020). During a geomagnetic storm the magnetopause can be significantly “eroded” (Wiltberger et al.2003; Le et al.2016), reducing the location of the standoff distance which may be within the geostationary orbit of 6.7 RE under extreme conditions (Shue et al.1998). In response to the storm, the polar-cap potential and the cross-tail current increase as energy is continually transferred from the solar wind into the geotail (Angelopoulos et al.2020). While it is curious that the electron density within the outer radiation belt can either increase or decrease under storm conditions (Reeves et al.2020), more relevant to establishing (nightside) or maintaining (dayside) the ionospheric D region is that the electrons within the outer belt can be pumped up to extremely high energies, including relativistic ones, by local wave activity and radial diffusion (Baker2019; Kanekal and Miyoshi2021). A fraction of these energetic electrons can be subsequently scattered into the atmospheric loss cone (Porazik et al.2014) by naturally occurring electromagnetic waves (Spjeldvik and Thorne1975; Gu et al.2020; Ripoll et al.2020; Aryan et al.2021) and precipitate into the lower ionosphere at mid-latitudes, where they collisionally ionise the neutral atmospheric constituents (Rodger et al.2007, 2010, 2012; Naidu et al.2020). The outer radiation belt can relax to its more quiescent state on timescales ranging from minutes (Turner et al.2013) to many days (Baker2019) following a significant geomagnetic storm. It should be noted, in conclusion, that although the evidence provided herein shows a clear association of magnetic storms with enhanced levels of electron precipitation, it is likely that even during geomagnetically quiet intervals electron precipitation persists but apparently at a greatly reduced rate (Mironova et al.2021).

As previously discussed, Appleton and contemporaries (Colwell and Friend1936) in the early 1900s used short-distance VLF radio-wave transmissions to probe and study the D region, whereas the efforts of Wait and colleagues (Wait and Spies1964) were to facilitate global VLF communications within the EIWG. The early work to quantify storm-related D-region impacts involved the general technique of HF radio-wave absorption (Eccles et al.2005; Pederick and Cervera2014; Scotto and Settimi2014; Siskind et al.2017). In this regard, Lauter and Knuth (1967) found that the after-effects of geomagnetic storms on the absorption of 245 kHz radio waves in the mid-latitude D region could persist for more than 10 d. Supporting in situ rocket data (Dickinson and Bennett1978) revealed that the electron density on the days following an “intense” geomagnetic storm could be 4 to 10 times the normal daytime density and that these measurements were well correlated with changes in HF radio-wave absorption. An interesting finding by Satori (1991) was the countering effect of a Forbush decrease in the density of the D region at mid-latitudes. A Forbush decrease refers to the measured reduction in the galactic cosmic ray background due to an Earth-passing CME (Forbush1954; Raghav et al.2020; Janvier et al.2021). Cosmic rays are an important D-region ionisation source at night and a minor contributor during the day (Moler1960). According to Satori (1991), the reduction in the galactic cosmic ray flux associated with a CME counters the increased D-region ionisation from precipitating radiation belt electrons during a geomagnetic storm.

Again, the focus of this report is on the geomagnetic storm-related impacts on VLF TRGCP radio-wave propagation and the information on the D region that can be gleaned from this approach. In this regard, Belrose and Thomas (1968) reported that mid-latitude VLF amplitudes during a geomagnetic storm were unaffected, whereas phase measurements showed rapid fluctuations with residual effects lasting several days following the storm. Muraoka (1979) found that the prevalence of these phase anomalies was dependent on the strength of the magnetic storm. Upon further examination, as reported by Rodger et al. (2007), both signal amplitude and phase variations in VLF TRGCP transmissions were found to be sensitive to electron precipitation events during geomagnetic storms. Choudhury et al. (2015) found that receiver position electron density was the main controlling factor in their storm time metric of the sunrise DLPT depth in VLF TRGCP radio-wave propagation. This metric plus the MDP parameter (after Nwankwo et al.2016) were used by Naidu et al. (2020) to ascertain, again, that geomagnetic-storm-related impacts on the D region were due to energetic electron precipitation. Recently, Kerrache et al. (2021) clarified the role of lightning-induced electron precipitation (LEP) events (Voss et al.1998; Blake et al.2001; Inan et al.2010) in the pitch angle scattering of outer radiation belt electrons during geomagnetic storms and their impact on VLF TRGCP transmissions within the mid-latitude EIWG.

The mechanisms that affect the low-latitude D region in response to geomagnetic storms have not been extensively studied and are still relatively unknown (Araki1974; Kleimenova et al.2004; Kumar and Kumar2014; Kumar et al.2015; Maurya et al.2018). It is recognised that the low-latitude E- and F-region ionospheres can be affected by storm-induced prompt penetration electric fields (PPEFs) (Tsurutani et al.2008; Timocin2022) and disturbance dynamo electric fields (DDEFs) (Blanc and Richmond1980; Fejer et al.1983; Scherliess and Fejer1997) plus related substorm effects (Sastri2006; Chakraborty et al.2015; Hui et al.2017). Generally speaking, PPEFs refer to the immediate and sustained fields resulting from the active response of the magnetosphere to an externally imposed forcing function such as an interplanetary shock (Nava et al.2016), whereas DDEFs result from the delayed equatorward motion of the thermosphere in response to auroral heating and chemistry (Zesta and Oliveira2019; Robinson and Zanetti2021). Transient substorms affect the dynamics of the low-latitude ionosphere as the magnetotail attempts to accommodate an increased storm time reconnection rate (Hajra2021). However, these storm time perturbations do not appear to affect the lower ionosphere except under the most extreme situations. For example, in a limited study of seven moderate (Dst <-50) to intense (Dst <-100) geomagnetic storms, Kumar and Kumar (2014) found that only the intense storm of 16 December 2006, with Dst =-145, had a clear measurable effect on VLF TRGCP transmissions at low latitudes. These findings are consistent with early studies (Araki1974) on the effects of large storms on transequatorial VLF propagation. Other case studies of intense (Kumar et al.2015; Maurya et al.2018) to super (Dst <-200, Wu et al.2016) geomagnetic storms validate this D-region response, although the specific mechanisms for coupling the low-latitude D-region ionosphere to the higher latitudes in terms of gravity waves or chemical processes have not yet been confirmed.

We have provided an extensive and, hopefully, well-documented background on the use of VLF radio waves to probe the ionospheric D region, with an emphasis on the impacts of geomagnetic storms. Monitoring VLF radio-wave propagation within the EIWG along TRGCPs is a convenient and cost-effective technique for determining how space weather affects this lowest traditional ionospheric density layer. While the D region is an important enabling element for VLF communications, it is also an interface linking the middle atmospheric regions to the upper ionosphere. We have stressed the role of the D region within an M–I perspective of a driven system during geomagnetic storms, perhaps at the expense of expanding on the impacts on and the feedback from the neutral atmosphere. While much of the basic physics and chemistry of the D region is well understood, how geomagnetic storms affect this region of space remains an active area of research. Our goal within this paper is to contribute in some small way to the existing body of knowledge concerning the D region as it pertains to VLF radio-wave propagation. Therefore, in this study we combine the observed diurnal VLF amplitude variations in the D region with standard measurements of the E and F regions to perform a diagnostic investigation of coupled geomagnetic storm effects in order to understand the observed storm-induced variations in VLF propagation based on the state and responses of the ionosphere.

2 Data and method

We obtained the VLF amplitude data for the DHO-A118 and GQD-A118 propagation paths received at the A118 SID monitoring station in Muret, southern France (lat 43.53 N, long 1.39 E). The locations of the transmitters (GQD – 22.1 kHz, lat 54.73 N, long 2.88 W – and DHO – 23.4 kHz, lat 53.08 N, long 7.61 W) and the receiver (A118) are shown in Fig. 1. The DHO-A118 and GQD-A118 TRGCPs are 1169 and 1316 km, respectively. Supporting ancillary data include solar X-ray flux, solar wind speed (Vsw) and particle density (PD), planetary geomagnetic Ap and the Dst index. These data adequately described in Nwankwo et al. (2016) and references therein.

Figure 1VLF signal propagation paths (DHO-A118 and GQD-A118) used in the study. Map adopted from the A118 channel; details at (last access: 18 May 2022).

We consider the TRGCP VLF data for the extended intervals 16 through 31 September 2011 and 22 October through 5 November 2011 which included the geomagnetic storm sub-intervals of interest on 17 September, 26–27 September, 25 October and 1 November 2011. The MBSR and MASS indices, as previously described, are the respective 2 to 4 h (dependent on signal quality or season, e.g. winter/summer) mean VLF signal amplitudes measured before local sunrise and after local sunset, whereas the MDP index is the peak signal amplitude measured near local noon. The frequency and power for the DHO and GQD transmitters can be considered constants in our analysis, whereas the received signal strength at A118 is subject to the path-integrated effect of attenuation due to the structure of the D region within the EIWG and the impact of atmospheric absorption. Figure 2 includes several examples of the diurnal variation in VLF signal strength for the DHO-A118 TRGCP. Times are referenced to Universal Time (UT), which is essentially local time of day at the A118 receiver site. The top trace is a rather well-behaved sample where the MBSR, MASS and MDP indices can be readily determined. Also shown in this top trace are the signal effects of the sunrise and sunset terminators, denoted as SRT and SST, respectively, which transit the TRGCP and are generally understood in terms of VLF-mode theory (Wait1960). However, not all diurnal traces are so well behaved. The middle example in Fig. 2 shows the transient effect of a solar flare on the VLF amplitude, which in this case would be classified as a SID and, more specifically, as a sudden enhancement of the signal strength of atmospherics (SEA) in accordance with Lincoln (1961). The lower trace includes an example of a “pseudo-SRT”, which illustrates the complicating effect of terminator-related modal interference along the TRGCP (Chand and Kumar2017). Relatedly, the evening-side SST signatures in the middle and lower traces are rather ambiguous, which, again, may be due to modal interference. The gist of these comments is to highlight the deleterious nature of VLF radio-wave propagation as a complicating factor in deriving the MBSR, MASS, MDP, SRT and SST indices. In previous work (Nwankwo et al.2016, 2020, 2022), we have discussed how these factors can be effectively addressed and the somewhat subjective nature of the determination.

The examples in Fig. 2 clearly illustrate that VLF TRGCP propagation within the EIWG is affected by the dynamics of the ionosphere and, to a lesser extent, the thermosphere. The “classic” quantifications of the VLF reflection height, H and sharpness, b (after Wait and Spies1964), are useful tools for assessing the local D-region response to solar-induced perturbations, for example SIDs, EPP, solar eclipses and geomagnetic storms as well as non-solar phenomena such as gravity waves and lightning (Silber and Price2017). For the TRGCP dataset at hand the received VLF amplitude reflects the integrated state of the traversed ionosphere along the signal path. In this regard we believe that the MBSR, MASS, MDP, SRT and SST indices provide a simpler yet more comprehensive set of metrics for assessing VLF TRGCP radio-wave propagation, particularly in the aftermath of a geomagnetic storm. By comparing and contrasting the changes in these indices for pre-storm versus storm time intervals, we hope to facilitate the development of a more comprehensive D-region specification that can be integrated into an empirical ionospheric model such as the IRI.

Figure 2Diurnal VLF signal amplitude signatures (from the DHO-A118 propagation path) showing analysed signal metrics.


Figure 3(a) Diurnal VLF amplitude for DHO-A118 and (b) GQD-A118 propagation paths, (c) daily variation in X-ray flux output, (d) solar wind speed (Vsw), (e) solar particle density (PD), (f) disturbance storm time (Dst), (g) planetary Ap, and (h) auroral electrojet (AE) indices during 16–30 September 2011.


3 Results and discussion

3.1 Analysis of VLF amplitude variations during intervals of geomagnetic storms

Figure 3 is a composite plot including the diurnal variation in VLF amplitudes for the (a) DHO-A118 and (b) GQD-A118 propagation paths plus the daily variation in the (c) solar X-ray flux output, (d) solar wind speed (Vsw), (e) solar wind PD, (f) Dst, (g) planetary geomagnetic Ap and (h) auroral electrojet (AE) indices during 16–30 September 2011. Four storm conditions occurred during this extended period, including an isolated storm of moderate intensity on 17 September (Dst =-60) and consecutive storms on 26 September (Dst =-101), 27 September (Dst =-88) and 28 September (Dst =-62), presumably driven by the significant increase in Vsw and PD on 17 September and 26 September (Fig. 3a–f). However, the main reference storms are those of 17 September and 26 September. The variation of the AE (especially between 26 September and 29 September) appears to be consistent with storm-related high-intensity, long-duration continuous AE activity events (HILDCAAs) during which “fresh energy was presumably injected” into the magnetosphere (Tsurutani et al.2011). A notable drop occurred in the DHO-A118 VLF signal level on 26 September around midday, coincident with a relatively intense storm with Dst of 101 (Fig. 3a). This scenario (signal strength decrease) has been associated with storm-induced variations in energetic electron precipitation flux (Kikuchi and Evans1983; Peter et al.2006). During a geomagnetic storm, the current system in the ionosphere and the energetic particles precipitating into the ionosphere can influence the density and distribution of density in the atmosphere (NOAA2016). The characterised metrics (e.g. MBSR, MDP, MASS, SST and SRT) of the VLF signal amplitude make it easier to study the behaviour of the signal during the storms. We therefore monitor trends of the signal variation in the analysis to follow for possible identification of storm-induced variation in the signal due to lower ionosphere responses.

Figure 4 shows daily mean fluctuation of Dst and AE and variations in the VLF midday signal amplitude peak (MDP), mean signal amplitude before local sunrise (MBSR), mean signal amplitude after sunset (MASS), sunrise terminator (SRT) and sunset terminator (SST) for the (a) DHO-A118 and (b) GQD-A118 propagation paths from 16 to 30 September 2011. In the GQD-A118 propagation path (Fig. 4a), we observed a dipping of the MDP on 17 September (extending to 20 September) as well as dipping of the MASS on 17 September but an increase in the MBSR, SRT and SST. Following the recurrent storms between 26 and 28 September, we observed dipping of the MDP on 26 September (extending to 29 September). The slight increase in the signal (MDP) on 28 September appears to be due to the significant flare activity (three C-class and one M-class), suggesting an increase in both the instantaneous and background X-ray flux output that usually results in an increase in signal amplitude (as depicted in Fig. 2b). High flare activity can “overshadow” the signal's response to geomagnetic storms when the event coincides with storm time (Nwankwo et al.2016). There is also a significant dipping of all the signal metrics (MDP, MBSR, MASS, SRT and SST) on 27 September. We note dipping of the MBSR on the days following the main (reference) storms on 18 and 27 September. Since the events occurred after dawn (around midday), the post-storm ionospheric effects are expected well into the day following the storm. This trend (post-storm day signal dip) suggests that the signals dipped in response to continuous driving of the ionosphere on the days following the events. However, such a response also depends on the characteristics of the signal propagation path. In the DHO-A118 propagation path, dipping of the MDP, MBSR, SRT and SST occurred on 17 September, while the MDP, MASS and SST also decreased on 26 September. The MASS and SRT maintained the pre-storm day values of 16 and 25 September, respectively. While the MBSR increased slightly on 26 September (main storm day), there is a significant dipping of the signal following the recurrent storm of 27 September.

Figure 4Daily deviations of (a) Dst, (b) AE, (c) variations in the peak value of midday signal amplitude (MDP), (d) mean signal amplitude before local sunrise (MBSR), (e) mean signal amplitude after sunset (MASS), (f) variation in sunrise terminator (SRT), and (g) sunset terminator (SST) for the GQD-A118 (left column) and DHO-A118 (right column) propagation paths from 16 to 30 September 2011.


Figure 5 depicts the diurnal VLF amplitude for the (a) DHO-A118 and (b) GQD-A118 propagation paths, daily variation in (c) X-ray flux output, (d) Vsw, (e) PD, (f) Dst, (g) Ap and (h) AE indices from 22 October to 5 November 2011. This period was associated with a severe storm with the main phase on 25 October (Dst =-132) and consecutive storms on 1 November (Dst =-71) and 2 November (Dst =-57), presumably driven by the highly variable Vsw and PD (Fig. 5d–e). It has been shown that the capability of a given value of the solar wind electric field (SWEF) to create a Dst disturbance or geo-efficiency is enhanced by high solar wind density (Weigel2010; Tsurutani et al.2011). Variation of the AE between 30 October and 3 November also appears to be consistent with a HILDCAA event (Fig. 5h). The DHO-A118 VLF signal level on 25 October around midday also showed a visible reduction following the intense storm condition with a Dst of 132 (Fig. 5a). VLF signal data for the GQD-A118 propagation path are not available from 12:00 on 25 October to 18:00 on 26 October (Fig. 5b).

Figure 5(a) Diurnal VLF amplitude for the DHO-A118 and (b) GQD-A118 propagation paths, (c) daily variation in X-ray flux output, (d) solar wind speed (Vsw), (e) solar particle density (PD), (f) Dst, (g) planetary Ap, and (h) AE indices from 22 October to 5 November 2011.


Figure 6 plots the daily deviations of Dst and AE and variations in the MDP, MBSR, MASS, SRT and SST for the (a) DHO-A118 and (b) GQD-A118 propagation paths from 22 October to 5 November 2011. Although data for the GQD-A118 propagation path during 25 and 26 October were inadequate for the present analysis, we did observe a dipping of the MBSR on the main storm day (25 October). Dipping of the MDP, MASS and SST occurred on 1 November and those of MBSR, MASS and SRT on 2 November, following the consecutive storms. In the DHO-A118 propagation path, we observed dipping of the MDP, MBSR, MASS and SRT on 25 October, dipping of the MDP, MBSR, MASS and SST on 1 November and dipping of the MBSR and SRT on 2 November. Similar to the first case (Figs. 4 and 5), we note the numerous flare events on 2 November (up to seven C-class and one M-class) that may have induced a spike in the MDP on the day in both the GQD-A118 and DHO-A118 propagation paths. Although dipping of the MDP signal (following storm events) has shown considerable consistency across the cases presented so far, the MBSR and MASS (in particular) appear to be influenced by storm occurrence time and the high variability or fluctuation of the dusk-to-dawn ionosphere (and signal) (Nwankwo et al.2016). However, presenting a consistency across a substantial number of cases is vital to the conclusion of this work. We, therefore, analysed 15 additional storm cases between September 2011 and October 2012 in order to obtain a statistically significant set of observations. The 15 storm cases are presented in Table 1.

Figure 6Daily deviations of (a) Dst and (b) AE and (c) variations in MDP, (d) MBSR, (e) MASS, (f) SRT, and (g) SST for the GQD-A118 (left column) and DHO-A118 (right column) propagation paths from 22 October to 5 November 2011.


Table 1 Summary of the 15 analysed geomagnetic storm events.

Download Print Version | Download XLSX

In Fig. 7, we show Dst deviations (σDst) and trends in variation of the MDP, MBSR, MASS, SRT and SST signals on the day before the storm (blue bar), the storm day (red bar) and after the storm day (brown bar) for the 15 selected storm cases in the (a) GQD-A118 and (b) DHO-A118 propagation paths. σDst is the measure or extent of daily fluctuation in measured values of Dst. We recognised the 3 consecutive days as day before event (BE), during event (DE) and after event (AEv). A missing index indicates an absence of data. It should be noted however (for this analysis) that these events were separate events and not continuous events. In the GQD-A118 propagation path, about 8 of 12 MDPs, 10 of 13 MBSRs, 7 of 12 MASSs, 3 of 12 SRTs and 5 of 12 SSTs showed dipping features, while 12 of 15 MDPs, 9 of 15 MBSRs, 10 of 15 MASSs, 5 of 15 SRTs and 7 of 15 SSTs showed dipping in the DHO-A118 propagation path. These values correspond to 67 %, 77 %, 58 %, 25 % and 42 % dipping in the GQD-A118 propagation path and 80 %, 60 %, 67 %, 33 % and 47 % dipping in the DHO-A118 propagation path. The signal levels along with the percentage dip are presented in Table 2. The MDP signals (in both propagation paths) generally show a tendency for dipping following geomagnetic storm conditions. However, we also observe a few cases of increase in the MDP during separate events in each of the propagation paths (e.g. events 4 and 7 in GQD-A118 and 9 in DHO-A118) as well as coincident increases occurring in both propagation paths during the same event (e.g. events 3 and 12). While the probable reason for this coincidence is suggestive of factors such as propagation characteristics and/or X-ray flux-induced spikes in amplitude (e.g. significant X-ray output during events 3 and 7 in Fig. 8c), further investigation into why this characteristic exists will be pursued. In the meantime, we analysed variations in X-ray flux output and geomagnetic indices during events 3 and 12 to better interpret the prevailing ionospheric conditions at the time.

Figure 7Dst deviation (or fluctuation) and variations in MDP, MBSR, MASS, SRT and SST signals 1 d before, during and after each of the 15 events for the GQD-A118 and DHO-A118 propagation paths. Note that each Dst bar represents the deviation (σ) corresponding to the VLF amplitude before, during and after the selected events (storm) as listed in Table 1 and shown in Fig. 8. The events are selected (for the intervals shown in Fig. 8) and are not continuous.


Table 2 Summary of the trend in dipping of the signals' metrics during 15 geomagnetic storm cases in the (a) DHO-A118 and GQD-A118 propagation paths.

Download Print Version | Download XLSX

In Fig. 8, we present the diurnal VLF amplitudes for the (a) DHO-A118 and (b) GQD-A118 propagation paths, daily variation in (c) X-ray flux output, (d) Vsw, (e) PD and (f) Dst indices for a day before and after each of the 15 storms. Data showed (Fig. 8c, f) the occurrence of an M-class flare in association with the storm on 22–23 January 2012 (event 3 on 21 January), both events having near-corresponding peaks. This scenario suggests an enhancement of both the instantaneous and background X-ray flux output (as stated earlier) that may have caused an increase (or spike) in the signal level and thus probably overshadowed geomagnetic effects on the signal. While this explanation may be argued for events 1 (25–27 September 2011) and 6 (6–8 March 2012), it should be noted that the flare events started well before the storms and continued until the storm time (in each case), suggesting an established increase in the overall background X-ray before the storms. Hence, it is possible for a storm-induced dipping to manifest under such a condition. However, further investigation is encouraged, which is beyond the scope of this work. For event 12 (from 15 to 17 July 2012), we observed that the peak of the storm (that commenced by midnight on 16 July) was on 17 July (recognised as AEv). Therefore, any geomagnetic influence on the signal (e.g. dipping) is expected on 17 July (or after) and not 16 July, hence the dipping of the AEv signal (on 17 September) in the DHO-A118 propagation path.

Figure 8Diurnal VLF amplitude for the (a) DHO-A118 and (b) GQD-A118 propagation paths, daily variation in (c) X-ray flux output, (d) Vsw, (e) PD, and (f) Dst indices for a day before and after each of the 15 storms.


Figure 9 shows Dst deviation (fluctuation) and 2 d mean variations of MDP, MBSR, MASS, SRT and SST signals before, during and after each event for the (a) GQD-A118 and (b) DHO-A118 propagation paths. This analysis is important for corroborating the result presented in Fig. 7, because the data selection criteria differ from those of Fig. 7 in some ways. While BE, DE and AEv represent data for 3 consecutive days with reference to the event's day (DE) in the former analysis (presented in Fig. 7), each abbreviation (BE, DE or AEv) represents a relatively quiescent 2 d mean amplitude before and after DE (but not necessarily in succession to/after DE). However, it should be noted that, due to the data averaging (2 d), a “pronounced” increase or dipping in the signals (comparable to those in the former analysis, Fig. 7) is not expected. Another important data selection criterion for this analysis is a relatively geomagnetic quiet day BE and AEv with respect to DE.

Figure 9Dst deviation and 2 d mean variations of MDP, MBSR, MASS, SRT and SST signals before, during and after each event for the GQD-A118 and DHO-A118 propagation paths.


In the GQD-A118 propagation path, 7 of 12 MDPs, 7 of 13 MBSRs, 7 of 12 MASSs, 6 of 12 SRTs and 3 of 12 SSTs showed dipping following the storms, while 10 of 15 MDPs, 11 of 15 MBSRs, 11 of 15 MASSs, 6 of 14 SRTs and 6 of 15 SSTs showed dipping in the DHO-A118 propagation path. These values correspond to, respectively, 58 %, 54 %, 58 %, 50 % and 25 % dipping in the GQD-A118 propagation path and 67 %, 73 %, 73 %, 43 % and 40 % dipping in the DHO-A118 propagation path. The signal levels along with the percentage dip of the signals are presented in Table 3. In general, the trend of variation of the signal metrics reflected the prevailing space weather-coupled effects in the lower ionosphere. The MDP signal appears to be more responsive to geomagnetic perturbations (58 % and 67 % (67 % and 80 % for 1 d mean analysis) in the respective propagation paths shown in Figs. 7 and 9) than other signal metrics. However, the 2 d mean analysis showed improvement in MBSR and MASS (73 %) of the DHO-A118 propagation path, thereby reenforcing the responsiveness of this propagation path to geomagnetic storm impacts. While the mechanism of VLF amplitude (and/or MDP) response to flare-induced SIDs is well understood (see Sect. 1.3 of Nwankwo et al.2016), the mechanism of the transient response of the signal to geomagnetic storming is not well developed. We speculate that the dipping response of the MDP may be related to (i) a positive storm effect, which affects, albeit to a small extent, the attenuation of the VLF radio waves (Fagundes et al.2016), (ii) adjustment of the D layer to storm-driven energy input (McPherron1979), (iii) precipitation of energetic electrons (Rodger et al.2010, 2012; Naidu et al.2020) and/or (iv) charge exchange between surrounding ionospheric regions (Tatsuta et al.2015). Since the processes are gradual (when compared with the SID scenario) and originate from the magnetosphere, it is not likely that a remote spike would occur in the diurnal signal (as observed during the flare condition). One other likely reason for the observed dipping characteristic is the modification of the chemistry of the D region by storm precipitation of SEPs (Turunen et al.2009), since chemistry (mainly NO) controls the quiescent D region (Siskind et al.2017).

Nwankwo et al. (2016) noted the existence of pseudo SRT and/or SST sometimes exhibited by diurnal VLF signal (see Fig. 2c) as a drawback in SRT and SST analysis. This anomaly is due to the secondary destructive interference pattern in signals and/or occurrence of solar flares during sunrise/sunset (Sandip Chakrabarti, personal communication, 2016). The authors concluded in their study that the post-storm SRT and SST variations do not appear to have a well-defined trend associated with the storm effect based on the approach used in the analysis (Nwankwo et al.2016). In this work, we considered the “first” SRT and SST values (in the event of a pseudo-terminator) during analysis of the signal metrics. A rise in SRT and SST amplitude under geomagnetic storm conditions appears to occur more than otherwise in both propagation paths. We found a respective dipping of 50 % and 25 % (43 % and 40 %) of the SRT and SST. Storm-induced disturbances may not have a significant influence on these metrics (SRT and SST), and since the sunrise and sunset signatures relate to mode conversion in the VLF propagation path, it might imply that the D-region density is not a significant contributor to this effect. It is important to note that, of the two propagation paths used in this study, the DHO-A118 signal appears to be more sensitive to geomagnetic storm-induced magnetosphere–ionospheric dynamics. We do not expect a “perfect” consistency in signal trend and variations across all cases, because the individual effects of solar and other forcing mechanisms (including those of lithospheric and atmospheric sources) on the ionosphere are difficult to estimate (Kutiev2013; Nwankwo et al.2016). This scenario can also cause non-linear coupling processes and consequent significant fluctuations in radio signals.

Table 3 Summary of trends in 2 d mean signal dipping following 15 geomagnetic storm cases in the DHO-A118 and GQD-A118 propagation paths.

Download Print Version | Download XLSX

3.2 Investigating the state of the ionosphere over the propagation paths of the VLF signals

Here, we study the state of the ionosphere over the two VLF propagation paths using the virtual heights (h'E, h'F1 and h'F2) and critical frequencies (foE, foF1 and foF2) of the E and F regions obtained from two ionosonde stations near the GQD and DHO transmitters. Although we made an effort to obtain data from stations near each transmitter/receiver and at the mid-point, we found no ionosonde station at the mid-point, and the nearest station to the receiver (Tortosa) has no data for the period/intervals under study. However, to make up for this dearth of data, we will complement the analysis with the results in the extended study that utilised the GNSS data in the region (e.g. Nwankwo et al.2022). Details of the ionosonde stations used in this study are provided in Table 4. We treat Chilton station as nearest to the GQD transmitter and Juliusruh station as nearest to the DHO transmitter. Tortosa station is closest to the A118 transmitter but has no data for the intervals. We obtained and calculated the daytime (08:00–15:00) mean values and standard deviations (σ) of the parameters and then analysed them for the storms of interest (on 17 and 26 September and 1 November) within the intervals 16–19 and 25–28 September and 29 October to 2 November 2011. We exclude analysis of the 25 October storm because the data for this interval are inadequate.

Table 4 Ionosonde stations near the VLF transmitters and receiver and/or propagation paths.

Download Print Version | Download XLSX

Figure 10 shows the daily mean and standard deviation (SD or σ) of foF2, foF1, foEs, foE, h'F2, h'F, h'Es and h'E from 16 to 17 September 2011 for Chilton and Juliusruh stations. We compared the pre-storm day (blue broken line) values with the storm day (red broken line) values. At Chilton station (near the GQD transmitter), results show significant increase and/or fluctuation (increase in SD) in foF2 and a decrease (with significant fluctuation) in foF1 on the storm day, 17 September. The height of the E and F regions (h'F2, h'F, h'Es and h'E) significantly increased following the storm. A similar pattern of variations was observed at Juliusruh station. foF2 and foF1 increased (and/or fluctuated) significantly, as well as h'F2, h'F, h'Es and h'E. Values of foEs decreased, while foE remained unaffected in both stations. Also, there appear to be a sustained post-storm increase and/or fluctuations of the parameters on 18 September, suggesting a continuous driving of the ionosphere by the storm.

Figure 10Daytime mean variations and standard deviation (SD) of foF2, foF1, foEs, foE, h'F2, h'F, h'Es and h'E from 16 to 17 September 2011 for Chilton and Juliusruh stations. The blue broken line represents the pre-storm day values, while the red broken line represents the storm day values of the ionospheric parameters.


Figure 11 shows the daytime mean variations and SD of foF2, foF1, foEs, foE, h'F2, h'F, h'Es and h'E from 25 to 28 September 2011 for Chilton and Juliusruh stations. The storm during this interval (on 26 September) was well developed (with Dst up to 101 nT) and larger than the 17 September event. The result of this analysis shows a slight increase in foF2 and foF1 but a decrease in foEs and foE for Chilton station. The height of F2 (h'F2) decreased (by 6.90 km), while those of F, Es and E increased on the storm day, 26 September. Near the DHO transmitter (Juliusruh station), there is an anti-correlated variation in the critical frequencies of the E and F regions: a depression of foF2 and foF1 but an increase in foEs and foE (when compared with the scenario at Chilton station). The height of the F2, F and E regions increased by 47.89, 16.08 and 9.14 km, respectively (which are so far the largest increase in the parameters), while the height of the Es region decreased by 0.16 km (see Table 5).

Figure 11Daytime mean variation and SD of foF2, foF1, foEs, foE, h'F2, h'F, h'Es and h'E from 25 to 28 September 2011 for Chilton and Juliusruh stations.


Figure 12 shows the daytime mean variation and SD of foF2, foF1, foEs, foE, h'F2, h'F, h'Es and h'E from 29 October to 2 November 2011 for Chilton and Juliusruh stations. This interval is of interest because of the fluctuation in geophysical parameters during the days preceding the storm. We include this interval to investigate the coupling effect of this extended period (30–31 October) of geomagnetic disturbances preceding the storm on 1 November. It appears that energy began building up in the magnetosphere–ionosphere system after the first significant spike in Vsw around 07:00 on 29 October (and subsequent increase on 30 October) and PD around 10:00 on 30 October until around 10:00 on 1 November, when the storm was triggered following a sudden increase in Vsw and southward turning of the Bz (Nwankwo et al.2020). Here, we compare the parameter level on the relatively quiet day (29 October) with those of the storm day (on 1 November), since the 2 days preceding the storm were significantly disturbed. The result shows significant increases in foF2 and foF1 at Chilton station. Like the 17 September storm scenario, values of foEs decreased, while foE remained unaffected for this station. h'F2 decreased, while h'F, h'Es and h'E showed significant increase. At Juliusruh station only the critical frequency of the F2 region increases, while those of F1, Es and E decreased. However, the increase and/or fluctuation of the parameters were significant (in most cases) during the disturbed days (30 and 31 October) preceding the storm, suggesting responses of the E and F ionosphere regions (coupled to the D region) before the storm commencement (as a result of increased geomagnetic activity on the days). We present a summary of the storm day variations in h'F2, h'F, h'Es and h'E for the two stations in Table 5.

Figure 12Daytime mean variation and SD of foF2, foF1, foEs, foE, h'F2, h'F, h'Es and h'E from 29 October to 2 November 2011 for Chilton and Juliusruh stations.


Figure 13Daytime variation in VLF amplitude (red line plot) for the DHO-A118 propagation path, together with VTEC values obtained from HERT (black line), EUSK (blue line), OPMT (green line) and ESCO (brown line) stations across some locations in Europe during 16–19 and 25–28 September, 24–27 October and 29 October–1 November 2011 (from Nwankwo et al.2022).


Table 5 Observed increase (or decrease) in h'F2, h'F, h'Es and h'E during the storms on 17 and 25 September and 1 November 2011.

Download Print Version | Download XLSX

In summary, foF2, foF1, h'F2, h'F, h'Es and h'E generally showed significant increases and/or fluctuations near both transmitters (GQD and DHO) during the geomagnetic storms, whereas foEs and foE mostly decreased (but slightly increased in fewer cases) or remained unaffected. It appears that the observed storm-induced increases and fluctuations were largely sustained or further enhanced on the day (or days) following the event (post storm day), suggesting a continuous driving of the ionosphere by the storms and/or a substorm effect. Although the analysis for the 1 November storm scenario showed weak correlation, variations of the parameters reflected the coupled responses of the ionosphere to energy build-up ahead of storm commencement. Nwankwo and Chakrabarti (2018) reported significant depression and fluctuations of foF2 following significant geomagnetic disturbances and/or storms at high and mid latitudes and distortion in the quasi-periodic pattern of the parameter. Their inference was, however, based on a preliminary analysis from the result of a single ionosonde station. From this comparatively detailed analysis, it is clear that the reported depression of foF2 may occur during some (isolated) storms and locations, and should, therefore, not be treated as a global response. Negative storm effects (in which foF2 decreases) have also been reported (e.g. Blanch et al.2013; Kane2005). In this analysis, the largest increase in h'F2, h'F, h'Es and h'E occurred in Juliusruh near the DHO transmitter (see Table 5). The ionosonde observations indicated that fluctuations in the reference heights appeared to be the dominant response of the E and F regions to geomagnetic storms, whereas attenuation of the VLF radio-wave signal strength was responsive to the storm-induced dynamics in the ionospheric D region. This observation is instructive in that the observed large ionosonde increase and the large amplitude decreases in the DHO-A118 propagation path signal may be related to coupled effects between the ionospheric regions but are also suggestive of strong storm responses (more intense) around/near the DHO transmitter or DHO-A118 propagation path. This result is in agreement with the recent findings reported in Nwankwo et al. (2022). Their study combined observed VLF amplitude variations with TEC/VTEC data obtained from multiple GNSS stations, including Euskirchen in Germany (EUSK), Hailsham in the UK (HERT), Paris in France (OPMT) and Naut Aran in Spain (ESCO), to investigate the ionospheric response to storms over some signal propagation paths during the same events. They showed and reported a simultaneous decrease in VLF amplitude and enhancement of electron density profiles near the DHO transmitter. In Fig. 13 we show the daytime variation in VLF amplitude (red line plot) for the DHO-A118 propagation path, together with VTEC values obtained from HERT (black line), EUSK (blue line), OPMT (green line) and ESCO (brown line) stations across some locations in Europe during 16–19 and 25–28 September, 24–27 October and 29 October–1 November 2011. HERT is closest to the GQD transmitter (about 508.12 km), EUSK is closest to the DHO transmitter (about 279.99 km), while ESCO is nearest to the receiver (about 90.47 km).

One clear observation (from Fig. 13) is the strong dipping (or reduction) of the daytime VLF amplitude and the simultaneous increase in VTEC values on the storm days in the DHO-A118 propagation path. We note the large increase in VTEC values for ESCO station located near DHO, especially during the 17 September and 25 October storms. This feature is in agreement with the findings of Choudhury et al. (2015), who reported that the receiver position electron density is the main factor influencing the VLF signal at ionospheric sunrise time during long-duration geomagnetic storms. It is also worth noting that the ancillary information of the timing, classification and location of associated solar flares, CMEs, SPEs, and timings for the SSCs showed that the strong storm intervals during which large dipping or decrease in DHO signal level occurred were associated with SPEs (see Table 5 in Nwankwo et al.2022). The results of this effort that combined the diagnostics of the D, E and F regions (to probe the geostorm effects in the lower ionosphere) demonstrate that, despite the tenuousness of the coupling between the dayside upper- and lower-ionospheric regions, the adjoining regions of E and F play significant roles in driving the storm-induced dynamics of the D region and the associated observed responses of VLF radio waves within the context of solar–terrestrial coupling.

4 Conclusions

In this work we performed a diagnostic study of geomagnetic storm-induced disturbances that were coupled to the lower ionosphere in the mid-latitude D region using propagation characteristics of VLF radio signals. We characterised the diurnal signal into five metrics (i.e. MBSR, MDP, MASS, SRT and SST) and monitored the trend in variations of the signal metrics for up to 20 storms between September 2011 and October 2012. The goal of the analysis was to understand deviations in the signal that were attributable to the storms. Up to five (5) storms and their effects on the signals were studied in detail, followed by statistical analysis of 15 other cases. Our results showed that the MDP exhibited characteristic dipping in about 67 % and 80 % of the cases in the GQD-A118 and DHO-A118 propagation paths, respectively. The MBSR showed respective dipping of about 77 % and 60 %, while the MASS dipped by 58 % and 67 %. Conversely, the SRT and SST showed respective dipping of 25 % and 33 % and 42 % and 47 %, favouring the rise of the signals following storms. The MDP consistently showed strong responses to the storms than the other metrics (followed by the MBSR and the MASS). Among other possible reasons outlined in this paper, we speculate that the responses were related to positive storm effects, resulting in an attenuation of the VLF radio waves. Of the two propagation paths examined in this study, we observed stronger dipping of the VLF amplitude of the DHO-A118 propagation path during the storms. To understand the state of the ionosphere over the propagation paths and examine how the upper ionosphere (E and F regions) might have affected VLF transmissions within the EIWG, we further analysed virtual heights (h'E, h'F1 and h'F2) and critical frequencies (foE, foF1 and foF2) of the E and F regions (from ionosonde stations near the GQD and DHO transmitters). The results of this analysis showed a significant increase and/or fluctuation in the height of the E and F regions (h'F2, h'F, h'Es and h'E) near both transmitters during the geomagnetic storms, with the largest increase occurring in Juliusruh (Germany) station near the DHO transmitter. This scenario suggest a strong storm response over the region, possibly leading to the large dipping of VLF amplitude for the DHO-A118 propagation path. The ionosonde observations show that fluctuations in the reference heights appear to be the dominant responses of the E and F regions to geomagnetic storms, whereas dipping of the VLF radio waves reflects storm-induced dynamics in the ionospheric D region. Our findings demonstrate that ionospheric E and F regions play significant roles in driving the storm-induced dynamics of the D region and the associated observed responses of VLF radio waves despite the tenuousness of the coupling between the dayside upper- and lower-ionospheric regions.

Data availability

The VLF data which formed the basis of this report were acquired from the A118 SID monitoring station website at (A1182022). Supporting ionosonde data were obtained from the UK Solar System Data Centre at (UKSSDC2022). The geomagnetic indices for Ap, AE and Dst were sourced from the German Research Centre for Geosciences at (GFZ2022) and the World Data Center for Geomagnetism at (WDCG2022). Solar-interplanetary data for the X-ray flux, Vsw and PD were secured from the National Aeronautics and Space Administration website at (NASA2022).

Author contributions

VUJN conceived the idea of the study, designed the methodology, data processing and analysis and coordinated the interpretation and discussion of the results and writing the paper. WD led the writing and editing of the paper and interpretation of solar-geophysical data and phenomena. SKC assisted in conceiving the study and contributed to the discussion of results. MPA contributed in data analysis work and discussion of results. OO, JOF, PIA, OEO, OEO and FVF contributed in discussion and interpretation of results and validation of paper content.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Victor U. J. Nwankwo thanks the World Academy of Sciences (Italy) and the S. N. Bose National Centre for Basic Sciences (India) for the award of a postgraduate research fellowship during which a portion of this work was performed. The authors are grateful for the contributions of Lionel Loudet (A118) for maintaining the A118 SID station and providing access to the VLF data. We appreciate the assistance of Matthew Wild (UKSSDC) for acquiring the ionosonde data. The solar-geophysical data were kindly provided by GFZ, NASA, and WDCG. Victor U. J. Nwankwo thanks Bruce Tsurutani and Jan Lastovicka for useful discussions and clarification of certain aspects related to the geomagnetic indices and geophysical processes.

Review statement

This paper was edited by Keisuke Hosokawa and reviewed by Sergey Sokolov and one anonymous referee.


A118: SID Monitoring Station,, last access: 15 May 2022. a

Abd Rashid, M. M., Ismail, M., Hasbie, A. M., Salut, M. M., and Abdullah, M.: VLF observation of D-region disturbances associated with solar flares at UKM Selangor Malaysia, IEEE International Conference on Space Science and Communication (IconSpace), 249–252,, 2013. a

Ahrens, C. D. and Henson, R.: Meteorology Today: An Introduction to Weather, Climate and the Environment, Cengage Learning, Boston, MA, 736, ISBN-13: 978-0357452073, 2021. a

Akasofu, S.: Relationship Between Geomagnetic Storms and Auroral/Magnetospheric Substorms: Early Studies, Frontiers Astron. Space Sci, 7, 16,, 2020. a

Akasofu, S.-I.: The development of the auroral substorm, Planet. Space Sci., 12, 273–282,, 1964. a

Akasofu, S.-I.: A Review of the Current Understanding in the Study of Geomagnetic Storms, Int. J. Earth Sci. Geophys., 4, 13,, 2018. a

Andersson, M. E., Verronen, P. T., Marsh, D. R., Paivarinta, S.-M., and Plane, J. M. C.: WACCM-D-Improved modeling of nitric acid and active chlorine during energetic particle precipitation, J. Geophys. Res.-Atmos., 121, 10328–10341,, 2016. a

Angelopoulos, V., Artemyev, A., Phan, T. D., and Miyashita, Y.: Near-Earth Magnetotail Reconnection Powers Space Storms, Nat. Phys., 16, 317–321,, 2020. a

Appleton, E. V.: The Existence of more than one Ionised Layer in the Upper Atmosphere, Nature, 120, 1476–4687,, 1927. a

Appleton, E. V.: Meeting for discussion on the ionosphere, P. Roy. Soc. Lond. A, 141, 697–721,, 1933. a

Appleton, E. and Barnett, M.: Local Reflection of Wireless Waves from the Upper Atmosphere, Nature, 115, 333–334,, 1925a. a

Appleton, E. V. and Barnett, M. A. F.: On some direct evidence for downward atmospheric reflection of electric rays, P. Roy. Soc. Lond. A., 109, 621–641,, 1925b. a, b

Appleton, E. V. and Barnett, M. A. F.: On wireless interference phenomena between ground waves and waves deviated by the upper atmosphere, P. Roy. Soc. Lond. A, 113, 450–458,, 1926. a

Appleton, E. V. and Naismith, R.: Some further measurements of upper atmospheric ionization, P. Roy. Soc. Lond. A, 150, 685–708,, 1935. a

Araki, T.: Anomalous Phase Changes of Trans equatorial VLF Radio Waves during Geomagnetic Storms, J. Geophys. Res., 79, 4811–4813,, 1974. a, b

Aryan, H., Bortnik, J., Meredith, N. P., Horne, R. B., Sibeck, D. G., and Balikhin, M. A.: Multi-parameter chorus and plasmaspheric hiss wave models, J. Geophys. Res.-Space, 126, e2020JA028403,, 2021. a

Aubry, M. P., Russell, C. T., and Kivelson, M. G.: Inward motion of the magnetopause before a substorm, J. Geophys. Res., 75, 7018–7031,, 1970. a

Baker, D. N.: Effects of the Sun on the Earth's environment, J. Atmos. Sol.-Terr. Phys., 62, 1669–1681,, 2000. a

Baker, D. N., Hoxie, V., Zhao, H., Jaynes, A. N., Kanekal, S., Li, X., and Elkington, S.: Multiyear measurements of radiation belt electrons: Acceleration, transport, and loss, J. Geophys. Res.-Space, 124, 588–2602,, 2019. a, b, c

Banks, P. M. and Kockarts, G.: Aeronomy, Academic Press Inc, NY, USA, eBook ISBN 978-1-48326-0-068, 1973. a

Barr, R., Jones, D. L., and Rodger, C. J.: ELF and VLF radio waves, J. Atmos. Sol.-Terr. Phys., 62, 1689–1718,, 2000. a

Bates, D. R. and Massey, H. S. W.: The basic reactions in the upper atmosphere, P. Roy. Soc. Lond. A., 187, 261–296,, 1946. a

Bates, D. R. and Massey, H. S. W.: The basic reactions in the upper atmosphere II. The theory of recombination in the ionized layers, P. Roy. Soc. Lond. A, 192, 1028,, 1947. a

Belehaki, A., James, S., Hapgood, M., Ventouras, S., Galkin, I., Lembesis, A., Tsagouri, I., Charisi, A., Spogli, L., Berdermann, J., and Häggström, I.: The ESPAS e-infrastructure: Access to data from near-Earth space, Adv. Space Res., 58, 1177–1200,, 2016. a, b

Belrose, J. S. and Thomas, L.: Ionization changes in the middle latitude D-region associated with geomagnetic storms, J. Atmos. Sol.-Terr. Phys., 30, 1397–1413,, 1968. a

Bennington, T. W.: Radio Waves and the Ionosphere, Nature, 154, 413,, 1944. a

Benson, R. F.: Four Decades of Space-Borne Radio Sounding, URSI Radio Science Bulletin, 2010, 24–44, 2010. a

Betz, H. D., Schmidt, K., and Oettinger, W. P.: LINET–An international VLF/LF lightning detection network in Europe, in: Lightning: principles, instruments and applications, Springer, Dordrecht, 115–140,, 2009. a

Beynon, W. J. G.: The physics of the ionosphere, Sci. Prog., 57, 415–433, (last access: 18 May 2022), 1969. a

Bibl, K.: Evolution of the ionosonde, Ann. Geophys., 41,, 1998. a

Bilitza, D.: Electron density in the D-region as given by the International Reference Ionosphere, in: International Reference Ionosphere IRI 79, UAG-82, World Data Center A for Solar-Terrestrial Physics, edited by: Lincoln, J. V. and Conkright, R.O., 7–10, (last access: 2 Decmeber 2021), 1981. a

Bilitza, D.: International Reference Ionosphere 1990, NSSDC, Report 90-22, Greenbelt, MD, 160, (last access: 2 Demceber 2021), 1990. a

Bilitza, D.: The E- and D-region in IRI, Adv. Space Res., 21, 871–874,, 1998. a, b

Bilitza, D.: International Reference Ionosphere 2000, Radio Sci., 36, 261–275,, 2001. a

Bilitza, D.: IRI the International Standard for the Ionosphere, Adv. Radio Sci., 16, 1–11,, 2018. a

Bilitza, D. and Reinisch, B. W.: International Reference Ionosphere 2007: Improvements and new parameters, Adv. Space Res., 42, 599–609,, 2008. a

Blake, J. B., Inan, U. S., Walt, W., Bell, T. F., Bortnik, J., Chenette, D. L., and Christian, H. J.: Lightning-induced energetic electron flux enhancements in the drift loss cone, J. Geophys. Res., 106, 29733–29744,, 2001. a

Blanc, M.: Magnetosphere-Ionosphere Coupling, Comput. Phys. Commun., 49, 103–118,, 1988. a

Blanc, M. and Richmond, A.: The ionospheric disturbance dynamo, J. Geophys. Res., 85, 1669–1686,, 1980. a

Blanch, E., Marsal, S., Segarra, A., Torta, J. M., Altadill, D., and Curto, J. J.: Space weather effects on Earth's environment associated to the 24–25 October 2011 geomagnetic storm, Space Weather, 11, 153–168,, 2013. a

Bonde, R. E. F., Lopez, R. E., and Wang, J. Y.: The effect of IMF fluctuations on the subsolar magnetopause position: A study using a global MHD model, J. Geophys. Res.-Space, 123 2598–2604,, 2018. a

Borovsky, J. E. and Denton, M. H.: Differences between CME-driven storms and CIR-driven storms, J Geophys. Res., 111, A07S08,, 2006. a

Borovsky, J. E. and Shprits, Y. Y.: Is the Dst index sufficient to define all geospace storms? J. Geophys. Res.-Space, 122, 11543–11547,, 2017. a

Breit, G., and Tuve, M. A.: A radio method of estimating the height of the conducting layer, Nature, 116, p. 357,, 1925. a

Budden, K. D.: I. The Propagation of a Radio-Atmospheric, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 42, 1–19,, 1951. a

Budden, K. G.: The propagation of very low frequency radio waves to great distances, Philos. Mag., 44, 504–513,, 1953. a

Budden, K. G.: The “Waveguide Mode” Theory of the Propagation of Very-Low-Frequency Radio Waves, Proc. IRE, 45, 772–774,, 1957. a

Burch, J. L.: Magnetosphere-Ionosphere Coupling, Past to Future, in: Magnetosphere‐Ionosphere Coupling in the Solar System, Geophysical Monograph Series, edited by: Chappell, C. R., Schunk, R. W., Banks, P. M., Burch, J. L., and Thorne, R. M., ISBN: 9781119066774, 1–17,, 2016. a

Cassak, P. A.: Inside the Black Box: Magnetic Reconnection and the Magnetospheric Multiscale Mission, Space Weather, 14, 186–197,, 2016. a

Chakrabarti, S. K., Sasmal, S., and Chakrabarti, S.: Ionospheric anomaly due to seismic activities – Part 2: Evidence from D-layer preparation and disappearance times, Nat. Hazards Earth Syst. Sci., 10, 1751–1757,, 2010. a

Chakraborty, M., Kumar, S., De, B. K. and Guha, A.: Effects of geomagnetic storm on low latitude ionospheric total electron content: A case study from Indian sector, J. Earth Syst. Sci., 124, 1115–1126,, 2015. a

Chand, A. E. and Kumar, S.: VLF modal interference distance and nighttime D region VLF reflection height for west-east and east-west propagation paths to Fiji, Radio Sci., 52, 1004–1015,, 2017. a

Chandra, R., Gopalswamy, N., Mäkelä, P., Xie, H., Yashiro, S., Akiyama, S., Uddin, W., Srivastava, A. K., Joshi, N. C., Jain, R., Awasthi, A. K., Manoharan, P. K., Mahalakshmi, K., Dwivedi, V. C., Choudhary, D. P., and Nitta, N. V.: Solar energetic particle events during the rise phases of solar cycles 23 and 24, Adv. Space Res., 52, 2102–2111,, 2013. a

Chapman, S.: The absorption and dissociative or ionizing effect of monochromatic radiation in an atmosphere on a rotating earth, Proc. Phys. Soc., 43, 26–45,, 1931. a

Chapman, S. and Ferraro, V.: A New Theory of Magnetic Storms, Nature, 126, 129–130,, 1930. a

Chapman, J. H. and Warren, E. S.: Topside sounding of the Earth's ionosphere, Space Sci. Rev., 8, 846–865,, 1968. a

Chilton, C. J., Crombie, D. D., and Jean, A. G.: Phase Variations in V.L.F. Propagation (Chapter 19), in: Propagation of Radio Waves at Frequencies Below 300 kc/s: Proceedings of the Seventh meeting of the AGARD Ionospheric Research Committee, Munich 1962, edited by: Blackband, W. T., AGARDograph, Elsevier, Amsterdam, Netherlands, 74, 257–290,, 1964. a

Choi, Y., Moon, Y. J., Choi, S., Baek, J. H., Kim, S. S., Cho, K. S., and Choe, G. S.: Statistical Analysis of the Relationships among Coronal Holes, Corotating Interaction Regions, and Geomagnetic Storms, Sol. Phys., 254, 311–323,, 2009. a

Choudhury, A., De, B. K., Guha, A., and Roy, R.: Long-duration geomagnetic storm effects on the D region of the ionosphere: Some case studies using VLF signal, J. Geophys. Res.-Space, 120, 778–787,, 2015. a, b, c

Clilverd, M. A., Seppala, A., Rodger, C. J., Thomson, N. R., Lichtenberger, J., and Steinbach, P.: Temporal variability of the descent of high-altitude NOX inferred from ionospheric data, J. Geophys. Res.-Space, 112, A09307,, 2007. a

Clilverd, M. A., Rodger, C. J., Thomson, N. R., Brundell, J. B., Ulich, T., Lichtenberger, J., Cobbett, N., Collier, A. B., Menk, F. W., Seppälä, A., Verronen, P. T., and Turunen, E.: Remote sensing space weather events: Antarctic-Arctic Radiation-belt (Dynamic) Deposition-VLF Atmospheric Research Konsortium network, Space Weather, 7, S04001,, 2009. a

Clilverd, M. A., Rodger, C. J., Neal, J. J., and Cresswell-Moorcock, K.: Remote sensing space weather events through ionospheric radio: The AARDDVARK network, 2014 XXXIth URSI General Assembly and Scientific Symposium (URSI GASS), 1–1,, 2014. a

Colwell, R. and Friend, A.: The D Region of the Ionosphere, Nature, 137, p. 782., 1936. a, b

Crombie, D. D.: Phase and Time Variations in VLF Propagation Over Long Distances, J. Res. NBS Radio Science, 68D, 1223–1224, 1964. a

Crombie, D. D.: Further Observations of Sunrise and Sunset Fading of Very-Low-Frequency Signals, Radio Sci., 1, 47–51,, 1966. a

Cummer, S. A., Inan, U. S., and Bell, T. F.: Ionospheric D region remote sensing using VLF radio atmospherics, Radio Sci., 33, 1781–1792,, 1998. a

Danilov, A. D. and Smirnova, N. V.: Improving the 75 to 300 km ion composition model of the IRI, Adv. Space Res., 15, 171–177,, 1995. a

Davies, K. and Hartmann, G. K.: Studying the ionosphere with the Global Positioning System, Radio Sci., 32, 1695–1703,, 1997. a

Dellinger, J. H.: Sudden ionospheric disturbances, Terr. Magn. Atmos. Electr., 42, 49–53,, 1937. a

Dickinson, P. H. G. and Bennett, F. D. G: Diurnal variations in the D-region during a storm after-effect, J. Atmos. Terr. Phys., 40, 549–558,, 1978. a, b

Dierckxsens, M., Tziotziou, K., Dalla, S., Patsou, I., Marsh, M. S., Crosby, N. B., Malandraki, O., and Tsiropoula, G.: Relationship between Solar Energetic Particles and Properties of Flares and CMEs: Statistical Analysis of Solar Cycle 23 Events, Sol. Phys., 290 841–874,, 2015. a

Dougherty, J. P. and Farley, D. T.: A theory of incoherent scattering of radio waves by a plasma, P. Roy. Soc. Lond. A, 259, 79–99,, 1961. a

Dougherty, J. P. and Farley, D. T.: A theory of incoherent scattering of radio waves by a plasma: 3. Scattering in a partly ionized gas, J. Geophys. Res., 68, 5473–5486,, 1963. a

Dungey, J. W.: Interplanetary Magnetic Field and the Auroral Zones, Phys. Rev. Lett., 6, 47–48,, 1961. a

Eccles, J. V., Hunsucker, R. D., Rice, D., and Sojka, J. J.: Space weather effects on midlatitude HF propagation paths: Observations and a data-driven D region model, Space Weather, 3, S01002,, 2005. a, b

Evans, J. V.: Millstone Hill Thomson scatter results for 1965, Tech. Rep. 474, Lincoln Lab., Mass. Inst. of Technol., Cambridge, 8 December, (last access: 18 May 2022), 1969a. a

Evans, J. V.: Theory and practice of ionosphere study by Thomson scatter radar, Proc. IEEE, 57, 496–530,, 1969b. a

Fagundes, P. R., Cardoso, F. A., Fejer, B. G., Venkatesh, K., Ribeiro, B. A. G., and Pillat, V. G.: Positive and negative GPS‐TEC ionospheric storm effects during the extreme space weather event of March 2015 over the Brazilian sector, J. Geophys. Res.-Space, 121, 5613–5625,, 2016. a

Fairfield, D. H., Average and unusual locations of the Earth's magnetopause and bow shock, J. Geophys. Res., 76, 6700–6716,, 1971. a

Farley, D. T., Dougherty, J. P., and Barron, D. W.: A Theory of Incoherent Scattering of Radio Waves by a Plasma II. Scattering in a Magnetic Field, P. Roy. Soc. London A, 263, 238–258,, 1961. a

Fejer, B. G., Larsen, M. F., and Farley, D. T.: Equatorial disturbance dynamo electric fields, Geophys. Res. Lett., 10, 537–540,, 1983. a

Forbush, S. E: World-wide cosmic ray variations, 1937–1952, J. Geophys. Res., 59, 525–542,, 1954. a

Friedrich, M. and Torkar, K. M.: An empirical model of the non-auroral D Region, Radio Sci., 27, 945–953,, 1992. a

Fuller-Rowell, T. J., Codrescu, M. V., Moffett, R. J., and Quegan, S.: Response of the thermosphere and ionosphere to geomagnetic storms, J. Geophys. Res., 99, 3893–3914,, 1994. a

Ganushkina, N. Y., Liemohn, M. W., and Dubyagin, S.: Current systems in the Earth's magnetosphere, Rev. Geophys., 56, 309–332,, 2018. a

Gardiner, G. W.: Origin of the Term Ionosphere, Nature, 224, 1096,, 1969. a

Gauss, C. F.: Allgemeine Theorie des Erdmagnetismus, in: Resultate aus den Beobachtungen des magnetischen Vereins im Jahre 1838, edited by: Gauss, C. F. and Weber, W., Weidmannsche Buchhandlung, Leipzig, 1–57, (last access: 18 May 2022), 1838. a

George, H., Kilpua, E., Osmane, A., Asikainen, T., Kalliokoski, M. M. H., Rodger, C. J., Dubyagin, S., and Palmroth, M.: Outer Van Allen belt trapped and precipitating electron flux responses to two interplanetary magnetic clouds of opposite polarity, Ann. Geophys., 38, 931–951,, 2020. a, b

GFZ: German Research Center for Geosciences, Potsdam (Germany),, last access: 1 June 2022. a

Glassmeier, K.-H. and Tsurutani, B. T.: Carl Friedrich Gauss – General Theory of Terrestrial Magnetism – a revised translation of the German text, Hist. Geo Space. Sci., 5, 11–62,, 2014. a

Gonzalez, W. D., Joselyn, J. A., Kamide, Y., Kroehl, H. W., Rostoker, G., Tsurutani, B. T., and Vasyliunas, V. M.: What is a geomagnetic storm?, J. Geophys. Res.-Space, 99, 5771–5792,, 1994. a, b

Gonzalez, W. D., Tsurutani, B. T., and Clúa de Gonzalez, A. L.: Interplanetary origin of geomagnetic storms, Space Sci. Rev., 88, 529–562,, 1999. a, b

Gopalswamy N., Chapter 2 – Extreme Solar Eruptions and their Space Weather Consequences, in: Extreme Events in Geospace, edited by: Buzulukova, N., Elsevier Publishing Co., Amsterdam, Netherlands, 37–63,, 2018. a

Grafe, A., Lauter, E.-A., Nikutowski, B., and Wagner, C. U.: Precipitation of Energetic Electrons into the Mid-Latitude Ionosphere After Geomagnetic Storms, edited by: Rycroft, M. J., COSPAR Colloquia Series, Pergamon, Oxford, U.K., 20, 157–162,, 1980. a

Greenwald, R. A.: History of the Super Dual Auroral Radar Network (SuperDARN)-I: pre-SuperDARN developments in high frequency radar technology for ionospheric research and selected scientific results, Hist. Geo Space. Sci., 12, 77–93,, 2021. a

Greenwald R. A., Baker, K. B., Dudeney, J. R., Pinnock, M., Jones, T. B., Thomas, E. C., and Yamagishi, H.: Darn/superdarn, Space Sci. Rev., 71, 761–796,, 1995. a

Greer, K. R., Immel, T., and Ridley, A.: On the variation in the ionospheric response to geomagnetic storms with time of onset, J. Geophys. Res.-Space, 122, 4512–4525,, 2017. a

Gross, N. C. and Cohen, M. B.: VLF remote sensing of the D region ionosphere using neural networks, J. Geophys. Res.-Space, 125, e2019JA027135,, 2020. a

Gu, T. T. and Xu, H. L.: Mode Interferences of VLF Waves in the Presence of an Anisotropic Terrestrial Waveguide, in: Electromagnetic Propagation and Waveguides in Photonics and Microwave Engineering, edited by: Steglich, P., 23, InTechOpen Limited, London,, 2020. a

Gu, X., Xia, S., Fu, S., Xiang, Z., Ni, B., Guo, J., and Cao, X.: Dynamic Responses of Radiation Belt Electron Fluxes to Magnetic Storms and their Correlations with Magnetospheric Plasma Wave Activities, Astrophys. J., 891, 127,, 2020. a

Guerrero, A., Palacios, J., Rodríguez-Bouza, M., Rodríguez-Bilbao, I., Aran, A., Cid, C., Herraiz, M., Saiz, E., Rodrigues-Caderot, G., and Cerrato, Y.: Storm and substorm causes and effects at midlatitude location for the St. Patrick's 2013 and 2015 events, J. Geophys. Res.-Space, 122, 9994–10011,, 2017. a

Hajra, R.: September 2017 Space-Weather Events: A Study on Magnetic Reconnection and Geoeffectiveness, Sol. Phys., 296, 50,, 2021. a

Hayes, L. A., O'Hara, O. S. D., Murray, S. A., and Gallagher, P. T.: Solar Flare Effects on the Earth's Lower Ionosphere, Sol. Phys., 296, 157,, 2021. a

Heaviside, O.: Telegraphy, in: Encyclopaedia Britannica, 214, Edinburgh and London, Adam and Charles Black, 1902. a

Heelis, R. A. and Maute, A.: Challenges to understanding the Earth's ionosphere and thermosphere, J. Geophys. Res.-Space, 125, e2019JA027497,, 2020. a

Hegde, S., Bobra, M. G., and Scherrer, P. H.: Classifying Signatures of Sudden Ionospheric Disturbances, Research Notes of the AAS, 2, 162,, 2018. a

Horne, R. B., Lam, M. M., and Green, J. C.: Energetic electron precipitation from the outer radiation belt during geomagnetic storms, Geophys. Res. Lett., 36, L19104,, 2009. a

Hui, D., Chakrabarty, D., Sekar, R., Reeves, G. D., Yoshikawa, A., and Shiokawa, K.: Contribution of storm time substorms to the prompt electric field disturbances in the equatorial ionosphere, J. Geophys. Res.-Space, 122, 5568–5578,, 2017. a

Hunsucker, R. D.: Auroral and polar-cap ionospheric effects on radio propagation, IEEE Trans. Antennas Propagation, 40, 818–828,, 1992. a

Immel, T. J. and Mannucci, A. J.: Ionospheric redistribution during geomagnetic storms, J. Geophys. Res.-Space, 118, 7928–7939,, 2013. a

Inan, U. S., Cummer, S. A., and Marshall, R. A.: A survey of ELF and VLF research on lightning-ionosphere interactions and causative discharges, J. Geophys. Res., 115, A00E36,, 2010. a, b

Jackson, J. E.: Alouette-ISIS Program Summary, NSSDC/WDC-A-R and S 86-09, 94, (last access: 18 May 2022), 1986. a

Janvier, M., Demoulin, P., Guo, J., Dasso, S., Regnault, F., Topsi-Moutesidou, S., Gutierrez, C., and Perri, B.: The Two-step Forbush Decrease: A Tale of Two Substructures Modulating Galactic Cosmic Rays within Coronal Mass Ejections, Astrophys. J., 922, 216,, 2021. a

Johnson, C. Y.: Ionospheric composition and density from 90 to 1200 kilometers at solar minimum, J. Geophys. Res., 71, 330–332,, 1966. a

Jordanova, V. K., Ilie, R., and Chen, M. W.: Introduction and historical background, in: Ring Current Investigations (The Quest for Space Weather Prediction), edited by: Jordanova, V. K., Ilie, R., and Chen, M. W., Elsevier, Amsterdam, ISBN: 9780128155714, 1–13,, 2020. a

Kane, R. P.: Ionospheric foF2 anomalies during some intense geomagnetic storms, Ann. Geophys., 23, 2487–2499,, 2005. a

Kelley, M. C.: The Earth's Ionosphere, Plasma Physics and Electrodynamics, Academic Press, Cambridge, Massachusetts, U.S.A., 576 pp., ISBN: 978012088425, 2009. a

Kennelly, A.: On the Elevation of the Electrically-Conducting Strata of the Earth's Atmosphere, Elec. World and Engr., 39, p. 473, 1902. a

Kanekal, S. and Miyoshi, Y.: Dynamics of the terrestrial radiation belts: a review of recent results during the VarSITI (Variability of the Sun and Its Terrestrial Impact) era, 2014–2018, Prog. Earth Planet Sci., 8, 22,, 2021. a

Kerrache, F., Amor, S. N., and Kumar, S.: Ionospheric D region disturbances due to FAC and LEP associated with three severe geomagnetic storms as observed by VLF signals, J. Geophys. Res.-Space, 126, e2020JA027838,, 2021. a

Kikuchi, T. and Evans, D. S.: Quantitative study of substorm-associated VLF phase anomalies and precipitating energetic electrons on November 13, 1979, J. Geophys. Res., 88, 871–880,, 1983. a

Kilfoyle, B. and Jacka, F.: Geomagnetic L Coordinates, Nature, 220, 773–775,, 1968. a

Kim, R. S., Cho, K. H., Kim, Y. D., Park, Y. J., Moon, Y. J., Yi, Y., Lee, J., Wang, H., Song, H., and Dryer, M.: CME Earthward Direction as an Important Geoeffectiveness Indicator, Astrophys. J., 677, 1378–1384,, 2008. a

Kleimenova, N. G., Kozyreva, O. V., Rozhnoy, A. A., and Soloveva, M. S.: Variations in the VLF signal parameters on the Australia-Kamchatka radio path during magnetic storms, Geomagn. Aeron. 44, 385–393, 2004. a

Kovács, T., Plane, J. M. C., Feng, W., Nagy, T., Chipperfield, M. P., Verronen, P. T., Andersson, M. E., Newnham, D. A., Clilverd, M. A., and Marsh, D. R.: D-region ion–neutral coupled chemistry (Sodankylä Ion Chemistry, SIC) within the Whole Atmosphere Community Climate Model (WACCM 4) – WACCM-SIC and WACCM-rSIC, Geosci. Model Dev., 9, 3123–3136,, 2016. a

Kulyamin, D. V. and Dymnikov, V. P.: Numerical modelling of coupled neutral atmospheric general circulation and ionosphere D region, Russ. J. Numer. Anal. M., 31, 159–171,, 2016. a

Kumar, A. and Kumar, S.: Space weather effects on the low latitude D-region ionosphere during solar minimum, Earth Planet. Space, 66, 76,, 2014. a, b, c

Kumar, A. and Kumar, S.: Ionospheric D region parameters obtained using VLF measurements in the South Pacific region, J. Geophys. Res.-Space, 125, e2019JA027536,, 2020. a

Kumar, S., Kumar, A., Menk, F., Maurya, A. K., Singh, R., and Veenadhari, B.: Response of the low-latitude D region ionosphere to extreme space weather event of 14–16 December 2006, J. Geophys. Res.-Space, 120, 788–799,, 2015. a, b, c

Kutiev, I., Tsagouri, I., Perrone, L., Pancheva, D., Mukhtarov, P., Mikhailov, A., Lastovicka, J., Jakowski, N., Buresova, D., Blanch, E., Andonov, B., Altadill, D., Magdaleno, S., Parisi, M., and Torta, J. M.: Solar activity impact on the Earth's upper atmosphere, J. Space Weather Space Clim., 3,, 2013. a

Lanzagorta, M.: Underwater Communications (Synthesis Lectures on Communications), Morgan and Claypool Publishers, San Rafael, California (USA), 129,, 2012. a

Lauter, E. A. and Knuth, R. H.: Precipitation of high energy particles into the upper atmosphere at medium latitudes after magnetic storms, J. Atmos. Sol.-Terr. Phys., 29, 411–417,, 1967. a

Lastovicka, J.: Effects of geomagnetic storms in the lower ionosphere, middle atmosphere and troposphere, J. Atmos. Sol.-Terr. Phys., 58, 831–843,, 1996. a, b

Laughlin, L. K., Turner, N. E., and Mitchell, E. J. J.: Geoeffectiveness of CIR and CME Events: Factors Contributing to Their Differences, Southeast. Assoc. Res. Astron., 2, 19–22, 2008. a

Le, H., Liu, L., Ren, Z., Chen, Y., Zhang, H., and Wan, W.: A modeling study of global ionospheric and thermospheric responses to extreme solar flare, J. Geophys. Res.-Space, 121, 832–840,, 2008. a

Lincoln, J. V.: Geomagnetic and solar data, J. Geophys. Res., 66, 1561–1563,, 1961. a

Lincoln, J. V.: The listing of sudden ionospheric disturbances, Planet. Space Sci., 12, 419–434,, 1964. a

Liu, S. L. and Li, L. W.: Study on Relationship between Southward IMF Events and Geomagnetic Storms, Chinese J. Geophys., 45, 301–310,, 2002. a

Lodge, O.: Mr. Marconi's Results in Day and Night Wireless Telegraphy, Nature, 66, 1705,, 1902. a

Lynn, K. J. W.: Some differences in diurnal phase and amplitude variations for VLF signals, J. Atmos. Terr. Phys., 40, 145–150,, 1978. a

Machol, J., Snow, M., Woodraska, D., Woods, T., Viereck, R., and Coddington. O.: An improved lyman-alpha composite, Earth Space Sci., 6, 2263–2272,, 2019. a

Mangla and Yadov: Ionospheric Layers: D, E, F, F1, F2, Regions, (last access: 30 November 2021), 2011. a, b

Mannucci, A. J., Wilson, B. D., Yuan, D. N., Ho, C. H., Lindqwister, U. J., and Runge, T. F.: A global mapping technique for GPS-derived ionospheric total electron content measurements, Radio Sci., 33, 565–582,, 1998. a

Mannucci, A. J., Ao, C. O., and Williamson, W.: GNSS Radio Occultation, in: Position, Navigation, and Timing Technologies in the 21st Century, edited by: Morton, Y. T. J., Diggelen, F., Spilker, J. J., Parkinson, B. W., Lo, S., and Gao, G., chap. 33, 971–1013,, 2020. a

Marconi, G.:, Syntonic Wireless Telegraphy, J. Soc. Arts, 49, 505–520, 1901. a

Marr, G. V.: The penetration of solar radiation into the atmosphere, P. Roy. Soc. Lond. A, 288, 531–539,, 1965. a

Maurya, A. K., Venkatesham, K., Kumar, S., Singh, R., Tiwari, P., and Singh, A. K.: Effects of St. Patrick's Day geomagnetic storm of March 2015 and of June 2015 on low-equatorial D region ionosphere, J. Geophys. Res.-Space, 123, 6836–6850,, 2018. a, b

Mayaud, P. N: Derivation, Meaning, and Use of Geomagnetic Indices, Geophys. Monogr. Ser., 22, ISBN: 9780875900223,, 1980. a

McCormick, J. C. and Morris, M. B.: D region Ionospheric Imaging Using VLF/LF Broadband Sferics, Forward Modeling, and Tomography, 25th International Lightning Detection Conference and 7th International Lightning Meteorology Conference, 12–15 March 2018, Ft. Lauderdale, FL, 2018. a

McIlwain, C. E.: Coordinates for mapping the distribution of magnetically trapped particles, J. Geophys. Res., 66, 3681–3691,, 1961. a

McPherron, R. L.: Magnetospheric substorms, Rev. Geophys., 17, 657–681,, 1979. a, b, c

McPherron, R., Weygand, J., and Hsu, T.-S.: Response of the Earth's magnetosphere to changes in the solar wind, J. Atmos. Sol.-Terr. Phys., 70, 303–315,, 2008. a

McRae, W. M. and Thomson, N. R.: VLF phase and amplitude: daytime ionospheric parameters, J. Atmos. Sol.-Terr. Phys., 62, 609–618,, 2000. a

McRae, W. M. and Thomson, N. R.: Solar flare induced ionospheric D-region enhancements from VLF phase and amplitude observations, J. Atmos. Sol.-Terr. Phys., 66, 77–87,, 2004. a

Mironova, I., Sinnhuber, M., Bazilevskaya, G., Clilverd, M., Funke, B., Makhmutov, V., Rozanov, E., Santee, M. L., and Sukhodolov, T.: Exceptional middle latitude electron precipitation detected by balloon observations: implications for atmospheric composition, Atmos. Chem. Phys. Discuss. [preprint],, in review, 2021. a

Mitra, W. B.: Ionospheric effects of solar flares, D. Reidel Publishing Company, Dordrecht, Holland, ISBN 978-94-010-2233-0, 1974. a, b, c

Mittal, N., Gupta, A., Negi, P. S., and Narain, U.: On Some Properties of SEP Effective CMEs, International Scholarly Research Notice, 2011, 727140,, 2011. a

Miyoshi, Y., Kurita, S., Oyama, S.-I., Ogawa, Y., Saito, S., Shinohara, I., Kero, A., Turunen, E., Verronen, P. T., Kasahara, S., Yokota, S., Mitani, T., Takashima, T., Higashio, N., Kasahara, Y., Matsuda, S., Tsuchiya, F., Kumamoto, A., Matsuoka, A., Hori, T., Keika, K., Shoji, M., Teramoto, M., Imajo, S., Jun, C., and Nakamura, S.: Penetration of MeV electrons into the mesosphere accompanying pulsating aurorae, Sci. Rep., 11, 13724,, 2021. a

Moler, W. F.: VLF propagation effects of a D-region layer produced by cosmic rays, J. Geophys. Res., 65, 1459–1468,, 1960. a

Moore, R. K.: Radio communication in the sea, IEEE Spectrum, 4, 42–51,, 1967. a

Moral, A. C., Eyiguler, E. C. K., and Kaymaz, Z.: Sudden Ionospheric Disturbances and their detection over Istanbul, 2013 6th International Conference on Recent Advances in Space Technologies (RAST), Istanbul, Turkey, 12–14 June 2013, 765–768,, 2013. a

Muraoka, Y.: Lower ionospheric disturbances observed in long-distance VLF transmission at middle latitude, J. Atmos. Terr. Phys., 41, 1031–1042,, 1979. a

Naidu, P. P., Madhavilatha, T., and Devi, M. I.: Influence of geomagnetic storms on the mid latitude D and F2 regions, Ann. Geophys., 63, GM214,, 2020. a, b, c

NASA: National Space and Aeronautics Administration, Goddard Space Flight Center, Greenbelt, Maryland, USA,, last access: 15 May 2022. a

Nava, B., Rodriguez-Zuluaga, J., Alazo-Cuartas, K., Kashcheyev, A., Migoya-Orue, Y., Radicella, S. M., Amory-Mazaudier, C., and Fleury, R.: Middle- and low-latitude ionosphere response to 2015 St. Patrick's Day geomagnetic storm, J. Geophys. Res.-Space, 121, 3421–3438, 2016. a

Neal, J. J., Rodger, C. J., and Green, J. C.: Empirical determination of solar proton access to the atmosphere: Impact on polar flight paths, Space Weather, 11, 420–433,, 2013. a

Neal, J. J., Rodger, C. J., Thomson, N. R., Clilverd, M. A., Raita, T., and Ulich, T.: Long-term determination of energetic electron precipitation into the atmosphere from AARDDVARK subionospheric VLF observations, J. Geophys. Res.-Space, 120, 2194–2211,, 2015. a

Nicolet, M. and Aikin, A. C.: The formation of the D region of the ionosphere, J. Geophys. Res., 65, 1469–1483,, 1960. a, b

Nina, A., Nico, G., Mitrović, S. T., Čadež, V. M., Milošević, I. R., Radovanović, M., and Popović, L. Č.: Quiet Ionospheric D-Region (QIonDR) Model Based on VLF/LF Observations, Remote Sensing, 13,, 2021. a

NOAA: Geomagnetic Storms, NOAA Space Weather Prediction Center, (last access: 18 May 2022), 2016. a

Nunn, D., Clilverd, M. A., Rodger, C. J., and Thomson, N. R.: The impact of PMSE and NLC particles on VLF propagation, Ann. Geophys., 22, 1563–1574,, 2004. a

Nwankwo, V. U. J. and Chakrabarti, S. K.: Effects of space weather on the ionosphere and LEO satellites' orbital trajectory in equatorial, low and middle latitude, Adv. Space Res., 61, 1880–1889,, 2018. a

Nwankwo, V. U. J., Chakrabarti, S. K., and Ogunmodimu, O.: Probing geomagnetic storm-driven magnetosphere-ionosphere dynamics in D-region via propagation characteristics of very low frequency radio signals, J. Atmos. Sol.-Terr. Phys., 145, 154–169,, 2016. a, b, c, d, e, f, g, h, i, j, k, l, m, n

Nwankwo, V. U. J., Chakrabarti, S. K., Sasmal, S., Denig, W., Ajakaiye, M. P., Akinsola, T., Adeyanju, M., Anekwe, P., Iluore, K., Olatunji, M., Bhowmick, D., Fatokun, J., Ayoola, M. A., Soneye, O. O., and Ajamu, J.: Radio aeronomy in Nigeria: First results from very low frequency (VLF) radio waves receiving station at Anchor University, Lagos, 2020 IEEE-ICMCECS, Lagos, Nigeria, 1–7,, 2020. a, b

Nwankwo, V. U. J., Raulin, J.-P., Correia, E., Denig, W., Folarin, O., Ogunmodimu, O., and De Oliveira, R. R.: Investigation of ionosphere response to geomagnetic storms over the propagation paths of very low frequency radio waves, Radio Sci., 2021RS0007331, submitted, 2022. a, b, c, d, e

Palit, S., Basak, T., Mondal, S. K., Pal, S., and Chakrabarti, S. K.: Modeling of very low frequency (VLF) radio wave signal profile due to solar flares using the GEANT4 Monte Carlo simulation coupled with ionospheric chemistry, Atmos. Chem. Phys., 13, 9159–9168,, 2013. a

Pavlov, A. V.: Ion Chemistry of the Ionosphere at E- and F-Region Altitudes: A Review, Surv. Geophys., 33, 1133–1172,, 2012. a

Pederick, L. H. and Cervera M. A.: Semiempirical Model for Ionospheric Absorption based on the NRLMSISE‐00 atmospheric model, Radio Sci., 49, 81–93,, 2014. a

Pedersen, A.: Time, height, and latitude distribution of D layers in the subauroral zone and their relation to geomagnetic activity and aurora, J. Geophys. Res., 67, 2685–2694,, 1962. a

Peter, W. B., Chevalier, M. W., and Inan, U. S.: Perturbations of mid-latitude sub-ionospheric VLF signals associated with lower ionospheric disturbances during major geomagnetic storms, J. Geophys. Res., 111, AO3301,, 2006. a, b

Pierce, E. T.: Sferics, Conference Proceedings, NAS-NRC Atmospheric Exploration by Remote Probes, Vol. 2, NTRS 19720017733, 25, 1969. a

Pierce, J. A.: The Diurnal Carrier-Phase Variation of a 16-Kilocycle Transatlantic Signal, Proc. IRE, 43, 584–588,, 1955. a

Poole, I.: Radio Waves and the Ionosphere, in QST (Calling All Stations), November 1999, monthly publication of the American Radio Relay League (ARRL), 3 pp., (last access: 13 December 2021), 1999. a

Potemra, T. A., Zmuda, A. J., Shaw, B. W., and Have, C.R.: VLF phase disturbances, HF absorption, and solar protons in the PCA events of 1967, Radio Sci., 5, 1137–1145,, 1970. a

Porazik, P., Johnson, J. R., Kaganovich, I., and Sanchez, E.: Modification of the loss cone for energetic particles, Geophys. Res. Lett., 41, 8107–8113,, 2014. a

Prol, F. S., Kodikara, T., Hoque, M. M., and Borries, C.: Global-scale ionospheric tomography during the March 17, 2015 geomagnetic storm, Space Weather, 19, e2021SW002889,, 2021. a

Prolss, G. W.: Physics of the Earth's space environment, Springer Berlin Heidelberg, Germany, ISBN: 978-3-540-21426-7, 159–208,, 2004. a

Quan, L., Cai, B., Hu, X., Xu, Q., and Li, L.: Study of ionospheric D region changes during solar flares using MF radar measurements, Adv. Space Res., 67, 715–721,, 2021. a

Raghav, A., Shaikh, Z., Misal, D., Rajan, G., Mishra, W., Kasthurirangan, S., Bhaskar, A., Bijewar, N., Johri, A., and Vichare, G.: Exploring the common origins of the Forbush decrease phenomenon caused by the interplanetary counterpart of coronal mass ejections or corotating interaction regions, Phys. Rev. D, 101, 062003,, 2020. a

Rawer, K.: Introduction to IRI 1979, in International Reference Ionosphere IRI 97, UAG-82, World Data Center A for Solar-Terrestrial Physics, edited by: Lincoln, J. V. and Conkright, R. O., 1–6, (last access: 2 December 2021), 1981. a

Rawer, K., Bilitza, D., and Ramakrishnan, S.: Goals and status of the International Reference Ionosphere, Rev. Geophys., 16, 177–181,, 1978. a

Reeves, G. D. and Daglis, I. A.: Geospace Magnetic Storms and the Van Allen Radiation Belts, in: Waves, Particles, and Storms in Geospace: A Complex Interplay, chap. 3, edited by: Balasis, G., Daglis, I. A., and Mann, I. R.,, 2016. a

Reeves, G. D., Spence, H. E., Henderson, M. G., Morley, S. K., Friedel, R. H., Funsten, H. O., Baker, D. N., Kanekal, S. G., Blake, J. B., Fennell, J. F., Claudepierre, S. G., Thorne, R. M., Turner, D. L., Kletzing, C. A., Kurth, W. S., Larsen, B. A., and Niehof, J. T.: Electron acceleration in the heart of the Van Allen radiation belts, Science, 341, 6149,, 2013. a

Reinisch, B. W. and Xueqin, H.: Automatic calculation of electron density profiles from digital ionograms: 3. Processing of bottomside ionograms, Radio Sci., 18, 477–492,, 1983. a

Rishbeth, H.: Physics and chemistry of the ionosphere, Contemporary Physics, 14, 229–249,, 1973. a, b

Ries, G.: Results Concerning the Sunrise Effect of VLF Signals Propagated Over Long Paths, Radio Sci., 2, 531–538,, 1967. a

Ripoll, J.-F., Denton, M., Loridan, V., Santolík, O., Malaspina, D., Hartley, D. P., Cunningham, G. S., Reeves, G., Thaller, S., Turner, D. L., Fennell, J. F., Drozdov, A. Y., Cervantes Villa, J. S., Shprits, Y. Y., Chu, X., Hospodarsky, G., Kurth, W. S., Kletzing, C. A., Wygant, J., Henderson, M. G., and Ukhorshiy, A. Y.: How whistler mode hiss waves and the plasmasphere drive the quiet decay of radiation belts electrons following a geomagnetic storm, J. Phys. Conf. Ser., 1623, paper 012005, 14th Int. Conf. on Numerical Modeling of Space Plasma Flows: ASTRONUM-2019 1-5 July 2019, Paris, France, 2020. a

Robinson, R. M. and Zanetti, L. J.: Auroral energy flux and Joule heating derived from global maps of field-aligned currents, Geophys. Res. Lett., 48, e2020GL091527,, 2021. a

Robinson, R. M., van Eyken, A., and Farley, D.: Fiftieth Anniversary of the First Incoherent Scatter Radar Experiment, Eos T. AGU, 90, p. 267,, 2009. a

Rodger, C. J., Clilverd, M. A., Thomson, N. R., Gamble, R. J., Seppälä, A., Turunen, E., Meredith, N. P., Parrot, M., Sauvaud, J.-A., and Berthelier, J. J.: Radiation belt electron precipitation into the atmosphere: Recovery from a geomagnetic storm, J. Geophys. Res., 112, A11307,, 2007. a, b

Rodger, C. J., Clilverd, M. A., Seppälä, A., Thomson, N. R., Gamble, R. J., Parrot, M., Sauvaud, J.-A., and Ulich, T.: Radiation belt electron precipitation due to geomagnetic storms: Significance to middle atmosphere ozone chemistry, J. Geophys. Res., 115, A11320,, 2010. a, b

Rodger, C. J., Clilverd, M. A., Kavanagh, A. J., Watt, C. E. J., Verronen, P. T., and Raita, T.: Contrasting the responses of three different ground-based instruments to energetic electron precipitation, Radio Sci., 47, RS2021,, 2012. a, b

Rodger, C. J., Hendry, A. T., Clilverd, M. A., Kletzing, C. A., Brundell, J. B., and Reeves, G. D.: High‐resolution in situ observations of electron precipitation‐causing EMIC waves, Geophys. Res. Lett., 42, 9633–9641,, 2015. a

Rogers, N. and Honary, F.: D-region HF absorption models incorporating real-time riometer measurements, XXXIth URSI General Assembly and Scientific Symposium (URSI GASS), 2014, 1–2,, 2014. a, b

Rogers, N. C., Kero, A., Honary, F., Verronen, P. T., Warrington, E. M., and Danskin, D. W.: Improving the twilight model for polar cap absorption nowcasts, Space Weather, 14, 950–972,, 2016. a

Rose, D. C. and Ziauddin, S.: The Polar Cap Absorption Effect, Space Sci. Rev., 1, 115–134,, 1962. a

Rostoker, G., Akasofu, S. I., Foster, J., Greenwald, R., Kamide, Y., Kawasaki, K., Lui, A., McPherron, R., and Russell, C.: Magnetospheric substorms—definition and signatures, J. Geophys. Res., 85, 1663–1668,, 1980. a

Rozhnoi, A., Solovieva, M., Fedun, V., Gallagher, P., McCauley, J., Boudjada, M. Y., Shelyag, S., and Eichelberger, H. U.: Strong influence of solar X-ray flares on low-frequency electromagnetic signals in middle latitudes, Ann. Geophys., 37, 843–850,, 2019. a

Russell, C. T.: The Magnetosphere, Annu. Rev. Earth Pl. Sc., 19, 169–182,, 1991. a

Russell, C. T., McPherron, R. L., and Burton, R. K.: On the cause of geomagnetic storms, J. Geophys. Res., 79, 1105–1109,, 1974. a

Samanes, J. E., Raulin, J., Macotela, E. L., and Guevara Day, W. R.: Estimating the VLF modal interference distance using the South America VLF Network (SAVNET), Radio Sci., 50, 122–129,, 2015. a

Samsonov, A. A., Bogdanova, Y. V., Branduardi-Raymont, G., Sibeck, D. G., and Toth, G.: Is the relation between the solar wind dynamic pressure and the magnetopause standoff distance so straightforward?, Geophys. Res. Lett., 47, e2019GL086474,, 2020. a

Sasmal, S. and Chakrabarti, S. K.: Ionosperic anomaly due to seismic activities – Part 1: Calibration of the VLF signal of VTX 18.2 KHz station from Kolkata and deviation during seismic events, Nat. Hazards Earth Syst. Sci., 9, 1403–1408,, 2009. a

Sastri, J. H.: Effect of magnetic storms and substorms on the low- latitude/ equatorial ionosphere, in: Solar Influence on the Heliosphere and Earth's Environment: Recent Progress and Prospects, edited by: Gopalswamy, N. and Bhattacharyya, A., Proc. ILWS Workshop in Goa, India, 19–24 February 2006, (last access: 18 January 2022), 2006. a

Satori, G.: Combined ionospheric effect due to Forbush decreases and magnetospheric high energy particles at mid-latitudes, J. Atmos. Terr. Phys., 53, 325–332,, 1991. a, b

Sauer, H. H. and Wilkinson, D. C.: Global mapping of ionospheric HF/VHF radio wave absorption due to solar energetic protons, Space Weather, 6, S12002,, 2008. a, b, c

Sharma, A. K. and More, C.: Diurnal Variation of VLF Radio Wave Signal Strength at 19.8 and 24 kHz Received at Khatav India (16 46'N, 75 53'E), J. Space Sci. Tech., 6, 15–34,, 2017. a, b

Scherliess, L. and Fejer, B. G.: Storm time dependence of equatorial disturbance dynamo zonal electric fields, J. Geophys. Res.-Space, 102, 24037–24046,, 1997. a

Schunk, R. W.: Handbook of Ionospheric Models, Report, Aeronomic Models of the Ionosphere, Solar-Terrestrial Energy Program (STEP), Working Group 3.6, 301, 1996. a

Schunk, R. W.: Guide to Reference and Standard Ionosphere, ANSI/AIAA G-034-1998, American National Standards Institute, 67, 1999. a

Schunk, R. and Nagy, A.: Ionospheres: Physics, Plasma Physics, and Chemistry, 2nd Edn., Cambridge Atmospheric and Space Science Series, Cambridge University Press,, 2009. a

Schunk, R. W., Scherliess, L., Sojka, J. J., Thompson, D. C., Anderson, D. A., Codrescue, M., Minter, C., Fuller-Rowell, T. J., Heelis, R. A., and Howe, B. M.: Global Assimilation of Ionospheric Measurements (GAIM), Radio Sci., 39, RS1S0,, 2004. a

Scotto, C. and Settimi, A.: The calculation of ionospheric absorption with modern computers, Adv. Space Res., 54, 1642–1650,, 2014. a

Sechrist, C. F.: Comparisons of techniques for measurement of D‐region electron densities, Radio Sci., 9, 137–149,, 1974. a

Seppälä, A., Clilverd, M. A., Beharrell, M. J., Rodger, C. J., Verronen, P. T., Andersson, M. E., and Newnham, D. A.: Substorm-induced energetic electron precipitation: Impact on atmospheric chemistry, Geophys. Res. Lett., 42, 8172–8176,, 2015. a

Shue, J. H., Chao, J. K., Fu, H. C., Russell, C. T., Song, P., Khurana, K. K., and Singer, H. J.: A new functional form to study the solar wind control of the magnetopause size and shape, J. Geophys. Res.-Space, 102, 9497–9511,, 1997. a

Shue, J.-H., Song, P., Russell, C. T., Steinberg, J. T., Chao, J. K., Zastenker, G., Vaisberg, O. L., Kokubun, S., Singer, H. J., Detman, T. R., and Kawano, H.: Magnetopause location under extreme solar wind conditions, J. Geophys. Res.-Space Phy., 103, 17691–17700,, 1998. a

Silber, I. and Price, C.: On the use of VLF narrowband measurements to study the lower ionosphere and the mesosphere-lower thermosphere, Surv. Geophys., 38, 407–441,, 2017. a

Singh, R., Verma, U. P., and Singh A. K.: Exploring Middle Atmosphere (D-Region) by Very Low Frequency (VLF) Waves, International Journal of Latest Technology in Engineering, Management and Applied Science (IJLTEMAS), V, VI, 51–55, 2016. a

Siskind, D. E., Zawdie, K., Sassi, F., Drob, D., and Friedrich, M.: Global modeling of the low and mid latitude ionospheric D and lower E regions and implications for HF radio wave absorption, Space Weather, 15, 115–130,, 2017. a, b, c, d

Snay, R. A. and Soler, T.: Continuously Operating Reference Station (CORS): History, Applications, and Future Enhancements, J. Surv. Eng., 134, 95–104,, 2008. a

Spence, H. E.: The what, where, when, and why of magnetospheric substorm triggers, Eos T. AGU, 77, 81–86,, 1996. a

Spies, K. P. and Wait, J. R.: Mode calculations for VLF propagation in the earth-ionosphere waveguide, NBS Technical Note 114, 116, Call Number: QC100.U5753 no.114 1961, (last access: 10 April 2020), 1961. a

Spjeldvik, W. N. and Thorne, R. M.: A simplified D-region model and its application to magnetic storm after-effects, J. Atmos. Terr. Phys., 37, 1313–1325,, 1975. a, b

Stoker, P. H.: Energetic Electron Power Flux Deposition at Sanae (L=4.0) from Riometer Recording, J. Geophys. Res., 98, 19111–19116,, 1993. 

Stamper, R., Davis, C., and Bradford, V.: RAL Low-Cost Ionosonde System, paper SM21A-0361, AGU Fall Meeting, 5–9 December 2005, San Francisco, CA U.S.A., 2005. a

Soni, S. L., Yadav, M. L., Gupta, R. S., and Verma, P. L.: Exhaustive study of three-time periods of solar activity due to single active regions: sunspot, flare, CME, and geo-effectiveness characteristics, Astrophys. Space Sci., 365, 189,, 2020. a

Sun, K., Cui, W., and Chen, C.: Review of Underwater Sensing Technologies and Applications, Sensors, 21, 7849,, 2021. a

Suvorova, A. V. and Dmitriev, A. V.: Magnetopause inflation under radial IMF: Comparison of models, Earth Space Sci., 2, 107–114,, 2015. a

Takefu, M.: Bragg scattering of radio waves by ionospheric wavelike irregularities, J. Geomag. Geoelect., 41, 647–672, 1989. a

Tatsuta, K., Hobara, Y., Pal, S., and Balikhin, M.: Sub-ionospheric VLF signal anomaly due to geomagnetic storms: a statistical study, Ann. Geophys., 33, 1457–1467,, 2015. a, b

Taylor, W. L.: Daytime Attenuation Rates in the Very Low Frequency Band Using Atmospherics, J. Res. NBS, 64D, 349–355, 1960. a

Thomson, N. R., Clilverd, M. A., and McRae, W. M.: Nighttime ionospheric D-region parameters from VLF phase and amplitude, J. Geophys. Res., 112, A07304,, 2007. a, b

Thomson, N. R.: Experimental daytime VLF ionospheric parameters, J. Atmos. Terr. Phys., 55, 173–184,, 1993. a

Thomson, N. R. and Clilverd, M. A.: Solar flare induced ionospheric D-region enhancements from VLF amplitude observations, J. Atmos. Sol.-Terr. Phys., 63, 1729–1737,, 2001. a, b

Thomson, N. R. and McRae, W. M.: Nighttime ionospheric D region: Equatorial and nonequatorial, J. Geophys. Res., 114, A08305,, 2009. a

Timocin, E.: Swarm satellite observations of the effect of prompt penetration electric fields (PPEFs) on plasma density around noon and midnight side of low latitudes during the 07–08 September 2017 geomagnetic storm, Adv. Space Res., 69, 1335–1343,, 2022. a

Tsurutani, B. T., Gonzalez, W. D., Gonzalez, A. L. C., Tang, F., Arballo, J. K., and Okada, M.: Interplanetary origin of geomagnetic activity in the declining phase of the solar cycle, J. Geophys. Res., 100, 21717–21733,, 1995. a

Tsurutani, B. T., Gonzalez, W. D., Lakhina, G. S., and Alex, S.: The extreme magnetic storm of 1-2 September 1859, J. Geophys. Res.-Space, 108, 1268,, 2003. a

Tsurutani, B. T., Gonzalez, W. D., Gonzalez, A. L. C., Guarnieri, F. L., Gopalswamy, N., Grande, M., Kamide, Y., Kasahara, Y., Lu, G., Mann, I., McPherron, R. L., Soraas, F., and Vasyliunas, V. M.: Corotating solar wind streams and recurrent geomagnetic activity: A review, J. Geophys. Res., 111, A07S01,, 2006. a

Tsurutani, B. T., Verkhoglyadova, O. P., Mannucci, A. J., Saito, A., Araki, T., Yumoto, K., Tsuda, T., Abdu, M. A., Sobral, J. H. A., Gonzalez, W. D., McCreadie, H., Lakhina, G. S., and Vasyliūnas, V. M.: Prompt penetration electric fields (PPEFs) and their ionospheric effects during the great magnetic storm of 30–31 October 2003, J. Geophys. Res., 113, A05311,, 2008. a

Tsurutani, B. T., Echer, E., Guarnieri, F. L., and Gonzalez, W. D.: The properties of two solar wind high speed streams and related geomagnetic activity during the declining phase of solar cycle 23, J. Atmos. Sol.-Terr. Phys., 73, 164,, 2011. a, b

Turner, N. E., Mitchell, E. J., Knipp, D. J., and Emery, B. A.: Energetics of Magnetic Storms Driven by Corotating Interaction Regions: A Study of Geoeffectiveness, in: Recurrent Magnetic Storms: Corotating Solar Wind Streams, edited by: Tsurutani, B., McPherron, R., Lu, G., Sobral, J. H. A., and Gopalswamy, N.,, 2006. a

Turner, D. L., Angelopoulos, V., Li, W., Hartinger, M. D., Usanova, M., Mann, I. R., and Shprits, Y.: On the storm‐time evolution of relativistic electron phase space density in Earth's outer radiation belt, J. Geophys. Res.-Space, 118, 2196–2212, 2013. a

Turunen, E., Matveinen, H., Tolvanen, J., and Ranta, H.: D-region ion chemistry model, in STEP Handbook of Ionospheric Models, edited by: Schunk, R. W., SCOSTEP Secretariat, Boulder, Colo., reference STEP Handbook, 1–25, 1996. a

Turunen, E., Verronen, P. T., Seppälä, A., Rodger, C. J., Clilverd, M. A., Tamminen, J., Enell, C.-F., and Ulich, V.: Impact of different energies of precipitating particles on NOx generation in the middle and upper atmosphere during geomagnetic storms, J. Atmos. Sol.-Terr. Phys., 71, 1176–1189,, 2009. a, b, c

Turunen, E., Kero, A., Verronen, P. T., Miyoshi, Y., Oyama, S. I., and Saito, S.: Mesospheric ozone destruction by high-energy electron precipitation associated with pulsating aurora, J. Geophys. Res.-Atmos., 121, 11852–11861,, 2016. a

Tuve, M. A. and Breit, G.: Note on a radio method of estimating the height of the conducting layer, Terrestrial Magnetism and Atmospheric Electricity, 30, 15–16, 1925. a

UKSSDC: U.K. Solar System Data Centre, Rutherford Appleton Laboratory, Oxfordshire (UK),, last access: 15 May 2022. a

van Allen, J. A., Ludwig, G. H., Ray, E. C., and McIlwain, C. E.: Observation of High Intensity Radiation by Satellites 1958 Alpha and Gamma, Jet Propulsion, September 1958, American Rocket Society, Inc., 5,, 1958. a

Verbanac G., Vršnak, B., Živković, S., Hojsak, T., Veronig, A. M., and Temmer, M.: Solar wind high-speed streams and related geomagnetic activity in the declining phase of solar cycle 23, Astron. Astrophys., 533, 6,, 2011. a

Veronig, A., Temmer, M., Hanslmeier, A., Itruba, W., and Messerotti, M.: Temporal aspects and frequency distributions of solar soft X-ray flares, Astron. Astrophys., 382, 1070–1080,, 2002. a

Verronen, P. T., Seppala, A., Clilverd, M. A., Rodger, C. J., Kyrola, E., Enell, C. F., Ulich, T., and Turunen, E.: Diurnal variation of ozone depletion during the October–November 2003 solar proton events, J. Geophys. Res.-Space, 110, A09S32,, 2005. a

Verronen, P. T., Andersson, M. E., Marsh, D. R., Kovacs, T., and Plane, J. M. C.: WACCM-D—Whole Atmosphere Community Climate Model with D-region ion chemistry, J. Adv. Model. Earth Sy., 8, 954–975,, 2016. a, b

Voss, H. D., Walt, M., Imhof, W. L., Mobilia, J., and Inan, U. S.: Satellite observations of lightning-induced electron precipitation, J. Geophys. Res., 103, 11725–11744,, 1998. a

Waheed-uz-Zaman, M. and Yousufzai, M. A. K.: Design and Construction of Very Low Frequency Antenna, J. Basic Appl. Sci., 7, 141–145, 2011. a

Wait, J. R.: Terrestrial Propagation of Very-Low-Frequency Radio Waves, J. Res. Nat Bureau Standards, 64D, 153–204,, 1960. a, b

Wait, J. R.: A New Approach to the Mode Theory of VLF Propagation, J. Res. NBS: D. Radio Propagation, 65D, 37–46,, 1961. a

Wait, J. R.: A Note on Diurnal Phase Changes of Very‐Low‐Frequency Waves for Long Paths, J. Geophys. Res., 68, 338–340,, 1963. a

Wait, J. R.: Two-Dimensional Treatment of Mode Theory of the Propagation of VLF Radio Waves, J. Res. NBS, 68D, 1, 81–93,, 1964. a

Wait, J. R.: 1968 Mode conversion and refraction effects in the Earth‐ionosphere waveguide for VLF radio waves, J. Geophys. Res., 73, 3537–3548,, 1964. a

Wait, J. R.: Electromagnetic waves in stratified media, Pergamon Press, Oxford, ISBN: 978-0-08-006636-3, 609 pp.,, 1970. a

Wait, J. R. and Spies, K. P.: Characteristics of the Earth-ionosphere wave-guide for VLF radio waves, NBS Tech. Note 300, 1964. a, b, c, d, e

Walker, D.: Phase steps and amplitude fading of VLF signals at dawn and dusk, Radio Sci., 68D, 1435–1443,, 1965. a

WDCG (World Data Center for Geomagnetism): Data Analysis Center for Geomagnetism and Space Magnetism, Kyoto University (Japan),, last access: 18 May 2022.  a

Weigel, R. S.: Solar wind density influence on geomagnetic storm intensity, J. Geophys. Res., 115, A09201,, 2010. a

Wiltberger, M., Lopez, R. E., and Lyon, L. G.: Magnetopause erosion: A global view from MHD simulation, J. Geophys. Res., 108, 1235,, 2003. a

Wu, C. C., Liou, K., Lepping, R. P., Hutting, L., Plunkett, S., Howard, R. A., and Socker, D.: The first super geomagnetic storm of solar cycle 24: The St. Patrick's day event (17 March 2015), Earth Planet Space, 68, 151,, 2016. a

Yokoyama, E. and Tanimura, I.: Some Long-Distance Transmission Phenomena of Low-Frequency Waves, Proc. IRE, 21, 263–270,, 1933. a

Youssef, M.: On the relation between the CMEs and the solar flares, NRIAG J. Astron. Geophys., 1, 172–178,, 2012. a

Yue, X., Schreiner, W. S., Pedatella, N., Anthes, R. A., Mannucci, A. J., Straus, P. R., and Liu, J. Y.: Space Weather Observations by GNSS Radio Occultation: From FORMOSAT-3/COSMIC to FORMOSAT-7/COSMIC-2, Space Weather, 12, 616–621,, 2014. a

Zawedde, A. E., Nesse Tyssøy, H., Stadsnes, J., and Sandanger, M. I.: The impact of energetic particle precipitation on mesospheric OH – Variability of the sources and the background atmosphere, J. Geophys. Res.-Spac, 123, 5764–5789,, 2018. a, b

Zesta, E. and Oliveira, D. M.: Thermospheric heating and cooling times during geomagnetic storms, including extreme events, Geophys. Res. Lett., 46, 12739–12746,, 2019. a

Short summary
We combined the observed diurnal VLF amplitude variation in the D region with standard measurements of the E and F regions to perform a diagnostic investigation of coupled geomagnetic storm effects in order to understand the observed storm-induced variations in VLF narrowband based on state and responses of the ionosphere. The dayside VLF amplitude showed a tendency for attenuation following geomagnetic storms, and the h’E and h’F variations confirmed strong storm response over the signal paths.