# Labs Landscape Evolution Modeling With Child Part 3: Difference between revisions

(9 intermediate revisions by the same user not shown) | |||

Line 1: | Line 1: | ||

==Landscape Evolution Modeling with CHILD== | ==Landscape Evolution Modeling with CHILD== | ||

< | <big>'''Part 3: Erosion and Transport by Running Water, Multiple Grain Sizes, and the Ten Commandments of Landscape Evolution Modeling'''</big><br> | ||

< | <br> | ||

< | <big>Gregory E. Tucker, University of Colorado Boulder, and Stephen T. Lancaster, Oregon State University</big><br> | ||

<br> | |||

These short course notes prepared for SIESD 2012: Future Earth: Interaction of Climate and Earth-surface Processes, University of Minnesota, MN, August 2012. <br> | |||

Notes and exercises updated for WMT by Stephanie Higgins, University of Colorado Boulder, August 2014 | |||

<br> | |||

<br> | |||

<div class="center">'''Before beginning these exercises, please download the CHILD visualization tools here: [[:File:CHILDVisTools.tar.gz]].'''</div> | |||

<br> | <br> | ||

<div class="center"><big>This material is Part 3 of a three-part introduction to landscape evolution modeling with CHILD in WMT. | <div class="center"><big>This material is Part 3 of a three-part introduction to landscape evolution modeling with CHILD in WMT. | ||

Part 1: [[Labs Landscape Evolution Modeling With Child Part 1]] | Part 1: [[Labs Landscape Evolution Modeling With Child Part 1]] | ||

Line 21: | Line 25: | ||

<p>where <span class="math">Φ <sub><em>c</em></sub></span> is, again, a threshold below which detachment is ineffective. Stream power per unit width turns out to be proportional to <span class="math"><em>τ</em><sup>3 / 2</sup></span>, so the two erosion formulas are closely related <span class="citation">(Whipple and Tucker 1999)</span>. In the following example, we will use the unit stream power formula with <span class="math">Φ <sub><em>c</em></sub> = 0</span>.</p> | <p>where <span class="math">Φ <sub><em>c</em></sub></span> is, again, a threshold below which detachment is ineffective. Stream power per unit width turns out to be proportional to <span class="math"><em>τ</em><sup>3 / 2</sup></span>, so the two erosion formulas are closely related <span class="citation">(Whipple and Tucker 1999)</span>. In the following example, we will use the unit stream power formula with <span class="math">Φ <sub><em>c</em></sub> = 0</span>.</p> | ||

<h2 id="exercise-6-detachment-limited-hills-and-mountains" class="unnumbered"><em>Exercise 6: Detachment-Limited Hills and Mountains</em></h2> | <h2 id="exercise-6-detachment-limited-hills-and-mountains" class="unnumbered"><em>Exercise 6: Detachment-Limited Hills and Mountains</em></h2> | ||

< | <p>>> In WMT, open the model titled "dlim." Save and submit the job.</p> | ||

<p>>> The 3 m.y. run should take about 20 seconds. Download the results and move them to the folder ''childFolderLocation''∖ChildExercises∖dlim. For ''childFolderLocation'', use the path to the Child visualization tools package that you downloaded.</p> | |||

<p>>> In Matlab, navigate to the dlim folder. Type:</p> | |||

<ul> | |||

<span><p><li>figure(1), clf, colormap jet</li></p></span> | |||

<span><p><li>cmovie( 'dlim', 31, 3e4, 3e4, 1e3, 500 );</li></p></span> | |||

<span><p><li>figure(2), clf</li></p></span> | |||

<span><p><li>csa( 'dlim', 31 ); % Shows slope-area graph</li></p></span> | |||

</ul> | |||

<span><p>Notice that the landscape has come close to a state of equilibrium between erosion and relative uplift. The resulting terrain has about 200 m of relief over a 30 km half-width mountain range---more Appalachian than Himalayan. Notice that the log-log slope-area graph shows a straight line, indicating a power-law relationship. This is exactly to be expected, and we can predict the plot slope and intercept analytically. Finally, note the points on the upper left of the graph. These “first order” cells, at about 2500 m<math>^2</math> contributing area, have slopes less than 10%. They represent embedded channels, not hillslopes, which are too small to resolve at this grid spacing.</p></span> | |||

<span><p>Now, what happens when we increase the relative uplift rate?</p></span> | |||

<p>>> In WMT, open the model titled "dlimC1." Save and submit the job.</p> | |||

<p>>> Download the results and move them to the folder ''childFolderLocation''∖ChildExercises∖dlim. For ''childFolderLocation'', use the path to the Child visualization tools package that you downloaded.</p> | |||

<span><p>This run starts off where the previous one ended, but with a 10x higher rate of relative uplift.</p></span> | |||

<p>>> In Matlab, type:</p> | |||

<ul> | |||

<span><p><li>figure(1)</li></p></span> | |||

<span><p><li>cmovie( 'dlimC1', 31, 3e4, 3e4, 1e4, 5000 ); % 10x vertical scale</li></p></span> | |||

<span><p><li>figure(2)</li></p></span> | |||

<span><p><li>hold on, csa( 'dlimC1', 31, 'r.' ); hold off</li></p></span> | |||

</ul> | |||

<span><p>Because we are using a slope-linear detachment law, a 10x increase in relative uplift rate leads to a 10x increase in relief. Notice that the points have shifted upward by a factor of 10 on the slope-area graph.</p></span> | |||

<span><p>We still do not see any hillslopes, because the scale of landscape dissection is too fine for the model to resolve. </p></span> | |||

<h2 id="exercise-7-zooming-in-to-the-hillslopes" class="unnumbered"><em>Exercise 7: Zooming in to the Hillslopes</em></h2> | <h2 id="exercise-7-zooming-in-to-the-hillslopes" class="unnumbered"><em>Exercise 7: Zooming in to the Hillslopes</em></h2> | ||

< | Next, we will “zoom in” by repeating the dlim run but with a twenty-fold decrease in domain size and model cell size. | ||

<p>>> In WMT, open the model titled "dlim_small." Save and submit the job.</p> | |||

<p>>> Download the results and move them to the folder ''childFolderLocation''∖ChildExercises∖dlim. For ''childFolderLocation'', use the path to the Child visualization tools package that you downloaded.</p> | |||

This run is identical to dlim but with a domain of 1.5 by 1.5km and <math>\sim</math>25m wide cells, instead of 30x30km and <math>\sim</math>500m cells. | |||

In Matlab, type: | |||

<ul> | |||

<span><p><li>figure(1)</li></p></span> | |||

<span><p><li>cmovie( 'dlim\_small', 31, 1.5e3, 1.5e3, 500, 200 );</li></p></span> | |||

