Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Rare earth, major and trace elements in Jurassic manganese shales of the Northern Calcareous Alps: hydrothermal versus hydrogenous origin of stratiform manganese deposits

Summary ?The rare earth, major and trace element geochemistry of Jurassic deep marine manganese shales allow insight into their environment of deposition. We present data of 24 samples from the Northern Calcareous Alps (Eastern Alps), collected from Late Jurassic strata of the Tennengebirge and from Early Jurassic strata of the Karwendel Mts. Whereas major and trace element geochemistry is controlled by the detrital input, rare earth element data and the relationship between total organic carbon and (pyritic) sulfur can be used to constrain the origin of the manganese accumulation. The data of this study are compatible with a hydrogenous precipitation of manganese by strongly varying redox-conditions at the sediment–water interface of a slope basin. Received December 21, 2001; revised version accepted February 18, 2002     

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.     

Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

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