Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Structurally controlled gold and sulfosalt mineralization: the Altenberg example, Salzburg Province, Austria

Summary ?In the south-eastern Altenbergkar–Silbereck area in the eastern Tauern window (Lungau, Salzburg) structurally controlled precious-metal (Au–Ag) mineralization is hosted in marbles of the Permo(?)-Mesozoic Silbereck Formation and in the underlying Variscan Central gneiss. During the Alpine otogeny both lithologies were affected by ductile deformation (shearing, D1; folding, D2/D3) and subsequent brittle deformation (tension gashes, D4; normal faulting, D5) related to the uplift and exhumation of the Tauern window. Mineralization is controlled by brittle D4 structures. NE–SW trending steeply dipping tension gashes of the “Tauerngoldgang” type occur within the Central gneiss. Three different marble-hosted ore types following fracture systems as well as foliation and bedding planes can be distinguished: 1) metasomatic replacement ores, 2) ores in tension gashes and 3) ores in talc-bearing structures, often containing high-grade gold and silver mineralization (native gold in association with Ag–Pb–Bi–sulfosalts). Four stages of mineralization can be distinguished which occur in all ore types: arsenopyrite–pyrite–pyrrhotite (first stage), Au–(Ag–Pb–Bi–sulfosalts) (second stage), base-metal sulfides and tetrahedrite–tennantite (third stage) and Ag-rich galena (fourth stage). Preliminary fluid inclusion data indicate temperatures of ore formation well above 300 °C (346 °C mean) for the second stage within the Central gneiss and temperatures between 310 and 230 °C for the second and third stages in the marble. Received October 12, 2001; revised version accepted September 5, 2002 Published online March 10, 2003     

Richtmyer-Meshkov Experiments on the Omega Laser

Observations of the interstellar medium reveal a dynamic realm permeated by shocks. These shocks are generated on a large range of scales by galactic rotation, supernovae, stellar winds, and other processes. Whenever a shock encounters a density interface, Richtmyer-Meshkov instabilities may develop. Perturbations along the interface grow, leading to structure formation and material mixing. An understanding of the evolution of Richtmyer-Meshkov instabilities is essential for understanding galactic structure, molecular cloud morphology, and the early stages of star formation. An ongoing experimental campaign studies Richtmyer-Meshkov mixing in a convergent, compressible, miscible plasma at the Omega laser facility. Cylindrical targets, consisting of a low density foam core and an aluminum shell covered by an epoxy ablator, are directly driven by fifty laser beams. The aluminum shell is machined to produce different perturbation spectra. Surface types include unperturbed (smooth), single-mode sinusoids, multi-mode (rough), and multi-mode with particular modes accentuated (specified-rough). Experimental results are compared to theory and numerical simulations.     

Biological and granulometric controls on sedimentary organic matter of an intertidal mudflat

A variety of measures of organic matter concentration and quality were made on samples collected from the top few mm of intertidal mudflat sediment over the course of a year, in order to assess the relative importance of biological and sedimentological influences on sedimentary organic matter. Winter and summer were times of relatively fine-grained sediment accumulation, caused by biological deposition or stabilization processes and resulting in higher organic matter concentrations. Stable carbon isotope and Br:C ratios indicated a planktonic source of bulk organic matter. Ratios of organic carbon to specific surface area of the sediments were consistent with an organic monolayer coverage of sediment grains. Correction for changing grain size during the year showed no change in the organic concentration per unit surface area, in spite of organic matter inputs by in situ primary production, buildup of heterotroph biomass and mucus coatings, and biodeposition of organic-rich seston. There were also no indications of changes in bulk organic quality, measured as hydrolyzable carbohydrates and amino acids, in response to these biological processes. It is concluded that biological processes on a seasonal time scale affect the bulk organic matter of these sediments via a modulation of grain size rather than creation or decay of organic matter.     

Effect of metamorphism and partial melting of host rocks on zircons

Abstract Zircons have been studied from different layers of migmatites (from Arvika, western Sweden and Nelaug, southern Norway) and from a paragneiss (from Arvika) associated with one of the migmatites. The main purpose of the investigation is to establish whether or not information about zircons can help in the elucidation of the parentage and rock-forming processes of migmatites.
The elongation ratio of zircons from all layers is small and characteristic of sedimentary zircons. Further, the absence of characteristic colours and the growth trends of the zircons (indicated by the reduced major axes) observed in the various samples both support a sedimentary parentage for these rocks. The zircons of all layers exhibit secondary growth (overgrowth, outgrowth and multiple growth) due to metamorphism. Compared with the zircons from the paragneiss, those of the migmatite layers are more clouded and less rounded, some of them becoming opaque or even skeletal; this is especially true of the zircons from the leucosomes. These observations indicate an alteration of the original sedimentary zircons in the migmatite, especially in the leucosomes, in response to the migmatization process, previously interpreted as partial melting.     

The nature of Europa's dark non-ice surface material: Spatially-resolved high spectral resolution spectroscopy from the Keck telescope

We present new 1.45-1.75 μm spectra of Europa's dark non-ice material with a spectral resolution (λ/δλ) of 1200, obtained by using adaptive optics on the Keck telescope to spatially separate the spectrum of the non-ice material from that of the surrounding ice-rich regions. Despite the great increase in spectral resolution over the previous best spectra of the non-ice material, taken with Galileo's near-infrared mapping spectrometer (NIMS) with λ/δλ=66, no new fine-scale spectral structure is revealed. The smoothness of the spectra is inconsistent with available laboratory spectra of crystalline hydrated salts at Europa temperatures, but is more consistent with various combinations of non-crystalline hydrated salts and/or hydrated sulfuric acid, as have been matched to the lower-resolution NIMS spectra.     

The NAIAD experiment for WIMP searches at Boulby mine and recent results

The NAIAD experiment (NaI Advanced Detector) for weakly interacting massive particle (WIMP) dark matter searches at Boulby mine (UK) is described. The detector consists of an array of encapsulated and unencapsulated NaI(Tl) crystals with high light yield. Six crystals are collecting data at present. Data accumulated by four of them (10.6 kg × year exposure) have been used to set upper limits on the WIMP–nucleon spin-independent and WIMP–proton spin-dependent cross-sections. Pulse shape analysis has been applied to discriminate between nuclear recoils, as may be caused by WIMP interactions, and electron recoils due to gamma background. Various calibrations of crystals are presented.     

A statistical evaluation of anomalous redshift claims

Claims that ordinary spiral galaxies and some classes of QSO show periodicity in their redshift distributions are investigated using recent high-precision data and rigorous statistical procedures. The claims are broadly upheld. The periodicites are strong and easily seen by eye in the datasets. Observational, reduction or statistical artefacts do not seem capable of accounting for them. This revised version was published online in August 2006 with corrections to the Cover Date.     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.