A new parallel PM code for very large-scale cosmological simulations

We have developed a parallel Particle–Particle, Particle–Mesh (P³M) simulation code for the Cray T3E parallel supercomputer that is well suited to studying the time evolution of systems of particles interacting via gravity and gas forces in cosmological contexts. The parallel code is based upon the public-domain serial Adaptive P³M-SPH (http://coho.astro.uwo.ca/pub/hydra/hydra.html) code of Couchman et al. (1995)[ApJ, 452, 797]. The algorithm resolves gravitational forces into a long-range component computed by discretizing the mass distribution and solving Poisson's equation on a grid using an FFT convolution method, and a short-range component computed by direct force summation for sufficiently close particle pairs. The code consists primarily of a particle–particle computation parallelized by domain decomposition over blocks of neighbour-cells, a more regular mesh calculation distributed in planes along one dimension, and several transformations between the two distributions. The load balancing of the P³M code is static, since this greatly aids the ongoing implementation of parallel adaptive refinements of the particle and mesh systems. Great care was taken throughout to make optimal use of the available memory, so that a version of the current implementation has been used to simulate systems of up to 10⁹ particles with a 1024³ mesh for the long-range force computation. These are the largest Cosmological N-body simulations of which we are aware. We discuss these memory optimizations as well as those motivated by computational performance. Performance results are very encouraging, and, even without refinements, the code has been used effectively for simulations in which the particle distribution becomes highly clustered as well as for other non-uniform systems of astrophysical interest.     

Petrographic investigation of shatter cone melt films recovered from MEMIN impact experiments in sandstone and iSALE modeling of their formation boundary conditions

Shatter cones are diagnostic for the recognition of meteorite impact craters. They are unambiguously identifiable in the field and the only macroscopic shock deformation feature. However, the physical boundary conditions and exact formation mechanism(s) are still a subject of debate. Melt films found on shatter cone surfaces may allow the constraint of pressure–temperature conditions during or immediately after their formation. Within the framework of the MEMIN research group, we recovered 24 shatter cone fragments from the ejecta of hypervelocity impact experiments. Here, we focus on silicate melt films (now quenched to glass) found on shatter cone surfaces formed in experiments with 20–80 cm sized sandstone targets, impacted by aluminum and iron meteorite projectiles of 5 and 12 mm diameter at velocities of 7.0 and 4.6 km s⁻¹, respectively. The recovered shatter cone fragments vary in size from 1.2 to 9.3 mm. They show slightly curved, striated surfaces, and conical geometries with apical angles of 36°–52°. The fragments were recovered from experiments with peak pressures ranging from 46 to 86 GPa, and emanated from a zone within 0.38 crater radii. Based on iSale modeling and petrographic investigations, the shatter coned material experienced low bulk shock pressures of 0.5–5 GPa, whereas deformation shows a steep increase toward the shatter cone surface leading to localized melting of the rock, resulting in both vesicular as well as polished melt textures visible under the SEM. Subjacent to the melt films are zones of fragmentation and brittle shear, indicating movement away from the shatter cone apex of the rock that surrounds the cone. Smearing and extension of the melt film indicates subsequent movement in opposite direction to the comminuted and brecciated shear zone. We believe the documented shear textures and the adjacent smooth melt films can be related to frictional melting, whereas the overlying highly vesiculated melt layer could indicate rapid pressure release. From the observation of melting and mixing of quartz, phyllosilicates, and rutile in this overlying texture, we infer high, but very localized postshock temperatures exceeding 2000 °C. The melted upper part of the shatter cone surface cross-cuts the fragmented lower section, and is accompanied by PDFs developed in quartz parallel to the {112} plane. Based on the overprinting textures and documented shock effects, we hypothesize shatter cones start to form during shock loading and remain an active fracture surface until pressure release during unloading and infer that shatter cone surfaces are mixed mode I/II fracture surfaces.     

Einfluß der TTC-Konzentration auf die Dehydrogenasenaktivität von Mikroorganismen

Die aus dem Sediment eines Gewässers isolierte Mischpopulation von Bakterien wurde auf Nährbouillon kultiviert sowie über unterschiedlich lange Zeit an Glucose und Phenol adaptiert. Die zur Bestimmung der Dehydrogenasenaktivität optimale Konzentration des TTC lag für die verschiedenen Kulturen zwischen 0,04 und 0,2 Volumenprozent der Bakterien-Suspension bzw. zwischen 0,1 bzw. 1 mg TTC/mg TS · h bezogen auf die getrocknete Bakterienbiomasse bei einstündiger Inkubation. Die Unterschiede in der Formazanbildung wie in der toxischen Wirkung des TTC sind abhängig von der Struktur der einzelnen Bakterienzelle wie von der Struktur der Belebtschlammflocken und dem ,,Alter“ der Kultur, weniger abhängig vom Substrat.     

Postglacial relative sea-level change and stratigraphy of raised coastal basins on Kola Peninsula, northwest Russia

A relative sea-level curve for the Holocene is constructed for Polyarny on the Kola Peninsula, northwest Russia. The curve is based on 18 radiocarbon dates of isolation contacts, identified from lithological and diatomological criteria, in nine lake basins situated between 12 and 57 m a.s.l. Most of the lakes show a conformable, regressive I–II–III (marine–transitional–freshwater) facies succession, indicating a postglacial history comprising an early (10,000–9000 radiocarbon years BP) phase of rapid, glacio-isostatically induced emergence (5 cm year⁻¹) and a later phase (after 7000 years BP,) having a moderate rate of emergence (<0.5 cm year⁻¹). Three lakes together record a phase of very low rate of emergence or slight sea-level rise at a level of 27 m a.s.l., between 8500 and 7000 years BP, which correlates with the regional Tapes transgression. Pollen stratigraphy in the highest lake shows that the area was deglaciated before the Younger Dryas and that previously reconstructed Younger Dryas glacier margins along the north Kola coast lie too far north.     

Using the anchoring device of a comet lander to determine surface mechanical properties

Owing to the low surface gravity of the Rosetta target comet 46P/Wirtanen, a means of anchoring the Rosetta Lander to the cometary surface will be necessary. This task can be accomplished by firing an anchor into the cometary soil immediately after touchdown to prevent a rebound of the spacecraft from the surface or subsequent ejection by other forces, and to allow for mechanical activities (drilling, etc.) at the landing site.

The rationale for anchoring is examined, based on estimates of the main forces likely to act on the spacecraft after landing. We report on the development of an anchoring device using a pyrotechnic gas generator as a power source and an instrumented anchor.

In addition to the anchoring function, which is the primary purpose of this system, the integration of acceleration and temperature sensors into the tip offers the possibility to determine some important material properties of the cometary surface layer. The accelerometer is designed to measure the deceleration history of the projectile and is thus expected to give information on how the material properties (in particular strength) change within the penetrated layer(s), while the temperature sensor will measure temperature variations at the depth at which the anchor finally comes to rest. As the mechanical properties of the material are not known, it is difficult to predict the final depth of the anchor with any great certainty, but it may well be greater than that reached by any other of the lander's instruments.

The instrumented anchor will be part of the MUPUS experiment, selected to form part of the Rosetta Lander payload. We report on results of laboratory simulations of anchor penetration performed at the Institut für Weltraumforschung, Graz, and compare these with models of projectile penetration. The value of the results expected from the penetrometry experiment in the context of an improved understanding of cometary processes is discussed.     

Tektonische Untersuchungen im Gebiet der Eisenspatgrube „Pfannenberger Einigkeit“ im Siegerland

Ohne Zusammenfassung     

Review of Acquisition and Signal Processing Methods for Electromagnetic Noise Reduction and Retrieval of Surface Nuclear Magnetic Resonance Parameters

Surveys in Geophysics - Surface nuclear magnetic resonance (sNMR) is an electromagnetic hydrogeophysical method directly sensitive to liquid phase water in the upper $$\approx $$ 100 m of the...     

A high-quality monthly total cloud amount dataset for Australia

A high-quality monthly total cloud amount dataset for 165 stations has been developed for monitoring and assessing long-term trends in cloud cover over Australia. The dataset is based on visual 9 a.m. and 3 p.m. observations of total cloud amount, with most records starting around 1957. The quality control process involved examination of historical station metadata, together with an objective statistical test comparing candidate and reference cloud series. Individual cloud series were also compared against rainfall and diurnal temperature range series from the same site, and individual cloud series from neighboring sites. Adjustments for inhomogeneities caused by relocations and changes in observers were applied, as well as adjustments for biases caused by the shift to daylight saving time in the summer months. Analysis of these data reveals that the Australian mean annual total cloud amount is characterised by high year-to-year variability and shows a weak, statistically non-significant increase over the 1957–2007 period. A more pronounced, but also non-significant, decrease from 1977 to 2007 is evident. A strong positive correlation is found between all-Australian averages of cloud amount and rainfall, while a strong negative correlation is found between mean cloud amount and diurnal temperature range. Patterns of annual and seasonal trends in cloud amount are in general agreement with rainfall changes across Australia, however the high-quality cloud network is too coarse to fully capture topographic influences. Nevertheless, the broadscale consistency between patterns of cloud and rainfall variations indicates that the new total cloud amount dataset is able to adequately describe the broadscale patterns of change over Australia. Favourable simple comparisons between surface and satellite measures of cloudiness suggest that satellites may ultimately provide the means for monitoring long-term changes in cloud over Australia. However, due to the relative shortness and homogeneity problems of the satellite record, a robust network of surface cloud observations will be required for many years to come.