Discriminating between pore‐filling load and bed‐structure load: a new porosity‐based method,exemplified for the river Rhine

Sediments contained in the river bed do not necessarily contribute to morphological change. The finest part of the sediment mixture often fills the pores between the larger grains and can be removed without causing a drop in bed level. The discrimination between pore‐filling load and bed‐structure load, therefore, is of practical importance for morphological predictions. In this study, a new method is proposed to estimate the cut‐off grain size that forms the boundary between pore‐filling load and bed‐structure load. The method evaluates the pore structure of the river bed geometrically. Only detailed grain‐size distributions of the river bed are required as input to the method. A preliminary validation shows that the calculated porosity and cut‐off size values agree well with experimental data. Application of the new cut‐off size method to the river Rhine demonstrates that the estimated cut‐off size decreases in a downstream direction from about 2 to 0·05 mm, covariant with the downstream fining of bed sediments. Grain size fractions that are pore‐filling load in the upstream part of the river thus gradually become bed‐structure load in the downstream part. The estimated (mass) percentage of pore‐filling load in the river bed ranges from 0% in areas with a unimodal river bed, to about 22% in reaches with a bimodal sand‐gravel bed. The estimated bed porosity varies between 0·15 and 0·35, which is considerably less than the often‐used standard value of 0·40. The predicted cut‐off size between pore‐filling load and bed‐structure load (D_c,p) is fundamentally different from the cut‐off size between wash‐load and bed‐material load (D_c,w), irrespective of the method used to determine D_c,p or D_c,w. D_c,w values are in the order of 10^?1 mm and mainly dependent on the flow characteristics, whereas D_c,p values are generally much larger (about 10⁰ mm in gravel‐bed rivers) and dependent on the bed composition. Knowledge of D_c,w is important for the prediction of the total sediment transport in a river (including suspended fines that do not interact with the bed), whereas knowledge of D_c,p helps to improve morphological predictions, especially if spatial variations in D_c,p are taken into account. An alternative to using a spatially variable value of D_c,p in morphological models is to use a spatially variable bed porosity, which can also be predicted with the new method. In addition to the morphological benefits, the new method also has sedimentological applications. The possibility to determine quickly whether a sediment mixture is clast‐supported or matrix‐supported may help to better understand downstream fining trends, sediment entrainment thresholds and variations in hydraulic conductivity.     

Trace element distribution among rock-forming minerals from metamorphosed to partially molten basic igneous rocks in a contact aureole (Fuerteventura, Canaries)

Low pressure partial melting of basanitic and ankaramitic dykes gave rise to unusual, zebra-like migmatites, in the contact aureole of a layered pyroxenite–gabbro intrusion, in the root zone of an ocean island (Basal Complex, Fuerteventura, Canary Islands). These migmatites are characterised by a dense network of closely spaced, millimetre-wide leucocratic segregations. Their mineralogy consists of plagioclase (An_32–36), diopside, biotite, oxides (magnetite, ilmenite), +/− amphibole, dominated by plagioclase in the leucosome and diopside in the melanosome. The melanosome is almost completely recrystallised, with the preservation of large, relict igneous diopside phenocrysts in dyke centres. Comparison of whole-rock and mineral major- and trace-element data allowed us to assess the redistribution of elements between different mineral phases and generations during contact metamorphism and partial melting.

Dykes within and outside the thermal aureole behaved like closed chemical systems. Nevertheless, Zr, Hf, Y and REEs were internally redistributed, as deduced by comparing the trace element contents of the various diopside generations. Neocrystallised diopside – in the melanosome, leucosome and as epitaxial phenocryst rims – from the migmatite zone, are all enriched in Zr, Hf, Y and REEs compared to relict phenocrysts. This has been assigned to the liberation of trace elements on the breakdown of enriched primary minerals, kaersutite and sphene, on entering the thermal aureole. Major and trace element compositions of minerals in migmatite melanosomes and leucosomes are almost identical, pointing to a syn- or post-solidus reequilibration on the cooling of the migmatite terrain i.e. mineral–melt equilibria were reset to mineral–mineral equilibria.     

Partial crustal melting beneath the Betic Cordillera (SE Spain): The case study of Mar Menor volcanic suite

