Simulations of deep pencil-beam redshift surveys

We create mock pencil-beam redshift surveys from very large cosmological N -body simulations of two cold dark matter (CDM) cosmogonies, an Einstein–de Sitter model ( τ CDM) and a flat model with Ω₀=0.3 and a cosmological constant (ΛCDM). We use these to assess the significance of the apparent periodicity discovered by Broadhurst et al. Simulation particles are tagged as 'galaxies' so as to reproduce observed present-day correlations. They are then identified along the past light-cones of hypothetical observers to create mock catalogues with the geometry and the distance distribution of the Broadhurst et al. data. We produce 1936 (2625) quasi-independent catalogues from our τ CDM (ΛCDM) simulation. A couple of large clumps in a catalogue can produce a high peak at low wavenumbers in the corresponding one-dimensional power spectrum, without any apparent large-scale periodicity in the original redshift histogram. Although the simulated redshift histograms frequently display regularly spaced clumps, the spacing of these clumps varies between catalogues and there is no 'preferred' period over our many realizations. We find only a 0.72 (0.49) per cent chance that the highest peak in the power spectrum of a τ CDM (ΛCDM) catalogue has a peak-to-noise ratio higher than that in the Broadhurst et al. data. None of the simulated catalogues with such high peaks shows coherently spaced clumps with a significance as high as that of the real data. We conclude that in CDM universes, the regularity on a scale of ∼130  h ⁻¹ Mpc observed by Broadhurst et al. has a priori probability well below 10⁻³.     

Clustering of galaxy clusters in cold dark matter universes

We use very large cosmological N -body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ( τ CDM) and a flat low density model with a cosmological constant (ΛCDM). Our simulations each use 10⁹ particles to follow the mass distribution within cubes of side 2  h ⁻¹ Gpc ( τ CDM) and 3  h ⁻¹ Gpc (ΛCDM) with a force resolution better than 10⁻⁴ of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from     for samples with mean separation     to     for samples with     The lower value here corresponds to τ CDM and the upper to ΛCDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10 per cent using the simple model of Sheth, Mo & Tormen. This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the automated plate measurement (APM) cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for ΛCDM. The observed power spectrum (for a cluster sample of mean separation     ) lies significantly above the predictions of both models.     

Nonlinear convection in an imposed horizontal magnetic field

Abstract

Nonlinear two-dimensional magnetoconvection, with a Boussinesq fluid driven across the field-lines, is taken as a model for giant-cell convection in the sun and late-type stars. A series of numerical experiments shows the sensitivity of the horizontal scale of convection to the applied field and to the Rayleigh number R. Overstable oscillations occur in cells as broad as they are deep, but increasing R leads to steady motions of much greater wavelength. Purely geometrical effects can cause oscillation: this work implies that strong horizontal field will in general lead to time-dependent convection.     

Structure of an active foreland fold and thrust belt,Papua New Guinea

Geological structure of the active foreland fold and thrust belt of Papua New Guinea has been interpreted using high-quality seismic-reflection data. Three en échelon anticlines, the Strickland, Cecilia and Wai Asi, are located along the frontal margin of the Papuan Fold Belt. All three are foreland-vergent and cut by hinterland-dipping thrust faults that sole into a common detachment beneath the Oligocene to Miocene Darai Limestone. Two of the anticlines are linked by a right-lateral transfer zone. Folding occurs primarily in the upper 2000 m of strata, which consist of Darai Limestone overlain by Miocene to Quaternary siliciclastic sedimentary rocks. Beneath the Darai Limestone lies the less-competent shaly Ieru Formation, which exhibits disharmonic folding and variable bed thickness. Seismic-reflection data clearly show that the Plio-Pleistocene upper Era Beds are deformed to the same extent as the underlying Darai Limestone, demonstrating that most of the observed deformation has occurred during the Late Pliocene and Pleistocene.     

Evolution of non-linear α 2-dynamos and taylor's constraint

A key non-linear mechanism in a strong-field geodynamo is that a finite amplitude magnetic field drives a flow through the Lorentz force in the momentum equation and this flow feeds back on the field-generation process in the magnetic induction equation, equilibrating the field. We make use of a simpler non-linear?α?²-dynamo to investigate this mechanism in a rapidly rotating fluid spherical shell. Neglecting inertia, we use a pseudo-spectral time-stepping procedure to solve the induction equation and the momentum equation with no-slip velocity boundary conditions for a finitely conducting inner core and an insulating mantle. We present calculations for Ekman numbers (E) in the range 2.5× 10^?3 to 5.0× 10^?5, for?α?=α ₀cos?θ?sin?π?(r?r_i ) (which vanishes on both inner and outer boundaries). Solutions are steady except at lower E and higher values of?α?₀. Then they are periodic with a reversing field and a characteristic rapid increase then equally rapid decrease in magnetic energy. We have investigated the mechanism for this and shown the influence of Taylor's constraint. We comment on the application of our findings to numerical hydrodynamic dynamos.     

Logical data modelling for geographical applications

Abstract.

A formal, yet practical, GeoRelational Data Model (GRDM) is presented for the logical database design phase of the development of spatial information systems. Geographic applications are viewed in the context of information systems development. The generic needs of modelling spatial data are analyzed; it is concluded that they are not served satisfactorily by existing data models, so specifications of modelling tools for spatial application design are given. GRDM provides a set of representational constructs (relations and layers for the logical schema; virtual layers, object classes and spatial constraints for the user views) on top of well-established models. It constitutes part of a full, easily automated application design methodology. Extensive examples demonstrate the relevance, and ease-of-use of the platform-independent GRDM.     

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.     

Evolution of a small intraplate basaltic lava field: Jerrabattgulla Creek,upper Shoalhaven River catchment,southeast New South Wales

The Jerrabattgulla Creek basalts are in the upper catchment of the Shoalhaven River of southeastern New South Wales. The basalts erupted into a narrow, north-draining valley and modified the local drainage system, re-routing the paleo-Jerrabattgulla Creek, preserving a series of sub-basaltic quartzose gravels with silcretes in the paleovalley. The paleovalley indicates that a north-flowing drainage existed in this place in the Miocene. The high-relief, narrow valley has preserved a volcanic stratigraphy allowing the magmatic evolution of this small lava field to be determined. The lavas have a large compositional range from olivine nephelinite through to quartz tholeiite, which is unusual in such a small lava field. They represent three distinct magma batches, most likely from an amphibole–apatite metasomatised sub-continental lithospheric mantle, and underwent fractional crystallisation of olivine, clinopyroxene and plagioclase and assimilated upper crust. The lava field underwent temporal change from dominantly alkaline, to mixed alkaline and subalkaline, to dominantly alkaline magmatism over the course of its evolution.     

Transient events in the COBE database simultaneous with gamma ray bursts

We are investigating the COBE DMR data at instances of known -ray bursts (GRBs) when any of the six DMR horn directions was coincident with the direction of a burst. The BATSE instrument on board GRO has detected 207 bursts during the eight-month period of overlap corresponding to the current release of COBE data. The odds of a GRB occurring within the DMR field of view are near one coincidence per year. Here we report on one such serendipitous observation in 1991, GRB 911226, for which a detailed analysis is currently in progress.