Source author record

Adrian L. Melott

Adrian L. Melott appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
20topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2020arXiv

From Cosmic Explosions to Terrestrial Fires?: A Reply

Deschamps and Mottez (hereafter DM) argue that the Gauss-Matuyama terrestrial magnetic field reversal may have left a vanishing main dipole moment to the field for a time of order 10,000 years. They say this may have allowed an enhanced cosmic ray flux, boosting the effect we proposed in Melott and Thomas (2019). We point out that the bulk of the cosmic ray flux from a nearby supernova should be too energetic, up to a million times more energetic than the limits of deflection by the terrestrial magnetic field. In fact, only those highly energetic ones will directly reach the troposphere, relevant for cloud-to-ground lightning. From Cosmic Explosions to Terrestrial Fires?: A Discussion. F. Deschamps and F. Mottez. J. Geology 128, online ahead of print. (2020) From Cosmic Explosions to Terrestrial Fires?: A Reply A.L. Melott and B.C. Thomas. J. Geology 128, online ahead of print. (2020) From cosmic explosions to terrestrial fires? (A.L. Melott and B.C. Thomas) Journal of Geology, 127, 475-481 10.1086/703418 (2019) [arXiv:1903.01501]

preprint2020arXiv

Supernova Triggers for End-Devonian Extinctions

The Late Devonian was a protracted period of low speciation resulting in biodiversity decline, culminating in extinction events near the Devonian-Carboniferous boundary. Recent evidence indicates that the final extinction event may have coincided with a dramatic drop in stratospheric ozone, possibly due to a global temperature rise. Here we study an alternative possible cause for the postulated ozone drop: a nearby supernova explosion that could inflict damage by accelerating cosmic rays that can deliver ionizing radiation for up to $\sim 100$ kyr. We therefore propose that the end-Devonian extinctions were triggered by supernova explosions at $\sim 20$ pc, somewhat beyond the "kill distance" that would have precipitated a full mass extinction. Such nearby supernovae are likely due to core-collapses of massive stars; these are concentrated in the thin Galactic disk where the Sun resides. Detecting either of the long-lived radioisotopes Sm-146 or Pu-244 in one or more end-Devonian extinction strata would confirm a supernova origin, point to the core-collapse explosion of a massive star, and probe supernova nucleosythesis. Other possible tests of the supernova hypothesis are discussed.

preprint2016arXiv

Atmospheric ionization by high-fluence, hard spectrum solar proton events and their probable appearance in the ice core archive

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.

preprint2016arXiv

Reply to comment by K.A. Duderstadt et al. on "Atmospheric ionization by high-fluence, hard spectrum solar proton events and their probable appearance in the ice core archive"

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.

preprint2015arXiv

Comment on 'Investigations into the impact of astronomical phenomena on the terrestrial biosphere and climate' (arXiv:1505.07856 [astro-ph.EP]) by Fabo Feng

This work by Feng and papers which published its conclusions do not cite nor do they deal with objections by the author published in 2013-2014. There are many fundamental problems. We will summarize here the principal problems, as published by Melott and Bambach which render irrelevant most of the work presented by Feng.

preprint2014arXiv

Biological radiation dose from secondary particles in a Milky Way gamma ray burst

Gamma ray bursts (GBRs) are a class of highly energetic explosions emitting radiation in a very short timescale of a few seconds and with a very narrow opening angle. Although, all GRBs observed so far are extragalactic in origin, there is a high probability of a GRB of galactic origin beaming towards the Earth in the past ~ 0.5 Gyr. Such an intense burst of gamma rays would ionize the atmosphere and deplete the ozone layer. With depleted ozone, there will be an increased flux of solar UVB on the Earth's surface with harmful biological effects. In addition to the atmospheric damage, secondary particles produced by gamma ray-induced showers will reach the surface. Amongst all secondary particles, muons dominate the ground-level secondary particle flux (99% of the total number of particles) and are potentially of biological significance. Using the Monte Carlo simulation code CORSIKA, we modeled the air showers produced by gamma ray primaries up to 100 GeV. We found that the number of muons produced by hypothetical galactic GRBs significantly increases the total muon flux. However, since the muon production efficiency is extremely low for photon energy below 100 GeV, and because GRBs radiate strongly for only a very short time, we find that the biological radiation dose from secondary muons is negligible. The main mechanism of biological damage from GRBs is through solar UVB irradiation from the loss of ozone in the upper atmosphere.

