Aplicación de las pruebas estadísticas de discordancia y significancia en la comparación del vulcanismo dacítico de la parte central de Cinturón Volcánico Mexicano

Application of discordancy and significance statistical tests for the comparison of dacitic volcanism from the central part of the Mexican Volcanic Belt

Lorena Díaz-González¹ and René Cruz-Huicochea²

¹ Facultad de Ciencias, Universidad Autónoma de Estado de Morelos. E-mail: ldg@uaem.mx.

Recepción: 29-10-2012
Aceptación: 29-10-2013

Resumen

Nuestro objetivo es presentar una metodología estadística, junto con dos nuevos programas (DODESSYS y UDASYS). Para esta tarea compilamos una base de datos de 249 muestras de rocas dacitas provenientes de cuatro regiones del cinturón volcánico mexicano (MVB): volcanes monogenéticos de la Sierra de Chichinautzin y el Valle de México, estratovolcán Nevado de Toluca, estratovolcán Iztaccíhuatl y estratovolcán Popocatépetl. Las pruebas estadísticas de discordancia y significancia (ANOVA –ANalysis Of Variance–, F de Fisher y t de Student) fueron aplicadas al 99% de nivel de confianza. Se calculó la estadística final para 98 parámetros, incluyendo óxidos mayores, elementos de tierras raras, elementos traza y parámetros adicionales, tales como parámetros de relaciones logarítmicas usados en nuevos diagramas de discriminación tectónica. Estos parámetros fueron tratados como muestras estadísticas univariadas y fueron clasificados en cuatro regiones del MVB. Las pruebas estadísticas de discordancia detectaron datos discordantes en 124 (aproximadamente en el 35%) muestras estadísticas univariadas. La prueba ANOVA mostró diferencias significativas entre todos los grupos en 32 parámetros. Las similitudes y diferencias entre los parámetros de relaciones logarítmicas pueden ser útiles en el futuro para proponer diagramas de discriminación tectónica a partir de una base de datos representativa.

Palabras clave: ANOVA, F de Fisher, t de Student, datos discordantes, datos geoquímicos, manejo estadístico de datos composicionales.

Abstract

Keywords: ANOVA, Fisher's F, Student's t, discordant outliers, geochemical data, statistical handling of compositional data.

Introduction

Recently, a new computer programs has been developed, UDASYS (Univariate Data Analysis SYStem) [1]. UDASYS is freely available from any of the authors to any scientist interested in correctly processing experimental data. This program, written in Java [2], provides statistical tools pertaining to both robust and outlierbased methods for univariate data. UDASYS also incorporates an updated version of the DODESSYS software [3]. Whereas DODESSYS allowed the application of discordancy tests ([3] for more details on these tests see Table S1 in online Supplementary Material) for statistical sample sizes up to 1000, all discordancy tests can now be applied to statistical sample sizes as large as 30000. Computer programs to enable the application of discordancy tests were practically nonexistent as documented by Barnett and Lewis (1994). Later about 12 years ago, a computer program (SIPVADE) was published by Verma et al. (1998), but it is now outdated for several reasons. The most important among them are that SIPVADE uses old, less precise and sometimes even inaccurate critical values then available in the literature (Barnett and Lewis 1994; Verma 2005) and relies on linear interpolation of these values when for a given statistical sample size n, the corresponding critical values were not tabulated. Both of these aspects have been shown to cause errors in the final statistical inferences. More importantly, unlike all available software to date (e.g., [4]), UDASYS allows a highly efficient use of significance tests of Fisher's F, Student's t, and ANOVA.

