Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

The unresolved mystery of dust particle swarms
within the magnetosphere

Max Sommer1,2
( 1University of Stuttgart, Germany
2University of Cambridge, UK
13 May 2024
)
Abstract

Early-generation in-situ dust detectors in near-Earth space have reported the occurrence of clusters of sub-micron dust particles that seemed unrelated to human spaceflight activities. In particular, data from the impact ionization detector onboard the HEOS-2 satellite indicate that such swarms of particles occur throughout the Earth’s magnetosphere up to altitudes of \qty60000—far beyond regions typically used by spacecraft. Further account of high-altitude clusters has since been given by the GEO-deployed GORID detector, however, explanations for the latter have so far only been sought in GEO spaceflight activity.

This perspective piece reviews dust cluster detections in near-Earth space, emphasizing the natural swarm creation mechanism conjectured to explain the HEOS-2 data—that is, the electrostatic disruption of meteoroids. Highlighting this mechanism offers a novel viewpoint on more recent near-Earth dust measurements. We further show that the impact clusters observed by both HEOS-2 and GORID are correlated with increased geomagnetic activity. This consistent correlation supports the notion that both sets of observations stem from the same underlying phenomenon and aligns with the hypothesis of the electrostatic breakup origin. We conclude that the nature of these peculiar swarms remains highly uncertain, advocating for their concerted investigation by forthcoming dust science endeavours, such as the JAXA/DLR DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT mission.

1 Introduction

Since the first reliable and highly sensitive in-situ dust detectors were deployed in the 1970s, the clustering of dust particles, that is, the occurrence of multiple impacts on timescales of minutes or hours, has been a commonly observed phenomenon in near-Earth space. In particular, it was the impact ionization detectors onboard the satellites Prospero and HEOS-2 that gave the first conclusive account of such clusters (Bedford and Sayers, 1973; Hoffmann et al., 1975a). This detector type registers the plasma generated in hypervelocity impacts of micron and sub-micron particles, and posed a major improvement over previous designs, such as the microphone-type sensors, which proved to be highly susceptible to the detection of spurious events (McDonnell, 1978). The detected clusters were interpreted as spatially confined swarms of micrometeoroids, encountered by the spacecraft. Such formations of particles would be expected to rapidly disperse, and thus, an active mechanism that continuously created these clusters would be required. This led Fechtig and Hemenway (1976) to propose that the clusters are generated by the fragmentation of larger meteoroids, for instance, through phase-change-induced or electrostatically-driven disruption triggered upon entering the ionosphere. Analysing the HEOS-2 dust counter data, which covered a broad range of distances from Earth due to the satellites highly elliptic orbit, Fechtig et al. (1979) (in the following referred to as F79) conclude that rather confined particle swarms, presumably generated by the electrostatic breakup of fluffy meteoroids, occur throughout the Earth’s magnetosphere, that is, up to altitudes of roughly \qty60000. Further account of clustering was given by several other dust counter instruments, in particular in low Earth orbit (LEO), such as those onboard the Long Duration Exposure Facility (LDEF, launched 1984) and other satellites (see Section 2.1), as well as at altitudes beyond LEO by the Munich Dust Counter (MDC, 1990) and the Geostationary Orbit Impact Detector (GORID, 1996).

Some clusters detections during these later missions could be attributed to human spaceflight activities, in particular, to specific firings of solid rocket motors (SRMs) (e.g., Schobert and Paul, 1997). SRMs produce large amounts of solid combustion products in the form of micron-sized aluminium oxide particles, which, if used for geostationary orbit (GEO) insertions or retrograde LEO deorbiting manoeuvres, can create circumterrestrial streams of micro-debris with lifetimes of up to months (Bunte, 2003; Stabroth et al., 2008). As the use of orbital SRMs ramped up in the late 1970s and 1980s (to enable utilization of the GEO, see e.g., McDowel, 1997; Wegener et al., 2004), they gained attention for their potential to cause hazard to spacecraft through their generated micro-debris trails (Mueller and Kessler, 1985; Akiba and Inatani, 1990)—amidst the growing awareness of the space debris problem in general (Kessler and Su, 1985; Kessler, 1991). In fact, the analysis of impacts on retrieved spacecraft surfaces showed that, in LEO, the debris dust flux dominates over the natural dust flux (Laurance and Brownlee, 1986; Graham et al., 2001; Hörz et al., 2002).

Following this increased scrutiny of dust sources from human activity, clustered impacts recorded by the GEO-deployed GORID were generally attributed to micro-debris (Drolshagen et al., 2001a, b), even though only few such events could be directly linked to SRM firings. It appears that earlier hypotheses regarding natural particle swarms occurring throughout the magnetosphere, as indicated by HEOS-2 data—remarkably, recorded prior to the extensive utilization of the GEO—were not considered in the analysis of the GORID data. However, aspects of HEOS-2 and GORID-detected clusters suggest that both studies may have observed the same phenomenon. In this perspective, we first review the measurements of particle clusters in near-Earth space in general in Section 2. In Section 3, we then reassess the observations of HEOS-2 and GORID of clusters observed beyond LEO (i.e., those by HEOS-2 and GORID), under the premise that both instruments have observed the ‘magnetospheric swarms’ postulated by F79. We also provide an outlook for the future investigation of this phenomenon by the upcoming DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT mission.

2 Measurements and interpretations of near-Earth dust clusters

2.1 At LEO altitudes

The multitude of dust-counter-type instruments deployed in LEO gathered extensive evidence for clustered impacts of submicron and micron-sized particles at such altitudes. The first such account was given by the impact ionization detector onboard the Prospero satellite, which registered large variations in the daily impact rate, suggesting that around 64% of impacts occurred in clusters (Bedford et al., 1975a, b). In search for explanations, there was cautious speculation that meteoroids grazing the Earth’s atmosphere could break up and shed a swarm of sub-micron fragments, which would then be intercepted by the spacecraft (Bedford et al., 1975a). (This mechanism was later ruled out for swarms observed at higher altitudes, see Section 2.2.1.) Similarly, the capacitor-type dust detector onboard Explorer 46 reported variations in the daily impact rate, which were attributed to the presence of submicron particles within known meteor showers (Singer and Stanley, 1980)—contrary to the theoretical expectation of micro-particle-depleted streams due to the action of radiation pressure (e.g., Dohnanyi, 1972).

