Rare earth element fractionation in authigenic illite–smectite from Late Permian clastic rocks, Bowen Basin, Australia: implications for physico-chemical environments of fluids during illitization

REE analyses were performed on authigenic illitic clay. minerals from Late Permian mudrocks, sandstones and bentonites from the Bowen Basin (Australia). The mixed-layer illite-smectite exhibit REE patterns with an obvious fractionation of the HREE from the LREE and MREE, which is an apparent function of degree of illitization reaction. The highly illitic (R greater than or equal to 3) illite-smectite from the northern Bowen Basin show a depletion of LREE relative to the less illitic (R=0 and 1) clays. In contrast, an enrichment of HREE for the illite-rich clays relative to less. illitic clays is evident for the southern Bowen Basin samples. The North American Shale Composite-normalized (La/Lu)(sn) ratios show negative correlations with the illite content in illite-smectite and positive correlations with the delta(18)O values of the clays for both the northern and southern Bowen Basin samples. These correlations indicate that the increasing depletion of LREE in hydrothermal fluids is a function of increasing water/rock ratios in the northern Bowen Basin. Good negative correlations between (La/Lu)(sn) ratios and illite content in illite-smectite from the southern Bowen Basin suggest the involvement of fluids with higher alkalinity and higher pH in low water/ rock ratio conditions. Increasing HREE enrichment with delta(18)O decrease indicates the effect of increasing temperature at low water/rock ratios in the southern Bowen Basin. Results of the present study confirm the conclusions of some earlier studies suggesting that REE in illitic clay minerals are mobile and fractionated during illitization and that this fact should be considered in studies of sedimentary processes and in identifying provenance. Moreover, our results show that REE systematic of illitic clay minerals can be applied as an useful technique to gain information about physico-chemical conditions during thermal and fluid flow events in certain sedimentary basins. (C) 2003 Elsevier Science B.V. All rights reserved.

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.

Chemical characterization of earth's most ancient clastic metasediments from the Isua Greenstone Belt, southern West Greenland

New and published major and trace element abundances of elastic metasediments (mainly garnet-biotite-plagioclase schists) from the similar to 3.8 Ga Isua Greenstone Belt (IGB), southern West Greenland, are used in an attempt to identify the compositional characteristics of the protoliths of these sediments. Compositionally, the metasediments are heterogeneous with enrichment of LREE (La/Sm-chord = 1.1-3.9) and variable enrichment and depletion of HREE (Gd/Yb-chord = 0.8-4.3). Chondrite-normalized Eu is also variable, spanning a range from relative Eu depletion to enrichment (Eu/Eu* = 0.6-1.3). A series of geochemical and geological criteria provides conclusive evidence for a sedimentary origin, in disagreement with some previous studies that questioned the presence of genuine elastic metasediments. In particular, trace element systematics of IGB metasediments show strong resemblance to other well-documented Archaean clastic sediments, and are consistent with a provenance consisting of ultramafic, malic and felsic igneous rocks. Two schists, identified as metasomatized mafic igneous rocks from petrographic and field evidence, show distinct compositional differences to the metasediments. Major element systematics document incipient-to-moderate source weathering in the majority of metasediments, while signs of secondary K-addition are rare. Detailed inspection of Eu/Eu*, Fe2O3 and CIW (chemical index of weathering) relationships reveals that elevated iron contents (when compared to averages for continental crust) and strong relative enrichment in Eu may be due to precipitation of marine Fe-oxyhydroxides during deposition or diagenesis on the seafloor. Some of the IGB metasediments have yielded anomalous Nd-142 and W-182 isotopic compositions that were respectively interpreted in terms of early mantle differentiation processes and the presence of a meteorite component. Alternatively, W and possibly Nd isotopes could have been affected by thermal neutron capture on the Hadean surface. The latter process was tested in this study by analysis of Sm isotope compositions, which serve as an effective monitor for neutron capture effects. As no anomalous variation from terrestrial values was detected, we infer that isotope systematics (including W-182 and Nd-142) of IGB metasediments were not affected by neutron capture, but reflect decay of radioactive parent isotopes. Copyright (c) 2005 Elsevier Ltd.

Rare earth element and yttrium variability in South East Queensland waterways

We present data for the rare earth elements and yttrium (REY) in the National Research Council of Canada natural river water reference material SLRS-4 and 19 natural river waters from small catchments in South-East Queensland, Australia, by a direct ICP-MS method. The 0.22 mu m filtered river water samples show a large degree of variability in both the REY concentration, e.g., La varies from 13 to 1157 ppt, and shape of the alluvial-sediment-normalised REY patterns with different samples displaying light, middle or heavy rare earth enrichment. In addition, a spatial study was undertaken along the freshwater section of Beerburrum Creek, which demonstrates that similar to 75% of the total REYs in this waterway are removed prior to estuarine mixing without evidence of fractionation.

Direct quantification of rare earth element concentrations in natural waters by ICP-MS

A direct quadrupole ICP-MS technique has been developed for the analysis of the rare earth elements and yttrium in natural waters. The method has been validated by comparison of the results obtained for the river water reference material SLRS-4 with literature values. The detection limit of the technique was investigated by analysis of serial dilutions of SLRS-4 and revealed that single elements can be quantified at single-digit fg/g concentrations. A coherent normalised rare earth pattern was retained at concentrations two orders of magnitude below natural concentrations for SLRS-4, demonstrating the excellent inter-element accuracy and precision of the method. The technique was applied to the analysis of a diluted mid-salinity estuarine sample, which also displayed a coherent normalised rare earth element pattern, yielding the expected distinctive marine characteristics. (c) 2006 Published by Elsevier Ltd.

Implementation of particle-scale rotation in the 3-d Lattice Solid Model

In this study, 3-D Lattice Solid Model (LSMearth or LSM) was extended by introducing particle-scale rotation. In the new model, for each 3-D particle, we introduce six degrees of freedom: Three for translational motion, and three for orientation. Six kinds of relative motions are permitted between two neighboring particles, and six interactions are transferred, i.e., radial, two shearing forces, twisting and two bending torques. By using quaternion algebra, relative rotation between two particles is decomposed into two sequence-independent rotations such that all interactions due to the relative motions between interactive rigid bodies can be uniquely decided. After incorporating this mechanism and introducing bond breaking under torsion and bending into the LSM, several tests on 2-D and 3-D rock failure under uni-axial compression are carried out. Compared with the simulations without the single particle rotational mechanism, the new simulation results match more closely experimental results of rock fracture and hence, are encouraging. Since more parameters are introduced, an approach for choosing the new parameters is presented.

The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.