preprint2014arXiv

Cosmic Rays and Terrestrial Life: A Brief Review

"The investigation into the possible effects of cosmic rays on living organisms will also offer great interest." - Victor F. Hess, Nobel Lecture, December 12, 1936 High-energy radiation bursts are commonplace in our Universe. From nearby solar flares to distant gamma ray bursts, a variety of physical processes accelerate charged particles to a wide range of energies, which subsequently reach the Earth. Such particles contribute to a number of physical processes occurring in the Earth system. A large fraction of the energy of charged particles gets deposited in the atmosphere, ionizing the atmosphere, causing changes in its chemistry and affecting the global electric circuit. Remaining secondary particles contribute to the background dose of cosmic rays on the surface and parts of the subsurface region. Life has evolved over the past ~ 3 billion years in presence of this background radiation, which itself has varied considerably during the period. As demonstrated by the Miller-Urey experiment, lightning plays a very important role in the formation of complex organic molecules, which are the building blocks of more complex structures forming life. There is growing evidence of increase in the lightning rate with increasing flux of charged particles. Is there a connection between enhanced rate of cosmic rays and the origin of life? Cosmic ray secondaries are also known to damage DNA and cause mutations, leading to cancer and other diseases. It is now possible to compute radiation doses from secondary particles, in particular muons and neutrons. Have the variations in cosmic ray flux affected the evolution of life on earth? We describe the mechanisms of cosmic rays affecting terrestrial life and review the potential implications of the variation of high-energy astrophysical radiation on the history of life on earth.

preprint2013arXiv

Analysis of periodicity of extinction using the 2012 geological time scale

Analysis of two independent data sets with increased taxonomic resolution (genera rather than families) using the revised 2012 time scale reveals that an extinction periodicity first detected by Raup and Sepkoski (1984) for only the post-Paleozoic actually runs through the entire Phanerozoic. Although there is not a local peak of extinction every 27 million years, an excess of the fraction of genus extinction by interval follows a 27 million year timing interval and differs from a random distribution at the p ~ 0.02 level. A 27 million year periodicity in the spectrum of interval lengths no longer appears, removing the question of a possible artifact arising from it. Using a method originally developed in Bambach (2006) we identify 19 intervals of marked extinction intensity, including mass extinctions, spanning the last 470 million years (and with another six present in the Cambrian) and find that 10 of the 19 lie within 3 Myr of the maxima in the spacing of the 27 Myr periodicity, which differs from a random distribution at the p = 0.004 level. These 19 intervals of marked extinction intensity also preferentially occur during decreasing diversity phases of a well-known 62 Myr periodicity in diversity (16 of 19, p = 0.002). Both periodicities appear to enhance the likelihood of increased severity of extinction, but the cause of neither periodicity is known. Variations in the strength of the many suggested causes of extinction coupled to the variation in combined effect of the two different periodicities as they shift in and out of phase is surely one of the reasons that definitive comparative study of the causes of major extinction events is so elusive.

preprint2013arXiv

Cosmogenic nuclide enhancement via deposition from long-period comets as a test of the Younger Dryas impact hypothesis

We explore the idea that detectable excursions in 26Al may arise from direct deposition by any bolide, and excursions in 14C and 10Be abundances in the atmosphere may result from long-period comet impacts. This is very different from the usual processes of production by cosmic rays within Earths atmosphere. Long-period comets experience greatly increased cosmic ray flux beyond the protection of the suns magnetic field. We report the computed amount of 14C, 10Be, and 26Al present on long-period comets as a function of comet mass. We find that the amount of nuclide mass on large long-period comets entering the Earths atmosphere may be sufficient for creating anomalies in the records of 14C and 10Be from past impacts. In particular, the estimated mass of the proposed Younger Dryas comet is consistent with its having deposited sufficient isotopes to account for recorded 14C and 10Be increases at that time. The 26Al/10Be ratio is much larger in extraterrestrial objects than in the atmosphere, and so, we note that measuring this ratio in ice cores is a suitable definitive test for the Younger Dryas impact hypothesis, even if the hypothetical bolide is not a long-period comet and/or did not contribute to the 14C and 10Be increases.