This work illustrates the application of statistical discordancy and significance tests using geochemical data. A geochemical database of major-elements in rocks from the Mexican volcanic belt (MVB) was established long ago by Pal et al. [5]. These authors used their database to objectively characterise for the first time the nature of volcanism in the MVB. This work was later extended by including more analyses of MVB rocks in this database which permitted to highlight the complexity of magmas in the MVB (e.g., [6]). Mean and standard deviation estimates of compositional data were presented by these authors, but this was done without the application of discordancy tests [7]. Similarly in local geochemical studies from this volcanic province (MVB), these two statistical parameters for laboratory analytical data were specifically reported by Verma [8-10] and Verma et al. [11]. Other researchers have used mean and standard deviation estimates for geochemical interpretation [12].

In this work geochemical data are compiled for dacitic rocks from four nearby areas of the MVB. The geochemical parameters are compared through the significance tests such as Fisher's F and Student's t [13] without and with the application of discordancy tests [14-17]. The results highlight the importance of these statistical tests in geosciences.

We searched the published geoscience literature for specific applications of discordancy and significance tests and found that it is not a common practice to apply them in geoscientific studies. Below we list some the reports found that made use of either one of these statistical methodologies.

Rice and Church [18] presented a statistical study on the variability in grain size of sediment from two confluent rivers in northeastern British Columbia (Canada). They stated that it was not appropriate to apply ANOVA test because the statistical samples did not show a normal distribution and their variances were unequal. However, the validity of the first condition can be checked by discordancy tests, whereas the second condition (equal variances) is not a requisite for ANOVA. They applied tests, such as Brown-Forsythe and chi-square, for comparing statistical sample means when sample variances are unequal.

Takano et al. [19] made a statistical comparison of inter-laboratory analytical data of fluid samples from crater lake of Maly Semiachik volcano, located in the central part of the Eastern Volcanic Belt of Kamchatka (Japan), obtained from eight different institutions. They used different analytical techniques (atomic absorption spectrometry, atomic emission spectrometry, mass spectrometry, ion chromatography, high performance liquid chromatography, colorimetry, and titrimetry) to compare the measured isotopic data coming from elements such as hydrogen, sulfur, and oxygen. Their comparison consisted of simply calculating the central tendency (mean) and dispersion (coefficient of variation) parameters for each one of these techniques. Experience shows that it would have been worthwhile to apply the discordancy and significance tests for such inter-laboratory evaluations as suggested earlier by several authors [20-21].

The correct statistic application, such as the work we are reporting, tries to promote the evolution of geochemistry towards geochemometrics, where statistics is an essential part of experimental data treatment. In general, in the area of geochemistry is not customary to apply a correct statistics methodology in the processing of databases. For example, Takano et al. [19] assessed the statistic differences in their experimental databases, but failed to apply the methodology based on significance tests and discordancy. However, recently some authors applied successfully this complete methodology in processing geochemical data [17, 23]. Particularly, the discordance tests have been applied in a diversity of scientific and engineering fields, including some branches of earth sciences such as determination of Nernst distribution coefficients [24]; quality control through reference materials [14-16, 23]; geothermal research [25-27]; geochemistry [12, 15, 17]; volcanoes studies [28, 29]; pollution studies [30]; petroleum research [31]; soil research [32]; proteomics research [33]. Also, sensitivity and uncertainty analysis is another important statistical application [34-38].

Method

Database and procedures

Geochemical data for 249 Neogene dacitic rock samples from four closely located areas of the MVB were compiled. The literature sources were as follows: [9, 21, 39-60]. Data are identifiqued as group numbers Gr1 to Gr4 (see locations of these regions on a map presented in Figure 1): Region 1 (Gr1)–diverse locations of the Sierra de Chichinautzin and the southern of the Valle de México (monogenetic volcanoes); Region 2 (Gr2)–the Nevado de Toluca stratovolcano; Region 3 (Gr3)–the Iztaccíhuatl stratovolcano, and Region 4 (Gr4)–the Popocatépetl stratovolcano.

TAS (Total Alkalis vs Silica) diagram was generated by IgRocs sofware [61]; see Figure 2. Geochemical data are concentrated in classification area for dacite rocks.

