A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

Rate of growth of the preserved North American continental crust: Evidence from Hf and O isotopes in Mississippi detrital zircons

Detrital zircons from the Mississippi River have been analyzed for U-Th-Pb, Lu-Hf and O isotopes to constrain the rate of growth of the preserved North American continental crust. One hundred and forty two concordant zircon U/Pb dates on grains mounted in epoxy, obtained by Excimer laser ablation ICP-MS method, resolved six major periods of zircon crystallization: 0-0.25, 0.3-0.6, 0.95-1.25, 1.3-1.5, 1.65-1.95 and 2.5-3.0 Ga. These age ranges match the ages of the recognized lithotectonic units of the North American continent in the hinterland of the Mississippi River. Ninety-six zircons mounted on tape, which show no age zonation and were within 7.5% of concordance, were selected to represent the six U/Pb age time intervals and analyzed for Lu-Hf and O isotope by laser ablation MC-ICP-MS and SHRIMP II, respectively. The δ¹⁸O values of the zircons show a small step increase in the maximum δ¹⁸O values at the Archean-Proterozoic boundary from 7.5‰ in the Archean to 9.5‰, and rarely 13‰, in the Proterozoic and Phanerozoic. However, the average value of δ¹⁸O in zircons changes little with time, showing that the increase in the maximum δ¹⁸O values between 2.5 and 2.0 Ga, which can be attributed to an increase in the sediment content of the source regions of younger granitoids, is largely balanced by an increase in zircons with anomalously low δ¹⁸O, which can be attributed to hydrothermally altered crust in the granitoid source region.εHf_i values for the zircons range from 13.1 to −26.9. Zircons derived from juvenile crust, which we define as having mantle δ¹⁸O (4.5-6.5‰) and lying within error of the Hf depleted mantle growth curve, are rare or absent in the Mississippi basin. The overwhelming majority of zircons crystallized from melted pre-existing continental crust, or mantle-derived magmas that were contaminated by continental crust. The average time difference between primitive crust formation and remelting for each of the recognized lithotectonic time intervals, which is defined as crustal incubation time in this study, is 890 ± 460 Myr. There is also a suggestion that the crustal incubation time increases with decreasing age in the Mississippi basin, which is consistent with the declining role of radioactive heat production in the lower crust with time.The average Hf model age (1.94 Ga), weighted by fraction of zircons in the river load is in reasonable agreement with the Nd model age (1.7 Ga) for the Mississippi River. However, if the zircons are weighted by the area of North America covered by the six recognized periods of zircon crystallization the average model age is 2.35 Ga, which compares favorably with an area weighted Nd model age of 2.36 Ga. Our preferred approach is to use the measured O isotope values to constrain variations in the ¹⁷⁶Lu/¹⁷⁷Hf ratio of the granitic source region from which the zircons crystallized, making the assumption that zircons with mantle-like O isotopic ratios have higher ¹⁷⁶Lu/¹⁷⁷Hf than zircons with higher O isotope values. This method gives an average Hf model age of 2.53 Ga, which is 180 Myr older than the constant ¹⁷⁶Lu/¹⁷⁷Hf calculation.The area weighted zircon Hf model ages show two distinct periods of crust formation for the North American continent, 1.6-2.2 and 2.9-3.4 Ga. At least 50% of the preserved North American continental crust was extracted from the mantle by 2.9 Ga and 90% by 1.6 Ga. Two similar periods of crustal growth are also recognized in Gondwana (Hawkesworth C. J. and Kemp A. I. S. (2006) Using hafnium and oxygen isotopes in zircons to unravel the record of crustal evolution. Chem. Geol.226, 144-162.), suggesting that these may be periods of global continental crustal growth. However, we stress that more data from other continents are required before the hypothesis of episodic global continental growth can be accepted with confidence.     

Examining potential genetic links between Jurassic porphyry Cu–Au ± Mo and epithermal Au ± Ag mineralization in the Toodoggone district of North-Central British Columbia,Canada

