Satellite Measured Ionospheric Magnetic Field Variations over Natural Hazards Sites

: Processes and threats related to natural hazards play an important role in the evolution of the Earth and in human history. The purpose of this study is to investigate magnetic ﬁeld variations measured at low Earth orbit (LEO) altitudes possibly associated with earthquakes, volcanic eruptions, and artiﬁcial outbursts. We focus on two missions with well equipped magnetometer packages, the China Seismo-Electromagnetic Satellite (CSES) and ESA’s three spacecraft Swarm ﬂeet. After a natural hazards survey in the context of this satellites, and consideration of external magnetospheric and solar inﬂuences, together with spacecraft interferences, wavelet analysed spatio-temporal patterns in ionospheric magnetic ﬁeld variations related to atmospheric waves are examined in detail. We provide assessment of the links between speciﬁc lithospheric or near surface sources and ionospheric magnetic ﬁeld measurements. For some of the diverse events the achieved statistical results show a change in the pattern between pre- and post-event periods, we show there is an increase in the ﬂuctuations for the higher frequency (smaller scales) components. Our results are relevant to studies which establish a link between space based magnetic ﬁeld measurements and natural hazards.


Introduction
The broad context of this study is addressing the issues of space based magnetic field measurements in the frame of natural hazards investigations. The impacts on society of, e.g., volcanic eruptions or earthquakes could be severe and systematic action is of vital importance. The purpose and focal point of this work is on the advanced magnetic field observation capabilities of two space missions. These are the China Seismo-Electromagnetic Satellite (CSES; or Zhangheng-1, named after the Chinese polymath), a China National Space Administration (CNSA) and Italian Space Agency (ASI) collaboration, and ESA's Swarm mission, and how a synergy of different instruments can achieve a global space based picture of the mentioned threats. The satellites are equipped with a rich and strong portfolio of ambitious field-and plasma-packages for a wide range of scientific investigations, but the main aim of this work is a characterisation of natural lithospheric events primarily with the magnetic field in the key role, i.e., a single dominant physical parameter. In recent complementary efforts, a broad spectrum of satellite based observations had already an important contribution to the overall characterisation of these phenomena with the potential of a variety of societal benefits and applications [1].

Magnetic Field Measurements at Low Earth Orbit Altitudes
The scalar-and vector-magnetic field measurements are acquired from two low Earth orbit (LEO) missions around 500 km altitude with their respective instruments operating in the ionospheric F-region, CSES and the three Swarm satellites.

CSES
The main objective of CSES is the study of ionospheric responses to lithospheric processes, i.e., electromagnetic-and plasma-phenomena related to natural hazards, in particular earthquakes (EQ) with possibly precursor phenomena. The three axes stabilised spacecraft with eight scientific instruments was launched on 2 February 2018, into a sun synchronous polar orbit with an altitude of ∼507 km and an inclination of 97.4 deg [2]. One of the instrument packages is the high precision magnetometer (HPM) consisting of two fluxgate magnetometer (FGM; FGM1 and FGM2) and a scalar coupled dark state magnetometer (CDSM). The HPM data processing pipeline includes corrections, interference mitigation, synchronisation, and calibration steps for the nominally sampled 60 Hz FGM vectors and 30 Hz CDSM values [3]. The science data cover the latitude range from −65 deg to +65 deg with a temporal resolution of 1 s. The CDSM instrument is a scalar magnetometer based on two-photon spectroscopy and has a noise level of ∼50 pT/sqrt(Hz) at 1 Hz, an accuracy of 0.19 nT (1σ), and a total magnetic field measurement range from 10 3 -10 5 nT [4]. The HPM instruments are parallel aligned, from the satellite surface outwards in the layout FGM1-FGM2-CDSM, and mounted along a 4.7 m nadir pointing boom in order to minimise the spacecraft magnetic disturbances. In-orbit results with a performance analysis have been reported by [5], the CSES data description, and data service are detailed in [6].

