Yb+Ta-Rb discriminant diagrams for syn-collision granites (syn-COLG), volcanic arc granites (VAG), within plate granites (WPG) and ocean ridge granites (ORG).
| Reference | Pearce, Julian A., Harris, Nigel B. W., Tindle, Andrew G. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks, 1984 [ref] |
| Data | bulk |
| Input | Yb_ppm, Ta_ppm, Rb_ppm |
require('ggplot2');
classifydiag.granites_ybta_rb_pearce <- function (wtp_data, ...)
{
ds <- selectNames(wtp_data,c('Name','Yb_ppm','Ta_ppm','Rb_ppm'))
pts <- data.frame(Name=ds$Name,
X = ds$Yb_ppm + ds$Ta_ppm,
Y = ds$Rb_ppm);
labels <- c("Yb + Ta", "Rb");
lines <- data.frame(X1 = c( 6, 200, 0, 40),
Y1 = c( 0, 450, 105, 2000),
X2 = c( 6, 6, 6, 6),
Y2 = c(200, 8, 200, 200));
marks <- data.frame(Label=c("VAG", "ORG", "syn-COLG", "WPG"),
X=c( 1.05, 40, 1.15, 80),
Y=c( 5 , 7, 1000, 700));
ret <- c();
ret$type <- "diagram";
ret$points <- pts;
ret$labels <- labels;
ret$marks <- marks;
ret$lines <- lines;
ret$limits <- c(0.1,200,1,2000);
ret$scale <- c('log','log');
return(ret);
}