The statistical central tendency (mean) and dispersion (standard deviation) parameters were calculated for several conventional variables, which were 11 major oxides (adjusted values) from (SiO₂)_Adj to (P₂O₅)_Adj, selected normative minerals, Mg-number (or Mg-value), and 6 other indices detailed by [61], followed by 14 rare earth elements from La to Lu, and 22 trace elements from Ba to Zr. In addition to these conventional chemical data, 30 additional parameters were computed and evaluated. These include two ratio parameters defined by Verma [62] called Nb-anomaly with respect to Ba and La and Ta-anomaly with respect to Ba and La, as well as 28 log-ratio parameters of elements used in new multi-dimensional tectonic discrimination diagrams [63-65].

Figure 3 shows the flow diagram of statistical methodology applied in this work. Conventionally, significant test are applied without prior application of discordancy tests. However, because these tests should be applied to normally distributed statistical samples, data for each variable from all individual groups (Gr1-Gr4) were first processed for discordant outliers by single-outlier type discordancy tests (see Table S1 in [1]) at 99% confidence level, and the discordant outlier-free groups were evaluated from the two-sided ANOVA-test and t-test at 99% confidence level (see Geological implications in [1] for more details on application of two-sided version of significant tests). The statistical parameters of mean and standard deviation were simply calculated from the discordant outlier-free individual groups.

We note that, ANOVA test can only be applied to three or more groups or statistical samples [7], therefore this significant test was applied to the data from each group (Gr1-Gr4). The application of ANOVA would result in any of the following: (i) no statistically significant differences among the four regions (Gr1-Gr4); (ii) one –e.g., Gr1– of the four regions showing a statistically significant difference as compared to the other three regions –e.g., Gr2- Gr4–; and (iii) statistically significant differences among the four regions (Gr1- Gr4), which will have to be resolved by Fisher's F and Student's t significance tests. When ANOVA detects significant differences among the four regions, the data should be processed thorough the combination of Fisher's F and Student's t tests, which are applicable to only two groups at a time [63, 64]. The Fisher's F test compares the two variances and could result in either the two variances are equal or the two are different. Depending on the result of the F test, the appropriate version of the t test should be applied. The Fisher's F and Student's t tests were applied to each one of the combinations Gr1-Gr2, Gr1-Gr3, Gr1-Gr4, Gr2-Gr3, Gr2-Gr4 and Gr3-Gr4.

Resultados

Identification and separation of discordant outliers

Geochemical data for a total of 96 variables o parameters from each group (Gr1-Gr4) were processed in this work. Single-outlier type discordancy tests at a very strict 99% confidence level were then applied to individual groups, outlying observations were separated, and statistical parameters were calculated from discordant outlier-free data. These statistical parameters are reported in Table 1; the first column gives the name of the chemical or ratio parameter, the next columns gives statistical parameters such as statistical sample size (n), mean and standard deviation from all individual groups (Gr1-Gr4); i.e. columns 2-4 show stastistical parameters from Sierra de Chichinautzin and Valle de Mexico monogenetic volcanoes. The second column gives the final statistical sample size (n) after discordant outlier detection and separation, the third column reports the mean, and the fourth one provided the standard deviation. The number of discordant outliers is represented by a symbol as superscript: α –one–; β –two–; γ –three–; δ –four –; £ –five–; ζ –seven–; η –eight–; λ –ten–. For the total of 350 statiscal samples processed in this work, 124 (35%) samples showed discordant outliers.

Application of ANOVA, t and F tests after elimination of outliers

