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     

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.     

Evidence for the loss of snow-deposited MSA to the interstitial gaseous phase in central Antarctic firn

We have examined several MSA (methanesulfonic acid) records from the upper 200 m of the Antarctic ice sheet and in particular the new Dome F profile. At all the four sites studied, concentration profiles exhibit similar patterns as a function of depth. They suggest that snow metamorphism and solid phase migration are responsible for a marked release of gaseous MSA to interstitial firn air as well as probably to the free atmosphere, in particular at extremely low accumulation sites. Snow acidity can also modify MSA concentration. It is proposed that, below the upper few metres where the communication with the free atmosphere is possible, gaseous MSA may remain in the firn layers and be entrapped later in air bubbles at pore close-off, i.e. when firn is transformed into ice. Chemical measurements on the firn core do not take into account the MSA released to the gaseous phase, but this fraction is measurable in ice samples. In spite of these alterations occurring in the firn layers, relative changes of the atmospheric MSA concentration in the past are probably still there deep within the Antarctic ice sheet. However, for glacial periods, different processes have to be considered in relation to modified aerosol properties.     

Tephrochronology and tephrostratigraphy of two Pleistocene continental fossiliferous successions from central Italy

The stratigraphical context of two Middle Pleistocene fossiliferous palaeosols from Central Italy (Abruzzo and Tuscany) have been studied. Small mammals and molluscs occur in both palaeosols, which are covered by tephra layers that were analysed using an interdisciplinary approach. Application of fission‐track dating to apatites separated from the Case Picconetto tephra (Pescara, Abruzzo), yielded an age of 0.48 ± 0.04 Ma, indistinguishable from those previously determined for the Campani Quarry (Lower Valdarno, Tuscany) (0.46 ± 0.05 Ma and 0.48 ± 0.05 Ma). Geochemical and petrographic investigations indicate that these tephra originated from different volcanoes, the Alban Hills Volcanic Complex and the Vico Volcano (Latium) respectively. Small mammal and mollusc assemblages indicate different palaeoclimatic and palaeoenvironmental conditions for the Case Picconetto and Campani Quarry palaeosols. Warm and humid conditions can be inferred for the Campani Quarry site, whereas open and cold conditions can be inferred for Case Picconetto. On the basis of faunal data, fission‐track dates and attribution of tephra to specific volcanic eruptions, we suggest a correlation of these faunas with marine oxygen isotope stage 14 (Case Picconetto) and with marine oxygen isotope stage 11 (Campani Quarry), respectively. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Middle-ear stiffness of the bottlenose dolphin tursiops truncatus

Previous research on the cetacean auditory system has consisted mostly of behavioral studies on a limited number of species. Little quantitative physiologic data exists on cetacean hearing. The frequency range of hearing varies greatly across different mammalian species. Differences among species correlate with differences in the middle-ear transfer function. Middle-ear transfer functions depend on the mechanical stiffness of the middle ear and the cochlear input impedance. The purpose of this study was to measure the middle-ear stiffness for the bottlenose dolphin (Tursiops truncatus), a species specialized for underwater high-frequency hearing and echolocation. Middle-ear stiffness was measured with a force probe that applied a known displacement to the stapes and measured the restoring force. The average middle-ear stiffness in ten dolphin ears was 1.37 N//spl mu/m, which is considerably higher than that reported for most terrestrial mammals. The relationship between middle-ear stiffness and low-frequency hearing cutoff in Tursiops was shown to be comparable to that of terrestrial mammals.     

A major sulphur isotope event at c. 510 Ma: a possible anoxia–extinction–volcanism connection during the Early–Middle Cambrian transition?

A new approach to constraining seawater δ³⁴S and sulphate concentration using francolite‐bound sulphate reveals an abrupt increase in δ³⁴S to +50‰ around the Early–Middle Cambrian boundary. Such high δ³⁴S values are best explained by increased rates of pyrite burial due to ocean anoxia coupled with an increased sensitivity of the ocean sulphate reservoir to perturbations due to low sulphate concentrations of 500–700 μgL^?1. We argue that the spread of anoxic waters at this time was partly the result of greenhouse warming related to the eruption of the Kalkarindji Large Igneous Province of northern Australia and that it triggered the collapse of early metazoan reef ecosystems during the latest Early Cambrian. Mass extinctions of the last 260 Myr have all coincided with enhanced volcanic activity, while several are also associated with positive shifts in seawater δ³⁴S. Extending this correlation back in time further implicates volcanically induced climate change as a major determining factor in biosphere evolution. Terra Nova, 18, 257–263, 2006     

Chemical alteration of tephra in the depositional environment: theoretical stability modelling

The study of the chemical stability of vitreous material in aqueous media is well‐established. There has to date been little consideration of the implications of variations in the chemical durability of tephra in Quaternary tephrochronology. Chemical alteration can take the form of cationic leaching from the matrix, or complete destruction of the silica network, either of which could constrain the ability to chemically identify distal tephra. Here we apply established models of vitreous durability to the published chemical analyses of a large number of Icelandic tephras in order to predict their relative durabilities under equivalent conditions. This suggests that some important tephras have relatively poor chemical stability, and that rhyolitic tephras are, in general, more stable than basaltic. We conclude that tephras should be expected to show predictable differential chemical stability in the post‐depositional environment. Copyright © 2003 John Wiley & Sons, Ltd.     

A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. Copyright © 2003 John Wiley & Sons, Ltd.