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 randomized trial comparing hormone replacement therapy (HRT) and HRT plus calcitriol in the treatment of postmenopausal osteoporosis with vertebral fractures: Benefit of the combination on total body and hip density

We report a prospective, randomized, multi-center, open-label 2-year trial of 81 postmenopausal women aged 53-79 years with at least one minimal-trauma vertebral fracture (VF) and low (T-score below 2) lumbar bone mineral density (BMD). Group HRT received piperazine estrone sulfate (PES) 0.625 - 1.25 mg/d +/- medroxyprogesterone acetate (MPA) 2.5 - 5 mg/d,- group HRT/D received HRT plus calcitriol 0.25 mug bd. All with a baseline dietary calcium (Ca) of < I g/d received Ca carbonate 0.6 g nocte. Final data were on 66 - 70 patients. On HRT/D, significant (P < 0.001) BNID increases from baseline by DXA were at total body - head, trochanter, Ward's, total hip, inter-trochanter and femoral shaft (% group mean Delta 4.2, 6.1, 9.3. 3.7. 3.3 and 3.3%, respectively). On HRT, at these significant Deltas were restricted to the trochanter and sites. si Wards. Significant advantages of HRT/D over HRT were in BMD of total body (- head), total hip and trochanter (all P = 0.01). The differences in mean Delta at these sites were 1.3, 2.6 and 3.9%. At the following, both groups Improved significantly -lumbar spine (AP and lateral), forearm shaft and ultradistal tibia/fibula. The weightbearing, site - specific benefits of the combination associated with significant suppression of parathyroid hormone-suggest a beneficial effect on cortical bone. Suppression of bone turnover was significantly greater on HRT/D (serum osteocalcin P = 0.024 and urinary hydroxyproline/creatinine ratio P = 0.035). There was no significant difference in the number of patients who developed fresh VFs during the trial (HRT 8/36, 22%; HRT/D 4/34, 12% - intention to treat); likewise in the number who developed incident nonvertebral fractures. This Is the first study comparing the 2 treatments in a fracture population. The results indicate a significant benefit of calcitriol combined with HRT on total body BMD and on BNID at the hip, the major site of osteoporotic fracture.

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.

Effects of settler size and density on early post-settlement survival of Ciona intestinalis in the field

Effects of variation in larval quality on post-metamorphic performance in marine invertebrates are increasingly apparent. Recently, it has been shown that variation in offspring size can also strongly affect post-settlement survival, but variation in environmental conditions can mediate this effect. The quality of habitat into which marine invertebrate larvae settle can vary markedly, and 1 influence on quality is the number of conspecifics present. We tested the effects of settler size and settler density on early (1 wk after settlement) post-settlement survival in the field for the solitary ascidian Ciona intestinalis. Larger settlers survived better than smaller settlers, within and among groups of siblings. Increases in the density of settlers decreased survival, but the density-dependent effects were much stronger for smaller settlers. We suggest that larger settlers are better able to cope with intra-specific competition because they have greater energetic reserves or a greater capacity to feed than smaller settlers.

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.

Use of bulk density profiles from X-radiography to examine structural crust models

Exponential and sigmoidal functions have been suggested to describe the bulk density profiles of crusts. The present work aims to evaluate these conceptual models using high resolution X-radiography. Repacked seedbeds from two soil materials, air-dried or prewetted by capillary rise, were subjected to simulated rain, which resulted in three types of structural crusts, namely, slaking, infilling, and coalescing. Bulk density distributions with depth were generated using high-resolution (70 mum), calibrated X-ray images of slices from the resin-impregnated crusted seedbeds. The bulk density decreased progressively with depth, which supports the suggestion that a crust should be considered as a nonuniform layer. For the slaking and the coalescing crusts, the exponential function underestimated the strong change in bulk density across the morphologically defined transition between the crust and the underlying material; the sigmoidal function provided a better description. Neither of these crust models effectively described the shape of the bulk density profiles through the whole seedbed. Below the infilling and slaking crusts, bulk density increased linearly with depth as a result of slumping. In the coalescing crusted seedbed, the whole seedbed uniformly collapsed and most of the bulk density change within the crust could be ascribed to slumping (0.33 g cm(-3)) rather than to crusting (0.12 g cm(-3)). Finally, (i) X-radiography appears as a unique tool to generate high resolution bulk density profiles and (ii) in structural crusts, bulk density profiles could be modeled using the existing exponential and sigmoidal crusting models, provided a slumping model would be coupled.

Low prenatal vitamin D alters morphology, cell density and gene expression in adult rat brain: Correlations with schizophrenia