The Neogene Volcanic Province (NVP) within the Betic Cordillera (SE Spain) consists of three main metapelitic enclave suites (from SW to NE: El Hoyazo, Mazarrón and Mar Menor). Since the NVP represents a singular place in the world where crustal enclaves were immediately quenched after melting, their microstructures provide a “photograph” of the conditions at depth just after the moment of the melting.

The thermobarometric information provided by the different microstructural assemblages has been integrated with the geophysical and geodynamical published data into a model of the petrologic evolution of the Mar Menor enclaves. They were equilibrated at 2–3 kbar, 850–900 °C, and followed a sequence of heating melt producing reactions. A local cooling event evidenced by minor melt crystallization preceded the eruption.

The lower crustal studies presented in this work contribute to the knowledge of: (i) the partial melting event beneath the Mar Menor volcanic suite through a petrologic detailed study of the enclaves; (ii) how the microstructures of fast cooled anatectic rocks play an important role in tracing the magma evolution in a chamber up to the eruption, and how they can be used as pseudothermobarometers; (iii) the past and current evolution of the Alborán Domain (Betic Cordillera) and Mediterranean Sea, and how the base of a metapelitic crust has melted within an active geodynamic setting.     

Cenozoic Eucla Basin and associated palaeovalleys,southern Australia — Climatic and tectonic influences on landscape evolution,sedimentation and heavy mineral accumulation

The Eucla Basin including the vast Nullarbor Plain lies on the margins of the Yilgarn, Musgrave and Gawler cratons in southern Australia and owes its distinctive landscape to a unique set of interactions between eustatic, climatic and tectonic processes over the last ~ 50 Ma. Understanding of the history of the basin and the palaeovalleys that drained from the surrounding cratons are important because they contain major mineral deposits, and the sediments derived from them contain remobilised gold, uranium, and heavy minerals. In particular, a remarkably preserved palaeoshoreline sequence along the north-eastern margin of the Eucla Basin is highly prospective for heavy mineral placer deposits. The record of marine, marginal marine, estuarine, fluvial and lacustrine environments, as constrained mainly by an extensive borehole dataset, reflects major depositional events during the Palaeocene–Early Eocene, Middle–Late Eocene, Oligocene–Early Miocene, Middle Miocene–Early Pliocene and Pliocene–Quaternary. These events reflect the key role of eustatic sea-level variation which, during highstands, inundated the craton margins, flooding palaeovalleys to up to 400 km inboard of the present coastline. However, a systematic eastward migration of the depocentre across the Eucla Basin during the Neogene, together with apparent flow reversals in a number of palaeovalley systems draining the Gawler Craton, suggest that the Eucla Basin has also been subject to differential vertical movements, expressed as a west-side up, east-side down tilting of ~ 100–200 m. This differential movement forms part of a broader north-down–southwest-up dynamic topographic tilting of the Australian continent associated with relatively fast (6–7 cm/yr) northward plate motion since fast spreading commenced in the Southern Ocean at ~ 43 Ma. We suggest that the evolving dynamic topography field has played a key role in facilitating development of placer deposits, largely through multistage, eastward reworking of near-shore sequences during highstand transgressive cycles on a progressively tilting platform under the influence of persistent westerly weather systems.     

The Anarak, Jandaq and Posht-e-Badam metamorphic complexes in central Iran: New geological data, relationships and tectonic implications

The Anarak, Jandaq and Posht-e-Badam metamorphic complexes occupy the NW part of the Central-East Iranian Microcontinent and are juxtaposed with the Great Kavir block and Sanandaj-Sirjan zone. Our recent findings redefine the origin of these complexes, so far attributed to the Precambrian–Early Paleozoic orogenic episodes, and now directly related to the tectonic evolution of the Paleo-Tethys Ocean. This tectonic evolution was initiated by Late Ordovician–Early Devonian rifting events and terminated in the Triassic by the Eocimmerian collision event due to the docking of the Cimmerian blocks with the Asiatic Turan block.

