Occurrence Probability of Large Solar Energetic Particle Events: Assessment from Data on Cosmogenic Radionuclides in Lunar Rocks

We revisited assessments of the occurrence probability distribution of large events in solar energetic particles (SEP), based on measurements of cosmogenic radionuclides in lunar rocks. We present a combined cumulative occurrence probability distribution of SEP events based on three timescales: directly measured SEP fluences for the past 60 years; estimates based on the terrestrial cosmogenic radionuclides ¹⁰Be and ¹⁴C for the multi-millennial (Holocene) timescale; and cosmogenic radionuclides measured in lunar rocks on a timescale of up to 1 Myr. These three timescales yield a consistent distribution. The data suggest a strong roll-over of the occurrence probability, so that SEP events with a proton fluence with energy >?30 MeV greater than 10¹¹ (protons cm^?2?yr^?1) are not expected on a Myr timescale.     

Nitrate in Polar Ice: A New Tracer of Solar Variability

Knowledge of the long-term variability of solar activity is of both astrophysical and geoscientific interest. Reconstructions of solar activity over multiple millennia are traditionally based on cosmogenic isotopes ¹⁴C or ¹⁰Be measured in natural terrestrial archives, but the two isotopes exhibit significant differences on millennial time scales, so that our knowledge of solar activity at this time scale remains somewhat uncertain. Here we present a new potential proxy of solar activity on the centennial-millennial time scale, based on a chemical tracer, viz. nitrate content in an ice core drilled at Talos Dome (Antarctica). We argue that this location is optimal for preserving the solar signal in the nitrate content during the Holocene. By using the firn core from the same location we show that the 11-year and Gleissberg cycles are present with the variability of 10??C?25?% in nitrate content in the pre-industrial epoch. This is consistent with the results of independent efforts of modeling HNO₃ and NO _y in Antarctic near surface air. However, meteorological noise on the interannual scale makes it impossible to resolve individual solar cycles. Based on different processes of formation and transport compared to cosmogenic isotopes, it provides new, independent insight into long-term solar activity and helps resolve the uncertainties related to cosmogenic isotopes as diagnostics of solar activity.     

Observations of solar gamma ray continuum between 360 keV and 7 MeV on August 4, 1972

Measurements were made of the time-averaged gamma ray energy loss spectrum in the energy range 360 keV to 7 MeV by the gamma ray detector on the OSO-7 satellite during the 3B flare on August 4, 1972. The differential photon spectrum unfolded from this spectrum after subtracting the background spectrum and contributions from gamma ray lines is best described by a power law with spectral index of 3.4±0.3 between 360–700 keV and by an exponential law of the form exp (-E/E ₀) with E ₀ = 1.0±0.1 MeV above 700 keV. It is suggested that this spectrum is due to nonthermal electron bremsstrahlung from a population of electrons, with a strong break in the spectrum at 2 MeV. Since the observational data indicates that the matter number density must be n _H ? 5 × 10¹⁰ cm^-3 in the production region, the number of electrons above 100 keV required to explain the results is ?2 × 10³⁴.     

Long term variations in solar flare activity

A statistical analysis of the contemporary (1954-1975) solar flare particle events has been made for the parametersF (integrated, proton fluence in cm^-2 in an event with kinetic energy above 10 MeV) andR ₀ (the characteristic rigidity). These data are compared with the long-term averaged values determined from stable- and radio-nuclide measurements of lunar samples. The analysis shows that the ancient solar flare proton spectrum was harder (higher R₀ values) compared to that observed in contemporary flares. A similar analysis can not be made for the mean long-term averaged flux (ˉJ, cm^-2 S^-1), since the contemporary averages suffer from an uncertainty due to the statistics of a single event. However, the average flux estimates for time durations 〈T〉 exceeding 10³ yr, are free from such uncertainties. The long-term averaged ˉJ values obtained over different time scales (10⁴ - 10⁶ yr) suggest a possible periodic variation in solar flare activity, with enhanced flux level during the last 10⁵ yr. The available data rule out the occurrence of giant flares, with proton fluence exceeding 10¹⁵ cm^-2 during the last million years.     