<span><p><li>figure(2)</li></p></span> | |||

<span><p><li>hold on, csa( 'dlim\_small', 31, 'g.' ); hold off</li></p></span> | |||

</ul> | |||

Note how the hillslopes become evident in the topography. In the slope-area plot, the points seem to continue the trend of the coarser-scale run, but somewhat shifted upward. Can you guess why they are shifted upward? (The answer is subtle, and lies hidden in dlim_small2). | |||

<h2 id="exercise-8-knickzones-and-transient-response" class="unnumbered"><em>Exercise 8: Knickzones and Transient Response</em></h2> | <h2 id="exercise-8-knickzones-and-transient-response" class="unnumbered"><em>Exercise 8: Knickzones and Transient Response</em></h2> | ||

< | For the next exercise, we return to our earlier dlimC1 run and plot a representative stream profile at different times, to look at how the profile responds to the increased rate of relative uplift. | ||

<p>>> In Matlab, type: </p> | |||

<ul> | |||

<span><p><li>figure(1), clf</li></p></span> | |||

<span><p><li>[d,h,x,y] = cstrmproseries( 'dlimC1', 10, 15000, 29000 );</li></p></span> | |||

</ul> | |||

<span><p>This command traces the stream profile starting from <span><em>x=15</em></span> km, <span><em>y=29</em></span> km. It will plot the first 10 profiles.</p></span> | |||

<ul> | |||

<span><p><li>figure(2), clf, plot( x, y )</li></p></span> | |||

</ul> | |||

<span><p>This shows the horizontal trace of the stream course.</p></span> | |||

<span><p>During the period of transient response, the stream profile shows a pronounced convexity, or knickzone, along the profile. The knickzone marches upstream through time. This pattern is characteristic of the “stream power” erosion law, which is actually a form of wave equation.</p></span> | |||

<h2 id="transport-limited-models">Transport-Limited Models</h2> | <h2 id="transport-limited-models">Transport-Limited Models</h2> | ||

<p>We next explore the dynamics of landscapes and networks with transport-limited models. One caution as we do so: we will assume that channel width is independent of grain size, slope, etc.</p> | <p>We next explore the dynamics of landscapes and networks with transport-limited models. One caution as we do so: we will assume that channel width is independent of grain size, slope, etc.</p> | ||

<h2 id="exercise-9-a-pile-of-fine-sand" class="unnumbered"><em>Exercise 9: A Pile of Fine Sand</em></h2> | <h2 id="exercise-9-a-pile-of-fine-sand" class="unnumbered"><em>Exercise 9: A Pile of Fine Sand</em></h2> | ||

< | <p>>> In WMT, open the model titled "tlim1." Save and submit the job.</p> | ||

<p>>> The 1 m.y. run should take about 2 minutes. Download the results and move them to the folder ''childFolderLocation''∖ChildExercises∖tlim. For ''childFolderLocation'', use the path to the Child visualization tools package that you downloaded.</p> | |||

<p>>> In Matlab, type:</p> | |||

<ul> | |||

<span><p><li>cd('../tlim')</li></p></span> | |||

<span><p><li>figure(1), clf</li></p></span> | |||

<span><p><li>cmovie( 'tlim1', 21, 3e4, 3e4, 40, 10 );</li></p></span> | |||

<span><p><li>figure(2), clf</li></p></span> | |||

<span><p><li>csa( 'tlim1', 21 ); axis([1e-1 1e3 1e-4 1e-3])</li></p></span> | |||

</ul> | |||

<span><p>In this run, we are effectively assuming that 0.1 mm sand moves as bed-load, according to a Meyer-Peter and Mueller-like transport formula. The landscape takes on an effectively uniform and very shallow gradient, on the order of <math>3\times 10^{-4}</math>.</p></span> | |||

<h2 id="exercise-10-a-pile-of-cobbles" class="unnumbered"><em>Exercise 10: A Pile of Cobbles</em></h2> | <h2 id="exercise-10-a-pile-of-cobbles" class="unnumbered"><em>Exercise 10: A Pile of Cobbles</em></h2> | ||

< | Now let's try the same experiment with 5cm cobbles. | ||

<p>>> In WMT, open the model titled "tlim2" Save and submit the job.</p> | |||

<p>>> The 3 m.y. run should take about 2-3 minutes. Download the results and move them to the folder ''childFolderLocation''∖ChildExercises∖tlim. For ''childFolderLocation'', use the path to the Child visualization tools package that you downloaded.</p> | |||

<p>>> In Matlab, type:</p> | |||

<ul> | |||

<span><p><li>figure(1), clf</li></p></span> | |||

<span><p><li>cmovie( 'tlim2', 31, 3e4, 3e4, 1000, 300 );</li></p></span> | |||

<span><p><li>figure(2)</li></p></span> | |||

<span><p><li>hold on, csa( 'tlim2', 31, 'r.' ); hold off</li></p></span> | |||

<span><p><li>axis([1e-1 1e3 1e-4 1e-1])</li></p></span> | |||

</ul> | |||

<span><p>Lesson: grain size matters! </p></span> | |||

<span><p>But let's remember the caveat that channel width matters too, and we haven't taken that into account with these simple runs. Also, Nicole Gasparini's work (Gasparini et al., 1999; Gasparini et al., 2004) tells us that channel concavity is less sensitive to grain size when there is a mixture of sizes available to the river. </p></span> | |||

<span><p>Optional exercise: Reconfigure tlim2 to have a higher uplift rate. Save the model with a new name, run and download the results. Use the Matlab script <tt>cstrmproseries</tt> to plot fluvial profiles undergoing transient response. How do these compare with the detachment-limited model?</p></span> | |||

<h2 id="hybrid-model-combining-detachment-and-transport">Hybrid Model: Combining Detachment and Transport</h2> | <h2 id="hybrid-model-combining-detachment-and-transport">Hybrid Model: Combining Detachment and Transport</h2> | ||

<p>Next, we’ll look at a more complex situation with simultaneous erosion and sedimentation, and simultaneous detachment-limited and transport-limited behavior. In this case, we use a fluvial model in which erosion rate can be limited either by transport capacity or by detachment capacity, depending on their relative magnitudes:</p> | <p>Next, we’ll look at a more complex situation with simultaneous erosion and sedimentation, and simultaneous detachment-limited and transport-limited behavior. In this case, we use a fluvial model in which erosion rate can be limited either by transport capacity or by detachment capacity, depending on their relative magnitudes:</p> | ||

Line 53: | Line 124: | ||

\end{cases}</math> (19)</center></span><br /></p> | \end{cases}</math> (19)</center></span><br /></p> | ||

