High-performance direct gravitational N-body simulations on graphics processing units

We present the results of gravitational direct N-body simulations using the commercial graphics processing units (GPU) NVIDIA Quadro FX1400 and GeForce 8800GTX, and compare the results with GRAPE-6Af special purpose hardware. The force evaluation of the N-body problem was implemented in Cg using the GPU directly to speed-up the calculations. The integration of the equations of motions were, running on the host computer, implemented in C using the 4th order predictor–corrector Hermite integrator with block time steps. We find that for a large number of particles (N  10⁴) modern graphics processing units offer an attractive low cost alternative to GRAPE special purpose hardware. A modern GPU continues to give a relatively flat scaling with the number of particles, comparable to that of the GRAPE. The GRAPE is designed to reach double precision, whereas the GPU is intrinsically single-precision. For relatively large time steps, the total energy of the N-body system was conserved better than to one in 10⁶ on the GPU, which is impressive given the single-precision nature of the GPU. For the same time steps, the GRAPE gave somewhat more accurate results, by about an order of magnitude. However, smaller time steps allowed more energy accuracy on the grape, around 10⁻¹¹, whereas for the GPU machine precision saturates around 10⁻⁶ For N  10⁶ the GeForce 8800GTX was about 20 times faster than the host computer. Though still about a factor of a few slower than GRAPE, modern GPUs outperform GRAPE in their low cost, long mean time between failure and the much larger onboard memory; the GRAPE-6Af holds at most 256k particles whereas the GeForce 8800GTX can hold 9 million particles in memory.     

Impact of dynamic feedbacks between sedimentation,sea-level rise,and biomass production on near-surface marsh stratigraphy and carbon accumulation

Salt marshes accrete both organic and inorganic sediments. Here we present analytical and numerical models of salt marsh sedimentation that, in addition to capturing inorganic processes, explicitly account for above- and belowground organic processes including root growth and decay of organic carbon. The analytical model is used to examine the bias introduced by organic processes into proxy records of sedimentation, namely ¹³⁷Cs and ²¹⁰Pb. We find that accretion rates estimated using ²¹⁰Pb will be less than accretion rates estimated using the ¹³⁷Cs peak in steadily accreting marshes if (1) carbon decay is significant and (2) data for ²¹⁰Pb extend below the ¹³⁷Cs peak. The numerical model expands upon the analytical model by including belowground processes such as compaction and root growth, and by explicitly tracking the evolution of aboveground biomass and its effect on sedimentation rates. Using the numerical model we explore how marsh stratigraphy responds to sediment supply and the rate of sea-level rise. It is calibrated and tested using an extensive data set of both marsh stratigraphy and measurements of vegetation dynamics in a Spartina alterniflora marsh in South Carolina, USA. We find that carbon accumulation in marshes is nonlinearly related to both the supply of inorganic sediment and the rate of sea-level rise; carbon accumulation increases with sea-level rise until sea-level rise reaches a critical rate that drowns the marsh vegetation and halts carbon accumulation. The model predicts that changes in carbon storage resulting from changing sediment supply or sea-level rise are strongly dependent on the background sediment supply: if inorganic sediment supply is reduced in an already sediment poor marsh the storage of organic carbon will increase to a far greater extent than in a sediment-rich marsh, provided that the rate of sea-level rise does not exceed a threshold. These results imply that altering sediment supply to estuaries (e.g., by damming upstream rivers or altering littoral sediment transport) could lead to significant changes in the carbon budgets of coastal salt marshes.     

Environmental correlates of large-scale spatial variation in the δ13C of marine animals

Carbon stable isotopes can be used to trace the sources of energy supporting food chains and to estimate the contribution of different sources to a consumer's diet. However, the δ¹³C signature of a consumer is not sufficient to infer source without an appropriate isotopic baseline, because there is no way to determine if differences in consumer δ¹³C reflect source changes or baseline variation. Describing isotopic baselines is a considerable challenge when applying stable isotope techniques at large spatial scales and/or to interconnected food chains in open marine environments. One approach is to use filter-feeding consumers to integrate the high frequency and small-scale variation in the isotopic signature of phytoplankton and provide a surrogate baseline, but it can be difficult to sample a single consumer species at large spatial scales owing to rarity and/or discontinuous distribution. Here, we use the isotopic signature of a widely distributed filter-feeder (the queen scallop Aequipecten opercularis) in the north-eastern Atlantic to develop a model linking base δ¹³C to environmental variables. Remarkably, a single variable model based on bottom temperature has good predictive power and predicts scallop δ¹³C with mean error of only 0.6‰ (3%). When the model was used to predict an isotopic baseline in parts of the overall study region where scallop were not consistently sampled, the model accounted for 76% and 79% of the large-scale spatial variability (10¹–10⁴ km) of the δ¹³C of two fish species (dab Limanda limanda and whiting Merlangus merlangius) and 44% of the δ¹³C variability in a mixed fish community. The results show that source studies would be significantly biased if a single baseline were applied to food webs at larger scales. Further, when baseline δ¹³C cannot be directly measured, a calculated baseline value can eliminate a large proportion of the unexplained variation in δ¹³C at higher trophic levels.     

The Monitor project: JW 380 – a 0.26-, 0.15-M⊙, pre-main-sequence eclipsing binary in the Orion nebula cluster

We present a general recipe for constructing N -body realizations of galaxies comprising near spherical and disc components. First, an exact spherical distribution function for the spheroids (halo and bulge) is determined, such that it is in equilibrium with the gravitational monopole of the disc components. Second, an N -body realization of this model is adapted to the full disc potential by growing the latter adiabatically from its monopole. Finally, the disc is sampled with particles drawn from an appropriate distribution function, avoiding local-Maxwellian approximations. We performed test simulations and find that the halo and bulge radial density profile very closely match their target model, while they become slightly oblate due to the added disc gravity. Our findings suggest that vertical thickening of the initially thin disc is caused predominantly by spiral and bar instabilities, which also result in a radial re-distribution of matter, rather than scattering off interloping massive halo particles.     

Fine‐scale diet of the Australian sea lion (Neophoca cinerea) using DNA‐based analysis of faeces

We applied DNA‐based faecal analysis to determine the diet of female Australian sea lions (n = 12) from two breeding colonies in South Australia. DNA dietary components of fish and cephalopods were amplified using the polymerase chain reaction and mitochondrial DNA primers targeting the short (~100 base pair) section of the 16S gene region. Prey diversity was determined by sequencing ~50 amplicons generated from clone libraries developed for each individual. Faecal DNA was also combined and cloned from multiple individuals at each colony and fish diversity determined. Diets varied between individuals and sites. Overall, DNA analysis identified a broad diversity of prey comprising 23 fish and five cephalopod taxa, including many species not previously described as prey of the Australian sea lion. Labridae (wrasse), Monacanthidae (leatherjackets) and Mullidae (goat fish) were important fish prey taxa. Commonly identified cephalopods were Octopodidae (octopus), Loliginidae (calamary squid) and Sepiidae (cuttlefish). Comparisons of fish prey diversity determined by pooling faecal DNA from several samples provided a reasonable but incomplete resemblance (55–71%) to the total fish diversity identified across individual diets at each site. Interpretation of diet based on the recovery of prey hard‐parts identified one cephalopod beak (Octopus sp.) and one fish otolith (Parapriacanthus elongatus). The present study highlights the value of DNA‐based analyses and their capabilities to enhance information of trophic interactions.