Introduction to Special issue on the Symposium on the Application of Neural Networks to the Earth Sciences

Papers in this issue of Natural Resources Research are from the “Symposium on the Application of Neural Networks to the Earth Sciences,” held 20–21 August 2002 at NASA Moffet Field, Mountain View, California. The Symposium represents the Seventh International Symposium on Mineral Exploration (ISME-02). It was sponsored by the Mining and Materials Processing Institute of Japan (MMIJ), the US Geological Survey, the Circum-Pacific Council, and NASA. The ISME symposia have been held every two years in order to bring together scientists actively working on diverse quantitative methods applied to the earth sciences. Although the title, International Symposium on Mineral Exploration, suggests exclusive focus on mineral exploration, interests and presentations always have been wide-ranging—talks presented at this symposium are no exception.

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.

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.

Towards a Modeling Environment for Earth Systems Simulations

The developments of models in Earth Sciences, e.g. for earthquake prediction and for the simulation of mantel convection, are fare from being finalized. Therefore there is a need for a modelling environment that allows scientist to implement and test new models in an easy but flexible way. After been verified, the models should be easy to apply within its scope, typically by setting input parameters through a GUI or web services. It should be possible to link certain parameters to external data sources, such as databases and other simulation codes. Moreover, as typically large-scale meshes have to be used to achieve appropriate resolutions, the computational efficiency of the underlying numerical methods is important. Conceptional this leads to a software system with three major layers: the application layer, the mathematical layer, and the numerical algorithm layer. The latter is implemented as a C/C++ library to solve a basic, computational intensive linear problem, such as a linear partial differential equation. The mathematical layer allows the model developer to define his model and to implement high level solution algorithms (e.g. Newton-Raphson scheme, Crank-Nicholson scheme) or choose these algorithms form an algorithm library. The kernels of the model are generic, typically linear, solvers provided through the numerical algorithm layer. Finally, to provide an easy-to-use application environment, a web interface is (semi-automatically) built to edit the XML input file for the modelling code. In the talk, we will discuss the advantages and disadvantages of this concept in more details. We will also present the modelling environment escript which is a prototype implementation toward such a software system in Python (see www.python.org). Key components of escript are the Data class and the PDE class. Objects of the Data class allow generating, holding, accessing, and manipulating data, in such a way that the actual, in the particular context best, representation is transparent to the user. They are also the key to establish connections with external data sources. PDE class objects are describing (linear) partial differential equation objects to be solved by a numerical library. The current implementation of escript has been linked to the finite element code Finley to solve general linear partial differential equations. We will give a few simple examples which will illustrate the usage escript. Moreover, we show the usage of escript together with Finley for the modelling of interacting fault systems and for the simulation of mantel convection.

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.

A rare earth element study of complex zircons from early Archaean Amitsoq gneisses, Godthåbsfjord, south-west Greenland

We present high spatial resolution ion-microprobe rare earth element (REE) data for discrete growth phases of complex polyphase zircons from early Archaean Amitsoq gneisses, outer Godthabsfjord, SW Greenland. In Matsuda diagrams, the two major growth phases, >3.8 Ga cores and ca. 3.65 Ga rims, have steep positive slopes from La to Lu, prominent positive Ce anomalies and negative Eu anomalies that are consistent with growth in a melt. Exceptions to this are non-cathodolurnmescent zircon developed between the cores and rims, sometimes truncating zoning in the cores, and late Archaean prismatic tip overgrowths, both of which exhibit flatter light REE (LREE) patterns and have small or no Eu anomaly, which we interpret as the result of metamorphism and/or small-degree, isolated partial melting. Our data support previous interpretations that the ca. 3.65 Ga zircon phase was generated in a melt, with the >3.8 Ga phase representing either original protolith zircons in a large degree partial melt or inherited zircons in an introduced magma. Regardless which of these two interpretations is correct for these, and similar, rocks in the outer GodthAbsfjord, the 3.65 Ga event will have profoundly affected isotopic systems and obscured beyond recognition any earlier igneous features such as cross-cutting relationships, which may only be assigned a minimum 3.65 Ga age. (C) 2003 Elsevier Science B.V. All rights reserved.

A refined solution to Earth's hidden niobium: implications for evolution of continental crust and mode of core formation

