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     

The evolution of binary fractions in globular clusters

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code ( startrack ) and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing three- and four-body interactions ( fewbody ). We find that the combination of stellar evolution and dynamical interactions (binary–single and binary–binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster such as 47 Tuc, assuming an initial binary fraction of 100 per cent, is only ∼ 5–10 per cent. We show that this is in good agreement with recent Hubble Space Telescope observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modelling of globular clusters, suggesting that 'realistic models' should incorporate much larger initial binary fractions than has usually been the case in the past.     

The depolarization properties of powerful extragalactic radio sources as a function of cosmic epoch

We use the observed polarization properties of a sample of 26 powerful radio galaxies and radio-loud quasars to constrain the conditions in the Faraday screens local to the sources. We adopt the cosmological redshift, low-frequency radio luminosity and physical size of the large-scale radio structures as our 'fundamental' parameters. We find no correlation of the radio spectral index with any of the fundamental parameters. The observed rotation measure is also independent of these parameters, suggesting that most of the Faraday rotation occurs in the Galactic foreground. The difference between the rotation measures of the two lobes of an individual source, as well as the dispersion of the rotation measure, shows significant correlations with the source redshift, but not with the radio luminosity or source size. This is evidence that the small-scale structure observed in the rotation measure is caused by a Faraday screen local to the sources. The observed asymmetries between the lobes of our sources show no significant trends with each other or other source properties. Finally, we show that the commonly used model for the depolarization of synchrotron radio emission by foreground Faraday screens is inconsistent with our observations. We apply alternative models to our data and show that they require a strong increase of the dispersion of the rotation measure inside the Faraday screens with cosmological redshift. Correcting our observations with these models for redshift effects, we find a strong correlation of the depolarization measure with redshift and a significantly weaker correlation with radio luminosity. We do not find any (anti-)correlation of depolarization measure with source size. All our results are consistent with a decrease in the order of the magnetic field structure of the Faraday screen local to the sources for increasing cosmological redshift.     

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.     

Pc1-2 Discrete regular daytime pulsation bursts at high latitudes

A family of related Pc1-2 (0.2–10 s) discrete daytime geomagnetic pulsations is presented using pulsation data obtained at Davis, Antarctica, a typical polar-cap station. The morphological properties of IPRP and Pclb pulsation regimes, which maximize in amplitude and frequency of occurrence under the projection of the polar cusp, are examined. Furthermore, two other variations of discrete pulsation bursts yet to be named are also presented, viz IPFP (Intervals of Pulsations with Falling Period) and IPAP (Intervals of Pulsations with Alternating Period) which are observed on rare occasions. It is also suggested that the Pc1b (0.2–5 s) should be extended to incorporate Pc2b (5–10 s) which from the results in this paper are physically the same phenomenon and could be collectively classified as IPCP (Intervals of Pulsations with Constant Period).     

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.