Improving subduction interface implementation in dynamic numerical models

Abstract

Numerical subduction models often implement an entrained weak layer (WL) to facilitate decoupling of the slab and upper plate. This approach is attractive in its simplicity, and can provide stable, asymmetric subduction systems that persist for many tens of millions of years. In this study we undertake a methodological analysis of the WL approach, and use these insights to guide improvements to the implementation. The issue that primarily motivates the study is the emergence of significant spatial and temporal thickness variations within the WL. We show that these variations are mainly the response to volumetric flux gradients, caused by the change in boundary conditions as the WL material e..