Hässeldala – a key site for Last Termination climate events in northern Europe

The Last Termination (19 000–11 000 a BP) with its rapid and distinct climate shifts provides a perfect laboratory to study the nature and regional impact of climate variability. The sedimentary succession from the ancient lake at Hässeldala Port in southern Sweden with its distinct Lateglacial/early Holocene stratigraphy (>14.1–9.5 cal. ka BP) is one of the few chronologically well‐constrained, multi‐proxy sites in Europe that capture a variety of local and regional climatic and environmental signals. Here we present Hässeldala's multi‐proxy records (lithology, geochemistry, pollen, diatoms, chironomids, biomarkers, hydrogen isotopes) in a refined age model and place the observed changes in lake status, catchment vegetation, summer temperatures and hydroclimate in a wider regional context. Reconstructed mean July temperatures increased between c. 14.1 and c. 13.1 cal. ka BP and subsequently declined. This latter cooling coincided with drier hydroclimatic conditions that were probably associated with a freshening of the Nordic Seas and started a few hundred years before the onset of Greenland Stadial 1 (c. 12.9 cal. ka BP). Our proxies suggest a further shift towards colder and drier conditions as late as c. 12.7 cal. ka BP, which was followed by the establishment of a stadial climate regime (c. 12.5–11.8 cal. ka BP). The onset of warmer and wetter conditions preceded the Holocene warming over Greenland by c. 200 years. Hässeldala's proxies thus highlight the complexity of environmental and hydrological responses across abrupt climate transitions in northern Europe.     

Flocculation settling characteristics of mud: sand mixtures

When natural muds become mixed with sandy sediments in estuaries, it has a direct effect on the flocculation process and resultant sediment transport regime. Much research has been completed on the erosion and consolidation of mud/sand mixtures, but very little is known quantitatively about how mixed sediments interact whilst in suspension, particularly in terms of flocculation. This paper presents the settling velocity findings from a recent laboratory study which examined the flocculation dynamics for three different mud/sand mixtures at different concentrations (0.2–5 g.l^?1) and turbulent shear stresses (0.06–0.9 Pa) in a mini-annular flume. The low intrusive video-based Laboratory Spectral Flocculation Characteristics instrument was used to determine floc/aggregate properties (e.g., size, settling velocity, density and mass) for each population. Settling data was assessed in terms of macrofloc (>160 μm) and microfloc (<160 μm) settling parameters: Ws_macro and Ws_micro, respectively. For pure muds, the macroflocs are regarded as the most dominant contributors to the total depositional flux. The parameterised settling data indicates that by adding more sand to a mud/sand mixture, the fall velocity of the macrofloc fraction slows and the settling velocity of microflocs quickens. Generally, a mainly sandy suspension comprising 25% mud and 75% sand (25M:75S), will produce resultant Ws_macro which are slower than Ws_micro. The quickest Ws_micro appears to consistently occur at a higher level of turbulent shear stress (τ?～?0.6 Pa) than both the macrofloc and microfloc fractions from suspensions of pure natural muds. Flocculation within a more cohesively dominant muddy-sand suspension (i.e., 75M:25S) produced macroflocs which fell at similar speeds (±10%) to pure mud suspensions at both low (200 mg l^?1) and intermediate (1 g?l^?1) concentrations at all shear stress increments. Also, low sand content suspensions produced Ws_macro values that were faster than the Ws_micro rates. In summary, the experimental results of the macrofloc and microfloc settling velocities have demonstrated that flocculation is an extremely important factor with regards to the depositional behaviour of mud/sand mixtures, and these factors must be considered when modelling mixed sediment transport in the estuarine or marine environment.     

Multi-Collector SIMS Determination of Trace Lanthanides in Zircon

A method is presented for the determination of rare earth elements (REE) in zircon (ZrSiO₄) using a multi-collector equipped ion-microprobe operating at moderately high mass resolution (M/ΔM ˜ 3900). The low abundance light REE (La-Eu) were measured simultaneously in six ion counting electron multipliers, reducing analytical time to less than half of that in a monocollection routine. Effective mass filtering, moderate energy filtering and careful set up of detector parameters in order to minimise background counts, yielded a highly coherent set of analyses from the 91500 zircon reference sample.     

Detrital zircon ages from southern Norway – implications for the Proterozoic evolution of the southwestern Baltic Shield

