a depth sequence for freshwater ΔRs?

13 views
Skip to first unread message

Erik Marsh

unread,
Mar 26, 2026, 9:41:50 AM (7 days ago) Mar 26
to OxCal
Hi all,
I'm looking at lake core data and the freshwater ΔR increases with depth. I want to make a depth sequence for these ΔRs.

I have a few anchor points (defined ΔRs at set depths), but I want to let ΔRs float at other depths, just like in a P_sequence. The code below shows the logic, but of course Delta_R doesn't run in this context. It seems like I might need to use a general parameter, perhaps with the operator Pois? I haven't worked much with these variables.

Any ideas welcome, thanks!
Erik

P_Sequence(1,2,U(-2,2))
{
Boundary();
Delta_R("R5",U(400,600){ z=500; };
Delta_R("R4",U(300,500)){ z=400; };
Delta_R("known R3",300,30){ z=300; };
Delta_R("R2",U(100,300)){ z=200; };
Delta_R("R1",U(0,200)){ z=100; };
Boundary();
};
Delta_R("=R4");
R_Date(1000,50);

Christopher Ramsey

unread,
Mar 26, 2026, 11:36:39 AM (7 days ago) Mar 26
to OxCal group
Dear Erik

P_Sequence does work with Number. Note however that P_Sequence like Sequence goes from low to high and Delta_R needs to have a distribution so the following will run:

Plot()
{
P_Sequence(1,2,U(-2,2))
{
Boundary();
Number("R1",U(0,200))
{
z=100;
};
Number("R2",U(100,300))
{
z=200;
};
Number("known R3",N(300,30))
{
z=300;
};
Number("=R4")
{
z=400;
};
Number("R5",U(400,600))
{
z=500;
};
Boundary();
};
Delta_R("R4",U(300,500));
R_Date(1000,50);
};

Whether it makes sense or not in the context is a different matter!

Best wishes

Christopher
> --
> You received this message because you are subscribed to the Google Groups "OxCal" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to oxcal+un...@googlegroups.com.
> To view this discussion visit https://groups.google.com/d/msgid/oxcal/b32283b1-c1b5-45aa-b747-c6f445a2d787n%40googlegroups.com.

Erik Marsh

unread,
Mar 26, 2026, 8:59:08 PM (6 days ago) Mar 26
to OxCal
Thanks Christopher, this works well. In this case, I think we will need a new model of lake levels and carbon inputs (above my pay grade) before this code can be built into the age model. But it's great to see how it can be done with cross-referenced Delta_R, Number and P_Sequence. And this is not the only lake where FRE increases with depth; it may be useful to others. –Erik
Reply all
Reply to author
Forward
0 new messages