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.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Determination of mechanical properties of PECVD silicon nitride thin films for tunable MEMS Fabry-Pérot optical filters

This paper reports an investigation on techniques for determining elastic modulus and intrinsic stress gradient in plasma-enhanced chemical vapor deposition (PECVD) silicon nitride thin films. The elastic property of the silicon nitride thin films was determined using the nanoindentation method on silicon nitride/silicon bilayer systems. A simple empirical formula was developed to deconvolute the film elastic modulus. The intrinsic stress gradient in the films was determined by using micrometric cantilever beams, cross-membrane structures and mechanical simulation. The deflections of the silicon nitride thin film cantilever beams and cross-membranes caused by in-thickness stress gradients were measured using optical interference microscopy. Finite-element beam models were built to compute the deflection induced by the stress gradient. Matching the deflection computed under a given gradient with that measured experimentally on fabricated samples allows the stress gradient of the PECVD silicon nitride thin films introduced from the fabrication process to be evaluated.

Characterization of mechanical properties of silicon nitride thin films for MEMS devices by nanoindentation

An experimental investigation of mechanical properties of thin films using nanoindentation was reported. Silicon nitride thin films with different thicknesses were deposited using plasma enhanced chemical vapor deposition (PECVD) on Si substrate. Nanoindentation was used to measure their elastic modulus and hardness. The results indicated that for a film/substrate bilayer system, the measured mechanical properties are significantly affected by the substrate properties. Empirical formulas were proposed for deconvoluting the film properties from the measured bilayer properties.

A continuum model for periodic two-dimensional block structures

A continuum model for regular block structures is derived by replacing the difference quotients of the discrete equations by corresponding differential quotients. The homogenization procedure leads to an anisotropic Cosserat Continuum. For elastic block interactions the dispersion relations of the discrete and the continuous models are derived and compared. Yield criteria for block tilting and sliding are formulated. An extension of the theory for large deformation is proposed. (C) 1997 by John Wiley & Sons, Ltd.

Analysis and implementation of a new constant acceleration subcycling algorithm

An algorithm for explicit integration of structural dynamics problems with multiple time steps is proposed that averages accelerations to obtain subcycle states at a nodal interface between regions integrated with different time steps. With integer time step ratios, the resulting subcycle updates at the interface sum to give the same effect as a central difference update over a major cycle. The algorithm is shown to have good accuracy, and stability properties in linear elastic analysis similar to those of constant velocity subcycling algorithms. The implementation of a generalised form of the algorithm with non-integer time step ratios is presented. (C) 1997 by John Wiley & Sons, Ltd.

Finite element analysis of mechanical properties in discontinuously reinforced metal matrix composites with ultrafine microstructure

This investigation focused on the finite element analyses of elastic and plastic properties of aluminium/alumina composite materials with ultrafine microstructure. The commonly used unit cell model was used to predict the elastic properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the plastic properties of the composites and the associated strengthening mechanism within the metal matrix material were investigated. The grain size of the matrix material was found to be an important factor influencing the mechanical properties of the composites studied. (C) 1997 Elsevier Science S.A.

Proton transport model in the ionosphere .1. Multistream approach of the transport equations

The suprathermal particles, electrons and protons, coming from the magnetosphere and precipitating into the high-latitude atmosphere are an energy source of the Earth's ionosphere. They interact with ambient thermal gas through inelastic and elastic collisions. The physical quantities perturbed by these precipitations, such as the heating rate, the electron production rate, or the emission intensities, can be provided in solving the kinetic stationary Boltzmann equation. This equation yields particle fluxes as a function of altitude, energy, and pitch angle. While this equation has been solved through different ways for the electron transport and fully tested, the proton transport is more complicated. Because of charge-changing reactions, the latter is a set of two-coupled transport equations that must be solved: one for protons and the other for H atoms. We present here a new approach that solves the multistream proton/hydrogen transport equations encompassing the collision angular redistributions and the magnetic mirroring effect. In order to validate our model we discuss the energy conservation and we compare with another model under the same inputs and with rocket observations. The influence of the angular redistributions is discussed in a forthcoming paper.

Folding processes and solitary waves in structural geology

Evolution of localized folding patterns in layered elastic and visco-elastic materials is reviewed in the context of compressed geological systems. The thin strut or plate embedded in a visco-elastic medium is used as an archetypal example to describe localized buckles which, in contrast to those from earlier formulations, appear in the absence of triggering imperfections. Structural and material effects are surveyed and important nonlinear characteristics are identified. A brief review of possible methods of analysis is conducted.