On the correlation between the solar gamma–ray line emission,radio bursts and proton fluxes in the interplanetary space

It is shown that the correlation takes place between the 4–7 MeV gamma–ray line flare fluence F_4–7 and the intensity of the > 10 MeV proton flux in the interplanetary space as well as between F_4–7 and the peak flux density of microwave bursts. Besides, the energy spectral index of protons displays the definite dependence from parameters of the radio burst frequency spectrum. These testify that: a) there is a close physical association between the acceleration of electrons and protons in flares; b) protons, giving gamma–ray lines, and ones, registered in the interplanetary space, belong to the same population.     

EGRET observations of > 30 MeV emission from the brightest bursts detected by BATSE

EGRET on the Compton Gamma Ray Observatory has detected 5 gammaray bursts above 30 MeV. The sub-MeV emission, as detected by BATSE, for these 5 bursts has the largest fluence,F, and peak intensity,I, of any of the BATSE detected bursts within the EGRET field of view. The BATSE data reported in the 2B catalog and the EGRET exposure,E, are combined to select additional bursts with largeF ×E andI ×E. The EGRET data for these bright BATSE bursts are searched for prompt, as well as delayed, emission above 30 MeV. The average properties of the >30 MeV emission are obtained by adding the EGRET data from the 5 bursts. On average the fluence is greater than 15% of the fluence detected by BATSE below an MeV, and the average spectrum is flatter than the spectrum from 1-30 MeV.     

Prompt acceleration of ≳ 30 MeV per nucleon ions in solar flares

Observation of prompt γ-rays in solar flares requires that ions be accelerated to >30 MeV nucl^-1 in ? 2 s. A model for prompt acceleration is developed. The energy release is assumed to occur in a flaring loop with the energy release region being ? 10⁴ km in dimensions and with an Alfvén speed υ _A ? 3 × 10³ km s^-1. The acceleration is assumed to occur in two steps. The second-step acceleration from ? ? _T = 1/2m _p υ_A ² nucl^-1 to ? 30 MeV nucl^-1 is attributed to stochastic acceleration by hydromagnetic turbulence which is found to be fast enough under conditions which are not extreme. Main emphasis is placed on the first step, called preacceleration, to ? _T ? 100 keV nucl^-1. Preacceleration mechanisms which involve accelerating a small fraction of ions from the tail of a Maxwellian distribution are unacceptable because they would lead to enormous abundance anomalies. Preacceleration is attributed either to localized heating of ions to ? 10⁹ K or to acceleration by potential electric fields. The latter mechanism is favoured and some theoretical ideas are outlined based on observations of reconnection in the Earth's magnetotail. Whether energetic ions are prompt, delayed or unobservable depends only on the rate at which the stochastic acceleration proceeds. The second-step acceleration of electrons, invoked to account for a harder microwave component, is predicted to be slower by a factor ? 3 than for ? 30 MeV nucl^-1 ions.     

Peak-Size Distributions of Proton Fluxes and Associated Soft X-Ray Flares

A database combining information about solar proton enhancements (SPEs) near the Earth and soft X-ray flares (GOES measurements) has been used for the study of different correlations through the period from 1975 to May 2006. The emphasis of this work is on the treatment of peak-size distributions of SXR flares and SPEs. The frequency of SXR flares and solar proton events (>10 and >100 MeV, respectively) for the past three solar cycles has been found to follow mainly a power-law distribution over three to five orders of magnitude of fluxes, which is physically correct beyond the “sensitivity” problem with the smallest peak values. The absence of significant spectral steepening in the domain of the highest peak values demonstrates that during the period considered, lasting 30 years, the limit of the highest flare’s energy release has not yet been achieved. The power-law exponents were found to be −2.19±0.04, −1.34±0.02, and −1.46±0.04, for the total SXR flare distribution and the total SPE distributions (for both E _P>10 MeV and E _P>100 MeV), respectively. For SPEs associated with flares located to the West of 20° W, the exponents are −1.22±0.05 (E _P>10 MeV) and −1.26±0.03 (E _P>100 MeV). The size distribution for corresponding flares follows a power law with a slope of −1.29±0.12. Thus, X-ray and proton fluxes produced in the same solar events have very similar distribution shapes. Moreover, the derived slopes are not incompatible with a linear dependence between X-ray flare power and proton fluxes near the Earth. A similar statistical relation is obtained independently from the direct comparison of the X-ray and proton fluxes. These all argue for a statistically significant relationship between X-ray and proton emissions.     