Conclusive evidence for the occurrence for compact particle clusters in LEO was eventually given by the high-time-resolution dust counters onboard LDEF, which occasionally recorded bursts of impacts with instantaneous rates up to four orders of magnitude higher than the mean rate (Oliver et al., 1995). Some of these bursts were seen repeatedly at certain points along the spacecraft orbit, indicating intersections with distinct circumterrestrial dust streams, thus interpreted to stem from human activity (Cooke et al., 1995). Two of those repeating clustering events could then indeed be linked to the exhaust dust streams of specific SRM firings (Schobert and Paul, 1997; Stabroth et al., 2007). However, the origin of most of the detected clusters remained unresolved.

Further confirmation that the majority of dust impacts in LEO occurs in clusters, some of which representing circumterrestrial streams, has since been given by the dust-counter-type instruments SPADUS (Tuzzolino et al., 2005), DEBIE-1 & 2 (Schwanethal et al., 2005; Menicucci et al., 2013), and SODAD-1 & 2 (Durin and Mandeville, 2009; Durin et al., 2016, 2022). One of the streams detected by SPADUS could be attributed to a known SRM firing (Neish et al., 2004; Bunte and Drolshagen, 2005), whereas another one appeared to be associated with the explosion of a launch vehicle upper stage (Tuzzolino et al., 2001).

2.2 Beyond LEO

While the occurrence of clusters in the densely-populated LEO region may plausibly be the result of human activity in a number of ways (i.e., collisions and explosions of satellites or debris, orbital SRM firings, substances released by crewed spacecraft, etc.), the occurrence of clusters at higher altitudes is more difficult to explain. Nonetheless, evidence for their existence has been gathered by a few in-situ dust detectors, which are reviewed in the following.

2.2.1 HEOS-2

The first and most extensive account of particle clusters at altitudes beyond LEO was given by the dust counter onboard HEOS-2. The HEOS-2 spacecraft, launched in 1972, carried an impact ionization detector on a highly eccentric orbit with varying perigee in a range of \qtyrange[range-units=single,range-phrase=–]3503000 and apogee of \qty240000km, as illustrated in Figure 1. Due to that orbit, HEOS-2 could sample the dust flux over a wide range of altitudes, from an almost ‘pristine’ interplanetary background flux around its apogee, all the way to the near-Earth dust flux at only a few \qty1000 altitude. Besides an interplanetary sporadic flux described by random impacts, the HEOS-2 data indicated the presence of two types of cluster phenomena: (1) the ‘swarms’, which were characterized by short intervals between consecutive impacts of up to 15 minutes, and (2) the ‘groups’, with longer intervals of up to 12 hours (Hoffmann et al., 1975a, b; Fechtig et al., 1979).

Refer to caption
Figure 1: Illustration of the HEOS-2 orbit through the plasma and field regions of the Earth’s magnetosphere. Reproduced from Fechtig et al. (1979). (Due to copyright restrictions this figure is only included in the journal version of this article.)111Journal version of this article: https://doi.org/10.1098/rsta.2023.0370

The swarms were exclusively encountered at distances of up to 10 Earth radii (R_⊕), which is about the extent of the Earth’s magnetosphere dipole field (about twice the GEO altitude). These ‘magnetospheric swarms’ generated an averaged particle flux rate about one order of magnitude above the sporadic background. Examining different sensor viewing directions perpendicular to the Sun—namely, the Earth apex, antapex, ecliptic north, and south—HEOS-2 also revealed that the swarms occurred anisotropically, with a slight preference for the Earth apex direction (F79). The derived masses of the individual particles ranged mostly within \qtyrangee-14e-12,222 Corresponding to particle radii of about \qtyrange[range-units=single,range-phrase=–]100500, or β𝛽\betaitalic_β factors of in the order of 0.5–1. similar to the largest of the hyperbolic \textbeta-meteoroids irradiating from the Sun (e.g., Wehry and Mann, 1999), which however the HEOS-2 sensor was insensitive to, due to its pointing. Impact speeds of the swarm particles derived by rise time analysis of the impact charge signal were found to be in the order of \qtyrange[range-units=single,range-phrase=–]416\per.

F79 show that the swarms could be caused by the breakup of fluffy interplanetary meteoroids that fragment due to inner, repulsive electrostatic forces (a mechanism first described by Öpik, 1956), resulting from being charged up to high potentials while traversing the Earth’s magnetosphere. From the number of registered impacts and the total encounter duration of each detected swarm, they estimate the total numbers of particles per swarm to fall in the range from 1014 to 1019rangeE14E19{10}^{14}{10}^{19}start_ARG start_ARG end_ARG start_ARG ⁢ end_ARG start_ARG power start_ARG 10 end_ARG start_ARG 14 end_ARG end_ARG end_ARG to start_ARG start_ARG end_ARG start_ARG ⁢ end_ARG start_ARG power start_ARG 10 end_ARG start_ARG 19 end_ARG end_ARG end_ARG, corresponding to individual parent meteoroid masses between \qty10 and \qty1000\kilo (F79), with a geometric mean of \qty5.2\kilo. (Note that the largest observed swarm, corresponding to a progenitor of nearly \qty1000\kilo, appears to be divided into four sub-swarms, suggesting the breakup of multiple smaller objects in close temporal proximity.) Such meteoroids should be observable as visible bright meteors, also called ‘fireballs’, when entering the atmosphere. The meteoroids producing the Ceplecha type III fireballs, which show high ablation ability and are thus thought to stem from low-bulk-density meteoroids (Ceplecha, 1977; Ceplecha et al., 1998), are considered by F79 the most likely candidates for the swarm progenitors. These low-density-type meteors, however, are associated with cometary source orbits that generally produce high atmospheric entry velocities (\qtyrange[range-phrase=–,range-units=single]2070\per), which is incompatible with the rather low-velocity swarms. F79 acknowledge a lack of low-velocity type III fireballs, which could represent the sought-after swarm progenitors, but note that this apparent deficiency could potentially be explained by selection effects in meteor observations as well as swarm production. Based on the swarms’ directionality, they rule out creation via aero-fragmentation (the speculated origin of the Prospero clusters, see Section 2.1) or a relation to meteor showers (Dohnanyi and Fechtig, 1977; Fechtig et al., 1979).

