On the overabundance of light rare earth elements in terrestrial zircons and its implication for Earth's earliest magmatic differentiation

We present whole-rock and zircon rare earth element (REE) data from two early Archaean gneisses (3.81 Ga and 3.64 Ga) from the Itsaq gneiss complex, south-west Greenland. Both gneisses represent extremely rare examples of unaltered, fresh and relatively undeformed igneous rocks of such antiquity. Cathodoluminescence imaging of their zircons indicates a single crystallisation episode with no evidence for either later metamorphic and/or anatectic reworking or inheritance of earlier grains. Uniform, single-population U/Pb age data confirm the structural simplicity of these zircons. One sample, a 3.64 Ga granodioritic gneiss from the Gothabsfjord, yields a chondrite-normalised REE pattern with a positive slope from La to Lu as well as substantial positive Ce and slight negative Eu anomalies, features generally considered to be typical of igneous zircon. In contrast, the second sample, a 3.81 Ga tonalite from south of the Isua Greenstone Belt, has variable but generally much higher light REE abundances, with similar middle to heavy REE. Calculation of zircon/melt distribution coefficients (D-REE(zircon/melt)) from each sample yields markedly different values for the trivalent REE (i.e. Ce and Eu omitted) and simple application of one set of D-REE(zircon/melt) to model the melt composition for the other sample yields concentrations that are in error by up to two orders of magnitude for the light REE (La-Nd). The observed light REE overabundance in the 3.81 Ga tonalite is a commonly observed feature in terrestrial zircons for which a number of explanations ranging from lattice strain to disequilibrium crystallisation have been proposed and are further investigated herein. Regardless of the cause of light REE overabundance, our study shows that simple application of zircon/melt distribution coefficients is not an unambiguous method for ascertaining original melt composition. In this context, recent studies that use REE data to claim that > 4.3 Ga Hadean detrital zircons originally crystallised from an evolved magma, in turn suggesting the operation of geological processes in the early Earth analogous to those of the present day (e.g. subduction and melting of hydrated oceanic crust), must be regarded with caution. Indeed, comparison of terrestrial Hadean and > 3.9 Ga lunar highland zircons shows remarkable similarities in the light REE, even though subduction processes that have been used to explain the terrestrial zircons have never operated on the Moon. (C) 2002 Elsevier Science B.V. All rights reserved.

Upper Devonian biostratigraphy of the Iberian Pyrite Belt, southwest Spain Part One: Miospores

The Iberian Pyrite Belt (IPB), which forms part of the Variscan orogenic massif, is renowned for the magnitude and extent of its massive sulfide mineralization. The stratigraphic record of the IPB consists of Upper Palaeozoic sedimentary and igneous rocks. In ascending order, these comprise the thick Phyllite-Quartzite Group attributed to the Middle and Upper Devonian and characterized by shales and quartzites with conglomeratic and carbonate intercalations towards the top; the appreciably thinner Volcano-Sedimentary Complex, a heterogeneous uppermost Devonian-Mississippian unit embodying diverse volcanic, subvolcanic, and sedimentary rocks that host the massive sulfide deposits; and the shaly and sandy, turbiditic Culm Group (Carboniferous). This entire succession was folded and faulted during the Asturian phase of the Variscan Orogeny that gave rise to a thin-skinned type structure. The present study constitutes a detailed blostratigraphic investigation of palynologically productive samples representative of the Phyllite-Quartzite Group and the basal (anoxic) portion of the Volcano-Sedimentary Complex. These were collected from surface and mine exposures variously located in the Spanish part of the IPB; out of 282 samples processed, 117 proved to be productive palynologically. The aim of this project is to provide comprehensive palynostratigraphic data applicable to precise dating and correlation of the IPB's stratigraphic succession (i.e., of the two sampled lithostratigraphic units), which has hitherto been investigated biostratigraphically on a relatively localized basis. The results are incorporated in two successive parts. The first of these, i. e., the present paper, focuses on the systematic analysis of the terrestrial (miospore) component of the palynological assemblages. The second part, devoted to the marine, organic-walled microphytoplankton (acritarchs and prasinophytes), will evaluate the stratigraphic significance of the IPB palynofloras and their application to elucidating the geological history of the region. In the systematic-descriptive section, which occupies the bulk of this paper, 55 species of trilete miospores are described and are allocated among 34 genera, two of which (Cristicavatispora and Epigruspora) are newly instituted herein. The majority of the species are either positively identifiable or closely affiliable with previously named species. The nine newly established species are as follows: Camptozonotriletes confertus, Indotriradites diversispinosus, Cristicavatispora dispersa (type species), Epigruspora regularis (type species), Ancyrospora? implicata, Endosporites tuberosus, Rugospora explicata, Spelaeotriletes plicatus, and Teichertospora iberica.

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.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

