Solar Cycle-Modulated Deformation of the Earth–Ionosphere Cavity

The Earth–ionosphere cavity resonator is occupied primarily by the electromagnetic radiation of lightning below 100 Hz. The phenomenon is known as Schumann resonances (SR). SR intensity is an excellent indicator of lightning activity and its distribution on global scales. However, long-term measurements from high latitude SR stations revealed a pronounced in-phase solar cycle modulation of SR intensity seemingly contradicting optical observations of lightning from satellite, which do not show any significant solar cycle variation in the intensity and spatial distribution of lightning activity on the global scale. The solar cycle-modulated local deformation of the Earth–ionosphere cavity by the ionization of energetic electron precipitation (EEP) has been suggested as a possible phenomenon that may account for the observed long-term modulation of SR intensity. Precipitating electrons in the energy range of 1–300 keV can affect the Earth–ionosphere cavity resonator in the altitude range of about 70–110 km and modify the SR intensities. However, until now there was no direct evidence documented in the literature supporting this suggestion. In this paper we present long-term SR intensity records from eight stations, each equipped with a pair of induction coil magnetometers: five high latitude (|lat| > 60°), two mid-high latitude (50° < |lat| < 60°) and one low latitude (|lat| < 30°). These long-term, ground-based SR intensity records are compared on the annual and interannual timescales with the fluxes of precipitating 30–300 keV medium energy electrons provided by the POES NOAA-15 satellite and on the daily timescale with electron precipitation events identified using a SuperDARN radar in Antarctica. The long-term variation of the Earth–ionosphere waveguide’s effective height, as inferred from its cutoff frequency, is independently analyzed based on spectra recorded by the DEMETER satellite. It is shown that to account for all our observations one needs to consider both the effect of solar X-rays and EEP which modify the quality factor of the cavity and deform it dominantly over low- and high latitudes, respectively. Our results suggest that SR measurements should be considered as an alternative tool for collecting information about and thus monitoring changes in the ionization state of the lower ionosphere associated with EEP.

The Earth-ionosphere cavity resonator is occupied primarily by the electromagnetic radiation of lightning below 100 Hz. The phenomenon is known as Schumann resonances (SR). SR intensity is an excellent indicator of lightning activity and its distribution on global scales. However, long-term measurements from high latitude SR stations revealed a pronounced in-phase solar cycle modulation of SR intensity seemingly contradicting optical observations of lightning from satellite, which do not show any significant solar cycle variation in the intensity and spatial distribution of lightning activity on the global scale. The solar cycle-modulated local deformation of the Earth-ionosphere cavity by the ionization of energetic electron precipitation (EEP) has been suggested as a possible phenomenon that may account for the observed long-term modulation of SR intensity. Precipitating electrons in the energy range of 1-300 keV can affect the Earth-ionosphere cavity resonator in the altitude range of about 70-110 km and modify the SR intensities. However, until now there was no direct evidence documented in the literature supporting this suggestion. In this paper we present longterm SR intensity records from eight stations, each equipped with a pair of induction coil magnetometers: five high latitude (|lat| > 60°), two mid-high latitude (50°< |lat| < 60°) and one low latitude (|lat| < 30°). These long-term, ground-based SR intensity records are compared on the annual and interannual timescales with the fluxes of precipitating 30-300 keV medium energy electrons provided by the POES NOAA-15 satellite and on the daily timescale with electron precipitation events identified using a SuperDARN radar in Antarctica. The long-term variation of the Earth-ionosphere waveguide's effective height, as inferred from its cutoff frequency, is independently analyzed based on spectra recorded by the DEMETER satellite. It is shown that to account for all our observations one needs to consider both the effect of solar X-rays and EEP which modify the quality factor of the cavity and deform it dominantly over low-and high latitudes, respectively. Our results

