Viscous Model


Viscous Model



This is a link to the ipython notebook *.ipynb

Geometry and model parm summary

Set up a vertical channel and made the lower crust very weak (low visc) and hot (asth temp) so that it would be sure to detach.





AGU abstract

dynamic topography



Add To Cart

Vertical & horizontal detachment

Vertical channel is set to vc=1.95 (so that is .05 wide or 132 km wide!)


No plastic rheology. Only viscosity variation with depth.

# Ideal continental crust: Jelly sandwich
AsthViscosity        =  0.01
lowerMantleViscosity =  100.0
UCrustViscosity      =  100.0
LCrustViscosity      =  0.01
LithViscosity        =  10.0
CUCViscosity         =  100.0
CLCViscosity         =  100.0
CLithViscosity       =  100.0


Temp of lower crust set to 1. Top = 0 base =1. 

I used a conductive lithosphere as a temperature initial condition defined by a variable heatflux between the craton (45mW/m2) and cordillera (60 mW/m2).  The Cordillera Lithosphere is set to be 1300C at 120km, while the craton is set to be 1300C at 200km.


Density of the lith is greater than the asthenosphere. It wants to sink. Temperature dependent density is applied like in Currie et al.

To do:

  1. Implement Plastic rheology
  2. Implement strain weakening like shear stress example
  3. Figure out how to extract topography from models (see below)
  4. Investigate if need air layer in models
  5. Figure out min model parms to product delam and dripping like Adams stuff.
  6. Write up equations using mac (for the formulas)

Very helpful papers

  1. Hardebol et al 2012
  2. Petersen et al 2016
  3. Sharples et al 2015