<h2 id="exercise-11-erosion-and-deposition-together-at-last" class="unnumbered"><em>Exercise 11: Erosion and Deposition, Together at Last</em></h2> | <h2 id="exercise-11-erosion-and-deposition-together-at-last" class="unnumbered"><em>Exercise 11: Erosion and Deposition, Together at Last</em></h2> | ||

< | <p>>> In WMT, open the model titled "erodep1.” Save and submit the job.</p> | ||

<p>>> The 1 m.y. run should take about 5 minutes. Download the results and move them to the folder ''childFolderLocation''∖ChildExercises∖hyrbid. For ''childFolderLocation'', use the path to the Child visualization tools package that you downloaded.</p> | |||

<p>>> In Matlab, navigate to the hybrid folder. Type:</p> | |||

<ul> | |||

<span><p><li>figure(1), clf</li></p></span> | |||

<span><p><li>cmovie( 'erodep1', 21, 6e4, 6e4, 4000 );</li></p></span> | |||

</ul> | |||

<span><p>Here we have a block rising at 1 mm/yr and an adjacent block subsiding at 0.25 mm/yr. Uplift and subsidence shut down after 500 ky. The subsiding block forms a large lake that gradually fills in with fan-deltas.</p></span> | |||

<h2 id="other-sediment-flux-dependent-fluvial-models">Other Sediment-Flux-Dependent Fluvial Models</h2> | <h2 id="other-sediment-flux-dependent-fluvial-models">Other Sediment-Flux-Dependent Fluvial Models</h2> | ||

<p>We won’t take the time to address some of the other models, including</p> | <p>We won’t take the time to address some of the other models, including</p> | ||

Line 97: | Line 174: | ||

