Duderstadt et al. [2016b] comment that the Melott et al. [2016] study of nitrate formation by solar proton events (SPEs) and comparison with the ice core archive is "fundamentally flawed," because it does not include pre-existing HNO3 in the stratosphere. We show that they exaggerate both the enhancement predicted by our findings and pre-industrial HNO3 levels in their model, and fail to prove this assertion. Our feasibility study matched expected SPE nitrate production with ground truth measurements. It is not clear that their approach is more realistic and absence of a detailed mechanism does not disprove our results. Models can be no better than the information they are provided and in this case there continue to be significant unknowns and uncertainties, especially in the role of polar stratospheric clouds (PSCs) and possible interactions with cosmic rays that constitute lower boundary conditions. Duderstadt et al. [2014; 2016a] used incomplete, poorly-constrained and incorrect initial and boundary conditions, and they continue to advocate on the basis of uncertain results. Meanwhile, Smart et al. [2014] identified a series of ice core nitrate spikes that have since been confirmed in 10Be by McCracken and Beer [2015]. Melott et al. [2016] computationally reproduced the ionization profile of the only major balloon measurements to date. We show that our calculated nitrate enhancement is consistent with measured results, given current levels of uncertainty, and that extreme SPEs can potentially produce occasional nitrate spikes with hundreds of percent increases. Instead of repeating old arguments to dismiss nitrates as proxies of SPEs, it is past time for a dedicated, fine-resolution, multi-parameter, replicate ice core field campaign to resolve this debate.

Solar energetic particles ionize the atmosphere, leading to production of nitrogen oxides. It has been suggested that some such events are visible as layers of nitrate in ice cores, yielding archives of energetic, high fluence solar proton events (SPEs). There has been controversy, due to slowness of transport for these species down from the upper stratosphere; past numerical simulations based on an analytic calculation have shown very little ionization below the mid stratosphere. These simulations suffer from deficiencies: they consider only soft SPEs and narrow energy ranges; spectral fits are poorly chosen; with few exceptions secondary particles in air showers are ignored. Using improved simulations that follow development of the proton-induced air shower, we find consistency with recent experiments showing substantial excess ionization down to 5 km. We compute nitrate available from the 23 February 1956 SPE, which had a high fluence, hard spectrum, and well-resolved associated nitrate peak in a Greenland ice core. For the first time, we find this event can account for ice core data with timely (~ 2 months) transport downward between 46 km and the surface, thus indicating an archive of high fluence, hard spectrum SPE covering the last several millennia. We discuss interpretations of this result, as well as the lack of a clearly-defined nitrate spike associated with the soft-spectrum 3-4 August 1972 SPE. We suggest that hard-spectrum SPEs, especially in the 6 months of polar winter, are detectable in ice cores, and that more work needs to be done to investigate this.

Ice cores are archives of climate change and possibly large solar proton events (SPEs). Wolff et al. (2012) used a single event, a nitrate peak in the GISP2-H core, which McCracken et al. (2001a) time associated with the poorly quantified 1859 Carrington event, to discredit SPE-produced, impulsive nitrate deposition in polar ice. This is not the ideal test case. We critique the Wolff et al. analysis and demonstrate that the data they used cannot detect impulsive nitrate events because of resolution limitations. We suggest re-examination of the top of the Greenland ice sheet at key intervals over the last two millennia with attention to fine resolution and replicate sampling of multiple species. This will allow further insight into polar depositional processes on a sub-seasonal scale, including atmospheric sources, transport mechanisms to the ice sheet, post-depositional interactions, and a potential SPE association.

Firestone (2014) asserted evidence for numerous (23) nearby (d<300 pc) supernovae within the Middle and Late Pleistocene. If true, this would have strong implications for the irradiation of the Earth; at this rate, mass extinction level events due to supernovae would be more frequent than 100 Myr. However, there are numerous errors in the application of past research. The paper overestimates likely nitrate and 14C production from moderately nearby supernovae by about four orders of magnitude. Moreover, the results are based on wrongly selected (obsolete) nitrate and 14C datasets. The use of correct and up-to-date datasets does not confirm the claimed results. The claims in the paper are invalidated.