ANOVA test determined that 32 variables did not show statistically significant differences among all groups, hence they were combined; e.g., (Na₂O)_Adj, (K₂O)_Adj, or_Norm, ab_Norm, an_Norm , La, Pr, Nd, Sm, Eu, ln((Na₂O)_Adj/Si), ln((K₂O)_Adj/Si), ln((P₂O₅)_Adj/Si), ln(Nb/Yb), ln(Th/Yb), ln(Y/Yb) and ln(Zr/Yb). ANOVA also identified a discordant group (Gr2, Gr3 and Gr4 in 3, 17 and 12 variables, respectively) in 32 variables; e.g., the Gr2 group was identified as discordant group in (P₂O₅)_Adj variable, therefore, Gr2 group was separated and Gr1, Gr2 and Gr4 groups were combined. Finally, ANOVA determined statistically significant differences among the four regions in 32 elements, e.g., all groups from (TiO₂)_Adj major element were identified as discordant groups and were not combined. Fisher's F and Student's t tests were applied to these 32 variables.

Application of discordancy tests after combining data (significance tests)

Single-outlier type discordancy tests were applied to the combined groups, outlying observations were separated, and statistical parameters were calculated from discordant outlier-free data (see Table 2 in appendix). These discordant outliers were rejected (or separated) and final statistical were calculated and shown in Table 2. Discordant outliers were represented by a symbol as superscript: α –one–; β –two–; γ –three–; δ –four –; £ –five–; ζ –seven–; η –eight–; λ –ten–.

In this work, we have shown a statistical procedure to decipher mean compositions and uncertainty estimates including various regions. For this, geochemical data are compiled for 249 Neogene dacitic rock samples from the four MVB regions.

All single-outlier type discordancy tests and significance (ANOVA –ANalysis Of Variance–, Fishers' F and Student's t) statistical tests were applied at the very strict 99% confidence level. These statistical tests were applied to each one of the 98 geochemical parameters, which were major oxides, selected normative minerals, rare earth, trace elements, two ratio parameters called Nb-anomaly and Ta-anomaly, as well as 28 log-ratio parameters of elements used in new multi-dimensional tectonic discrimination diagrams.

All geochemical parameters were treated as univariate statistical samples. Final statistical parameters were calculated from discordant outlier-free data. We suggest that the final mean compositions could be used to compare statistically the geochemical data for the same type of igneous rocks, i.e., dacite type, sampled around the world.

Furthermore, significance statistical tests determined significant differences and similarities among various geochemical parameters from the four MVB regions. Particularly, the similarities and differences among the log-ratios parameters could be useful to propose new diagrams to discriminate tectonic settings, with a more representative database.

Acknowledgements

L.D-G. acknowledges PROMEP support to the project "Estadística Computacional para el tratamiento de datos experimentales" (PROMEP/103-5/10/7332). R.C-H. acknowledges CONACyT-Mexico for a scholarship to carry out his master studies.

Referencias

[2] Deitel P J and Deitel H M, Java: How to Program, 7 Ed., Prentice Hall, Java: How to Program, 2006.         [ Links ]

[3] Verma S P and Díaz-González L, Application of the discordant outlier detection and separation system in the geosciences, International Geology Review, Vol. 54, No. 5, 2012, pp. 593-614.         [ Links ]

[4] StatSoft I. STATISTICA (data analysis software system), statsoft.com, 2004.         [ Links ]

[5] Pal S, et al., Magma characterization of the Mexican volcanic belt (Mexico), Bulletin of Volcanology, Vol. 41, No. 4, 1978, pp. 379-389.         [ Links ]

[6] Verma S P, Aguilar-Y-Vargas, V.H., Bulk chemical composition of magmas in the Mexican Volcanic Belt (Mexico) and inapplicability of generalized arc-models, Chemie der Erde, Vol. 48, 1988, pp. 203-221.         [ Links ]

[7] Barnett V and Lewis T, Outliers in statistical data, J. Wiley, 1994.         [ Links ]

[8] Verma S P, Alkali and alkaline earth element geochemistry of Los Humeros Caldera, Puebla, Mexico, Journal of Volcanology and Geothermal Research, Vol. 20, No. 1-2, 1984, pp. 21-40.         [ Links ]

