DynEarthSol2D: An efficient unstructured finite element method to study long‐term tectonic deformation

[1] Many tectonic problems require to treat the lithosphere as a compressible elastic material, which can also flow viscously or break in a brittle fashion depending on the stress level applied and the temperature conditions. We present a flexible methodology to address the resulting complex material response, which imposes severe challenges on the discretization and rheological models used. This robust, adaptive, two-dimensional, finite element method solves the momentum balance and the heat equation in Lagrangian form using unstructured meshes. An implementation of this methodology is released to the public with the publication of this paper and is named DynEarthSol2D (available at http://bitbucket.org/tan2/dynearthsol2). The solver uses contingent mesh adaptivity in places where shear strain is focused (localization) and a conservative mapping assisted by marker particles to preserve phase and facies boundaries during remeshing. We detail the solver and verify it in a number of benchmark problems against analytic and numerical solutions from the literature. These results allow us to verify and validate our software framework and show its improved performance by an order of magnitude compared against an earlier implementation of the Fast Lagrangian Analysis of Continua algorithm.

[1]  Robert B. Trumbull,et al.  Numerical models of crustal scale convection and partial melting beneath the Altiplano–Puna plateau , 2002 .

[2]  C. Thieulot,et al.  DOUAR: A new three-dimensional creeping flow numerical model for the solution of geological problems , 2008 .

[3]  John A. Trangenstein,et al.  A second-order Godunov algorithm for two-dimensional solid mechanics , 1994 .

[4]  Marie-Pierre Doin,et al.  A comparison of methods for the modeling of thermochemical convection , 1997 .

[5]  A. Poliakov,et al.  Factors controlling normal fault offset in an ideal brittle layer , 2000 .

[6]  Eunseo Choi,et al.  Thermomechanics of mid-ocean ridge segmentation , 2008 .

[7]  David A. Yuen,et al.  Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems , 2007 .

[8]  Vladimir Lyakhovsky,et al.  A rheological model of a fractured solid , 1993 .

[9]  Jean Braun,et al.  Dynamical Lagrangian Remeshing (DLR): A new algorithm for solving large strain deformation problems and its application to fault-propagation folding , 1994 .

[10]  R Stephenson,et al.  Small-Scale Mantle Convection Produces Stratigraphic Sequences in Sedimentary Basins , 2010, Science.

[11]  P Cundall,et al.  A MICROCOMPUTER PROGRAM FOR MODELLING LARGE-STRAIN PLASTICITY PROBLEMS. PROCEEDINGS OF THE SIXTH INTERNATIONAL CONFERENCE ON NUMERICAL METHODS IN GEOMECHANICS, 11-15 APRIL 1988, INNSBRUCK, AUSTRIA. VOLUMES 1 - 3 , 1988 .

[12]  M. Handy,et al.  The solid‐state flow of polymineralic rocks , 1990 .

[13]  Stuart S. Nagel.  A Microcomputer Program for Decision-making , 1991 .

[14]  Phillip Colella,et al.  A higher-order Godunov method for modeling finite deformation in elastic-plastic solids , 1991 .

[15]  C. Scholz.  The Mechanics of Earthquakes and Faulting , 1990 .

[16]  Peter B. Flemings,et al.  Large‐scale stratigraphic architecture, eustatic variation, and unsteady tectonism: A theoretical evaluation , 1991 .

[17]  Magali I. Billen,et al.  Rheologic controls on the dynamics of slab detachment , 2009 .

[18]  C. Beaumont,et al.  Asymmetric lithospheric extension: The role of frictional plastic strain softening inferred from numerical experiments , 2002 .

[19]  Louis Moresi,et al.  Timing of formation of Beta Regio and its geodynamical implications , 2003 .

[20]  W. R. Buck,et al.  Flexural rotation of normal faults , 1988 .

[21]  B. Parsons,et al.  The motion of crustal blocks driven by flow of the lower lithosphere and implications for slip rates of continental strike-slip faults , 1998, Nature.

[22]  Jun Korenaga,et al.  Thermal cracking and the deep hydration of oceanic lithosphere: A key to the generation of plate tectonics? , 2007 .

[23]  Jian Lin,et al.  A continuum mechanics model for normal faulting using a strain-rate softening rheology: implications for thermal and rheological controls on continental and oceanic rifting , 2002 .

[24]  Vladimir Lyakhovsky,et al.  Deformation and seismicity associated with continental rift zones propagating toward continental margins , 2012 .

[25]  Christopher H. Scholz,et al.  The brittle-plastic transition and the depth of seismic faulting , 1988 .

[26]  Hans Muhlhaus,et al.  A Lagrangian integration point finite element method for large deformation modeling of viscoelastic geomaterials , 2003 .

[27]  Katia Mocellin,et al.  A new efficient explicit formulation for linear tetrahedral elements non‐sensitive to volumetric locking for infinitesimal elasticity and inelasticity , 2009 .

[28]  Ted Belytschko,et al.  On the dynamic effects of explicit FEM in sheet metal forming analysis , 1998 .

[29]  Yuri Y. Podladchikov,et al.  FRACTAL PLASTIC SHEAR BANDS , 1994 .

[30]  Ryota Hino,et al.  Mantle wedge structure in the Miyagi Prefecture forearc region, central northeastern Japan arc, and its relation to corner‐flow pattern and interplate coupling , 2011 .

[31]  Laurent G. J. Montési,et al.  Spacing of faults at the scale of the lithosphere and localization instability: 1. Theory , 2003 .

[32]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[33]  B. Gruzinskaya,et al.  What drives orogeny in the Andes , 2005 .

[34]  David H. Sharp,et al.  A conservative Eulerian formulation of the equations for elastic flow , 1988 .

[35]  David H. Sharp,et al.  A conservative formulation for plasticity , 1992 .

[36]  P. A. Cundall,et al.  An Explicit Inertial Method for the Simulation of Viscoelastic Flow: An Evaluation of Elastic Effects on Diapiric Flow in Two- and Three- Layers Models , 1993 .

[37]  Giorgio Ranalli,et al.  Rheology of the lithosphere in space and time , 1997, Geological Society, London, Special Publications.

[38]  S. V. Sobolev,et al.  SLIM3D: A tool for three-dimensional thermomechanical modeling of lithospheric deformation with elasto-visco-plastic rheology , 2008 .

[39]  Paul J. Tackley,et al.  Self‐consistent generation of tectonic plates in time‐dependent, three‐dimensional mantle convection simulations , 2000 .

[40]  Eduardo N. Dvorkin,et al.  An Eulerian finite element formulation for modelling stationary finite strain elastic deformation processes , 2005 .

[41]  Phillip Colella,et al.  A high-order Eulerian Godunov method for elastic-plastic flow in solids , 2001 .

[42]  Victor M. Calo,et al.  A finite strain Eulerian formulation for compressible and nearly incompressible hyperelasticity using high‐order B‐spline finite elements , 2012 .

[43]  Simon M. Peacock,et al.  Serpentinization of the forearc mantle , 2003 .

[44]  Boris J. P. Kaus,et al.  A stabilization algorithm for geodynamic numerical simulations with a free surface , 2010 .

[45]  Vincent Godard,et al.  Incorporating metamorphism in geodynamic models: the mass conservation problem , 2011 .

[46]  Cian R. Wilson,et al.  Fluidity: A fully unstructured anisotropic adaptive mesh computational modeling framework for geodynamics , 2011 .

[47]  Ali Karrech,et al.  Thermal‐elastic stresses and the criticality of the continental crust , 2012 .

[48]  Arie P. van den Berg,et al.  Latent heat effects of the major mantle phase transitions on low-angle subduction , 2001 .

[49]  Arie P. van den Berg,et al.  The influence of rheological weakening and yield stress on the interaction of slabs with the 670 km discontinuity , 2002 .

[50]  Marcin Dabrowski,et al.  MILAMIN: MATLAB‐based finite element method solver for large problems , 2008 .

[51]  Patrick E. Farrell,et al.  Conservative interpolation between volume meshes by local Galerkin projection , 2011 .

[52]  Andrew J. Dombard,et al.  A test of the validity of yield strength envelopes with an elastoviscoplastic finite element model , 2000 .

[53]  P. Fullsack.  An arbitrary Lagrangian-Eulerian formulation for creeping flows and its application in tectonic models , 1995 .

[54]  P. Cundall.  Numerical experiments on localization in frictional materials , 1989 .

[55]  A. Poliakov,et al.  Mechanics of Stretching Elastic‐Plastic‐Viscous Layers: Applications to Slow‐Spreading Mid‐Ocean Ridges , 2013 .

[56]  W. Buck,et al.  Effect of lithospheric thickness on the formation of high- and low-angle normal faults , 1993 .

[57]  A. B. WATTS,et al.  Isostasy and Flexure of the Lithosphere , 2001 .

[58]  Gabriele Morra,et al.  A model comparison study of large-scale mantle–lithosphere dynamics driven by subduction , 2008 .

[59]  Boris J. P. Kaus,et al.  Factors that control the angle of shear bands in geodynamic numerical models of brittle deformation , 2010 .

[60]  C. Scholz,et al.  Dilatancy in the fracture of crystalline rocks , 1966 .

[61]  Taras Gerya,et al.  Discretization errors and free surface stabilization in the finite difference and marker‐in‐cell method for applied geodynamics: A numerical study , 2011 .

[62]  P. Colella,et al.  A conservative three-dimensional Eulerian method for coupled solid-fluid shock capturing , 2002 .

[63]  D. Culler,et al.  Comparison of methods , 2000 .

[64]  Luc L. Lavier,et al.  A mechanism to thin the continental lithosphere at magma-poor margins , 2006, Nature.

[65]  B. Evans,et al.  Strength of the lithosphere: Constraints imposed by laboratory experiments , 1995 .

[66]  Anthony Watts,et al.  Lithospheric strength and its relationship to the elastic and seismogenic layer thickness , 2003 .

[67]  Eh Tan,et al.  The role of frictional strength on plate coupling at the subduction interface , 2012 .

[68]  Eh Tan,et al.  The role of frictional strength on plate coupling at the subduction interface , 2012 .

[69]  G. Suhr,et al.  Stacked gabbro units and intervening mantle: A detailed look at a section of IODP Leg 305, Hole U1309D , 2008 .