Gamma-Ray Bursts (GRBs) are likely to have made a number of significant impacts on the Earth during the last billion years. We have used a two-dimensional atmospheric model to investigate the effects on the Earth's atmosphere of GRBs delivering a range of fluences, at various latitudes, at the equinoxes and solstices, and at different times of day. We have estimated DNA damage levels caused by increased solar UVB radiation, reduction in solar visible light due to $\mathrm{NO_2}$ opacity, and deposition of nitrates through rainout of $\mathrm{HNO_3}$. For the ``typical'' nearest burst in the last billion years, we find globally averaged ozone depletion up to 38%. Localized depletion reaches as much as 74%. Significant global depletion (at least 10%) persists up to about 7 years after the burst. Our results depend strongly on time of year and latitude over which the burst occurs. We find DNA damage of up to 16 times the normal annual global average, well above lethal levels for simple life forms such as phytoplankton. The greatest damage occurs at low to mid latitudes. We find reductions in visible sunlight of a few percent, primarily in the polar regions. Nitrate deposition similar to or slightly greater than that currently caused by lightning is also observed, lasting several years. We discuss how these results support the hypothesis that the Late Ordovician mass extinction may have been initiated by a GRB.

Based on cosmological rates, it is probable that at least once in the last Gy the Earth has been irradiated by a gamma-ray burst in our Galaxy from within 2 kpc. Using a two-dimensional atmospheric model we have performed the first computation of the effects upon the Earth's atmosphere of one such impulsive event. A ten second burst delivering 100 kJ/m^2 to the Earth penetrates to the stratosphere and results in globally averaged ozone depletion of 35%, with depletion reaching 55% at some latitudes. Significant global depletion persists for over 5 years after the burst. This depletion would have dramatic implications for life since a 50% decrease in ozone column density results in approximately three times the normal UVB flux. Widespread extinctions are likely, based on extrapolation from UVB sensitivity of modern organisms. Additional effects include a shot of nitrate fertilizer and NO2 opacity in the visible providing a cooling perturbation to the climate over a similar timescale. These results lend support to the hypothesis that a GRB may have initiated the late Ordovician mass extinction (Melott et al. 2004).

Gamma-ray bursts (hereafter GRB) produce a flux of radiation detectable across the observable Universe, and at least some of them are associated with galaxies. A GRB within our own Ggalaxy could do considerable damage to the Earth's biosphere; rate estimates suggest that a dangerously near GRB should occur on average two or more times per billion years. At least five times in the history of life, the Earth experienced mass extinctions that eliminated a large percentage of the biota. Many possible causes have been documented, and GRB may also have contributed. The late Ordovician mass extinction approximately 440 million years ago may be at least partly the result of a GRB. A special feature of GRB in terms of terrestrial effects is a nearly impulsive energy input of order 10 s. Due to expected severe depletion of the ozone layer, intense solar ultraviolet radiation would result from a nearby GRB, and some of the patterns of extinction and survivorship at this time may be attributable to elevated levels of UV radiation reaching the Earth. In addition a GRB could trigger the global cooling which occurs at the end of the Ordovician period that follows an interval of relatively warm climate. Intense rapid cooling and glaciation at that time, previously identified as the probable cause of this mass extinction, may have resulted from a GRB.

Estimates made in the 1970's indicated that a supernova occurring within tens of parsecs of Earth could have significant effects on the ozone layer. Since that time, improved tools for detailed modeling of atmospheric chemistry have been developed to calculate ozone depletion, and advances have been made in theoretical modeling of supernovae and of the resultant gamma-ray spectra. In addition, one now has better knowledge of the occurrence rate of supernovae in the galaxy, and of the spatial distribution of progenitors to core-collapse supernovae. We report here the results of two-dimensional atmospheric model calculations that take as input the spectral energy distribution of a supernova, adopting various distances from Earth and various latitude impact angles. In separate simulations we calculate the ozone depletion due to both gamma-rays and cosmic rays. We find that for the combined ozone depletion roughly to double the ``biologically active'' UV flux received at the surface of the Earth, the supernova must occur at <8 pc. Based on the latest data, the time-averaged galactic rate of core-collapse supernovae occurring within 8 pc is ~1.5/Gyr. In comparing our calculated ozone depletions with those of previous studies, we find them to be significantly less severe than found by Ruderman (1974), and consistent with Whitten et al. (1976). In summary, given the amplitude of the effect, the rate of nearby supernovae, and the ~Gyr time scale for multicellular organisms on Earth, this particular pathway for mass extinctions may be less important than previously thought.