[9] Verma S P, Geochemistry of evolved magmas and their relationship to subduction-unrelated mafic volcanism at the volcanic front of the central Mexican Volcanic Belt, Journal of Volcanology and Geothermal Research, Vol. 93, No. 1-2, 1999, pp. 151-171.         [ Links ]

[10] Verma S P, Geochemical evidence for a lithospheric source for magmas from Los Humeros caldera, Puebla, Mexico, Chemical Geology, Vol. 164, No. 1-2, 2000, pp. 35-60.         [ Links ]

[11] Verma S P, Carrasco-N'uñez G, and Milán M, Geology and geochemistry of Amealco Caldera, Qro., Mexico, Journal of Volcanology and Geothermal Research, Vol. 47, No. 1-2, 1991, pp. 105-127.         [ Links ]

[12] Madhavaraju J L, Y. I., Geochemistry of the Dalmiapuram Formation of the Uttatur Group (Early Cretaceous), Cauvery basin, southeastern India: Implications on provenance and paleo-redox conditions, Revista Mexicana de Ciencias Geológicas, Vol. 26, No. 2, 2009, pp. 380-394.         [ Links ]

[13] Miller J N and Miller J C, Statistics and chemometrics for analytical chemistry, Fifth Ed., Pearson Prentice Hall, 2005.         [ Links ]

[14] Marroquín-Guerra S G, Velasco-Tapia F, and Díaz-González L, Evaluación estadística de Materiales de Referencia Geoquímica del Centre de Recherches Pétrographiques et Géochimiques (Francia) aplicando un esquema de detección y eliminación de valores desviados, Revista mexicana de ciencias geológicas, Vol. 26, 2009, pp. 530-542.         [ Links ]

[15] Pandarinath K, Evaluation of geochemical sedimentary reference materials of the Geological Survey of Japan (GSJ) by an objective outlier rejection statistical method, Revista Mexicana de Ciencias Geológicas, Vol. 26, No. 3, 2009b, pp. 638-646.         [ Links ]

[16] González-Ramírez R, Díaz-González, L., Verma, S.P., Eficiencia relativa de 15 pruebas de discordancia con 33 variantes aplicadas al procesamiento de datos geoquímicos, Revista Mexicana de Ciencias Geológicas, Vol. 26, No. 2, 2009, pp. 501-515.         [ Links ]

[17] Armstrong-Altrin J S, Provenance of sands from Cazones, Acapulco, and Bahía Kino beaches, México, Revista Mexicana de Ciencias Geológicas, Vol. 26, No. 3, 2009, pp. 764-782.         [ Links ]

[18] Rice S and Church M, Grain size along two gravel-bed rivers: statistical variation, spatial pattern and sedimentary links, Earth Surface Processes and Landforms, Vol. 23, No. 4, 1998, pp. 345-363.         [ Links ]

[19] Takano B, Fazlullin S M, and Delmelle P, Analytical laboratory comparison of major and minor constituents in an active crater lake, Journal of Volcanology and Geothermal Research, Vol. 97, No. 1-4, 2000, pp. 497-508.         [ Links ]

[20] Guevara M, Verma, S.P., Velasco-Tapia, F., Evaluation of GSJ intrusive rocks JG1, JG2, JG3, JG1a, and Jgb1 by an objective outlier rejection statistical procedure, Revista Mexicana de Ciencias Geológicas, Vol. 18, No. 1, 2001, pp. 74-88.         [ Links ]

[21] Velasco-Tapia F, Verma, S.P., Magmatic Processes at the volcanic front of Central Mexican Volcanic Belt: Sierra de Chichinautzin volcanic field (Mexico), Turkish Journal of Earth Sciences, No. in press, 2011.         [ Links ]

[22] Wani H and Mondal M E A, Evaluation of provenance, tectonic setting, and paleoredox conditions of the Mesoproterozoic-Neoproterozoic basins of the Bastar craton, Central Indian Shield: Using petrography of sandstones and geochemistry of shales, Lithosphere, Vol. 3, No. 2, 2011, pp. 143-154.         [ Links ]