Swarm
Swarm is an ESA mission (launched on 22 September 2013) consisting of three identical satellites in near-polar, circular orbits with the goal to study the Earth's magnetic field in all aspects. The implemented orbital constellation are two satellites closely spaced side-by-side (i.e., 1.4 deg longitudinal separation and only a few seconds differential delay in orbit) at an altitude of <470 km (Alpha, Charlie) and the third one (Bravo) a bit higher at <520 km (altitude decrease during life time). The Level-1b magnetic field products are derived from two instruments, the absolute scalar magnetometer (ASM) and the vector field magnetometer (VFM) with 1 Hz data rates for both instruments (low rate) and in addition 50 Hz (high rate) for the VFMs. A detailed description, including nomenclature, ephemeris-and plasma-products together with derived Level-2 products is given by [7].

Magnetic Field Models and Multi-Experiment Observations
The generation of geomagnetic field models including secular variations is an important application of space based measurements to study long term behaviour in terms of month and timescales beyond, e.g., [8]. In case of distinct natural hazards we rely on observations immediately before and after the events. The time frame could be from a few weeks, days, up to hours-in any case as much as possible steady conditions-and also depends on the time of the preparation phase of an, e.g., earthquake. This implies, at first glance, that carrying out and comparing such models are not fully useful, narrow preand post-spatio-temporal high-quality samplings in multi experiment settings and complementary measurements are preferable [9,10]. The 1 Hz magnetic field data from CSES HPM permit the analysis of large scale structures, whereas the Swarm VFM 50 Hz high rate values enable the study of, e.g., fine structured details during boundary crossings in the ionosphere above target objects on the ground.

Lithospheric Atmospheric Ionospheric Coupling
Important scientific objectives are the couplings and energy carriers from solid earth sources up to ionospheric heights, the so-called lithospheric atmospheric ionospheric coupling (LAIC). This is in particular relevant if physically based short-term earthquake forecasting in the context of precursory signals of natural hazards shall be achieved, see [11] and references therein. Joint collaborative efforts in a multi-parameter approach shall facilitate early risk mitigation of this most of the time cascade of events and protect society and infrastructure. Simões et al. (2012) [12] gave a summary of these couplings in terms of electromagnetic waves, refs. [13][14][15] investigated precursory phenomena, the latter with Swarm data in respect of, e.g., precursor time versus main shock magnitude [16]. Marchetti et al. (2020) [17] conducted co-seismic investigations with Swarm magnetic field data, ref. [18] used magnetic field and electron density measurements in the analysis of large earthquakes.

Earth-Ionosphere Waveguide
The Earth-Ionosphere Waveguide extends from the surface (lower boundary) up to the lower ionosphere (upper boundary), which is the D-layer during daytime and the Elayer during night. The wave propagation between narrowband very low frequency (VLF; f = 3-30 kHz) and low frequency (LF; f = 30-300 kHz) transmitter and receiver stations is confined in this cavity, e.g., [19]. In a network configuration with several VLF/LF propagation paths it is possible to sample larger areas, e.g., over Southern Europe [20]. The upper waveguide boundary could be externally driven, e.g., through solar X-ray flares [21], but significant internal drivers can be acoustic waves (AW) and atmospheric gravity waves (AGW) [22]. A thorough and systematic understanding of phenomena up to the lower ionosphere is a basic requirement for investigations of the upper atmosphere (thermosphere) with decisive in-situ measurements. An important task is to systematically quantify uncertainties in the datasets and finally validate and verify the processes in accurate and predictive ways.

