Scripting High Performance Earth Systems Simulations on the SGI Altix 3700

A Python language scripting interface for optimised numerical software written in C/C++ has been designed to facilitate the rapid development of 3D parallel simulations on the Altix 3700. A recent advance focusing on the extensions implemented within a high performance finite element kernel module is described. Performance issues, measurements and results on the Altix 3700 are presented

The Influence Of Elasticity, Temperature And Fracture On Large Scale Geological Flow

Realistic simulations of earth processes such as faulting, shearing, magma flow, subduction and convection often require the consideration of non-Newtonian Effects such as elasticity and power law creep. As the deformations involved in geological deformation are often large the constitutive relationships must maintain certain geometric terms to ensure that the tensor properties of the model are conserved. A model with such properties is termed as objective. There are a wide range of objective, visco-elasto-plastic models to choose from. The main structural difference between these models consists in the choice of the objective stress rate, e.g. Jaumann, Oldroyd, Truesdell - rates (see Kolymbas and Herle, 2003, for a recent discussion). In this paper we give an outline of a thermo-visco-elastic-plastic model including a discussion of numerical aspects such as the derivation of a consistent incremental form. The viscous part of the deformation involves a combination of both Newtonian and power law creep. Plastic deformations are described by means of a standard Prandtl-Reuss flow rule combined with a von Mises yield criterion. In planetary scale flow modeling the yield criterion is required as a stress limiter during episodic events e.g. in connection with the initiation of subduction. The temperature sensitivity of the viscous deformation is considered by means of an Arrhenius relation involving a pressure dependent reference (melting) temperature. The salient features of the model are first explored by means of analytical and numerical solutions of a simple shear problem for an infinite strip with fixed and prescribed shear velocities on the bottom and top of the layer respectively

Multicycle dynamics of fault systems and static and dynamic triggering of earthquakes

Dynamic simulations of rupture propagation and multiple earthquake cycles for varying fault geometries are presented. We investigate the role of both dynamic and static stress changes on earthquake triggering. Dynamic stress triggering of earthquakes is caused by the passage of seismic waves, whereas static stress triggering is due to net slippage on a fault resulting from an earthquake. Static stress changes represented by a Coulomb failure function and its relationship to seismicity rate change is a relatively well-known mechanism, whereas the physical origin of dynamic triggering remains one of the least understood aspects of earthquake nucleation. We investigate these mechanisms by analysing seismicity patterns with varying fault separation, geometry and with and without dynamic triggering present

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

Universal and wide shear zones in granular bulk flow

We present experiments on slow granular flows in a modified (split-bottomed) Couette geometry in which wide and tunable shear zones are created away from the sidewalls. For increasing layer heights, the zones grow wider (apparently without bound) and evolve towards the inner cylinder according to a simple, particle-independent scaling law. After rescaling, the velocity profiles across the zones fall onto a universal master curve given by an error function. We study the shear zones also inside the material as function of both their local height and the total layer height.Comment: Minor corrections, accepted for PRL (4 pages, 6 figures

An anisotropic viscous representation of Mohr-Coulomb failure for use is modelling coupled mantle-continent dynamics

In mantle convection models it has become common to make use of a modified (pressure sensitive, Boussinesq) von Mises yield criterion to limit the maximum stress the lithosphere can support. This approach allows the viscous, cool thermal boundary layer to deform in a relatively plate-like mode even in a fully Eulerian representation. In large-scale models with embedded continental crust where the mobile boundary layer represents the oceanic lithosphere, the von Mises yield criterion for the oceans ensures that the continents experience a realistic broad-scale stress regime. In detailed models of crustal deformation it is, however, more appropriate to choose a Mohr-Coulomb yield criterion based upon the idea that frictional slip occurs on whichever one of many randomly oriented planes happens to be favorably oriented with respect to the stress field. As coupled crust/mantle models become more sophisticated it is important to be able to use whichever failure model is appropriate to a given part of the system. We have therefore developed a way to represent Mohr-Coulomb failure within a code which is suited to mantle convection problems coupled to large-scale crustal deformation. Our approach uses an orthotropic viscous rheology (a different viscosity for pure shear to that for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation - neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane which is formally equivalent in any symmetric formulation. It is also important to allow strain-weakening of the material. The material should remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-Integration-point finite element code and will show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion, and comparisons between mantle convection models using the von Mises versus the Mohr-Coulomb yield criteria. The formulation itself is general and applies to 2D and 3D problems, although it is somewhat more complicated to identify the slip plane in 3D