INTRODUCTION
The Earth's surface and the lower ionosphere together form an electromagnetic cavity resonator with resonance frequencies in the extremely low frequency (ELF, 3 Hz-3 kHz) band at ∼8, ∼14, ∼20, ∼26 Hz etc.. These resonances are named after Schumann (1952) who described them theoretically in 1952, 8 years before full frequency spectra were actually observed by Balser and Wagner (1960). Schumann resonances (SR) are excited primarily by lightning-radiated electromagnetic waves and the approximate number of 30-100 lightning strokes per second worldwide (Christian et al., 2003) maintains the resonance field quasi-continuously. It follows that SR are excellent indicators of lightning activity and distribution on global scales (see e.g., Williams et al., 2021). The possible connection between atmospheric electricity (including SR) and biological systems is another topic with substantial scientific interest (e.g., Elhalel et al., 2019;Fdez-Arroyabe et al., 2020;Price et al., 2021;Savoska et al., 2021;Sukhov et al., 2021). Until the end of the 1980's SR measurements were sparse and due to the restricted computational and data storage capacities long-term continuous recording was not really feasible at that time. In the early 90s' owing to the rapid development of large computational capacities and motivated by a fundamental Science article (Williams, 1992) many SR stations started to operate on a continuous basis worldwide (e.g., Price and Melnikov, 2004;Neska et al., 2007;Ondrášková et al., 2007;Bór et al., 2020) which gave a new impetus for SR-related research.
Much published material has appeared describing the standard (recurrent) daily and seasonal variations of SR spectral parameters (modal frequency, intensity and sharpness of maxima) at different stations (e.g., Sátori, 1996;Nickolaenko and Hayakawa, 2002;Price and Melnikov, 2004;Roldugin et al., 2004a;Ondrášková et al., 2007;Yatsevich et al., 2008;Toledo-Redondo et al., 2010;Zhou et al., 2013;Nickolaenko and Hayakawa, 2014;Fornieles-Callejón et al., 2015;Manu et al., 2015;Musur and Beggan, 2019;Tatsis et al., 2020). These standard variations in SR are mainly the consequence of the daily and seasonal variations of global lightning activity which fundamentally affect the observed SR spectral parameters (Sátori, 1996). On the other hand, standard variations in SR are also connected to the effect of the day-night terminator, i.e., to changes in the shape of the Earth-ionosphere cavity Sátori et al., 2007;Prácser et al., 2021). SR measurements demonstrated that the characteristic variations of global lightning activity are highly recurrent on a day-to-day and year-to-year basis (see e.g., Sátori, 1996). The latter was confirmed by optical observations of lightning from space as well (Williams et al., 2000;Christian et al., 2003).
Occasional variations in SR have also been documented, which are partly connected to transient extra-terrestrial phenomena that can modify the shape of the Earth-ionosphere cavity and thereby to "untune" its resonance frequencies. It is important to note that as a consequence of any change in the shape of the Earth-ionosphere cavity all the resonance frequencies shift in the same direction at all stations on Earth, but whether they increase or decrease depends on the actual modification of the conductivity profile in the upper atmosphere (Zhou and Qiao, 2015;Sátori et al., 2016;Kudintseva et al., 2018). Solar X-ray bursts (Roldugin et al., 2004b;Sátori et al., 2005;Dyrda et al., 2015;Sátori et al., 2016;Shvets et al., 2017) and solar proton events (Schlegel and Füllekrug, 1999;Roldugin et al., 2003;Singh et al., 2014) are known to cause such changes in frequencies (of both signs). In contrast, occasional variations in SR intensity with extra-terrestrial origin have been associated only with solar proton events so far (Schlegel and Füllekrug, 1999;Roldugin et al., 2003), and recently, indication for periods with increased SR intensity during geomagnetic storms was demonstrated as well (Salinas et al., 2016;Pazos et al., 2019). We note that variations in SR were also reported connected to seismic activity (e.g., Christofilakis et al., 2019;Galuk et al., 2019;Florios et al., 2020;Hayakawa et al., 2020).
SR measurements have also revealed the long-term solar cycle modulation of SR frequencies (e.g., Kulak et al., 2003a;Kulak et al., 2003b;Sátori et al., 2005;Ondrášková et al., 2011;Nickolaenko et al., 2015;Musur and Beggan, 2019;Koloskov et al., 2020) which has been attributed to the more than two orders of magnitude difference in solar X-ray flux between the minimum and maximum of solar activity (Sátori et al., 2005;Williams and Sátori, 2007). This effect is inherently non-uniform as only the sunlit portion of the lower ionosphere is affected by X-rays (and EUV radiation). It is to be noted that the long-term modulation of SR frequencies can be observed at all SR stations around the globe just as in the case of transient extra-terrestrial phenomena (X-ray bursts, solar proton events) discussed above.
In 2015 long-term SR intensity records from the Ukrainian Antarctic station "Akademik Vernadsky" were published which clearly showed a pronounced solar cycle modulation (up to 60%) of SR intensity (Nickolaenko et al., 2015), just as in an earlier paper by Füllekrug et al. (2002) using data from the Antarctic station Arrival Heights. Recently, Koloskov et al. (2020) reported that SR intensity records at the Vernadsky station are still in phase with the solar cycle and that a new SR station began operation at Svalbard in the Arctic in 2013 that also shows the solar cycle modulation of SR intensity. The origin of this effect is a challenging question as satellite observations do not show a solar cycle variation in the intensity and distribution of lightning activity (which is recognized as the primary origin of SR intensity variations) (Christian et al., 2003;Williams et al., 2014;Williams, 2016). Furthermore, as highlighted above, there is little clear evidence reported in the literature so far regarding the direct influence of extra-terrestrial processes on SR intensity over shorter timescales which could form the basis of the longerterm modulation from the solar cycle. In studies by Williams et al. (2014) and Sátori et al. (2016) the authors argued that energetic electron precipitation (EEP) may account for the observed interannual SR intensity modulation by affecting the ionospheric height locally over the observer. As for the physical basis the authors suggested (following Sentman and Fraser, 1991) that the energy flux through the waveguide, parallel with the Earth's surface, must be conserved and when the local height is reduced by EEP, the magnetic field and attendant intensity must increase locally (a scenario illustrated in Figure 9). However, until now no direct evidence has been presented in the literature supporting the suggested EEP-effect on SR either on the solar-cycle or on shorter timescales.
Energetic particles are trapped in the magnetosphere in radiation belts (RBs). They are heated to relativistic energies (from several tens of keV up to the order of MeV) by various wave particle interactions (WPI), field line reconnection processes and large-scale electric fields (see e.g., Baker et al., 2018). Particle acceleration, as well as losses from pitch angle scattering into the atmosphere, can be intensified during geomagnetic storm periods in space weather events (Reeves et al., 2003). The characteristic timescale of single precipitation loss events reflects the duration of the causative source mechanism, and generally falls in the subsecond to hours range.
The expected effect of electron precipitation on SR depends in part on the characteristic penetration depth of the electrons into the atmosphere (Sátori et al., 2016) and the main factor determining the penetration depth of precipitating electrons is their energy (Rees, 1989), as illustrated in Figure 1. Model calculations (Artamonov et al., 2017;Mironova et al., 2019a) show that the upper ("magnetic") boundary of the Earth-ionosphere cavity resonator (near 90-110 km) (Sátori et al., 2005;2016) is strongly affected by precipitating electrons in the 1-30 keV energy range. However, EISCAT incoherent scatter radar measurements show that the electron density can increase dramatically down to ∼70 km during strongly ionizing EEP events (Belova et al., 2005;Miyoshi et al., 2015). This observation suggests that electrons in the energy range of 30-300 keV can influence the upper ("magnetic") boundary of the Earth-ionosphere cavity resonator as well (Figure 1). It is important to note that Figure 1 illustrates only necessary conditions that energetic electrons will lower the local ionospheric height. Sufficient conditions involve the orders-of-magnitude increase in ionizing flux of energetic electrons, as discussed earlier for X-radiation (Williams and Satori, 2007).
The paper by Toledo-Redondo et al. (2012) is a pioneering work on long-term quasi-global variations in the night-time height of the Earth-ionosphere cavity. Based on the extraction of the waveguide cutoff frequency in DEMETER satellite observations, the authors showed that the effective height of the waveguide was anti-correlated with the solar activity between 2006 and 2009.
A collection of simultaneously recording ELF stations is needed to characterize the deformation of the Earth-ionosphere cavity. In this paper our first objective is to extend the published set of long-term SR intensity records by presenting new results from the high latitude (|lat| > 60°) stations: Hornsund (Svalbard), Syowa (Antarctica) and Maitri (Antarctica), from the mid-high latitude (50°< |lat| < 60°) stations: Eskdalemuir (United Kingdom) and Belsk (Poland), and from the low latitude (|lat| < 30°) station: Shillong (India), as well as an extended set of measurements from the high latitude Ukranian stations in the Arctic (Sousy) and Antarctic (Vernadsky). SR observations are compared with long-term fluxes of precipitating 30-300 keV electrons provided by the National Oceanic and Atmospheric Administration (NOAA) Polar Operational Environmental Satellites (POES) and with electron precipitation events identified in Super Dual Auroral Radar Network (SuperDARN) measurements in Antarctica (Bland et al., 2019). An extended set of the global waveguide's effective height, as observed in the survey electric field ELF-VLF spectra recorded by the French DEMETER satellite (Toledo-Redondo et al., 2012), is also presented and analyzed.