preprint2013arXiv

Declining Volatility, a General Property of Disparate Systems: From Fossils, to Stocks, to the Stars

There may be structural principles pertaining to the general behavior of systems that lead to similarities in a variety of different contexts. Classic examples include the descriptive power of fractals, the importance of surface area to volume constraints, the universality of entropy in systems, and mathematical rules of growth and form. Documenting such overarching principles may represent a rejoinder to the Neodarwinian synthesis that emphasizes adaptation and competition. Instead, these principles could indicate the importance of constraint and structure on form and evolution. Here we document a potential example of a phenomenon suggesting congruent behavior of very different systems. We focus on the notion that universally there has been a tendency for more volatile entities to disappear from systems such that the net volatility in these systems tends to decline. We specifically focus on origination and extinction rates in the marine animal fossil record, the performance of stocks in the stock market, and the characters of stars and stellar systems. We consider the evidence that each is experiencing declining volatility, and also consider the broader significance of this.

preprint2013arXiv

Do Periodicities in Extinction -- with Possible Astronomical Connections -- Survive a Revision of the Geological Timescale?

A major revision of the geological timescale was published in 2012. We re-examine our past finding of a 27 Myr periodicity in marine extinction rates by re-assigning dates to the extinction data used previously. We find that the spectral power at this period is somewhat increased, and persists at a narrow bandwidth, which supports our previous contention that the Nemesis hypothesis is untenable as an explanation for the periodicity that was first noted by Raup and Sepkoski in the 1980s. We enumerate a number of problems in a recent study comparing extinction rates with time series models.

preprint2013arXiv

Terrestrial effects of possible astrophysical sources of an AD 774-775 increase in 14C production

We examine possible sources of a substantial increase in tree ring 14C measurements for the years AD 774-775. Contrary to claims regarding a coronal mass ejection (CME), the required CME energy is not several orders of magnitude greater than known solar events. We consider solar proton events (SPEs) with three different fluences and two different spectra. The data may be explained by an event with fluence about one order of magnitude beyond the October 1989 SPE. Two hard spectrum cases considered here result in moderate ozone depletion, so no mass extinction is implied, though we do predict increases in erythema and damage to plants from enhanced solar UV. We are able to rule out an event with a very soft spectrum that causes severe ozone depletion and subsequent biological impacts. Nitrate enhancements are consistent with their apparent absence in ice core data. The modern technological implications of such an event may be extreme, and considering recent confirmation of superflares on solar-type stars, this issue merits attention.

preprint2012arXiv

Causes of an AD 774-775 14C increase

Atmospheric 14C production is a potential window into the energy of solar proton and other cosmic ray events. It was previously concluded that results from AD 774-775 are orders of magnitude greater than known solar events. We find that the coronal mass ejection energy based on 14C production is much smaller than claimed, but still substantially larger than the maximum historical Carrington Event of 1859. Such an event would cause great damage to modern technology, and in view of recent confirmation of superflares on solar-type stars, this issue merits attention.

preprint2011arXiv

Astrophysical Ionizing Radiation and the Earth: A Brief Review and Census of Intermittent Intense Sources

Cosmic radiation backgrounds are a constraint on life, and their distribution will affect the Galactic Habitable Zone. Life on Earth has developed in the context of these backgrounds, and characterizing event rates will elaborate the important influences. This in turn can be a base for comparison with other potential life-bearing planets. In this review we estimate the intensities and rates of occurrence of many kinds of strong radiation bursts by astrophysical entities ranging from gamma-ray bursts at cosmological distances to the Sun itself. Many of these present potential hazards to the biosphere: on timescales long compared with human history, the probability of an event intense enough to disrupt life on the land surface or in the oceans becomes large. We enumerate the known sources of radiation and characterize their intensities at the Earth and rates or upper limits on these quantities. When possible, we estimate a "lethal interval", our best estimate of how often a major extinction-level event is probable given the current state of knowledge; we base these estimates on computed or expected depletion of stratospheric ozone. In general, moderate level events are dominated by the Sun, but the far more severe infrequent events are probably dominated by gamma-ray bursts and supernovae. We note for the first time that so-called "short-hard" gamma-ray bursts are a substantial threat, comparable in magnitude to supernovae and greater than that of the higher-luminosity long bursts considered in most past work. Given their precursors, short bursts may come with little or no warning.