Solar Activity during the Maunder Minimum: Comparison with the Dalton Minimum

The solar modulation potential has been reconstructed from data on the ¹⁰Be concentration in south and central Greenland ice over more than 500 last years. These two reconstructions, along with fourteen others obtained by various authors from data on the cosmogenic isotopes ¹⁴C and ¹⁰Be, have been investigated in the time interval 1630–1840 encompassing the Maunder and Dalton minima. The information contained in these sixteen paleoreconstructions has been generalized. The available data on the concentration of cosmogenic isotopes in terrestrial archives suggest that the solar activity in the first part of theMaunder minimum (1645–1680) was lower than that at the Dalton minimum (1792–1827), while in the second part (1680–1715) it was considerably lower. At the same time, at the beginning of theMaunder minimum (1645–1660) the solar activity could reach levels exceeding noticeably the estimates based on telescopic observations. Possible causes of these discrepancies and the directions of further research are discussed.     

Cosmogenic nuclides in core samples of the Chico L6 chondrite: Evidence for irradiation under high shielding

Abstract— Core samples were obtained from various locations of the ～ 105-kg Chico, NM, L6 chondrite in order to study the effects of large shielding on the production rates of cosmic-ray-produced nuclides. Relations between measured abundances of cosmogenic nuclides (¹⁰Be, ²⁶Al, and stable isotopes of He, Ne, and Ar) and the cosmogenic ²²Ne/²¹Ne ratio were determined and compared with recent model predictions of production rates. The measured ²²Ne/²¹Ne ratios (1.06-1.08) and significant variations observed in concentrations of cosmogenic ²¹Ne and ³He suggest an ～40-cm shielding gradient across Chico and irradiation within a large object (> 100-cm radius). Noble gas data indicate that Chico experienced greater shielding than chondrites Knyahinya or Keyes and similar to Jilin. Values of ¹⁰Be (average = 20.7 dpm/kg) and ²⁶Al (average = 71.1 dpm/kg) are nearly constant, however, and show no correlation with either ²²Ne/²¹Ne or ²¹Ne. Activities of ¹⁰Be and ²⁶Al suggest irradiation in a smaller object (～40–80 cm radius). The ²⁶Al activity and the ²⁶Al/¹⁰Be ratio (average value = 3.42) are both significantly larger than values for most other chondrites. These results could indicate a two-stage irradiation with t₁ ～ 104 Ma and t₂ ～ 4 Ma and a second-stage body the size of Knyahinya. The single stage, ¹⁰Be/²¹Ne exposure age for Chico is 65 Ma. The ²²Ne/²¹Ne ratio apparently becomes insensitive to shielding for objects the size of Chico. No substantial evidence exists for chondrites with ²²Ne/²¹Ne ratios significantly less than ～ 1.055.     

Fall,classification and cosmogenic records of the Sabrum (LL6) chondrite

