Albert de Montserrat Navarro, Jason Morgan, Robert Hall & Lloyd White

Background

The Earth has an amazingly flexible rheology, with brittle fault-like deformation predominating near the surface, and ductile deformation dominating in hotter deeper rocks.  Fault zones are thought to deform by brittle and ductile modes of deformation occurring at different spatial locations within the deforming lithosphere. Typically, numerical modelling studies of extension processes have used techniques that work well for idealizing brittle deformation (e.g. discrete element techniques), or ductile deformation (e.g. finite element techniques), but not both simultaneously.

BCs

Aims

This project will involve the development of a 2D numerical model that simultaneously includes both discrete element and finite element subdomains. To test this hybrid deformation code we will compare its predictions for the evolution of an extensional basinal system with previous analogue experiments and numerical models. Coding will be done using MATLAB and parallel-MATLAB, which allows us to combine these new deformation codetools with previous finite element codetools developed to model heat transport and melting processes. The resulting code will be used to model the evolution of extensional processes in one or more basins associated with Indonesian back-arc and within-arc regions that are currently being explored by the SE Asia Research Group.