[23] Pandarinath K, Clay minerals in SW Indian continental shelf sediment cores as indicators of provenance and palaeomonsoonal conditions: a statistical approach, International Geology Review, Vol. 51, No. 2, 2009, pp. 145-165.         [ Links ]

[24] Torres-Alvarado I S, et al., DC_BASE: a database system to manage Nernst distribution coefficients and its application to partial melting modeling, Computers & Geosciences, Vol. 29, No. 9, 2003, pp. 1191-1198.         [ Links ]

[25] Díaz-González L, Santoyo E, and Reyes-Reyes J, Tres nuevos geotermómetros mejorados de Na/K usando herramientas computacionales y geoquimiométricas: aplicación a la predicción de temperaturas de sistemas geotérmicos, Revista Mexicana de Ciencias Geológicas, Vol. 25, No. 3, 2008, pp. 465-482.         [ Links ]

[26] Gómez-Arias E, et al., Determinación de la viscosidad y su incertidumbre en fluidos de perforación usados en la construcción de pozos geotérmicos: aplicación en el campo de Los Humeros, Puebla, México, Revista mexicana de ciencias geológicas, Vol. 26, 2009, pp. 516-529.         [ Links ]

[27] Pandarinath K, Applicability of solute geothermometry for springs and wells of the Los Azufres and Las Tres Vírgenes geothermal fields, Mexico, International Geology Review, No. in press, 2011.         [ Links ]

[28] Sanci R, Panarello, H.O., and Ostera, H.A., Flujo de dióxido de carbono en el flanco oriental del Volcán Peteroa Andes del Sur, Revista Mexicana de Ciencias Geológicas, Vol. 27, No. 2, 2010, pp. 225-237.         [ Links ]

[29] Torres-Alvarado I S, Smith A D, and Castillo-Román J, Sr, Nd and Pb isotopic and geochemical constraints for the origin of magmas in Popocatépetl volcano (central Mexico) and their relationship with the adjacent volcanic fields, International Geology Review, Vol. 53, No. 1, 2011, pp. 84-115.         [ Links ]

[30] Heath E, et al., Second interlaboratory exercise on non-steroidal anti-inflammatory drug analysis in environmental aqueous samples, Talanta, Vol. 81, No. 4-5, 2010, pp. 1189-1196.         [ Links ]

[31] Salleh H S, Rosales, E., and Flores-de la Mota, I., Influence of different probability based models on oil prospect exploration decision making: a case from southern Mexico, Revista Mexicana de Ciencias Geológicas, Vol. 24, No. 3, 2007, pp. 306-317.         [ Links ]

[32] González-Márquez L C and Hansen A M, Adsorción y mineralización de atrazina y relación con parámetros de suelos del DR 063 Guasave, Sinaloa, Revista mexicana de ciencias geológicas, Vol. 26, 2009, pp. 587-599.         [ Links ]

[33] Viner R I, et al., Quantification of post-translationally modified peptides of bovine ?-crystallin using tandem mass tags and electron transfer dissociation, Journal of Proteomics, Vol. 72, No. 5, 2009, pp. 874-885.         [ Links ]

[34] Santoyo E, et al., Rheological property measurement of drilling fluids used in geothermal wells, Applied Thermal Engineering, Vol. 21, No. 3, 2001, pp. 283-302.         [ Links ]

[35] Santoyo E, et al., Convective heat-transfer coefficients of non-Newtonian geothermaldrilling fluids, Journal of Geochemical Exploration, Vol. 78-79, No. 0, 2003, pp. 249-255.         [ Links ]

[36] Santoyo E, Guevara M, and Verma S P, Determination of lanthanides in international geochemical reference materials by reversed-phase high-performance liquid chromatography using error propagation theory to estimate total analysis uncertainties, Journal of Chromatography A, Vol. 1118, No. 1, 2006, pp. 73-81.         [ Links ]

