Skip to content

Commit

Permalink
Last updates, weighting latitude
Browse files Browse the repository at this point in the history
  • Loading branch information
@tth030
tth030 committed Jan 10, 2022
1 parent 350f1d8 commit 9c0edc1
Show file tree
Hide file tree
Showing 9 changed files with 299 additions and 1,257 deletions.
1 change: 1 addition & 0 deletions .gitignore
View file
Original file line number Diff line number Diff line change
Expand Up @@ -133,3 +133,4 @@ data_figshare
.DS_Store
environment.yml_failed
data_figshare.zip
*.pdf
72 changes: 14 additions & 58 deletions isostasy_calib.ipynb
View file
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,7 @@
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"execution_count": null,
"id": "185d8870-6c3e-4b04-80fb-0accd59bd115",
"metadata": {},
"outputs": [],
Expand All @@ -18,47 +18,16 @@
},
{
"cell_type": "code",
"execution_count": 2,
"execution_count": null,
"id": "880ec5ee-c355-4da5-95e2-7de2475ae72d",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"# Column 1: Reference MOR ###################################################################################\n",
"Density and temperature profiles from 2-D geodynamic model M2 (temperature dependent only densities)\n",
"# Column 2: Reference Continental lithosphere ###############################################################\n",
" z bottom (km) = 600000.0\n",
" k asth. = 87.56\n",
" Grad slm = 0.4\n",
" Qbasal = 0.0195488281\n",
" Tbasal = 1520.0\n",
" Tp = 1280.0\n",
"Layer 6 T top = 1330.0 Q top = 19.55 (matSLM)\n",
"Layer 3 T top = 548.0 Q top = 19.55 (matLM1)\n",
"Layer 2 T top = 441.8 Q top = 28.61 (matLC)\n",
"Layer 1 T top = 295.3 Q top = 37.68 (matMC)\n",
"Layer 0 T top = 0.0 Q top = 51.28 (matUC)\n",
" Qsurface = 51.28 mW/m2\n",
"========= SUBSIDENCE AT THE RIDGE WITH FIXED WATER DEPTH ===========\n",
"Reference crustal density = 2860.0 kg/m3\n",
"Average crustal density = 2833.1 kg/m3\n",
"Average degree depletion CLM = -20.0 kg/m3\n",
"WATER DEPTH = 2750.0 m\n",
"SEA LEVEL = -395.7 m\n",
"dh = -3145.7 m\n"
]
}
],
"outputs": [],
"source": [
"# Calibration of the crustal density and degree of depletion of the continentlal lithospheric mantle\n",
"# and sensitivity on water depth/sea level\n",
"#\n",
"\n",
"water_depth_at_ridge = 2750.0e0 # meters = reference water depth at the ridge based on average MOR elevation\n",
"sea_level_ref = -400e0 # meters = reference sea level based on average elevation of continents\n",
"dz = 200e0 # meters\n",
"ref_depth = 125e3 # meters = compensation depth (base of the continental lithosphere)\n",
"\n",
Expand All @@ -74,7 +43,7 @@
"perplex_name = 'slm_sp2008' # composition of the continental lithospheric mantle\n",
"\n",
"\n",
"use_fixed_waterdepth = True # True means that sea level can change but that waterdepth at the ridge is fixed\n",
"use_fixed_waterdepth = False # True means that sea level can change but that waterdepth at the ridge is fixed\n",
"\n",
"\n",
"######## COMMENT/UNCOMMENT OR DEFINE A CASE #########################################################\n",
Expand Down Expand Up @@ -102,8 +71,8 @@
"printOut = True\n",
"use_perplex_table_mantle = False\n",
"allrefdensity_crust = [2860]\n",
"alllith_depletion = [-20]\n",
"allsea_level_ref = [-400]\n",
"alllith_depletion = [-13]\n",
"allsea_level_ref = [-200]\n",
"\n",
"############################################################################################\n",
"\n",
Expand Down Expand Up @@ -162,10 +131,12 @@
" for lith_depletion in alllith_depletion: \n",
" \n",
" # Update the book ------------------------\n",
" book_complith=update_book(book_complith,'refdensity_crust',refdensity_crust,['UC','MC','LC'])\n",
" book_complith=update_book(book_complith,'lith_depletion',lith_depletion,['LM1','LM2','LM3'])\n",
" book_complith=update_book(book_complith,'rho',refdensity_crust,['UC','MC','LC'])\n",
" book_complith=update_book(book_complith,'deltarho',lith_depletion,['LM1','LM2','LM3'])\n",
" # ----------------------------------------\n",
"\n",
" print(book_complith)\n",
" \n",
" if printOut:\n",
" print(\"# Column 2: Reference Continental lithosphere ###############################################################\")\n",
" lith2 = geodyn1d.lithosphere(complith,printOut=False)\n",
Expand Down Expand Up @@ -216,25 +187,10 @@
},
{
"cell_type": "code",
"execution_count": 4,
"execution_count": null,
"id": "82a66fcd-ac1c-42b7-bdbf-c27746ddde16",
"metadata": {},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "c921153b724b43deb171129ab4013e2e",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"outputs": [],
"source": [
"# Calibration of the crustal density and degree of depletion of the continentlal lithospheric mantle\n",
"# and sensitivity on water depth/sea level\n",
Expand All @@ -252,8 +208,8 @@
"\n",
"# Calibration of the crustal density and degree of depletion of the continentlal lithospheric mantle\n",
"#\n",
"#refdensity_mantle = \"SP2008_speed3_0.55_bl400_waterload\" # fixed sea level\n",
"refdensity_mantle = \"SP2008_speed3_0.55_waterdepth2750_waterload\" # fixed water depth the ridge\n",
"refdensity_mantle = \"SP2008_speed3_0.55_bl400_waterload\" # fixed sea level\n",
"#refdensity_mantle = \"SP2008_speed3_0.55_waterdepth2750_waterload\" # fixed water depth the ridge\n",
"\n",
"\n",
"# Sensitivity on water depth/sea level for temperature dependent only densities\n",
Expand Down
133 changes: 10 additions & 123 deletions isostasy_calib2.ipynb
View file
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,7 @@
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"execution_count": null,
"id": "166d6aa2-7d95-4802-80e6-3ec96e5bca53",
"metadata": {},
"outputs": [],
Expand All @@ -13,7 +13,7 @@
},
{
"cell_type": "code",
"execution_count": 2,
"execution_count": null,
"id": "bd867c0d-6056-422b-840e-8864fcb81b7d",
"metadata": {},
"outputs": [],
Expand All @@ -28,39 +28,10 @@
},
{
"cell_type": "code",
"execution_count": 3,
"execution_count": null,
"id": "828a9b1c-3166-4221-b179-d3e20f8e201e",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" z bottom (km) = 600000.0\n",
" k asth. = 87.56\n",
" Grad slm = 0.4\n",
" Qbasal = 0.0195488281\n",
" Tbasal = 1520.0\n",
" Tp = 1280.0\n",
"Layer 6 T top = 1330.0 Q top = 19.55 (matSLM)\n",
"Layer 3 T top = 548.0 Q top = 19.55 (matLM1)\n",
"Layer 2 T top = 441.8 Q top = 28.61 (matLC)\n",
"Layer 1 T top = 295.3 Q top = 37.68 (matMC)\n",
"Layer 0 T top = 0.0 Q top = 51.28 (matUC)\n",
" Qsurface = 51.28 mW/m2\n",
"Pref = 3.79540e+09 Pa with compensation depth of 125.0 km\n",
"========= SUBSIDENCE AT THE REFERENCE MOR ===========\n",
"WATER DEPTH = 2743.614290913309 m\n",
"SEA LEVEL = -400.0 m\n",
"Reference dh = -3143.6 m\n",
"========= SUBSIDENCE AT THE RIDGE MOR2 ===========\n",
"correction mantle = -11.2 kg/m3\n",
"WATER DEPTH = 2742.908277926956 m\n",
"SEA LEVEL = -400.0 m\n",
"diff with ref = 0.7 m\n"
]
}
],
"outputs": [],
"source": [
"# Calibration of the reference mantle density with temperature dependent only densities (model M2)\n",
"#\n",
Expand Down Expand Up @@ -161,53 +132,10 @@
},
{
"cell_type": "code",
"execution_count": 4,
"execution_count": null,
"id": "84611c37-f8d9-4b48-9e02-053f000d41a7",
"metadata": {},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "094ca40ee5f4478dae8605419bf8810a",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "9d87d3c878bc40d196648ddf1cfa98fa",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "ac64aabcedc14815812af3e1be2be7c5",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"outputs": [],
"source": [
"# Display temperature, pressure, and density profiles of the two columns:\n",
"# MOR: reference MOR column from model M2\n",
Expand Down Expand Up @@ -261,28 +189,10 @@
},
{
"cell_type": "code",
"execution_count": 5,
"execution_count": null,
"id": "97a4ebeb-d6af-410a-90dc-c584093fec14",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Reference density mantle: 3250.0 kg/m3\n",
"Reference density mantle: 3260.0 kg/m3\n",
"Reference density mantle: 3270.0 kg/m3\n",
"Reference density mantle: 3280.0 kg/m3\n",
"Reference density mantle: 3290.0 kg/m3\n",
"Reference density mantle: 3300.0 kg/m3\n",
"Reference density mantle: 3310.0 kg/m3\n",
"Reference density mantle: 3320.0 kg/m3\n",
"Reference density mantle: 3330.0 kg/m3\n",
"Reference density mantle: 3340.0 kg/m3\n",
"Done !!\n"
]
}
],
"outputs": [],
"source": [
"# Calibration of the reference mantle density with temperature dependent only densities\n",
"#\n",
Expand Down Expand Up @@ -391,33 +301,10 @@
},
{
"cell_type": "code",
"execution_count": 6,
"execution_count": null,
"id": "49d10926-50a9-4936-907b-722baef58af1",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Calibration:\n",
"The reference mantle density is 3311.6524564678734 kg/m3 to fit the water depth at the MOR (with fixed sea level at -400 m)\n"
]
},
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "9f382472d1d448a9bcccd7a24d510d51",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"outputs": [],
"source": [
"# Calibration of the reference mantle density with temperature dependent only densities\n",
"# Plot the solution from STEP 2 (Figure 10d)\n",
Expand Down
Loading

0 comments on commit 9c0edc1

Please sign in to comment.