Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

Role of compaction versus aggregate disruption on slumping and shrinking of repacked hardsetting seedbeds

Slumping of hardsetting seedbeds upon wetting is likely to determine the shrinking and development of strength on drying. Different processes have been invoked, including aggregate disruption, material relocation, and compaction. To gain a better understanding of the role played by compaction compared with aggregate disruption in seedbed slumping and shrinking, mechanical analysis was combined with previous morphogenetical description. The global structural behavior of repacked seedbeds of a hardsetting sandy loam soil was studied after wetting and again after subsequent drying. Bulk density was measured in 5-mm-depth increments using gamma attenuation, and water content was determined at 10-mm-depth increments. Various wetting conditions were used to simulate a range of climatic and management conditions, including flood irrigation, furrow irrigation of a formed seedbed, drip irrigation, and rainfall. Aggregate coalescence under overburden pressure played the main role in slumping, even though microcracking enhanced coalescence. Most of the slumping occurred at calculated effective stress > 1.1 kPa. Intense aggregate breakdown at the top of seedbeds under fast wetting led to slight slumping because the resulting clogging of the initial interaggregate packing voids was balanced, in part, by the increase in microporosity resulting from aggregate disruption. However, aggregate coalescence induced by overburden pressure developing at the seedbed bottom often resulted in a strong decrease in total porosity. The effect of rainfall kinetic energy on crust bulk density was strong compared with the effect of fast wetting (bulk density increase of about 0.07 Mg m(-3) and 0.03 Mg m(-3), respectively) and could be ascribed to compaction rather than to aggregate breakdown. Shrinking on drying was related to the continuity of the microstructure resulting from wetting rather than to the intensity of slumping. Aggregate breakdown led to more shrinking than did aggregate coalescence.

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

Micromorphological study of slumping in a hardsetting seedbed under various wetting conditions

Slumping of hardsetting seedbeds upon wetting has not been extensively studied despite the likelihood that it determines the physical properties after drying. Slumping results from processes similar to those involved in crusting except that overburden pressure can dominate rather than rainfall kinetic energy. Only a few studies have dealt with the morphological description of slumping. To simulate different climatic and management conditions, repacked seedbeds of a hardsetting sandy-loam soil were subjected to a range of wetting conditions, e.g. capillary rise, immersion, and rainfall simulation. Slumping processes were characterized using qualitative and quantitative micromorphological observations of polished blocks and thin sections from resin-impregnated samples. A morphogenetical framework was proposed to help description of the complex associations of processes which can lead to structural collapse (crusting and slumping) on wetting. Three main stages were considered, i.e. aggregate disruption or abrasion, relocation of the released material, and compaction. In the hardsetting material studied here, structural collapse under slow wetting occurred at the bottom of cores due to aggregate coalescence under overburden pressure. Coalescence required aggregate cohesion being reduced by microcracking; therefore, it differed from the coalescence previously described in unstable silty loam soils where microcracking was not necessary for aggregates to coalesce. Macroporosity decreased most strongly under fast wetting due to physical dispersion and aggregate breakdown. Under simulated rainfall, compaction by raindrops could not be distinguish from aggregate breakdown. The role of overburden pressure and of rainfall kinetic energy remains to be stated; new data are required including measurement of total porosity in the initial, wet, and dry states. (C) 2003 Elsevier B.V. All rights reserved.

The Mesozoic and Cenozoic motion of Adria (central Mediterranean): a review of constraints and limitations

This paper presents kinematic analysis on the motion of Adria, which is the continental mass that bridges Africa and Europe in the central Mediterranean. Palaeomagnetic data show a general coherence between the motion of Adria and Africa since the Late Paleozoic. This mutual motion, for the period from 120 Ma and the present, is verified by comparing inferred palaeolatitudes from relatively stable parts of Adria (Apulia, Gargano, Istria, and the Southern Alps) and the Hyblean Plateau, with latitudinal changes that are calculated from the motion of Africa with respect to hotspots. Additional constraints on the motion of Adria are provided from the Late Paleozoic-Early Mesozoic passive margin of Adria in the Ionian Sea. The seismic structure of the floor of the Ionian Sea resembles the structure of the oceanic crust in marginal back-arc basins, suggesting that it formed as a small ocean basin. Furthermore, the Ionian lithosphere in the Calabrian arc has been subjected to rapid rollback, which commonly occurs only when the subducting slab is made of oceanic lithosphere. This oceanic domain marks the Pennian-Triassic to Jurassic plate boundary between Adria and Africa, suggesting that a small amount of independent motion between Adria and Africa took place at that time. Since the Jurassic, Adria and Africa have shared a relatively coherent motion path. (C) 2004 Lavoisier SAS. All rights reserved.

