Code and data accompanying the paper entitled "Seismic fault slip behaviour predicted from internal microphysical processes"

Improving the basis for earthquake hazard assessemt relies in part on structural observation of natural fault zones, on laboratory experiments, and on theoretical developments. While quantitative models that reproduce the laboratory data on frictional slip on faults is the minimum necessary to simulate their mechanical behavior in nature, convincing models must ultimately also account for all key fault zone observations. Taking a fault in carbonate rock as an example, and using fundamental data on the microscale processes that lead to fault friction, this work takes a first step towards simulating the dynamic evolution and self-organization of internal fault zone (micro)structure, processes and properties during a seismic slip event. The results show that the model captures many key observations on seismic fault behaviour, paving the way for still further improvements in developing a physics-based understanding of earthquake rupture in future. The present dataset includes the raw codes of four cases that are implemented in the finite element model (FEM) package Comsol (version 3.5a), and the numerical data produced by the codes (Results I to IV).