Skip to main content


  • Access
  • Open access


      • Send article to Kindle

        To send this article to your Kindle, first ensure is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

        Note you can select to send to either the or variations. ‘’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

        Find out more about the Kindle Personal Document Service.

        Relationship between solar activity and Δ14C peaks in AD 775, AD 994, and 660 BC
        Available formats
        Send article to Dropbox

        To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your Dropbox account. Find out more about sending content to Dropbox.

        Relationship between solar activity and Δ14C peaks in AD 775, AD 994, and 660 BC
        Available formats
        Send article to Google Drive

        To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your Google Drive account. Find out more about sending content to Google Drive.

        Relationship between solar activity and Δ14C peaks in AD 775, AD 994, and 660 BC
        Available formats
Export citation



Since the AD 775 and AD 994 Δ14C peak (henceforth M12) was first measured by Miyake et al. (2012, 2013), several possible production mechanisms for these spike have been suggested, but the work of Mekhaldi et al. (2015) shows that a very soft energy spectrum was involved, implying that a strong solar energetic particle (SEP) event (or series of events) was responsible. Here we present Δ14C values from AD 721–820 Sequoiadendron giganteum annual tree-ring samples from Sequoia National Park in California, USA, together with Δ14C in German oak from 650–670 BC. The AD 721–820 measurements confirm that a sharp Δ14C peak exists at AD 775, with a peak height of approximately 15‰ and show that this spike was preceded by several decades of rapidly decreasing Δ14C. A sharp peak is also present at 660 BC, with a peak height of about 10‰, and published data (Reimer et al. 2013) indicate that it too was preceded by a multi-decadal Δ14C decrease, suggesting that solar activity was very strong just prior to both Δ14C peaks and may be causally related. During periods of strong solar activity there is increased probability for coronal mass ejection (CME) events that can subject the Earth’s atmosphere to high fluencies of solar energetic particles (SEPs). Periods of high solar activity (such as one in October–November 2003) can also often include many large, fast CMEs increasing the probability of geomagnetic storms. In this paper we suggest that the combination of large SEP events and elevated geomagnetic activity can lead to enhanced production of 14C and other cosmogenic isotopes by increasing the area of the atmosphere that is irradiated by high solar energetic particles.


The strongest peak of Δ14C, corresponding to the event of AD 775 (called henceforth M12), and another peak in AD 993/4 were measured by Miyake (Miyake et al. 2012, 2013) using 14C in annually resolved Japanese cedar tree-ring sequences, and M12 has been confirmed by several studies (Usoskin et al. 2013; Jull et al. 2014; Güttler et al. 2015; Rakowski et al. 2015; Büntgen et al. 2016) including three that showed peaks with slightly different starting times. Data from high-latitude Siberian and Altai (Jull et al. 2014; Büntgen et al. 2016) record the AD 775 peak but show a Δ14C rise beginning one year earlier, and Δ14C data in New Zealand kauri (Güttler et al. 2015) indicate a peak delayed by half a year.

Numerous causes of M12 have been suggested and examined (Thomas et al. 2013), including a supernova (Miyake et al. 2012), a gamma ray burst (GRB) (Hambaryan et al. 2013; Pavlov et al. 2013), and a solar proton event, also known as a solar energetic particle (SEP) event (Melott et al. 2012; Miyake et al. 2012; Thomas et al. 2013; Usoskin et al. 2013; Jull et al. 2014; Mekhaldi et al. 2015.) However, all of these hypotheses have some weak points. If M12 was caused by a supernova, a candidate supernova should exist, but observations spanning several hundred years before AD 775 record no such event. Supernovae emit both high-energy particles and gamma rays, but high-energy particles from supernova have sufficient energy spread that extended arrival times spanning at least 10 yr would be expected (Gehrels 2003; Dee et al. 2017) whereas M12 shows a sharp 1–2-yr rise. Gamma rays from supernova and gamma ray bursts have a very short duration time, implying that a short 14C peak could occur, but M12 concentrations of the cosmogenic isotope 10Be measured in Antarctic ice are much too high to have been produced by a GRB (Usoskin et al. 2013). Since no evidence exists for any Δ14C increase due to an historical SEP event (Jull et al. 2014), any solar phenomenon associated with the M12 peak must have been extremely intense and/or was a much more prolific producer of 14C than those in the instrumental record. However, comparison of M12 10Be and 36Cl data from polar ice cores (Mekhaldi et al. 2015) show a larger percentage increase over the long term background for 36Cl than 10Be. This result is consistent with production from a flux of relatively soft spectrum of incoming particles (Webber et al. 2007) and hence supports a solar rather than galactic cosmic ray origin for M12.

