A sedimentary record from a deep Quaternary valley in the southern Lillebælt area, Denmark: Eemian and Early Weichselian lithology and chronology at Mommark

A unique sequence of Late Saalian, Eemian and Early Weichselian strata is exposed in a coastal outcrop at Mommark in the western Baltic. The sedimentary facies and faunas reflect palaeoenvironmental changes from an initial freshwater lake followed by marine transgression and interglacial deposition in a palaeo-Baltic sea. The upper part of the Eemian marine record indicates regression followed by lacustrine sedimentation and deposition of Early Weichselian aeolian sediments, which are truncated by an erosional unconformity overlain by a till bed. The lower and middle parts of the sequence have previously been correlated with the European glacial-interglacial stratigraphy on the basis of pollen analysis, while the upper part has been dated for the present study using optically stimulated luminescence (OSL) of samples from the aeolian and glacial deposits. A similar complete glacial-interglacial-glacial succession has not previously been recorded from this area. The Mommark sequence of conformable strata has been subjected to lateral compression, evidenced by folding and low-angle reverse faults. Seismic records from the adjacent waters in the western Baltic reveal a system of buried Quaternary valleys in the area. It is suggested that the interglacial deposition took place in a basin within one of these valleys and that a slab constituting the Mommark sequence, originating from the margin of a valley, has been glaciotectonically displaced northwestwards to the present location.     

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.     

Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses

A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 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.     

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     

Mezen Bay—a key area for understanding Weichselian glaciations in northern Russia

Sediment successions in coastal cliffs around Mezen Bay, southeastern White Sea, record an unusually detailed history of former glaciations, interstadial marine and fluvial events from the Weichselian. A regional glaciation model for the Weichselian is based on new data from the Mezen Bay area and previously published data from adjacent areas. Following the Mikulinian (Eemian) interglacial a shelf‐centred glaciation in the Kara Sea is reflected in proglacial conditions at 100–90 ka. A local ice‐cap over the Timan ridge existed between 75 and 65 ka. Renewed glaciation in the Kara Sea spread southwestwards around 60 ka only, interrupted by a marine inundation, before it advanced to its maximum position at about 55–50 ka. After a prolonged ice‐free period, the Scandinavian ice‐sheet invaded the area from the west and terminated east of Mezen Bay about 17 ka. The previously published evidence of a large ice‐dammed lake in the central Arkhangelsk region, Lake Komi, finds no support in this study. 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.