preprint2011arXiv

Biological implications of high-energy cosmic ray induced muon flux in the extragalactic shock model

A ~ 62 My periodicity in fossil biodiversity has been observed in independent studies of paleontology databases over ~0.5Gy. The period and phase of this biodiversity cycle coincides with the oscillation of our solar system normal to the galactic disk with an amplitude ~70 parsecs and a period ~64 My. Our Galaxy is falling toward the Virgo cluster, forming a galactic shock at the north end of our galaxy due to this motion, capable of accelerating particles and exposing our galaxy's northern side to a higher flux of cosmic rays. These high-energy particles strike the Earth's atmosphere initiating extensive air showers, ionizing the atmosphere by producing charged secondary particles. Secondary particles such as muons produced as a result of nuclear interactions are able to reach the ground and enhance the biological radiation dose. Using a Monte Carlo simulation package CORSIKA, we compute the biological dose resulting from enhanced muon exposure from cosmic rays and discuss their implications for terrestrial biodiversity variations.

preprint2011arXiv

Comment on "N-body simulations of gravitational dynamics" by Dehnen and Read

It is shown that the historical summary of the growth in size of N-body simulations as measured by particle number in this review is missing some key milestones. Size matters, because particle number with appropriate force smoothing is a key method to suppress unwanted discreteness, so that the initial conditions and equations of motion are appropriate to growth by gravitational instability in a Poisson-Vlasov system appropriate to a Universe with dark matter. Published strong constraints on what can be done are not included in the review.

preprint2011arXiv

Modeling high-energy cosmic ray induced terrestrial muon flux: A lookup table

On geological timescales, the Earth is likely to be exposed to an increased flux of high energy cosmic rays (HECRs) from astrophysical sources such as nearby supernovae, gamma ray bursts or by galactic shocks. Typical cosmic ray energies may be much higher than the ~ 1 GeV flux which normally dominates. These high-energy particles strike the Earth's atmosphere initiating an extensive air shower. As the air shower propagates deeper, it ionizes the atmosphere by producing charged secondary particles. Secondary particles such as muons and thermal neutrons produced as a result of nuclear interactions are able to reach the ground, enhancing the radiation dose. Muons contribute 85% to the radiation dose from cosmic rays. This enhanced dose could be potentially harmful to the biosphere. This mechanism has been discussed extensively in literature but has never been quantified. Here, we have developed a lookup table that can be used to quantify this effect by modeling terrestrial muon flux from any arbitrary cosmic ray spectra with 10 GeV - 1 PeV primaries. This will enable us to compute the radiation dose on terrestrial planetary surfaces from a number of astrophysical sources.

preprint2010arXiv

A ubiquitous ~62 Myr periodic fluctuation superimposed on general trends in fossil biodiversity: II, Evolutionary dynamics associated with periodic fluctuation in marine diversity

We investigate evolutionary dynamics related to periodicity fossil biodiversity. Coherent periodic fluctuation in origination/extinction of marine genera that survive <45 million years is the source of an observed ~62 million year periodicity analyzed in Paper I. We also show that the evolutionary dynamics of "long-lived" genera (those that survive >45 million years) do not participate in the periodic fluctuation in diversity and differ from those of "short-lived" genera. The difference between the evolutionary dynamics of these 2 genera classes indicates that the periodic pattern is not an artifact of variation in quality of the geologic record. The interplay of these two previously undifferentiated systems, together with the secular increase in abundance of "long-lived" genera, is probably the source of heretofore unexplained differences in evolutionary dynamics between the Paleozoic and post-Paleozoic as reported by others. Testing for cycles similar to the 62 Myr cycle in fossil biodiversity superimposed on the long-term trends of the Phanerozoic as described in Paper I, we find a significant (but weaker) signal in sedimentary rock packages, particularly carbonates, which suggests a connection. The presence of a periodic pattern in evolutionary dynamics of the vulnerable "short-lived" component of marine fauna demonstrates that a long-term periodic fluctuation in environmental conditions capable of affecting evolution in the marine realm characterizes our planet. Coincidence in timing is more consistent with a common cause than sampling bias. A previously identified set of mass extinctions preferentially occur during the declining phase of the 62 Myr periodicity, supporting the idea that the periodicity relates to variation in biotically important stresses. Further work should focus on finding links to physical phenomena that might reveal the causal system or systems.