This study introduces new M12 related Δ14C data from sequoia tree rings in Sequoia National Park together with German oak Δ14C values from 650–670 BC and published data (Reimer et al. 2013) that span a Δ14C peak at 660 BC. Following an analysis and comparisons of our data with the findings of earlier work, the possible causes and processes associated with the M12 peak are discussed.


Dendrochronological single-year tree-ring samples for AD 721–820 were prepared in the Laboratory of Tree Ring Research (LTRR) at the University of Arizona from sequoia (Sequoiadendron giganteum) in Sequoia National Park, California, USA. Sequoia is an evergreen, common in coastal and montane regions of northern California and is among the tallest living trees on Earth. The tree used for this research had rings large enough to ensure that the tested samples contained only material from the precise indicated year, and each ring was subsampled perpendicular to the ring direction to contain an even representation of the year’s growth. The samples were treated by acid-base-acid and burned to produce CO2, which was reduced to graphite and measured at Keck-Carbon Cycle AMS in University of California, Irvine, USA. Additional aliquots of some samples were prepared and measured independently, especially for the decades preceding the M12 spike.

The German oak samples from Oberhaid in Bavaria, Germany, were also prepared as single-year samples at the Institute of Botany, University of Hohenheim, Stuttgart, Germany, and acid-base-acid treated before combustion and subsequent graphitization at Keck-Carbon Cycle AMS.


M12 Peak Analysis

The M12 Δ14C data for the sequoia are presented and compared with results of Miyake et al. (2012) in Figure 1, while Figure 2 shows six published datasets related to M12 plus results from the present work. Since the Northern Hemisphere (NH) growing season is generally from April to September, the appearance time of M12 in NH datasets is recorded here as AD 775.5 and if a Southern Hemisphere (SH) tree ring showing M12 grew from October AD 774 to April AD 775, the appearance time of M12 is recorded as AD 775.0 (Güttler et al. 2015).

Figure 1 Δ14C (‰) between AD 721 and 820 in tree rings from sequoia in Sequoia National Park, California, USA, compared to Δ14C in Japanese cedar (Miyake et al. 2012). The California data include numerous duplicates prepared from the same rings: the A–E designations indicate the time order of treatment and measurement for different batches.

Figure 2 Δ14C in Siberia (Larch, Jull et al. 2014), California, USA (bristlecone pine, Jull et al. 2014), Japan (cedar, Miyake et al. 2012), Germany (oak, Usoskin et al. 2013), New Zealand (kauri, Güttler et al. 2015), Poland (Quercus, Rakowski et al. 2015) plus average values for California sequoia from this study.

The M12s in five NH datasets appear in AD 775.5, but Siberian larch from 67°31′N (Jull et al. 2014) shows a rise of about 6‰ (< two standard deviations) one year early in AD 774.5, with additional increase of 9‰ in AD 775.5 and 2‰ in AD 776.5. The kauri dataset from New Zealand shows M12 appearing 6 months later than in the NH (Güttler et al. 2015), consistent with differences in the NH and SH growing seasons.

