Solar flare proton spectrum averaged over the last 5000 years

The solar flare differential energy spectrum and integral proton flux above 10 MeV averaged over the last several thousand years are derived from thermoluminescence (TL) data. The dose-rate depth profile is obtained using a new TL equilibrium technique which is independent of the TL decay kinetics. The dose-rate depth profile produced by solar flare protons with a differential energy spectrum of the form dJ/dE = KE to the minus gamma power is calculated for arbitrary gamma. The best fit to the TL data in rock 14310 is obtained for gamma = 2.3 plus or minus 0.2 and an omnidirectional (4 pi) integral flux above 10 MeV of 40 to 80 prot/sq cm/sec. The TL half-life is determined to be 2000 yr. These results are compared to those for Na-22 (half-life of. 2.6 yr) and Al-26 (half-life of 740,000 yr) obtained by Wahlen et al. (1972) and Rancitelli et al. (1972), and it is concluded that the spectral shape and flux of protons in the interval from 25 to 100 MeV is the same within experimental errors when averaged over these three very different time periods.