The Toodoggone district comprises Upper Triassic to Lower Jurassic Hazelton Group Toodoggone Formation volcanic and sedimentary rocks, which unconformably overlie submarine island-arc volcanic and sedimentary rocks of the Lower Permian Asitka Group and Middle Triassic Takla Group, some of which are intruded by Upper Triassic to Lower Jurassic plutons and dikes of the Black Lake suite. Although plutonism occurred episodically from ca. 218 to 191 Ma, the largest porphyry Cu–Au ± Mo systems formed from ca. 202 to 197 Ma, with minor mineralization occurring from ca. 197 to 194 Ma. Porphyry-style mineralization is hosted by small-volume (<1 km³), single-phase, porphyritic igneous stocks or dikes that have high-K calc-alkaline compositions and are comparable with volcanic-arc granites. The Fin porphyry Cu–Au–Mo deposit is anomalous in that it is 16 m.y. older than any other porphyry Cu–Au ± Mo occurrence in the district and has lower REEs. All porphyry systems are spatially restricted to exposed Asitka and Takla Group basement rocks, and rarely, the lowest member of the Hazelton Group (i.e., the ca. 201 Ma Duncan Member). The basement rocks to intrusions are best exposed in the southern half of the district, where high rates of erosion and uplift have resulted in their preferential exposure. In contrast, low- and high-sulfidation epithermal systems are more numerous in the northern half of the district, where the overlying Hazelton Group rocks dominate exposures. Cogenetic porphyry systems might also exist in the northern areas; however, if they are present, they are likely to be buried deeply beneath Hazelton Group rocks. High-sulfidation epithermal systems formed at ca. 201 to 182 Ma, whereas low-sulfidation systems were active at ca. 192 to 162 Ma. Amongst the studied epithermal systems, the Baker low-sulfidation epithermal deposit displays the strongest demonstrable genetic link with magmatic fluids; fluid inclusion studies demonstrate that its ore fluids were hot (>468°C), saline, and deposited metals at deep crustal depths (>2 km). Sulfur, C, O, and Pb isotope data confirm the involvement of a magmatic fluid, but also suggest that the ore fluid interacted with Asitka and Takla Group country rocks prior to metal deposition. In contrast, in the Shasta, Lawyers, and Griz-Sickle low-sulfidation epithermal systems, there is no clear association with magmatic fluids. Instead, their fluid inclusion data indicate the involvement of low-temperature (175 to 335°C), low-salinity (1 to 11 equiv. wt.% NaCl) fluids that deposited metals at shallow depths (<850 m). Their isotope (i.e., O, H, Pb) data suggest interaction between meteoric and/or metamorphic ore fluids with basement country rocks.     

Composition and origin of the Çaldağ oxide nickel laterite,W. Turkey

The Çalda? nickel laterite deposit located in the Aegean region of W. Turkey contains a reserve of 33 million tons of Ni ore with an average grade of 1.14% Ni. The deposit is developed on an ophiolitic serpentinite body which was obducted onto Triassic dolomites in the Late Cretaceous. The deposit weathering profile is both laterally and vertically variable. A limonite zone, which is the main ore horizon, is located at the base of the profile. A hematite horizon is located above the limonite, which in the south of the deposit is capped by Eocene freshwater limestones and in the north by a siliceous horizon. The deposit is unusual in lacking a significant saprolite zone with little development of Ni-silicates. The boundary between the limonite zone and serpentinite below is sharp with a marked decrease in concentrations of MgO from 13 to 1 wt.% over a distance of 2 mm representing the ‘Mg discontinuity’. Ni concentrations within goethite, the main ore mineral, reach a maximum of ~3 wt.% near the base of the limonite zone. Silica concentrations are high throughout most of the laterite with up to 80 wt.% silica in the upper portion of some profiles. The combination of a serpentinite protolith and a high water table at Çalda?, in association with an aggressive weathering environment in a tropical climate, resulted in the formation of an oxide-dominated deposit. The precipitation of silica may coincide with a change in climate with silica precipitation linked to an increase in seasonality. Additional variations within profile morphology are attributed to transportation during and after laterite development as a result of faulting, pocket type laterite formation and slumping, each of which produces a contrasting set of textural and geochemical features.     

Adsorption kinetics of CO2, CH4, and their equimolar mixture on coal from the Black Warrior Basin,West-Central Alabama

Laboratory experiments were conducted to investigate the adsorption kinetic behavior of pure and mixed gases (CO_2, CH₄, approximately equimolar CO₂ + CH₄ mixtures, and He) on a coal sample obtained from the Black Warrior Basin at the Littleton Mine (Twin Pine Coal Company), Jefferson County, west-central Alabama. The sample was from the Mary Lee coal zone of the Pottsville Formation (Lower Pennsylvanian). Experiments with three size fractions (45–150 µm, 1–2 mm, and 5–10 mm) of crushed coal were performed at 40 °C and 35 °C over a pressure range of 1.4–6.9 MPa to simulate coalbed methane reservoir conditions in the Black Warrior Basin and provide data relevant for enhanced coalbed methane recovery operations. The following key observations were made: (1) CO₂ adsorption on both dry and water-saturated coal is much more rapid than CH₄ adsorption; (2) water saturation decreases the rates of CO₂ and CH₄ adsorption on coal surfaces, but it appears to have minimal effects on the final magnitude of CO₂ or CH₄ adsorption if the coal is not previously exposed to CO₂; (3) retention of adsorbed CO₂ on coal surfaces is significant even with extreme pressure cycling; and (4) adsorption is significantly faster for the 45–150 μm size fraction compared to the two coarser fractions.     

Continental shelf nurseries and recruitment variability in American plaice and yellowtail flounder on the Grand Bank: insights into stock resiliency