For a more detailed comparison (Table 1) we define baseline and intensity values for M12. The baseline is defined as the 4-yr Δ14C average for AD 770.5–773.5 in the NH (the AD 774.5 point is not used, due to the possible early Δ14C rise in AD 774.5 in the Siberian record) and that for AD 771.0–774.0 in SH. The NH intensity value is defined as the 5-yr average for AD 776.5–780.5, (the AD 775.5 point is not used since some records show Δ14C still rising in AD 775.5) and the average for AD 777.0–781.0 is used for the SH intensity value. The baseline and intensity values in the present work are higher by 2.7±1.4‰ and 3.1±1.2‰ than those of Miyake et al. (2012) respectively, similar to most of the other NH data: whether this represents a real geographic difference or simply an interlaboratory offset is presently unknown. The baseline and intensity values from New Zealand are ~6‰ lower than the European and U.S. data, consistent with measured NH-SH offsets over the last two millennia (Hogg et al. 2013; Reimer et al. 2013). Given the uncertainties for the M12 heights, the magnitude of the M12 spike is essentially identical for six of the seven records, and although the M12 height for the bristlecone pine record is low, it is nevertheless within two standard deviations of the mean of the other six datasets.

Table 1 Baseline and intensity M12 Δ14C values plus heights (=intensity-baseline) for datasets shown in Figure 2 (see text for details).

*Average of data from two laboratories (Usoskin et al. 2013)

The early (774 AD) rise in the Siberian and Altai ∆14C data is difficult to reconcile with 14C results from the 1960s, when high-yield Soviet atmospheric nuclear tests in Novaya Zemlya (~74°N) in September–November 1961 and again in the latter half of 1962 produced very high burdens of excess 14C in the NH stratosphere (Feely et al. 1966). After each major injection, measurements of atmospheric 14CO2 in air samples from NH ground stations (Nydal and Lovseth 1983; Levin et al. 1985) showed a sharp rise beginning in the early spring of the following year and continuing through the entire NH growing season. By 1963 when the largest 14C increases were observed, a large network of sampling stations extending from 78°N to 9°N was in place (Nydal and Lovseth 1983) and recorded an essentially synchronous ∆14C rise for all stations. Since no anomalous early high-latitude tropospheric 14CO2 signal in a bomb peak was observed, it seems unlikely that any of the mechanisms put forward to explain the M12 peak from a high latitude like Siberia’s could generate an equivalent early 14C spike in biomass. Hence, confirmations are required from measurements at other high-latitude sites.

As Figure 2 shows, the M12 peak is delayed in New Zealand kauri, but the delay is consistent with the differences in timing of the NH and SH growing seasons. Furthermore, whereas the maximum SH tropospheric 14C increase during the bomb spike was just 70% of mid-latitude NH values (Manning et al. 1990), reflecting the fact that close to 100% of the input was generated within the NH, the M12 amplitude in the kauri is ~90% of the mean of the NH stations (Table 1), indicating that significant 14C production occurred in both hemispheres.

Large-scale stratospheric air transport is dominated by the so-called Brewer-Dobson circulation (Butchart 2014) characterized by injection of tropospheric air into the stratosphere in the tropics and subsequently by movement towards the winter pole and subsidence at high latitudes, driven primarily by seasonally varying large-scale eddy motions that act to reduce the mean zonal velocity of air parcels and thus induce a poleward drift (Holton et al. 1995). Descending air in the winter polar vortex carrying a high burden of stratospheric tracers is released near the tropopause when the vortex breaks up in the spring, and extratropical injection of stratospheric air into the troposphere peaks in mid-to-late summer (Stohl et al. 2003). Hence in both hemispheres, tropospheric 14CO2 concentrations in the years following the bomb spike showed seasonal cycles, with ∆14C increasing through the spring and summer and reaching a maximum near the end of the growing season (Nydal and Lovseth 1983; Levin et al. 1985; Manning et al. 1990).

