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.

Diplacanthid acanthodians from the Aztec Siltstone (late Middle Devonian) of southern Victoria Land, Antartica

One articulated and several partial, semi-articulated specimens of acanthodians were collected in 1970 from the freshwater deposits of the Aztec Siltstone (Middle Devonian; Givetian), Portal Mountain, southern Victoria Land, Antarctica, during a Victoria University of Wellington Antarctic Expedition. The Portal Mountain fish fauna, preserved in a finely laminated, non-calcareous siltstone, includes acanthodians, palaeoniscoids, and bothriolepid placoderms. The articulated acanthodian specimens are the most complete fossil fish remains documented so far from the Aztec assemblage, which is the most diverse fossil vertebrate fauna known from Antarctica. They are described as a new taxon, Milesacanthus antarctica gen. et sp. nov., which is assigned to the family Diplacanthidae. Its fin spines show some similarities to spine fragments named Byssacanthoides debenhami from glacial moraine at Granite Harbour, Antarctica, and much larger spines named Antarctonchus glacialis from outcrops of the Aztec Siltstone in the Boomerang Range, southern Victoria Land. Both of these are reviewed, and retained as form taxa for isolated spines. Various isolated remains of fin spines and scales are described from Portal Mountain and Mount Crean (Lashly Range), and referred to Milesacanthus antarctica gen. et sp. nov. The histology of spines and scales is documented for the first time, and compared with acanthodian material from the Devonian of Australia and Europe. Distinctive fin spines from Mount Crean are provisionally assigned to Culmacanthus antarctica Young, 1989b. Several features on the most complete of the new fish specimens - in particular, the apparent lack of an enlarged cheek plate - suggest a revision of the diagnosis for the Diplacanthidae.

Geochemistry of late Archaean stromatolites from Zimbabwe: evidence for microbial life in restricted epicontinental seas

We report comprehensive trace element and Sr-isotope data for microbial carbonates from the Archaean Mushandike limestone, Masvingo Greenstone Belt, Zimbabwe. The stromatolites have very coherent REE + Y patterns and share the essential shale-normalised characteristics of other Archaean marine precipitates (positive La and Gd anomalies, absence of a negative Cc anomaly and a strongly superchondritic Y/Ho ratio). Mixing models constrain the maximum amount of shale contamination to 0.25-1% and calculated detritus-free carbonate REE + Y systematics require precipitation from seawater. In terms of light-REE over heavy-REE depletion, however, the studied samples are very different from all other known Archaean marine precipitates. In shale-normalised plots, the Mushandike samples yield a negative slope. A very restricted, regional input source of the dissolved load is indicated because normalisation with locally occurring tonalite gneiss REE + Y data yields a pattern closely resembling typical shale-normalised Archaean marine chemical sediments. The disappearance of a negative Eu anomaly when patterns are normalised with local tonalite gneiss strengthens this interpretation. Sr-isotope ratios are strongly correlated with trace element contents and ratios, which explains the modest scatter in Sr-isotope ratios as representing (minor) clastic contamination. Importantly, even the least contaminated samples have very radiogenic initial Sr-87/Sr-86 ratios (0.7184) implying Sr input from an ancient high Rb/Sr source, such as the early Archaean gneisses of south-central Zimbabwe. A local ancient (3.5-3.8 Ga) source is also indicated by previously published Pb-isotope datasets for the Mushandike stromatolites. This is entirely compatible with the occurrence of 3.7-3.8 Ga zircons in quartzites and metapelites from comparably old greenstone belts within less than 150 km of the studied locality. Comparison of the Pb-isotope ratios of the Mushandike stromatolites with 2.7 and 2.6 Ga old stromatolites from the neighbouring, Belingwe Greenstone Belt demonstrates differences in initial isotope composition that relate to the extent of exchange with the open ocean. The development of numerous basins on old continental crust, with water masses variably restricted from the open ocean. suggests a lack of strong vertical topography on this late Archaean craton. (C) 2004 Elsevier B.V. All rights reserved.

Volcanic resurfacing and the early terrestrial crust: Zircon U-Pb and REE constraints from the Isua Greenstone Belt, southern West Greenland