A sensitive vorticity gauge using rotated porphyroblasts, and its application to rocks adjacent to the Alpine Fault, New Zealand

Variable aspect ratio porphyroblasts deformed in non-coaxial flow. and internally containing rotated relicts of an external foliation, can be used to characterise plane strain flow regimes. The distribution obtained by plotting the orientation of the long axis of such grains, classified by aspect ratio, against the orientation of the internal foliation is potentially a sensitive gauge of both the bulk shear strain (as previously suggested) and kinematic vorticity number. We illustrate the method using rotated biotite porphyroblasts in the Alpine Schist: a sequence of mid-crustal rocks that have been ramped to the surface along the Alpine Fault. a major transpressional plate boundary. Results indicate that, at distances greater than or equal to similar to1 km from the fault, the rocks have undergone a combination of irrotational fattening and dextral-oblique, normal-sense shear, with a bulk shear strain of similar to0.6 and kinematic vorticity number of similar to0.2. The vorticity analysis is compatible with estimates of strongly oblate bulk strain of similar to 75% maximum shortening. Dextral-reverse transpressional flow characterises higher strain S-tectonite mylonite within similar to1 km of the Alpine Fault. These relationships provide insight into the kinematics of flow and distribution of strain in the hangingwall of the Alpine Fault and place constraints on numerical mechanical models for the exhumation of these mid-crustal rocks. (C) 2001 Elsevier Science Ltd. All rights reserved.

Petrology and geochemistry of the Bandas del Sur Formation, Las Canadas Edifice, Tenerife (Canary Islands)

The Bandas del Sur Formation preserves a Quaternary extra-caldera record of central phonolitic explosive volcanism of the Las Canadas volcano at Tenerife. Volcanic rocks are bimodal in composition, being predominantly phonolitic pyroclastic deposits, several eruptions of which resulted in summit caldera collapse, alkali basaltic lavas erupted from many fissures around the flanks. For the pyroclastic deposits, there is a broad range of pumice glass compositions from phonotephrite to phonolite. The phonolite pyroclastic deposits are also characterized by a diverse, 7-8-phase phenocryst assemblage (alkali feldspar + biotite + sodian diopside + titanomagnetite + ilmenite + nosean-hauyne + titanite + apatite) with alkali feldspar dominant, in contrast to interbedded phonolite lavas that typically have lower phenocryst contents and lack hydrous phases. Petrological and geochemical data are consistent with fractional crystallization (involving the observed phenocryst assemblages) as the dominant process in the development of phonolite magmas. New stratigraphically constrained data indicate that petrological and geochemical differences exist between pyroclastic deposits of the last two explosive cycles of phonolitic volcanism. Cycle 2 (0.85-0.57 Ma) pyroclastic fall deposits commonly show a cryptic compositional zonation indicating that several eruptions tapped chemically, and probably thermally stratified magma systems. Evidence for magma mixing is most widespread in the pyroclastic deposits of Cycle 3 (0.37-0.17 Ma), which includes the presence of reversely and normally zoned phenocrysts, quenched mafic glass blebs in pumice, banded pumice, and bimodal to polymodal phenocryst compositional populations. Syn-eruptive mixing events involved mostly phonolite and tephriphonolite magmas, whereas a pre-eruptive mixing event involving basaltic magma is recorded in several banded pumice-bearing ignimbrites of Cycle 3. The periodic addition and mixing of basaltic magma ultimately may have triggered several eruptions. Recharge and underplating by basaltic magma is interpreted to have elevated sulphur contents (occurring as an exsolved gas phase) in the capping phonolitic magma reservoir. This promoted nosean-hauyne crystallization over nepheline, elevated SO3 contents in apatite, and possibly resulted in large, climatologically important SO2 emissions.