The M12 radiocarbon (14C) was produced simultaneously in both hemispheres according to globally measured data and the events occurred in both hemispheres are undoubtedly confirmed by 10Be data from Greenland and Antarctic ice cores (Sigl et al. 2015). Based on this situation, we infer that this occurred in the early months of 775 AD, at the end of the northern winter. The timing was early enough (and sufficient production was generated low in the stratosphere) that the NH stations received high-14C air when cross-tropopause exchange occurred in the NH spring and summer of 775 AD, and hence the signal was recorded in NH wood produced in mid-year. However, the event took place in the end of SH summer so that little or no 14C was injected into troposphere in that hemisphere during the 774/775 AD austral growing season, and the signal only appeared in the New Zealand kauri in the 775/776 AD ring, almost a year after the production event.

The 660 BC Peak and Solar Activity

Figure 3 shows Δ14C data covering the period from 550 to 750 BC including single year European oak results from this study. A sharp peak with a rise time of 3–4 yr and amplitude of 10‰ is present at 660 BC. Although the rise time of the 660 BC peak is longer than that of M12 and the height of the peak is lower, the peaks in several of the M12 records actually take 2 yr to reach their maximum value (Figure 2), and we suggest that the two peaks were probably produced by similar processes and more measurement of the 660 BC peak are needed for confirmation.

Figure 3 Single-year Δ14C data in German oak spanning the 660 BC peak (this work) are shown together with single-year and decadal data from IntCal13 (Reimer et al. 2013).

The ∆14C data in Figure 1 show a decreasing trend (~0.2‰/yr) from 721 to 774 AD and distinct periodic variations of an 11-yr solar cycle [the Schwabe cycle, average variation is a few per mil (‰)] are seen on both sides of the M12 peak (Güttler et al. 2015). This suggests that an interval of several decades prior to M12 showed elevated solar activity, characterized by a strong solar cycle and low 14C production due to enhanced heliomagnetic shielding against galactic cosmic rays. The decadal ∆14C data from 725 to 675 BC in Figure 3 also show a decreasing trend, suggesting that the 660 BC peak was also preceded by a period of strong solar activity. Any downturn in Δ14C immediately preceding the single-year AD 994 peak, which Miyake et al. (2013) suggest is also similar to M12, was very small and short-lived at best, but we note that this peak did occur within an interval of generally low ∆14C. Hence although additional research and data are needed, there is some evidence suggesting that these short-lived Δ14C peaks are correlated with intervals of strong solar activity, and the simplest explanation is that all three have a common cause.

Solar Proton Events and Coronal Mass Ejections

Why could a relationship between strong solar activity and the Δ14C peaks exist? Mekhaldi et al. (2015) have pointed out that the high 36Cl/10Be ratios for the M12 and AD994 peaks in polar ice are consistent with production by relatively low energy particles (tens to hundreds of MeV) and hence point to SEPs rather than galactic cosmic rays (GeV energies) as the source. SEP events occur more frequently during strong solar activity periods (Kurt et al. 2004), and the probability for very intense SEP events increases concurrently. A hypothesized explanation for M12 and similar ∆14C peaks is that the Earth was irradiated by particles from one or more SEP events involving extremely intense particle fluxes, preferably characterized by unusually harder energy spectra than those of most SEP events. The lower energy SEP component is deflected by Earth’s magnetic field and cannot interact with the atmosphere except at the poles, but the high energy tail of the energy spectrum may have contained a sufficiently high fluence of protons to generate the observed burden of cosmogenic isotopes. Recently, superflares from Sun-like stars were observed with a thousand times the energy release of the largest flares on our Sun, although the life stage at which Sun-like stars make superflares and differences between the Sun and Sun-like stars are still under debate (Maehara et al. 2012, 2015).

This is a popular theory (Miyake et al. 2012; Usoskin et al. 2013; Jull et al. 2014; Mekhaldi et al. 2015) but a potential problem is that it requires a SEP fluence at least five times larger than for any event recorded in the era of ground-based or satellite observations (Mekhaldi et al. 2015). Therefore, we consider the fact that large SEP events are often also associated with large geomagnetic storms (Li et al. 2006), which can have the effect of enlarging the area of the atmosphere that is irradiated by solar energetic particles (Leske et al. 2001).