Abstract— The petrographic and chemical characteristics of a fresh Indian meteorite fall at Sabrum are described. Its mean mineral composition is defined by olivine (Fa_31.4), orthopyroxene (Fs_25.1,Wo_2.0), clinopyroxene (Wo₄₅En_45.6Fs_9.4) and plagioclase (An_10.6Ab_83.6Or_5.8). The meteorite shows moderate shock features, which indicate that it belongs to the S4 category. Based on mineralogical and chemical criteria the meteorite is classified as an LL6 brecciated veined chondrite. Several cosmogenic radioisotopes (⁴⁶Sc, ⁷Be, ⁵⁴Mn, ²²Na and ²⁶Al), noble gas (He, Ne, Ar, Kr and Xe), nitrogen isotopes, and particle tracks density have been measured. Concentrations of cosmogenic ²¹Ne and ³⁸Ar indicate that its cosmic‐ray exposure age is 24.8 Ma. Small amounts of trapped Kr and Xe, consistent with petrologic class 5/6, are present. The track density in olivines is found to be (1.3 ± 0.3) × 10⁶/cm². Activities of most of the short‐lived isotopes are lower than those expected from solar cycle variation. ²²Na/²⁶Al (1.12 ± 0.02) is found to be significantly anomalous, being ?25% lower than expected from the Climax neutron monitor data. These results indicate that the cosmic‐ray flux during the terminal segment of the meteoroid orbit was low. The activities of ²⁶Al and ⁶⁰Co and the track density indicate small meteoroid size with a radius ?15 cm.     

An high resolution FDIRC for the measurement of cosmic-ray isotopic abundances

Measurements of the relative abundance of cosmic isotopes and of the energy dependence of their fluxes may clarify our present understanding on the confinement time of charged cosmic rays in the Galaxy. Experimental studies of these propagation clocks have been carried out by balloon and space missions at energies of a few 100 MeV/amu by means of detection techniques based on multiple dE/dx sampling, coupled with a measurement of the energy released in a thick absorber. At larger energies, the isotopic separation of light nuclei (as, for instance,⁹Be/¹⁰Be) can be achieved by combining a precise measurement of the particle’s rigidity with an high resolution determination of its velocity, via the observation of the Cherenkov effect in a radiator.In this paper, we propose the introduction - for the first time in a space experiment - of the DIRC technique (Detection of Internal Reflected Cherenkov light) for the identification of cosmic-ray isotopes. This type of detector has been successfully used in electron-positron colliders for particle identification and in particular for π-K separation. While for particles with unit charge the light yield is a limiting factor, in the case of a nucleus of charge Z the larger photostatistics (due to the Z² dependence of Cherenkov light emission) is the key to reach an adequate angular resolution to provide a mass discrimination for isotopes of astrophysical interest. We report on the early development phase of a DIRC prototype with a focussing scheme (FDIRC) to collect the Cherenkov light onto a detector plane instrumented with a Silicon PhotoMultiplier (SiPM) array.     

The spectra and light curves of two gamma-ray bursts

During the last half of 1977 the UCSD/MIT Hard X-Ray and Low Energy Gamma-Ray Experiment of HEAO-1 observed two of the three gamma-ray bursts detected by at least three satellites. The first of these bursts (20 October, 1977) had a fluence of (3.1±0.5)×10^–5 erg cm^–2 integrated over the energy range 0.135–2.05 MeV and over its duration of 38.7 s, placing it among the largest bursts observed. The second (10 November, 1977) had a fluence of (2.1±0.8)×10^–5 erg cm^–2 integrated over the energy range 0.125–3 MeV and over its duration of 2.8 s. The light curves of both bursts exhibit time fluctuations down to the limiting time resolution of the detectors (0.1 s). The spectrum of the 20 October, 1977 burst can be fitted with a power law (index –1.93±0.16), which is harder than other reported gamma-ray burst spectral fits. This burst was detected up to 2.05 MeV, and approximately half of its energy was emitted at photon energies above 0.5 MeV. The spectrum of the 10 November, 1977 burst is softer (index –2.4±0.7) and is similar to the spectrum of the 27 April, 1972 burst.Paper presented at the Symposium on Cosmic Gamma-Ray Bursts held at Toulouse, France, 26–29 November, 1979.     

Cosmogenic activation of a natural tellurium target

