library(“moments”, lib.loc=“C:/Program Files/R/R-3.2.2/library”) surface<- read.csv(“C:/Users/lschifma/Desktop/surface_kernel.csv”)

density(surface$log_Kunsat,na.rm=T)->Kunsat density(surface$log_nswc,na.rm=T)->nswc density(surface$log_ssurgo,na.rm=T)->ssurgo density(surface$log_Texture,na.rm=T)->Texture density(surface$log_SSC,na.rm=T)->SSC density(surface$log_adjSSC,na.rm=T)->adjSSC density(surface$log_SSCBD,na.rm=T)->SSCBD density(surface$log_adjSSCBD,na.rm=T)->adjSSCBD

png(filename=“Surface_Infil_2panel2.png”,width=8,height=9,units=“in”,res=550)

par(mfrow=c(1,2),oma=c(1,0,0,0),tcl=-0.4) par(mar=c(4,4,1,1)) plot(density(surface$log_Kunsat,na.rm=T),main=“ ”,xlab=“ ”, ylab= “ ”,lty=0, lwd=0,xlim=c(-3,3.),ylim=c(0,2.2)) polygon(nswc,lwd=2, border=(col=rgb(31/255, 120/255, 180/255)),col=rgb(31/255, 120/255, 180/255,0.9)) polygon(ssurgo, lwd=2,border=(col=rgb(166/255, 206/255, 227/255,1)),col=rgb(166/255, 206/255, 227/255,0.8)) polygon(Kunsat,lty=1,lwd=2, border=T,col=“black”,density=10) legend(0.5,2.3,cex=0.9, x.intersp=0.5,y.intersp=0.7, legend=c(“Measured”,“NSWC”,“SSURGO”), lty = c(1,1,1),lwd=2, bty='n', col=c(“black”,“#1F78B4” ,“#A6CEE3”)) axis(1,labels=F) axis(2,labels=F) mtext(expression(italic(“Density”)), side=2, outer=F, line=3,at=1)

par(mar=c(4,2,1,1)) plot(density(surface$log_Kunsat,na.rm=T),main=“ ”,yaxt=“n”,ylab=“ ”,xlab=“ ”,lty=0, lwd=0,xlim=c(-3,3.),ylim=c(0,2.2)) polygon(adjSSC, lwd=2,border=(col=rgb(178/255, 223/255, 138/255,1)),col=rgb(178/255, 223/255, 138/255,0.9)) polygon(SSC,lwd=2,border=(col=rgb(51/255,160/255,44/255,1)),col=rgb(51/255,160/255,44/255,0.8)) polygon(adjSSCBD, lwd=2,border=(col=rgb(251/255, 154/255, 153/255,1)),col=rgb(251/255, 154/255, 153/255,0.7)) polygon(SSCBD, lwd=2,border=(col=rgb(227/255, 26/255, 28/255,1)),col=rgb(227/255, 26/255, 28/255,0.6)) polygon(Texture, lwd=2,border=(col=rgb(253/255, 191/255, 30/255,1)), col=rgb(253/255, 191/255, 50/255,0.5)) polygon(Kunsat,lty=1,lwd=2, border=T,col=“black”,density=10) legend(-0.5,2.3,cex=0.9, x.intersp=0.7,y.intersp=0.7, legend=c(“Measured”,“Texture Class”,“SSC”, “Adjusted SSC”, “SSC+BD”,“Adjusted SSC+BD”), lty = c(1,1,1,1,1,1),lwd=2, bty='n', col=c(“black” ,“#FDBF1E”,“#33A02C” , “#B2DF8A” ,“#E31A1C”,“#FB9A99”)) axis(1,labels=F) axis(2,labels=F) text(-0.52, 1.25, srt=0, expression(italic(“loam”)),cex=0.75,col=“#ca9818”) arrows(-0.47, 1.2,-0.8, 1, angle=20,length=0.15,code = 2, col = “gray22”, lty = 1,lwd = 1) text(0.46, 1.5,srt=0, expression(italic(“sandy loam”)),cex=0.75,col=“#ca9818”)

mtext(expression(italic(“Log Infiltration (cm/hr)”)), side=1, outer=F, line=3,at=-4)

dev.off()

skewness(surface$Kunsat,na.rm=T) skewness(surface$nswc,na.rm=T) skewness(surface$ssurgo,na.rm=T) skewness(surface$Texture,na.rm=T) skewness(surface$SSC,na.rm=T) skewness(surface$adjSSC,na.rm=T) skewness(surface$SSCBD,na.rm=T) skewness(surface$adjSSCBD,na.rm=T)

kurtosis(surface$Kunsat,na.rm=T) kurtosis(surface$nswc,na.rm=T) kurtosis(surface$ssurgo,na.rm=T) kurtosis(surface$Texture,na.rm=T) kurtosis(surface$SSC,na.rm=T) kurtosis(surface$adjSSC,na.rm=T) kurtosis(surface$SSCBD,na.rm=T) kurtosis(surface$adjSSCBD,na.rm=T)

####################################################################

subsurface_kernel<- read.csv(“C:/Users/lschifma/Desktop/subsurface_kernel.csv”) View(subsurface)

density(subsurface_kernel$log_Ksat,na.rm=T)->Sub_Ksat density(subsurface_kernel$log_Texture,na.rm=T)->Sub_Texture density(subsurface_kernel$log_SSC,na.rm=T)->Sub_SSC density(subsurface_kernel$log_adjSSC,na.rm=T)->Sub_adjSSC

plot(Sub_Ksat) plot(Sub_Texture) plot(Sub_SSC) plot(Sub_adjSSC)

png(filename=“Subsurface_Drain_1panel2.png”,width=8,height=9,units=“in”,res=550)

par(mfrow=c(1,1),oma=c(1,0,0,0),tcl=-0.5) par(mar=c(4,4,1,1)) plot(density(subsurface_kernel$log_Ksat,na.rm=T),main=“ ”,xlab=“ ”,ylab=“ ”,lty=0, lwd=0,xlim=c(-3,3.),ylim=c(0,2.2)) polygon(Sub_adjSSC, lwd=2,border=(col=rgb(178/255, 223/255, 138/255,1)),col=rgb(178/255, 223/255, 138/255,0.7)) polygon(Sub_SSC,lwd=2,border=(col=rgb(51/255,160/255,44/255,1)),col=rgb(51/255,160/255,44/255,0.7)) polygon(Sub_Texture, lwd=2,border=(col=rgb(253/255, 191/255, 30/255,1)), col=rgb(253/255, 191/255, 50/255,0.4)) polygon(Sub_Ksat,lty=1,lwd=2, border=T,col=“black”,density=10) axis(1,labels=F) axis(2,labels=F) legend(cex=0.9, x.intersp=0.5,y.intersp=0.7,“topright”, legend=c(“Measured”,“Texture Class”,“SSC”, “Adjusted SSC”), lty = c(1, 1,1,1,1,1,1,1),lwd=2, bty='n', col=c(“black” ,“#FDBF1E”,“#33A02C” , “#B2DF8A”)) text(1.55, 0.3, expression(italic(“sand”)),cex=0.75,col=“#ca9818”) text(1.0, 0.6, expression(italic(“loamy sand”)),cex=0.75,col=“#ca9818”) text(0.65, 1, expression(italic(“sandy loam”)),cex=0.75,col=“#ca9818”) text(0.2, 1.4, expression(italic(“loam\nsandy clay loam”)),cex=0.75,col=“#ca9818”) mtext(expression(italic(“Log Drainage (cm/hr)”)), side=1, outer=F, line=3,at=0) mtext(expression(italic(“Density”)), side=2, outer=F, line=3,at=1)

dev.off()