An Anisotropic Cosmic-Ray Enhancement Event on 07-June-2015: A Possible Origin

A usual event, called anisotropic cosmic-ray enhancement (ACRE), was observed as a small increase (≤5%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}${\leq}\,5\%$\end{document}) in the count rates of polar neutron monitors during 12 – 19 UT on 07 June 2015. The enhancement was highly anisotropic, as detected only by neutron monitors with asymptotic directions in the southwest quadrant in geocentric solar ecliptic (GSE) coordinates. The estimated rigidity of the corresponding particles is ≤1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}${\leq}\,1$\end{document} GV. No associated detectable increase was found in the space-borne data from the Geostationary Operational Environmental Satellite (GOES), the Energetic and Relativistic Nuclei and Electron (ERNE) on board the Solar and Heliospheric Observatory (SOHO), or the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) instruments, whose sensitivity was not sufficient to detect the event. No solar energetic particles were present during that time interval. The heliospheric conditions were slightly disturbed, so that the interplanetary magnetic field strength gradually increased during the event, followed by an increase of the solar wind speed after the event. It is proposed that the event was related to a crossing of the boundary layer between two regions with different heliospheric parameters, with a strong gradient of low-rigidity (<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}${<}\,1$\end{document} GV) particles. It was apparently similar to another cosmic-ray enhancement (e.g., on 22 June 2015) that is thought to have been caused by the local anisotropy of Forbush decreases, with the difference that in our case, the interplanetary disturbance was not observed at Earth, but passed by southward for this event.


Introduction
Variations in galactic cosmic-ray (GCR) flux near Earth are ultimately driven by solar variability on different timescales (Vainio et al., 2009). The temporal evolution of the cosmic-ray flux near Earth is traditionally and continuously monitored by the global network of neutron monitors (NMs), which is in continuous operation since the early 1950s (e.g. Simpson, 2000). The interannual variability of GCRs is driven by the 11-year solar cycle in the solar magnetic field that modulates them in the heliosphere (Potgieter, 2013). This is observed as a slow roughly 11-year variability of the NM count rate with a magnitude of about 25% at polar NMs (e.g. Usoskin, Bazilevskaya, and Kovaltsov, 2011), with alternating flat and sharp tops of the cycles representing the 22-year magnetic polarity cycle (Jokipii and Levy, 1977). On shorter timescales, the flux of GCRs is sporadically suppressed by interplanetary transient events (shocks, corotating regions, magnetic clouds, etc.), which cause reductions in GCR intensity, known as Forbush decreases (e.g. Lockwood, 1971;Cane, 2000;Ahluwalia and Fikani, 2007). In addition to GCR variability, sporadic enhancements of NM count rates are caused by solar energetic particle (SEP) events when a bulk of particles can be accelerated during solar eruptive events, flares, and/or coronal mass ejections (CMEs) to energies of hundreds of MeV up to several GeV (e.g. Desai and Giacalone, 2016;Klein and Dalla, 2017). Such events, when recorded by the ground-based NMs, are called ground-level enhancements (GLEs) and typically occur several times per solar cycle (see the International GLE database at http://gle.oulu.fi). Sometimes, the count rates of NMs can slightly increase as the result of a transient reduction of the geomagnetic shielding (Kudela, Bučík, and Bobík, 2008;Mohanty et al., 2016), but this is only related to detectors located at low or medium latitudes. Moreover, slight temporal enhancements in NM count rates can be observed during Forbush decreases (e.g. Cane, 2000). Cosmic-ray flux was enhanced on 22 June 2015, which was observed at high energies (several tens of GeV) by the global muon detector network during the active phase of a strong Forbush decrease and was speculated to be caused by a local anisotropic effect near the heliospheric current sheet in highly disturbed conditions (Munakata et al., 2018). Another apparently similar cosmic-ray enhancement was observed on 18 January 2017, whose origin is not known yet (Mangeard et al., 2017).
Here we report an event of highly anisotropic cosmic-ray enhancement (ACRE) that was first registered by the Dome C standard and bare NMs, conventionally denoted as DOMC and DOMB, respectively, located at the Central Antarctic Plateau Usoskin et al., 2015) during 7 June 2015, between roughly 12 and 19 UT. In contrast to the event on 22 June 2015 mentioned above, it was not accompanied by a Forbush decrease or other heliospheric disturbances recorded on Earth and only affected the lowest detectable energies below 1 GeV. This was a moderate increase in NM count rate that was not associated with a SEP flux, it was very anisotropic and lasted for several hours. The event was a subject of an informal discussion during the 34th International Cosmic Ray Conference in the Hague (2015), but has not been studied in detail since then. We present a detailed analysis of the event as recorded by the worldwide NM network and space-borne instruments, considering also the corresponding interplanetary conditions, and conclude that it was caused by focused scattering of GCRs on an interplanetary disturbance that passed Earth by.