The oldest known bona fide succession of elastic metasediments Occurs in the Isua Greenstone Belt. SW Greenland and consists of a variety of mica schists and rare metaconglomerates. The metasediments are in direct contact with a felsic metavolcanic lithology that has previously been dated to 3.71 Ga. Based on trace element geochemical data for 30 metasediments, we selected the six samples with highest Zr concentrations for zircon extraction. These samples all yielded very few or no zircon, Those extracted from mica schists yielded ion probe U/Pb ages between 3.70 and 3,71 Ga. One metaconglomerate sample yielded just a single zircon of 3.74 Ga age. The mica schist hosted zircons have U/Pb ages. Th/U ratios, REE patterns and Eu anomalies indistinguishable from zircon in the adjacent 3.71 Ga felsic metavolcanic unit. Trace element modelling requires the bulk of material in the metasediments to be derived from variably weathered mafic lithologies but some metasediments contain substantial contribution from more evolved source lithologies. The paucity of zircon in the mica schists is thus explained by incorporation of material from largely zircon-free volcanic lithologies. The absence of older zircon in the mica schists and the preponderance of mafic source material imply intense, mainly basaltic resurfacing of the early Earth. The implications of this process are discussed, Thermal considerations suggest that horizontal growth of Hadean crust by addition of mafic ultramafic lavas must have triggered self-reorganisation of the protocrust by remelting. Reworking oft Hadean crust may have been aided by burial of hydrated (weathered) metabasalt due to semi-continuous addition of new voluminous basalt Outpouring,;, This process Causes a bias towards eruption of Zr-saturated partial melts at the surface with O-isotope corn posit ion,, potentially different from the mantle. The oldest zircons hosted in sediments would have been buried to substantial depth or formed in plutons that crystallised at some depth from which it took hundreds of millions of years for them to be exhumed and incorporated into much younger sediments. (C) 2005 Elsevier B.V.All rights reserved.

Integrated Pb- and S-isotope investigation of sulphide minerals from the early Archaean of southwest Greenland

We present high-spatial resolution secondary ion mass spectrometry (SIMS) measurements of Pb and S isotopes in sulphides from early Archaean samples at two localities in southwest Greenland. Secondary pyrite from a 3.71 Ga sample of magnetite-quartz banded iron formation in the Isua Greenstone Belt, which has previously yielded unradiogenic Pb consistent with its ancient origin, contains sulphur with a mass independently fractionated (MIF) isotope signature (Delta(33)S =+3.3 parts per thousand). This reflects the secondary mineralization of remobilized sedimentary S carrying a component modified by photochemical reactions in the early Archaean atmosphere. It further represents one of the most extreme positive excursions so far known from the early Archaean rock record. Sulphides from a quartz-pyroxene rock and an ultramafic boudin from the island of Akilia, in the Godth (a) over circle bsfjord, have heterogeneous and generally radiogenic Pb isotopic compositions that we interpret to represent partial re-equilibration of Pb between the sulphides and whole rocks during tectonothermal events at 3.6, 2.7 and 1.6 Ga. Both these samples have Delta(33)S=0 (within analytical error) and therefore show no evidence for MIF sulphur. These data are consistent with previous interpretations that the rock cannot be proven to have a sedimentary origin. Our study illustrates that SIMS S-isotope measurements in ancient rocks can be used to elucidate early atmospheric parameters because of the ability to obtain combined S and Pb-isotope data, but caution must be applied when using such data to infer protolith. When information from geological context, petrography and chronology (i.e. by Pb isotopes) is combined and fully evaluated, Delta(33)S signatures from sulphides and their geological significance can be interpreted with a higher degree of confidence. (c) 2005 Elsevier B.V All rights reserved.

Multiscale coupling and multiphysics approaches in earth sciences : applications

Geoscientists are confronted with the challenge of assessing nonlinear phenomena that result from multiphysics coupling across multiple scales from the quantum level to the scale of the earth and from femtoseconds to the 4.5 Ga of history of our planet. We neglect in this review electromagnetic modelling of the processes in the Earth’s core, and focus on four types of couplings that underpin fundamental instabilities in the Earth. These are thermal (T), hydraulic (H), mechanical (M) and chemical (C) processes which are driven and controlled by the transfer of heat to the Earth’s surface. Instabilities appear as faults, folds, compaction bands, shear/fault zones, plate boundaries and convective patterns. Convective patterns emerge from buoyancy overcoming viscous drag at a critical Rayleigh number. All other processes emerge from non-conservative thermodynamic forces with a critical critical dissipative source term, which can be characterised by the modified Gruntfest number Gr. These dissipative processes reach a quasi-steady state when, at maximum dissipation, THMC diffusion (Fourier, Darcy, Biot, Fick) balance the source term. The emerging steady state dissipative patterns are defined by the respective diffusion length scales. These length scales provide a fundamental thermodynamic yardstick for measuring instabilities in the Earth. The implementation of a fully coupled THMC multiscale theoretical framework into an applied workflow is still in its early stages. This is largely owing to the four fundamentally different lengths of the THMC diffusion yardsticks spanning micro-metre to tens of kilometres compounded by the additional necessity to consider microstructure information in the formulation of enriched continua for THMC feedback simulations (i.e., micro-structure enriched continuum formulation). Another challenge is to consider the important factor time which implies that the geomaterial often is very far away from initial yield and flowing on a time scale that cannot be accessed in the laboratory. This leads to the requirement of adopting a thermodynamic framework in conjunction with flow theories of plasticity. This framework allows, unlike consistency plasticity, the description of both solid mechanical and fluid dynamic instabilities. In the applications we show the similarity of THMC feedback patterns across scales such as brittle and ductile folds and faults. A particular interesting case is discussed in detail, where out of the fluid dynamic solution, ductile compaction bands appear which are akin and can be confused with their brittle siblings. The main difference is that they require the factor time and also a much lower driving forces to emerge. These low stress solutions cannot be obtained on short laboratory time scales and they are therefore much more likely to appear in nature than in the laboratory. We finish with a multiscale description of a seminal structure in the Swiss Alps, the Glarus thrust, which puzzled geologists for more than 100 years. Along the Glarus thrust, a km-scale package of rocks (nappe) has been pushed 40 km over its footwall as a solid rock body. The thrust itself is a m-wide ductile shear zone, while in turn the centre of the thrust shows a mm-cm wide central slip zone experiencing periodic extreme deformation akin to a stick-slip event. The m-wide creeping zone is consistent with the THM feedback length scale of solid mechanics, while the ultralocalised central slip zones is most likely a fluid dynamic instability.