The “Variscan accretionary complex” is a new name we proposed for the most widely distributed metamorphic rocks connected to the Anarak and Jandaq complexes. This accretionary complex exposed from SW of Jandaq to the Anarak and Kabudan areas is a thick and fine grain siliciclastic sequence accompanied by marginal-sea ophiolitic remnants, including gabbro-basalts with a supra-subduction-geochemical signature. New ⁴⁰Ar/³⁹Ar ages are obtained as 333–320 Ma for the metamorphism of this sequence under greenschist to amphibolite facies. Moreover, the limy intercalations in the volcano-sedimentary part of this complex in Godar-e-Siah yielded Upper Devonian–Tournaisian conodonts. The northeastern part of this complex in the Jandaq area was intruded by 215 ± 15 Ma arc to collisional granite and pegmatites dated by ID-TIMS and its metamorphic rocks are characterized by some ⁴⁰Ar/³⁹Ar radiometric ages of 163–156 Ma.

The “Variscan” accretionary complex was northwardly accreted to the Airekan granitic terrane dated at 549 ± 15 Ma. Later, from the Late Carboniferous to Triassic, huge amounts of oceanic material were accreted to its southern side and penetrated by several seamounts such as the Anarak and Kabudan. This new period of accretion is supported by the 280–230 Ma ⁴⁰Ar/³⁹Ar ages for the Anarak mild high-pressure metamorphic rocks and a 262 Ma U–Pb age for the trondhjemite–rhyolite association of that area. The Triassic Bayazeh flysch filled the foreland basin during the final closure of the Paleo-Tethys Ocean and was partly deposited and/or thrusted onto the Cimmerian Yazd block.

The Paleo-Tethys magmatic arc products have been well-preserved in the Late Devonian–Carboniferous Godar-e-Siah intra-arc deposits and the Triassic Nakhlak fore-arc succession. On the passive margin of the Cimmerian block, in the Yazd region, the nearly continuous Upper Paleozoic platform-type deposition was totally interrupted during the Middle to Late Triassic. Local erosion, down to Lower Paleozoic levels, may be related to flexural bulge erosion. The platform was finally unconformably covered by Liassic continental molassic deposits of the Shemshak.

One of the extensional periods related to Neo-Tethyan back-arc rifting in Late Cretaceous time finally separated parts of the Eocimmerian collisional domain from the Eurasian Turan domain. The opening and closing of this new ocean, characterized by the Nain and Sabzevar ophiolitic mélanges, finally transported the Anarak–Jandaq composite terrane to Central Iran, accompanied by large scale rotation of the Central-East Iranian Microcontinent (CEIM). Due to many similarities between the Posht-e-Badam metamorphic complex and the Anarak–Jandaq composite terrane, the former could be part of the latter, if it was transported further south during Tertiary time.     

Petrography and textural development of inorganic and biogenic lithotypes in a relict barite tufa deposit at Flybye Springs, NT, Canada

A relict mound of Holocene barite (BaSO₄) tufa underlies the Flybye Springs, a small, barium‐rich, cold sulphur spring system in the Northwest Territories of Canada. The tufa is composed of relatively pure barite with ≤0·34 wt% Ca²⁺ and ≤0·77 wt% Sr²⁺. The mound is made up of coated bubble, raft, undulatory sheet, stromatolitic, coated grain and detrital conglomerate barite tufa. Although previously unreported in barite, these lithotypes are akin to facies found in many carbonate spring deposits. Raft and ooid‐coated grain tufa was formed via ‘inorganic’ barite precipitation in spring water ponds and tributaries where rapid oxidation of sulphide to sulphate established barite supersaturation. Undulatory sheet tufa may have formed by the reaction of dissolved barium with sulphate derived from the oxidation of extracellular polysaccharide‐rich colloidal sulphur films floating in oxygenated, barite‐saturated spring water ponds. Coated bubble, oncoid‐coated grain and stromatolitic tufa with filamentous microfossils was formed in close association with sulphur‐tolerant microbes inhabiting dysoxic and oxygenated spring water tributaries and ponds. Adsorption of dissolved barium to microbial extracellular polysaccharide probably facilitated the development of these ‘biogenic’ lithotypes. Detrital conglomerate tufa was formed by barite cementation of microdetrital tufa, allochthonous lithoclasts and organic detritus, including caribou hair. Biogenic textures, organic artefacts and microfossils in the Flybye barite tufa have survived diagenetic aggradational recrystallization and precipitation of secondary cements, indicating the potential for palaeoecological information to be preserved in barite in the geological record. Similarities between the Flybye barite tufa and carbonate spring deposits demonstrate that analogous textures can develop in chemical sedimentary systems with distinct mineralogy, biology and physiochemistry.     