In 1994, a directed fishing moratorium was declared for Grand Bank American plaice (Hippoglossoides platessoides) and yellowtail flounder (Limanda ferruginea) stocks because both stocks showed severe declines in abundance from heavy exploitation during the mid 1980s and early 1990s. Four years later, the fishery for yellowtail re-opened while the plaice stock has shown little recovery and the moratorium is still in effect. To assess the possible causes of the differences in recovery between species, we examined the spatial structure and environmental characteristics of the continental shelf nursery habitats of plaice and yellowtail, and their relationship to recruitment variability and overall population size. Depth plays a major influential role determining the spatial pattern and the abundance of juveniles of both species and in the case of plaice the spatial structure of the adult population also determines the amount of nursery area utilised by juveniles. Recruitment variability was higher in plaice than in yellowtail. We found year class synchrony in both species indicating that common environmental conditions and/or biological processes are affecting recruitment in a similar manner. Density-dependent regulation appears to be more severe in yellowtail and this should contribute to a more stable population when compared to plaice. These results are discussed in terms of resiliency of both stocks to over-exploitation.     

Relating millimeter-scale turbulence to meter-scale subtidal erosion and accretion across the fringe of a coastal mangrove forest

Within a wave-exposed mangrove forest, novel field observations are presented, comparing millimeter-scale turbulent water velocity fluctuations with contemporaneous subtidal bed elevation changes. High-resolution velocity and bed level measurements were collected from the unvegetated mudflat, at the mangrove forest fringe, and within the forest interior over multiple tidal cycles (flood–ebb) during a 2-week period. Measurements demonstrated that the spatial variability in vegetation density is a control on sediment transport at sub-meter scales. Scour around single and dense clusters of pneumatophores was predicted by a standard hydraulic engineering equation for wave-induced scour around regular cylinders, when the cylinder diameter in the equations was replaced with the representative diameter of the dense pneumatophore clusters. Waves were dissipated as they propagated into the forest, but dissipation at infragravity periods (> 30 s) was observed to be less than dissipation at shorter periods (< 30 s), consistent with the predictions of a simple model. Cross-wavelet analysis revealed that infragravity-frequency fluctuations in the bed level were occasionally coherent with velocity, possibly indicating scour upstream of dense pneumatophore patches when infragravity waves reinforced tidal currents. Consequently, infragravity waves were a likely driver of sediment transport within the mangrove forest. Near-bed turbulent kinetic energy, estimated from the turbulent dissipation rate, was also correlated with bed level changes. Specifically, within the mangrove forest and over the unvegetated mudflat, high-energy events were associated with erosion or near-zero bed level change, whereas low-energy events were associated with accretion. In contrast, no single relationship between bed level changes and mean current velocity was applicable across both vegetated and unvegetated regions. These observations support the theory that sediment mobilization scales with turbulent energy, rather than mean velocity, a distinction that becomes important when vegetation controls the development of turbulence.     

The statistical properties of Λ cold dark matter halo formation

We present a comparison of the statistical properties of dark matter halo merger trees extracted from the Millennium Simulation with Extended Press–Schechter (EPS) formalism and the related galform Monte Carlo method for generating ensembles of merger trees. The volume, mass resolution and output frequency make the Millennium Simulation a unique resource for the study of the hierarchical growth of structure. We construct the merger trees of present-day friends-of-friends groups and calculate a variety of statistics that quantify the masses of their progenitors as a function of redshift, accretion rates, and the redshift distribution of their most recent major merger. We also look in the forward direction and quantify the present-day mass distribution of haloes into which high-redshift progenitors of a specific mass become incorporated. We find that the EPS formalism and its Monte Carlo extension capture the qualitative behaviour of all these statistics, but as redshift increases they systematically underestimate the masses of the most massive progenitors. This shortcoming is worst for the Monte Carlo algorithm. We present a fitting function to a scaled version of the progenitor mass distribution and show how it can be used to make more accurate predictions of both progenitor and final halo mass distributions.     

Southern high latitude climate variability in the Late Cretaceous greenhouse world

A palynological study of oil exploration wells in the Gippsland Basin southeastern Australia has provided a record of southern high latitude climate variability for the last 12 million years of the Cretaceous greenhouse world. During this time, the vegetation was dominated by a cool to temperate flora of Podocarpaceae, Proteaceae and Nothofagidites spp. at a latitude of 60°S. Milankovitch forced cyclic alternations from drier to wetter climatic periods caused vegetation variability from 72 to 77 Ma. This climate change was probably related to the waxing and waning of ephemeral (100 ky) small ice sheets in Antarctica during times of insolation minima and maxima. Drying and cooling after 72 Ma culminated from 68 to 66 Ma, mirroring trends in global δ¹⁸O data. Quantitative palynofloral analyses have the potential to provide realistic proxies for small-scale climate variability in the predominantly ice-free Late Cretaceous.