Challenges in Earth sciences: the 21st century

The 21st century poses many challenges for global sustainability. Among them, most importantly, the human race will encounter scarcity of raw materials and conventional energy resources. And, India may have to take the brunt of these problems as it is going to be the most populated region of the world with concomitant increase in energy demand and requirement of other resources. India will be the testing ground for introducing newer ways of green technology and innovative principles of resource management and utilization. With the vagaries of potential climate change gathering clouds in the background, Earth sciences will have a special and predominant role in guiding the society in prioritizing our resource discovery, utilization and their consumption and the upkeep of environment. On the fundamental level, Earth sciences are going through a most exciting phase of development as a born-again science. Technological breakthroughs including the satellite-based observations augur well for gaining new insights into Earth processes. A set of exciting fundamental problems that are globally identified will set the stage for an exhilarating period of new discoveries. Improvements in numerical and computer-based techniques will assist in modelling of Earth processes to unprecedented levels. India will have to take special effort in improving the existing experimentation facilities in the Earth science departments of the country, and also the general level of Earth science education to meet the global standards. This article presents an Earth science vision for the 21st century in an Indian context.

Best practices in prediction for decision-making : lessons from the atmospheric and earth sciences

Predictions that result from scientific research hold great appeal for decision-makers who are grappling with complex and controversial environmental issues, by promising to enhance their ability to determine a need for and outcomes of alternative decisions. A problem exists in that decision-makers and scientists in the public and private sectors solicit, produce, and use such predictions with little understanding of their accuracy or utility, and often without systematic evaluation or mechanisms of accountability. In order to contribute to a more effective role for ecological science in support of decision-making, this paper discusses three ``best practices'' for quantitative ecosystem modeling and prediction gleaned from research on modeling, prediction, and decision-making in the atmospheric and earth sciences. The lessons are distilled from a series of case studies and placed into the specific context of examples from ecological science.

Multiscale coupling and multiphysics approaches in earth sciences : theory

Measuring Earth material behaviour on time scales of millions of years transcends our current capability in the laboratory. We review an alternative path considering multiscale and multiphysics approaches with quantitative structure-property relationships. This approach allows a sound basis to incorporate physical principles such as chemistry, thermodynamics, diffusion and geometry-energy relations into simulations and data assimilation on the vast range of length and time scales encountered in the Earth. We identify key length scales for Earth systems processes and find a substantial scale separation between chemical, hydrous and thermal diffusion. We propose that this allows a simplified two-scale analysis where the outputs from the micro-scale model can be used as inputs for meso-scale simulations, which then in turn becomes the micro-model for the next scale up. We present two fundamental theoretical approaches to link the scales through asymptotic homogenisation from a macroscopic thermodynamic view and percolation renormalisation from a microscopic, statistical mechanics view.

Entropic bounds for multi-scale and multi-physics coupling in earth sciences

The ability to understand and predict how thermal, hydrological,mechanical and chemical (THMC) processes interact is fundamental to many research initiatives and industrial applications. We present (1) a new Thermal– Hydrological–Mechanical–Chemical (THMC) coupling formulation, based on non-equilibrium thermodynamics; (2) show how THMC feedback is incorporated in the thermodynamic approach; (3) suggest a unifying thermodynamic framework for multi-scaling; and (4) formulate a new rationale for assessing upper and lower bounds of dissipation for THMC processes. The technique is based on deducing time and length scales suitable for separating processes using a macroscopic finite time thermodynamic approach. We show that if the time and length scales are suitably chosen, the calculation of entropic bounds can be used to describe three different types of material and process uncertainties: geometric uncertainties,stemming from the microstructure; process uncertainty, stemming from the correct derivation of the constitutive behavior; and uncertainties in time evolution, stemming from the path dependence of the time integration of the irreversible entropy production. Although the approach is specifically formulated here for THMC coupling we suggest that it has a much broader applicability. In a general sense it consists of finding the entropic bounds of the dissipation defined by the product of thermodynamic force times thermodynamic flux which in material sciences corresponds to generalized stress and generalized strain rates, respectively.