preprint2010arXiv

A ubiquitous ~62-Myr periodic fluctuation superimposed on general trends in fossil biodiversity. I. Documentation

We use Fourier analysis and related techniques to investigate the question of periodicities in fossil biodiversity. These techniques are able to identify cycles superimposed on the long-term trends of the Phanerozoic. We review prior results and analyze data previously reduced and published. Joint time-series analysis of various reductions of the Sepkoski Data, Paleobiology Database, and Fossil Record 2 indicate the same periodicity in biodiversity of marine animals at 62 Myr. We have not found this periodicity in the terrestrial fossil record. We have found that the signal strength decreases with time because of the accumulation of apparently "resistant" long-lived genera. The existence of a 62-Myr periodicity despite very different treatment of systematic error, particularly sampling-strength biases, in all three major databases strongly argues for its reality in the fossil record.

preprint2010arXiv

Can periodicity in low altitude cloud cover be induced by cosmic ray variability in the extragalactic shock model?

Variation in high energy cosmic rays (HECRs) has been proposed to explain a 62 My periodicity in terrestrial fossil biodiversity. It has been suggested that the infall of our galaxy toward the Virgo cluster could generate an extragalactic shock, accelerating charged particles and exposing the earth to a flux of high energy cosmic rays (HECRs). The oscillation of the Sun perpendicular to the galactic plane could induce 62 My periodicity in the HECR flux on the Earth, with a magnitude much higher than the Galactic cosmic ray change we see in a solar cycle. This mechanism could potentially explain the observed 62 My periodicity in terrestrial biodiversity over the past 500 My. In addition to direct effects on life from secondaries, HECRs induced air showers ionize the atmosphere leading to changes in atmospheric chemistry and microphysical processes that can lead to cloud formation including low altitude cloud cover. An increase in ionization changes the global electric circuit which could enhance the formation of cloud condensation nuclei (CCN) through microphysical processes such as electroscavenging and ion mediated nucleation, leading to an increase in the cloud cover. This could increase the albedo and reduce the solar flux reaching the ground, reducing the global temperature. Using an existing model, we have calculated the enhancement in atmospheric ionization at low altitudes resulting from exposure to HECRs. We use a conservative model to estimate the change in low altitude cloud cover from this increased ionization.

preprint2010arXiv

Cometary airbursts and atmospheric chemistry: Tunguska and a candidate Younger Dryas event

We find agreement between models of atmospheric chemistry changes from ionization for the 1908 Tunguska airburst event and nitrate enhancement in GISP2H and GISP2 ice cores, plus an unexplained ammonium spike. We then consider a candidate cometary impact at the Younger Dryas onset (YD). The large estimated NO_x production and O_3 depletion are beyond accurate extrapolation, but the ice core peak is much lower, possibly because of insufficient sampling resolution. Ammonium and nitrate spikes have been attributed to biomass burning at YD onset in both GRIP and GISP2 ice cores. A similar result is well-resolved in Tunguska ice core data, but that forest fire was far too small to account for this. Direct input of ammonia from a comet into the atmosphere is adequate for YD ice core data, but not Tunguska data. An analog of the Haber process with hydrogen contributed by cometary or surface water, atmospheric nitrogen, high pressures, and possibly catalytic iron from a comet could in principle produce ammonia, accounting for the peaks in both data sets.

preprint2010arXiv

Lookup tables to compute high energy cosmic ray induced atmospheric ionization and changes in atmospheric chemistry