The more spread-out ‘groups’, on the other hand, were detected at all distances from Earth (i.e., up to the satellite’s apogee of around \qty240000), yet crucially, they appeared almost exclusively when the Moon was within the sensor’s field of view (Hoffmann et al., 1975a; Fechtig et al., 1979). Backtracing of the group particles’ trajectories further substantiated the idea of a lunar origin (Hoffmann et al., 1975b). Consequently, the authors interpreted these groups as originating from occasional larger meteoroid impacts on the Moon that would create significant amounts of ejecta escaping lunar gravity, subsequently roaming through the Earth-Moon system in somewhat dispersed formations of particles. By modelling the motion of escaping ejecta clouds through the Earth-Moon system, it was further concluded that their expected spatial dimensions are consistent with HEOS-2’s encounter times of the groups, and that lunar impactors in the order of \qty1 could suffice to create the observed group number densities (Dohnanyi, 1977). Attempts to correlate the occurrence of groups with lunar impact events sensed by Apollo-deployed seismometers were inconclusive, due to the much higher number of seismic events compared to group detections (F79). In this paper, we will not consider these lunar-origin clusters further, and instead focus on the near-Earth swarms observed by HEOS-2.

2.3 Munich Dust Counter

The findings of HEOS-2 are reinforced by the data of the impact-ionization-type Munich Dust Counter (MDC) onboard Hiten, launched in 1990. As HEOS-2, the Hiten spacecraft was placed on a highly eccentric orbit with initial perigee of a few \qty1000 and apogee of \qty300000. MDC, similarly to the HEOS-2 detector, reported a ten-fold increase of flux in regions near Earth (closer than \qty100000) compared to interplanetary space and an anisotropy in the flow direction, with a preference for the Earth apex direction (Iglseder et al., 1993a, b). In addition, strong clustering with instantaneous particle rates of up to five orders of magnitudes over the mean flux was observed. Unfortunately, however, no further details about the MDC-detected clusters have been published.

2.4 GORID

The Geostationary Orbit Impact Detector (GORID) was a refurbished engineering model unit essentially identical to the impact ionization dust detectors flown on the interplanetary probes Ulysses and Galileo (Grün et al., 1992a, b). GORID was mounted on the Russian telecom satellite Express-2, which was placed in GEO in 1996 (Drolshagen et al., 1999). GORID data indicated a significant excess of clustered impacts, which were characterized by separations in time from minutes to an hour, dominating over random impacts by a factor of about four (Drolshagen et al., 2001a; Graps et al., 2007). Similarly to the case of LDEF, few clusters were observed recurrently at certain points along the orbit, and could be dynamically linked to the exhaust dust streams of specific SRM firings (Bunte, 2003; Bunte and Drolshagen, 2005). The vast majority of detected clusters, however, were non-recurrent and could not be attributed to any specific source. This led Bunte and Drolshagen (2005) to conjecture an unknown micro-debris-cloud-generating mechanism present at GEO altitudes. In that sense, it was speculated that the breakup of larger SRM slag particles via electrostatic fragmentation could act as such a mechanism (Graps et al., 2007).333 SRM slag particles are larger (up to cm-sized) aluminium-oxide clumps ejected at the end of an SRM burn with low relative velocity (Jackson et al., 1997). Fragments of slag produced in GEO insertions (consequentially exhibiting near-GEO orbital parameters) would have too-low velocities relative to GORID (few \qty100\per) to be efficiently detected via impact ionization. However, GTO insertions would produce slag with velocities relative to GORID of \qtyrange[range-units=single,range-phrase=–]1.52\per.

Ultimately, all GORID-detected clusters were attributed to debris-related phenomena (Drolshagen et al., 2001a; Graps et al., 2007) with no consideration of the possibility that the clusters are of natural origin. One reason for this is that the impact velocities inferred from charge signal rise times indicate most impactors to have slow speeds (\qty¡5\per), which would be compatible with relative speeds of debris in GTO-like orbits. However, it is noted that these impact speeds had been found to be unreliable for a number of reasons. In particular, impacts on the insides of the detector’s side walls instead on the impact target would generate events with abnormally long rise times, and thus appear as slow, even though they are not (see also Stübig, 2002; Altobelli et al., 2004; Willis et al., 2004, 2005). Later analysis of the GORID data therefore ignores the rise-time-derived speeds altogether (Graps et al., 2007). By geometrical consideration, one might expect that around half of the impacts inside the detector occur on the side walls (assuming an isotropic flow).444 The effective solid angle ΩeffsubscriptΩeff\Omega_{\mathrm{eff}}roman_Ω start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT of GORID including its inner side walls can be approximated by an unobstructed plate detector with the same nominal sensitive area (all particles entering the instrument are detected), which has Ωeff=πsrsubscriptΩeff𝜋sr\Omega_{\mathrm{eff}}\,=\,\pi\,\mathrm{sr}roman_Ω start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT = italic_π roman_sr. The ΩeffsubscriptΩeff\Omega_{\mathrm{eff}}roman_Ω start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT of the GORID impact target is \qty1.45sr (Grün et al., 1992a). Thus, the ΩeffsubscriptΩeff\Omega_{\mathrm{eff}}roman_Ω start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT of just the inner walls (\qty1.69sr, calculated as the remainder of the two) is comparable to that of the impact target.

