Diffusive-advective evolution of porewater chemistry in a marine sediment column with changing seafloor chemistry linked to climate and ice sheet grounding line behavior.
Based on diffusion-advection model framework of Neuhaus & Tulaczyk 2023, translated into Julia from the original MATLAB for speed 🚀
The model is wrapped in a Markov chain Monte Carlo algorithm. Given a climate record, the MCMC inverts porewater chemistry data for posterior chains of the onset of glaciation at the marine site, the climate sensitivity of the system, and the basal boundary conditions of the sediment column.
To install CorePore.jl on your own computer, just type ] into the Julia REPL to enter the built-in package manager and then type:
add https://github.com/grahamedwards/CorePore.jland hit enter.
After installing, just type using CorePore to use.
The primary functionality of this package is hosted in the function porewatermetropolis. Please read the documentation for its full functionality and the underlying forward models. The example below (also available as notebook example.ipynb) will give you a quick start:
using CorePore
p = Proposal(5320., 1e-4,1e-3,3.5, 4.2, 1000, deepbonney()...)
jumpsize = Proposal(20., .1e-4, .1e-3, .1, .1,10, 10, 1.)
climate, k = LR04(), Constants()
priors = ProposalPriors(climate, k)
chains = porewatermetropolis(p, jumpsize, priors, andrill2a(), climate, k; burnin=1000, chainsteps=1000, k=Constants(), seawater=mcmurdosound(), onlychloride=true)