Late Quaternary contourites and glaciomarine sedimentation in the Fram Strait

Two sites in the eastern Fram Strait, the Vestnesa Ridge and the Yermak Plateau, have been surveyed and sampled providing a depositional record over the last glacial‐interglacial cycle. The Fram Strait is the only deep‐water connection from the Arctic Ocean to the North Atlantic and contains a marine sediment record of both high latitude thermohaline flow and ice sheet interaction. On the Vestnesa Ridge, the western Svalbard margin, a sediment drift was identified in 1226 m of water. Gravity and multicores from the crest of the drift recovered turbidites and contourites. ¹⁴C dating indicates an age range of 8287 to 26 900 years BP (Early Holocene to Late Weichselian). The Yermak Plateau is characterized by slope sediments in 961 m of water. Gravity and multicores recovered contourites and hemipelagites. ¹⁴C ages were between 8615 and 46 437 years BP (Early Holocene to mid‐Weichselian). Downcore dinoflagellate cyst analyses from both sites provide a record of changing surface water conditions since the mid‐Weichselian, suggesting variable sea ice extent, productivity and polynyas present even during the Last Glacial Maximum. Four layers of ice‐rafted debris were also identified and correlated within the cores. These events occurred ca at 9, 24 to 25, 26 to 27 and 43 ka, asynchronous with Heinrich layers in the wider north‐east Atlantic and here interpreted as reflecting instability in the Svalbard/Barents Ice sheet and the northward advection of warm Atlantic water during the Late Weichselian. The activity of the ancestral West Spitsbergen Current is interpreted using mean sortable silt records from the cores. On the Vestnesa Ridge drift the modern mass accumulation rate, calculated using excess ²¹⁰Pb, is 0·076 g cm^?2 year^?1. On the Yermak Plateau slope the modern mass accumulation rate is 0·053 g cm^?2 year^?1.     

Fluid and source magma evolution of the Questa porphyry Mo deposit,New Mexico,USA

Combined fluid inclusion microthermometry and microanalysis by laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) are used to constrain the hydrothermal processes forming a typical Climax-type porphyry Mo deposit. Molybdenum mineralisation at Questa occurred in two superimposed hydrothermal stages, a magmatic-hydrothermal breccia and later stockwork veining. In both stages, texturally earliest fluids were single-phase, of low salinity (~7 wt.% NaCl_equiv.) and intermediate-density. Upon decompression to ~300 bar, they boiled off a vapour phase, leaving behind a residual brine (up to 45 wt.% NaCl_equiv) at temperatures of ~420°C. The highest average Mo concentrations in this hot brine were ~500 μg/g, exceeding the Mo content of the intermediate-density input fluid by about an order of magnitude and reflecting pre-concentration of Mo by fluid phase separation prior to MoS₂ deposition from the brine. Molybdenum concentrations in brine inclusions, then, decrease down to 5 μg/g, recording Mo precipitation in response to cooling of the saline liquid to ~360°C. Molybdenite precipitation from a dense, residual and probably sulphide-depleted brine is proposed to explain the tabular shape of the ore body and the absence of Cu-Fe sulphides in contrast to the more common Cu-Mo deposits related to porphyry stocks. Cesium and Rb concentrations in the single-phase fluids of the breccia range from 2 to 8 and from 40 to 65 μg/g, respectively. In the stockwork veins, Cs and Rb concentrations are significantly higher (45–90 and 110–230 μg/g, respectively). Because Cs and Rb are incompatible and hydrothermally non-reactive elements, the systematic increase in their concentration requires two distinct pulses of fluid exsolution from a progressively more fractionated magma. By contrast, major element and ore metal concentrations of these two fluid pulses remain essentially constant. Mass balance calculations using fluid chemical data from LA-ICPMS suggest that at least 25 km³ of melt and 7 Gt of deep input fluid were necessary to provide the amount of Mo contained in the stockwork vein stage alone. While the absolute amounts of fluid and melt are uncertain, the well-constrained element ratios in the fluids together with empirical fluid/melt partition coefficients derived from the inclusion analyses suggest a high water content of the source melt of ~10%. In line with other circumstantial evidence, these results suggest that initial fluid exsolution may have occurred at a confining pressure exceeding 5 kbar. The source of the molybdenum-mineralising fluids probably was a particularly large magma chamber that crystallised and fractionated in the lower crust or at mid-crustal level, well below the shallow intrusions immediately underlying Questa and other porphyry molybdenum deposits. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Age and depositional environment of the Draa Sfar massive sulfide deposit, Morocco