A variety of events such as gamma-ray bursts and supernovae may expose the Earth to an increased flux of high-energy cosmic rays, with potentially important effects on the biosphere. Existing atmospheric chemistry software does not have the capability of incorporating the effects of substantial cosmic ray flux above 10 GeV . An atmospheric code, the NASA-Goddard Space Flight Center two-dimensional (latitude, altitude) time-dependent atmospheric model (NGSFC), is used to study atmospheric chemistry changes. Using CORSIKA, we have created tables that can be used to compute high energy cosmic ray (10 GeV - 1 PeV) induced atmospheric ionization and also, with the use of the NGSFC code, can be used to simulate the resulting atmospheric chemistry changes. We discuss the tables, their uses, weaknesses, and strengths.

preprint2010arXiv

Nemesis Reconsidered

The hypothesis of a companion object (Nemesis) orbiting the Sun was motivated by the claim of a terrestrial extinction periodicity, thought to be mediated by comet showers. The orbit of a distant companion to the Sun is expected to be perturbed by the Galactic tidal field and encounters with passing stars, which will induce variation in the period. We examine the evidence for the previously proposed periodicity, using two modern, greatly improved paleontological datasets of fossil biodiversity. We find that there is a narrow peak at 27 My in the cross-spectrum of extinction intensity time series between these independent datasets. This periodicity extends over a time period nearly twice that for which it was originally noted. An excess of extinction events are associated with this periodicity at 99% confidence. In this sense we confirm the originally noted feature in the time series for extinction. However, we find that it displays extremely regular timing for about 0.5 Gy. The regularity of the timing compared with earlier calculations of orbital perturbation would seem to exclude the Nemesis hypothesis as a causal factor.

preprint2009arXiv

Testing the link between terrestrial climate change and Galactic spiral arm transit

We re-examine past suggestions of a close link between terrestrial climate change and the Sun's transit of spiral arms in its path through the Milky Way galaxy. These links produced concrete fits, deriving the unknown spiral pattern speed from terrestrial climate correlations. We test these fits against new data on spiral structure based on CO data that does not make simplifying assumptions about symmetry and circular rotation. If we compare the times of these transits to changes in the climate of Earth, not only do the claimed correlations disappear, but also we find that they cannot be resurrected for any reasonable pattern speed.

preprint2008arXiv

Long-term cycles in the history of life: Periodic biodiversity in the Paleobiology Database

Time series analysis of fossil biodiversity of marine invertebrates in the Paleobiology Database (PBDB) shows a significant periodicity at approximately 63 My, in agreement with previous analyses based on the Sepkoski database. I discuss how this result did not appear in a previous analysis of the PBDB. The existence of the 63 My periodicity, despite very different treatment of systematic error in both PBDB and Sepkoski databases strongly argues for consideration of its reality in the fossil record. Cross-spectral analysis of the two datasets finds that a 62 My periodicity coincides in phase by 1.6 My, equivalent to better than the errors in either measurement. Consequently, the two data sets not only contain the same strong periodicity, but its peaks and valleys closely correspond in time. Two other spectral peaks appear in the PBDB analysis, but appear to be artifacts associated with detrending and with the increased interval length. Sampling-standardization procedures implemented by the PBDB collaboration suggest that the signal is not an artifact of sampling bias. Further work should focus on finding the cause of the 62 My periodicity.

preprint2006arXiv

Terrestrial Consequences of Spectral and Temporal Variability in Ionizing Photon Events

Gamma-Ray Bursts (GRBs) directed at Earth from within a few kpc may have damaged the biosphere, primarily though changes in atmospheric chemistry which admit greatly increased Solar UV. However, GRBs are highly variable in spectrum and duration. Recent observations indicate that short (~0.1 s) burst GRBs, which have harder spectra, may be sufficiently abundant at low redshift that they may offer an additional significant effect. A much longer timescale is associated with shock breakout luminosity observed in the soft X-ray (~10^3 s) and UV (~10^5 s) emission, and radioactive decay gamma-ray line radiation emitted during the light curve phase of supernovae (~10^7 s). Here we generalize our atmospheric computations to include a broad range of peak photon energies and investigate the effect of burst duration while holding total fluence and other parameters constant. The results can be used to estimate the probable impact of various kinds of ionizing events (such as short GRBs, X-ray flashes, supernovae) upon the terrestrial atmosphere. We find that the ultimate intensity of atmospheric effects varies only slightly with burst duration from 10^-1 s to 10^8 s. Therefore, the effect of many astrophysical events causing atmospheric ionization can be approximated without including time development. Detailed modeling requires specification of the season and latitude of the event. Harder photon spectra produce greater atmospheric effects for spectra with peaks up to about 20 MeV, because of greater penetration into the stratosphere.