An ion-microprobe (SIMS) U-Pb zircon dating study on four samples of Precambrian metasediments from the high-grade Bamble Sector, southern Norway, gives the first information on the timing of discrete crust-forming events in the SW part of the Baltic Shield. Recent Nd and Pb studies have indicated that the sources of the clastic metasediments in this area have crustal histories extending back to 1.7 to 2.1 Ga, although there is no record of rocks older than 1.6 Ga in southern Norway. The analysed metasediments are from a sequence of intercalated, centimetre to 10-metre wide units of quartzites, semi-metapelites, metapelites and mafic granulites. The zircons can be grouped in two morphological populations: (1) long prismatic; (2) rounded, often flattened. The BSE images reveal that both populations consist of oscillatory zoned, rounded and corroded cores (detrital grains of magmatic origin), surrounded by homogeneous rims (metamorphic overgrowths). The detrital zircons have ²⁰⁷Pb/²⁰⁶Pb ages between 1367 and 1939 Ma, with frequency maxima in the range 1.85 to 1.70 Ga and 1.60 to 1.50 Ga. There is no correlation between crystal habit and age of the zircon. One resorbed, inner zircon core in a detrital grain is strongly discordant and gives a composite inner core-magmatic outer core ²⁰⁷Pb/²⁰⁶Pb age of 2383 Ma. Two discrete, unzoned zircons have ²⁰⁷Pb/²⁰⁶Pb ages of 1122 and 1133 Ma, representing zircon growth during the Sveconorwegian high-grade metamorphism. Also the μm wide overgrowths, embayments in the detrital cores and apparent “inner cores” which represent secondary metamorphic zircon growth in deep embayments in detrital grains, are of Sveconorwegian age. The composite-detrital-metamorphic zircon analyses give generally discordant ²⁰⁶Pb/²³⁸U versus ²⁰⁷Pb/²³⁵U ratios and maximum ²⁰⁷Pb/²⁰⁶Pb ages of 1438 Ma. These data demonstrate the existence of a protocrust of 1.7 to 2.0 Ga in the southwestern part of the Baltic Shield, implying a break in the overall westward younging trend of the Precambrian crust, inferred from the southeastern part of the Baltic Shield. Received: 8 April 1997 / Accepted: 14 July 1997     

4D analogue modelling of transtensional pull-apart basins

Scaled sandbox models were used to investigate the 4D evolution of pull-apart basins formed above underlapping releasing stepovers in both pure strike-slip and transtensional basement fault systems. Serial sectioning and 3D volume reconstruction permitted analysis of the full 3D fault geometries. Results show that very different pull-apart basins are developed in transtension compared to pure strike-slip. Both types of models produced elongate, sigmoidal to rhomboidal pull-apart systems, but the transtensional pull-apart basins were significantly wider and uniquely developed a basin margin of en-echelon oblique-extensional faults. Dual, opposing depocentres formed in the transtensional model whereas a single, central depocentre formed in pure strike-slip. In transtension, a distinct narrow graben system formed above the principal displacement zones (PDZs). Cross-basin fault systems that linked the offset PDZs formed earlier in the transtensional models.Sequential model runs to higher PDZ displacements allowed the progressive evolution of the fault systems to be evaluated. In cross-section, transtensional pull-aparts initiated as asymmetric grabens bounded by planar oblique-extensional faults. With increasing displacement on the PDZs, basin subsidence caused these faults to become concave-upwards and lower in dip angle due to fault block collapse towards the interior of the basin. In addition, strain partitioning caused fault slip to become either predominantly extensional or strike-slip. The models compare closely with the geometries of natural pull-apart basins including the southern Dead Sea fault system and the Vienna Basin, Austria.     

An integrated zircon geochronological and geochemical investigation into the Miocene plutonic evolution of the Cyclades, Aegean Sea, Greece: part 2—geochemistry