¹³⁰Te is one of the candidates for the search for neutrinoless double beta decay. It is currently planned to be used in two experiments: CUORE and SNO+. In the CUORE experiment TeO₂ crystals cooled at cryogenic temperatures will be used. In the SNO+ experiment ^natTe will be deployed up to 0.3% loading in the liquid scintillator volume. A possible background for the signal searched for, are the high Q-value, long-lived isotopes, produced by cosmogenic neutron and proton spallation reaction on the target material. A total of 18 isotopes with Q-value larger than 2 MeV and T_1/2 > 20 days have been identified as potential backgrounds. In addition low Q-value, high rate isotopes can be problematic due to pile-up effects, specially in liquid scintillator based detectors. Production rates have been calculated using the ACTIVIA program, the TENDL library, and the cosmogenic neutron and proton flux parametrization at sea level from Armstrong and Gehrels for both long and short lived isotopes. The obtained values for the cross sections are compared with the existing experimental data and calculations. Good agreement has been generally found. The results have been applied to the SNO+ experiment for one year of exposure at sea level. Two possible cases have been considered: a two years of cooling down period deep underground, or a first purification on surface and 6 months of cooling down deep underground. Deep underground activation at the SNOLAB location has been considered.     

Long-Term Modulation of Galactic Cosmic Rays in the Heliosphere by Curvature of the Interplanetary Magnetic Field

Long-term changes of the GCR intensity in the heliosphere were analysed making use of records of cosmogenic isotopes abundances (¹⁰Be and ¹⁴C) in natural archives and neutron monitor data. These data contain information about short-term, 11-year, 22-year, centennial and two-centennial cycles of solar activity. A mechanism to explain the long-term modulation of GCR due to changes in the curvature of the interplanetary magnetic field, resulting from changes in solar activity, is suggested.     

The X-ray source SLX 1732-304 in the globular cluster Terzan 1: The spectral states and an X-ray burst

ART-P/Granat observations of the X-ray burster SLX 1732-304 in the globular cluster Terzan 1 are presented. The X-ray (3–20 keV) fluxes from the source differed by more than a factor of 4 during the observing sessions on September 8 (F _x ? 6.95 × 10^?10 erg cm^?2 s^?1) and October 6, 1990 (F _x ? 1.64 × 10^?10 erg cm^?2 s^?1). The intensity variations of SLX 1732-304 were apparently accompanied by variations in its hardness: whereas the source in its high state had the spectrum with a distinct exponential cutoff typical of bright low-mass X-ray binaries, its low-state spectrum could be satisfactorily described by a simple power law with a photon index α?1.7. During the ART-P observation on September 8, a type I X-ray burst was detected from SLX 1732-304.     

Core-collapse supernovae and the Younger Dryas/terminal Rancholabrean extinctions

Early predictions that some supernovae release large quantities of prompt high energy photons are now corroborated by optical identification of core-collapse supernovae associated with extragalactic GRBS (beamed γ-ray bursts) and XRFS (beamed or un-beamed X-ray flashes). Given the in-galaxy supernova frequency and GRB and XRF recurrence statistics, significant Earth-incident events during the past several million years very likely occurred and nearby events should have affected the Earth and other planetary atmospheres, including terrestrial surface solar UV, the Earth’s climate, and its ecology. The Younger Dryas Stadial (∼12,900 to 11,550 calendar yr BP) began with sharply cooler temperatures in the Earth’s northern hemisphere, regional drought, paleoecological evidence compatible with increased UV, and abrupt increases in cosmogenic ¹⁴C and ¹⁰Be in ice and marine cores and tree rings. In North America, stratigraphic and faunal sequences indicate that a major pulse of mammalian extinctions (at least 23-31 genera) began very close to 12,830 calendar yr BP and was sudden: deposits one century younger are devoid of diverse extinct fauna remains. A 10 s beamed GRB within 2 kpc of the Earth delivers 100 kJ m⁻² fluence to the Earth’s atmosphere, where it causes spallation and catalytic reactions depleting 35-50% O³, and producing excess NO_x species (which favor cooling, drought, and surface fertility), ¹⁴C, and ¹⁰Be. An un-beamed, 10⁵⁰ erg hard photon impulse at ∼250 pc produces similar terrestrial atmospheric effects. A well-characterized massive star supernova, the unusually close Vela event (d = 250 ± 30 pc; total energy of 1-2 × 10⁵¹ erg; age constrained from remnant nebula shock velocities considerations at 13,000-16,000 yr and from the pulsar characteristic age at ∼11,400 yr) may have initiated the Younger Dryas climate change, and caused the extinction of the terminal Rancholabrean fauna.     