preprint2005arXiv

Gamma-Ray Bursts and the Earth: Exploration of Atmospheric, Biological, Climatic and Biogeochemical Effects

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.

preprint2004arXiv

Numerical study of the cosmological velocity field as a function of density

We report on a new study of the velocity distribution in N-body simulations. We investigate the center-of-mass and internal kinetic energies of coarsening cells as a function of time, cell size and cell mass. By using self-similar cosmological models, we are able to derive theoretical predictions for comparison and to assess the influence of finite-size and resolution effects. The most interesting result is the discovery of a polytropic-like relationship between the average velocity dispersion (internal kinetic energy) and the mass density in an intermediate range of densities, K ~ rho^{2-eta}. The exponent eta measures the deviations from the virial prediction, eta_virial=0. For self-similar models, eta depends only on the spectral index of the initial power spectrum. We also study CDM models and confirm a previous result that the same polytropic-like dependence exists (astro-ph/0103313), with a time and coarsening length dependent eta. The dependence K(rho) is an important input for a recently proposed theoretical model of cosmological structure formation which improves over the standard dust model (pressureless fluid) by regularizing the density singularities.

preprint1999arXiv

Disentangling the Cosmic Web I: Morphology of Isodensity Contours

We apply Minkowski functionals and various derived measures to decipher the morphological properties of large-scale structure seen in simulations of gravitational evolution. Minkowski functionals of isodensity contours serve as tools to test global properties of the density field. Furthermore, we identify coherent objects at various threshold levels and calculate their partial Minkowski functionals. We propose a set of two derived dimensionless quantities, planarity and filamentarity, which reduce the morphological information in a simple and intuitive way. Several simulations of the gravitational evolution of initial power-law spectra provide a framework for systematic tests of our method.

preprint1996arXiv

Testing Tree-Level Perturbation Theory for Large-Scale Structure with the Local Lagrangian Approximation

We test tree-level perturbation theory for Gaussian initial conditions with power spectra $P(k)\propto k^n$ by comparing the probability distribution function (PDF) for the density predicted by the Local Lagrangian Approximation (LLA) with the results of numerical gravitational clustering simulations. Our results indicate that our approximation correctly reproduces the evolved density PDF for $-3 \leq n \leq-1$ power spectra up to the weakly nonlinear regime, while it shows marginal agreement for power indices n=0 and +1 in the linear regime and poor agreement beyond this point. This suggests that tree-level perturbation theory (as realized in the Local Lagrangian Approximation) can accurately predict the density distribution function for $-3 \leq n \leq -1$ but fails for $n \ge 0$.

preprint1996arXiv

The Ellipticity and Orientation of Clusters of Galaxies from N-Body Experiments

In this study we use simulations of 128$^3$ particles to study the ellipticity and orientation of clusters of galaxies in N-body simulations of differing power-law initial spectra (P(k) \propto k^n ,n = +1, 0, -1, -2$), and density parameters ($Ω_0 = 0.2$ to 1.0). Furthermore, unlike most theoretical studies we mimic most observers by removing all particles which lie at distances greater than 2 1/h Mpc from the cluster center of mass. We computed the axial ratio and the principal axes using the inertia tensor of each cluster. The mean ellipticity of clusters increases strongly with increasing $n$. We also find that clusters tend to become more spherical at smaller radii. We compared the orientation of a cluster to the orientation of neighboring clusters as a function of distance (correlation). In addition, we considered whether a cluster's major axis tends to lie along the line connecting it to a neighboring cluster, as a function of distance (alignment). Both alignments and correlations were computed in three dimensions and in projection to mimic observational surveys. Our results show that significant alignments exist for all spectra at small separations ($D < 15 h^{-1}$ Mpc) but drops off at larger distance in a strongly $n-$dependent way.