SR Data
In this paper we present long-term induction coil measurements from eight different SR stations which observe the variation of the horizontal magnetic field in the lowest part of the ELF band and analyze monthly average magnetic intensities of the first SR mode. Figure 2 shows the location of the SR stations as well as the great circle paths (GCPs) of the wave propagation directions in which the coils are most sensitive. A common FIGURE 1 | Ionization rates versus altitude for monoenergetic electron flux based on the models described in Artamonov et al., 2017 andMironova et al., 2019a. Frontiers in Earth Science | www.frontiersin.org August 2021 | Volume 9 | Article 689127 map is shown for the Hornsund and Sousy stations in Figure 2 as these two stations are in close proximity (within 150 km, both at Svalbard) compared to the wavelength of ELF waves (of the order of 10 Mm). Induction coils are usually aligned with the local geographical meridian and perpendicular to it, except at the Syowa station where they are oriented along the geomagnetic north-south and east-west directions ( Figure 2F).
The Institute of Geophysics of the Polish Academy of Sciences established SR measurements at the Polish Polar Station Hornsund at Svalbard (HRN; 77.0°N, 15.6°E; L 14.0; Figure 2A) and at the Central Geophysical Observatory in Belsk, Poland (BEL; 51.8°N, 20.8°E; L 2.2; Figure 2C) in 2004 and 2005, respectively Neska et al., 2019). From the HRN station measured data were processed from September 2004 to October 2020. The main data gaps within this time period are: 06/200606/ -11/200706/ , 07-11/200806/ , 03/200906/ , 09/200906/ , 07-08/201106/ , 09/201506/ -03/201606/ and 01-09/2017  Raw time series measurements were bandpass-filtered and 10min estimates of average power spectral density (PSD) with a frequency resolution of ∼0.2 Hz were generated by applying Welch's method (Welch, 1967). Welch's method estimates the PSD by dividing the data into overlapping segments, determining the PSD of each segment and averaging them. In order to minimize the aliasing effect of regional lightning activity and exceptionally intense lightning strokes generating Q-bursts (Guha et al., 2017) the PSD of data segments containing spikes greater than 40 pT in absolute value were omitted before averaging. From the obtained 10-min average PSD the intensity of the first SR mode was determined by means of the weighted average method (Nickolaenko et al., 2015). The peak frequency of the first mode (f 1 ) was calculated from the discrete PSD as: and the PSD of the nearest frequency value to this peak frequency was assigned as the intensity of the first SR mode. SR measurements can easily be affected (contaminated) by different local noise processes (with natural and artificial origin), such as local lightning, wind, human activity, etc. (see e.g., Tatsis et al., 2021;Tritakis et al., 2021). In order to further improve the quality of the results a manual data sanitization step has been carried out before calculating the monthly averaged intensities. Our method relies on the evidence that the daily variation of SR intensity is usually highly similar within a month (Sátori, 1996). Therefore the intensity of the first SR mode was plotted as a function of UT time for all the days within the processed month and days with unusual daily intensity variation were excluded when calculating the monthly average. Frontiers in Earth Science | www.frontiersin.org August 2021 | Volume 9 | Article 689127 SR measurements are carried out by the British Geological Survey at Eskdalemuir Observatory (ESK; 55.3°N, 3.2°W; L 2.7; Figure 2B) near the Scottish Borders of the United Kingdom since September 2012 (Beggan and Musur, 2018;Musur and Beggan, 2019). For the present study processed data are available from June 2012 to December 2020 with a few gaps in the record. Details about the station and data processing can be found in Musur and Beggan (2019).
The Institute of Radio Astronomy (IRA) of the National Academy of Sciences of Ukraine established observations at two high latitude SR stations in the Arctic and in the Antarctic, respectively. The Arctic station is located at the Sousy facility of the Tromso Geophysical Observatory (the Arctic University of Norway, UiT) at Svalbard (SOU; 78.1°N, 16.0°E; L 16.1; Figure 2A) and began operation in September 2013 while the Antarctic station is located in the Western Antarctic at the Ukrainian station "Akademik Vernadsky" (VRN; 65.25°S, 64.25°W; L 2.6; Figure 2D) and began operation in March 2002. From the SOU station data are available from September 2013 to August 2020 with data gaps in 04/2014, 05/2014, 07/2016 while from the VRN station data are available from March 2002 to April 2020 with one long gap between 09/2009 and 03/2010. For more information about the stations and data processing we refer to the paper of Koloskov et al. (2020). Here we note only that data segments with distorted SR spectra were removed from the VRN and SOU databases before averaging. From the VRN station an extended dataset of the spectral intensity at 11 Hz (between the first and second resonance peaks) is also presented (Koloskov et al., 2020). Later in this study we show that a substantial portion of long-term SR intensity variations are related to changes in the quality factor (Q-factor) of the Earth-ionosphere cavity, which characterizes the rate of energy dissipation at a given resonance mode (Madden and Thompson, 1965;Kulak et al., 2003a). The investigation of spectral intensity between the first and second resonance peaks (near 11 Hz) is motivated by this observation as it is expected to have only a minor dependence on the Q-factor of the cavity. We note that recording at the frequency of about 11 Hz was originally proposed and implemented by Fraser-Smith et al. (1991a), Fraser-Smith et al. (1991b. The Indian Antarctic station Maitri (MAI; 70.8°S, 11.7°E; L 5.0; Figure 2E) is located in the Eastern Antarctic and SR data from the station are available from March 2010 to December 2019 with data gaps in 04/2018, 11-12/2018 (H EW ) and 11/2019. It is important to note that at the Maitri station the magnetic coils were oriented initially along the geomagnetic north-south and east-west directions. However, by the end of 2012 they were reoriented along the geographic main directions. In order to work with a homogeneous dataset and to be consistent with the other SR stations, SR intensity values from Maitri are shown only from January 2013 in the present study. For more information about the station we refer to the paper by Manu et al. (2015). Low latitude SR measurements are available from the Shillong station (SHI; 25.92°N, 91.88°E; L 1.15; Figure 2G) in India (Rawat et al., 2012) from January 2008 to December 2016 with data gaps in 07/2008, 05/2014 and 01-04/2015. Raw SR measurements from SHI were processed with the same technique as described for the HRN and BEL stations.
The induction coil magnetometers at the Japanese permanent research station Syowa (SYO; 69.0°S, 39.6°E; L 6.3; Figure 2F) were set up in February 2000 and processed data are available for the present study from January 2006 to December 2015 with data gaps in 03/2008, 07/2008, 06-07/2014 and 06-08/2015. At the SYO station the coils are aligned with the geomagnetic northsouth and east-west directions (D −48.5°in the year of installation). For more information about the station we refer to the paper by Sato and Fukunishi (2005). For the SYO data the same manual data sanitization process as for the HRN and BEL data has been applied.

Satellite Measurements
The flux of precipitating E > 30 keV, >100 keV, >300 keV electrons is provided by the T0 telescope of the Medium Energy Proton and Electron Detector (MEPED) instrument (as part of the Space Environment Monitor 2 instrument package) onboard the NOAA-15 satellite . Electrons that are detected by the T0 telescope (pointing outward along the local zenith) on the quasi-polar (98.70°), Sunsynchronous, low earth orbit (LEO) (∼807 km) of the NOAA-15 satellite are in the atmospheric loss cone and will enter the atmosphere . However, it is to be noted that the T0 telescope detects particles only near the center of the atmospheric loss cone and therefore underestimates the total flux of precipitating electrons van de Kamp et al., 2016;Mironova et al., 2019b).
To demonstrate the typical geographical configuration of the propagation paths (great circle paths) corresponding to the different SR stations and the electron precipitation-affected areas the distribution of precipitated 30-300 keV electrons (as a sum of the measured flux in the three corresponding channels of the T0 telescope) in 1 × 1 degree spatial resolution in 2005 (yielding the highest precipitation fluxes in the time interval investigated in this study) are displayed in the background in the subplots of Figure 2. As can be seen on the maps three main areas are affected primarily by electron precipitation: two high latitude, zonally continuous stripes in the Northern and Southern Hemispheres which we will refer to hereafter as precipitation belts, and the South Atlantic Magnetic Anomaly in South America. Five of the investigated SR stations (HRN, SOU, VRN, MAI, SYO) lie very close (within 1 Mm and << the wavelength of the ELF waves) to at least one of these areas. Furthermore, some propagation paths cross the South Atlantic Magnetic Anomaly. The great circle path corresponding to the H NS component at the SHI station needs to be highlighted as well which does not cross the precipitation belts at all. Furthermore, the South Atlantic Magnetic Anomaly is also very far from this measuring site. Therefore, we do not expect EEPrelated local SR intensity changes in this H NS component.