One striking feature of the GORID-observed clusters is the strong modulation of their occurrence with the local time of the satellite. Clusters were detected with a strong preference around local midnight, that is, when the satellite was on the Sun-opposing side of the Earth, and the spacecraft-body-fixed sensor was sensitive to a flux from the Earth apex direction (Bunte and Drolshagen, 2005; Graps et al., 2007). For the half of the orbit, when the sensor was pointing away from apex, virtually no clusters were detected. There is no obvious reason why debris would preferentially be located around local midnight, although Graps et al. (2007) note that an unknown process occurring only in the magnetotail region might facilitate debris cluster detections.

Another noteworthy feature is that the distribution of the impact-generated ion grid charges (‘Qisubscript𝑄iQ_{\mathrm{i}}italic_Q start_POSTSUBSCRIPT roman_i end_POSTSUBSCRIPT’) shows a greater tendency towards smaller charges for clustered impacts compared to random events (Graps et al., 2007). This trend could imply a steeper mass distribution, possibly indicating an increased prevalence of smaller particles within clusters.

Lastly, measurements of GORID’s charge-senstive entry grids indicated that cluster particles were highly negatively charged (Bunte et al., 2006; Graps et al., 2007), as opposed to the positive equilibrium potential particles are thought to assume in interplanetary space as well as inside the magnetosphere under average conditions (Horányi, 1996). Positive charges, however, occurred only among the randomly impacting particles, according to the GORID measurements.

3 Discussion

3.1 HEOS-2 and GORID clusters - a common origin?

Having reiterated the main findings about dust particle clusters in near-Earth space, we will now reassess them in context to each other. In particular, we will attempt to reconcile the high-altitude clusters observed by HEOS-2 and GORID, which have garnered notably different interpretations. The fact that both instruments reported an apparent anisotropy favouring the Earth apex, and typical cluster encounter durations of up to \qty90minutes suggests that they are related to each other, and thus, that they could stem from the same phenomenon. In the following, we will discuss several aspects that could support this notion—specifically, F79’s hypothesis of the magnetospheric swarms—and weigh in on potential swarm progenitor candidates, as well as other explanations for the observed clusters.

3.1.1 Anisotropy

Both instruments reported preferential detection of clusters when exposed to the Earth apex direction, although the observed anisotropy was arguably more pronounced for GORID than for HEOS-2 (compare, e.g., F79 Tab. 7 with Graps et al. (2007) Fig. 2). Interestingly, although the GORID data indicate a strong preference for cluster detection when exposed to the Earth apex, the effective sensitivity towards that direction remained relatively low. Due to GORID’s body-fixed mounting on the GEO satellite with a boresight pointing only \qty25 away from equatorial north, the minimum angle between boresight and the apex direction (which was achieved once per day, during local midnight on the geostationary orbit) varied between \qtyrange[range-phrase= and ]4288, depending on the time of the year. Thus, considering GORID’s angular sensitivity profile (Grün et al., 1992a), one would expect a distinct seasonal modulation of cluster detections, if the clusters were indeed coming preferentially from the apex direction. The maximum cluster occurrence rate should occur in autumn when the sensor’s exposure to the Earth apex direction is highest (during local midnight hours), whereas the minimum rate should occur in spring. To check this proposition, the mean flux of clustered events during each calendar month averaged over the mission duration is shown in Figure 2, alongside the minimum angle between sensor boresight and the apex direction. Although the highest rates are indeed observed in autumn, the minimum detection rate of clusters occurs in summer. Thus, a clear correlation of the cluster occurrence rate with the sensor’s sensitivity toward the Earth apex direction is not evident.

Refer to caption
Figure 2: Monthly mean impact rates recorded by GORID per calendar month averaged over the entire observing period (1997–2002). Error bars represent the year-to-year standard deviation. Based on data from Drolshagen (2006). Also shown is the minimum angle between the sensor boresight and the Earth apex direction (reached once per day at local midnight) over the course of the year.

A reduced modulation of rates with the exposure to the apex could be due to the sensor’s asserted sensitivity to wall impacts (Stübig, 2002; Altobelli et al., 2004), which considerably widens the angular sensitivity profile. On the other hand, the lack of a clear modulation throughout the year in combination with the observed strong modulation with local hour, indicates that it is not the sensor’s exposure to the apex direction that is the crucial factor for cluster detection, but rather the location of the spacecraft on the tail side of the magnetosphere. The tail side could indeed be preferential for the detection of the putative magnetospheric swarms, due to the proximity to the magnetosphere’s plasmasheet. The plasmasheet, with its elevated electron fluxes, had been suggested by F79 as another potential swarm genesis region, besides the field line regime connected to the auroral zones traversed by HEOS-2 (F79).

3.1.2 Geomagnetic activity

Considering that F79 identified high fluxes of energetic electrons (\qty¿10keV) in the surrounding plasma as a crucial factor for charging up meteoroids enough to cause their electrostatic breakup (F79), it is reasonable to assume that the occurrence of swarms is correlated with geomagnetic activity. Plasma motion and electric currents inside the magnetosphere are enhanced during geomagnetic storms, which are driven by solar wind disturbances, which in turn are produced by solar activity (Richardson et al., 2000). We may thus analyse the occurrence of clusters registered by HEOS-2 and GORID in the context of the solar cycle timing and the geomagnetic activity level.