The Western Alps from the Jurassic to Oligocene: spatio-temporal constraints and evolutionary reconstructions

Despite extensive research in the last 150 years, the regional tectonic reconstruction of the Western Alps has remained controversial. The curved orogenic belt consists of several ribbon-like continental terranes (Sesia/Austroalpine, Internal Crystalline Massifs, Brianconnais), which are separated by two or more ophiolitic sutures (Piemonte, Valais, Antrona?, Lanzo/ Canavese?). High-pressure (HP) metamorphism of each terrane occurred during distinct orogenic episodes: at similar to65 Ma in the Sesia/Austroalpine, at similar to45 Ma in the Piemonte zone and at similar to35 Ma in the Internal Crystalline Massifs. It is suggested that these events reflect individual accretionary episodes, which together with kinematic indicators and the speed and direction of plate motions, provide constraints for the discussed reconstruction model. The model involves a prolonged orogenic history that took place during relative convergence of Europe and Adria (here considered as a promontory of the African plate). The first accretionary event involved the Sesia/Austroalpine terrane. Final closure of the Piemonte Ocean occurred during the Eocene (similar to45 Ma) and involved ultra-high-pressure (UHP) metamorphism of the Piemonte oceanic crust. Incorporation of the Brianconnais terrane in the accretionary wedge occurred thereafter, possibly during or after subduction of the Valais Ocean in the late Eocene (45-35 Ma). This subduction was terminated at ca. 35 Ma, when the Internal Crystalline Massifs (i.e. the assumed internal parts of the Brianconnais terrane) were buried into great depths and underwent HP and UHP metamorphism. (C) 2004 Elsevier B.V. All rights reserved.

Continental extension: From core complexes to rigid block faulting

Extension of overthickened continental crust is commonly characterized by an early core complex stage of extension followed by a later stage of crustal-scale rigid block faulting. These two stages are clearly recognized during the extensional destruction of the Alpine orogen in northeast Corsica, where rigid block faulting overprinting core complex formation eventually led to crustal separation and the formation of a new oceanic backarc basin (the Ligurian Sea). Here we investigate the geodynamic evolution of continental extension by using a novel, fully coupled thermomechanical numerical model of the continental crust. We consider that the dynamic evolution is governed by fault weakening, which is generated by the evolution of the natural-state variables (i.e., pressure, deviatoric stress, temperature, and strain rate) and their associated energy fluxes. Our results show the appearance of a detachment layer that controls the initial separation of the brittle crust on characteristic listric faults, and a core complex formation that is exhuming strongly deformed rocks of the detachment zone and relatively undeformed crustal cores. This process is followed by a transitional period, characterized by an apparent tectonic quiescence, in which deformation is not localized and energy stored in the upper crust is transferred downward and causes self-organized mobilization of the lower crust. Eventually, the entire crust ruptures on major crosscutting faults, shifting the tectonic regime from core complex formation to wholesale rigid block faulting.

Hepatogenous photosensitization associated with liver copper accumulation in buffalos

Four 1.5-year-old, male, Murrah buffalos were maintained during eight months without direct solar exposure during a study of copper toxicosis. Four days after solar exposure, all buffalos presented clinical manifestations consistent with acute photosensitization, including anorexia, apathy, loss of body weight, and generalized cutaneous lesions. Gross lesions were characterized by severe erythema, localized edema, fissures, tissue necrosis, gangrene and crust formation with h serous exudation. Liver copper concentration was evaluated, and cutaneous biopsies were taken when clinical signs were evident. The liver copper concentration before solar exposure was increased in all animals. Histopathologic examination of the skin revealed hepatogenous photosensitization characterized by parakeratotic hyperkeratosis, acantholysis, degeneration of squamous epithelial cells, epidermal necrosis with atrophy of sweat glands, and multifocal superficial and deep dermal edema. These findings suggest that asymptomatic accumulation of copper within the liver might have induced hepatic insufficiency thereby resulting in secondary photosensitization when these buffalos were exposed to sunlight. (C) 2009 Elsevier Ltd. All rights reserved.

