Pecube-D: Thermokinematic and Erosion Modeling Software for problems in Tectonics and Surface Processes
This repository submission provides the source code for the Pecube-D software developed by the workgroup and collaborators of Todd Ehlers. The source code is based on the original version of Pecube published by Jean Braun (Braun, 2003). It has been modified in several ways to simulate the predicted thermal history and thermochronometer age response for diverse problems in tectonic and surface process studies.
Pecube-D (where the -D represents ‘detrital’ thermochronometer ages) uses the finite element method to solve the 3D advection-diffusion equation with a crustal radiogenic heat source. The governing equations are solved for constant temperature lower and surface boundary conditions. The upper surface of the model (topography) is allowed to vary with time following the approach explained in Braun (2003). Modifications made to the original version of Pecube and included in this distribution include: 1) the transient calculation of detrial thermochronometer ages for user defined river or basin sample locations (e.g., Whipp et al., 2009, JGR); 2) Monte-Carlo and genetic inverse search algorithms for inversion of user-defined cooling ages for transient exhumation histories (e.g., Thiede and Ehlers, 2013, EPSL); 3) inversion of bedrock thermochronometer ages for paleotopography and paleotopographic change (e.g., Olen et al., 2012, JGR-ES; Ehlers et al., 2006, Geology); 4) integration with the Move software to simulate thermochronometer ages as a function of spatially and temporally variable structural and kinematic fields (e.g., McQuarrie and Ehlers, 2015, Tectonics); 5) diverse kinematic scenarios to evaluate the thermal and thermochronometer response to Himalayan and syntaxis rock exhumation (e.g., Michel et al., 2018, Geology); and 6) modifications to how different thermochronometer ages are predicted (e.g., Ehlers, 2005, RMG; and also as described in the previous references).
Interested users are recommended to start with the included Pecube.in file, which provides an instruction manual via copious comments, as well as instructions for compilation of the source code. Compilation of the software requires GNU compilers and MPI. Model output includes the thermal history of samples exhumed to the surface, cooling ages across topography, the transient thermal field, and the 3D position history of exhumed samples. Visualization of results is left up to the user, although outputs are provided in the Tecplot data format. A singularity (*.sif) container is also provided. Finally, example input files are included for a thrust belt-related simulation (from Eizenhoefer et al., 2023, Tectonics).