Cosmic GRB 060428C detected in the field of view of the IBIS and SPI telescopes onboard the INTEGRAL observatory and its early afterglow

While analyzing the archival data of the INTEGRAL observatory, we detected and localized a cosmic gamma-ray burst recorded on April 28, 2006, by the IBIS/ISGRI and SPI telescopes in their fields of view. Since the burst was not revealed by the INTEGRAL burst alert system (IBAS), information about its coordinates was not distributed in time and no search for its afterglow was conducted. The burst was recorded by the KONUS/WIND and RHES SI satellites. Its 20–200-keV fluence was 2.3 × 10^?6 erg cm^?2, the peak flux was 3.6 × 10^?7 erg cm^?2 s^?1 (3.9 phot. cm^?2 s^?1). The burst had a complex multipeaked profile and stood out among typical bursts by an increase in its hardness with time. At the flux peak, the spectrum was characterized by a photon index α ? ?1.5 and a peak energy E _p ? 95 keV. The burst lasted for ～12 s, after which its afterglow decaying as a power law with an index γ ～ ?4.5 was observed at energies 15–45 keV. The spectral hardness decreased noticeably during the afterglow.     

Stable nickel isotopes and cosmogenic beryllium-10 and aluminum-26 in metallic spheroids from Meteor Crater,Arizona

Abstract— The isotopic abundances of Ni in 17 metallic spheroids from Meteor Crater, Arizona, were determined by inductively coupled plasma mass spectrometry (ICP-MS). Sixteen spheroids have normal isotopic abundances. A 17th shows a marginally detectable mass fractionation of 0.40 ± 0.14 %/AMU in favor of the heavier isotopes. The general absence of mass fractionation indicates that open system evaporation caused little loss of Ni. Variable activities of the cosmogenic radionuclides ¹⁰Be and ²⁶Al were measured by accelerator mass spectrometry in separate suites of spheroids. Activities of ²⁶Al in most samples and of ¹⁰Be in metal cores separated from spheroids indicate that they either (1) come from greater depths in the parent meteoroid than do hand specimens, or (2) lost Al and Be during the process of spheroid formation. One individual spheroid has ¹⁰Be and ²⁶Al activities comparable to those of bulk specimens. This result suggests that spheroid formation may occasionally include material from the outermost meter or so of the impactor. Relatively high activities of ¹⁰Be, ～3 dpm/kg, in the siliceous shells of Canyon Diablo spheroids very likely have a meteoric origin.     

Size Distributions of Solar Proton Events: Methodological and Physical Restrictions

Based on the new catalogue of solar proton events (SPEs) for the period of 1997?–?2009 (Solar Cycle 23) we revisit the long-studied problem of the event-size distributions in the context of those constructed for other solar-flare parameters. Recent results on the problem of size distributions of solar flares and proton events are briefly reviewed. Even a cursory acquaintance with this research field reveals a rather mixed and controversial picture. We concentrate on three main issues: i) SPE size distribution for \({>}\,10~\mbox{MeV}\) protons in Solar Cycle 23; ii) size distribution of \({>}\,1~\mbox{GV}\) proton events in 1942?–?2014; iii) variations of annual numbers for \({>}\,10~\mbox{MeV}\) proton events on long time scales (1955?–?2015). Different results are critically compared; most of the studies in this field are shown to suffer from vastly different input datasets as well as from insufficient knowledge of underlying physical processes in the SPEs under consideration. New studies in this field should be made on more distinct physical and methodological bases. It is important to note the evident similarity in size distributions of solar flares and superflares in Sun-like stars.