[DATASET] How do rift-related fault network distributions evolve?
Here we provide the Supplementary Material for the publication 'How do rift-related fault network distributions evolve? Quantitative comparisons between natural fault observations and 3D numerical models of continental extension' by Pan, Naliboff, Bell and Jackson in Tectonics.
(1) Models: these contain the five models (A-E) as an interpolated numpy grid. The grid is a 2D array given that we extract a 2D depth slice 5 km beneath the initial model grid. The original ASPECT parameter file is also included for each model.
(2) Scripts: contains the python scripts to generate fault statistics. All scripts use the interpolated numpy grid as a model input. Use the input.py file to dictate key parameter inputs (e.g. model extent) and run ‘fault_extraction_analysis_2D.py’ to execute the main script. Comments are provided in the script, but two jupyter notebooks provide an additional step-by-step context of the fault extraction workflow and reproduction of figures.
(3) Extracted fault statistics: are the generated outputs from the Python scripts. ‘Geometric_relationships’ is the main excel database generated from ‘fault_extraction_analysis_2D.py’ whereby each row represents the data for one single fault (or label in python) and contains coordinates, strain, length etc. statistics. Alternatively there is an option to generate an equivalent excel file, but whereby each row represents an x, y coordinate (with the associated label represented as a unique accessible index) – this option works better for plotting the results and the option to generate this is given in the input.py script and jupyter notebooks.
In the publication, we plot (a) overall strain summation of the study area, (b) fault accommodated strain and (c) OFD (i.e. the difference between a and b). The excel data for each of these are present as subfolders, and were extracted from the main python script, as well as additional scripts 'results_concat.py’ and ‘results_to_strain_time.py’ which process the excel files in ‘Geometric_relationships’. The calculation and plotting of OFD is shown in the jupyter notebook ‘OFF FAULT DEFORMATION calculations.ipynb’.