Acoustic Waves and Atmospheric Gravity Waves
The approach to delineate irregularities in the upper atmosphere with atmospheric gravity waves (AGW) including propagation effects along the channel from bottom to top is based on seminal work by [23]. Due to seismic events like earthquakes, volcanic eruptions or tsunamis, acoustic-and atmospheric gravity waves (internal waves; continuously stratified medium) are able to propagate upwards into the ionosphere. Therefore, they could generate oscillations in the ionospheric plasma. The governing equations of motion for adiabatic conditions including inertia-, gravity-, and pressure gradient forces are, see [24], adiabaticity: where v corresponds to the neutral gas velocity, ρ is the density, p the pressure, g the gravitational acceleration, and C s the constant speed of sound. The dispersion relation, assuming harmonic processes, small perturbations, and linearised equations is: If the wave numbers k x and k z = k z + i/2H are real, acoustic waves are determined with an angular frequency ω > ω a = γg 2C s and internal atmospheric gravity waves with ω < ω g = (γ−1) 1/2 g C s . Here γ corresponds to the ratio of the specific heat, ω a is called the acoustic low frequency cut-off and ω g the high frequency cut-off for internal gravity waves, the so-called Brunt-Väisälä frequency (buoyancy frequency). For upward propagation, because of the imaginary part of k z , the amplitudes of these waves grow as exp(z/2H), where H = C 2 s γg is the atmospheric scale height. Kinetic energy can be conserved while the exponential density losses are compensated at the same time, which means that the energy flux remains constant.

Data Processing Method
For a source event with specific date and geographical coordinate a virtual spatial box of usually a few degrees in lat/long is selected. If line of sight projected footprints from satellites are within this area, the magnetic field values from these orbit segments (transits of the area) are considered and selected for further processing. In general, these segments are crossings before and after the event, the extend of the time interval is selectable (usually a few days). All field components for each segment are individually detrended, this corresponds to a high-pass filtering with low corner frequency of f c = 0.02 × f Ny , where f Ny is the Nyquist frequency. A wavelet transform (Morse wavelet, symmetry γ = 3, time-bandwidth product P 2 = 60) converts the components (time series) into a time-scale decomposition. With the absolute values of this representation a binning into frequency bands (scales) is conducted. For Swarm the statistics of the fluctuations are computed over the time dimension of the scalogram for each crossing, i.e., over the time shift for each scale of the wavelet transform (not over different crossings). In case of CSES HPM the segments before the event are merged to a single entity, a second entity are the merged segments after the event. This is necessary to get enough measurements for the further processing, i.e., a higher sample rate is preferable. For CSES HPM the step size is in log-scale according to ln( f n ) = −0.7651 -0.069315×n, where n is the number of steps shown in the ranges n = {1 -50, 1 -60, 1 -63}, for Swarm A/B/C VFM the step size is ln( f n ) = 3.1469 -0.069315 × n, shown are two ranges n = {1 -40, 1 -61}. The significance of the values are higher at the lower end of the range (the higher frequencies), because the fluctuations are usually weak for the upper range but depending also on the strength of the source event. The statistics of the segment bands are compared. For the evaluation and graphical representation (box plots) median values and the 25/75 percentiles are used, outliers are individually drawn, and the boxes/points are colour and shape coded. Segments before (after) the event are in blue (red) colour with fading intensity for larger time differences to the event time, i.e., focus is on crossings in the temporal vicinity of the event, notches indicate the medians. The widths of the boxes are decreasing with time, the furthest crossing in the past (future) has maximum (minimum) elongation. The first crossing after the event is visualised with a red bar, the length marks the percentiles range, and a circle symbol for the median is used. A Mann-Whitney U test with significance level 5 percent for the earthquakes investigated with CSES FGM1 is carried out, flags (yellow stars) at the top of the panels indicate the shift of the medians.
The method is designed to be sensitive to detect magnetic field variations at LEO altitudes which potentially originate from natural hazards. The measured signals in the ionosphere shall carry the imprints of (and close to) lithospheric processes. The conjecture is that the transport of energy and momentum from bottom to top takes place through atmospheric waves.