The geochemistry of late Archaean microbial carbonate: Implications for ocean chemistry and continental erosion history

Trace element concentrations and combined Sr- and Nd-isotope compositions were determined on stromatolitic carbonates (microbialites) from the 2.52 Ga Campbellrand carbonate platform (South Africa). Shale-normalised rare earth element and yttrium patterns of the ancient samples are similar to those of modern seawater in having positive La and Y anomalies and in being depleted in light rare earth elements. In contrast to modem seawater (and microbialite proxies), the 2.52 Ga samples lack a negative Ce anomaly but possess a positive Eu anomaly. These latter trace element characteristics are interpreted to reflect anoxic deep ocean waters where, unlike today, hydrothermal Fe input was not oxidised, and scavenged and rare earth elements were not coprecipitated with Fe-oxyhydroxides. The persistence of a positive Eu anomaly in relatively shallow Campbellrand platform waters indicates a dramatic reversal from hydrothermally dominated (Archaean) to continental erosion-dominated (Phanerozoic) rare earth element flux ratio. The dominant hydrothermal input is also expressed in the initial Sr- and Nd-isotope ratios. There is collinear variation in Sr-Nd systematics, which range from primitive values (Sr-87/Sr-86 of 0.702386 and epsilon (Nd) of +2.1) to more evolved crustal ratios. Mixing calculations show that the range in trace element ratios (e.g., Y/Ho) and initial isotope ratios is not a result of contamination by trapped sediment, but that the chemical band isotopic variation reflects carbonate deposition in an environment where different water masses mixed. Calculated Nd flux ratios yield a hydrothermal input into the 2.52 Ga oceans one order of magnitude larger than continental input. Such a change in flux ratio most likely required substantially reduced continental inputs, which could, in turn, reflect a plate tectonic causation (e.g., reduced topography or expansion of epicontinental seas). Copyright (C) 2001 Elsevier Science Ltd.

Tungsten isotope evidence from similar to 3.8-Gyr metamorphosed sediments for early meteorite bombardment of the Earth

The 'Late Heavy Bombardment' was a phase in the impact history of the Moon that occurred 3.8-4.0 Gyr ago, when the lunar basins with known dates were formed(1,2). But no record of this event has yet been reported from the few surviving rocks of this age on the Earth. Here we report tungsten isotope anomalies, based on the Hf-182-W-182 system (half-life of 9 Myr), in metamorphosed sedimentary rocks from the 3.7-3.8-Gyr-old Isua greenstone belt of West Greenland and closely related rocks from northern Labrador, Canada. As it is difficult to conceive of a mechanism by which tungsten isotope heterogeneities could have been preserved in the Earth's dynamic crust-mantle environment from a time when short-lived Hf-182 was still present, we conclude that the metamorphosed sediments contain a component derived from meteorites.

New W-isotope evidence for rapid terrestrial accretion and very early core formation

The short-lived Hf-182-W-182-isotope system is an ideal clock to trace core formation and accretion processes of planets. Planetary accretion and metal/silicate fractionation chronologies are calculated relative to the chondritic Hf-182-W-182-isotope evolution. Here, we report new high-precision W-isotope data for the carbonaceous chondrite Allende that are much less radiogenic than previously reported and are in good agreement with published internal Hf-W chronometry of enstatite chondrites. If the W-isotope composition of terrestrial rocks, representing the bulk silicate Earth, is homogeneous and 2.24 epsilon(182W) units more radiogenic than that of the bulk Earth, metal/silicate differentiation of the Earth occurred very early. The new W-isotope data constrain the mean time of terrestrial core formation to 34 million years after the start of solar system accretion. Early terrestrial core formation implies rapid terrestrial accretion, thus permitting formation of the Moon by giant impact while Hf-182 was still alive. This could explain why lunar W-isotopes are more radiogenic than the terrestrial value. Copyright (C) 2002 Elsevier Science Ltd.

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.

