Evolution of the cosmological density distribution function

We present a new calculation for the evolution of the one-point probability distribution function (PDF) of the cosmological density field based on an exact statistical treatment. Using the Chapman–Kolmogorov equation and second-order Eulerian perturbation theory we propagate the initial density distribution into the non-linear regime. Our calculations yield the moment generating function, allowing a straightforward derivation of the skewness of the PDF to second order. We find a new dependence on the initial perturbation spectrum. We compare our results with other approximations to the one-point PDF, and with N -body simulations. We find that our distribution accurately models the evolution of the one-point PDF of dark matter.     

Crystal scale anatomy of a dying supervolcano: an isotope and geochronology study of individual phenocrysts from voluminous rhyolites of the Yellowstone caldera

A voluminous (>600 km³) and long-lived (~520–75 ka) phase of rhyolitic eruptions followed collapse of the Yellowstone caldera 640 ka. Whether these eruptions represent a dying cycle, or the growth of a new magma chamber, remains an important question. We use new U–Th zircon ages and δ¹⁸O values determined by ion microprobe, and sanidine Pb isotope ratios determined by laser ablation, to investigate the genesis of voluminous post-caldera rhyolites. The oldest post-caldera rhyolites, erupted between ~520 and 470 ka, exhibit extreme age and oxygen isotopic heterogeneity, requiring derivation from individual parcels of low-δ¹⁸O melts. We find a progressive increase in zircon homogeneity for rhyolite eruptions from ~260 to 75 ka, with homogeneous low-δ¹⁸O zircon values of 2.7–2.8‰ that are in equilibrium with low-δ¹⁸O host melts for the majority of the youngest eruptions. New sanidine Pb isotope data define separate arrays for post-caldera rhyolites and the caldera-forming tuffs that preceded them, indicating that they were not sourced from a mushy Lava Creek Tuff batholith that remained after caldera collapse. Rather, our new age and isotopic data indicate that the post-caldera rhyolites were generated by remelting of a variety of intracaldera source rocks, consisting of pre-Lava Creek Tuff volcanic and plutonic rocks and earlier erupted post-Lava Creek Tuff rhyolites. Batch assembly of low-δ¹⁸O melts starting at ~260 ka resulted in progressive homogenization, followed by differentiation and cooling up until the last rhyolite eruption ~75 ka, a trend that we interpret to be characteristic of a dying magma reservoir beneath the Yellowstone caldera.     

Pacific Basin tsunami hazards associated with mass flows in the Aleutian arc of Alaska

We analyze mass-flow tsunami generation for selected areas within the Aleutian arc of Alaska using results from numerical simulation of hypothetical but plausible mass-flow sources such as submarine landslides and volcanic debris avalanches. The Aleutian arc consists of a chain of volcanic mountains, volcanic islands, and submarine canyons, surrounded by a low-relief continental shelf above about 1000–2000 m water depth. Parts of the arc are fragmented into a series of fault-bounded blocks, tens to hundreds of kilometers in length, and separated from one another by distinctive fault-controlled canyons that are roughly normal to the arc axis. The canyons are natural regions for the accumulation and conveyance of sediment derived from glacial and volcanic processes. The volcanic islands in the region include a number of historically active volcanoes and some possess geological evidence for large-scale sector collapse into the sea. Large scale mass-flow deposits have not been mapped on the seafloor south of the Aleutian Islands, in part because most of the area has never been examined at the resolution required to identify such features, and in part because of the complex nature of erosional and depositional processes. Extensive submarine landslide deposits and debris flows are known on the north side of the arc and are common in similar settings elsewhere and thus they likely exist on the trench slope south of the Aleutian Islands. Because the Aleutian arc is surrounded by deep, open ocean, mass flows of unconsolidated debris that originate either as submarine landslides or as volcanic debris avalanches entering the sea may be potential tsunami sources.To test this hypothesis we present a series of numerical simulations of submarine mass-flow initiated tsunamis from eight different source areas. We consider four submarine mass flows originating in submarine canyons and four flows that evolve from submarine landslides on the trench slope. The flows have lengths that range from 40 to 80 km, maximum thicknesses of 400–800 m, and maximum widths of 10–40 km. We also evaluate tsunami generation by volcanic debris avalanches associated with flank collapse, at four locations (Makushin, Cleveland, Seguam and Yunaska SW volcanoes), which represent large to moderate sized events in this region. We calculate tsunami sources using the numerical model TOPICS and simulate wave propagation across the Pacific using a spherical Boussinesq model, which is a modified version of the public domain code FUNWAVE. Our numerical simulations indicate that geologically plausible mass flows originating in the North Pacific near the Aleutian Islands can indeed generate large local tsunamis as well as large transoceanic tsunamis. These waves may be several meters in elevation at distal locations, such as Japan, Hawaii, and along the North and South American coastlines where they would constitute significant hazards.     

