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.

Inheritance of early Archaean Pb-isotope variability from long-lived Hadean protocrust

Comparison of initial Pb-isotope signatures of several early Archaean (3.65-3.82 Ga) lithologies (orthogneisses and metasediments) and minerals (feldspar and galena) documents the existence of substantial isotopic heterogeneity in the early Archaean, particularly in the Pb-207/Pb-204 ratio. The magnitude of isotopic variability at 3.82-3.65 Ga requires source separation between 4.3 and 4.1 Ga, depending on the extent of U/Pb fractionation possible in the early Earth. The isotopic heterogeneity could reflect the coexistence of enriched and depleted mantle domains or the separation of a terrestrial protocrust with a U-238/Pb-204 (mu) that was ca. 20-30% higher than coeval mantle. We prefer this latter explanation because the high-p signature is most evident in metasediments (that formed at the Earth's surface). This interpretation is strengthened by the fact that no straightforward mantle model can be constructed for these high-mu lithologies without violating bulk silicate Earth constraints. The Pb-isotope evidence for a long-lived protocrust complements similar Hf-isotope data from the Earth's oldest zircons, which also require an origin from an enriched (low Lu/Hf) environment. A model is developed in which greater than or equal to3.8-Ga tonalite and monzodiorite gneiss precursors (for one of which we provide zircon U-Pb data) are not mantle-derived but formed by remelting or differentiation of ancient (ca. 4.3 Ga) basaltic crust which had evolved with a higher U/Pb ratio than coeval mantle in the absence of the subduction process. With the initiation of terrestrial subduction at, we propose, ca. 3.75 Ga, most of the greater than or equal to3.8-Ga basaltic shell (and its differentiation products) was recycled into the mantle, because of the lack of a stabilising mantle lithosphere. We argue that the key event for preservation of all greater than or equal to3.8-Ga terrestrial crust was the intrusion of voluminous granitoids immediately after establishment of global subduction because of complementary creation of a lithospheric keel. Furthermore, we argue that preservation of !3.8-Ga material (in situ rocks and zircons) globally is restricted to cratons with a high U/Pb source character (North Atlantic, Slave, Zimbabwe, Yilgarn, and Wyoming), and that the Pb-isotope systematics of these provinces are ultimately explained by reworking of material that was derived from ca. 4.3 Ga (i.e. Hadean) basaltic crust.

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.

40Ar/39Ar dating of Mn-iron ore of the Sapecado Mine, Qadrilatero Ferrifero, Brazil

The SW region of Amazonian craton presents policyclic evolution between 1.80-1.00 Ga and is comprised of the Rio Negro-Juruena, Rondoniana and Sunsas Provinces. The evolution of this region has being characterized by four orogens: Alto Jauru (1.79-1.74 Ga), Cachoeirinha (1.58-1.52 Ga), Suíte Santa Helena (1.45-1.42 Ga) e Sunsas/Aguapeí (1.0-0.9 Ga). The Alto Jauru orogen consists of TTG gneissic associations, greenstone sequences and intrusive granitoids origined in volcanic arc setting. Eight 40Ar/39Ar step-heating analyses were carried out in minerals (biotiteand hornblende) to investigate the thermal history and crustal evolution of this region. From the Alto Jauru orogen was sampled the gneiss banded and two biotite grains provide large dispersion of apparent ages, suggesting heterogenity in reservoir of the argon. Apparent age diagram yielded integrated ages of 1472 ± 6 Ma, interpreted as minimum ages of regional cooling episode. Three analyses of hornblende present ages varing from 1310 to 1400 Ma, possibly because smaller grain size become more susceptible to argon loss. 40Ar/39Ar step-heating methodology applied on biotite of pyroclastic tuff (U-Pb age about 1758 ± 7 Ma) presented integrated age of 1507 ± 7 Ma. The results found for this terrane demonstrated a geochronological correlation with metamorphic process linked Cachoeirinha orogen. Biotite and hornblende grains separates from granite and a tonalite origined during Cachoeirinha orogen were analyzed and the apparent age diagrams indicated well-defined plateau ages of 1520-1540 Ma. Biotite grains from a granitic sample were analized, and integrated ages about 1526 ± 2 Ma were obtained due argon loss in the initial steps. Thermochronologic history of SW region Amazonian craton is coherent with regional policyclic events and 40Ar/39Ar ages here presented probably correspond to regional cooling period of Cachoeirinha orogen.

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.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

Microstructural evolution of the Serido Belt, NE-Brazil: the effect of two tectonic events on development of c-axis preferred orientation in quartz

The polyphase evolution of the Serido Belt (NE-Brazil) includes D, crust formation at 2.3-2.1 Ga, D-2 thrust tectonics at 1.9 Ga and crustal reworking by D-3 strike-slip shear zones at 600 Ma. Microstructural investigations within mylonites associated with D-2 and D-3 events were used to constrain the tectono-thermal evolution of the belt. D-2 shear zones commenced at deeper crustal levels and high amphibolite facies conditions (600-650 degreesC) through grain boundary migration, subgrain rotation and operation of quartz Q-prism slip. Continued shearing and exhumation of the terrain forced the re-equilibration of high-T fabrics and the switching of slip systems from (c)-prism to positive and negative (a)-rhombs. During D-3, enhancement of ductility by dissipation of heat that came from syn-D-3 granites developed wide belts of amphibolite facies mylonites. Continued shearing, uplift and cooling of the region induced D-3 shear zones to act in ductile-brittle regimes, marked by fracturing and development of thinner belts of greenschist facies mylonites. During this event, switching from (a)-prism to a basal slip indicates a thermal path from 600 to 350 degreesC. Therefore, microstructures and quartz c-axis fabrics in polydeformed rocks from the Serido Belt preserve the record of two major events, which includes contrasting deformation mechanisms and thermal paths. (C) 2003 Elsevier Ltd. All rights reserved.

Changes in organic matter production and accumulation as a mechanism for isotopic evolution in the Mesoproterozoic ocean

Mesoproterozoic marine successions worldwide record a shift in average delta(13)C values from 0 to +3.5parts per thousand, with the latter value evident in successions younger than 1250 Ma. New carbon isotope data from the similar to 1300 to 1270 Ma Dismal Lakes Group, Arctic Canada, provide further insight into this fundamental transition. Data reveal that the shift to higher VC values was gradual and marked by occasional excursions to values less than 0 parts per thousand. When compared to records from older and younger marine successions, it is evident that the difference between isotopic minima and maxima increased with time, indicating that the marine system evolved to become isotopically more variable. We interpret these patterns to record an increase in the crustal inventory of organic carbon, reflecting eukaryotic diversification and a change in the locus of organic carbon burial to include anoxic deep marine sites where preservation potential was high. We speculate that the release of O-2 to Earth's surface environments associated with increased organic carbon storage induced irreversible changes in the Mesoproterozoic biosphere, presaging the more extreme environmental and evolutionary developments of the Neoproterozoic.

Geochronology of weathering and landscape evolution, Dugald River valley, NW Queensland, Australia

40Ar/39Ar laser incremental heating analyses of individual grains of supergene jarosite, alunite, and cryptomelane from weathering profiles in the Dugald River area, Queensland, Australia, show a strong positive correlation between a sample’s age and its elevation. We analyzed 125 grains extracted from 35 hand specimens collected from weathering profiles at 11 sites located at 3 distinct elevations. The highest elevation profile hosts the oldest supergene minerals, whereas progressively younger samples occur at lower positions in the landscape. The highest elevation sampling sites (three sites), located on top of an elongated mesa (255 to 275 m elevation), yield ages in the 16 to 12 Ma range. Samples from an intermediate elevation site (225 to 230 m elevation) yield ages in the 6 to 4 Ma range. Samples collected at the lowest elevation sites (200 to 220 m elevation) yield ages in the 2.2 to 0.8 Ma interval. Grains of supergene alunite, jarosite, and cryptomelane analyzed from individual single hand specimens yield reproducible results, confirming the suitability of these minerals to 40Ar/39Ar geochronology. Multiple samples collected from the same site also yield reproducible results, indicating that the ages measured are true precipitation ages for the samples analyzed. Different sites, up to 3 km apart, sampled from weathering profiles at the same elevation again yield reproducible results. The consistency of results confirms that 40Ar/39Ar geochronology of supergene jarosite, alunite, and cryptomelane yields ages of formation of weathering profiles, providing a reliable numerical basis for differentiating and correlating these profiles. The age versus elevation relationship obtained suggest that the stepped landscapes in the Dugald River area record a progressive downward migration of a relatively flat weathering front. The steps in the landscape result from differential erosion of previously weathered bedrock displaying different susceptibility to weathering and contrasting resistance to erosion. Combined, the age versus elevation relationships measured yield a weathering rate of 3.8 m. Myr−1 (for the past 15 Ma) if a descending subhorizontal weathering front is assumed. The results also permit the calculation of the erosion rate of the more easily weathered and eroded lithologies, assuming an initially flat landscape as proposed in models of episodic landscape development. The average erosion rate for the past 15 Ma is 3.3 m. Myr−1, consistent with erosion rates obtained by cosmogenic isotope studies in the region.

Thermal effects in the evolution of initially layered mantle material

A simplified model for anisotropic mantle convection based on a novel class of rheologies, originally developed for folding instabilities in multilayered rock (MUHLHAUS et al., 2002), is extended ¨ through the introduction of a thermal anisotropy dependent on the local layering. To examine the effect of the thermal anisotropy on the evolution of mantle material, a parallel implementation of this model was undertaken using the Escript modelling toolkit and the Finley finite-element computational kernel (DAVIES et al., 2004). For the cases studied, there appears too little if any effect. For comparative purposes, the effects of anisotropic shear viscosity and the introduced thermal anisotropy are also presented. These results contribute to the characterization of viscous anisotropic mantle convection subject to variation in thermal conductivities and shear viscosities.