preprint1994arXiv

Biased Power Spectrum and Bispectrum for an Ensemble of Three-Dimensional Scale Free Numerical Simulations

We examine the effect of a threshold bias on the power spectrum and the bispectrum in an ensemble of numerical simulations (Gaussian initial perturbations with power law spectra P(k) \sim k^n, n=+1, 0, -1, -2) and compare our results with theoretical predictions. Our simulations are evolved sufficiently that on the scale where we apply the threshold the rms fluctuation has developed significantly into the nonlinear regime. Thus, predictions based on perturbation theory do not necessarily apply. Nevertheless, we find our results for the power spectrum, biased power simply amplified by a numerical factor, follow predicted trends, far beyond the regime where perturbation theory is expected to be valid. We find that the biased bispectrum continues to follow the so-called hierarchical form, with reduced three-point amplitude Q \approx 1 in the strongly nonlinear regime, independent of initial spectrum. In the quasi-linear perturbative regime the three-point amplitude depends on configuration shape, a behavior that is found to give useful information about the amount of bias without information about the unbiased matter distribution.

preprint1994arXiv

Hierarchical Pancaking: Why the Zel'dovich Approximation Describes Coherent Large-Scale Structure in N-Body Simulations of Gravitational Clustering

To explain the rich structure of voids, clusters, sheets, and filaments apparent in the Universe, we present evidence for the convergence of the two classic approaches to gravitational clustering, the ``pancake'' and ``hierarchical'' pictures. We compare these two models by looking at agreement between individual structures -- the ``pancakes'' which are characteristic of the Zel'dovich Approximation (ZA) and also appear in hierarchical N-body simulations. We find that we can predict the orientation and position of N-body simulation objects rather well, with decreasing accuracy for increasing large-$k$ (small scale) power in the initial conditions. We examined an N-body simulation with initial power spectrum $P(k) \propto k^3$, and found that a modified version of ZA based on the smoothed initial potential worked well in this extreme hierarchical case, implying that even here very low-amplitude long waves dominate over local clumps (although we can see the beginning of the breakdown expected for $k^4$). In this case the correlation length of the initial potential is extremely small initially, but grows considerably as the simulation evolves. We show that the nonlinear gravitational potential strongly resembles the smoothed initial potential. This explains why ZA with smoothed initial conditions reproduces large-scale structure so well, and probably why our Universe has a coherent large-scale structure.

preprint1993arXiv

Decaying Neutrinos in Galaxy Clusters

Davidsen et al. (1991) have argued that the failure to detect uv photons from the dark matter DM) in cluster A665 excludes the decaying neutrino hypothesis. Sciama et al. (1993) argued that because of high central concentration the DM in that cluster must be baryonic. We study the DM profile in clusters of galaxies simulated using the Harrison--Zel'dovich spectrum of density fluctuations, and an amplitude previously derived from numerical simulations (Melott 1984b; Anninos et al. 1991) and in agreement with microwave background fluctuations (Smoot et al. 1992). We find that with this amplitude normalization cluster neutrino DM densities are comparable to observed cluster DM values. We conclude that given this normalization, the cluster DM should be at least largely composed of neutrinos. The constraint of Davidsen et al. can be somewhat weakened by the presence of baryonic DM; but it cannot be eliminated given our assumptions.

preprint1993arXiv

Testing the Frozen-Flow Approximation

We investigate the accuracy of the frozen--flow approximation (FFA), recently proposed by Matarrese \etal (1992), for following the nonlinear evolution of cosmological density fluctuations under gravitational instability. We compare a number of statistics between results of the FFA and nbody simulations, including those used by Melott, Pellman \& Shandarin (1993) to test the Zel'dovich approximation. The FFA performs reasonably well in a statistical sense, e.g. in reproducing the counts--in--cell distribution, at small scales, but it does poorly in the crosscorrelation with nbody which means it is generally not moving mass to the right place, especially in models with high small--scale power.