The 1990 May 24 solar cosmic-ray event

This paper presents an integrated analysis of GOES 6, 7 and neutron monitor observations of solar cosmic-ray event following the 1990 May 24 solar flare. We have used a model which includes particle injection at the Sun and at the interplanetary shock front and particle propagation through the interplanetary medium. The model does not attempt to simulate the physical processes of coronal transport and shock acceleration, therefore the injections at the Sun and at the shock are represented by source functions in the particle transport equation. By fitting anisotropy and angle-average intensity profiles of high-energy (>30 MeV) protons as derived from the model to the ones observed by neutron monitors and at GOES 6 and 7, we have determined the parameters of particle transport, the injection rate and spectrum at the source. We have made a direct fit of uncorrected GOES data with both primary and secondary proton channels taken into account.The 1990 May 24–26 energetic proton event had a double-peaked temporal structure at energies 100 MeV. The Moreton (shock) wave nearby the flare core was seen clearly before the first injection of accelerated particles into the interplanetary medium. Some (correlated with this shock) acceleration mechanism which operates in the solar corona at a height up to one solar radius is regarded as a source of the first (prompt) increase in GOES and neutron monitor counting rates. The proton injection spectrum during this increase is found to be hard (spectral index 1.6) at lower energies ( 30 MeV) with a rapid steepening above 300 MeV. Large values of the mean free path ( 1.8 AU for 1 GV protons in the vicinity of the Earth) led to a high anisotropy of arriving protons. The second (delayed) proton increase was presumably produced by acceleration/injection of particles by an interplanetary shock wave at height of 10 solar radii. Our analysis of the 1990 May 24–26 event is in favour of the general idea that a number of components of energetic particles may be produced while the flare process develops towards larger spatial/temporal scales.Visiting Associate from St. Petersburg State Technical University, St. Petersburg 195251, Russia.     

Atmospheric refraction of solar neutrons during the event of 24 May 1990