The spin evolution of nascent neutron stars

The loss of angular momentum owing to unstable r-modes in hot young neutron stars has been proposed as a mechanism for achieving the spin rates inferred for young pulsars. One factor that could have a significant effect on the action of the r-mode instability is fallback of supernova remnant material. The associated accretion torque could potentially counteract any gravitational-wave-induced spin-down, and accretion heating could affect the viscous damping rates and hence the instability. We discuss the effects of various external agents on the r-mode instability scenario within a simple model of supernova fallback on to a hot young magnetized neutron star. We find that the outcome depends strongly on the strength of the magnetic field of the star. Our model is capable of generating spin rates for young neutron stars that accord well with initial spin rates inferred from pulsar observations. The combined action of r-mode instability and fallback appears to cause the spin rates of neutron stars born with very different spin rates to converge, on a time-scale of approximately 1 year. The results suggest that stars with magnetic fields ≤10¹³ G could emit a detectable gravitational wave signal for perhaps several years after the supernova event. Stars with higher fields (magnetars) are unlikely to emit a detectable gravitational wave signal via the r-mode instability. The model also suggests that the r-mode instability could be extremely effective in preventing young neutron stars from going dynamically unstable to the bar-mode.     

The energy dependence of burst oscillations from the accreting millisecond pulsar XTE J1814−338

The nature of the asymmetry that gives rise to Type I X-ray burst oscillations on accreting neutron stars remains a matter of debate. Of particular interest is whether the burst oscillation mechanism differs between the bursting millisecond pulsars and the non-pulsing systems. One means to diagnose this is to study the energy dependence of the burst oscillations: here we present an analysis of oscillations from 28 bursts observed during the 2003 outburst of the accreting millisecond pulsar XTE J1814−338. We find that the fractional amplitude of the burst oscillations falls with energy, in contrast to the behaviour found by Muno et al. in the burst oscillations from a set of non-pulsing systems. The drop with energy mirrors that seen in the accretion-powered pulsations; in this respect XTE J1814−338 behaves like the other accreting millisecond pulsars. The burst oscillations show no evidence for either hard or soft lags, in contrast to the persistent pulsations, which show soft lags of up to 50 μs. The fall in amplitude with energy is inconsistent with current surface-mode and simple hotspot models of burst oscillations. We discuss improvements to the models and uncertainties in the physics that might resolve these issues.     

Geochemical mapping using stream sediments in west-central Nigeria: Implications for environmental studies and mineral exploration in West Africa

This paper provides an overview of regional geochemical mapping using stream sediments from central and south-western Nigeria. A total of 1569 stream sediment samples were collected and 54 major and trace elements determined by ICP-MS and Au, Pd and Pt by fire assay. Multivariate statistical techniques (e.g., correlation analysis and principal factor analysis) were used to explore the data, following appropriate data transformation, to understand the data structure, investigate underlying processes controlling spatial geochemical variability and identify element associations. Major geochemical variations are controlled by source geology and provenance, as well as chemical weathering and winnowing processes, more subtle variations are a result of land use and contamination from anthropogenic activity.