[37] Santoyo E, et al., Separation and quantification of lanthanides in synthetic standards by capillary electrophoresis: A new experimental evidence of the systematic "odd-even" pattern observed in sensitivities and detection limits, Journal of Chromatography A, Vol. 1149, No. 1, 2007, pp. 12-19.         [ Links ]

[38] Espinosa-Paredes G, et al., Mass flow rate sensitivity and uncertainty analysis in natural circulation boiling water reactor core from Monte Carlo simulations, Nuclear Engineering and Design, Vol. 240, No. 5, 2010, pp. 1050-1062.         [ Links ]

[39] Gunn B M, Mooser, F., Geochemistry of the volcanics of central Mexico, Bulletin Volcanologique, Vol. 34, No. 2, 1971, pp. 577-613.         [ Links ]

[40] Negendank J F W, Volcanics of the Valley of Mexico. Description of some Mexican volcanic rocks with special consideration of the opaques. Part I: petrography of the volcanics, Neues Jahrbuch für Mineralogie-Abhandlungen, Vol. 116, No. 3, 1972, pp. 308-320.         [ Links ]

[41] Bloomfield K, A late-Quaternary monogenetic volcano field in central Mexico, Geologische Rundschau, Vol. 64, No. 1, 1975, pp. 476-497.         [ Links ]

[42] Pérez R. J, Pal, S., Terrell, D.J., Urrutia F., J., López M., M., Preliminary report on the analysis of some "in-house" geochemical reference samples from Mexico, Geofísica Internacional, 18, Vol. 2, 1979, pp. 197-209.         [ Links ]

[43] Robin C, Le Volcan Popocatepetl (Mexique): structure, evolution pétrologique et risques, Bulletin of Volcanology, Vol. 47, No. 1, 1984, pp. 1-23.         [ Links ]

[44] Boudal. Pétrologie d'un grand volcan andésitique mexicain: le Popocatépetl, Univ. Clermont-Ferrand II., 1985, pp. 140.         [ Links ]

[45] Nixon G T, Petrology of the Younger Andesites and Dacites of Iztaccíhuatl Volcano, Mexico: I. Disequilibrium Phenocryst Assemblages as Indicators of Magma Chamber Processes, Journal of Petrology, Vol. 29, No. 2, 1988, pp. 213-264.         [ Links ]

[46] Martin del Pozzo A L. Geoquímica y paleomagnetismo de la Sierra Chichinautzin, U.N.A.M., 1989.         [ Links ]

[47] Swinamer R T. The geomorphology, petrography, geochemistry and petrogenesis of the volcanic rocks in the Sierra del Chichinautzin, Mexico, Queen's University., 1989.         [ Links ]

[48] Nixon G T, The geology of Iztaccíhualtl volcano and adjacent areas of the Sierra Nevada and Valley of Mexico, Geological Society of America, Vol. Special paper, 1993, pp. 219, 1-58.         [ Links ]

[49] Larocque A C L, Stimac J A, and Siebe C, Metal-residence sites in lavas and tuffs from Volcán Popocatépetl, Mexico: implications for metal mobility in the environment, Environmental Geology, Vol. 33, No. 2, 1998, pp. 197-208.         [ Links ]

[50] Siebe C, Schaaf P, and Urrutia-Fucugauchi J, Mammoth bones embedded in a late Pleistocene lahar from Popocatépetl volcano, near Tocuila, central México, Geological Society of America Bulletin, Vol. 111, No. 10, 1999, pp. 1550-1562.         [ Links ]

[51] Wallace P J and Carmichael I S E, Quaternary volcanism near the Valley of Mexico: implications for subduction zone magmatism and the effects of crustal thickness variations on primitive magma compositions, Contributions to Mineralogy and Petrology, Vol. 135, No. 4, 1999, pp. 291-314.         [ Links ]

