Near-Earth Cosmic Ray Decreases Associated with Remote Coronal Mass Ejections
Total Page:16
File Type:pdf, Size:1020Kb
The Astrophysical Journal, 801:5 (8pp), 2015 March 1 doi:10.1088/0004-637X/801/1/5 C 2015. The American Astronomical Society. All rights reserved. NEAR-EARTH COSMIC RAY DECREASES ASSOCIATED WITH REMOTE CORONAL MASS EJECTIONS S. R. Thomas1,2,M.J.Owens2, M. Lockwood2, L. Barnard2, and C. J. Scott2 1 Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Surrey, RH5 6NT, UK 2 Department of Meteorology, University of Reading, Reading, RG6 6BB, UK; [email protected] Received 2014 October 13; accepted 2014 December 25; published 2015 February 23 ABSTRACT Galactic cosmic ray (GCR) flux is modulated by both particle drift patterns and solar wind structures on a range of timescales. Over solar cycles, GCR flux varies as a function of the total open solar magnetic flux and the latitudinal extent of the heliospheric current sheet. Over hours, drops of a few percent in near-Earth GCR flux (Forbush decreases, FDs) are well known to be associated with the near-Earth passage of solar wind structures resulting from corotating interaction regions (CIRs) and transient coronal mass ejections (CMEs). We report on four FDs seen at ground-based neutron monitors which cannot be immediately associated with significant structures in the local solar wind. Similarly, there are significant near-Earth structures which do not produce any corresponding GCR variation. Three of the FDs are during the STEREO era, enabling in situ and remote observations from three well-separated heliospheric locations. Extremely large CMEs passed the STEREO-A spacecraft, which was behind the West limb of the Sun, approximately 2–3 days before each near-Earth FD. Solar wind simulations suggest that the CMEs combined with pre-existing CIRs, enhancing the pre-existing barriers to GCR propagation. Thus these observations provide strong evidence for the modulation of GCR flux by remote solar wind structures. Key words: cosmic rays – Sun: coronal mass ejections (CMEs) – Sun: heliosphere – Sun: magnetic fields 1. INTRODUCTION commonly referred to as a “Forbush decrease” (FD). Due to the recent “weak” solar cycle (e.g., Lockwood et al. 2012), there When a galactic cosmic ray (GCR) proton enters the atmo- have been fewer fast/large CMEs than is normally expected for sphere, it interacts with atmospheric particles and produces a this phase of the solar cycle (Webb & Howard 2012), and thus cascade of secondary particles, such as neutrons. Some of these we can expect a smaller number of FDs in this period. neutrons with high enough energies are able to penetrate through A number of studies have investigated how FDs vary with to the surface can be detected by neutron monitors, located at a the properties of the driving solar wind structure (e.g., Blanco number of stations at a variety of latitudes and longitudes. Thus et al. 2013 and references therein). However, not all FDs are neutron count rates serve as an indirect proxy for the top-of-the- associated with an obvious signature of a CME or CIR seen atmosphere GCR flux, once the dependence of the atmospheric in near-Earth space. Indeed, Cane et al. (1993) found FDs that absorption of neutrons on the atmospheric pressure is consid- appear to have arisen from remote connections with shocks, ered (e.g., Raubenheimer & Stoker 1974). The standard pressure primarily driven by CMEs propagating from toward the east correction to neutron monitor counts is applied to all neutron limb of the Sun. This study provides strong evidence for such monitor stations used in this study. modulation of GCRs by remote solar wind structures. In Sec- The flux of GCRs is modulated within the heliosphere on a tion 2 we examine a particularly large example of an FD with- wide range of timescales. On decadal timescales, GCRs are out a significant structure in near-Earth space. In Section 3 we influenced by differing energetic particle drift patterns and present further events from a period with greater spatial sam- long-term variations in the heliospheric magnetic field (HMF) pling of the heliosphere from the STEREO spacecraft. Finally, in magnitude and direction (e.g., Jokipii et al. 1977; Thomas et al. Section 4 we use solar wind simulations to provide interpreta- 2014a). Over timescales of hours to days, solar wind structures tion of these observations in terms of the global structure of the associated with coronal mass ejections (CMEs; Cane 2000)or inner heliosphere. corotating interaction regions (CIRs; e.g., El Borie et al. 1998; 2. A “PHANTOM” FORBUSH DECREASE Thomas et al. 2014b) affect GCR flux at Earth. The passage of a CME past an observing spacecraft is frequently associated Figure 1 shows a large FD, observed in neutron monitor with an enhancement in the magnetic field magnitude and a count rates on the 2002 November 16. The top row shows variation in the magnetic field direction away from typical the percentage change in neutron count rates relative to the Parker spiral orientations (Lepping et al. 1990). If the speed 20 day mean for McMurdo (latitude 77◦.9 south; longitude of the CME exceeds that of the surrounding solar wind by 166◦.6 east; rigidity cut-off 0.3 GV), Newark (latitude 39◦.7 north; greater than the fast magnetosonic wave speed, it will drive a longitude 75◦.7 west; rigidity cut-off 2.4 GV), Oulu (latitude shock front, compressing and heating the ambient solar wind 65◦.1 north; longitude 25◦.5 east; rigidity cut-off 0.81 GV), and ahead of it in a “sheath” region bounded by the shock and CME Alma-Ata (latitude 43◦.1 north; longitude 76◦.6 east; rigidity leading edge. As GCRs can be scattered by inhomogeneities cut-off 6.7 GV) in red, blue, green, and orange respectively. and discontinuities in the HMF (e.g., Giaclone & Jokipii 1999), The estimated rigidity cut-off values are accessed from the the GCR flux at Earth is strongly modulated by the arrival of an respective neutron monitor’s Web site, e.g., the Bartol Institute Earth-impacting CME, in particular fast CMEs. This can result website for McMurdo and Newark. However, the rigidity cut- in GCR flux reducing suddenly over a number of hours, before off can be influenced substantially by geomagnetic storms and it recovers to pre-CME levels over a number of days (Forbush planetary indices such as the Kp index must be considered 1937;Forbush1954). Such a temporal variation in GCR flux is (Smart et al. 2006). 1 The Astrophysical Journal, 801:5 (8pp), 2015 March 1 Thomas et al. Figure 1. Forbush decrease on 2002 November 17. The top panel shows neutron monitor counts from McMurdo (red), Newark (blue), Oulu (green), and Alma-Ata (orange). The second panel shows the hourly variability in the near-Earth HMF, δB . The near-Earth magnetic field strength, |B|, radial solar wind speed, |Vx|,and the proton density, np, are displayed in the next three panels. The green horizontal line give the 1σ ,2σ ,and3σ occurrences of enhancements of each variable from OMNI2 with respect to the 2001–2003 mean. During the period shown in Figure 1, each of the neutron potentially GCR-modulating structures in the solar wind, we monitors show very similar variations, despite the large spread analyzed the probability distribution functions of all solar wind of rigidity cut-offs of the stations. There is a slight spread parameters, using data from 2001 through 2003. From these in the onset timings of the FD due to the longitudes of the probability density functions, we calculate the threshold values stations essentially providing different viewing angles into the corresponding to 1σ ,2σ , and 3σ above the mean, which are heliosphere. Also, it is found from NOAA data that at the time of shown as horizontal green lines in Figures 1–3. the FD the Kp index did not increase above a value of 3 ruling out For each station, the neutron count rates reduce from ap- influence on neutron monitor counts from geomagnetic activity. proximately 2%–3% above the background to 4%–5% below The second panel shows a measure of the variability in the it. This represents one of the largest FDs within the neutron HMF (δB ) calculated from the hourly OMNI-2 data set (King monitor data series. Such a decrease in neutron count rates is & Papitashvili 2005). Here δB is the modulus of the differences expected to be due to a sizable GCR barrier in near-Earth space, between consecutive HMF vectors, i.e., δB =|Bi − B i −1|. typically resulting from a large, fast Earth-directed CME or a The third panel shows HMF magnitude (|B|), the fourth the particularly strong Earth-impacting CIR. Although a CME in radial solar wind speed (|Vx|), and the final shows the proton near-Earth space identified at this time (Richardson & Cane density (np). In order to objectively assess the presence of 2010) primarily on the basis of a reduced proton temperature, 2 The Astrophysical Journal, 801:5 (8pp), 2015 March 1 Thomas et al. Figure 2. Forbush decrease on 2012 May 30. The top panel shows neutron monitor counts from McMurdo (red), Newark (blue), Oulu (green), and Alma-Ata (orange). The second shows the hourly magnetic field variability, δB , in near-Earth space (black), at STEREO-A (red) and STEREO-B (blue). STEREO data are scaled so that the probability density functions match OMNI2 data, to allow direct comparisons of significance. The remaining three panels show (unscaled) magnetic field strength, |B|, radial solar wind speed, |Vx| and the proton density, np. The green horizontal lines give the 1σ ,2σ ,and3σ from OMNI2 with respect to the 2011–2013 mean. The schematic in the bottom-right gives the approximate locations of the measurement spacecraft.