We revise the published neutron monitor raw data for the increase caused by the solar neutron event of the 24 May 1990. With these data we calculate the attenuation length, , of solar neutrons in the Earth's atmosphere assuming either a minimum path as given by the spread of elastically scattered neutrons, or using the minimum mass path estimated by Smart, Shea, and O'Bren (1995) due to an atmospheric refraction effect. In both cases reduces to a value around 100 g cm^–2, which is more in accordance with data on neutron cross-sections (Shibata, 1994). These two phenomenological calculations suggest that solar neutrons do not propagate in straight lines in the atmosphere. The previous estimate of the attenuation length, =208 g cm^–2, was calculated assuming straight-ahead transport (Smart, Shea, and O'Bren, 1995). Dorman, Valdes-Galicia, and Dorman (1999) performed a numerical simulation and an analytical approximation to the problem of solar neutron scattering and attenuation in the Earth's atmosphere. These solutions incorporate the refraction effect as a natural consequence of the greater absorption experienced by neutrons scattered to large zenith angles. They are able to reproduce the normalised observed counting rates of neutron monitors for this event.     

Cosmic-Ray Variations During the Two Greatest Bursts of Solar Activity in the 23rd Solar Cycle

During two extreme bursts of solar activity in March–April 2001 and October–November 2003, the ground-based neutron monitor network recorded a series of outstanding events distinguished by their magnitude and unusual peculiarities. The important changes that lead to increased activity initiated not with the sunspot appearance, but with the large-scale solar magnetic field reconfiguration. A series of strong and moderate magnetic storms and powerful proton events (including ground-level enhancements, GLE) were registered during these periods. The largest and most productive in the 23rd solar cycle, active region 486, generated a significant series of solar flares among which the 4 November 2003 flare (X28/3B) was the most powerful X-ray solar event ever observed. The fastest arrival of the interplanetary disturbance from the Sun (after August 1972) and the highest solar wind velocity and IMF intensity were recorded during these events. Within 1 week, three GLEs of solar cosmic rays were registered by the neutron monitor network (28 and 29 October and 2 November 2003). In this work, we perform a tentative analysis of a number of the effects seen in cosmic rays during these two periods, using the neutron monitor network and other relevant data.     

The 1991 March 22 flare: Possible anisotropy of high-energy neutral emission

We made a parameter fit to the Haleakala neutron monitor counting rate during the 1991 March 22 solar flare (Pyle and Simpson, 1991) using the time profiles of -rays at 0.42–80 MeV obtained with the GRANAT satellite (Vilmeret al., 1994) and the microwave data from Owens Valley Radio Observatory. We use a two-component neutron injection function to find that either an impulsive injection or the impulsive-plus-prolonged neutron injection is possible. In both cases, the number of > 300 MeV neutrons emitted towards the Earth is estimated as 2 × 10²⁷ sr^–1, which is less than that of the 1990 May 24 flare by an order of magnitude.We tested if such a big difference in neutron number detected on the Earth can be accounted for solely by their different positions on the solar disk. For the estimation of the degree of anisotropy of high-energy secondary emission, we made use of macroscopic parameters of the flare active region, in particular, the vector magnetogram data from the Big Bear Solar Observatory. In our result, the anisotropy factor for the neutral emissions of the 1991 March 22 flare is only 1 – 10, which is rather small compared with previous theoretical predictions for a disk flare. Such a moderate anisotropy is due to the relatively large inclination angles of the magnetic fields at the footpoints of the flaring loop where accelerated particles are trapped. We thus concluded that the smaller number of neutrons of the 1991 March 22 flare would be not only due to its location on the disk, but also due to fewer protons accelerated during this event as compared with the 1990 May 24 limb event. For a more precise determination of the anisotropy factor in a flare, we need a detailed spectrum of electron bremsstrahlung in 0.1 – 10 MeV and the fluence of -ray emission from the ⁰-decay.Visting Associate from St. Petersburg State Technical University, St. Petersburg, 195251, Russia.     

A joint analysis of high-energy neutrons and neutron-decay protons from a flare

A joint analysis of neutron monitor and GOES data is performed to study the production of high-energy neutrons at the Sun. The main objects of the research are the spectrum of >50 MeV neutrons and a possible spectrum of primary (interacting) protons which produced those neutrons during the major 1990 May 24 solar flare. Different possible scenarios of the neutron production are presented. The high magnitude of the 1990 May 24 neutron event provided an opportunity to detect neutron decay protons of higher energies than ever before. We compare predictions of the proposed models of neutron production with the observations of protons on board GOES 6 and 7. It is shown that the precursor in high-energy GOES channels observed during 20:55–21:09 UT can be naturally explained as originating from decay of neutrons in the interplanetary medium. The ratio of counting rates observed in different GOES channels can ensure the selection of the model parameters.The set of experimental data can be explained in the framework of a scenario which assumes the existence of two components of interacting protons in the flare. A hard spectrum component (the first component) generates neutrons during a short time while the interaction of the second (soft spectrum) component lasts longer. Alternative scenarios are found to be of lesser likelihood. The intensity-time profile of neutron - decay protons as predicted in the framework of the two-component exponential model of neutron production (Kocharov et al., 1994a) is in an agreement with the proton profiles observed on board GOES. We compare the deduced characteristics of interacting high-energy protons with the characteristics of protons escaping into the interplanetary medium. It is shown that, in the 100–1000 MeV range, the spectrum of the second component of interacting protons was close to the spectrum of the prompt component of interplanetary protons. However, it is most likely that, at 300 MeV, the interacting proton spectrum was slightly softer than the spectrum of interplanetary protons. An analysis of gamma-ray emission is required to deduce the spectrum of interacting protons below 100 MeV and above 1 GeV.     

Neutron and electromagnetic emissions during the 1990 May 24 solar flare

In this paper, we are primarily concerned with the solar neutron emission during the 1990 May 24 flare, utilizing the counting rate of the Climax neutron monitor and the time profiles of hard X-rays and γ-rays obtained with the GRANAT satellite (Pelaezet al., 1992; Talonet al., 1993; Terekhovet al., 1993). We compare the derived neutron injection function with macroscopic parameters of the flare region as obtained from theHα and microwave observations made at the Big Bear Solar Observatory and the Owens Valley Radio Observatory, respectively. Our results are summarized as follows: (1) to explain the neutron monitor counting rate and 57.5–110 MeV and 2.2 MeV γ-ray time profiles, we consider a two-component neutron injection function,Q(E, t), with the form $$Q(E,t) = N_f {\text{ exp[}} - E/E_f - t/T_f ] + N_s {\text{ exp[}} - E/E_s - t/T_s ],$$ whereN _f(s),E _f(s), andT _f(s) denote number, energy, and decay time of the fast (slow) injection component, respectively. By comparing the calculated neutron counting rate with the observations from the Climax neutron monitor we derive the best-fit parameters asT _f ≈ 20 s,E _f ≈ 310 MeV,T _s ≈ 260 s,E _s ≈ 80 MeV, andN _f (E > 100 MeV)/N _s (E > 100 MeV) ≈ 0.2. (2) From the Hα observations, we find a relatively small loop of length ≈ 2 × 10⁴ km, which may be regarded as the source for the fast-decaying component of γ-rays (57.5–110 MeV) and for the fast component of neutron emission. From microwave visibility and the microwave total power spectrum we postulate the presence of a rather big loop (≈ 2 × 10⁵ km), which we regard as being responsible for the slow-decaying component of the high-energy emission. We show how the neutron and γ-ray emission data can be explained in terms of the macroscopic parameters derived from the Hα and microwave observations. (3) The Hα observations also reveal the presence of a fast mode MHD shock (the Moreton wave) which precedes the microwave peak by 20–30 s and the peak of γ-ray intensity by 40–50 s. From this relative timing and the single-pulsed time profiles of both radiations, we can attribute the whole event as due to a prompt acceleration of both electrons and protons by the shock and subsequent deceleration of the trapped particles while they propagate inside the magnetic loops.     

Effective Rigidity of a Polar Neutron Monitor for Recording Ground-Level Enhancements

The “effective” rigidity of a neutron monitor for a ground-level enhancement (GLE) event is defined so that the event-integrated fluence of solar energetic protons with rigidity above it is directly proportional to the integral intensity of the GLE as recorded by a polar neutron monitor, within a wide range of solar energetic-proton spectra. This provides a direct way to assess the integral fluence of a GLE event based solely on neutron-monitor data. The effective rigidity/energy was found to be 1.13?–?1.42 GV (550?–?800 MeV). A small model-dependent, systematic uncertainty in the value of the effective rigidity is caused by uncertainties in the low-energy range of the neutron-monitor yield function, which requires more detailed computations of the latter.     

A search for solar neutrons during solar flares

The upper limit on the solar neutron flux from 1–20 MeV has been measured, by a neutron detector on the OGO-6 satellite, to be less than 5 × 10^–2 n cm^–2 s^–1 at the 95% confidence level for several flares including two flares of importance 3B and a solar proton event of importance 3B. The measurements are consistent with the models proposed by Lingenfelter (1969) and by Lingenfelter and Ramaty (1967) for solar neutron production during solar flares. The implied upper limit on the flux of 2.2 MeV solar gamma rays is about the same as the 2.2 MeV flux observed by Chupp et al. (1973).     

Solar Energetic Particle Events with Protons Above 500 MeV Between 1995 and 2015 Measured with SOHO/EPHIN

The Sun is an effective particle accelerator that produces solar energetic particle (SEP) events, during which particles of up to several GeVs can be observed. These events, when they are observed at Earth with the neutron monitor network, are called ground-level enhancements (GLEs). Although these events with their high-energy component have been investigated for several decades, a clear relation between the spectral shape of the SEPs outside the Earth’s magnetosphere and the increase in neutron monitor count rate has yet to be established. Hence, an analysis of these events is of interest for the space weather and for the solar event community.In this article, SEP events with protons accelerated to above 500 MeV were identified using data obtained with the Electron Proton Helium Instrument (EPHIN) onboard the Solar and Heliospheric Observatory (SOHO) between 1995 and 2015. For a statistical analysis, onset times were determined for the events and the proton energy spectra were derived and fitted with a power law.As a result, we present a list of 42 SEP events with protons accelerated to above 500 MeV measured with the EPHIN instrument onboard SOHO. The statistical analysis based on the fitted spectral slopes and absolute intensities is discussed, with special emphasis on whether an event has been observed as a GLE. Furthermore, we are able to determine that the derived intensity at 500 MeV and the observed increase in neutron monitor count rate are correlated for a subset of events.     

On forecasting solar proton events using a ground-based neutron monitor

Results of applying the method developed for early warning about arrival of ∼10–100-MeV proton fluxes to the Earth after powerful eruptive events on the Sun, which uses the real time observable data received by the global network of ground-based neutron monitors (Mavromichalaki et al., 2009), are discussed. The retrospective analysis and comparison to the 2001–2006 observations indicate that more than 50% of solar proton events were omitted in such a forecasting method. For higher reliability, it is necessary to use additional data on the state of solar and heliospheric activity.     

Particle trapping and acceleration during the August 1972 event

Several features of the August 1972 events are studied using neutron monitor data together with solar wind streamlines calculated on the basis of an approximate kinematic approach. Examination of the evolution of these streamlines shows that the streamline which passes from the Earth undergoes dramatic changes during the main phase of these events. In a few hours this streamline, which was estimated with HEOS-2 solar wind velocity data, was decreased (i.e., compressed) to a total radial extend of 0.2 AU (at the beginning of 5 August), although its initial length was 1 AU. An exact MHD time-dependent solution by Dryer et al. (1978a) gives similar results.The relative cosmic ray increase (3–7 UT, 5 August), immediately after the deep F.d., is attributed to trapping and acceleration of particles between two shock waves. Similar acceleration was found by Pomerantz and Duggal (1974) and Levy et al. (1976) for another cosmic ray increase during this event.The extremely large solar wind velocities during the main phase of the event are not only due to the large energy of the flare but also to the fact that the ambient solar wind was already almost empty because of the sweeping action of previous shock waves.Work performed partly at the University of Athens and partly at Imperial College of Science and Technology.     

THE RELATIONSHIP BETWEEN SOLAR FLARE GAMMA-RAY EMISSION AND NEUTRON PRODUCTION

We have examined six solar neutron events measured by satellite instruments and/or neutron monitors (NM) to understand the relationship between the intensity–time profiles of the -ray lines, the pion-related -rays, and the neutron production. In all six events the solar neutron production was clearly time-extended. We find that neutron emission as detected by NMs most closely follows the emission of pion-related -rays, whereas lower energy neutron production may follow that of nuclear -ray line emissions. Although this distinction is not unexpected, it is safe to say that the 2.223 MeV -ray line from neutron capture on hydrogen is a poor measure of the neutron production at energies >200 MeV. During the three events on 1982, June 3, 1990, May 24 and 1991, June 4 solar neutrons with energies greater than 200 MeV were recorded by NMs. The NM increases on 1982, June 3 and 1990, May 24 can be modeled using the time profile of the pion-related -rays. For the 1991, June 4 event the NM signal was small but lasted for 60 min and the high-energy -ray data available to us are insufficient to conclude unambiguously that the high-energy neutron production followed the pion-related -rays. In the other three events on 1991, June 9, 11, and 15 solar neutrons with energies 10–100 MeV were observed by the COMPTEL -ray instrument on the Compton Gamma Ray Observatory. The duration of the low-energy neutron production on 1991, June 9 corresponded clearly to the high-energy and not to the low-energy -ray emission.     

Non-diffusive particle pulse transport – Application to an anisotropic solar GLE

A theory of the transport of an anisotropic pulse of charged particles injected into the interplanetary magnetic field is applied to an anisotropic ground level event on 24 May 1990. For this event the kinetic regime is considered when the mean free path is comparable with the distance from particle source. Both the source angular particle distribution and the angular dependence of a detector response are included. The theoretically predicted temporal profiles are compared with the particle intensity records measured by several neutron monitors with different asymptotic directions.     

First Analysis of Ground-Level Enhancement (GLE) 72 on 10 September 2017: Spectral and Anisotropy Characteristics

Using data obtained with neutron monitors and space-borne instruments, we analyzed the second ground-level enhancement (GLE) of Solar Cycle 24, namely the event of 10 September 2017 (GLE 72), and derived the spectral and angular characteristics of associated GLE particles. We employed a new neutron-monitor yield function and a recently proposed model based on an optimization procedure. The method consists of simulating particle propagation in a model magnetosphere in order to derive the cutoff rigidity and neutron-monitor asymptotic directions. Subsequently, the rigidity spectrum and anisotropy of GLE particles are obtained in their dynamical evolution during the event on the basis of an inverse-problem solution. The derived angular distribution and spectra are discussed briefly.     

A Yohkoh search for “black-light flares”

Calculations which predict that a phenomenon analogous to stellar negative pre-flares could also exist on the Sun were published by Hénouxet al. (1990), and Aboudarhamet al., (1990), who showed that at the beginning of a solar white-light flare (WLF) event an electron beam can cause a transient darkening before the WLF emission starts, under certain conditions. They named this event a black light flare (BLF). Such a BLF event should appear as diffuse dark patches lasting for about 20 seconds preceding the WLF emission, which would coincide with intense and impulsive hard X-ray bursts. The BLF location would be at (or in the vicinity of) the forthcoming bright patches. Their predicted contrast depends on the position of the flare on the solar disc and on the wavelength band of the observation.TheYohkoh satellite provided white-light data from the aspect camera of the SXT instrument (Tsunetaet al., 1991), at 431 nm and with a typical image interval of 10–12 s. We have studied nine white-light flares observed with this instrument, with X-ray class larger than M6. We have found a few interesting episodes, but no unambiguous example of the predicted BLF event. This study, although the best survey to date, was not ideal from the observational point of view. We therefore encourage further searches. Successful observations of this phenomenon on the Sun would greatly strengthen our knowledge of the lower solar atmosphere and its effects on solar luminosity variations.     

Impacts on Cosmic-Ray Intensity Observed During Geomagnetic Disturbances

Geomagnetic disturbances are the results of interplanetary causes such as high-speed streamers (HSSs), interplanetary coronal mass ejections (ICMEs), corotating interaction regions (CIRs), and magnetic clouds. During different forms of geomagnetic disturbances, we observed changes in the count rate at neutron monitors that are kept at various locations. We studied the count rates measured by neutron monitors at four stations at various latitudes during different categories of geomagnetic events and compared them. We analysed five events: a geomagnetically quiet event, a non-storm high-intensity long-duration continuous AE activity (HILDCAA) event, a storm-preceded HILDCAA event, a geomagnetic substorm event, and a geomagnetic moderate storm event. We based our analysis on geomagnetic indices, solar wind parameters, and interplanetary magnetic field (IMF) parameters. We found that the strength of the modulation was least during the quiet event and highest during the storm-preceded HILDCAA. By analysing the cause of these geomagnetic disturbances, we related each decrease in the neutron monitor data with the corresponding solar cause. For the ICME-driven storm, we observed a decrease in neutron monitor data ranging from 6% to 12% in all stations. On the other hand, we observed a decrease ranging from 2% to 5% for the HSS-driven storm. For the non-storm HILDCAA, we observed a decrease in neutron monitor data of about 1% to 1.5%. For the quiet event, the neutron monitor data fluctuated such that there was no overall decrease in all stations.     

Interplanetary disturbances produced by a simulated solar flare and equatorially-fluctuating heliospheric current sheet

A recently developed nonplanar, time-dependent magnetohydrodynamic (MHD) model (Wuet al., 1983) was used to study the interplanetary disturbances produced by a compound event in the heliosphere. That is, a steady-state interplanetary medium is first disturbed by a simulated equatorially-fluctuating current sheet. After a few days (100 hr), the disturbed interplanetary medium is again perturbed by a solar-flare-generated shock wave. Attention is directed toward the differences that are caused by the presence of the equatorially-fluctuating (warped) current sheet.     

The OSO-1 solar neutron experiment

BF₃ counters on OSO-1 were used to look for solar neutrons by trying to observe a diurnal variation in count rate. No effect was observed and an upper limit was placed on the solar neutron flux at the earth of J _n < 2 × 10^–3 neut/cm²/sec of 10 keV < E _n < 10 meV for the period March to May 1962. No proton-producing flares occurred during this time, so the most obvious source of solar neutrons could not be studied. The emulsion experiment of Apparao et al., flown during this period, which seemed to indicate solar neutrons has probably been misinterpreted.     

Phase differences between irradiance and velocity low degree solar acoustic modes revisited

Since 1984, simultaneous observations of irradiance and velocity solar acoustic modes, have been carried out by several authors in order to measure the phase difference between irradiance and velocity modes. Following the earliest observations with stratospheric balloon (Frolich and van Der Raay, 1984), a two ground-based stations (Tenerife and Baja California) were established (Jimenez et al, 1990) obtaining coherence results in the frequency range from 2.5 mHz to 4.3 mHz. These phase differences between irradiance and velocity solar acoustic modes are interpreted in terms of the non-adiabatic behaviour of the solar atmosphere. In 1988 the IPHIR (Frolich et al, 1988) instrument flown on the PHOBOS-2 mission to Mars and measured the solar irradiance during 150 consecutive days. The best velocity observations obtained in Tenerife for this period were compared with IPHIR data to compute the phase differences (Schrijver et al, 1991). The final conclusion is that good agreement is attained between space quadsi-space and ground observations which yield a phase diffenrece of about -125 degrees in the frequency range 2.5 mHz to 4.2 mHz, with a slight increase suggested by the data running up to 4.6 mHz.     

Effect of solar heliospheric parameters on different components of daily variation in cosmic rays

In the present work the data of three different neutron monitoring stations, Deep River, Tokyo and Inuvik located at different geomagnetic cutoff rigidities and altitudes has been harmonically analysed for the period 1980–1993, 1980–1990 and 1981–1993 respectively to investigate for a comparative study of diurnal, semi-diurnal and tri-diurnal anisotropies in cosmic ray (CR) intensity in connection with the change in IMF Bz component and solar wind velocity on 60 quietest days. It is observed that the amplitude of first harmonic is highly anti-correlated to the solar wind velocity during the period of high-speed solar wind stream (HSSWS) epoch on quiet days for three neutron monitor stations at different geomagnetic rigidity thresholds. During quiet days the amplitude of all the three harmonics significantly deviates on the onset of HSSWS epoch, whereas the direction of the anisotropy of all the three harmonics remains time invariant at three different cut off rigidity stations. The amplitude as well as the direction of anisotropy of all the three harmonics does not have time variation characteristics associated with Bz component of IMF on geo-magnetically most quiet days.