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.

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.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Estimation of crystalline phase present in the glucose crystal-solution mixture by water activity measurement

The amount of crystalline fraction present in monohydrate glucose crystal-solution mixture up to 110% crystal in relation to solution (crystal:solution=110:100) was determined by water activity measurement. It was found that the water activity had a strong linear correlation (R-2=0.994) with the amount of glucose present above saturation. Difference in the water activities of the crystal-solution mixture (a(w1)) and the supersaturated solution (a(w2)) by re-dissolving the crystalline fraction allowed calculation of the amount of crystalline phase present (DeltaG) in the mixture by an equation DeltaG=846.97(a(w1)-a(w2)). Other methods such as Raoult's, Norrish and Money-Born equations were also tested for the prediction of water activity of supersaturated glucose solution. (C) 2003 Swiss Society of Food Science and Technology. Published by Elsevier Science Ltd. All rights reserved.

Seedling responses of four Australian tree species to toxic concentrations of manganese in solution culture

Little is known about the responses of Australian plants to excess metal, including Mn. It is important to remedy this lack of information so that knowledgeable decisions can be made about managing Mn contaminated sites where inhabited by Australian vegetation. Acacia holosericea, Melaleuca leucadendra, Eucalyptus crebra and Eucalyptus camaldulensis were grown in dilute solution culture for 10 weeks. The seedlings ( 42 days old) were exposed to six Mn treatments viz., 1, 8, 32, 128, 512 and 2048 muM. The order of tolerance to toxic concentrations of Mn was A. holosericea congruent to = E. crebra < M. leucadendra < E. camaldulensis, the critical external concentrations being approximately 5.1, 5.0, 21 and 330 muM, respectively. The critical tissue Mn concentrations for the youngest fully expanded leaf and total shoots were, respectively, 265 and 215 mug g(-1) DM for A. holosericea, 445 and 495 mug g(-1) DM for M. leucadendra, 495 and 710 mug g(-1) DM for E. crebra and 7230 and 6510 mug g(-1) DM for E. camaldulensis. The high tolerance of E. camaldulensis ( as opposed to the sensitivity of E. crebra) to excess Mn raises concern about fauna feeding on the plant and is consistent with hypotheses suggesting the Eucalyptus subgenus Symphomyrtus is particularly tolerant of stress, including excess Mn. The results from this paper provide the first comprehensive combination of growth responses, critical external concentrations, critical tissue concentrations and plant toxicity symptoms for three important Australian genera, viz., Eucalyptus, Acacia and Melaleuca, for use in the management of Mn toxic sites.

Integrability and exact solution of an electronic model with long range interactions

We present an electronic model with long range interactions. Through the quantum inverse scattering method, integrability of the model is established using a one-parameter family of typical irreducible representations of gl(211). The eigenvalues of the conserved operators are derived in terms of the Bethe ansatz, from which the energy eigenvalues of the Hamiltonian are obtained.

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.

Solution structure of χ-conopeptide MrIA, a modulator of the human norepinephrine transporter

The chi-conopeptides MrIA and MrIB are 13-residue peptides with two disulfide bonds that inhibit human and rat norepinephrine transporter systems and are of significant interest for the design of novel drugs involved in pain treatment. In the current study we have determined the solution structure of MrIA using NMR spectroscopy. The major element of secondary structure is a hairpin with the two strands connected by an inverse gamma-turn. The residues primarily involved in activity have previously been shown to be located in the turn region (Sharpe, I. A.; Palant, E.: Schroder, C. L; Kaye, D. M.; Adams, D. I.; Alewood, P. F.; Lewis, R. J. J Biol Client 2003, 278, 40317-40323), which appears to be more flexible than the beta-strands based on disorder in the ensemble of calculated structures. Analogues of MrIA with N-terminal truncations indicate that the N-terminal residues play a role in defining a stable conformation and the native disulfide connectivity. In particular, noncovalent interactions between Val3 and Hypl2 are likely to be involved in maintaining a stable conformation. The N-terminus also affects activity, as a single N-terminal deletion introduced additional pharmacology at rat vas deferens, while deleting the first two amino acids reduced chi-conopeptide potency. This article was originally published online as an accepted preprint. The Published Online date corresponds to the preprint version. You can request a copy of the preprint by entailing the Biopolymers editorial office at biopolymers@wiley.com (c) 2005 Wiley Periodicals, Inc.

Numerical simulation of dendrite growth during solidification

A comprehensive probabilistic model for simulating dendrite morphology and investigating dendritic growth kinetics during solidification has been developed, based on a modified Cellular Automaton (mCA) for microscopic modeling of nucleation, growth of crystals and solute diffusion. The mCA model numerically calculated solute redistribution both in the solid and liquid phases, the curvature of dendrite tips and the growth anisotropy. This modeling takes account of thermal, curvature and solute diffusion effects. Therefore, it can simulate microstructure formation both on the scale of the dendrite tip length. This model was then applied for simulating dendritic solidification of an Al-7%Si alloy. Both directional and equiaxed dendritic growth has been performed to investigate the growth anisotropy and cooling rate on dendrite morphology. Furthermore, the competitive growth and selection of dendritic crystals have also investigated.

A numerical study of effect of grain boundaries on elastic and plastic properties in nanocomposite materials

Nanocomposite materials have received considerable attention in recent years due to their novel properties. Grain boundaries are considered to play an important role in nanostructured materials. This work focuses on the finite element analysis of the effect of grain boundaries on the overall mechanical properties of aluminium/alumina composites. A grain boundary is incorporated into the commonly used unit cell model to investigate its effect on material properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the ''effective'' plastic property of the grain boundary is estimated. In addition, the strengthening mechanism is also discussed based on the Estrin-Mecking model.

Determination of the solution structures of conantokin-G and conantokin-T by CD and NMR spectroscopy

Conantokin-G and conantokin-T are two paralytic polypeptide toxins originally isolated from the venom of the fish-hunting cone snails of the genus Conus. Conantokin-G and conantokin-T are the only naturally occurring peptidic compounds which possess N-methyl-D-aspartate receptor antagonist activity, produced by a selective non-competitive antagonism of polyamine responses, They are also structurally unusual in that they contain a disproportionately large number of acid labile post-translational gamma-carboxyglutamic acid (Gla) residues, Although no precise structural information has previously been published for these peptides, early spectroscopic measurements have indicated that both conantokin-G and conantokin-T form alpha-helical structures, although there is some debate whether the presence of calcium ions is required for these peptides to adopt this fold, We now report a detailed structural study of synthetic conantokin-G and conantokin-T in a range of solution conditions using CD and H-1 NMR spec troscopy. The three-dimensional structures of conantokin-T and conantokin-G were calculated from H-1 NMR-derived distance and dihedral restraints. Both conantokins were found to contain a mixture of alpha- and 3(10) helix, that give rise to curved and straight helical conformers. Conantokin-G requires the presence of divalent cations (Zn2+, Ca2+, Cu2+, Or Mg2+) to form a stable iv-helix, while conantokin-T adopts a stable alpha-helical structure in aqueous conditions, in the presence or absence of divalent cations (Zn2+, Ca2+, Cu2+, Or Mg2+).