The Draa Sfar mineralization consists of two main stratabound orebodies, Sidi M’Barek and Tazacourt, located north and south of the Tensift River (“Oued Tessift”), respectively. Each orebody is comprised by at least two massive sulfide lenses. The hosting rocks are predominantly black shales, although minor rhyolitic rocks are also present in the footwall to the southern orebody. Shales, rhyolitic volcanic rocks, and massive sulfides are all included into the Sarhlef Series, which is recognized as one of the main stratigraphic units of the Moroccan Variscan Meseta. Hydrothermal activity related with an anomalous thermal gradient, together with a high sedimentation rate in a tectonically driven pull-apart marine basin, favored the accumulation of organic-rich mud in the deepest parts of the basin and the sedimentary environment suitable for massive sulfide deposition and preservation. This took place by replacement of the hosting unlithified wet mud below the sediment–water interface. Geochemical data suggest a sedimentary environment characterized by oxic water column and anoxic sediment pile with the redox boundary below the sediment–water interface. The low oxygen availability within the sediment pile inhibited oxidation and pyritization of pyrrhotite. Biostratigraphic analysis, based on the palynological content of the hosting black shales, restricts the age of the sulfides to the Asbian substage (mid-Mississippian). This age is consistent with earlier geochronological constraints.     

Detrital zircon provenance of Neoproterozoic to Cenozoic deposits in Iran: Implications for chronostratigraphy and collisional tectonics

Ion-microprobe U–Pb analyses of 589 detrital zircon grains from 14 sandstones of the Alborz mountains, Zagros mountains, and central Iranian plateau provide an initial framework for understanding the Neoproterozoic to Cenozoic provenance history of Iran. The results place improved chronological constraints on the age of earliest sediment accumulation during Neoproterozoic–Cambrian time, the timing of the Mesozoic Iran–Eurasia collision and Cenozoic Arabia–Eurasia collision, and the contribution of various sediment sources of Gondwanan and Eurasian affinity during opening and closure of the Paleotethys and Neotethys oceans. The zircon age populations suggest that deposition of the extensive ~ 1 km-thick clastic sequence at the base of the cover succession commenced in latest Neoproterozoic and terminated by Middle Cambrian time. Comparison of the geochronological data with detrital zircon ages for northern Gondwana reveals that sediment principally derived from the East African orogen covered a vast region encompassing northern Africa and the Middle East. Although most previous studies propose a simple passive-margin setting for Paleozoic Iran, detrital zircon age spectra indicate Late Devonian–Early Permian and Cambrian–Ordovician magmatism. These data suggest that Iran was affiliated with Eurasian magmatic arcs or that rift-related magmatic activity during opening of Paleotethys and Neotethys was more pronounced than thought along the northern Gondwanan passive-margin. For a Triassic–Jurassic clastic overlap assemblage (Shemshak Formation) in the Alborz mountains, U–Pb zircon ages provide chronostratigraphic age control requiring collision of Iran with Eurasia by late Carnian–early Norian time (220–210 Ma). Finally, Cenozoic strata yield abundant zircons of Eocene age, consistent with derivation from arc magmatic rocks related to late-stage subduction and/or breakoff of the Neotethys slab. Together with the timing of foreland basin sedimentation in the Zagros, these detrital zircon ages help bracket the onset of the Arabia–Eurasia collision in Iran between middle Eocene and late Oligocene time.