Here we suggest a possible solution to this difficulty, based on the association of SEPs with coronal mass ejection events (CMEs). A CME occurs when magnetic disturbances on the Sun lead to the violent ejection of a large volume of plasma into interplanetary space, in the form of an expanding thick shell that can extend over angular ranges of 120º or more, anchored at the source region on the Sun. Since the velocities of fast CMEs can reach ~2500–3000 km/s, far greater than the ~350–750 km/s speed of the solar wind, a shockwave forms in front of the CME and the shocked solar wind plasma and the CME itself both contain embedded magnetic fields that interact and reconnect with the geomagnetic field. Fast CMEs are almost always accompanied by solar flares, but typically most of the associated SEPs are not from the flare itself, but are solar wind and suprathermal ions (mainly protons) and electrons accelerated to tens to thousands of MeV by the expanding CME shockwave (Gosling 1993; Reames 1999). Since protons at these energies outrun the shock, SEP arrival at the Earth typically begins 1–2 days prior to the interaction of the shockwave and the CME with the geomagnetic field. In the largest events a reservoir forms in the inner heliosphere (Reames et al. 1997) and elevated particle fluxes may persist for a week or more.

The fastest CMEs can reach Earth in a day or less and if the flare site is located near the central meridian the geomagnetic storm can often commence while SEP intensities are near their maximum. Examples include the August 4, 1972, and July 14, 2000, events, among the largest of the space era. In addition, Gopalswamy et al. (2004) and Li et al. (2012) have shown that many CMEs involving especially high fluences of SEPs are twin events where a second CME erupts within 24 hours of the first. The 14C production from the two events would add nearly linearly (Melott et al. 2012). Additionally the passage of the initial shockwave leaves a region in its wake that for several hours contains high densities of suprathermal particles and elevated levels of solar wind turbulence. If a second shock arrives before the suprathermal decay, it had been suggested that particles are apparently accelerated with greater efficiency, producing a hard SEP flux that will generate unusually large numbers of secondary neutrons and correspondingly high levels of 14C. Gopalswamy et al. (2004) find that the second of these two events is typically more intense by a median value of ~7. Furthermore, if the first CME and/or the shocked plasma in front of it contain strong embedded southward-pointing magnetic fields, these southward fields reconnect with the north-directed fields on the front side of the magnetosphere, releasing a large amount of energy. The resulting magnetic storm acts to reduce the geomagnetic field, so that when SEPs accelerated by the second shock reach the neighborhood of the Earth, they may encounter greatly reduced geomagnetic cutoffs and hence can reach an abnormally large fraction of the atmosphere (Leske et al. 2001; Kress et al. 2010).

To our knowledge, this hypothesis of a synergy between unusually intense and high energy SEPs and reduced geomagnetic shielding has not been considered in the derivation of particle fluences required to explain the M12 peak, but elements of it are certainly consistent with CME and SEP observations. Twin CMEs occur relatively frequently (Li et al. 2012; Ding et al. 2013) and multiple large events within a few days (often from the same active region) are common. Gopalswamy et al. (2004) reported that most of the large SEP events originating from CME-driven shock acceleration in 1997 to early 2002 were preceded by an earlier CME-driven shock event originating from the same solar active region within 24 hours earlier. For example, 80 CMEs, numerous intense flares, and five large SEP events occurred over a three-week period in October–November 2003, the four largest from a single, very active region on the Sun (Gopalswamy et al. 2005; Mewaldt et al. 2005). Furthermore, over several periods of some hours within this sequence, the geomagnetic cutoff latitude for 19–27 MeV protons was reduced from 65° to 50–55° (Kress et al. 2010), corresponding approximately to doubling of the fraction of the atmosphere exposed to those particles. The Carrington Event of 1859 involved major magnetic storms on August 28–29 and September 2–3, with auroras visible at geomagnetic latitudes as low as 25° and 18°, respectively (Green and Boardsen 2006), suggesting that two closely spaced CMEs interacted strongly with the geomagnetic field. As noted by Jull et al. (2014), there is no evidence for this event in the annual 14C record, which at first sight appears to challenge the hypothesis outlined above, but we have no information on the intensity, spectrum, and duration of the SEP event and the conditions for efficient particle acceleration leading to a very hard energy spectrum may not have been met (Li et al. 2006).


