A Matlab function “isostaticmassflux1D” is included as part of the supplemental materials to the manuscript: Conder, J.A., Oceanic isostasy as a trigger for the rift to drift transition, submitted to Geology, 2021.
The function numerically calculates the asthenosphere mass flux needed to maintain isostasy with differential thermal subsidence. Basic theory is laid out in Conder, Non-Pratt component of oceanic isostasy, Lithosphere, 2012. Written in Matlab ver. R2021a (Mac).
>> U = isostaticmassflux1D(x,t)
>> U = isostaticmassflux1D(x,t,water)
The function requires two inputs: x (spatial vector in km across axis) and t (thermal ages in Myr). A third optional input, water, is a binary vector denoting which parts of the line are above water (0) and which are below (1). Default without the input is all submerged.
Output mass flux, U, has units km3/Myr/km
August 10, 2021