Gridplates global sediment circulation model
Source code
here,
brief model description here.
- May 16. Bug fixes of the previous revisions, and tuned the physical erosion rate to match the thickness of ocean sediments (after changing the porosity which screwed it up).
- April 26. Revamped the surface geomorphology and the sediment bookkeeping to align the model with the weathering model of Lechuga-Crespo et al 2020 Global and Planetary Change 192 (2020) 103226. New plots with river concentrations and fluxes of multiple solutes.
- April 11. Minor bug fix over the last version -- added pore water to the density of continental sediments, and increased the orogenic erosion rate a bit.
- April 10. Revamped the sediment bookkeeping to allow variable porosity. The goal was to fill up Laurentia with lots of CaCO3 during the early high sea level period. Isostatic subsidence is very sensitive to porosity via bulk density. In the oceans, trying to fill up the margins required a higher porosity of the sediments. Added a "weathering age" calculation of the land sediment. Found that the elevation of land-sedimented grid points next to orogenic areas requires a fairly high land diffusion coefficient, otherwise the soil at the edge of the valley is higher than the mountain.
- March 15, Worked on ocean sedimentation, to try to create more continental rise at the margins. There is a new Laurentia elevation transect animation, and a slice-by-slice comparison of the model end state with observed elevation and sediment thickness. Concentrated sedimentation near the margin by decreasing the ocean sediment transport coefficient by 2 orders of magnitude.
- March 10, Increased Paleozoic sea level rise to 400 meters, which increased the thickness of the CaCO3 on continents to 1.5 km, which is still not enough but proof of concept. Came up with a scheme to only allow ocean-to-continent grid point transitions while limiting the 'land-spreading' in the last simulation.
- March 5, eliminated land sediment transport of CaCO3, turned off continent to ocean transitions. Now the sea floor looks less lumpy near the coasts but the continents get all raggedy because grid points flip to continent randomly at the edges and never flip back.
- February 20, fixed isostacy again to get mid ocean ridge depths right. Next goal is to get more CaCO3 on continents and keep it there.
- February 16, sea level now changes as instructed.
- February 15, fixed problems in isostacy but sea level change was not being applied.
- February 13, turned off subaereal CaCO3 dissolution to see if that's where it was disappearing to.
- February 12, Fixed a problem with subaereal sediment dissolution going below sea level. Doh!
- February 9, Fixed the clay bug and some other problems. Weathering and erosion fluxes are probably wrong, need to be checked.
- February 7, First results
with altitude tracking the Scotese reconstruction (realigned to the Merdith plate locations), and ice
sheets, parameterized as locations where climate model results make it -10 deg. C or colder. There are
problems with the chemical vs physical weathering in this simulation
- January 24, First results
with altitude tracking the Scotese reconstruction (realigned to the Merdith plate locations), and ice
sheets, parameterized as locations where climate model results make it -10 deg. C or colder. There was a bug, no clay production on land, so everywhere is flooded. n
- December 11, The continental crust now thins after orography to minimize secular trend in crust thickness. Land and ocean sediment loads now approximately correct.
- December 9, CaCO3 burial at a specified rate, trying to increase land sediment thickness by not letting the crust get too thick
- December 6, Ca-conserving run, choked for CaCO3
Additional animations: Runoff from CLIMBER-X using CO2 = 280 ppm with interpolated slices, 1120 ppm. Temperature fields from CLIMBER-X, using pCO2 of 280 ppm and 1120 ppm. Comparison of my simple runoff generator vs. Baum et al here, animation of plateID changes (red means information has to migrate from one plate file to another) here.