[52] García-Palomo A, Macías, J.L., Arce, J.L., Capra, L., Garduño, V.H., Espíndola, J.M.. Geology of Nevado de Toluca volcano and surrounding areas, Central Mexico, Geological Society of America, 2002, pp. 26.         [ Links ]

[53] Arce J L, Macías J L, and Vázquez-Selem L, The 10.5 ka Plinian eruption of Nevado de Toluca volcano, Mexico: Stratigraphy and hazard implications, Geological Society of America Bulletin, Vol. 115, No. 2, 2003, pp. 230-248.         [ Links ]

[54] Martínez-Serrano R G, et al., Sr, Nd and Pb isotope and geochemical data from the Quaternary Nevado de Toluca volcano, a source of recent adakitic magmatism, and the Tenango Volcanic Field, Mexico, Journal of Volcanology and Geothermal Research, Vol. 138, No. 1-2, 2004, pp. 77-110.         [ Links ]

[55] Siebe C, et al., Geochemistry, Sr-Nd isotope composition, and tectonic setting of Holocene Pelado, Guespalapa and Chichinautzin scoria cones, south of Mexico City, Journal of Volcanology and Geothermal Research, Vol. 130, No. 3-4, 2004, pp. 197-226.         [ Links ]

[56] Schaaf P, et al., Geochemical Evidence for Mantle Origin and Crustal Processes in Volcanic Rocks from Popocatépetl and Surrounding Monogenetic Volcanoes, Central Mexico, Journal of Petrology, Vol. 46, No. 6, 2005, pp. 1243-1282.         [ Links ]

[57] Meriggi L, Macías, J.L., Tommasini, S., Capra, L., Conticelli, S., Heterogeneous magmas of the Quaternary Sierra Chichinautzin volcanic field (central Mexico): the role of an amphibole-bearing mantle and magmatic evolution processes, Revista Mexicana de Ciencias Geológicas, Vol. 25, 2008, pp. 197-216.         [ Links ]

[58] Straub S M and Martin-Del Pozzo A L, The significance of phenocryst diversity in tephra from recent eruptions at Popocatepetl volcano (central Mexico), Contributions to Mineralogy and Petrology, Vol. 140, No. 4, 2001, pp. 487-510.         [ Links ]

[59] Straub S M, et al., Evidence from high-Ni olivines for a hybridized peridotite/pyroxenite source for orogenic andesites from the central Mexican Volcanic Belt, Geochem. Geophys. Geosyst., Vol. 9, No. 3, 2008.         [ Links ]

[60] Velasco-Tapia F, Guevara, M., Verma, S.P., Evaluation of concentration data in geochemical reference materials, Chemie der Erde, Vol. 61, No. 1, 2001, pp. 69-91.         [ Links ]

[61] Verma, S.P., Rivera-Gómez, M.A., Computer programs for the classification and nomenclature of igneous rocks: Episodes, 2013 (in press).         [ Links ]

[62] Verma, S.P., 2009, Continental rift setting for the central part of the Mexican Volcanic Belt: a statistical approach: Open Geology Journal, v. 3, p. 8-29.         [ Links ]

[63] Verma, S.P., Pandarinath, K., Verma, S.K., and Agrawal, S., 2013, Fifteen new discriminant-function-based multi-dimensional robust diagrams for acid rocks and their application to Precambrian rocks: Lithos, v. 168-169, p. 113-123.         [ Links ]

[64] Verma, S.P., and Verma, S.K., 2013, First fifteen probability-based multi-dimensional discrimination diagrams for intermediate magmas and their robustness against post-emplacement compositional changes and petrogenetic processes: Turkish Journal of Earth Sciences, in press.

[65] Verma, S.P., 2013. Application of 50 multi-dimensional discrimination diagrams and significance tests: deciphering compositional similarities and differences between Hawaiian and Icelandic volcanism International Geology Review, 2013. Vol. 00, No. 00, 1–20. http://dx.doi.org/10.1080/00206814.2013.788239.         [ Links ]