Гранат-клинопироксеновый термометр Д. Накамура, 2009

Это старая версия документа!

Термометр с учетом неидеальности Fe-Mg обмена в гранате.

Подробнее

Информация о публикации	D. Nakamura, A new formulation of garnet–clinopyroxene geothermometer based on accumulation and statistical analysis of a large experimental data set, J. metamorphic Geol., 2009, 27, 495–508
Ссылка	nakamura2009.pdf
Калибровка	Экспериментальные данные. 333 пары в графитовых капсулах
Точность	±74 ºС
Диапазон значений	800-1820 ºС
Обновляет сенсоры	grt_cpx_krogh_ravna

Код метода

source('recalc_grt_8cat_norm.r');
source('recalc_cpx_4cat_norm_jd.r');

sensors.t.grt_cpx_nakamura <- function(Grt, Cpx, P, usePFU = FALSE, ...){

    if (!usePFU) Grt <- recalc.grt_8cat_norm(Grt);
    if (!usePFU) Cpx <- recalc.cpx_4cat_norm_jd(Cpx);
    

    X_Grt_Ca <- Grt$Ca_pfu / (Grt$Mg_pfu + Grt$Fe_pfu + Grt$Ca_pfu);
    X_Grt_Fe <- Grt$Fe_pfu / (Grt$Mg_pfu + Grt$Fe_pfu + Grt$Ca_pfu);
    X_Grt_Mg <- Grt$Mg_pfu / (Grt$Mg_pfu + Grt$Fe_pfu + Grt$Ca_pfu);
    
    X_Cpx_Fe <- Cpx$Fe_p2_pfu / (Cpx$Al_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + Cpx$Mg_pfu);
    X_Cpx_Mg <- Cpx$Mg_pfu / (Cpx$Al_pfu + Cpx$Fe_p2_pfu + Cpx$Fe_p3_pfu + Cpx$Mg_pfu);
    
    Kd <- (Grt$Fe_pfu / Grt$Mg_pfu) / (Cpx$Fe_p2_pfu / Cpx$Mg_pfu);
    
    X_Prp <- Grt$Mg_pfu / (Grt$Fe_pfu + Grt$Mn_pfu + Grt$Mg_pfu + Grt$Ca_pfu);
    X_Alm <- Grt$Fe_pfu / (Grt$Fe_pfu + Grt$Mn_pfu + Grt$Mg_pfu + Grt$Ca_pfu);
    X_Sps <- Grt$Mn_pfu / (Grt$Fe_pfu + Grt$Mn_pfu + Grt$Mg_pfu + Grt$Ca_pfu);
    X_Grs <- Grt$Ca_pfu / (Grt$Fe_pfu + Grt$Mn_pfu + Grt$Mg_pfu + Grt$Ca_pfu);
    
    A <- 0.5*X_Grs*(X_Prp - X_Alm - X_Sps);
    B <- 0.5*X_Grs*(X_Prp - X_Alm + X_Sps);
    C <- 0.5*(X_Grs + X_Sps)*(X_Prp - X_Alm);
    
    
    Tcalc <- (  2784 + 14.52 * P +
               (2610 + 1.44*P)*(2*X_Grs*X_Prp - A) +
               (1183 + 6.98*P)*(X_Grs^2 - A) -
               105*(2*X_Grs*X_Alm + B) +
               (814.6 + 3.61*P)*(X_Grs^2 + B) -
               (254.6 + 8.42*P)*(2*X_Prp*X_Alm - X_Alm^2 + C) -
               83.6*(X_Prp^2 - 2*X_Prp*X_Alm + C) + 
               1388*X_Sps - 
               462*(X_Cpx_Mg - X_Cpx_Fe)  )  /
             ( log(Kd) + 1.431 + 0.695*(2*X_Grs*X_Prp + X_Grs^2 - 2*A) +
               0.203 * (X_Grs^2 - 2*X_Grs*X_Alm) +
               0.922 * X_Sps ) -
              273;
    
    return (data.frame(Name=paste0(Grt$Name,'-',Cpx$Name),P = P,A = A*100, B=B*100,C=C*100, T_calc_Nakamura = Tcalc));
}

sensors.t.grt_cpx_nakamura.args <- function() {
    ars <- list();
    ars$minerals = c('Grt','Cpx');
    ars$ranges = c('P');
    return(ars);
}

sensors.t.grt_cpx_nakamura.test <- function () {
  
  Grt <- data.frame(Name='Grt',SiO2_wtp=37.50,TiO2_wtp=0.03, Al2O3_wtp=20.21, FeO_wtp=26.3, MnO_wtp=2.94, MgO_wtp=3.25, CaO_wtp=9.45);
  Cpx <- data.frame(Name='Cpx',SiO2_wtp=50.95,TiO2_wtp=0.23, Al2O3_wtp=2.42, FeO_wtp=12.35, MnO_wtp=0.33, MgO_wtp=10.93, CaO_wtp=21.39, Na2O=0.62);

  result <- sensors.t.grt_cpx_krogh_ravna (Grt, Cpx, 10);
  
  expect_equal(result, 711, tolerance = 1);
}