3D lithological structure in a steady state model drives divide migration E. Graf 1 , S. Mudd 1 , F. Kober 2 , A. Landgraf 2 , A. Ludwig 2 1 School of Geosciences, University of Edinburgh, UK 2 National Cooperative for the Disposal of Radioactive Waste (Nagra), Switzerland vEGU2021 – April 2021
24
Embed
3D lithological structure in a steady state model drives ...
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
3D lithological structure in a steady
state model drives divide migration
E. Graf1, S. Mudd1, F. Kober2, A. Landgraf2, A. Ludwig2
1 School of Geosciences, University of Edinburgh, UK2 National Cooperative for the Disposal of Radioactive Waste (Nagra), Switzerland
vEGU2021 – April 2021
In order to confidently model scenarios of future
topography, we want to account for variations in
lithology as well as potential changes in drainage
patterns (e.g. through river capture).
We present a model capable of incorporating 3D
geological data and adapting to alternative
drainage axes.
We demonstrate the model using a site in the
Swiss Jura Mountains. Values of erodibility K are
calibrated for different lithological units, and the
model is then run for selected incision and
alternative drainage scenarios.
Motivation
Right: Location of the study area in northern Switzerland (figure modified from Yanites et al., 2017). Black rectangle delineates the extent of the
study area according to the maps on subsequent slides.
The MuddPILE (Parsimonious Integrated Landscape Evolution)
Model (Mudd, 2017):
• calculates local relief using steady state solutions of the
stream power incision model (where gradient is related to
drainage area via a concavity index and a steepness index)
• quantifies hillslope relief using a very simple critical slope
gradient where hillslope angles are set to a critical value on
pixels that have a small drainage area
• allows drainage divides to migrate to minimize sharp breaks
in relief across them
For the 3D lithological structure, we use a 3D geological model
of the study area by Gmünder et al. (2013). See litho-
stratigraphic scheme at end of display.
Model overview
Right: model domain, draining into the fixed Rhine-Aare channel (blue).
A given combination of K values that achieves good fit is not necessarily the “correct” combination,
since multiple combinations can produce similar results in terms of goodness of fit (equifinality).
We pick the K combination resulting in the highest p-value and assign these values to the
corresponding erodibility classes in the 3D geological model (see end of display for detailed scheme,
including assignment of hillslope gradients).
Table: K combinations resulting in the 5 highest p-values. Shaded in red is the best-fit combination, which is used in the subsequent simulations
with 3D lithology.
Applications: Incision scenarios
3D UNIFORM
Final model topography using 3D lithology and uniform lithology, respectively, with 350 m main channel incision in both cases. White boxes indicate
characteristic areas of divide migration. Left: 3D lithology using the best-fit combination of K values from the Monte Carlo simulations. Right: Uniform
lithology, using K = 3.5E-06 and Sc = 0.21, which corresponds to Group 2 (Tertiary Sediments) in the best-fit combination. E = 0.0001 (implying
duration of incision = 3.5 My); m = 0.45; n = 1. Note that the uplift component is just indirectly simulated via the cumulative main channel incision,
partly resulting in hypothetical model elevations near or even below sea level.
Applications: Incision scenarios
3D UNIFORM
Lithologies exposed in the final model landscape using 3D lithology and uniform lithology, respectively, with 350 m main channel incision in
both cases. Left: 3D lithology using the best-fit combination of K values from the Monte Carlo simulations. Right: Uniform lithology, using K =
3.5E-06 and Sc = 0.21, which corresponds to Group 2 (Tertiary Sediments) in the best-fit combination. E = 0.0001 (implying duration of
incision = 3.5 My); m = 0.45; n = 1. See end of display for full legend of the 3D geological model..
Applications: Alternative drainage
We next simulate the incision scenario in
combination with the reactivation of a paleo-
channel, thus assessing the effect of lateral
changes of the main channel axis on local relief
and drainage divide migration.
Left: Base level channel for the alternative drainage scenario (blue).
The original course of the Aare is plotted in red. Digital terrain model:
Map legend, modified after Isler et al. (1984) – all units
Unit ID corresponds to violin plot (for bedrock units). Q = Quaternary unit; B = bedrock unit.
Map legend, mod. after Isler et al. (1984) – Quat. separation
Unit ID corresponds to violin plot. QD = Quaternary domain; B = bedrock unit; EC = erodibility class, corresponding to grouping of units for the Monte Carlo simulations. See map on slide 6 for spatial separation of QD1 and QD2.
Legend 3D Geol. Model (mod. from Gmünder et al. 2013)