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.

Mantle convection modeling with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling

The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.

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.

The effect of energy feedbacks on continental strength

The classical strength profile of continents(1,2) is derived from a quasi-static view of their rheological response to stress-one that does not consider dynamic interactions between brittle and ductile layers. Such interactions result in complexities of failure in the brittle-ductile transition and the need to couple energy to understand strain localization. Here we investigate continental deformation by solving the fully coupled energy, momentum and continuum equations. We show that this approach produces unexpected feedback processes, leading to a significantly weaker dynamic strength evolution. In our model, stress localization focused on the brittle-ductile transition leads to the spontaneous development of mid-crustal detachment faults immediately above the strongest crustal layer. We also find that an additional decoupling layer forms between the lower crust and mantle. Our results explain the development of decoupling layers that are observed to accommodate hundreds of kilometres of horizontal motions during continental deformation.

Magma Flow Instabilities in a Volcanic Conduit: Implications for Long-Period Seismicity

Silicic volcanic eruptions are typically accompanied by repetitive Long-Period (LP) seismicity that originates from a small region of the upper conduit. These signals have the capability to advance eruption prediction, since they commonly precede a change in the eruption vigour. Shear bands forming along the conduit wall, where the shear stresses are highest, have been linked to providing the seismic trigger. However, existing computational models are unable to generate shear bands at the depths where the LP signals originate using simple magma strength models. Presented here is a model in which the magma strength is determined from a constitutive relationship dependent upon crystallinity and pressure. This results in a depth-dependent magma strength, analogous to planetary lithospheres. Hence, in shallow highly-crystalline regions a macroscopically discontinuous brittle type of deformation will prevail, whilst in deeper crystal-poor regions there will be a macroscopically continuous plastic deformation mechanism. This will result in a depth where the brittle-ductile transition occurs, and here shear bands disconnected from the free-surface may develop. We utilize the Finite Element Method and use axi-symmetric coordinates to model magma flow as a viscoplastic material, simulating quasi-static shear bands along the walls of a volcanic conduit. Model results constrained to the Soufrière Hills Volcano, Montserrat, show the generation of two types of shear bands: upper-conduit shear bands that form between the free-surface to a few 100 metres below it and discrete shear bands that form at the depths where LP seismicity is measured to occur corresponding to the brittle-ductile transition and the plastic shear region. It is beyond the limitation of the model to simulate a seismic event, although the modelled viscosity within the discrete shear bands suggests a failure and healing cycle time that supports the observed LP seismicity repeat times. However, due to the paucity of data and large parameter space available these results can only be considered to be qualitative rather than quantitative at this stage.

The role of mobile belts for the longevity of deep cratonic lithosphere: The crumple zone model

Many Archean cratons are surrounded by Proterozoic mobile belts that have experienced episodes of tectonic re-activation over their lifetimes. This suggests that mobile belt lithosphere may be associated with long lived, inherited weakness. It is proposed that the proximity of this weakness can increase the longevity of deep Archean lithosphere by buffering Archean cratons from mantle derived stresses. The physical plausibility of this idea is explored through numerical simulations of mantle convection that include continents and allow for material rheologies that model the combined brittle and ductile behavior of the lithosphere. Within the simulations, the longevity of deep cratonic lithosphere does increase if it is buffered by mobile belts that can fail at relatively low stress levels.

Surface characterization of 6H-SiC (0001) substrates in indentation and abrasive machining

Surface characterization of 6H-SiC (0001) substrates in indentation and abrasive machining was carried out to investigate microfracture, residual damage, and surface roughness associated with material removal and surface generation. Brittle versus plastic deformation was studied using Vickers indention and nano-indentation. To characterize the abrasive machining response, the 6H-SiC (0001) substrates were ground using diamond wheels with grit sizes of 25, 15 and 7 mum, and then polished with diamond suspensions of 3 and 0.05 mum. It is found that in indentation, there was a scale effect for brittle versus plastic deformation in 6H-SiC substrates. Also, in grinding, the scales of fracture and surface roughness of the substrates decreased with a decrease in diamond grit size. However, in polishing, a reduction in grit size of diamond suspensions gave no significant improvement in surface roughness. Furthermore, the results showed that fracture-free 6H-SiC (0001) surfaces were generated in polishing with the existence of the residual crystal defects, which were associated with the origin of defects in single crystal growth. (C) 2003 Elsevier Ltd. 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.

Ductile fabrics in the zone of active oblique convergence near the Alpine Fault, New Zealand: identifying the neotectonic overprint

The mid-crustal Alpine Schist in central Southern Alps, New Zealand has been exhumed during the past similar to3 m.y. on the hanging wall of the oblique-slip Alpine Fault. These rocks underwent ductile deformation during their passage through the similar to 150-km-wide Pacific-Australia plate boundary zone. Likely to be Cretaceous in age, peak metamorphism predates the largely Pliocene and younger oblique convergence that continues to uplift the Southern Alps today. Late Cenozoic ductile deformation constructively reinforced a pre-existing fabric that was well oriented to accommodate a dextral-transpressive overprint. Quartz microstructures below a recently exhumed brittle-ductile transition zone reflect a late Cenozoic increment of ductile strain that was distributed across deeper levels of the Pacific Plate. Deformation was transpressive, including a dextral-normal shear component that bends and rotates a delaminated panel of Pacific Plate crust onto the oblique footwall ramp of the Alpine Fault. Progressive ductile shear in mylonites at the base of the Pacific Plate overprints earlier fabrics in a dextral-reverse sense, a deformation that accompanies translation of the schists up the Alpine Fault. Ductile shear along that structure affects not only the 12-km-thick section of Alpine mylonites, but is distributed across several kilometres of overlying nonmylonitic rocks. (C) 2001 Elsevier Science Ltd. All rights reserved.

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 load-unload response ratio and critical sensitivity in the lattice solid model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

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.

The plasticisation effect of glycerol and water on the gelatinisation of wheat starch

In the extrusion manufacture of starch-based thermoplastics, such as biodegradable packaging materials, glycerol is an effective additive as a plasticiser, that is, to diminish the brittle nature of the product and provide the desired extent of flexibility. However, the addition of glycerol may also affect the gelatinisation behaviour of the starch-water mixture, and hence the required processing conditions for producing a homogeneously gelatinised starch-based material. The effect of glycerol on the gelatinisation of wheat starch was studied using differential scanning calorimetry (DSC). Mixtures of starch, water and glycerol were investigated with a water content ranging from 12 - 40% and a glycerol concentration up to 75%. Dependent on composition, the enthalpy of gelatinisation ranged from 1.7 - 12.6 J/g (on a dry starch basis), while the onset and peak temperatures varied from 54 to 86 degreesC and 60 to 90 degreesC, respectively. As expected, water acted as a plasticiser in that the onset temperature for gelatinisation (TO) decreased with increasing moisture content. Glycerol, however, increased To. It is shown that the T-0 of starch-glycerol-water mixtures may be predicted on the basis of the effective moisture content of the starch fraction of these mixtures resulting from the relative speed of moisture absorption by glycerol and starch, respectively. Moisture sorption kinetics of wheat starch and glycerol in 100% relative humidity were determined and used to predict the preferential water absorption by glycerol in starch-glycerol-water mixtures and hence the resulting T-0 of the system.