The visual ecology of avian photoreceptors

The spectral sensitivities of avian retinal photoreceptors are examined with respect to microspectrophotometric measurements of single cells, spectrophotometric measurements of extracted or in vitro regenerated visual pigments, and molecular genetic analyses of visual pigment opsin protein sequences. Bird species from diverse orders are compared in relation to their evolution, their habitats and the multiplicity of visual tasks they must perform. Birds have five different types of visual pigment and seven different types of photo receptor-rods, double (uneven twin) cones and four types of single cone. The spectral locations of the wavelengths of maximum absorbance (lambda (max)) of the different visual pigments, and the spectral transmittance characteristics of the intraocular spectral filters (cone oil droplets) that also determine photoreceptor spectral sensitivity, vary according to both habitat and phylogenetic relatedness. The primary influence on avian retinal design appears to be the range of wavelengths available for vision, regardless of whether that range is determined by the spectral distribution of the natural illumination or the spectral transmittance of the ocular media (cornea, aqueous humour, lens, vitreous humour). Nevertheless, other variations in spectral sensitivity exist that reflect the variability and complexity of avian visual ecology. (C) 2001 Elsevier Science Ltd. All rights reserved.

The structures of Langmuir-Blodgett films of fatty acids and their salts

Recent advances in several experimental techniques have enabled detailed structural information to be obtained for floating (Langmuir) monolayers and Langmuir-Blodgett films. These techniques are described briefly and their application to the study of films of fatty acids and their salts is discussed. Floating monolayers on aqueous subphases have been shown to possess a complex polymorphism with phases whose structures may be compared to those of smectic mesophases. However, only those phases that exist at high surface pressures are normally used in Langmuir-Blodgett (LB) deposition. In single LB monolayers of fatty acids and fatty acid salts the acyl chains are in the all-cans conformation with their long axes normal to the substrate. The in-plane molecular packing is hexagonal with long-range bond orientational order and short-range positional order: known as the hexatic-B structure. This structure is found irrespective of the phase of the parent floating monolayer. The structures of multilayer LB films are similar to the structures of their bulk crystals, consisting of stacked bilayer lamellae. Each lamella is formed from two monolayers of fatty acid molecules or ions arranged head to head and held together by hydrogen bonding between pairs of acids or ionic bonding through the divalent cations. With acids the acyl chains are tilted with respect to the substrate normal and have a monoclinic structure, whereas the salts with divalent cations may have the chains normal to the substrate or tilted. The in-plane structures are usually centred rectangular with the chains in the trans conformation and packed in a herringbone pattern, Multilayer films of the acids show only a single-step order-disorder transition at the malting point, This temperature tends to rise as the number of layers increases. Complex changes occur when multilayer films of the salts are heated. Disorder of the chains begins at low temperatures but the arrangement of the head groups does not alter until the melting temperature is reached, Slow heating to a temperature just below the melting temperature gives, with some salts, a radical change in phase. The lamellar structure disappears and a new phase consisting of cylindrical rods lying parallel to the substrate surface and stacked in a hexagonal pattern is formed, In each rod the cations are aligned along the central axis surrounded by the disordered acyl chains. (C) 2001 Elsevier Science B,V. All rights reserved.

Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

The effect of heat treatment on the toughness, hardness and microstructure of low carbon white cast irons

The effect of destabilisation and subcritical heat treatment on the impact toughness, hardness, and the amount and mechanical stability of retained austenite in a low carbon white cast iron have been investigated. The experimental results show that the impact energy constantly increases when the destabilisation temperature is raised from 950 degreesC to 1200 degreesC. Although the hardness decreases, the heat-treated hardness is still greater than the as-cast state. After destabilisation treatment at 1130 degreesC, tempering at 200 to 250 degreesC for 3 hours leads to the highest impact toughness, and secondary hardening was observed when tempering over 400 degreesC. The amount of retained austenite increased with the increase in the destabilisation temperature, and the treatment significantly improves the mechanical stability of the retained austenite compared with the as-cast state. Tempering below 400 degreesC does not affect the amount of retained austenite and its mechanical stability. But the amount of retained austenite is dramatically reduced when tempered above 400 degreesC. The relationship between the mechanical properties and the microstructure changes was discussed. (C) 2001 Kluwer Academic Publishers.