Results
The observations indicate that after some energetic lithospheric natural and manmade events, magnetic disturbances at LEO altitudes can possibly be traced back to the source region. In the case studies magnetic field variations from satellite measurements are associated with eligible natural sources, among them earthquakes and volcanic erup-tions, and also secondly from artificial factors, in particular rocket launches and strong explosions, see Table 1. The manifold causal relationship via lithospheric atmospheric ionospheric coupling (LAIC) is difficult to establish, but atmospheric acoustic waves (AW) and gravity waves (AGW) can serve as almost line of sight channel to bridge the gap from solid earth to ionospheric heights. This conjecture is supported by theory, observations, and numerical investigations. Measurements indicate main wave activity and dissipation at somewhat lower altitudes, e.g., in case of tsunamis [25], but leaking of energy could lead to pronounced structures and boundary layers at LEO altitudes. If a satellite crosses these structures with sufficient orbital speed the AW or AGW induced pattern in magnetic variation can be considered in first approximation as stationary and sketch out the horizontal dimension along the track.

Volcanoes
Volcanic explosions and effusions, in particular Plinian eruptions with ejections of ashes into stratospheric heights, have the potential to devastate the surrounding community and local environment, significantly disturb and alter the atmosphere for longer time spans, in specific cases even years, with a marked global temperature decrease and a shift of weather pattern. Here, we consider immediate influences in close temporal proximity to reported and impending volcanic activities and investigate magnetic field variations at LEO altitudes in a spatial box of a few degrees in latitude and longitude around the volcano location. In addition, seismic activity prior to the eruption often relates to such cataclysmic processes. Shults et al. (2016) [26] and Heki et al. (2006) [27] reported ionospheric variations, with regard to volcanoes, based on total electron content (TEC) ground observations, [28] made a survey via measurements of the DEMETER microsatellite, and both [29,30] investigated case studies with the CSES HPM magnetic field package. Such conceptions and considerations can be extended beyond the Earth to the solar system [31].

Bárðarbunga Caldera Collapse and Effusive Eruptions
Bárðarbunga is one of Iceland's most active subglacial (under the Vatnajökull glacier) basaltic volcanoes. A specific event series-interconnected with magnitude M > 5 earthquakes and caldera collapse-lasted approximately 6 month, beginning from August 2014 with a seismic swarm to the end of seismic activity in February 2015. The massive lava fissure eruptions in the Holuhraun plain started on 31 August 2014, and were accompanied by huge blasts of SO 2 gas [32]. The geodynamic setting consists not only of a volcanic eruption and caldera collapse event, but is associated with high seismicity and caldera earthquakes over a prolonged period, hence a series of complex intertwined sequences. As a case study, Swarm A/B/C satellite magnetic field measurements are used to investigate these processes in a narrow confining spatial range with a box size of only ±1.5 deg in geographical lat/long around Bárðarbunga's location, and an extended time period of ±10 days centred around 28 August 2014. This time span is also split in a geomagnetic quiet period up to 26 August 2014 (except 19 August 2014), and, thereafter, a disturbed period caused by a geomagnetic storm. All pre-event box crossings are under quiet conditions, i.e., these segments can serve as a low disturbed baseline for post-event crossings. Basically during the magnetic storm phases external influences can have adverse effects, including unwanted spurious signals which interfere with AGW induced variations. The fluctuations of the magnetic field components are studied by wavelet analysis, the corresponding scalograms for one crossing path, together with an illustrative map, are shown in Figure 1.

Earthquakes
Seismic phenomena are a rich area of research and ubiquitous at substantial scales throughout the terrestrial body, in particular together with the propagation and coupling of seismic waves with the environment. We investigate four large earthquakes with moment magnitudes M ≥ 6.5 accessible with CSES and Swarm magnetic field data in a time period a few days around the event in order to get pre-and post-background behaviour and at best foreshock signatures, see Table 2. For common investigations and statistical analysis, according to the United States Geological Survey (USGS) database, 59 earthquakes occurred in the period 1 March 2018 to 31 July 2019. The so-called Pacific Ring of Fire, with its subduction zones, is a hotspot of volcanic and earthquake activity. Here we focus on the M6.5 earthquake near Severo-Kurilsk, Kuril Islands, Russia, where the data enables common CSES and Swarm investigations, see Figure 3. This provides a testbed at the Kuril-Kamchatka trench to compare high magnitude events at different times and hypocenters. CSES FGM1 shows an increase in fluctuations for all three components, i.e., the variations are higher after the earthquake occurrence. Swarm B had only one post-crossing. Swarm A VFM shows higher variations than Swarm C VFM, the three components have varying maxima in different frequency bands (wavelike transition between frequencies and scales).

