1. General Model Information
Name: TREMP BASIN STRATIGRAPHIC SIMULATION MODEL
Main medium: aquatic+terrstrial
Main subject: hydrology,other
Organization level: Landscape
Type of model: cellular automat, difference equations (3D)
Main application: research, education
Keywords: foreland basin, thrusting, alluvial fan, longitudinal river, sea level fluctuation, incised valleys
Utrecht University, Faculty of Earth Sciences, Sedimentology Group
TREMP is a numerical model coded in C used to simulate the development of the Eocene Tremp foreland
basin in the Spanish Pyrenees. It is based on the GOLEM (Tucker) landscape evolution modelling code,
which is extended with:
Length Scales and Resolution Constraints:
- 3D flexure: is based on a fourier transform method
- Thrusting: A velocity description of deformation is used to advect bedrock and sediment over the rectangular grid. The model is capable of translating the foreland basin fill (and the stratigraphy !) over a shallow detachment fault, creating thrust-sheet top basins.
- Bedrock collapse: using the Cullman Slope stability criterion.
- Bifurcating flow: used to simulate alluvial fan cones and deltas filling the basin.
- Grainsize sorting: is performed using a combination of perfect sorting and the streampower law, Qs~KfQmSn, were Kf is made function of grainsize.
- Marine sedimentation: slope dependent submarine collapse, water depth dependent carbonate deposition
- Automated output routines are directly called from Tremp C-code. To Matlab (visualization stratigraphy fence-diagrams and subsurface property voxel cubes) and Surfer scripts (landscapes)
The model spatial dimensions are 75 x 75 km, using rectangular grid cells of 100-500 m width.
Vertical stratigraphic resolution is between the 0-1.5 m, or 1000 yrs of equivalent sedimentation.
Depending on memory available on your PC, every gridcell holds about 1000 stratigraphic layers.
Time Scales and Resolution Constraints:
Computational time steps are 10-50 yrs and the maximum simulation time worked with ~ 3 Myr.
Numerical limitations and issues:
Available system memory to store 3D stratigraphy is a limiting factor.
3D stratigraphy (channel belt voxels, fence-diagrams, synthetic wells), basin relief,
drainage patterns, process-specific erosion and sedimentation rates.
II. Technical Information
Operating System(s): MS Windows2000/XP, Linux
Programming Language(s): C/C++
Required Input: starting surface, fault positions, tectonic,sealevel and rainfall history.
III. Mathematical Information
Key Physical Parameters: exponents and coefficients in the streampower transport law.
Key Output: 3D stratigraphy (channel belt voxels, fence-diagrams, synthetic wells), basin relief, drainage patterns, process-specific erosion and sedimentation rates.
Clevis, Quintijn. (2003) Three-dimensional modelling of thrust-controlled foreland basin stratigraphy.
PhD Thesis, Utrecht University, Faculty of Earth Sciences. Geologica Ultraiectina, 226, pp. 136.
Clevis,Quintijn, De Boer,Poppe L. & Wachter,Maarten (2003) Numerical modelling of drainage basin evolution and three-dimensional alluvial fan stratigraphy.
Sedimentary Geology, 163(1-2), 85-110. doi:10.1016/S0037-0738(03)00174-X
Clevis,Quintijn, De Boer,Poppe L. & Nijman,Wouter (2004) Differentiating the effect of episodic tectonism and eustatic sea-level fluctuations in foreland basins filled by alluvial fans and axial deltaic systems: insights from a three-dimensional stratigraphic forward model.
Sedimentology 51(4), 809-835. doi: 10.1111/j.1365-3091.2004.00652.x
Clevis,Quintijn, de Jager,Gerben, Nijman,Wouter & de Boer,Poppe L. (2004) Stratigraphic signatures of translation of thrust-sheet top basins over low- angle detachment faults.
Basin Research 16(2), 145-163. doi: 10.1046/j.1365-2117.2003.00226.x
V. Further information in the World-Wide-Web
VI. Additional remarks
Last review of this document by: : Mon Jan 20 12:56:29 2003
Status of the document: Contributed by Quintijn Clevis
last modified by
Joachim Benz Thu Jul 1 10:46:45 CEST 2004