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.     

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.     

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.     

Nanometer‐scale anatomy of entire Stardust tracks

Abstract– We have developed new sample preparation and analytical techniques tailored for entire aerogel tracks of Wild 2 sample analyses both on “carrot” and “bulbous” tracks. We have successfully ultramicrotomed an entire track along its axis while preserving its original shape. This innovation allowed us to examine the distribution of fragments along the entire track from the entrance hole all the way to the terminal particle. The crystalline silicates we measured have Mg‐rich compositions and O isotopic compositions in the range of meteoritic materials, implying that they originated in the inner solar system. The terminal particle of the carrot track is a ¹⁶O‐rich forsteritic grain that may have formed in a similar environment as Ca‐, Al‐rich inclusions and amoeboid olivine aggregates in primitive carbonaceous chondrites. The track also contains submicron‐sized diamond grains likely formed in the solar system. Complex aromatic hydrocarbons distributed along aerogel tracks and in terminal particles. These organics are likely cometary but affected by shock heating.     

The initial responses of hot liquid water released under low atmospheric pressures: Experimental insights

Experiments have been performed to simulate the shallow ascent and surface release of water and brines under low atmospheric pressure. Atmospheric pressure was treated as an independent variable and water temperature and vapor pressure were examined as a function of total pressure variation down to low pressures. The physical and thermal responses of water to reducing pressure were monitored with pressure transducers, temperature sensors and visible imaging. Data were obtained for pure water and for solutions with dissolved NaCl or CO₂. The experiments showed the pressure conditions under which the water remained liquid, underwent a rapid phase change to the gas state by boiling, and then solidified because of removal of latent heat. Liquid water is removed from phase equilibrium by decompression. Solid, liquid and gaseous water are present simultaneously, and not at the 611 Pa triple point, because dynamic interactions between the phases maintain unstable temperature gradients. After phase changes stop, the system reverts to equilibrium with its surroundings. Surface and shallow subsurface pressure conditions were simulated for Mars and the icy satellites of the outer Solar System. Freezing by evaporation in the absence of wind on Mars is shown to be unlikely for pure water at pressures greater than c. 670 Pa, and for saline solutions at pressures greater than c. 610 Pa. The physical nature of ice that forms depends on the salt content. Ice formed from saline water at pressures less than c. 610 Pa could be similar to terrestrial sea ice. Ice formed from pure water at pressures less than c. 100 Pa develops a low thermal conductivity and a ‘honeycomb’ structure created by sublimation. This ice could have a density as low as c. 450 kg m⁻³ and a thermal conductivity as low as 1.6 W m⁻¹ K⁻¹, and is highly reflective, more akin to snow than the clear ice from which it grew. The physical properties of ice formed from either pure or saline water at low pressures will act to reduce the surface temperature, and hence rate of sublimation, thereby prolonging the lifespan of any liquid water beneath.     

High outflow channels on Mars indicate Hesperian recharge at low latitudes and the presence of Canyon Lakes

A number of martian outflow channels were carved by discharges from large dilational fault zones. These channels were sourced by groundwater, not surface water, and when observed on high-standing plateaus they provide indicators of elevated paleo-groundwater levels. We identify three outflow channels of Hesperian age that issued from a 750-km-long fault zone extending from Candor Chasma to Ganges Chasma. Two of these channels, Allegheny Vallis and Walla Walla Vallis, have sources >2500 m above the topographic datum, too high to be explained by discharge from a global aquifer that was recharged solely in the south polar region. The indicated groundwater levels likely required regional sources of recharge at low latitudes. The floodwaters that erupted from Ophir Cavus to form Allegheny Vallis encountered two ridges that restricted the flow, forming temporary lakes. The flow probably breached or overtopped these obstructions quickly, catastrophically draining the lakes and carving several scablands. After the last obstacle had been breached, a single main channel formed that captured all subsequent flow. We performed hydrologic analyses of this intermediate phase of the flooding, prior to incision of the channel to its present depth. Using floodwater depths of 30-60 m, we calculated flow velocities of 6-15 m s⁻¹ and discharges in the range of . Locally higher flow velocities and discharges likely occurred when the transient lakes were drained. Variable erosion at the channel and scabland crossing of MOLA pass 10644 suggests that the upper 25-30 m may consist of poorly consolidated surface materials underlain by more cohesive bedrock. We infer that an ice-covered lake with a surface elevation >2500 m probably existed in eastern Candor Chasma because this canyon is intersected by the Ophir Catenae fault system from which Allegheny Vallis and Walla Walla Vallis originated. We introduce a new hydrology concept for Mars in which the groundwater system was augmented by recharge from canyon lakes that were formed when water was released by catastrophic melting of former ice sheets in Tharsis by effusions of flood basalts. This model could help to reconcile the expected presence of a thick cryosphere during the Hesperian with the abundant evidence for groundwater as a source for some of the circum-Chryse outflow channels.     

Chemostratigraphy of the Tamengo Formation (Corumbá Group,Brazil): A contribution to the calibration of the Ediacaran carbon-isotope curve

The Corumbá Group, cropping out in the southern Paraguay Belt in Brazil, is one of the most complete Ediacaran sedimentary archives of palaeogeographic, climatic, biogeochemical and biotic evolution in southwestern Gondwana. The unit hosts a rich fossil record, including acritarchs, vendotaenids (Vendotaenia, Eoholynia), soft-bodied metazoans (Corumbella) and skeletal fossils (Cloudina, Titanotheca). The Tamengo Formation, made up mainly of limestones and marls, provides a rich bio- and chemostratigraphic record. Several outcrops, formerly assigned to the Cuiabá Group, are here included in the Tamengo Formation on the basis of lithological and chemostratigraphical criteria. High-resolution carbon isotopic analyses are reported for the Tamengo Formation, showing (from base to top): (1) a positive δ¹³C excursion to +4‰ PDB above post-glacial negative values, (2) a negative excursion to −3.5‰ associated with a marked regression and subsequent transgression, (3) a positive excursion to +5.5‰, and (4) a plateau characterized by δ¹³C around +3‰. A U-Pb SHRIMP zircon age of an ash bed interbedded in the upper part of the δ¹³C positive plateau yielded 543 ± 3 Ma, which is considered as the depositional age ( Babinski et al., 2008a). The positive plateau in the upper Tamengo Formation and the preceding positive excursion are ubiquitous features in several successions worldwide, including the Nama Group (Namibia), the Dengying Formation (South China) and the Nafun and Ara groups (Oman). This plateau is constrained between 542 and 551 Ma, thus consistent with the age of the upper Tamengo Formation. The negative excursion of the lower Tamengo Formation may be correlated to the Shuram–Wonoka negative anomaly, although δ¹³C values do not fall beyond −3.5‰ in the Brazilian sections. Sedimentary breccias occur just beneath this negative excursion in the lower Tamengo Formation. One possible interpretation of the origin of these breccias is a glacioeustatic sea-level fall, but a tectonic interpretation cannot be completely ruled out.