Applications of accurate, high-precision Pb isotope ratio measurement by multi-collector ICP-MS

The isotope composition of Ph is difficult to determine accurately due to the lack of a stable normalisation ratio. Double and triple-spike addition techniques provide one solution and presently yield the most accurate measurements. A number of recent studies have claimed that improved accuracy and precision could also be achieved by multi-collector ICP-MS (MC-ICP-MS) Pb-isotope analysis using the addition of Tl of known isotope composition to Pb samples. In this paper, we verify whether the known isotope composition of Tl can be used for correction of mass discrimination of Pb with an extensive dataset for the NIST standard SRM 981, comparison of MC-ICP-MS with TIMS data, and comparison with three isochrons from different geological environments. When all our NIST SRM 981 data are normalised with one constant Tl-205/Tl-203 of 2.38869, the following averages and reproducibilities were obtained: Pb-207/Pb-206=0.91461+/-18; Pb-208/Ph-206 = 2.1674+/-7; and (PbPh)-Pb-206-Ph-204 = 16.941+/-6. These two sigma standard deviations of the mean correspond to 149, 330, and 374 ppm, respectively. Accuracies relative to triple-spike values are 149, 157, and 52 ppm, respectively, and thus well within uncertainties. The largest component of the uncertainties stems from the Ph data alone and is not caused by differential mass discrimination behaviour of Ph and Tl. In routine operation, variation of sample introduction memory and production of isobaric molecular interferences in the spectrometer's collision cell currently appear to be the ultimate limitation to better reproducibility. Comparative study of five different datasets from actual samples (bullets, international rock standards, carbonates, metamorphic minerals, and sulphide minerals) demonstrates that in most cases geological scatter of the sample exceeds the achieved analytical reproducibility. We observe good agreement between TIMS and MC-ICP-MS data for international rock standards but find that such comparison does not constitute the ultimate. test for the validity of the MC-ICP-MS technique. Two attempted isochrons resulted in geological scatter (in one case small) in excess of analytical reproducibility. However, in one case (leached Great Dyke sulphides) we obtained a true isochron (MSWD = 0.63) age of 2578.3 +/- 0.9 Ma, which is identical to and more precise than a recently published U-Pb zircon age (2579 3 Ma) for a Great Dyke websterite [Earth Planet. Sci. Lett. 180 (2000) 1-12]. Reproducibility of this age by means of an isochron we regard as a robust test of accuracy over a wide dynamic range. We show that reliable and accurate Pb-isotope data can be obtained by careful operation of second-generation MC-ICP magnetic sector mass spectrometers. (C) 2002 Elsevier Science B.V. All rights reserved.

Fluid-mobile trace element constraints on the role of slab melting and implications for Archaean crustal growth models

We use published and new trace element data to identify element ratios which discriminate between arc magmas from the supra-subduction zone mantle wedge and those formed by direct melting of subducted crust (i.e. adakites). The clearest distinction is obtained with those element ratios which are strongly fractionated during refertilisation of the depleted mantle wedge, ultimately reflecting slab dehydration. Hence, adakites have significantly lower Pb/Nd and B/Be but higher Nb/Ta than typical arc magmas and continental crust as a whole. Although Li and Be are also overenriched in continental crust, behaviour of Li/Yb and Be/Nd is more complex and these ratios do not provide unique signatures of slab melting. Archaean tonalite-trondhjemite-granodiorites (TTGs) strongly resemble ordinary mantle wedge-derived arc magmas in terms of fluid-mobile trace element content, implying that they-did not form by slab melting but that they originated from mantle which was hydrated and enriched in elements lost from slabs during prograde dehydration. We suggest that Archaean TTGs formed by extensive fractional crystallisation from a mafic precursor. It is widely claimed that the time between the creation and subduction of oceanic lithosphere was significantly shorter in the Archaean (i.e. 20 Ma) than it is today. This difference was seen as an attractive explanation for the presumed preponderance of adakitic magmas during the first half of Earth's history. However, when we consider the effects of a higher potential mantle temperature on the thickness of oceanic crust, it follows that the mean age of oceanic lithosphere has remained virtually constant. Formation of adakites has therefore always depended on local plate geometry and not on potential mantle temperature.