Rare earth element geochemistry of Late Devonian reefal carbonates, canning basin, Western Australia: Confirmation of a seawater REE proxy in ancient limestones

Rare earth element and yttrium (REE+Y) concentrations were determined in 49 Late Devonian reefal carbonates from the Lennard Shelf, Canning Basin, Western Australia. Shale-normalized (SN) REE+Y patterns of the Late Devonian samples display features consistent with the geochemistry of well-oxygenated, shallow seawater. A variety of different ancient limestone components, including microbialites, some skeletal carbonates (stromatoporoids), and cements, record seawater-like REE+Y signatures. Contamination associated with phosphate, Fe-oxides and shale was tested quantitatively, and can be discounted as the source of the REE+Y patterns. Co-occurring carbonate components that presumably precipitated from the same seawater have different relative REE concentrations, but consistent REE+Y patterns. Clean Devonian early marine cements (n = 3) display REE+Y signatures most like that of modern open ocean seawater and the highest Y/Ho ratios (e.g., 59) and greatest light REE (LREE) depletion (average Nd-SN/Yb-SN = 0.413, SD = 0.076). However, synsedimentary cements have the lowest REE concentrations (e.g., 405 ppb). Non-contaminated Devonian microbialite samples containing a mixture of the calcimicrobe Renalcis and micritic thrombolite aggregates in early marine cement (n = 11) have the highest relative REE concentrations of tested carbonates (average total REE = 11.3 ppm). Stromatoporoid skeletons, unlike modern corals, algae and molluscs, also contain well-developed, seawater-like REE patterns. Samples from an estuarine fringing reef have very different REE+Y patterns with LREE enrichment (Nd-SN/Yb-SN > 1), possibly reflecting inclusion of estuarine colloidal material that contained preferentially scavenged LREE from a nearby riverine input source. Hence, Devonian limestones provide a proxy for marine REE geochemistry and allow the differentiation of co-occurring water masses on the ancient Lennard Shelf. Although appropriate partition coefficients for quantification of Devonian seawater REE concentrations from out data are unknown, hypothetical Devonian Canning Basin seawater REE patterns were obtained with coefficients derived from modern natural proxies and experimental values. Resulting Devonian seawater patterns are slightly enriched in LREE compared to most modem seawaters and suggest higher overall REE concentrations, but are very similar to seawaters from regions with high terrigenous inputs. Our results suggest that most limestones should record important aspects of the REE geochemistry of the waters in which they precipitated, provided they are relatively free of terrigenous contamination and major diagenetic alteration from fluids with high, non-seawater-like REE contents. Hence, we expect that many other ancient limestones will serve as seawater REE proxies, and thereby provide information on paleoceanography, paleogeography and geochemical evolution of the oceans. Copyright (C) 2004 Elsevier Ltd.

Geological and trace element evidence for a marine sedimentary environment of deposition and biogenicity of 3.45 Ga stromatolitic carbonates in the Pilbara Craton, and support for a reducing Archaean ocean

Bedded carbonate rocks from the 3.45 Ga Warrawoona Group, Pilbara Craton, contain structures that have been regarded either as the oldest known stromatolites or as abiotic hydrothermal deposits. We present new field and petrological observations and high-precision REE + Y data from the carbonates in order to test the origin of the deposits. Trace element geochemistry from a number of laminated stromatolitic dolomite samples of the c. 3.40 Ga Strelley Pool Chert conclusively shows that they precipitated from anoxic seawater, probably in a very shallow environment consistent with previous sedimentological observations. Edge-wise conglomerates in troughs between stromatolites and widespread cross-stratification provide additional evidence of stromatolite construction, at least partly, from layers of particulate sediment, rather than solely from rigid crusts. Accumulation of particulate sediment on steep stromatolite sides in a high-energy environment suggests organic binding of the surface. Relative and absolute REE + Y contents are exactly comparable with Late Archaean microbial carbonates of widely agreed biological origin. Ankerite from a unit of bedded ankerite–chert couplets from near the top of the stratigraphically older (3.49 Ga) Dresser Formation, which immediately underlies wrinkly stromatolites with small, broad, low-amplitude domes, also precipitated from anoxic seawater. The REE + Y data of carbonates from the Strelley Pool Chert and Dresser Formation contrast strongly with those from siderite layers in a jasper–siderite–Fe-chlorite banded iron-formation from the base of the Panorama Formation (3.45 Ga), which is clearly hydrothermal in origin. The geochemical results, together with sedimentological data, strongly support: (1) deposition of Dresser Formation and Strelley Pool Chert carbonates from Archaean seawater, in part as particulate carbonate sediment; (2) biogenicity of the stromatolitic carbonates; (3) a reducing Archaean atmosphere; (4) ongoing extensive terrestrial erosion prior to ∼3.45 Ga.

Petrology and geochemistry of early cretaceous bimodal continental flood volcanism of the NW Etendeka, Namibia. Part 1: Introduction, mafic lavas and re-evaluation of mantle source components

