Discrete Element Method model data of biaxial shear deformation

Owing to their destructive potential, earthquakes receive considerable attention from laboratory studies. In friction experiments, stick-slips are studied as the laboratory equivalent of natural earthquakes, and numerous attempts have been made to simulate stick-slips numerically using the Discrete Element Method (DEM). However, while laboratory stick-slips commonly exhibit regular stress drops and recurrence times, stick-slips generated in DEM simulations are highly irregular. This discrepancy highlights a gap in our understanding of stick-slip mechanics, which propagates into our understanding of earthquakes. In this work, we show that regular stick-slips emerge in DEM when time-dependent compaction by pressure solution is considered. We further show that the stress drop and recurrence time of stick-slips is directly controlled by the kinetics of pressure solution. Since compaction is known to operate in faults, this mechanism for frictional instabilities directly relates to natural seismicity.
The zip-fle contains a Python script (render_figures.py) that is used to generate the data fgures as reported by Van den Ende & Niemeijer (2018), auxiliary script fles in the scripts directory, and the original model data in ASCII and HDF format in the data directory. The main Python script fle render_figures.py will read and process the original model data and generate the interactive data fgures. These fgures are automatically saved in PDF format. More information is given in Van den Ende & Niemeijer (2018) to which these data and scripts are supplementary material to.