Mineral chemistry, whole-rock compositions, and petrogenesis of leg 176 gabbros: Data and discussion

We report mineral chemistry, whole-rock major element compositions, and trace element analyses on Hole 735B samples drilled and selected during Leg 176. We discuss these data, together with Leg 176 shipboard data and Leg 118 sample data from the literature, in terms of primary igneous petrogenesis. Despite mineral compositional variation in a given sample, major constituent minerals in Hole 735B gabbroic rocks display good chemical equilibrium as shown by significant correlations among Mg# (= Mg/[Mg+Fe2+]) of olivine, clinopyroxene, and orthopyroxene and An (=Ca/[Ca+Na]) of plagioclase. This indicates that the mineral assemblages olivine + plagioclase in troctolite, plagioclase + clinopyroxene in gabbro, plagioclases + clinopyroxene + olivine in olivine gabbro, and plagioclase + clinopyroxene + olivine + orthopyroxene in gabbronorite, and so on, have all coprecipitated from their respective parental melts. Fe-Ti oxides (ilmenite and titanomagnetite), which are ubiquitous in most of these rocks, are not in chemical equilibrium with olivine, clinopyroxene, and plagioclase, but precipitated later at lower temperatures. Disseminated oxides in some samples may have precipitated from trapped Fe-Ti–rich melts. Oxides that concentrate along shear bands/zones may mark zones of melt coalescence/transport expelled from the cumulate sequence as a result of compaction or filter pressing. Bulk Hole 735B is of cumulate composition. The most primitive olivine, with Fo = 0.842, in Hole 735B suggests that the most primitive melt parental to Hole 735B lithologies must have Mg# ≤ 0.637, which is significantly less than Mg# = 0.714 of bulk Hole 735B.

Kinematics of oblique collision and ramping inferred from microstructures and strain in middle crustal rocks, central Southern Alps, New Zealand

The hanging wall of the Alpine Fault near Franz Josef Glacier has been exhumed during the past similar to2-3 m.y. providing a sample of the ductilely deformed middle crust of a modem obliquely convergent orogen. Presently exposed rocks of the Pacific Plate are inferred to have undergone several phases of ductile deformation as they moved westward above a mid-crustal detachment. Initially they were transpressed across the outboard part of the orogen, resulting in oblate fabrics with a down-dip stretch. Later, they encountered the Alpine Fault, experiencing an oblique-slip backshearing on vertical planes. This escalator-like deformation tilted and thinned the incoming crust onto that crustal-scale oblique ramp. This style of hanging wall deformation may affect only the most rapidly uplifting, central part of the Southern Alps because of the low flexural rigidity of the crust in that region and its displacement over a relatively sharp ramp-angle at depth. A 3D transpressive flow affected mylonites locally near the fault, but their shear direction remained parallel to plate motion, ruling out ductile 'extrusion' as an important process in this orogen. Outside the mylonite zone, late Cenozoic shortening is inferred to be modest (30-40%), as measured from deformation of younger biotite grains. Oblique collision is dominated by translation on the Alpine Fault, and rocks migrate rapidly through the deforming zone, preventing the accumulation of large finite strains. Transpression may play a minor role in oblique collision. (C) 2001 Elsevier Science Ltd. All rights reserved.

Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.