Coupled models of half-graben growth and river incision

The code.zip contains all the MATLAB scripts necessary to run the simulations presented in: Olive, J.-A., Malatesta, L.C., Behn, M.D., and Buck, W.R., Sensitivity of rift tectonics to global variability in the efficiency of river erosion, PNAS, 2022 (https://doi.org/10.1073/pnas.2115077119).

These simulations couple the geodynamic modeling code SiStER (https://github.com/jaolive/SiStER) with a landscape evolution model accounting for river incision and sedimentary infilling. The input files are set up to model slip on a half-graben master fault, ultimately leading to spontaneous strain re-localization.

The landscape evolution routines correspond to the files starting with "LEM-", and are based on a MATLAB implementation of the Fastscape algorithm (Braun & Willett, 2013) by B. Kaus and A. Popov. Each simulation corresponds to an input file in which the only parameters that are varied are the crustal radiogenic heat production H, the erodibility coefficient K, and whether sediment deposition is included or not. A summary of all simulation inputs is provided in Table S1 of Olive et al. (under review). Simulations can be launched in a MATLAB prompt by running the command:

>> SiStER_RUN('SiStER_Input_File_name')

where input file names match the simulation labels, e.g., SiStER_Input_File_half_graben_H60.m corresponds to simulation H60. Further details on how to run the code and visualize the outputs are accessible at https://github.com/jaolive/SiStER.