Viscous Model

dripping_channel_1_init.png
dripping_channel_1_init.png

Viscous Model

0.00

Script

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.

Movies

Viscosity

Temperature

Density

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!)

Viscosity/Rheology

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

Temperature

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

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
  4.