SuperDARN Measurements
A list of EEP events derived from the Syowa East SuperDARN radar was used to identify individual EEP-related anomalies in Frontiers in Earth Science | www.frontiersin.org August 2021 | Volume 9 | Article 689127 our SR records in June 2011. Although SuperDARN radars are designed primarily for detecting plasma structures in the E and F region ionospheres (Greenwald et al., 1995;Chisham et al., 2007), they can also be used to estimate HF radio wave attenuation in the ionospheric D region using an approach similar to riometry (Bland et al., 2018). In this study we use the EEP event list for the Syowa East SuperDARN radar described in Bland et al. (2019). This radar is located in Antarctica at 69.00°S, 39.58°E (geographic), so all of the June 2011 events were observed under night-time/twilight conditions.

Sunspot Number and AE Index
The traditional sunspot number is used to characterize the magnitude and length of the solar cycles and to identify the year(s) of the maximum and minimum solar activities. The auroral electrojet (AE) index is a measure of geomagnetic activity in the auroral zone (Davis and Sugiura, 1966) and is often applied as an indicator for EEP activity (e.g., Lam et al., 2010). These datasets are available from the OMNIWeb database (https://omniweb.gsfc.nasa.gov/form/dx1.html).

X-Ray Observations
The 0.1-0.8 nm X-ray measurements of the Space Environment Monitor (SEM) instrument onboard the GOES 10 and 15 satellites is used to characterize solar radiation-related changes in the Earth-ionosphere cavity (Sátori et al., 2005;Sátori et al., 2016). Monthly average flux values were determined from the 1min and 5-min resolution datafiles available from NOAA's corresponding database (https://satdat.ngdc.noaa.gov/sem/goes/ data/avg/).

Satellite-Based Observations of Waveguide Cutoff
Long-term variation of the upper boundary region of the Earth-ionosphere cavity resonator is also investigated in this study by using satellite-based wave recordings. The lower ionosphere forms the highly variable reflecting upper surface of the Earth-ionosphere waveguide (EIWG). Since far from the source, lightning impulses propagate in a guided manner, the effective height of the EIWG (reflecting altitudes) can be observed as suppressions in wave field strength at higher ELF, lower VLF frequencies around the cutoffs of distinct wave modes (see modeling details in Cummer, 2000;Ferencz et al., 2007). This effect is best recognized in case of the 1st guided mode at approximately 1.6 kHz, exhibiting often a lower intensity band in the spectrum of single impulses as well as in time averaged dynamic power spectra of the cavity background field (Toledo- Redondo et al., 2012). Due to the excitation of the lower ionosphere by waves in the EIWG with the above described character, the wide-band ELF-VLF wave records of LEO satellites also exhibit this specific spectral pattern. The electric field data of the French DEMETER satellite (Berthelier et al., 2006) with a polar LEO orbit has been utilized here to determine guided cutoffs, and thus to follow the change of the cavity's shape during approximately a quarter solar cycle period. Six years of night-time ELF-VLF spectra (whole years of 2005-2010) recorded at topside altitudes (660-710 km) between 65°geomagnetic latitudes have been analyzed. This uniquely rich, continuous ELF-VLF wave database is composed of power spectrum values under 20 kHz with 19.53 Hz resolution in frequencies and 2.048 s resolution in time (survey mode). Occurrences of wave power minima around the supposed cutoffs and their parameters (minimum frequency, sharpness, reliability factor) has been determined in consecutive spectrum vectors by applying a quadratic fit in the 1.4-1.95 kHz frequency range. The settings of the identifying algorithm have been fine-tuned by using a training set of several hundred of half orbits selected randomly from the whole recording base. In order to suppress noise and still maintain good spatial and temporal resolution, an averaged spectrum of eleven consecutive measurements was found optimal as input data. This step reliably finds the minima (if present) and yields approximately 1.5°spatial resolution in latitude. The effective height corresponding to the cutoff frequency is then determined by applying the condition for the transverse resonance (see e.g., Bliokh et al., 1977) of the waveguide (Toledo-Redondo et al., 2012): We present annual averages of the effective height for different latitudinal domains. Raw data consisting of satellite coordinates, times at the mid positions of the averaged recording set and the parameters of the identified effect were stored for later analysis. The investigated domains consist of various geomagnetic latitudes: low-(|mlat| < 20°), mid-(20°≤ |mlat| < 50°) and mid-high (50°≤ |mlat| ≤ 65°). Data points corresponding to the South Atlantic Magnetic Anomaly (defined as −50°≤ |lat| ≤ 0°and −90°≤ lon ≤ 40°) were excluded when calculating the yearly averages.

Results Based on SR Measurements
Multi-station Schumann resonance observations are used in this study to investigate long-term variations in SR intensity, selected year-to-year variations, and short-term variations attributable to single EEP events that affect both magnetic hemispheres. These topics are addressed in three separate sections below.  Figure 3A shows the monthly average values of the sunspot number, the AE index and the X-ray flux. The EEP flux of 30-300 keV electrons in daily (black) and monthly (red) time resolution are shown in Figure 3B. The EEP flux in daily time resolution is included to demonstrate its highly variable nature. The precipitation fluxes can vary by 3 orders of magnitude within a few days during geomagnetic storms . As can be seen in Figure 3A, the X-ray flux follows the time variation of the sunspot number in contrast to the AE index which has a distinct phase delay with respect to the solar maxima after 2002 and 2014/2015, and also shows extraordinary behavior in 2012 and in 2013 preceding the last solar cycle maximum in 2014. A good correlation is evident between the time series of the AE index and the 30-300 keV EEP fluxes as shown in Figures  3A,B. We note that a moderate annual variation is present in both parameters with Northern Hemisphere (NH) summer maxima as noted in earlier works (e.g., Suvorova, 2017;Lockwood et al., 2020). The SR intensity records in Figures 3C-J show a complex behavior. On one hand a seasonal variation is present in the data with the highest intensity values usually observed in NH summer months. This general behavior is in accordance with the NH summer annual maxima in the magnitude of global lightning activity (Sátori, 1996;Christian et al., 2003). We also note the presence of semi-annual variations (with spring and autumn peaks) at some SR stations (Sátori and Zieger, 1996). The absolute level of SR intensity may show up to an order of magnitude difference station-by-station. We attribute this observation primarily to the different source-observer geometries but it could also hint to problems with the absolute calibration of the stations. Calibration issues are the most likely explanation for the eye-catching differences between the SR intensity records of the SOU and HRN stations which should be almost the same based on the close proximity of the stations.

Long-Term Variations
On the other hand, the solar cycle modulation is evident at all SR stations and is in phase with the solar cycle (as indicated by the sunspot number). The solar cycle modulation of SR intensity is usually more pronounced in one of the magnetic field  (2005,2015,2016) or in the vicinity of the solar cycle maxima (2012,2013,2014), are in agreement with enhanced levels of EEP as indicated by the AE index and the 30-300 keV EEP fluxes. As the solar cycle modulation of SR intensity is clearly present at the two midhigh latitude stations (ESK and BEL) and at the one low latitude (SHI) station as well, the effect is not confined to high latitudes.
The magnitude of the annual variation in SR intensity tends to follow the solar cycle as well. This behavior is in keeping with the suggestion of the solar cycle modification of cavity height. Higher magnitudes can be identified in the years near the solar maxima and smaller ones near the solar minimum (see Figures 3C-J). This observation had been noted earlier for the total horizontal magnetic field by Koloskov et al. (2020). To take a closer look at this phenomenon we have removed the yearly averages from the VRN records for each (H NS , H EW ) field component and plotted the residuals (ΔI) in Figure 4. The solar cycle modulation of the annual SR intensity variation is different for the two field When calculating the yearly averages some missing values were determined with interpolation using the SR intensity values of the neighboring months (in case of 1 missing month) or based on the same months measured in the neighboring years (in case of >1 adjacent missing months). If neither of these two approaches was feasible for any year, that year was omitted. The long-term records of spectral intensity at 11 Hz are also displayed for each (H NS , H EW ) field component from the VRN station in the same format.
As we stated earlier in this section, solar radiation-related changes in the Earth-ionosphere cavity are expected to follow the solar cycle (as indicated by the sunspot number) while EEPrelated changes are expected to deviate from that cycle and to show a common behavior with the AE index. Next we investigate the normalized yearly average intensities displayed in Figure 5 in this aspect. The relative SR intensity values corresponding to the first resonance mode generally follow the different phases of the solar cycle. However, this statement is rather true for the field component exhibiting smaller long-term variability (H EW at

Year-to-Year Comparisons
As confirmed by Figure 5, it is not always the solar cycle maximum that yields the highest SR intensity values. This observation is in accordance with electron precipitation fluxes presented in Figure 3B. In Figure 6

Short-Term Variations Synchronous in Two Hemispheres
The well-organized motions of energetic electrons along the magnetic field lines suggest that EEP events and their associated cavity deformation may often be quasi-synchronous in the Northern and Southern Hemispheres. Figure 7 shows SR intensity records for two selected days from June 2011 (14 June and 23 June) together with independently identified EEP events based on Syowa East SuperDARN radar measurements from the event collection of Bland et al. (2019). To identify non-standard (not lightning-related) variations in the diurnal SR intensity records, quiet day diurnal SR intensity curves were determined for each field component of each station for June 2011. Quiet days were selected based on the AE index as days with AE values smaller than 400 nT. Figure 7 contains the average diurnal variation of these quiet days and in addition the corresponding standard deviations. In order to investigate whether EEP events are able to modify the global Q-factor of the Earth-ionosphere cavity resonator we applied the Lorentzian fitting technique (see e.g., Mushtak and Williams, 2002) on SR spectra from the Hornsund station. In Figure 7 the extracted Q-factors are displayed in the same format as the SR intensity records.   On June 14, 2011 the strongest EEP event within the day can be observed between 09:20 and 14:30 UT at Syowa with radio wave absorption values greater than 2 dB ( Figure 7A). The event starts with a substantial increase of the AE index. SR intensity clearly increases during the event not only in the H NS component at Syowa (so near the radar) but in the Northern Hemisphere under daytime conditions at the Hornsund (HRN) station as well where increased Q-factor values (by ∼10-30%) can be observed during the event. The relative increase of SR intensity is as large as 50-100% at Syowa during this event. On June 23, 2011 several strong (absorption >3 dB) EEP events occured between 00:00 and 13:00 UT and one other event can be observed between 21:00 and 23:30 UT. The imprint of these events can be identified in both field components at SYO, in the H NS component at HRN (with increased Q-factor values by ∼10-30%) as well as in the H EW component at Vernadsky (VRN), even though the VRN station lies in a considerably different longitudinal sector than Syowa. This documented synchronous behavior provides evidence for the large spatial extents of the EEP events (∼4 h magnetic local time extent).

Results Based on DEMETER Measurements
The effective height of the Earth-ionosphere waveguide as inferred from DEMETER satellite measurements of the waveguide cutoff frequency at night-time (see Section Satellite-Based Observations of Waveguide Cutoff for details) can be regarded as a quantity providing an independent view on the deformation of the Earth-ionosphere waveguide. Figure 8 shows the yearly average values of the effective height from 2005 to 2010 corresponding to low-, mid-and mid-high geomagnetic latitudes.
In general, the Earth-ionosphere waveguide height seems to be the lowest at mid-high latitudes while it is the highest at midlatitudes with average height values of 90.4 and 91.8 km, respectively. The effective height is also smaller at low latitudes (90.7 km) as compared to midlatitudes (91.8 km). Regarding the long-term trend in the data, the solar cycle variation is evident in all the investigated domains. Table 3 summarizes the relative height differences between 2005 and the solar minimum year of 2009 for low-, mid-and mid-high geomagnetic latitudes. The relative height difference is the greatest at mid-high latitudes (which represent the highest observed latitudes in this case) (2.2%) and the smallest at low latitudes (1.2%).

DISCUSSION
In this section, we suggest a general interpretation for the results obtained, before attempting to estimate the height changes of the Earth-ionosphere cavity associated with EEP. Finally, we describe the most important characteristics of our SR-and DEMETERbased results.

Solar X-Rays and EEP: The Two Main Effects to be Considered
In order to interpret correctly the obtained results it is first important to state that although SR intensity variations are generally connected with changes in the distribution and/or intensity of global lightning activity, satellite observations of lightning activity do not show any significant solar cycle variation on the global scale (Christian et al., 2003;Williams et al., 2014). This is in line with the observation that the total FIGURE 8 | Yearly averages of the Earth-ionosphere waveguide's effective height based on DEMETER observations. A similar solar cycle modulation is evident in each of three latitude ranges. energy output from the Sun varies on the order of 0.1% throughout the solar cycle (Fröhlich and Lean, 1998). Therefore, the solar cycle modulation of SR intensity needs to be accounted for by long-term changes in the Earth-ionosphere cavity, i.e., by changes in the propagation conditions of ELF waves. As the solar cycle variations in Galactic Cosmic Rays (GCR) is considered to have only a negligible effect (<1%) on the intensity SR (Sátori et al., 2005) hereafter we discuss the possible effects related to the Sun. Two fundamental mechanisms can be distinguished for the Sun-related impact on the near-Earth environment: one based on photons and the other on charged particles. The majority of solar photons arrive in the Earth's equatorial zone while charged particles precipitate in the upper atmosphere in the high latitude precipitation belts and in the South Atlantic Magnetic Anomaly (Figure 2). The long-term behavior of these two effects is markedly different on the solar cycle timescale. While the flux of solar photons generally follows the solar cycle, the precipitation of charged particles does not (Asikainen and Ruopsa, 2016) ( Figure 3). We documented both of these behaviors in longterm SR intensity records of different field components at different stations ( Figure 5). Therefore, we suggest that a combined effect of photons and charged particles needs to be considered to account for all our SR-based observations (as illustrated in Figure 9).
The long-term solar cycle modulation of the Earth-ionosphere cavity by solar X-rays had been first documented by Sátori et al. (2005) based on long-term SR frequency records. Both SR frequencies and Q-factors observed at different stations exhibited maxima during the two solar maxima in 1992 and 2002 and minima in 1996. The observed variations were interpreted on the basis of a uniform cavity and the presence of two characteristic ionospheric layers (Greifinger and Greifinger, 1978;Mushtak and Williams, 2002). The height and scale height of the lower (electric) layer were considered to be invariant over the solar cycle while for the upper (magnetic) layer the authors inferred a height decrease of 4.7 km and a scale height decrease from 6.9 to 4.2 km for the first SR mode from solar maximum to solar minimum. These changes result in an increase of the resonator's Q-factor by 12-15% (Sátori et al., 2005). The intensity of a damped simple harmonic oscillator is known to be proportional to the square of its Q-factor (Feynman et al., 1963), in agreement with the theoretical description of SR (Nelson, 1967, Eqs 2-24;Williams et al., 2006). It follows that a 12-15% change in the Q-factor of the Earth-ionosphere cavity can result in a 25-30% increase of the first SR mode's intensity. This estimate is in agreement with the relative SR intensity values corresponding to the H NS component at the SHI and VRN stations ( Figure 5). As the Q-factor is a general property of the Earth-ionosphere cavity this effect should contribute to the long-term SR intensity variation at all SR stations around the globe.
The phenomenon involving charged particles deposited into the atmosphere through various heliophysical and geomagnetic processes is called energetic particle precipitation (EPP). EPP covers a wide variety of charged particles but protons and electrons are the two main constituents to be considered with significantly high fluxes (Mironova et al., 2015). As it has been stated in earlier works (Williams et al., 2014;Sátori et al., 2016) proton events are too sporadic in time to yield the observed persistent long-term modulation of SR intensity. Therefore, we consider the precipitation of electrons hereafter which is usually referred to as energetic electron precipitation (EEP).
The expected effect of EEP on SR strongly depends on the characteristic penetration depth of the electrons into the atmosphere. SR are primarily responsive to two specific altitude regions where the ionospheric dissipation is predominant (Sátori et al., 2016). These are usually referred to as the electric (capacitive) and magnetic (inductive) heights (Madden and Thompson, 1965;Greifinger and Greifinger, 1978) to be found at around 50-60 km and 90-110 km altitudes, respectively (Mushtak and Williams, 2002;Greifinger FIGURE 9 | Illustration of the deformation of the Earth-ionosphere cavity at solar minimum (A) and solar maximum (B). At solar maximum X-rays and EEP reduce the cavity height more significantly (dominantly over low-and high latitudes, respectively) which results in enhanced SR intensity values beneath the deformed regions.
Frontiers in Earth Science | www.frontiersin.org August 2021 | Volume 9 | Article 689127 et al. Kulak and Mlynarczyk, 2013). The main factor determining the penetration depth of precipitating electrons is their energy (Figure 1). Model calculations show that the magnetic height can be affected by electrons in the 1-30 keV energy range while the electric height only by electrons with energy above 300 keV (Rees, 1989;Sátori et al., 2016;Artamonov et al., 2017;Mironova et al., 2019a). However, EISCAT incoherent scatter radar measurements show that in connection with EEP the isopycnic surfaces of electron density can occasionally descend by some tens of kilometers in the high latitude D-and E-region (Belova et al., 2005;Miyoshi et al., 2015). To account for such a large decrease in the SR magnetic height the effect of 30-300 keV electrons needs to be considered as well.
Generally, it is true that the lower the energy of electrons the higher the precipitated flux and the more continuous the precipitation in time . Therefore, to account for the long-term modulation of SR intensity, electrons in the 1-300 keV energy range (affecting the magnetic height of SR) are the best candidates. This energy range covers both auroral and radiation belt electrons (Mironova et al., 2015). EEP is expected to affect the height of the Earth-ionosphere cavity and the measured magnetic field locally where precipitation occurs consistent with the conservation of energy flux (Sentman and Fraser, 1991;Williams et al., 2014;Sátori et al., 2016). A permanent local decrease in the magnetic height is expected to modulate the amplitudes of all natural variations in SR intensity on all timescales as shown for the magnitudes of the annual variation in Figure 4. We also showed that EEP events are able to modify the Q-factor of the Earth-ionosphere cavity (Figure 7) which is in line with our observation that they can affect both hemispheres as well as large longitudinal regions. This effect may result in globally observable changes in SR intensity but this aspect needs to be investigated further in the future. It is possible that the long-term 12-15% change in the Q-factor of the Earth-ionosphere cavity documented by Sátori et al. (2005) may involve EEP-related changes as well.
Two important pieces of evidence were shown in this paper for an EEP-related deformation of the Earth-ionosphere cavity. First, exceptional increases of SR intensities were documented in certain years, as in 2005 and 2012 ( Figure 5), which are not solar maximum years, but are characterized by enhanced electron fluxes ( Figure 6). Based on satellite observations Asikainen and Ruopsa (2016) reported on enhanced electron precipitation at >30 and >100 keV fluxes related to Coronal Mass Ejections (CMEs) in 2005 and 2012 which is in agreement with our findings. Second, direct observation of EEP-related SR intensity enhancements were presented at widely separated SR stations on two selected days in June 2011 (Figure 7).

Estimating the Long-Term Ionospheric Height Changes Associated With EEP
In this section we estimate the long-term average ionospheric height changes associated with EEP for the H EW component at the VRN station yielding the highest intensity variation on the solar cycle timescale (40-60%) ( Table 4). If we accept that 25-30% of the observed SR intensity enhancement is attributable to a global change of the cavity's Q-factor (independently of whether this change is related solely to X-rays or to EEP as well) the remaining 15-30% in years with the largest enhancement can be considered to be related to the local waveguide deformation caused by EEP.
We estimate the EEP-related height changes of the waveguide from the intensity at 11 Hz (the first minimum between the first and second resonance peaks) which can be considered to be mostly unaffected by the Q-factor of the cavity. Greifinger et al. (2005) have made predictions for changes in the SR amplitude as a function of the ionospheric height above the source and above the observer in the form of symbolic equations. According to these equations the measured magnetic intensity at a given frequency is proportional to the reciprocal value of the squared local magnetic height. It follows that the 20-25% intensity enhancements documented at 11 Hz in 2005 and 2012 (which are in a good agreement with the resonance mode-based estimation of 15-30%) correspond to an average magnetic height decrease of ∼10% at VRN. This average decrease in the magnetic height is considerably smaller than the decrease of isopycnic surfaces during an EEP event as observed by the EISCAT radar, with enhanced ionization down to ∼70 km (e.g., Miyoshi et al., 2015). Although EEP events are very common especially near solar maximum and can result in SR intensity enhancements in order of 50-100% (as shown in Figure 7), they are not permanent and each event is often confined to a limited longitudinal range (Clilverd et al., 2007;Bland et al., 2021). Therefore, the smaller percentage average increase of SR intensity and the inferred decrease in the magnetic height is reasonable on the interannual timescale.

Interpretation of SR-Based Results
As can be seen in Figure 5 and Table 1 the solar cycle modulation of SR intensity can be different station-to-station depending on the station location and on the wave propagation path corresponding to the measured field component. A deformation in the Earth-ionosphere cavity is expected to affect SR intensity locally where the deformation occurs, whereas a change in its Q-factor is expected to affect SR intensity (I) globally, for all propagation paths that cross the affected area(s). As we stated earlier, we aimed to describe these effects by two simple formulas: and where Q and H m denote the Q-factor of the cavity and its local magnetic height, respectively. The solar cycle modulation of the H NS field component at the SHI and VRN stations can be well described by the X-ray-related global variation of the cavity's Q-factor and the H EW field component at the VRN station by the combined effect of the Q-factor variation and an EEP-related height decrease of ∼10%. At different SR stations the wave propagation directions corresponding to the perpendicular magnetic coils are differently affected by EEP. It seems to be generally true that the longer the propagation path crosses the precipitation-affected area(s) the larger the solar cycle modulation of SR intensity. It is usually the H NS coil, sensitive to the east-west propagation direction, which yields propagation paths with greater extents beneath the precipitation belt. This is the case for the high latitude SOU, HRN,MAI and SYO stations (see Figures 2A,E,F) where the solar cycle modulation is indeed larger in the H NS component as compared to the H EW component (Figures 3, 5). The H NS coil is also more responsive at the two mid-high latitude SR stations ESK and BEL in spite the fact that the GCPs corresponding to the H NS component at these two stations do not cross the northern precipitation belt according to Figures 2B,C. This result may indicate that the precipitation belts are more extended towards lower latitudes as shown by the POES data. Towards lower latitudes the flux of precipitating 30-300 keV electrons decreases and POES is known to underestimate flux values below ∼10 4 -10 5 electrons/cm 2 /sec/ster (van de Kamp et al., 2016;van de Kamp et al., 2018).
The magnetic height of SR is affected by electrons with energies lower than 30 keV (the lowest channel of POES satellites) and these electrons may form wider precipitation belts. At the VRN and SHI stations the SR response is more pronounced in the H EW coil, sensitive to the north-south propagation direction. At SHI the H EW is the one field component with a GCP crossing the precipitation belts ( Figure 2G) while at VRN the propagation path corresponding to the H EW component meets not only the southern precipitation belt but the South Atlantic Magnetic Anomaly as well ( Figure 2B). Besides the dominant response of the H EW component as compared to the H NS component this fact may account for the observation that the H EW component at the VRN station yields the highest response among all the investigated records. Another factor that should be important for this observation is that the longitudinal sector including America is shown to be affected by harder electron precipitation spectrum than the other parts of the globe which can be attributed to the tilted-offset dipole magnetic field of Earth (Clilverd et al., 2007). We note that 2012 shows the highest relative SR intensity values at the HRN and VRN stations ( Figure 5; Table 1) despite the observation that POES measurements show larger EEP fluxes in 2005 as compared to 2012 ( Figure 6). We can explain this observation by the fact that in 2012 both X-rays and EEP show enhanced values while 2005 is in the declining phase of the solar cycle characterized by a lower X-ray flux ( Figures 3A,B).
As we noted earlier in Section Results EEP is known to have seasonal dependence (Clilverd et al., 2007;Suvorova, 2017;Bland et al., 2019), i.e., June/July maxima in its occurrence rate can be identified in both hemispheres. We note that rocket measurements at polar stations also showed maximum electron density in June/July at the 80 km height of the ionospheric D-region in both hemispheres (Danilov and Vanina, 1999). We suggest that this seasonal dependence may contribute to the observed modulation of the magnitude of the annual SR intensity variation (Figure 4).

Interpretation of DEMETER-Based Results
DEMETER observations showed that the effective height of the Earth-ionosphere waveguide is generally the highest at mid geomagnetic latitudes, lower at low geomagnetic latitudes and the very lowest at the highest observed geomagnetic latitudes (in proximity to the precipitation belts) ( Table 3) which is in a general agreement with our SR-based findings that the Earth-ionosphere cavity is deformed by solar X-rays (dominantly over lower latitudes) and by EEP (dominantly over higher latitudes) (as illustrated in Figure 9). The deformation effect of EEP has a larger variability on the solar cycle timescale which results in a larger height reduction at higher latitudes as shown by the effective heights. This finding is in keeping with the high latitude location of the precipitation belts. Although the documented 2.2% height decrease between 2005 and 2009 corresponding to mid-high latitudes is considerably smaller than the ∼10% value inferred for the VRN station, one needs to keep in mind that 1) DEMETER data are restricted to geomagnetic latitudes lower than 65°, 2) the South Atlantic Magnetic Anomaly (which is in close vicinity to VRN) is excluded from the dataset, 3) electromagnetic waves with higher ELF, lower VLF frequencies are responsive at different altitude regions than the magnetic height of SR and 4) this is an integrated value for all longitudes while the longitudinal sector including America is shown to be affected by harder electron precipitation.
It is to be noted as well that the inferred effective height values were measured exclusively on the nightside by DEMETER at 22:30 LT where solar X-rays are absent. We can explain this seeming contradiction by assuming that lightning impulses may propagate in the dayside hemisphere of the EIWG before escaping the waveguide on the nightside (but close to the local sunset).