New high-precision niobium (Nb) and tantalum (Ta) concentration data are presented for early Archaean metabasalts, metabasaltic komatiites and their erosion products (mafic metapelites) from SW Greenland and the Acasta gneiss complex, Canada. Individual datasets consistently show sub-chondritic Nb/Ta ratios averaging 15.1+/-11.6. This finding is discussed with regard to two competing models for the solution of the Nb-deficit that characterises the accessible Earth. Firstly, we test whether Nb could have sequestered into the core due to its slightly siderophile (or chalcophile) character under very reducing conditions, as recently proposed from experimental evidence. We demonstrate that troilite inclusions of the Canyon Diablo iron meteorite have Nb and V concentrations in excess of typical chondrites but that the metal phase of the Grant, Toluca and Canyon Diablo iron meteorites do not have significant concentrations of these lithophile elements. We find that if the entire accessible Earth Nb-deficit were explained by Nb in the core, only ca. 17% of the mantle could be depleted and that by 3.7 Ga, continental crust would have already achieved ca. 50% of its present mass. Nb/Ta systematics of late Archaean metabasalts compiled from the literature would further require that by 2.5 Ga, 90% of the present mass of continental crust was already in existence. As an alternative to this explanation, we propose that the average Nb/Ta ratio (15.1+/-11.6) of Earth's oldest mafic rocks is a valid approximation for bulk silicate Earth. This would require that ca. 13% of the terrestrial Nb resided in the Ta-free core. Since the partitioning of Nb between silicate and metal melts depends largely on oxygen fugacity and pressure, this finding could mean that metal/silicate segregation did not occur at the base of a deep magma ocean or that the early mantle was slightly less reducing than generally assumed. A bulk silicate Earth Nb/Ta ratio of 15.1 allows for depletion of up to 40% of the total mantle. This could indicate that in addition to the upper mantle, a portion of the lower mantle is depleted also, or if only the upper mantle were depleted, an additional hidden high Nb/Ta reservoir must exist. Comparison of Nb/Ta systematics between early and late Archaean metabasalts supports the latter idea and indicates deeply subducted high Nb/Ta eclogite slabs could reside in the mantle transition zone or the lower mantle. Accumulation of such slabs appears to have commenced between 2.5 and 2.0 Ga. Regardless of these complexities of terrestrial Nb/Ta systematics, it is shown that the depleted mantle Nb/Th ratio is a very robust proxy for the amount of extracted continental crust, because the temporal evolution of this ratio is dominated by Th-loss to the continents and not Nb-retention in the mantle. We present a new parameterisation of the continental crust volume versus age curve that specifically explores the possibility of lithophile element loss to the core and storage of eclogite slabs in the transition zone. (C) 2003 Elsevier Science B.V. All rights reserved.

40Ar/39Ar constraints on a temporal link between gold mineralization, magmatism, and continental margin transtension in the Jiaodong Gold Province, eastern China

The Jiaodong gold province is the largest gold repository in China. Both mineralization and granitoid hosts are spatially related to the crustal-scale Tan-Lu strike-slip fault system, which developed along the Mesozoic continental margin in eastern China. A series of Ar-40/Ar-39 laser incremental heating analyses of hydrothermal sericite/muscovite from three major gold deposits (Jiaojia, Xincheng, and Wangershan) and igneous biotite from the granodiorite hosts were performed to establish a possible temporal link between gold mineralization, magmatism, and movement along the Tan-Lu fault zone. Magmatic biotite crystals yield well-defined and concordant plateau ages between 124.5+/-0.4 Ma and 124.0+/-0.4 Ma (2sigma), whereas sericite and muscovite samples (a total of 30 single separates) give reproducible plateau ages ranging from 121.0+/-0.4 Ma to 119.2+/-0.2 Ma (2sigma). An integration of our Ar-40/Ar-39 results with age data from other major gold deposits in Jiaodong demonstrates that widespread gold mineralization occurred contemporaneously during a 2-3-m.yr. period. Most gold deposits show intimate spatial associations with abundant mafic to intermediate dikes. The mafic dikes have K-Ar ages of 123.5-119.6 Ma, in excellent agreement with those of the gold deposits. These newly obtained Ar-40/Ar-39 ages, in combination with other independent geological and geochronological data on granodioritic intrusions (130-126 Ma), volcanic rocks (1243.6-114.7 Ma), and deformed rocks within strike-slip faults (132-120 Ma) in Jiaodong or adjacent areas, also support the idea that gold mineralization postdated the granodioritic magmatism but was contemporaneous with mafic magmatism and volcanism, all controlled by the transtensional motion along the Tan-Lu fault in the Early Cretaceous.

A redescription of Atopacanthus dentatus Hussakof and Bryant, 1918 (Acanthodii, Ischnacanthidae)

The holotype specimen of the ischnacanthid acanthodian Alopacanthus dentatus comprises most of the posterior half of a dentigerous jaw bone, rather than a fragment from the middle of the jaw as was previously believed. A new diagnosis and revised description for the taxon are based on the holotype from the Rhinestreet Shale (Frasnian) and other specimens front the North Evans and Genundewa limestones (early Frasnian) of the Genesee Formation (late Givetian-early Frasnian); all of: these dentigerous jaw bones were collected near Hamburg, New York state, U.S.A. An emended diagnosis for Atopacanthus clarifies the differences between it and other ischnacanthid genera.