Earthquake M7.3 Rio Caribe, Venezuela
The investigation is performed using CSES HPM (FGM1 and CDSM) and the data availability is from July 2018 to July 2019. In this period, USGS reported 21 earthquakes with magnitude M ≥ 7 and a wide spread in focal depth, ranging from shallow, e.g., 8 km in the case of the Ridgecrest earthquake sequence, up to deep seated earthquakes with a focal depth of 670.8 km in the case of Suva, Fiji. The Rio Caribe, Venezuela, EQ (USGS ID us1000gez7) had magnitude M7.3 with a focal depth of 146.8 km, the distribution of the magnetic field fluctuations are shown in Figure 4. CDSM shows wavelike features in the frequency bands, i.e., enhancement and attenuation in varying order. For FGM1 the fluctuations in the Y-component are larger than in the X-and Z-components, but for all three components the fluctuations over the whole frequency range are higher after the earthquake.  Figure 2, the geomagnetic activity is |Dst| ≤ 20 nT for the selected time span.

Earthquake M8.2 Ndoi Island, Fiji
The Ndoi Island, Fiji, EQ (USGS ID us1000gcii) had a magnitude M8.2 with a focal depth of 600 km, the distribution of the magnetic field fluctuations are shown in Figure 5. CDSM shows a wavelike behaviour along the frequency bands, the upper bands (lower frequencies) are dominated by the pre-earthquake crossings (higher blue bars according to the colour coding). The fluctuations for FGM1 are generally low, the X-and Z-components along the frequency bands are attenuated after the earthquake.  Figure 2, geomagnetic activity is |Dst| ≤ 37 nT for the selected time span.

Earthquake M7.1 Ridgecrest, USA
The Ridgecrest, USA, EQ (USGS ID ci38457511) had a magnitude M7.1 with a focal depth of 8 km and is the largest one from the corresponding seismic sequence. The distribution of the magnetic field fluctuations are shown in Figure 6. CDSM shows a wavelike behaviour, FGM1 a sinusoidal modulation along the bands.

Explosions
Among several natural and artificial phenomena, not only natural hazards could be responsible for acoustic waves and atmospheric gravity waves. Strong explosions and blasts, i.e., impulsive point sources, have a similar effect like powerful volcanic eruptions with an instantaneous energy release and could, therefore, generate gravity and acoustic waves too. In the recent past the devastating explosion from Beirut, Lebanon, on 4 August 2020, was detected with ionospheric total electron content (TEC) measurements from ground based Global Navigation Satellite System (GNSS) receivers [33].

Baumgarten Natural Gas Hub Technical Failure
A dangerous incident in Central Europe was a sudden gas release together with a serious fire in Baumgarten, Austria. On 12 December 2017, a technical failure in the natural gas hub effectively resulted in an explosion, which could be detected with a seismic array [34]. In order to investigate this precisely timed accident from LEO altitudes the box size has to be at least ±3.5 deg in lat/long to have a minimum of two satellites in the time span ±1 day around the event. There are small magnetic field variations in the traces, see Figure 7 (shown for Swarm A VFM after wavelet analysis), in particular an enhancement in the Y-component, but the time-and distance-gap between the satellite traces and the location is too large to assign these features to this event.  Figure 2, the geomagnetic activity is |Dst| ≤ 29 nT for the selected time span.

Tlahuelilpan Pipeline Explosion
Another man-made catastrophe was the gasoline pipeline explosion in the town of Tlahuelilpan in the Mexican state of Hidalgo on 18 January 2019. According to eyewitness reports a fuel spill ignited into a fireball with tragically loss of life for many people in the immediate vicinity. Such a shock wave at the ground has the ability to release energy into far distance. For CSES HPM the variations are shown in Figure 8. For most of the frequency bands CDSM has an enhancement, i.e., higher values after the explosion. The pre-and postvariations for FGM1 are similar, the Y-component values are larger, all three components have a slightly wavelike feature (redistribution between scales) along their particular frequency band axis.  Figure 2, the geomagnetic activity is |Dst| ≤ 15 nT for the selected time span.