Zircons from 14 compositionally variable granitic rocks were examined in detail using CL image-guided micro-analysis to unravel the complex magmatic history above the southward retreating Hellenic subduction zone system in the Aegean Sea. Previously published U–Pb ages document an episodic crystallisation history from 17 to 11?Ma, with peraluminous (S-type) granitic rocks systematically older than closely associated metaluminous (I-type) granitic rocks. Zircon O- and Hf isotopic data, combined with trace element compositions, are highly variable within and between individual samples, indicative of open-system behaviour involving mantle-derived melts and evolved supracrustal sources. Pronounced compositional and thermal fluctuations highlight the role of magma mixing and mingling, in accord with field observations, and incremental emplacement of distinct melt batches coupled with variable degrees of crustal assimilation. In the course of partial fusion, more fertile supracrustal sources dominated in the earlier stages of Aegean Miocene magmatism, consistent with systematically older crystallisation ages of peraluminous granitic rocks. Differences between zircon saturation and crystallisation temperatures (deduced from zircon Ti concentrations), along with multimodal crystallisation age spectra for individual plutons, highlight the complex and highly variable physico-compositional and thermal evolution of silicic magma systems. The transfer of heat and juvenile melts from the mantle varied probably in response to episodic rollback of the subducting lithospheric slab, as suggested by punctuated crystallisation age spectra within and among individual granitic plutons.     

The fate of zircon during UHT–UHP metamorphism: isotopic (U/Pb, δ18O,Hf) and trace element constraints

Garnet–clinopyroxene ultra‐high‐pressure (UHP) rocks from the northern Bohemian Massif contain zircon with micro‐diamond inclusions. Trace element concentrations, oxygen and hafnium isotopic composition and U–Pb age of distinct textural domains in zircon characterize their growth conditions and temporal evolution. Diamond‐bearing zircon mantle domains with relicts of oscillatory zoning have uniform Th/U ratios (~0.1–0.2), high‐Ti contents (110–190 ppm, corresponding to temperatures of at least 1100 °C), and some (two of 17 mantle analyses) preserve steep heavy rare earth element (HREE) patterns with Yb_N/Gd_N = 10–11, with a weak negative Eu anomaly. These signatures are consistent with crystallization from a melt under UHP/ultra‐high‐temperature (UHT) conditions. Some of the bright‐cathodoluminscence (CL) rims preserve Th/U and Ti values characteristic of the zircon mantles, but others show elevated Th/U ratios of ~0.3–0.4 and lower Ti contents (20–40 ppm; only 13 ppm in a rare low‐CL outer rim). As they feature flat HREE patterns and negative Eu anomalies and commonly make embayments and truncate the mantle zoning, we suggest that they have formed through recrystallization in the solid state during exhumation of the rock, when both garnet and plagioclase were stable. The three zircon domains, that is, cores, mantles and rims, yield U–Pb concordia ages of 340.9 ± 1.5, 340.3 ± 1.5 and 341.2 ± 3.4 Ma respectively. When linked to the previously reconstructed P–T path of the rock, the error limits of the zircon mantle and rim ages constrain the exhumation of the rocks from depth of ~140 km (UHP) to ~80 km (HP) to a minimum rate of 1.5 cm yr^?1. The zircon cores are heterogeneous in terms of Th/U ratio (below 0.1 but also above 0.2) and REE characteristics, and their εHf values scatter between ?15.7 and +4.8 with similar values for individual domains within a single zircon grain suggesting a very localized control on hafnium isotope composition on a grain scale. The non‐equilibrated εHf values as well as a large range of the Hf‐depleted mantle model ages possibly reflect the presence of a heterogeneous population of old zircon. Consequently, the uniform and young ²³⁸U/²⁰⁶Pb ages may represent (near‐)complete resetting of the U–Pb geochronometer during the UHP–UHT event at c. 340 Ma through dissolution–reprecipitation process. In contrast to Hf, the oxygen isotope composition of zircon is homogeneous, ranging between 7.8‰ and 9.6‰ VSMOW, reflecting a source containing upper crustal material and homogenization at UHP–UHT conditions. Our study documents that continental crust was subducted to mantle depths at c. 340 Ma during the Variscan orogeny and was subsequently very rapidly exhumed, implying that the sequence of events was faster than can be resolved by the secondary ion mass spectrometry technique.     

Behaviour of geochronometers and timing of metamorphic reactions during deformation at lower crustal conditions: phase equilibrium modelling and U–Pb dating of zircon,monazite, rutile and titanite from the Kalak Nappe Complex,northern Norway