The observation period of HEOS-2 (1972–1974) occurred well after the peak of solar cycle 20 (1969), and during the cycle’s decline (minimum in 1976). A cycle’s declining phase is typically not associated with major geomagnetic storms, due to the less frequent occurrence of coronal mass ejections (CMEs). The decline of solar cycle 20, however, was marked by a period of exceptionally high geomagnetic activity occurring in 1973–1975, caused by persistent high-speed solar wind streams (Gosling et al., 1977; Richardson et al., 2000; Richardson and Cane, 2012). High fluxes of energetic electrons, in particular, are associated with high-speed streams (Reeves et al., 2011; Boynton et al., 2013). Notably, HEOS-2 detected 13 of the 15 swarms in 1973 and 1974, even though this period corresponded to only 64% of the observation time. To test whether the occurrence of particle swarm events correlates with elevated geomagnetic activity, we can analyse the ‘Kp index’ preceding these events. The Kp index, derived from geomagnetic observatory data at a 3-hour interval, provides a standardized measure of global geomagnetic activity levels on a scale from 0 to 9, with higher values indicating increased disturbance, and has been recorded for more than 70 years (Matzka et al., 2021b). Our analysis shows that, compared to the overall mean Kp index of 2.45 (SD = 1.38) over the mission duration, the mean of the Kp indices directly preceding a HEOS-2 swarm detection is notably higher at 3.46 (SD = 1.19). This difference is statistically significant, as supported by the obtained p-value of 0.005 from a two-sample t-test, suggesting that there is only a 0.5% probability of observing such a difference in means if there was no correlation. These findings indicate a preferential association between particle swarm detections and heightened geomagnetic activity, as reflected by the Kp index. This is also illustrated in Figure 3 (top), which shows the timings of HEOS-2 swarm detections next to the Kp index (averaged for 5-day intervals).

The observation period of GORID (1997–2002) covers the rising phase and maximum of solar cycle 23 (min. in 1996 and max. in 2001). In the case of GORID, our dataset is limited to the monthly mean detection rate of clustered impacts (Drolshagen, 2006), rather than the precise timings of cluster detections. Therefore, we compare the cluster impact rate also to the monthly mean Kp index, as depicted in Figure 3 (bottom). This comparison reveals a correlation between the cluster detection rate and the geomagnetic activity level, quantified by the Pearson correlation coefficient r𝑟ritalic_r. While for the random event rate we obtain r=0.11𝑟0.11r\,=\,-0.11italic_r = - 0.11 (indicating no correlation), the cluster rate shows a moderate correlation with r= 0.57𝑟0.57r\,=\,0.57italic_r = 0.57.

HEOS-2: Refer to caption
GORID: Refer to caption
Figure 3: Cluster occurrence with respect to geomagnetic activity. (Top) Times of swarm detections by HEOS-2 over the entire observing period (obtained from F79), indicated by vertical blue lines. Also shown is the average Kp index for 5-day intervals, indicating the level of geomagnetic activity. (Bottom) Monthly mean detection rates recorded by GORID of random and clustered impacts over the entire observing period (obtained from Drolshagen (2006), considering only ‘class 3’ events, supposed to be real impacts). Also shown is the average Kp index for each month. Kp index data for both plots was obtained from Matzka et al. (2021a, b).

It is important to note that the electron fluxes might not always be directly correlated with the Kp index. In particular, the pumping of fluxes of relativistic electrons (\qty¿1MeV) by geomagnetic storms is known to be time-delayed and can cause elevated levels lasting for several days after a storm has subsided (Baker et al., 1998; Friedel et al., 2002; Boynton et al., 2013). Conditions facilitating swarm formation may thus persist beyond geomagnetic activity peaks. Nevertheless, the fact that both HEOS-2 and GORID detected clusters preferentially during periods of high geomagnetic activity further hints at a common phenomenology.

The notion that swarm creation is associated with high geomagnetic activity is also consistent with the findings of Graps and Grün (2000), who modelled the equilibrium potentials of grains within the magnetosphere for ‘quiet’ and ‘active’ conditions, as well as for different particle materials. They conclude that, while quite conditions usually generate positive potentials (\qty¡15), disturbed conditions can produce highly negative potentials (\qty¡-1000) for materials with low photoelectron yield.

3.1.3 Impact rates

We may also compare the incidence rates of clustered impacts between HEOS-2 and GORID. Table 1 summarizes the key geometrical characteristics of both sensors, the nominal sensitive area A0subscript𝐴0A_{0}italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, the maximum entry angle θmaxsubscript𝜃max\theta_{\mathrm{max}}italic_θ start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT, the effective solid angle ΩeffsubscriptΩeff\Omega_{\mathrm{eff}}roman_Ω start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT, and the geometric factor G𝐺Gitalic_G, which is the product of A0subscript𝐴0A_{0}italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and ΩeffsubscriptΩeff\Omega_{\mathrm{eff}}roman_Ω start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT, and represents the effective gathering power of an instrument with respect to an isotropic flux. These parameters allow us to relate the impact rates onto HEOS-2 and GORID to one another. For a collimated dust stream and a directly upstream-pointing sensor, only the nominal sensitive area A0subscript𝐴0A_{0}italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is relevant, such that GORID would be exposed to about 10 times as many impacts as the HEOS-2 detector. For an isotropic flow of dust (or a randomly varying instrument pointing) the geometric factor G𝐺Gitalic_G is the relevant characteristic, such that GORID would be exposed to about 15 times as many impacts as HEOS-2.

Table 1: Comparison of sensor geometries.

[b] Instrument A0subscript𝐴0A_{0}italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT θmaxsubscript𝜃max\theta_{\mathrm{max}}italic_θ start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT ΩeffsubscriptΩeff\Omega_{\mathrm{eff}}roman_Ω start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT G𝐺Gitalic_G cm2centimeter2{\mathrm{cm}}^{2}power start_ARG roman_cm end_ARG start_ARG 2 end_ARG °degree\mathrm{\SIUnitSymbolDegree}° srsteradian\mathrm{sr}roman_sr cm2 srtimescentimeter2steradian{\mathrm{cm}}^{2}\text{\,}\mathrm{sr}start_ARG power start_ARG roman_cm end_ARG start_ARG 2 end_ARG end_ARG start_ARG times end_ARG start_ARG roman_sr end_ARG GORID (Grün et al., 1992a) 1000100010001000 70707070 1.451.451.451.45 1450145014501450 GORID incl. walls (approx.) 1000100010001000 90909090 3.143.143.143.14 3140314031403140 DDA (Sommer, 2024) 300300300300 45454545 0.480.480.480.48 144144144144 HEOS-2 (Hoffmann et al., 1975b) 95.495.495.495.4 60606060 1.031.031.031.03 98989898

With HEOS-2, a total of 207 swarm particles (split among 15 swarms) have been identified during an accumulated observation time below \qty10RsubscriptRdirect-sum\textit{R}_{\oplus}R start_POSTSUBSCRIPT ⊕ end_POSTSUBSCRIPT of \qty70days (F79). This amounts to an average incidence rate of \qty2.9\per. On the other hand, GORID detected a total of 2477 clustered events during a total observation time of 1827 days (Drolshagen, 2006; Graps et al., 2007), amounting to a mean rate of \qty1.4\per. (This only represents ‘class 3’ events, i.e., those with the highest confidence of being real impacts.) Since the GORID sensor has a gathering power of about 15 times that of HEOS-2, it effectively detected a 30 times lower flux of clustered impacts than HEOS-2.555 This discrepancy increases by a factor of 2similar-toabsent2\sim 2∼ 2 if wall impacts are considered for GORID, due to a doubling of its geometrical factor. The discrepancy decreases by 33% if wall impacts are considered for both instruments. If the clusters detected by HEOS-2 and GORID stem for the most part from the same phenomenon (presumably, the magnetospheric swarms), the discrepancy may be explained by the following reasons:

  • The HEOS-2 sensor was more exposed to the anisotropic flow of swarms than GORID, which both instruments indicated to be coming preferentially from the Earth apex: HEOS-2 pointed straight towards the Earth apex for 40% of the observation time. GORID, on the other hand, was only marginally exposed to the apex direction, as discussed in Section 3.1.1.

  • The HEOS-2 sensor was more sensitive than GORID, and thus could detect more particles. The particle mass detection threshold of the HEOS-2 detector is indeed reported to be one order of magnitude lower than that of GORID at the same velocity: mmin,HEOS=\qty1.2e15subscript𝑚minHEOS\qty1.2𝑒15m_{\mathrm{min,HEOS}}\!=\!\qty{1.2e-15}{}italic_m start_POSTSUBSCRIPT roman_min , roman_HEOS end_POSTSUBSCRIPT = 1.2 italic_e - 15 and mmin,GORID=\qty1.5e14subscript𝑚minGORID\qty1.5𝑒14m_{\mathrm{min,GORID}}\!=\!\qty{1.5e-14}{}italic_m start_POSTSUBSCRIPT roman_min , roman_GORID end_POSTSUBSCRIPT = 1.5 italic_e - 14, both at vimp=\qty10\persubscript𝑣imp\qty10\perv_{\mathrm{imp}}\!=\!\qty{10}{\per}italic_v start_POSTSUBSCRIPT roman_imp end_POSTSUBSCRIPT = 10 (Dietzel et al., 1973; Göller and Grün, 1989). Particle masses derived from HEOS-2 measurements indicate that swarm particles are partially below the GORID threshold (F79).

  • More swarms occurred at the location of the HEOS-2 satellite, which flew through the pole regions of the magnetosphere, as opposed to GORID, which orbited in the equatorial plane. Or, more swarms occurred during the time of the HEOS-2 mission, potentially as a result of varying geomagnetic activity (see Section 3.1.2).

Given the above considerations, the rates as of GORID and HEOS-2 are arguably reconcilable, which upholds the notion that the clusters detected by both instruments could stem from the same phenomenon.

3.1.4 Swarm progenitor candidates

We have already mentioned the apparent deficiency of low-velocity, low-bulk-density meteors, whose progenitors could be the source of the swarms, yet which might not have been observed due to selection effects, as noted by F79. Especially considering the large masses derived by F79 for these progenitors of \qty10 to \qty1000\kilo (geom. mean of \qty5.2\kilo), the lack of their detection as meteors poses a conundrum. In addition, Mendis (1981, 1984) showed that—while conceivable for smaller grains—such large meteoroids would have to exhibit improbably small tensile strengths in order to electrostatically disrupt. Rather than complete disruption, they suggest that irregularly shaped meteoroids would be electrostatically ‘eroded’, as edges and protrusions are more readily chipped by the electrostatic tension (Hill and Mendis, 1981), thereby shedding the material to generate the swarms. As only a fraction of the bulk mass could be shed by this erosion, however, even larger meteoroids would be required to yield the swarm masses estimated by F79.

Similarly speculative, one might consider the possibility that it is not the disruption of a single massive progenitor that produces the swarms, but rather the quasi-simultaneous electrostatic breakup of numerous micrometeoroids within a certain volume of space. Given the asserted correlation of swarm detection with geomagnetic activity (see Section 3.1.2), it is conceivable that local intensity spikes of the plasma environment could cause the disruption of micron- and millimetre-sized particles within a whole region of the magnetosphere (for a review of turbulences of the magnetosphere, see, e.g., Zimbardo et al., 2010). Although only two of the 15 swarms detected by HEOS-2 exhibited a clear substructure (hinting at multiple simultaneous breakup events, F79), such a scenario would lift the high mass requirement for the swarm progenitors, so that other candidates could be considered.

In that regard, Jupiter family comets (JFC) have been found via dynamical modelling to be able to produce also micrometeoroids on prograde orbits with very low semi-major axis and aphelia near \qty1, which would generate low-velocity near-apex meteor radiants (Nesvorný et al., 2011). As cometary meteoroids, which are generally associated with low bulk densities, these could constitute the missing meteoroid subpopulation, that produce the magnetospheric swarms. Moreover, their preferential approach direction from the Earth apex would be consistent with the swarms’ apparent anisotropy. Quantitatively, their occurrence is expected to be only secondary to those from the likewise JFC-generated, higher-velocity helion/antihelion radiants (\qty∼30\per), as well as the yet-higher-velocity north/south apex radiants (\qty∼60\per) sustained by Halley-type and long-period comets (Nesvorný et al., 2011), which, however, might less readily produce swarms, due to their shorter residence times within the magnetosphere.

A population of grains that appear to exhibit similar dynamics has been detected in head-echo radar meteor with the Arecibo Observatory (AO). These slow meteors (\qty∼15\per) emerged when AO pointed near the Earth apex, alongside higher rates of fast meteors from the north/south apex radiants (Janches et al., 2003; Sulzer, 2004; Janches and Chau, 2005). The range of particle sizes constituting the AO dataset is stated to be \qtyrange[range-units=single,range-phrase=–]0.5100, yet the sensitivity of head-echo radar meteor observations to slow particles \qty¡30\per has been a matter of debate and is likely closer to the higher end of the stated range (Fentzke and Janches, 2008; Janches et al., 2014, 2015, 2017).

The dynamics of the swarms may hint at a connection to another group of interplanetary dust grains, namely, the \textalpha-meteoroids (Grün and Zook, 1980; Sommer, 2023). These particles move on highly eccentric orbits with aphelia near \qty1, where they thus exhibit low heliocentric velocities. The Earth effectively overtakes these slow-moving particles, such that they appear to be coming from around the Earth apex with v=subscript𝑣absentv_{\infty}\,=\,italic_v start_POSTSUBSCRIPT ∞ end_POSTSUBSCRIPT =\qtyrange[range-units=single,range-phrase=–]520\per. The \textalpha-meteoroids notably exhibit the dynamical properties demanded for swarm progenitors, that is, a preferred approach direction from the apex, as well as relative velocities below \qty20\per. Yet, these micron-sized particles are only about 1–2 orders of magnitude more massive than the swarm particles themselves.

3.1.5 Other explanations

At the time of the HEOS-2 mission, the phenomenon of SRM dust streams was not yet well-known, which is why they were not considered as explanation for the clusters in the interpretation of HEOS-2 data. As it was established that GEO-insertion SRM firings could create distinct dust formations in orbit around Earth, Fechtig (1984) acknowledged that they could be another explanation for the observed swarms. However, there are reasons why the HEOS-2-detected clusters are unlikely to be SRM exhaust dust formations, namely the altitude and the locations at which they occurred. The swarms were observed at all altitudes up to \qty60000. In SRM GEO insertions, the exhaust dust is expelled rather retrograde such that the maximum obtainable apogee of ejected particles is near the GEO altitude (about \qty36000) (Mueller and Kessler, 1985; Bunte, 2003). On the other hand, the orbit of HEOS-2 was oriented such that the spacecraft reached its far-out apogee near the ecliptic north direction, causing it to spend only very limited time near the equatorial region. (see Figure 1). The north polar region was in use at the time by Russian Molniya satellites. However, those were launched by the Molniya-M launcher, which used a liquid upper stage (notwithstanding the fact that SRM dust ejected in Molniya orbit insertions would fail to reach a stable orbit).

It is conceivable that orbital perturbations could have caused SRM dust formations to drift from their origin region to where the swarms were observed. However, it is questionable whether these displacements could take place in the lifetimes of SRM dust streams, which are limited by atmospheric drag (at most 1–2 months for particles \qty¡1 in size (Friesen et al., 1992; Bunte, 2003)). Moreover, electrodynamic perturbations on Earth-orbiting dust particles tend to dissipate their orbital energy, adding to the orbital decay due to atmospheric drag (Juhász and Horányi, 1997). Even if perturbations can meaningfully move GEO-insertion SRM dust to high altitude polar regions, their stochastic nature would arguably lead to the dispersion of the dust clouds.

Lastly, the Earth-apex anisotropy of the HEOS-2-detected swarms has no obvious explanation in the context of SRM dust, similarly to the preference of local midnight hours for cluster detections by GORID.

An entirely different explanation could be that the consistently detected clusters are not true particle impacts but rather noise events caused by an interaction of the impact ionization detectors with the Earth’s magnetosphere. Spurious detections caused by various environmental factors have been a notorious issue for dust-counter-type instruments, requiring sophisticated data reduction techniques (e.g., Poppe et al., 2011), or invalidating datasets outright (Nilsson, 1966). In order to distinguish true impacts from noise events, a coincidence criterion was used with HEOS-2 and GORID, requiring valid successive charge signals from different measurement channels. With HEOS-2, a two-way coincidence criterion (target and ion collector) was used, which for instance, classified all events registered during a historically powerful solar storm of August 1972 as noise (Dietzel et al., 1973). GORID on the other hand, employed a three-way coincidence criterion (target charge, ion grid, and channeltron) to identify true impacts with considerable confidence (Baguhl et al., 1993; Grün et al., 1995). (Note that GORID is identical to the Ulysses and Galileo detectors.) A significant amount of events classified by GORID as probably noise-induced appears to be linked to an electrostatic interaction with the spacecraft as well as to the operation of its plasma thrusters, as described in (Drolshagen et al., 2001a). Note that usually only the highest-confidence true-impact events have been used in the referred-to analyses of the GORID data. Nevertheless, it is at least conceivable that a so-far unknown interaction of the impact-ionization-type detectors with the magnetosphere’s turbulent plasma environment may cause the clustered occurrence of genuine-looking noise events, despite the scrutiny put in place. Given that electromagnetic effects are evidently conducive to the occurrence of noise events, the apparent correlation of cluster detections with geomagnetic activity (Section 3.1.2) is arguably consistent with that proposition.

3.2 Future investigation with DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT

JAXA’s DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT mission is a planned small body science mission to active asteroid 3200 Phaethon, scheduled for launch in 2025. It carries the DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT Dust Analyzer (DDA), a state-of-the-art dust detector that allows for the simultaneous analysis of particles’ dynamical (via charge-sensitive entry grids) and compositional (via an impact plasma mass spectrometer) information (Simolka and et al., 2024). In addition to studying the dust environment of Phaethon, DDA will conduct science operations during the entire nominal mission duration of 3 years, to study interplanetary and interstellar dust (Krüger et al., 2019). For the first two years of the mission, the spacecraft revolves in Earth-bound orbits, during which it gradually ascends from an initially GTO-like orbit all the way to the Moon’s orbit using its solar electric propulsion system, as shown in Figure 4. Eventually, DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT will conduct a series of lunar gravity assists to escape the Earth-Moon system and further advance towards Phaethon.

DDA will thus have extensive exposure to the near-Earth dust environment at all altitudes beyond LEO. DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT will pass through the magnetosphere at altitudes \qty¡10RsubscriptRdirect-sum\textit{R}_{\oplus}R start_POSTSUBSCRIPT ⊕ end_POSTSUBSCRIPT for roughly the first year of the mission, and then traverse the plasma sheet region for several more months after that. Unlike HEOS-2, DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT will remain within equatorial rather than the polar regions of the magnetosphere, at an orbital inclination of \qty∼31. The near-Earth mission phase of the mission will coincide with the solar cycle 25 maximum, which is expected for mid 2025.

Refer to caption
Figure 4: Trajectory of the initial phase of the DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT mission, during which the spacecraft will gradually raise its orbit towards the Moon. Phases of thrusting and coasting are indicated by thick and thin lines, respectively.

The DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT mission profile combined with the capabilities of the DDA instrument represent a unique opportunity to further investigate the clustering phenomena previously reported by HEOS-2 and GORID. By recording a valid mass spectrum of a particle’s impact plasma, DDA can unequivocally identify true particle impacts. It could thus rule out the possibility of the clusters being spurious noise events. The compositional analysis of swarm particles—should they be found—would then give definitive proof of their natural or artificial origin. Moreover, DDA’s charge-sensing entrance grids may register larger swarm particles (reported by HEOS-2 to have masses up to \qtye-11 (F79), corresponding to sizes of roughly \qty∼1), allowing for a direct measurement of their charge and velocity. Such measurements would undoubtably be crucial for understanding the swarms’ formation mechanism.

Similarly to how we related the incidence rates of HEOS-2 and GORID to each other, we may tentatively estimate the incidence rate of swarm particles onto DDA by converting the HEOS-2 and GORID rates (see Section 3.1.3) to DDA rates using the relation of their respective geometrical factors (see Table 1). This yields 133 swarm particle impacts per month onto DDA according to HEOS-2 rates and 4.2 impacts per months according to GORID rates (for the time within the magnetosphere). Identifying the reason(s) for this discrepancy (as speculated upon in Section 3.1.3) would be a key aspect of the investigation. It’s important to note that this simple calculation ignores the differences in spacecraft orbits and instrument pointings.

4 Conclusion

We have provided an overview of the various reports of clustered impact events given by dust-counter-type instruments in near-Earth space. The clustered impacts detected by the HEOS-2 and GORID impact ionization detectors within the magnetosphere and at GEO, respectively, represent a peculiar phenomenon, which have been interpreted as transient particle swarms created either by the breakup of meteoroids, or by GEO spaceflight activities. We have assessed the possibility that the swarms observed by both instruments are caused by the same phenomenon, which is indicated by commonalities regarding their encounter durations and their anisotropy. Remarkably, we have found that in both cases their occurrence was correlated with high geomagnetic activity, which is consistent with the notion that the swarms are caused by the electrostatic breakup of meteoroids proposed by F79, but could also be interpreted as in favour of the clusters being spurious, plasma-environment-induced noise events. While the GORID-detected clusters might also stem from the electrostatic breakup of SRM slag particles in GTO-like orbits (as speculated by Graps et al., 2007), this scenario can hardly explain the HEOS-2-detected swarms, which were observed far from the equatorial plane and at altitudes of up to \qty60000.

Lacking the proper dynamical and compositional characterization of dust particles—which both HEOS-2 and GORID were unable to provide—the nature of the swarms remains elusive. The upcoming DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT mission, however, will carry a sophisticated dust analyser instrument (the DDA), which will offer these capabilities. Considering the mission’s prolonged orbit raising phase through the magnetosphere, DESTINY+{}^{\scriptscriptstyle+}\mkern-2.5mustart_FLOATSUPERSCRIPT + end_FLOATSUPERSCRIPT/DDA poses a unique opportunity to finally unveil the origin of the putative particle swarms.

In light of the new insights presented here, it may be worth to revisit in particular the GORID data, and analyse the timings and locations of cluster detections more thoroughly with respect to geomagnetic activity and region within the magnetosphere. Exploring dust counter data retrieved from the magnetospheres of Jupiter and Saturn for similar swarms could provide additional avenue for investigating this phenomenon. While most observed clusters there were linked to dynamic nano-dust streams from the moons Io and Enceladus, there have also been reports of impact clusters of unresolved origin seen by the Cassini Cosmic Dust Analyzer (CDA) in orbit around Saturn (Hsu et al., 2011).

Finally, new insights about electrostatic effects in cosmic dust, in particular the ‘electrostatic lofting’ of grains from dusty surfaces suspended in a plasma or UV light (e.g., Wang et al., 2016), could help shed light on the formation mechanism of the magnetospheric swarms.

Acknowledgement

I thank Simon Green and Tony McDonnell for their assistance in acquiring valuable literature about the GORID results, as well as three anonymous reviewers for their constructive comments. Partial support by the German Aerospace Center (DLR, grant no. 50OO2101) is gratefully acknowledged.

References