New Δ14C data from AD 721–820 from sequoia tree rings of California, USA, are presented and compared with the result of Miyake et al. (2012). The calculated M12 height of 15.5±1.4‰ and the appearance time of AD 775 for M12 in our record are same to those of Miyake et al. (2012). Based on the amplitude and timing of the corresponding peaks in other NH records and in New Zealand kauri, 14C spike production was a global event and probably occurred in the late northern winter (January or February) of AD 775. Δ14C data from 550–750 BC from both the present study and Reimer et al. (2013) show that a sharp 10‰ peak with a rise time of 3–4 yr is present at 660 BC. Although the 660 BC rise time is longer than that of M12 and the height is lower, both peaks appear during periods of strong solar activity, and the simplest explanation for the 660 BC and AD 994 spikes and M12 is that all three have a common cause.

The high 36Cl/10Be ratios for the M12 and AD 994 peaks in polar ice (Mekhaldi et al. 2015) support a relatively soft energy spectrum, consistent with production by accelerated solar particles rather than galactic cosmic rays, but the SEP fluxes required are far larger than any of those in the instrumental SEP record. However, if a southward-pointing magnetic field associated with an intense CME event reconnects with the geomagnetic field, an unusually large fraction of the atmosphere may be exposed to particles accelerated efficiently by the shock wave associated with a closely following second CME (Li et al. 2012). The synergistic effects of reduced geomagnetic cutoffs and hard SEP spectra could lead to unusually effective production of 14C and other cosmogenic isotopes, and hence the calculated values for the required particle fluences may be too large. In addition, the largest events often occur during a period of elevated solar activity when elevated SEP intensities and geomagnetic storms can continue for weeks. Evidence for sharp reductions in geomagnetic cutoff values during CMEs and for multi-CME sequences leading to unusually efficient particle acceleration does exist, but determining whether this mechanism can provide a quantitative explanation for the observed 14C spikes will require detailed calculations in future studies.


This study was supported by a research project funded by the Ministry of Science, ICT and Future Planning of Korea. We are grateful for helpful and constructive comments from F Miyake and F Mekhaldi. The work at Caltech was supported by NASA under NNX13A66G, and subcontract 00008864 of NNX15AG09G, and by the NSF under grant AGS-1622487.


To view supplementary material for this article, please visit