CONCLUSIONS
In the present study long-term records of the first SR magnetic mode intensity from eight different stations were compared with the fluxes of precipitating medium energy electrons and with independently identified EEP events to confirm that the longterm solar cycle modulation of SR intensity (documented originally at high latitude stations) is caused by the EEPrelated local deformation of the Earth-ionosphere cavity. Although our results generally confirmed the presence of EEPrelated deformations in the Earth-ionosphere cavity, we also found that the EEP-effect alone cannot account for all our SRbased observations and that the combined effect of solar X-rays and EEP needs to be considered. We have identified four distinct factors that can play important role in shaping long-term SR intensity records: 1) X-ray related deformations of the cavity, 2) X-ray related changes in the Q-factor of the cavity, 3) EEP-related deformations of the cavity and 4) EEP-related changes in the Q-factor of the cavity. The exact interplay of these factors depends on the location of the SR station and the orientation of the magnetic coils, i.e., on the wave propagation path corresponding to the actual measurement. The deformation of the cavity is expected to affect SR intensity locally beneath the deformed area(s) whereas Q-factor variations are expected to introduce globally observable changes in SR intensity.
In connection to EEP we showed that there is a very good phase agreement between certain SR intensity records and the long-term variation of 30-300 keV EEP fluxes as provided by the POES NOAA-15 satellite which is most conspicuous in years not exactly coincident with the solar maximum, as in for example 2005 and 2012 when both the EEP flux and SR intensity showed exceptional values. We documented SR intensity and Q-factor enhancements connected to individual EEP events on the daily timescale in June 2011 by comparing the corresponding records with EEP events independently identified by Bland et al. (2019) in the measurements of the Syowa East SuperDARN radar. We showed that the effect of certain EEP events can be identified in SR intensity records from both the Southern and Northern Hemispheres (under different sunlit conditions), from different longitudinal sectors and the EEP-related relative enhancement of SR intensity can attain values as high as 50-100%. The solar cycle modulation of the magnitude of the annual SR intensity variation was demonstrated and efforts were made to quantify the ionospheric height changes associated with EEP. It has been shown that near the Antarctic station Vernadsky (yielding the largest long-term SR intensity variation) the upper (magnetic) characteristic height of the Earth-ionosphere cavity resonator decreased by ∼10% in the most conspicuous years of EEP activity. This paper is the first one to show that the solar cycle variation of SR intensity is consistently observable all around the globe and gives a qualitative explanation for this observation. Quantitative conclusions were made as well based on simple SR intensity-cavity height, SR intensity-Q-factor relations, which describe the identified effects reasonably well. To achieve a more detailed interpretation of our observationbased results it is necessary to apply a full model of the SR.
We hope that our work will motivate and assist such works in the future.
The effective height of the Earth-ionosphere waveguide as inferred from DEMETER measurements provided an independent view on the long-term deformation of the Earth-ionosphere cavity. We showed that the largest height of the cavity can be found at mid (geomagnetic) latitudes while the waveguide is depressed at low and high (geomagnetic) latitudes which is in a general agreement with our general conclusion that the Earth-ionosphere cavity is deformed by solar X-rays dominantly over lower latitudes and by EEP dominantly over higher latitudes. It has also been demonstrated that effective height of the waveguide varies on the solar cycle timescale at all (geomagnetic) latitudes and this effect is stronger at high (geomagnetic) latitudes where energetic electrons can enter along magnetic field lines.
We suggest that SR measurements may be suitable to collect information and thus to monitor changes in the lower ionosphere associated with EEP which would be highly valuable for investigating space weather processes and their impact on atmospheric chemistry (see e.g., Mironova et al., 2015;Duderstadt et al., 2021). Such measurements could yield quasicontinuous information on EEP in a global sense in contrast to the space and time limitations of satellite and radar observations. Finally, the value of high latitude stations for SR-based research should be emphasized. Due to their remote locations, these SR stations have usually much lower cultural noise contamination than stations at low-and mid-latitudes. And because they are also remote from thunderstorm activity they are not subjected to the interference caused by nearby lightning. These two properties make high latitude stations highly valuable for SR-connected research aims, such as the reconstruction of global lightning activity (see e.g., Bozóki et al., 2019;Prácser et al., 2019) by inversion methods. However, because of the extraterrestrial effects on the SR cavity at higher latitude as documented in this work, it is of vital importance to realize all possible masking effects when analyzing these data-sets in connection with lightning-related climate research.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
TB, GS, and EW led the interpretation of results and the writing of the manuscript. TB processed the SR data and prepared the figures. IM processed the POES data and contributed to the interpretation of EEP-related results. PS and ST-R processed the DEMETER data and contributed to the interpretation of DEMETER-based results. EB provided information on individual EEP events and contributed to the interpretation of the EEP-related results. AK, YY, and OB provided SR data from the Vernadsky and Sousy stations and contributed to the interpretation of SR-related results. MN provided SR data from the Hornsund and Belsk stations. AS and RR provided SR data from the Maitri and Shillong stations. MS provided SR data from the Syowa station. CB provided SR data from the Eskdalemuir station. YL collected and organized long-term SR records in the MIT SuperCloud database. RB applied the Lorentzian fitting technique on SR data.

FUNDING
The work of TB and GS was supported by the COST Action CA15211 and by the National Research, Development, and Innovation Office, Hungary-NKFIH, K115836. The contribution to this work by EW and YL has been supported by both the Grainger Foundation and by the US National Science Foundation (Grant # 6942679). The work of the Ukrainian team (AK, YY, and OB) was partially supported by MES of Ukraine projects 0121U112293 and NAS of Ukraine projects 0121U108110 and 0121U108635. EB is supported by the Research Council of Norway (grant #287427). IM work on processing, interpreting results of EEP data was supported by a grant from the Russian Science Foundation . The database of POES data was prepared by IM in the frame of the SPBU "Ozone Layer and Upper Atmosphere Research Laboratory" supported by the Ministry of Science and Higher Education of the Russian Federation under agreement 075-15-2021-583.