Neutron Monitors
We have analyzed data of all available NMs for the period around the event, viz. 6 through 8 June 2015. Data were obtained from the European Neutron Monitor Data Base (NMDB: Figure 1 Five-minute averaged count rates of the NMs (see Table 1) for the period of 6 through 8 June 2015. The time interval corresponding to the event is hatched.
http://nmdb.eu, Mavromichalaki et al., 2011). Data from DOMC and DOMB are also available at http://cosmicrays.oulu.fi. These data are collected in the International GLE database (http://gle.oulu.fi). The time profiles of the count rates of the NMs, which depict an increase, are shown in Figure 1. A smooth increase in NM count rates between 12 and 19 UT on 7 June 2015 is visible, and the maximum is observed around 15 UT. The characteristics of the NMs and the recorded event magnitude are gathered in Table 1. Figure 2 shows the same data as in Figure 1, but normalized to the maximum count rate during the peak of the event (15:00 UT on 07 June 2015). The strongest response of 5 -8% was recorded by DOMB ("bare" or lead-free neutron monitor) and DOMC (standard NM), both located at the Condordia station on the Central Antarctic Plateau. Other Antarctic NMs, including South Pole, McMurdo, and Terre Adelie, recorded a statistically significant  increase of 2 -4%. Thus, the significant signal was detected only by polar southern stations. In addition, two high-latitude northern-hemisphere NMs, Nain and Peawanuck, detected a marginal increase during the event that was statistically significant. We have also examined data from other available NMs for the period around the event and found no increases.  Physical location of an NM can be misleading because as a result of the bending of charged-particle trajectories in the Earth's magnetosphere, the directions of acceptance of charged particles by an NM may be quite far from its actual location (Smart, Shea, and Flückiger, 2000). For all the NMs we calculated the asymptotic viewing directions for the time of the event maximum (15:00 UT on 7 June 2015) by the standard method of particletrajectory backtracing described elsewhere (e.g. Mishev, Poluianov, and Usoskin, 2017). The asymptotic directions were computed using the PLANETOCOSMICS code (Desorgher et al., 2005) and applying the International Geomagnetic Reference Field (IGRF: Finlay et al., 2010), epoch 2015, and Tsyganenko-89 (Tsyganenko, 1989 as internal and external geomagnetic field models, respectively. The asymptotic directions (in geographical coordinates) for protons with rigidities from 0.7 -5 GV are shown for some NMs in Figure 3   for the stations with a marginal response, and red crosses show stations with no response. It is interesting to observe that the asymptotic directions of the northern stations NAIN and PWNK are more south than that of the South Pole station. The significant response to the event is limited to a relatively small solid angle covering the range of about 60 -180 • west longitude to the south (we emphasize that these are GSE, not geographical coordinates). It is important that this area is not related to the apparent direction of the interplanetary magnetic field (IMF), which is denoted as the black crossed circle for the time of the event maximum (according to data from the Advance Composition Explorer (ACE) satellite). There is no response in any other region, implying that the event was highly anisotropic, with the anisotropy direction being in the southwest sector and not aligned with the IMF direction.

Space-Borne Data
We have also checked available space-borne data for the same period.
Near-Earth data from the OMNIWeb database (https://omniweb.gsfc.nasa.gov/ow.html) imply the absence of any SEP event (top panels of Figure 4). The proton fluxes do not show any significant enhancement above the background variability during the time of the event in either energy range (> 10, 30 or 60 MeV). There is a small enhancement of > 10 MeV protons during the day of the event, but it is not different from the similar variability during other days.
Data from the Energetic and Relativistic Nuclei and Electron (ERNE: Torsti et al., 1995) detector on board the Solar and Heliospheric Observatory (SOHO) spacecraft are somewhat gappy during the period of interest, but fully consistent with the "no-response" background for the day of 7 June 2015 (Vainio, 2018, private communication). Since the instrument is located at the L1 point outside the magnetosphere, there is no need for geomagnetic correction.
We have also performed a careful search for a possible signal at the PAMELA detector (Adriani et al., 2014) for the date of the event. The PAMELA detector for studying cosmic rays flew during the period from June 2006 through January 2016 on board the Resurs-DK1 satellite on a polar orbit, which had a height of about 570 km in 2015. Since the full analysis of the data from the PAMELA spectrometer in this low-energy range would not provide sufficient precision during a polar flyby, we used the count rate of its neutron detector as an integral index of cosmic-ray flux. The neutron detector consists of 18 pairs of helium-filled proportional counters (18.5 mm diameter and 209 mm length) surrounded by a polyethylene moderator and a cadmium layer that shield it from the slow neutrons that are produced in the satellite body. The detector permanently monitors the background neutron flux along the orbit, recording hits between triggers of the main spectrometer (Picozza et al., 2007). The count rate of the background neutron detector channel has a strong latitudinal dependence that varies between about 200 and 1000 Hz in equatorial and polar regions, respectively. The count rate is strongly enhanced in the region of the South Atlantic Anomaly (SAA) because of precipitation of the trapped particles from the radiation belt. The orbit moves longitudinally by about 26 • at each rotation, leading to a variation of the count rates that is related to the tilt of the geomagnetic field. Using the known dependence of the neutron-detector count rate on the geomagnetic cutoff rigidity, we have estimated the expected count rate along the orbit, which works well outside the SAA region. We compared the actually recorded values with the expected count rates for the day of 7 June 2015 and found no statistically significant difference. This implies that no detectable excess of secondary neutrons has been found that would correspond to the ACRE. Thus, data from the neutron counter of the PAMELA instrument are also consistent with a "no-response" background.
However, since PAMELA is a low-orbiting detector that is always inside the geomagnetosphere, an analysis of the asymptotic directions should be performed. The asymptotic directions of PAMELA are shown in Figure 3 for the southernmost point of its orbit at 14:34 UT on 7 June 2015, near the peak of the event. The asymptotic direction lies outside the event region as identified with NMs, thus confirming the anisotropic nature of the event. We note that the sensitivity of the PAMELA neutron counter to this low-energy particles is about 1% above the background, and a signal at a level of several percentages, as observed by Antarctic NMs, would have been statistically detected there if the detector had crossed the event region.

Interplanetary Conditions
The near-Earth interplanetary conditions during the event are shown in Figure 4: the integral flux of protons with energies above 10, 30, and 60 MeV (top panel), plasma flow speed (central panel), and the IMF magnitude (bottom panel). The bottom panel of Figure 4 shows that the IMF strength started to grow from 5 nT shortly before the event, at about 07:15 UT, until about noon of the next day, when it reached about 20 nT and then suddenly dropped. This gradual increase was accompanied by a slow increase in solar wind speed from about 350 km s −1 before the event to about 690 km s −1 at noon of 8 June. This may indicate the presence of a moderate compression region that formed upstream of a high-speed wind stream. However, no severe change, such as an interplanetary shock, was observed. Accordingly, the heliospheric conditions were not greatly disturbed during the time of the event.
However, a large view shows that there was indeed a strong heliospheric disturbance that did not hit Earth, however, and thus was not measured in the near-Earth space in the ecliptic plane. A reconstruction (obtained from http://helioweather.net), using the ENLIL Model (Odstrčil, Smith, and Dryer, 1996) of the dynamical 3D heliospheric conditions is shown in Figure 5 for 14 UT of the day of the event. While the ecliptic plane was not much disturbed, there was a fast "ejecta" with a velocity of 600 -700 km s −1 . It passed south and west of the Earth, so that our planet was exactly at the boundary of the ejecta during the event.

Discussion and Speculation
We can summarize the experimental facts and speculate about the event as follows.
There was an increase in count rates of polar NMs, whose asymptotic directions lie in the southwest quadrant in GSE coordinates. The increase was strongest for the high-altitude nocutoff DOMB/DOMC station, moderate for the high-altitude low-cutoff SOPO/SOPB and the sea-level no-cutoff MCMD stations, and marginal for the low-cutoff sea-level stations NAIN and PWNK that are located in the northern hemisphere. This suggests that the effect was related to low-energy particles with a rigidity below 1 GV. No associated detectable increase was found in the space-borne data from the GOES, SOHO/ERNE, or PAMELA instruments. There are some hints, but they are not statistically significant. No SEP event occurred during the period in question.
The event took place on the background of relatively quiet near-Earth interplanetary conditions, without any strong disturbance as recorded by spacecraft in the ecliptic plane during the event. There was no Forbush decrease, in contrast to the event of 22 June 2015, nor any geomagnetic storm. However, a fast ejecta passed south and west of Earth, so that the planet was exactly at the boundary shear layer of the ejecta.
Based on an analysis of available data, we can speculate that the event was related to the fact that Earth brushed against a boundary layer between the high-speed stream and the slow solar wind region ( Figure 5). The IMF strength and solar wind speed both gradually increased during the ACRE, corresponding to a compression region that was very asymmetric with respect to the ecliptic plane. This led to a strong local anisotropy of low-rigidity (< 1 GV) cosmic-ray flux. This event was similar to that of 22 June 2015 (and possibly also to the event on 18 January 2017), with the difference that the interplanetary disturbance did not hit Earth, but passed nearby. A more systematic study of this type of events is pending. We note that another apparently similar ACRE took place on 26 August 2018, after submission of this manuscript, and its full analysis is in progress.