This code is used to generate age-depth models and pollen diagrams for a forthcoming publication:
David K. Wright, Sarah J. Ivory, Jago J. Birk, Jeong-Heon Choi, Benjamin Davies, Sabine Fiedler, Jacob Davis, Potiphar Kaliba, Jessica C. Thompson.”Palaeoenvironmental data track Late Quaternary anthropogenic impacts on vegetation and landscapes in Mzimba, northern Malawi” Submitted to Frontiers in Environmental Archaeology, May 2023.
This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.
rm(list = ls())
library(rbacon)
Loading required package: IntCal
library(rioja)
This is rioja 1.0-5
Bacon(core="AU3",d.min=0,d.max=67,accept.suggestions = TRUE,verbose=FALSE)
The run's files will be put in this folder: C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU3
Reading C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU3/AU3_15.bacon
Constant calibration curve.
IntCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_intcal20.14C
Marine20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_marine20.14C
SHCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_shcal20.14C
Added det: NA: 49.0+-964000.0 d=0 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69826: 49.0+-29.0 d=14 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69825: 267.0+-29.0 d=19 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69824: 964.0+-29.0 d=28 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: NA: 964.0+-964000000.0 d=67 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
BaconFixed: Bacon jumps model with fixed c's.
K=15, H=0, dim=17, Seed=161655, Dc=5.000000, c(0)=0.000000, c(K)=75.000000
twalk: 3740000 iterations to run, Sat Jun 3 07:54:09 2023
twalk thinning: 1 out of every 85 accepted iterations will be saved in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU3/AU3_15.out
twalk: Finished, 2.2% of moved pars per iteration (ratio 83855.705883/3740000). Output in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU3/AU3_15.out,
Sat Jun 3 07:54:13 2023
bacon: 184049 WarnBeyondLimits warnings:
bacon: WARNING: calibration attempted beyond MinYr= -96.000000 or MaxYr= 1000000.000000
bacon: burn in (initial iterations which will be removed): 17000
Eso es to...eso es to...eso es to...eso es toooodo amigos!
|
| | 0%
|
|= | 1%
|
|== | 3%
|
|=== | 4%
|
|==== | 6%
|
|===== | 7%
|
|====== | 9%
|
|======= | 10%
|
|======== | 12%
|
|========= | 13%
|
|========== | 15%
|
|=========== | 16%
|
|============ | 18%
|
|============= | 19%
|
|============== | 21%
|
|=============== | 22%
|
|================ | 24%
|
|================= | 25%
|
|================== | 27%
|
|=================== | 28%
|
|==================== | 30%
|
|===================== | 31%
|
|====================== | 33%
|
|======================= | 34%
|
|======================== | 36%
|
|========================= | 37%
|
|========================== | 39%
|
|=========================== | 40%
|
|============================ | 42%
|
|============================= | 43%
|
|============================== | 45%
|
|=============================== | 46%
|
|================================ | 48%
|
|================================= | 49%
|
|================================== | 51%
|
|=================================== | 52%
|
|==================================== | 54%
|
|===================================== | 55%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|============================================ | 66%
|
|============================================= | 67%
|
|============================================== | 69%
|
|=============================================== | 70%
|
|================================================ | 72%
|
|================================================= | 73%
|
|================================================== | 75%
|
|=================================================== | 76%
|
|==================================================== | 78%
|
|===================================================== | 79%
|
|====================================================== | 81%
|
|======================================================= | 82%
|
|======================================================== | 84%
|
|========================================================= | 85%
|
|========================================================== | 87%
|
|=========================================================== | 88%
|
|============================================================ | 90%
|
|============================================================= | 91%
|
|============================================================== | 93%
|
|=============================================================== | 94%
|
|================================================================ | 96%
|
|================================================================= | 97%
|
|================================================================== | 99%
|
|===================================================================| 100%
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|==== | 7%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|====== | 10%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|======== | 13%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|========== | 16%
|
|=========== | 16%
|
|=========== | 17%
|
|============ | 17%
|
|============ | 18%
|
|============ | 19%
|
|============= | 19%
|
|============= | 20%
|
|============== | 20%
|
|============== | 21%
|
|============== | 22%
|
|=============== | 22%
|
|=============== | 23%
|
|================ | 23%
|
|================ | 24%
|
|================ | 25%
|
|================= | 25%
|
|================= | 26%
|
|================== | 26%
|
|================== | 27%
|
|================== | 28%
|
|=================== | 28%
|
|=================== | 29%
|
|==================== | 29%
|
|==================== | 30%
|
|==================== | 31%
|
|===================== | 31%
|
|===================== | 32%
|
|====================== | 32%
|
|====================== | 33%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 35%
|
|======================== | 36%
|
|========================= | 37%
|
|========================= | 38%
|
|========================== | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|=========================== | 41%
|
|============================ | 41%
|
|============================ | 42%
|
|============================= | 43%
|
|============================= | 44%
|
|============================== | 44%
|
|============================== | 45%
|
|=============================== | 46%
|
|=============================== | 47%
|
|================================ | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================= | 50%
|
|================================== | 50%
|
|================================== | 51%
|
|=================================== | 52%
|
|=================================== | 53%
|
|==================================== | 53%
|
|==================================== | 54%
|
|===================================== | 55%
|
|===================================== | 56%
|
|====================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================= | 59%
|
|======================================== | 59%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================= | 62%
|
|========================================== | 62%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|=========================================== | 65%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 67%
|
|============================================= | 68%
|
|============================================== | 68%
|
|============================================== | 69%
|
|=============================================== | 69%
|
|=============================================== | 70%
|
|=============================================== | 71%
|
|================================================ | 71%
|
|================================================ | 72%
|
|================================================= | 72%
|
|================================================= | 73%
|
|================================================= | 74%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 75%
|
|=================================================== | 76%
|
|=================================================== | 77%
|
|==================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 78%
|
|===================================================== | 79%
|
|===================================================== | 80%
|
|====================================================== | 80%
|
|====================================================== | 81%
|
|======================================================= | 81%
|
|======================================================= | 82%
|
|======================================================= | 83%
|
|======================================================== | 83%
|
|======================================================== | 84%
|
|========================================================= | 84%
|
|========================================================= | 85%
|
|========================================================= | 86%
|
|========================================================== | 86%
|
|========================================================== | 87%
|
|=========================================================== | 87%
|
|=========================================================== | 88%
|
|=========================================================== | 89%
|
|============================================================ | 89%
|
|============================================================ | 90%
|
|============================================================= | 90%
|
|============================================================= | 91%
|
|============================================================= | 92%
|
|============================================================== | 92%
|
|============================================================== | 93%
|
|=============================================================== | 93%
|
|=============================================================== | 94%
|
|=============================================================== | 95%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 96%
|
|================================================================= | 97%
|
|================================================================= | 98%
|
|================================================================== | 98%
|
|================================================================== | 99%
|
|===================================================================| 99%
|
|===================================================================| 100%png
2
core<-rep("AU3",9)
pollenNo<-c("POL11","POL10","POL9","POL8","POL12","POL13","POL16","POL14","POL15")
depth<-c(3,14,19,28,34,40,46,55,67)
ageBP<-c()
ageMin<-c()
ageMax<-c()
for (i in depth){
ageBP<-append(ageBP,round(mean(Bacon.Age.d(i))))
ageMin<-append(ageMin,round(min(Bacon.Age.d(i))))
ageMax<-append(ageMax,round(max(Bacon.Age.d(i))))
}
AUG3data<-data.frame(core,pollenNo,depth,ageBP,ageMin,ageMax)
print(AUG3data)
write.csv(AUG3data,"AUG3data.csv",row.names = FALSE)
Bacon(core="AU7",d.min=0,d.max=77,accept.suggestions = TRUE,verbose=FALSE)
The run's files will be put in this folder: C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU7
Reading C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU7/AU7_17.bacon
Constant calibration curve.
IntCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_intcal20.14C
Marine20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_marine20.14C
SHCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_shcal20.14C
Added det: NA: 149.0+-2388000.0 d=0 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69829: 149.0+-29.0 d=15 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69526: 792.0+-27.0 d=32 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69527: 1233.0+-28.0 d=51 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69528: 927.0+-28.0 d=54 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69529: 2388.0+-37.0 d=77 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
BaconFixed: Bacon jumps model with fixed c's.
K=17, H=0, dim=19, Seed=814200, Dc=5.000000, c(0)=0.000000, c(K)=85.000000
twalk: 4180000 iterations to run, Sat Jun 3 07:54:25 2023
twalk thinning: 1 out of every 95 accepted iterations will be saved in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU7/AU7_17.out
twalk: Finished, 1.2% of moved pars per iteration (ratio 49705.526316/4180000). Output in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU7/AU7_17.out,
Sat Jun 3 07:54:31 2023
bacon: 185239 WarnBeyondLimits warnings:
bacon: WARNING: calibration attempted beyond MinYr= -73.000000 or MaxYr= 1000000.000000
bacon: burn in (initial iterations which will be removed): 19000
Eso es to...eso es to...eso es to...eso es toooodo amigos!
|
| | 0%
|
|= | 1%
|
|== | 3%
|
|=== | 4%
|
|=== | 5%
|
|==== | 6%
|
|===== | 8%
|
|====== | 9%
|
|======= | 10%
|
|======== | 12%
|
|========= | 13%
|
|========== | 14%
|
|========== | 16%
|
|=========== | 17%
|
|============ | 18%
|
|============= | 19%
|
|============== | 21%
|
|=============== | 22%
|
|================ | 23%
|
|================= | 25%
|
|================= | 26%
|
|================== | 27%
|
|=================== | 29%
|
|==================== | 30%
|
|===================== | 31%
|
|====================== | 32%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 36%
|
|========================= | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|============================ | 42%
|
|============================= | 43%
|
|============================== | 44%
|
|============================== | 45%
|
|=============================== | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================== | 51%
|
|=================================== | 52%
|
|==================================== | 53%
|
|===================================== | 55%
|
|===================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================== | 62%
|
|=========================================== | 64%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 68%
|
|============================================== | 69%
|
|=============================================== | 70%
|
|================================================ | 71%
|
|================================================= | 73%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 79%
|
|====================================================== | 81%
|
|======================================================= | 82%
|
|======================================================== | 83%
|
|========================================================= | 84%
|
|========================================================= | 86%
|
|========================================================== | 87%
|
|=========================================================== | 88%
|
|============================================================ | 90%
|
|============================================================= | 91%
|
|============================================================== | 92%
|
|=============================================================== | 94%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 97%
|
|================================================================== | 99%
|
|===================================================================| 100%
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|==== | 7%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|====== | 10%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|======== | 13%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|========== | 16%
|
|=========== | 16%
|
|=========== | 17%
|
|============ | 17%
|
|============ | 18%
|
|============ | 19%
|
|============= | 19%
|
|============= | 20%
|
|============== | 20%
|
|============== | 21%
|
|============== | 22%
|
|=============== | 22%
|
|=============== | 23%
|
|================ | 23%
|
|================ | 24%
|
|================ | 25%
|
|================= | 25%
|
|================= | 26%
|
|================== | 26%
|
|================== | 27%
|
|================== | 28%
|
|=================== | 28%
|
|=================== | 29%
|
|==================== | 29%
|
|==================== | 30%
|
|==================== | 31%
|
|===================== | 31%
|
|===================== | 32%
|
|====================== | 32%
|
|====================== | 33%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 35%
|
|======================== | 36%
|
|========================= | 37%
|
|========================= | 38%
|
|========================== | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|=========================== | 41%
|
|============================ | 41%
|
|============================ | 42%
|
|============================= | 43%
|
|============================= | 44%
|
|============================== | 44%
|
|============================== | 45%
|
|=============================== | 46%
|
|=============================== | 47%
|
|================================ | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================= | 50%
|
|================================== | 50%
|
|================================== | 51%
|
|=================================== | 52%
|
|=================================== | 53%
|
|==================================== | 53%
|
|==================================== | 54%
|
|===================================== | 55%
|
|===================================== | 56%
|
|====================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================= | 59%
|
|======================================== | 59%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================= | 62%
|
|========================================== | 62%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|=========================================== | 65%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 67%
|
|============================================= | 68%
|
|============================================== | 68%
|
|============================================== | 69%
|
|=============================================== | 69%
|
|=============================================== | 70%
|
|=============================================== | 71%
|
|================================================ | 71%
|
|================================================ | 72%
|
|================================================= | 72%
|
|================================================= | 73%
|
|================================================= | 74%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 75%
|
|=================================================== | 76%
|
|=================================================== | 77%
|
|==================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 78%
|
|===================================================== | 79%
|
|===================================================== | 80%
|
|====================================================== | 80%
|
|====================================================== | 81%
|
|======================================================= | 81%
|
|======================================================= | 82%
|
|======================================================= | 83%
|
|======================================================== | 83%
|
|======================================================== | 84%
|
|========================================================= | 84%
|
|========================================================= | 85%
|
|========================================================= | 86%
|
|========================================================== | 86%
|
|========================================================== | 87%
|
|=========================================================== | 87%
|
|=========================================================== | 88%
|
|=========================================================== | 89%
|
|============================================================ | 89%
|
|============================================================ | 90%
|
|============================================================= | 90%
|
|============================================================= | 91%
|
|============================================================= | 92%
|
|============================================================== | 92%
|
|============================================================== | 93%
|
|=============================================================== | 93%
|
|=============================================================== | 94%
|
|=============================================================== | 95%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 96%
|
|================================================================= | 97%
|
|================================================================= | 98%
|
|================================================================== | 98%
|
|================================================================== | 99%
|
|===================================================================| 99%
|
|===================================================================| 100%png
2
core<-rep("AU7",9)
pollenNo<-c("POL36","POL37","POL38","POL39","POL40","POL41","POL42","POL43","POL44")
depth<-c(6,15,26,32,42,51,54,66,77)
ageBP<-c()
ageMin<-c()
ageMax<-c()
for (i in depth){
ageBP<-append(ageBP,round(mean(Bacon.Age.d(i))))
ageMin<-append(ageMin,round(min(Bacon.Age.d(i))))
ageMax<-append(ageMax,round(max(Bacon.Age.d(i))))
}
AUG7data<-data.frame(core,pollenNo,depth,ageBP,ageMin,ageMax)
print(AUG7data)
write.csv(AUG7data,"AUG7data.csv",row.names = FALSE)
#Run age-depth model between 0 and max depth
Bacon(core="AU8",d.min=0,d.max=146,depth.unit = "cm",accept.suggestions = TRUE,verbose=FALSE)
The run's files will be put in this folder: C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU8
Reading C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU8/AU8_31.bacon
Constant calibration curve.
IntCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_intcal20.14C
Marine20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_marine20.14C
SHCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_shcal20.14C
Added det: NA: 1794.0+-15742000.0 d=0 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69524: 1794.0+-32.0 d=11 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69523: 7759.0+-124.0 d=59 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69828: 9690.0+-51.0 d=74 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69522: 15742.0+-59.0 d=86 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: NA: 15742.0+-15742000000.0 d=146 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
BaconFixed: Bacon jumps model with fixed c's.
K=31, H=0, dim=33, Seed=970277, Dc=5.000000, c(0)=0.000000, c(K)=155.000000
twalk: 7260000 iterations to run, Sat Jun 3 07:54:59 2023
twalk thinning: 1 out of every 165 accepted iterations will be saved in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU8/AU8_31.out
twalk: Finished, 0.5% of moved pars per iteration (ratio 38029.303031/7260000). Output in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU8/AU8_31.out,
Sat Jun 3 07:55:13 2023
bacon: 173536 WarnBeyondLimits warnings:
bacon: WARNING: calibration attempted beyond MinYr= -73.000000 or MaxYr= 1000000.000000
bacon: burn in (initial iterations which will be removed): 33000
Eso es to...eso es to...eso es to...eso es toooodo amigos!
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 3%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|====== | 10%
|
|======= | 10%
|
|======= | 11%
|
|======== | 12%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|=========== | 16%
|
|=========== | 17%
|
|============ | 18%
|
|============= | 19%
|
|============= | 20%
|
|============== | 21%
|
|=============== | 22%
|
|=============== | 23%
|
|================ | 23%
|
|================ | 24%
|
|================= | 25%
|
|================= | 26%
|
|================== | 27%
|
|=================== | 28%
|
|=================== | 29%
|
|==================== | 29%
|
|==================== | 30%
|
|===================== | 31%
|
|===================== | 32%
|
|====================== | 32%
|
|====================== | 33%
|
|====================== | 34%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 36%
|
|========================= | 37%
|
|========================= | 38%
|
|========================== | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|============================ | 41%
|
|============================ | 42%
|
|============================= | 43%
|
|============================= | 44%
|
|============================== | 45%
|
|=============================== | 46%
|
|=============================== | 47%
|
|================================ | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================== | 50%
|
|================================== | 51%
|
|=================================== | 52%
|
|=================================== | 53%
|
|==================================== | 53%
|
|==================================== | 54%
|
|===================================== | 55%
|
|====================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================= | 59%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================= | 62%
|
|========================================== | 62%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 66%
|
|============================================= | 67%
|
|============================================= | 68%
|
|============================================== | 68%
|
|============================================== | 69%
|
|=============================================== | 70%
|
|=============================================== | 71%
|
|================================================ | 71%
|
|================================================ | 72%
|
|================================================= | 73%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 76%
|
|=================================================== | 77%
|
|==================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 79%
|
|====================================================== | 80%
|
|====================================================== | 81%
|
|======================================================= | 82%
|
|======================================================== | 83%
|
|======================================================== | 84%
|
|========================================================= | 85%
|
|========================================================= | 86%
|
|========================================================== | 86%
|
|========================================================== | 87%
|
|=========================================================== | 88%
|
|============================================================ | 89%
|
|============================================================ | 90%
|
|============================================================= | 90%
|
|============================================================= | 91%
|
|============================================================= | 92%
|
|============================================================== | 92%
|
|============================================================== | 93%
|
|=============================================================== | 94%
|
|=============================================================== | 95%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 97%
|
|================================================================== | 98%
|
|================================================================== | 99%
|
|===================================================================| 99%
|
|===================================================================| 100%
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|==== | 7%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|====== | 10%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|======== | 13%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|========== | 16%
|
|=========== | 16%
|
|=========== | 17%
|
|============ | 17%
|
|============ | 18%
|
|============ | 19%
|
|============= | 19%
|
|============= | 20%
|
|============== | 20%
|
|============== | 21%
|
|============== | 22%
|
|=============== | 22%
|
|=============== | 23%
|
|================ | 23%
|
|================ | 24%
|
|================ | 25%
|
|================= | 25%
|
|================= | 26%
|
|================== | 26%
|
|================== | 27%
|
|================== | 28%
|
|=================== | 28%
|
|=================== | 29%
|
|==================== | 29%
|
|==================== | 30%
|
|==================== | 31%
|
|===================== | 31%
|
|===================== | 32%
|
|====================== | 32%
|
|====================== | 33%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 35%
|
|======================== | 36%
|
|========================= | 37%
|
|========================= | 38%
|
|========================== | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|=========================== | 41%
|
|============================ | 41%
|
|============================ | 42%
|
|============================= | 43%
|
|============================= | 44%
|
|============================== | 44%
|
|============================== | 45%
|
|=============================== | 46%
|
|=============================== | 47%
|
|================================ | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================= | 50%
|
|================================== | 50%
|
|================================== | 51%
|
|=================================== | 52%
|
|=================================== | 53%
|
|==================================== | 53%
|
|==================================== | 54%
|
|===================================== | 55%
|
|===================================== | 56%
|
|====================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================= | 59%
|
|======================================== | 59%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================= | 62%
|
|========================================== | 62%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|=========================================== | 65%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 67%
|
|============================================= | 68%
|
|============================================== | 68%
|
|============================================== | 69%
|
|=============================================== | 69%
|
|=============================================== | 70%
|
|=============================================== | 71%
|
|================================================ | 71%
|
|================================================ | 72%
|
|================================================= | 72%
|
|================================================= | 73%
|
|================================================= | 74%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 75%
|
|=================================================== | 76%
|
|=================================================== | 77%
|
|==================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 78%
|
|===================================================== | 79%
|
|===================================================== | 80%
|
|====================================================== | 80%
|
|====================================================== | 81%
|
|======================================================= | 81%
|
|======================================================= | 82%
|
|======================================================= | 83%
|
|======================================================== | 83%
|
|======================================================== | 84%
|
|========================================================= | 84%
|
|========================================================= | 85%
|
|========================================================= | 86%
|
|========================================================== | 86%
|
|========================================================== | 87%
|
|=========================================================== | 87%
|
|=========================================================== | 88%
|
|=========================================================== | 89%
|
|============================================================ | 89%
|
|============================================================ | 90%
|
|============================================================= | 90%
|
|============================================================= | 91%
|
|============================================================= | 92%
|
|============================================================== | 92%
|
|============================================================== | 93%
|
|=============================================================== | 93%
|
|=============================================================== | 94%
|
|=============================================================== | 95%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 96%
|
|================================================================= | 97%
|
|================================================================= | 98%
|
|================================================================== | 98%
|
|================================================================== | 99%
|
|===================================================================| 99%
|
|===================================================================| 100%png
2
#Create age-depth table
core<-rep("AU8",7)
pollenNo<-c("POL35","POL34","POL33","POL32","POL31","POL30","POL29")
depth<-c(11,45,59,74,86,110,146)
ageBP<-c()
ageMin<-c()
ageMax<-c()
for (i in depth){
ageBP<-append(ageBP,round(mean(Bacon.Age.d(i))))
ageMin<-append(ageMin,round(min(Bacon.Age.d(i))))
ageMax<-append(ageMax,round(max(Bacon.Age.d(i))))
}
AUG8data<-data.frame(core,pollenNo,depth,ageBP,ageMin,ageMax)
print(AUG8data)
write.csv(AUG8data,"AUG8data.csv",row.names = FALSE)
Notes: 1) Terrestrial pollen data are expressed as percentages of all terrestrial species (excl. aquatics). Aquatics (Cyperaceae and Ascolepis) are plotted as percentage of all species (terrestrial + aquatics). 2) In order to view microcharcoal concentrations in the same plot as percentages, these values were transformed using log10, and then multiplied by 10 to exaggerate (e.g., conc of 10000 -> transformed log10(10000) = 4 -> exaggerated 4 * 10 = 40). The x-axis tick mark labels were then manually replaced with the appropriate scientific notation (e.g., 40 = 10000 = 10^4). 3) Age-depth modeling was not used for composite terrace sequence; ages are estimated using adjacent OSL samples (included as “TERRACEdata.csv”). 4) The convention used for the sequences is “ka”, and several of the dates are less than 1000 years and are expressed as decimals (e.g., 0.4 ka for 400 years). The rioja package does not seem to plot ages as non-integer values, so these were plotted here in terms of integer years (e.g., 400 for 400 years), but y-axis tick mark labels were edited manually to accommodate the convention. Y-axis titles (“Years ka”) were left off as well and were added manually. 5) A red line indicating disturbance in AUG7 was added manually in the published figure. A function is included here to mark the location of the disturbance, but this will look slightly different.