U Büntgen , VS Myglan , FC Ljungqvist , M McCormick , N Di Cosmo , M Sigl , J Jungclaus , S Wagner , PJ Krusic , J Esper , JO Kaplan , MAC de Vaan , J Luterbacher , L Wacker , W Tegel , AV Kirdyanovet . 2016. Cooling and societal change during the late antique little ice age from 536 to around 660 AD. Nature Geoscience 9:231236. DOI: 10.1038/NGEO02652.
N Butchart . 2014. The Brewer-Dobson circulation. Reviews of Geophysics 52:157184.
M Dee , B Pope , D Miles , S Manning , F Miyake . 2017. Supernovae and single-year anomalies in the atmospheric radiocarbon record. Radiocarbon 59(2):293302. DOI: 10.1017/RDC.2016.50.
L Ding , Y Jiang , L Zhao , G Li . 2013. The “twin-CME” scenario and large solar energetic particle events in solar cycle 23. The Astrophysical Journal 763:30.
HW Feely , H Seitz , RJ Lagomarsino , PE Biscaye . 1966. Transport and fallout of stratospheric radioactive debris. Tellus 18(2):316328.
N Gehrels , CM Laird , CH Jackman , JK Cannizzo , BJ Mattson , W Chen . 2003. Ozone depletion from nearby supernovae. The Astrophysical Journal 585:11691176.
N Gopalswamy , S Yashiro , S Krucker , G Stenborg , RA Howard . 2004. Intensity variation of large solar energetic particle events associated with coronal mass ejections. Journal of Geophysical Research 109:A12105.
N Gopalswamy , S Yashiro , G Michalek , H Xie , RP Lepping , RA Howard . 2005. Solar source of the largest geomagnetic storm of cycle 23. Geophysical Research Letters 32:L12S09. DOI: 10.1029/2004GL021639.
JT Gosling . 1993. The solar flare myth. Journal of Geophysical Research 98:1893718949.
JL Green , S Boardsen . 2006. Duration and extent of the great auroral storm of 1859. Advances in Space Research 38:130135.
D Güttler , F Adolphi , J Beer , N Bleicher , G Boswijk , M Christl , A Hogg , J Palmer , C Vockenhuber , L Wacker , J Wunder . 2015. Rapid increase in cosmogenic 14C in AD 775 measured in New Zealand kauri trees indicates short-lived increase in 14C production spanning both hemispheres. Earth and Planetary Science Letters 411:290297.
VV Hambaryan , R Neuhäuser . 2013. Galactic short gamma-ray burst as cause for the 14C peak in AD 774/5. Monthly Notices of the Royal Astronomical Society 430:3236.
AG Hogg , Q Hua , PG Blackwell , N Mu , CE Buck , TP Guilderson , TJ Heaton , JG Palmer , PJ Reimer , RW Reimer , CSM Turney , SRH Zimmerman . 2013. SHCal13 Southern Hemisphere calibration, 0–50,000 years cal BP. Radiocarbon 55(4):18891903.
JR Holton , PH Haynes , ME McIntyre , AR Douglass , RB Rood , L Pfister . 1995. Stratosphere-troposphere exhange. Reviews of Gephysics 33(4):403439.
AJT Jull , IP Panyushkina , TE Lange , VV Kukarskih , VS Myglan , KJ Clark , MW Salzer , GS Burr , SW Leavitt . 2014. Excursions in the 14C record at AD 774–775 in tree rings from Russia and America. Geophysical Research Letters 41. DOI: 10.1002/2014GL059874.
BT Kress , CJ Mertens , M Wiltberger . 2010. Solar energetic particle cutoff variations during the 29–31 October 2003 geomagnetic storm. Space Weather 8:S05001. DOI: 10.1029/2009SW000488.
V Kurt , QA Belov , H Mavromichalaki , M Gerontidou . 2004. Statistical analysis of solar proton events. Annales Geophysicae 22:22552271.
RA Leske , RA Mewaldt , EC Stone , TT Rosenvinge . 2001. Observations of geomagnetic cutoff variations during solar energetic particle events and implications for the radiation environment at the space station. Journal of Geophysical Research 106:3001130022.
I Levin , B Kromer , H Schoch-Fischer , M Bruns , M Münnich , D Berdau , JC Vogel , KO Münnich . 1985. 25 years of tropospheric 14C observations in Central Europe. Radiocarbon 27(1):119.
G Li , R Moore , RA Mewaldt , L Zhao , AW Labrador . 2012. A twin-CME scenario for ground level enhancement events. Space Science Reviews 171:141160.
X Li , M Temerin , BT Tsurutani , S Alex . 2006. Modeling 1–2 September 1859 super magnetic storm. Advances in Space Research 38:273279.
H Maehara , T Shibayama , S Notsu , Y Notsu , T Nagao , S Kusaba , S Honda , D Nogami , K Shibata . 2012. Superflares on solar-type stars. Nature 485:478481.
H Maehara , T Shibayama , Y Notsu , S Notsu , S Honda , D Nogami , K Shibata . 2015. Statistical properties of superflares on solar-type stars based on 1-min cadence data. Earth, Planets and Space 67:59.
MR Manning , DC Lowe , WH Melhuish , RJ Sparks , G Wallace , CAM Brenninkmeijer , RC McGill . 1990. The use of radiocarbon measurements in atmospheric studies. Radiocarbon 32(1):3758.
F Mekhaldi , R Muscheler , F Adolphi , A Aldahan , J Beer , JR McConnell , G Possnert , M Sigl , A Svensson , H-A Synal , KC Welten , TE Woodruff . 2015. Multiradionuclide evidence for the solar origin of the cosmic-ray events of AD 774/5 and 993/4. Nature Communications 6:8611. DOI: 10.1038/ncomms9611.
AL Melott , BC Thomas . 2012. Causes of an AD 774–775 14C increase. Nature 491:E1E2.
RA Mewaldt , CMS Cohen , AW Labrador , RA Leske , GM Mason , MI Desai , MD Looper , JE Mazur , RS Selesnick , DK Haggerty . 2005. Proton, helium, and electron spectra during the large solar particle events of October–November 2003. Journal of Geophysical Research 110:A09S18. DOI: 10.1029/2005JA011038.
F Miyake , K Nagaya , K Masuda , T Nakamura . 2012. A signature of cosmic-ray increase in AD 774–775 from tree rings in Japan. Nature 486:240242.
F Miyake , K Masuda , T Nakamura . 2013. Another rapid event in the carbon-14 content of tree rings. Nature Communications 4:1748. DOI: 10.1038/ncomms2783.
R Nydal , K Lӧvseth . 1983. Tracing bomb 14C in the Atmosphere 1962–1980. Journal of Geophysical Research 88:36213642.
AK Pavlov , AV Blinov , AN Konstantinov , VM Ostryakov , GI Vasilyev , MA Vdovina , PA Volkovet . 2013. AD 775 pulse of cosmogenic radionuclides production as imprint of a galactic gamma-ray burst. MNRAS 435:28782884.
AZ Rakowski , M Krapiec , M Huels , J Pawlyta , A Dreves , J Meadows . 2015. Increase of radiocarbon concentration in tree rings from Kujawy (SE Poland) around AD 774–775. Nuclear Instruments and Methods in Physics Research B 361:564568.
DV Reames . 1999. Particle acceleration at the sun and in the heliosphere. Space Science Reviews 90:413491.
DV Reames , SW Kahler , CK Ng . 1997. Spatial and temporal invariance in the spectra of energetic particles in gradual solar events. The Astrophysical Journal 491:414420.
PJ Reimer , E Bard , A Bayliss , JW Beck , PG Blackwell , C Bronk Ramsey , CE Buck , H Cheng , RL Edwards , M Friedrich , PM Grootes , TP Guilderson , H Haflidason , I Hajdas , C Hatte , TJ Heaton , DL Hoffman , AG Hogg , KA Hughen , KF Kaiser , B Kromer , SW Manning , M Niu , RW Reimer , DA Richards , EM Scott , JR Southon , RA Staff , CSM Turney , J van der Plicht . 2013. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon 55(4):18691887.
M Sigl , M Winstrup , JR McConnell , et al. 2015. Timing and climate forcing of volcanic eruptions for the past 2,500 years. Nature 523:543549.
A Stohl , P Bonasoni , P Cristofanelli , et al. 2003. Stratosphere-troposphere exchange: A review, and what we have learned from STACCATO. Journal of Geophysical Research 108(D12). DOI: 10.1029/2002JD002490.
BC Thomas , AL Melott , KR Arkenberg , BR Snyder II. 2013. Terrestrial effects of possible astrophysical sources of an AD 774–775 increase in 14C production. Geophysical Research Letters 40:12371240.
IG Usoskin , B Kromer , F Ludlow , J Beer , M Friedrich , GA Kovaltsov , SK Solanki , L Wacker . 2013. The AD775 cosmic event revisited: the Sun is to blame. Astronomy & Astrophysics 552:L3. DOI: 10.1051/0004-6361/201321080.
WR Webber , PR Higbie , KG McCracken . 2007. Production of the cosmogenic isotopes 3H, 7Be, 10Be, and 36Cl in the Earth’s atmosphere by solar and galactic cosmic rays. Journal of Geophysical Research 112:A10106. DOI: 10.1029/2007JA012499.