This study investigates the behaviour of the geochronometers zircon, monazite, rutile and titanite in polyphase lower crustal rocks of the Kalak Nappe Complex, northern Norway. A pressure–temperature–time–deformation path is constructed by combining microstructural observations with P–T conditions derived from phase equilibrium modelling and U–Pb dating. The following tectonometamorphic evolution is deduced: A subvertical S1 fabric formed at ~730–775 °C and ~6.3–9.8 kbar, above the wet solidus in the sillimanite and kyanite stability fields. The event is dated at 702 ± 5 Ma by high‐U zircon in a leucosome. Monazite grains that grew in the S1 fabric show surprisingly little variation in chemical composition compared to a large spread in (concordant) U–Pb dates from c. 800 to 600 Ma. This age spread could either represent protracted growth of monazite during high‐grade metamorphism, or represent partially reset ages due to high‐T diffusion. Both cases imply that elevated temperatures of >600 °C persisted for over c. 200 Ma, indicating relatively static conditions at lower crustal levels for most of the Neoproterozoic. The S1 fabric was overprinted by a subhorizontal S2 fabric, which formed at ~600–660 °C and ~10–12 kbar. Rutile that originally grew during the S1‐forming event lost its Zr‐in‐rutile and U–Pb signatures during the S2‐forming event. It records Zr‐in‐rutile temperatures of 550–660 °C and Caledonian ages of 440–420 Ma. Titanite grew at the expense of rutile at slightly lower temperatures of ~550 °C during ongoing S2 deformation; U–Pb ages of c. 440–430 Ma date its crystallization, giving a minimum estimate for the age of Caledonian metamorphism and the duration of Caledonian shearing. This study shows that (i) monazite can have a large spread in U–Pb dates despite a homogeneous composition; (ii) rutile may lose its Zr‐in‐rutile and U–Pb signature during an amphibolite facies overprint; and (iii) titanite may record crystallization ages during retrograde shearing. Therefore, in order to correctly interpret U–Pb ages from different geochronometers in a polyphase deformation and reaction history, they are ideally combined with microstructural observations and phase equilibrium modelling to derive a complete P–T–t–d path.     

The settling dynamics of flocculating mud and sand mixtures: part 2—numerical modelling

Estuarine and coastal sediment transport is characterised by the transport of both sand-sized particles (of diameter greater than 63?μm) and muddy fine-grained sediments (silt, diameter less than 63?μm; clay, diameter less than 2?μm). These fractions are traditionally considered as non-cohesive and cohesive, respectively, because of the negligible physico-chemical attraction that occurs between sand grains. However, the flocculation of sediment particles is not only caused by physico-chemical attraction. Cohesivity of sediment is also caused by biology, in particular the sticky extra-cellular polymeric substances secreted by diatoms, and the effect of biology in binding sediment particles can be much larger than that of physico-chemical attraction. As demonstrated by Manning (2008) and further expanded in part 1 of this paper (Manning et al., submitted), the greater binding effect of biology allows sand particles to flocculate with mud. In many estuaries, both the sand and fine sediment fractions are transported in significant quantities. Many of the more common sediment transport modelling suites now have the capability to combine mud and sand transport. However, in all of these modelling approaches, the modelling of mixed sediment transport has still essentially separated the modelling of sand and mud fractions assuming that these different fractions do not interact except at the bed. However, the use of in situ video techniques has greatly enhanced the accuracy and reliability of settling velocity measurements and has led to a re-appraisal of this widely held assumption. Measurements of settling velocity in mixed sands presented by Manning et al. (2009) have shown strong evidence for the flocculation of mixed sediments, whilst the greater understanding of the role of biology in flocculation has identified mechanisms by which this mud-sand flocculation can occur. In the first part of this paper (Manning et al., submitted), the development of an empirical flocculation model is described which represents the interaction between sand and mud particles in the flocculation process. Measurements of the settling velocity of varying mud-sand mixtures are described, and empirical algorithms governing the variation of settling velocity with turbulence, suspended sediment concentration and mud-sand content are derived. The second part of this paper continues the theme of examination of the effects of mud-sand interaction on flocculation. A 1DV mixed transport model is developed and used to reproduce the vertical transport of mixed sediment fractions. The 1DV model is used to reproduce the measured settling velocities in the laboratory experiments described in the part 1 paper and also to reproduce measurements of concentration of mixed sediments in the Outer Thames. In both modelling exercises, the model is run using the algorithms developed in part 1 and repeated using an assumption of no interaction between mud and sand in the flocculation process. The results of the modelling show a significant improvement in the ability of the 1DV to reproduce the observed sediment behaviour when the empirical equations are used. This represents further strong evidence of the interaction between sand and mud in the flocculation process.