Statement of the study: Based on data from ecological and analytic epidemiological studies, we have proposed that low prenatal vitamin D is a candidate risk-modifying factor for schizophrenia. Previously, we demonstrated that low prenatal vitamin D adversely affected brain development in neonatal rats (Eyles et al, 2003). Here we examine the impact of both prenatal and early life hypovitaminosis D on various outcomes in the adult rat brain. Methods: Female Sprague-Dawley rats were made vitamin D deficient via the use of a special diet (Dyets CA) and lighting conditions that excluded UVB radiation. Animals were kept under these conditions for 6 weeks then mated with males kept under normal conditions. Vitamin deplete dams were kept under these conditions during pregnancy. Offspring from two test groups were examined. Offspring were either reared with dams repleted with vitamin D at birth or remained under deplete conditions till weaning. Both test groups were weaned under normal vitamin D conditions and remained so till testing at adulthood. We compared the brains of adult offspring kept under both test conditions with animals from control environments. Summary of results: We found a significant persistent dose-related increase in lateral ventricle volume and alterations in anterior cingulate and prefrontal cortical cell densities (consistent with the known prodifferentiation properties of this steroid). In both test groups we observed a reduced expression of NGF as well as a down-regulation of transcripts coding for GABAA alpha 4 receptor and two neuronal structural elements; MAP2 and Neurofilament L. Conclusion: These findings provide further evidence that vitamin D is involved in brain development. An increase in prefrontal cortical cell density, a reduction neuronal structural elements and persistent ventriculomegaly are all common anatomical findings in the brains of patients with schizophrenia. The specific reduction in transcripts for neuronal structural proteins but not GFAP is also in accordance with the proposal that frontal cortical architecture in schizophrenia reflects a reduction in connectivity rather than a reduction in glial processes(Goldman-Rakic and Selemon, 1997). These findings confirm the biological plausibility of early life hypovitaminosis D as a risk factor for schizophrenia.

Finite-temperature correlations and density profiles of an inhomogeneous interacting one-dimensional Bose gas

We calculate the density profiles and density correlation functions of the one-dimensional Bose gas in a harmonic trap, using the exact finite-temperature solutions for the uniform case, and applying a local density approximation. The results are valid for a trapping potential that is slowly varying relative to a correlation length. They allow a direct experimental test of the transition from the weak-coupling Gross-Pitaevskii regime to the strong-coupling, fermionic Tonks-Girardeau regime. We also calculate the average two-particle correlation which characterizes the bulk properties of the sample, and find that it can be well approximated by the value of the local pair correlation in the trap center.

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.

Modelling of species in hood fires by conditional moment closure

Carbon monoxide, the chief killer in fires, and other species are modelled for a series of enclosure fires. The conditions emulate building fires where CO is formed in the rich, turbulent, nonpremixed flame and is transported frozen to lean mixtures by the ceiling jet which is cooled by radiation and dilution. Conditional moment closure modelling is used and computational domain minimisation criteria are developed which reduce the computational cost of this method. The predictions give good agreement for CO and other species in the lean, quenched-gas stream, holding promise that this method may provide a practical means of modelling real, three-dimensional fire situations. (c) 2005 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

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.

Tidal influence on behaviour of a coastal aquifer adjacent to a low-relief estuary

The tidal influence on groundwater hydrodynamics, salt-water intrusion and submarine groundwater discharge from coastal/estuarine aquifers is poorly quantified for systems with a mildly sloping beach, in contrast to the case where a vertical beach face is assumed. We investigated the effect of beach slope for a coastal aquifer adjacent to a low-relief estuary, where industrial waste was emplaced over the aquifer. The waste was suspected to discharge leachate towards the estuary. Field observations at various locations showed that tidally induced groundwater head fluctuations were skewed temporally. Frequency analysis suggested that the fluctuation amplitudes decreased exponentially and the phase-tags increased Linearly for the primary tidal signals as they propagated inland. Salinisation zones were observed in the bottom part of the estuary and near the beach surface. Flow and transport processes in a cross-section perpendicular to the estuary were simulated using SEAWAT-2000, which is capable of depicting density-dependent flow and multi-species transport. The simulations showed that the modelled water table fluctuations were in good agreement with the monitored data. Further simulations were conducted to gain insight into the effects of beach slope. In particular the limiting case of a vertical beach face was considered. The simulations showed that density difference and tidal forcing drive a more complex hydrodynamic pattern for the mildly sloping beach than the vertical beach, as well as a profound asymmetry in tidally induced water table fluctuations and enhanced salt-water intrusion. The simulation results also indicated that contaminant transport from the aquifer to the estuary was affected by the tide, where for the mildly sloping beach, the tide tended to intensify the vertical mass exchange in the vicinity of the shorelines, (c) 2005 Elsevier B.V. All rights reserved.

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.