The bimodal NW Etendeka province is located at the continental end of the Tristan plume trace in coastal Namibia. It comprises a high-Ti (Khumib type) and three low-Ti basalt (Tafelberg, Kuidas and Esmeralda types) suites, with, at stratigraphically higher level, interstratified high-Ti latites (three units) and quartz latites (five units), and one low-Ti quartz latite. Khumib basalts are enriched in high field strength elements and light rare earth elements relative to low-Ti types and exhibit trace element affinities with Tristan da Cunha lavas. The unradiogenic Pb-206/Pb-204 ratios of Khumib basalts are distinctive, most plotting to the left of the 132 Ma Geochron, together with elevated Pb-207/Pb-204 ratios, and Sr-Nd isotopic compositions plotting in the lower Nd-143/Nd-144 part of mantle array (EM1-like). The low-Ti basalts have less coherent trace element patterns and variable, radiogenic initial Sr (similar to0.707-0.717) and Pb isotope compositions, implying crustal contamination. Four samples, however, have less radiogenic Pb and Sr that we suggest approximate their uncontaminated source. All basalt types, but particularly the low-Ti types, contain samples with trace element characteristics (e.g. Nb/Nb-*) suggesting metasediment input, considered source-related. Radiogenic isotope compositions of these samples require long-term isolation of the source in the mantle and depletions (relative to unmodified sediment) in certain elements (e.g. Cs, Pb, U), which are possibly subduction-related. A geodynamic model is proposed in which the emerging Tristan plume entrained subducted material in the Transition Zone region, and further entrained asthenosphere during plume head expansion. Mixing calculations suggest that the main features of the Etendeka basalt types can be explained without sub-continental lithospheric mantle input. Crustal contamination is evident in most low-Ti basalts, but is distinct from the incorporation of a metasedimentary source component at mantle depths.

Assigning dates to thin gneissic veins in high-grade metamorphic terranes: A cautionary tale from Akilia, southwest Greenland

A granodiorite from Akilia, southwest Greenland, previously suggested to date putative life-bearing rocks to greater than or equal to3.84 Ga, is re-investigated using whole-rock major and trace-element geochemistry, and detailed cathodoluminescence image-guided secondary ion mass spectrometer analyses of zircon U-Th-Pb and rare earth elements. Complex zircon internal structure reveals three episodes of zircon growth and/or recrystallization dated to c. 3.84 Ga, 3.62 Ga and 2.71 Ga. Rare earth element abundances imply a significant role for garnet in zircon generation at 3.62 Ga and 2.71 Ga. The 3.62 Ga event is interpreted as partial melting of a c. 3.84 Ga grey gneiss precursor at granulite facies with residual garnet. Migration of this 3.62 Ga magma (or melt-crystal mush) away from the melt source places a maximum age limit on any intrusive relationship. These early Archaean relationships have been complicated further by isotopic reworking in the 2.71 Ga event, which could have included a further episode of partial melting. This study highlights a general problem associated with dating thin gneissic veins in polyphase metamorphic terranes, where field relationships may be ambiguous and zircon inheritance can be expected.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

The geochemistry of late Archaean microbial carbonate: Implications for ocean chemistry and continental erosion history

Trace element concentrations and combined Sr- and Nd-isotope compositions were determined on stromatolitic carbonates (microbialites) from the 2.52 Ga Campbellrand carbonate platform (South Africa). Shale-normalised rare earth element and yttrium patterns of the ancient samples are similar to those of modern seawater in having positive La and Y anomalies and in being depleted in light rare earth elements. In contrast to modem seawater (and microbialite proxies), the 2.52 Ga samples lack a negative Ce anomaly but possess a positive Eu anomaly. These latter trace element characteristics are interpreted to reflect anoxic deep ocean waters where, unlike today, hydrothermal Fe input was not oxidised, and scavenged and rare earth elements were not coprecipitated with Fe-oxyhydroxides. The persistence of a positive Eu anomaly in relatively shallow Campbellrand platform waters indicates a dramatic reversal from hydrothermally dominated (Archaean) to continental erosion-dominated (Phanerozoic) rare earth element flux ratio. The dominant hydrothermal input is also expressed in the initial Sr- and Nd-isotope ratios. There is collinear variation in Sr-Nd systematics, which range from primitive values (Sr-87/Sr-86 of 0.702386 and epsilon (Nd) of +2.1) to more evolved crustal ratios. Mixing calculations show that the range in trace element ratios (e.g., Y/Ho) and initial isotope ratios is not a result of contamination by trapped sediment, but that the chemical band isotopic variation reflects carbonate deposition in an environment where different water masses mixed. Calculated Nd flux ratios yield a hydrothermal input into the 2.52 Ga oceans one order of magnitude larger than continental input. Such a change in flux ratio most likely required substantially reduced continental inputs, which could, in turn, reflect a plate tectonic causation (e.g., reduced topography or expansion of epicontinental seas). Copyright (C) 2001 Elsevier Science Ltd.

Fluid-mobile trace element constraints on the role of slab melting and implications for Archaean crustal growth models

We use published and new trace element data to identify element ratios which discriminate between arc magmas from the supra-subduction zone mantle wedge and those formed by direct melting of subducted crust (i.e. adakites). The clearest distinction is obtained with those element ratios which are strongly fractionated during refertilisation of the depleted mantle wedge, ultimately reflecting slab dehydration. Hence, adakites have significantly lower Pb/Nd and B/Be but higher Nb/Ta than typical arc magmas and continental crust as a whole. Although Li and Be are also overenriched in continental crust, behaviour of Li/Yb and Be/Nd is more complex and these ratios do not provide unique signatures of slab melting. Archaean tonalite-trondhjemite-granodiorites (TTGs) strongly resemble ordinary mantle wedge-derived arc magmas in terms of fluid-mobile trace element content, implying that they-did not form by slab melting but that they originated from mantle which was hydrated and enriched in elements lost from slabs during prograde dehydration. We suggest that Archaean TTGs formed by extensive fractional crystallisation from a mafic precursor. It is widely claimed that the time between the creation and subduction of oceanic lithosphere was significantly shorter in the Archaean (i.e. 20 Ma) than it is today. This difference was seen as an attractive explanation for the presumed preponderance of adakitic magmas during the first half of Earth's history. However, when we consider the effects of a higher potential mantle temperature on the thickness of oceanic crust, it follows that the mean age of oceanic lithosphere has remained virtually constant. Formation of adakites has therefore always depended on local plate geometry and not on potential mantle temperature.