Y-Nb 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 | Nb_ppm, Y_ppm |
require('ggplot2'); classifydiag.granites_y_nb_pearce <- function (wtp_data, ...) { ds <- selectNames(wtp_data,c('Name','Nb_ppm','Y_ppm')) pts <- data.frame(Name=ds$Name, X = ds$Y_ppm, Y = ds$Nb_ppm); labels <- c("Y", "Nb"); lines <- data.frame(X1 = c(40, 1000, 3, 500), Y1 = c(1, 95, 500, 250), X2 = c(50, 50, 50, 30), Y2 = c(11, 11, 11, 25)); marks <- data.frame(Label=c("ORG", "WPG", "VAG +\n syn-COLG"), X=c( 300, 20, 4), Y=c( 3, 90, 40)); ret <- c(); ret$type <- "diagram"; ret$points <- pts; ret$labels <- labels; ret$marks <- marks; ret$lines <- lines; ret$limits <- c(1,1000,1,2000); ret$scale <- c('log','log'); return(ret); }