Monte Carlo code for modelling the transient creep of different minerals

Postseismic relaxation from large earthquakes induce transient deformation of the solid Earth, particularly in the deeper part of the crust. The deformation of the upper and lower crust are mainly controlled by quartz and feldspar, respectively. The mechanical properties of quartz and feldspar at steady-state are well constrained and flow law parameters are known. However, the physical mechanism underlying transient creep is poorly understood and the corresponding flow law parameters are unknown. Here, we constrain a constitutive framework that captures transient creep and steady state creep consistently using the mechanical data from laboratory experiments. The constitutive framework represents a Burgers assembly with a thermally activated nonlinear stress versus strain-rate relationship for the dashpots. Using the Markov chain Monte Carlo (MCMC) method, we uniquely determine the flow law parameters for both quartz and feldspar. This repository contains the code and data that can be used to replicate the results