Rocket Launches
Acoustic shock waves, which are basically similar to processes in explosions, have the peculiarity to initiate ionospheric disturbances. Striking examples are rocket launches with wave generation in the exhaust wake during supersonic ascents, measurements and results are inter alia from ionospheric sounders and Global Positioning System (GPS) receiver networks [35]. In two case studies we investigate if the disturbances are detectable with magnetic field measurements at LEO altitudes.

Falcon 9 Launch
The statistical evaluation of Swarm A/B/C VFM magnetic field fluctuations for the SpaceX Falcon 9 launch on 2020-11-21T17:17:08 UTC (payload Sentinel-6 Michael Freilich satellite) from Vandenberg Air Force Base, California, USA, is shown in Figure 9. In general the components are similar, but a distinct attenuation of fluctuations over the entire frequency band is observed in the Y-component of Swarm C VFM, the Y-and Z-components of Swarm B VFM are clearly enhanced. In this particular event with the given parameters (box size lat/long ± 5.5 deg, time span ± 1 day) there are four crossings for each Swarm satellite. On the basis of the orbit constellation Swarm A/C are close in the temporal and spatial domains, which basically enables an intended dual satellite investigation. However, for this event the two post crossings have time gaps of ∆t = {+1.15} days, i.e., ionospheric disturbances triggered by the rocket launch are presumably dissipated, magnetic field fluctuations cannot be linked to the launch. Swarm B has one crossing before and three crossings after the SpaceX Falcon 9 lift-off, ∆t = {−0.40, +0.11, +0.59, +1.10} days, the measurements show distinct enhancements in the fluctuations for the first crossing after the launch. We contend that this magnetic field variation can be linked to ionospheric structures caused by atmospheric waves generated along the first part of the launch vehicle ascent trajectory. All three components for the first crossing after the launch (∆t = {+0.11} days) are enhanced (bars in red colour). The colour coding is the same as in Figure 2, the geomagnetic activity is |Dst| ≤ 38 nT for the selected time span.

Ariane 5 Launch
On 25 July 2018, an Ariane 5 ES rocket was successfully launched into orbit from the Guiana Space Centre, near Kourou, French Guiana, aboard with four navigation satellites from Galileo GNSS. We investigate this rocket ascent (VA-244, and further launches) from this site near the equator with tropical climate, because the defined sequences repeatedly occur and, therefore, enable the possibility of launch event comparisons several times per year. The fluctuation statistics for a Swarm A/B/C VFM scalogram investigation is shown in Figure 10. Swarm A VFM has an attenuation in the Y-and Z-components, Swarm B/C VFM measurements are similar in all components.  Figure 2, the geomagnetic activity is |Dst| ≤ 174 nT (start and main phase of a magnetic storm) for the selected time span.

Discussion
By applying novel satellite magnetic field techniques and approaches on existing state-of-the-art best practises it is possible to circumspectly answer some of the most prominent questions in natural hazards investigations. Magnetic field variations in the ionosphere can shed light on processes in the lithosphere, Earth's surface and above ground. In any configuration the orbital recurrence intervals are crucial points in order to get many crossing segments in the survey area of few degrees around a target. A constellation of satellites like ESA's Swarm mission has inherently advantages, the CSES-2 expansion, together with CSES, will double the revisiting frequency and is very encouraging. The highres 50 Hz sampling enables statistical analyses of each individual segment, either before or after an event. A 1 Hz temporal resolution, due to limited number of vectors per single crossing within a selected investigation area (box size), induces a merging of pre-and post-event periods (which could slightly differ), i.e., only two time spans are effectively available for the evaluation. An analysis to remove confounding effects is obligatory, among them spacecraft interferences and external sources like magnetospheric (mainly in high latitudes) and solar influences on the ionospheric F-region plasma. In that sense the practical experience of instrument teams and amassed evidence through real-time services in the field of near Earth space can provide important knowledge to detect spurious signals. Future research directions include a thorough characterisation of the AW-and AGW-channel into the upper atmosphere and the coupling mechanisms between neutral particles and the ambient plasma population. By combining a variety of multi-parameter methods to describe the phenomena more details will show up. An outlook is to transform these results with an all-encompassing perspective into an essential space based service to leverage natural hazards evaluation and risk assessment.