<li><p>If thou darest to predict future erosion, thou shalt heed the previous commandment ten times over (but thou mightest point out to skeptics that a process-based prediction is usually better than one based on pure guesswork, provided that commandment #2 is obeyed).</p></li> | <li><p>If thou darest to predict future erosion, thou shalt heed the previous commandment ten times over (but thou mightest point out to skeptics that a process-based prediction is usually better than one based on pure guesswork, provided that commandment #2 is obeyed).</p></li> | ||

</ol> | </ol> | ||

<h1 id="references">References</h1> | <h1 id="references">References</h1> | ||

Line 108: | Line 186: | ||

<p>Chase, Clement G. 1992. “Fluvial Landsculpting and the Fractal Dimension of Topography.” <em>Geomorphology</em> 5: 39–57.</p> | <p>Chase, Clement G. 1992. “Fluvial Landsculpting and the Fractal Dimension of Topography.” <em>Geomorphology</em> 5: 39–57.</p> | ||

<p>Clevis, Q., G. Tucker, S. Lancaster, A. Desitter, N. Gasparini, and G. Lock. 2006. “Geoarchaeological Simulation of Meandering River Deposits and Settlement Distributions: a Three-Dimensional Approach.” <em>Computers and Geosciences</em> 21 (8): 843–874.</p> | <p>Clevis, Q., G. Tucker, S. Lancaster, A. Desitter, N. Gasparini, and G. Lock. 2006. “Geoarchaeological Simulation of Meandering River Deposits and Settlement Distributions: a Three-Dimensional Approach.” <em>Computers and Geosciences</em> 21 (8): 843–874.</p> | ||

<p>Collins, DBG, RL Bras, and GE Tucker. 2004. “Modeling the Effects of Vegetation-Erosion Coupling on Landscape Evolution.” <em>Journal of Geophysical Research—Earth Surface</em> 109 (F3) ({AUG 3}). doi: | <p>Collins, DBG, RL Bras, and GE Tucker. 2004. “Modeling the Effects of Vegetation-Erosion Coupling on Landscape Evolution.” <em>Journal of Geophysical Research—Earth Surface</em> 109 (F3) ({AUG 3}). doi:10.1029/2003JF000028. | ||

<p>Coulthard, TJ, MJ Kirkby, and MG Macklin. 1996. “A Cellular Automaton Landscape Evolution Model.” <em>Proceedings of the First International Conference on GeoComputation</em>.</p> | <p>Coulthard, TJ, MJ Kirkby, and MG Macklin. 1996. “A Cellular Automaton Landscape Evolution Model.” <em>Proceedings of the First International Conference on GeoComputation</em>.</p> | ||

<p>Crave, A, and Ph Davy. 2001. “A Stochastic ’Precipiton’ Model for Simulating Erosion/Sedimentation Dynamics.” <em>Computers and Geosciences</em> 27: 815–827.</p> | <p>Crave, A, and Ph Davy. 2001. “A Stochastic ’Precipiton’ Model for Simulating Erosion/Sedimentation Dynamics.” <em>Computers and Geosciences</em> 27: 815–827.</p> | ||

<p>Crosby, Benjamin T., Kelin X. Whipple, Nicole M. Gasparini, and Cameron W. Wobus. 2007. “Formation of Fluvial Hanging Valleys: Theory and Simulation.” <em>JOURNAL OF GEOPHYSICAL RESEARCH-EARTH SURFACE</em> 112 (F3) ({AUG 9}). doi: | <p>Crosby, Benjamin T., Kelin X. Whipple, Nicole M. Gasparini, and Cameron W. Wobus. 2007. “Formation of Fluvial Hanging Valleys: Theory and Simulation.” <em>JOURNAL OF GEOPHYSICAL RESEARCH-EARTH SURFACE</em> 112 (F3) ({AUG 9}). doi:10.1029/2006JF000566.</p> | ||

<p>Culling, WEH. 1963. “Soil Creep and the Development of Hillside Slopes.” <em>The Journal of Geology</em>: 127–161.</p> | <p>Culling, WEH. 1963. “Soil Creep and the Development of Hillside Slopes.” <em>The Journal of Geology</em>: 127–161.</p> | ||

<p>Densmore, Alexander L, Michael A Ellis, and Robert S Anderson. 1998. “Landsliding and the Evolution of Normal-Fault-Bounded Mountains.” <em>Journal of Geophysical Research</em> 103: 15203–15219.</p> | <p>Densmore, Alexander L, Michael A Ellis, and Robert S Anderson. 1998. “Landsliding and the Evolution of Normal-Fault-Bounded Mountains.” <em>Journal of Geophysical Research</em> 103: 15203–15219.</p> | ||

Line 119: | Line 197: | ||

<p>Flores-Cervantes, Homero, Erkan Istanbulluoglu, and Rafael L Bras. 2006. “Development of Gullies on the Landscape: a Model of Headcut Retreat Resulting from Plunge Pool Erosion.” <em>Journal of Geophysical Research</em> 111: F01010.</p> | <p>Flores-Cervantes, Homero, Erkan Istanbulluoglu, and Rafael L Bras. 2006. “Development of Gullies on the Landscape: a Model of Headcut Retreat Resulting from Plunge Pool Erosion.” <em>Journal of Geophysical Research</em> 111: F01010.</p> | ||

<p>Garcia-Castellanos, D. 2002. “Interplay Between Lithospheric Flexure and River Transport in Foreland Basins.” <em>Basin Research</em> 14 (2): 89–104.</p> | <p>Garcia-Castellanos, D. 2002. “Interplay Between Lithospheric Flexure and River Transport in Foreland Basins.” <em>Basin Research</em> 14 (2): 89–104.</p> | ||

<p>Gasparini, N. M., K. X. Whipple, and R. L. Bras. 2007. “Predictions of Steady State and Transient Landscape Morphology Using Sediment-Flux-Dependent River Incision Models.” <em> | <p>Gasparini, N., G. Tucker, and R. Bras. 1999. "Downstream fining through selective particle sorting in an equilibrium drainage network." <em>Geology</em> 27(12), 1079.</p> | ||

<p>Gasparini, N., G. Tucker, and R. Bras. 2004. "Network-scale dynamics of grain-size sorting: Implications for downstream fining, stream-profile concavity, and drainage basin morphology. <em>Earth Surface Processes and Landforms</em> 29(4), 401-421.</p> | |||

<p>Gasparini, N. M., K. X. Whipple, and R. L. Bras. 2007. “Predictions of Steady State and Transient Landscape Morphology Using Sediment-Flux-Dependent River Incision Models.” <em>Journal of Geophysical Research Earth Surface</em> 112 (F3) ({MAY 19}). doi: 10.1029/2006JF000567.</p> | |||

<p>Gilbert, GK. 1877. “Report on the Geology of the Henry Mountains: US Geog. and Geol.” <em>Survey, Rocky Mtn. Region</em> 160.</p> | <p>Gilbert, GK. 1877. “Report on the Geology of the Henry Mountains: US Geog. and Geol.” <em>Survey, Rocky Mtn. Region</em> 160.</p> | ||

<p>Herman, F., and J. Braun. 2006. “Fluvial Response to Horizontal Shortening and Glaciations: a Study in the Southern Alps of New Zealand.” <em>Journal of Geophysical Research-Earth Surface</em> 111 (F1): F01008.</p> | <p>Herman, F., and J. Braun. 2006. “Fluvial Response to Horizontal Shortening and Glaciations: a Study in the Southern Alps of New Zealand.” <em>Journal of Geophysical Research-Earth Surface</em> 111 (F1): F01008.</p> | ||

Line 140: | Line 220: | ||

<p>Sólyom, P.B., and G.E. Tucker. 2004. “Effect of Limited Storm Duration on Landscape Evolution, Drainage Basin Geometry, and Hydrograph Shapes.” <em>Journal of Geophysical Research</em> 109: 13.</p> | <p>Sólyom, P.B., and G.E. Tucker. 2004. “Effect of Limited Storm Duration on Landscape Evolution, Drainage Basin Geometry, and Hydrograph Shapes.” <em>Journal of Geophysical Research</em> 109: 13.</p> | ||

<p>Tucker, G E, and R L Slingerland. 1994. “Erosional Dynamics, Flexural Isostasy, and Long-Lived Escarpments: a Numerical Modeling Study.” <em>Journal of Geophysical Research</em> 99: 12, 229–12, 243.</p> | <p>Tucker, G E, and R L Slingerland. 1994. “Erosional Dynamics, Flexural Isostasy, and Long-Lived Escarpments: a Numerical Modeling Study.” <em>Journal of Geophysical Research</em> 99: 12, 229–12, 243.</p> | ||

<p>Tucker, GE. 2004. “Drainage Basin Sensitivity to Tectonic and Climatic Forcing: Implications of a Stochastic Model for the Role of Entrainment and Erosion Thresholds.” <em>Earth Surface Processes and Landforms</em> 29 (2) ({FEB}): 185–205. doi: | <p>Tucker, GE. 2004. “Drainage Basin Sensitivity to Tectonic and Climatic Forcing: Implications of a Stochastic Model for the Role of Entrainment and Erosion Thresholds.” <em>Earth Surface Processes and Landforms</em> 29 (2) ({FEB}): 185–205. doi:10.1002/esp.1020.</p> | ||

<p>Tucker, Gregory E, and Rafael L Bras. 2000. “A Stochastic Approach to Modeling the Role of Rainfall Variability in Drainage Basin Evolution.” <em>Water Resources Research</em> 36 (7): 1953–1964.</p> | <p>Tucker, Gregory E, and Rafael L Bras. 2000. “A Stochastic Approach to Modeling the Role of Rainfall Variability in Drainage Basin Evolution.” <em>Water Resources Research</em> 36 (7): 1953–1964.</p> | ||

<p>Tucker, Gregory E, and Gregory R Hancock. 2010. “Modelling Landscape Evolution.” <em>Earth Surface Processes and Landforms</em> 46: 28–50.</p> | <p>Tucker, Gregory E, and Gregory R Hancock. 2010. “Modelling Landscape Evolution.” <em>Earth Surface Processes and Landforms</em> 46: 28–50.</p> |

## Latest revision as of 16:45, 24 September 2014

## Landscape Evolution Modeling with CHILD

**Part 3: Erosion and Transport by Running Water, Multiple Grain Sizes, and the Ten Commandments of Landscape Evolution Modeling**

Gregory E. Tucker, University of Colorado Boulder, and Stephen T. Lancaster, Oregon State University

These short course notes prepared for SIESD 2012: Future Earth: Interaction of Climate and Earth-surface Processes, University of Minnesota, MN, August 2012.

Notes and exercises updated for WMT by Stephanie Higgins, University of Colorado Boulder, August 2014

**Before beginning these exercises, please download the CHILD visualization tools here: File:CHILDVisTools.tar.gz.**

Part 1: Labs Landscape Evolution Modeling With Child Part 1

Part 2: Labs Landscape Evolution Modeling With Child Part 2# Erosion and Transport by Running Water

There are several competing models for erosion by channelized flow. Detachment-limited models assume that eroded material leaves the system without significant re-deposition and that lowering of channels is limited by the ability of the stream to detach material from the bed (Howard 1994; Whipple and Tucker 1999). Transport-limited models assume plentiful supply of loose sediment and that lowering of channels is limited by the stream’s capacity to transport sediment (Willgoose, Bras, and Rodriguez-Iturbe 1991; Whipple and Tucker 2002). In simple hybrid models, lowering may be limited either by excess transport capacity or by detachment rate, depending on local sediment supply and substrate resistance (Gregory E Tucker, Lancaster, Gasparini, and Bras 2001; Whipple and Tucker 2002). With the undercapacity concept, detachment rate depends on surplus transport capacity (Beaumont, Fullsack, and Hamilton 1992). In the saltation-abrasion model, detachment is driven by grain impacts and limited by sediment shielding (Gasparini, Whipple, and Bras, 2007; Whipple and Tucker, 2002.

## Detachment-Limited Models

On a cohesive or rock bed with a discontinuous or absent cover of loose sediment, detachment of particles from the bed may be driven primarily by hydraulic lift and drag (“plucking”). Most models assume that the rate of detachment (or more generally the capacity for detachment) depends on excess bed shear stress:

*D*

_{c}=

*K*

_{b}(

*τ*−

*τ*

_{c})

^{pb}, or alternatively,

*D*

_{c}=

*K*

_{b}(

*τ*

^{pb}−

*τ*

_{c}

^{pb}) (17)

where *τ* is local bed shear stress, *τ*_{c} is a threshold stress below which detachment is ineffective, *K*_{b} is a constant, and *p*_{b} is an exponent.

Bed shear stress fluctuates in space and time, but is often treated using the cross-sectional average, which in turn is based on a force balance between gravity and friction.

Some models assume that the detachment rate depends on stream power per unit width, *ω* = *ρ**g*(*Q* / *W*)*S*:

where Φ _{c} is, again, a threshold below which detachment is ineffective. Stream power per unit width turns out to be proportional to *τ*^{3 / 2}, so the two erosion formulas are closely related (Whipple and Tucker 1999). In the following example, we will use the unit stream power formula with Φ _{c} = 0.

*Exercise 6: Detachment-Limited Hills and Mountains*

>> In WMT, open the model titled "dlim." Save and submit the job.

>> The 3 m.y. run should take about 20 seconds. Download the results and move them to the folder *childFolderLocation*∖ChildExercises∖dlim. For *childFolderLocation*, use the path to the Child visualization tools package that you downloaded.

>> In Matlab, navigate to the dlim folder. Type:

- figure(1), clf, colormap jet
- cmovie( 'dlim', 31, 3e4, 3e4, 1e3, 500 );
- figure(2), clf
- csa( 'dlim', 31 ); % Shows slope-area graph

Notice that the landscape has come close to a state of equilibrium between erosion and relative uplift. The resulting terrain has about 200 m of relief over a 30 km half-width mountain range---more Appalachian than Himalayan. Notice that the log-log slope-area graph shows a straight line, indicating a power-law relationship. This is exactly to be expected, and we can predict the plot slope and intercept analytically. Finally, note the points on the upper left of the graph. These “first order” cells, at about 2500 m[math]\displaystyle{ ^2 }[/math] contributing area, have slopes less than 10%. They represent embedded channels, not hillslopes, which are too small to resolve at this grid spacing.

Now, what happens when we increase the relative uplift rate?

>> In WMT, open the model titled "dlimC1." Save and submit the job.

>> Download the results and move them to the folder *childFolderLocation*∖ChildExercises∖dlim. For *childFolderLocation*, use the path to the Child visualization tools package that you downloaded.

This run starts off where the previous one ended, but with a 10x higher rate of relative uplift.

>> In Matlab, type:

- figure(1)
- cmovie( 'dlimC1', 31, 3e4, 3e4, 1e4, 5000 ); % 10x vertical scale
- figure(2)
- hold on, csa( 'dlimC1', 31, 'r.' ); hold off

Because we are using a slope-linear detachment law, a 10x increase in relative uplift rate leads to a 10x increase in relief. Notice that the points have shifted upward by a factor of 10 on the slope-area graph.

We still do not see any hillslopes, because the scale of landscape dissection is too fine for the model to resolve.

*Exercise 7: Zooming in to the Hillslopes*

Next, we will “zoom in” by repeating the dlim run but with a twenty-fold decrease in domain size and model cell size.

>> In WMT, open the model titled "dlim_small." Save and submit the job.

>> Download the results and move them to the folder *childFolderLocation*∖ChildExercises∖dlim. For *childFolderLocation*, use the path to the Child visualization tools package that you downloaded.

This run is identical to dlim but with a domain of 1.5 by 1.5km and [math]\displaystyle{ \sim }[/math]25m wide cells, instead of 30x30km and [math]\displaystyle{ \sim }[/math]500m cells. In Matlab, type:

- figure(1)
- cmovie( 'dlim\_small', 31, 1.5e3, 1.5e3, 500, 200 );
- figure(2)
- hold on, csa( 'dlim\_small', 31, 'g.' ); hold off

Note how the hillslopes become evident in the topography. In the slope-area plot, the points seem to continue the trend of the coarser-scale run, but somewhat shifted upward. Can you guess why they are shifted upward? (The answer is subtle, and lies hidden in dlim_small2).

*Exercise 8: Knickzones and Transient Response*

For the next exercise, we return to our earlier dlimC1 run and plot a representative stream profile at different times, to look at how the profile responds to the increased rate of relative uplift.

>> In Matlab, type:

- figure(1), clf
- [d,h,x,y] = cstrmproseries( 'dlimC1', 10, 15000, 29000 );

This command traces the stream profile starting from *x=15* km, *y=29* km. It will plot the first 10 profiles.

- figure(2), clf, plot( x, y )

This shows the horizontal trace of the stream course.

During the period of transient response, the stream profile shows a pronounced convexity, or knickzone, along the profile. The knickzone marches upstream through time. This pattern is characteristic of the “stream power” erosion law, which is actually a form of wave equation.

## Transport-Limited Models

We next explore the dynamics of landscapes and networks with transport-limited models. One caution as we do so: we will assume that channel width is independent of grain size, slope, etc.

*Exercise 9: A Pile of Fine Sand*

>> In WMT, open the model titled "tlim1." Save and submit the job.

>> The 1 m.y. run should take about 2 minutes. Download the results and move them to the folder *childFolderLocation*∖ChildExercises∖tlim. For *childFolderLocation*, use the path to the Child visualization tools package that you downloaded.

>> In Matlab, type:

- cd('../tlim')
- figure(1), clf
- cmovie( 'tlim1', 21, 3e4, 3e4, 40, 10 );
- figure(2), clf
- csa( 'tlim1', 21 ); axis([1e-1 1e3 1e-4 1e-3])

In this run, we are effectively assuming that 0.1 mm sand moves as bed-load, according to a Meyer-Peter and Mueller-like transport formula. The landscape takes on an effectively uniform and very shallow gradient, on the order of [math]\displaystyle{ 3\times 10^{-4} }[/math].

*Exercise 10: A Pile of Cobbles*

Now let's try the same experiment with 5cm cobbles.

>> In WMT, open the model titled "tlim2" Save and submit the job.

>> The 3 m.y. run should take about 2-3 minutes. Download the results and move them to the folder *childFolderLocation*∖ChildExercises∖tlim. For *childFolderLocation*, use the path to the Child visualization tools package that you downloaded.

>> In Matlab, type:

- figure(1), clf
- cmovie( 'tlim2', 31, 3e4, 3e4, 1000, 300 );
- figure(2)
- hold on, csa( 'tlim2', 31, 'r.' ); hold off
- axis([1e-1 1e3 1e-4 1e-1])

Lesson: grain size matters!

But let's remember the caveat that channel width matters too, and we haven't taken that into account with these simple runs. Also, Nicole Gasparini's work (Gasparini et al., 1999; Gasparini et al., 2004) tells us that channel concavity is less sensitive to grain size when there is a mixture of sizes available to the river.

Optional exercise: Reconfigure tlim2 to have a higher uplift rate. Save the model with a new name, run and download the results. Use the Matlab script `cstrmproseries` to plot fluvial profiles undergoing transient response. How do these compare with the detachment-limited model?

## Hybrid Model: Combining Detachment and Transport

Next, we’ll look at a more complex situation with simultaneous erosion and sedimentation, and simultaneous detachment-limited and transport-limited behavior. In this case, we use a fluvial model in which erosion rate can be limited either by transport capacity or by detachment capacity, depending on their relative magnitudes:

*Exercise 11: Erosion and Deposition, Together at Last*

>> In WMT, open the model titled "erodep1.” Save and submit the job.

>> The 1 m.y. run should take about 5 minutes. Download the results and move them to the folder *childFolderLocation*∖ChildExercises∖hyrbid. For *childFolderLocation*, use the path to the Child visualization tools package that you downloaded.

>> In Matlab, navigate to the hybrid folder. Type:

- figure(1), clf
- cmovie( 'erodep1', 21, 6e4, 6e4, 4000 );

Here we have a block rising at 1 mm/yr and an adjacent block subsiding at 0.25 mm/yr. Uplift and subsidence shut down after 500 ky. The subsiding block forms a large lake that gradually fills in with fan-deltas.

## Other Sediment-Flux-Dependent Fluvial Models

We won’t take the time to address some of the other models, including

“Under-capacity” models (detachment rate depends on degree to which sediment flux falls below transport capacity), and

Saltation-abrasion models (detachment rate driven by particle impacts, and limited by alluvial shielding of bed)

Gasparini, Whipple, and Bras (2007) explore the behavior of these models with CHILD simulations.

# Multiple Grain Sizes

Although we won’t explore the effects of including multiple grain sizes of sediment in transport, grain size introduces some interesting issues, including:

Bed armoring and its impact on transport rates

Downstream fining

Abrasion and lithologic controls

# Exotica

Landscape evolution models include more than diffusion and stream-power models:

Stream meandering in the context of landscape evolution and valley stratigraphy (Clevis et al. 2006 a,b).

Vegetation, including both grass (Collins, Bras, and Tucker 2004; Istanbulluoglu and Bras 2005) and trees (Lancaster, Hayes, and Grant 2003)

Alternate forms of mass wasting, including landslides and debris flows (Densmore, Ellis, and Anderson 1998; Lancaster, Hayes, and Grant 2003; Istanbulluoglu et al. 2005)

Knickpoints, hanging valleys, and plunge pools (Flores-Cervantes, Istanbulluoglu, and Bras 2006; Crosby et al. 2007)

Glaciation (Herman and Braun 2006; Herman, Braun, and Dunlap 2007; Herman and Braun 2008)

# Forecasting or Speculation?

Some mathematical models in the physical sciences have such firm foundations that they can be relied upon to forecast the behavior of the natural world. For example, laws of motion of objects in a vacuum are absolutely reliable (as long as their speed is much less than that of light). The same can be said for numerical solutions to these equations, provided the solution is reasonably accurate. For these kinds of model, the verb “to model” means to calculate with high reliability what would happen under a particular set of initial and boundary conditions.

At the other end of the spectrum, we have mathematical models that are essentially tentative hypotheses. Such models are often based on intuition about a physical system, and represent a sort of educated guess about the quantitative relationships between things. For example, when Ahnert (1976) presented his inverse-exponential equation for regolith generation from bedrock, he was essentially expressing a conceptual hypothesis in mathematical terms. For these models-as-hypotheses, the phrase “to model” means to perform a quantitative “what if” experiment, asking the question: what kinds of pattern would I see if my hypothesis were correct? Comparing the prediction with observations provides a test of the hypothesis.

One can find many models that fall between these extremes. There are models that are based on well-known physics, but which are forced to use approximations of unknown accuracy in order to solve the governing equations. For example, climate models typically use simple parameterization schemes to represent convective mass and energy transport. Then too there are models that combine basic physical principles with elements of intuition, empiricism, and approximation. Arguably, many sediment-transport laws fall into this category: they are based on firm mechanical foundations (the force balance on a sediment grain) but also rely on strong approximations of factors like grain geometry, local flow velocity, and so on.

By now, it should be obvious that landscape evolution models also fall somewhere between the end-member cases of “model as truth” and “model as speculative hypothesis.” As we have seen throughout this course, there is a varying degree of experimental and observational support for the individual transport, weathering and erosion laws that go into a typical landscape model. In that sense, then, these models amount to more than just speculation. But equally there is still an element of speculation behind many of the process laws used in landscape models. Also, the process laws and algorithms represent a significant amount of upscaling in space and (especially) time. For example, the use of a steady precipitation rate as a proxy for the natural sequence of flows in a river channel represents a major approximation. For these reasons, we believe that three of the most important frontiers in landscape evolution research are (1) continuing to test individual process laws in the field and lab, (2) testing whole-landscape models using natural experiments, and (3) using mathematics, computation and experiments to study how the rates of various processes scale upward in time and space, and how these can be effectively parameterized.

# Ten Commandments of Landscape Evolution Modeling

Thou shalt not use a model without understanding the ingredients therein.

Be thou ever mindful of uncertainty.

Thou shalt use thy model to develop insight.

Thou shalt take delight when thy model surprises thee.

Thou shalt kick thy model hard, that it may notice thee (an injunction borrowed gratefully from the 10 Climate Modeling Commandments).

Thou shalt diagnose the reasons for thy model’s behavior.

Thou shalt conduct sensitivity experiments and “play around.”

Thou shalt use thy model to discover the necessary and sufficient conditions needed to explain thy target problem.

If thou darest use a model to calculate what happened in your field area in the past, thou shalt find a way to test and calibrate it first.

If thou darest to predict future erosion, thou shalt heed the previous commandment ten times over (but thou mightest point out to skeptics that a process-based prediction is usually better than one based on pure guesswork, provided that commandment #2 is obeyed).

# References

Ahnert, F. 1976. “Brief Description of a Comprehensive Three-Dimensional Process-Response Model of Landform Development.” *Zeitschrift FÜR Geomorfologie, Supplementband* 25: 29–49.

Ahnert, F. 1971. “Brief Description of a Comprehensive Three-Dimensional Process-Response Model of Landform Development.” *Zeitschrift Fur Geomorfologie, Supplementband* 25: 29–49.

Amos, Colin B, and Douglas W Burbank. 2007. “Channel Width Response to Differential Uplift.” *Journal of Geophysical Research: Earth Surface (2003–2012)* 112 (F2).

Attal, M., G. E. Tucker, A. C. Whittaker, P. A. Cowie, and G P Roberts. 2008. “Modeling Fluvial Incision and Transient Landscape Evolution: Influence of Dynamic Channel Adjustment.” *Journal of Geophysical Research* 113: F03013.

Beaumont, C., P. Fullsack, and J. Hamilton. 1992. “Erosional Control of Active Compressional Orogens.” *Thrust Tectonics* 99: 1–18.

Braun, Jean, and Malcom Sambridge. 1997. “Modelling Landscape Evolution on Geological Time Scales: a New Method Based on Irregular Spatial Discretization.” *Basin Research* 9: 27–52.

Carretier, S., and F. Lucazeau. 2005. “How Does Alluvial Sedimentation at Range Fronts Modify the Erosional Dynamics of Mountain Catchments?” *Basin Research* 17 (3): 361–381.

Chase, Clement G. 1992. “Fluvial Landsculpting and the Fractal Dimension of Topography.” *Geomorphology* 5: 39–57.

Clevis, Q., G. Tucker, S. Lancaster, A. Desitter, N. Gasparini, and G. Lock. 2006. “Geoarchaeological Simulation of Meandering River Deposits and Settlement Distributions: a Three-Dimensional Approach.” *Computers and Geosciences* 21 (8): 843–874.

Collins, DBG, RL Bras, and GE Tucker. 2004. “Modeling the Effects of Vegetation-Erosion Coupling on Landscape Evolution.” *Journal of Geophysical Research—Earth Surface* 109 (F3) ({AUG 3}). doi:10.1029/2003JF000028.

Coulthard, TJ, MJ Kirkby, and MG Macklin. 1996. “A Cellular Automaton Landscape Evolution Model.” *Proceedings of the First International Conference on GeoComputation*.

Crave, A, and Ph Davy. 2001. “A Stochastic ’Precipiton’ Model for Simulating Erosion/Sedimentation Dynamics.” *Computers and Geosciences* 27: 815–827.

Crosby, Benjamin T., Kelin X. Whipple, Nicole M. Gasparini, and Cameron W. Wobus. 2007. “Formation of Fluvial Hanging Valleys: Theory and Simulation.” *JOURNAL OF GEOPHYSICAL RESEARCH-EARTH SURFACE* 112 (F3) ({AUG 9}). doi:10.1029/2006JF000566.

Culling, WEH. 1963. “Soil Creep and the Development of Hillside Slopes.” *The Journal of Geology*: 127–161.

Densmore, Alexander L, Michael A Ellis, and Robert S Anderson. 1998. “Landsliding and the Evolution of Normal-Fault-Bounded Mountains.” *Journal of Geophysical Research* 103: 15203–15219.

DiBiase, R.A., and K.X. Whipple. 2011. “The Influence of Erosion Thresholds and Runoff Variability on the Relationships Among Topography, Climate, and Erosion Rate.” *Journal of Geophysical Research* 116 (F4): F04036.

Duvall, A., E. Kirby, and D. Burbank. 2004. “Tectonic and Lithologic Controls on Bedrock Channel Profiles and Processes in Coastal California.” *Journal of Geophysical Research* 109 (F3): F03002.

Fagherazzi, S., A.D. Howard, and P.L. Wiberg. 2002. “An Implicit Finite Difference Method for Drainage Basin Evolution.” *Water Resources Research* 38 (7): 21.

Flores-Cervantes, Homero, Erkan Istanbulluoglu, and Rafael L Bras. 2006. “Development of Gullies on the Landscape: a Model of Headcut Retreat Resulting from Plunge Pool Erosion.” *Journal of Geophysical Research* 111: F01010.

Garcia-Castellanos, D. 2002. “Interplay Between Lithospheric Flexure and River Transport in Foreland Basins.” *Basin Research* 14 (2): 89–104.

Gasparini, N., G. Tucker, and R. Bras. 1999. "Downstream fining through selective particle sorting in an equilibrium drainage network." *Geology* 27(12), 1079.

Gasparini, N., G. Tucker, and R. Bras. 2004. "Network-scale dynamics of grain-size sorting: Implications for downstream fining, stream-profile concavity, and drainage basin morphology. *Earth Surface Processes and Landforms* 29(4), 401-421.

Gasparini, N. M., K. X. Whipple, and R. L. Bras. 2007. “Predictions of Steady State and Transient Landscape Morphology Using Sediment-Flux-Dependent River Incision Models.” *Journal of Geophysical Research Earth Surface* 112 (F3) ({MAY 19}). doi: 10.1029/2006JF000567.

Gilbert, GK. 1877. “Report on the Geology of the Henry Mountains: US Geog. and Geol.” *Survey, Rocky Mtn. Region* 160.

Herman, F., and J. Braun. 2006. “Fluvial Response to Horizontal Shortening and Glaciations: a Study in the Southern Alps of New Zealand.” *Journal of Geophysical Research-Earth Surface* 111 (F1): F01008.

Herman, F., and J. Braun. 2008. “Evolution of the Glacial Landscape of the Southern Alps of New Zealand: Insights from a Glacial Erosion Model.” *Journal of Geophysical Research-Earth Surface* 113 (F2): F02009.

Herman, F., J. Braun, and W.J. Dunlap. 2007. “Tectonomorphic Scenarios in the Southern Alps of New Zealand.” *Journal of Geophysical Research-Solid Earth* 112 (B4): B04201.

Howard, A. 1971. “Simulation Model of Stream Capture.” *Geological Society of America Bulletin* 82 (5): 1355–1376.

Howard, Alan D. 1994. “A Detachment-Limited Model of Drainage Basin Evolution.” *Water Resources Research* 30 (7): 2261–2285.

Huang, X., and J.D. Niemann. 2006. “An Evaluation of the Geomorphically Effective Event for Fluvial Processes over Long Periods.” *Journal of Geophysical Research* 111 (F3): F03015.

Istanbulluoglu, Erkan, and Rafael L. Bras. 2005. “Vegetation-Modulated Landscape Evolution: Effects of Vegetation on Landscape Processes, Drainage Density, and Topography.” *Journal of Geophysical Research* 110: F02012.

Istanbulluoglu, Erkan, Rafael L. Bras, Homero Flores-Cervantes, and Gregory E. Tucker. 2005. “Implications of Bank Failures and Fluvial Erosion for Gully Development: Field Observations and Modeling.” *Journal of Geophysical Research* 110: F01014.

Kirkby, MJ. 1971. “Hillslope Process-Response Models Based on the Continuity Equation.” *Inst. Br. Geogr. Spec. Publ* 3: 15–30.

Lancaster, S.T., S.K. Hayes, and G.E. Grant. 2003. “Effects of Wood on Debris Flow Runout in Small Mountain Watersheds.” *Water Resources Research* 39 (6): 1168.

Lavé, J., and JP Avouac. 2001. “Fluvial Incision and Tectonic Uplift Across the Himalayas of Central Nepal.” *Journal of Geophysical Research* 106 (B11): 26561.

Mitas, L., and H. Mitasova. 1998. “Distributed Soil Erosion Simulation for Effective Erosion Prevention.” *Water Resources Research* 34 (3): 505–516.

Murray, A.B., and C. Paola. 1994. “A Cellular Model of Braided Rivers.” *Nature* 371 (6492): 54–57.

Perron, J.T. 2011. “Numerical Methods for Nonlinear Hillslope Transport Laws.” *Journal of Geophysical Research* 116 (F2): F02021.

Press, W.H., S.A. Teukolsky, W.T. Vetterling, and B.P. Flannery. 2007. *Numerical Recipes: the Art of Scientific Computing*. Cambridge Univ Pr.

Schoorl, JM, A. Veldkamp, and J. Bouma. 2002. “Modeling Water and Soil Redistribution in a Dynamic Landscape Context.” *Soil Science Society of America Journal* 66 (5): 1610.

Snyder, N.P., K.X. Whipple, G.E. Tucker, and D.J. Merritts. 2003. “Importance of a Stochastic Distribution of Floods and Erosion Thresholds in the Bedrock River Incision Problem.” *Journal of Geophysical Research* 108 (B2): 2117.

Sólyom, P.B., and G.E. Tucker. 2004. “Effect of Limited Storm Duration on Landscape Evolution, Drainage Basin Geometry, and Hydrograph Shapes.” *Journal of Geophysical Research* 109: 13.

Tucker, G E, and R L Slingerland. 1994. “Erosional Dynamics, Flexural Isostasy, and Long-Lived Escarpments: a Numerical Modeling Study.” *Journal of Geophysical Research* 99: 12, 229–12, 243.

Tucker, GE. 2004. “Drainage Basin Sensitivity to Tectonic and Climatic Forcing: Implications of a Stochastic Model for the Role of Entrainment and Erosion Thresholds.” *Earth Surface Processes and Landforms* 29 (2) ({FEB}): 185–205. doi:10.1002/esp.1020.

Tucker, Gregory E, and Rafael L Bras. 2000. “A Stochastic Approach to Modeling the Role of Rainfall Variability in Drainage Basin Evolution.” *Water Resources Research* 36 (7): 1953–1964.

Tucker, Gregory E, and Gregory R Hancock. 2010. “Modelling Landscape Evolution.” *Earth Surface Processes and Landforms* 46: 28–50.

Tucker, Gregory E, S T Lancaster, Nicole M Gasparini, and Rafael L Bras. 2001. “The Channel-Hillslope Integrated Landscape Development Model (CHILD).” In *Landscape Erosion and Evolution Modeling*, edited by R. S. Harmon and W. W. Doe, 349–388. Dordrecht: Kluwer Press.

Tucker, Gregory E, S T Lancaster, Nicole M Gasparini, Rafael L Bras, and S M Rybarczyk. 2001. “An Object-Oriented Framework for Hydrologic and Geomorphic Modeling Using Triangular Irregular Networks.” *Computers and Geosciences* 27: 959–973.

Turowski, Jens M, Dimitri Lague, and Niels Hovius. 2009. “Response of Bedrock Channel Width to Tectonic Forcing: Insights from a Numerical Model, Theoretical Considerations, and Comparison with Field Data.” *Journal of Geophysical Research: Earth Surface (2003–2012)* 114 (F3).

Whipple, K X, and G E Tucker. 1999. “Dynamics of the Stream-Power River Incision Model: Implications for Height Limits of Mountain Ranges, Landscape Response Timescales, and Research Needs.” *Journal of Geophysical Research* 104: 17661–17674.

Whipple, K. X., and G. E. Tucker. 2002. “Implications of Sediment-Flux-Dependent River Incision Models for Landscape Evolution.” *Journal of Geophysical Research* 107.

Whittaker, A.C., P.A. Cowie, M. Attal, G.E. Tucker, and G.P. Roberts. 2007. “Bedrock Channel Adjustment to Tectonic Forcing: Implications for Predicting River Incision Rates.” *Geology* 35 (2): 103.

Willgoose, Garry, Rafael L Bras, and Ignacio Rodriguez-Iturbe. 1991. “A Coupled Channel Network Growth and Hillslope Evolution Model, 1, Theory.” *Water Resources Research* 27 (7): 1671–1684.

Wobus, C.W., J.W. Kean, G.E. Tucker, and R.S. Anderson. 2008. “Modeling the Evolution of Channel Shape: Balancing Computational Efficiency with Hydraulic Fidelity.” *Journal of Geophysical Research* 113.

Wobus, C.W., G.E. Tucker, and R.S. Anderson. 2006. “Self-Formed Bedrock Channels.” *Geophys. Res. Lett* 33: 1–6.

Yanites, B.J., and G.E. Tucker. 2010. “Controls and Limits on Bedrock Channel Geometry.” *Journal of Geophysical Research* 115 (F4): F04019.

Yanites, B.J., G.E. Tucker, K.J. Mueller, Y.G. Chen, T. Wilcox, S.Y. Huang, and K.W. Shi. 2010. “Incision and Channel Morphology Across Active Structures Along the Peikang River, Central Taiwan: Implications for the Importance of Channel Width.” *Geological Society of America Bulletin* 122 (7-8): 1192.