Materials
The ESA Swarm mission data access and data policy is publicly available via the web page earth.esa.int/web/guest/swarm/data-access (accessed on 19 March 2021), L1b-data are used, in particular the magnetic field high-res 50 Hz files with their meta-data (stored in a single file per day per satellite A/B/C) under the "latest baseline folder". The CSES data access and data policy is by the link www.leos.ac.cn (accessed on 19 March 2021), a registration and completion of a request form is necessary, HPM L2-data in Hierarchical Data Format version 5 (HDF5) are mainly used; a second source are data (calibrated according the nominal standards, i.e., L2 pipeline) from intercalibration campaigns of the CSES HPM instrument suite (checked for spurious satellite disturbances). In both cases ordinary (ftp-)clients and unzip-software are sufficient to get the requested data products. The Swarm datasets are in the Common Data Format (CDF), ordinary software is used to read the data. The earthquake catalogue is from the United States Geological Survey (USGS) via the application programming interface earthquake.usgs.gov (accessed on 15 March 2021), a script handles the on-demand access (query string and result form). Information and activity status about volcanos are from VolcanoDiscovery web.volcanodiscovery.com (accessed on 22 March 2021) and from the Global Volcanism Program of the Smithsonian Institution volcano.si.edu (accessed on 22 March 2021). In both cases textual parts and descriptions are manually translated into machine readable sequences (mainly about eruption type and timing aspects). The events presented in this paper were checked with respect to geomagnetic disturbances. The calculated geomagnetic indices were made available by International Service of Geomagnetic Indices (ISGI) isgi.unistra.fr (accessed on 29 March 2021) collaborating institutes from data collected at magnetic observatories. The map is generated with the software GMT 5 [36]. Programmes (scripts and functions) are written with the software Matlab (version R2019a; "The MathWorks, Inc.; Natick, MA, USA") under Windows 10 ("Microsoft Corporation, Redmond, USA") operating system.

Conclusions
In most of the presented cases and events the statistical results from the wavelet analysis show a change in the pattern between pre-and post-event periods. Wavelike features are found for the China Seismo-Electromagnetic Satellite (CSES) Coupled Dark State Magnetometer (CDSM) scalar measurements. The analysis of the components from the vector field instruments, i.e., CSES fluxgate magnetometer (FGM1) and Swarm A/B/C Vector Field Magnetometer (VFM), show enhancements and attenuations along specific frequency bands. The spatio-temporal confinement to a few degrees in lat/long directs the focus to a limited area, which could be a part (diameter) of a complex atmospheric wave propagation channel from bottom to top, in a certain sense a cross section from a transmission tube which is traversed by satellites at Low Earth Orbit (LEO) altitudes. The synchronous behaviour of the Swarm A/C VFM components for the Bárðarbunga (Iceland) volcano sequence suggests a transit through the same structure, observational biases in causing this trend are unlikely. Enhancements, mainly for the earthquake events, see Figure 11 for CSES FGM1 case studies, indicate energy transfer to the plasma, the wavelike feature a redistribution between different scales and later on dissipation occurs to background levels. The small enhancements from the explosion events could be related to compressional signatures from acoustic waves, but overall the pattern are weak. For the rocket launches at least one component shows alterations, energy transfer into the plasma at LEO altitudes related to acoustic waves is a reasonable explanation.

Data Availability Statement:
Publicly available datasets were analyzed in this study. This data can be found and accessed as stated in the section "Materials".

Conflicts of Interest:
The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript: