El Metate es un volcán en escudo situado en el sector sur del Campo Volcánico de Michoacán-Guanajuato, uno de los dos campos volcánicos más grandes del Cinturón Volcánico Transmexicano. Su actividad tuvo lugar aproximadamente 4.700 ± 200 a B.P y produjo más de quince coladas de afinidad calcoalcalina que muestran diferente grado de diferenciación.
Las temperaturas calculadas mediante geotermómetros mineral-líquido para olivino, plagioclasa, y piroxenos muestran que la fase más temprana en cristalizar fue el olivino (1232–1198 °C), seguido de plagioclasa (1162–1126 °C), ortopiroxeno (1147–1027 °C) y clinopiroxeno (1147–1018 °C). Las estimaciones de presión sugieren que la cristalización comenzó a ~7 kbar y continuó hasta niveles superficiales. El contenido en agua del fundido durante la cristalización de la plagioclasa fue de ~1.6% en peso.
Las temperaturas calculadas a partir del contenido en Al de los anfíboles, indican que cristalizaron entre 995 y 922 °C, a una presión media de 3.5 kbar y con un contenido en H2O del fundido de entre 5.2% y 6.9%. Aunque estos valores estarían de acuerdo con que el anfíbol representa una fase tardía en el proceso de cristalización, el que estos cristales presenten siempre texturas de desequilibrio, que se observan también de forma ocasional en otras fases, sugiere que pueden representar xenocristales y/o que la cristalización de estas lavas ha tenido lugar en sistema abierto.
ABSTRACT
El Metate is a shield volcano located in the southern sector of the Michoacan-Guanajuato Volcanic Field, one of two largest monogenetic volcanic fields of the Transmexican Volcanic Belt. It was active c. 4.700 ± 200 years B.P and emitted about fifteen calcalkaline lava flows showing variable differentiation degrees.
Temperatures calculated from mineral-liquid geothermobarometers for olivine, plagioclase and pyroxene, suggest that olivine was the earliest fractionating phase (1232–1198 °C), followed by plagioclase (1162–1126 °C), orthopyroxene (1147–1027 °C) and clinopyroxene (1147–1018 °C). Pressure estimations indicate that crystallization started at ~7 kbar and progressed up to surface levels. Water contents in the melts during crystalliztion of plagioclase is estimated at ~1.6%.
Temperatures calculated on the basis of Al content in amphibole, provide a crystallization range between 995 and 922 °C, at an average pressure of 3.5 kbar and water contents between 5.2% and 6.9%. Although these values could agree with a scenario where amphibole represents a late crystallization phase along the previous fractionating sequence, the systematic presence of disequilibrium textures, which are also observed occasionally in other phases, suggest that other possibilities such as open-system crystallization cannot be discarded.
Para establecer la evolución de los fundidos desde su acumulación en la cámara magmática hasta su ascenso y extrusión, es necesario conocer la temperatura y la presión a la que se formaron las diferentes fases que han cristalizado en el magma. Una de las aproximaciones más frecuentes consiste en estimar dichos parámetros a partir de la composición de las fases minerales, mediante la aplicación de geotermobarómetros basados en equilibrios químicos mineral-mineral y mineral-líquido. En la actualidad la aplicación de esta metodología se ha facilitado mucho al contarse con un completo arsenal de calibraciones para las fases más comunes en sistemas ígneos y adaptadas a un amplio rango de condiciones de presión, temperatura y composición (ver Putirka, 2008).
En efecto, estudios recientes sobre geotermobarometría en magmas basálticos (ver p.ej. Costa et al., 2013; Dahren et al., 2012; Keiding & Sigmarsson, 2012) han confirmado que este tipo de cálculos permite obtener información detallada sobre la evolución y funcionamiento del sistema magmático durante la diferenciación, incluidos aquellos que tienen lugar en sistema abierto, como la asimilación o la mezcla de magmas, y que son esenciales para perfilar cualquier modelo petrogenético. De hecho, como han demostrado otro tipo de estudios, contar con información sobre la evolución de la paragénesis mineral permite cuantificar con mayor precisión los procesos de diferenciación (ver p.ej. Cebriá et al., 2011).
Por estas razones, y en el contexto de una investigación más amplia sobre la petrogénesis del volcanismo en el Campo Volcánico de Michoacán-Guanajuato (México) se ha realizado un estudio geotermobarométrico (P-T-H2O) a partir de la paragénesis mineral observada en un primer conjunto de 11 muestras representativas de las lavas emitidas por el volcán El Metate, como paso previo para establecer un modelo petrogenético cuantitativo que pueda ser aplicable a volcanes similares de la región. Las muestras fueron seleccionadas de modo que cubrieran la secuencia de lavas emitidas por el volcán. En este trabajo se detallan los resultados preliminares de dicha aproximación.
Contexto geológico
El Cinturón Volcánico Transmexicano (CVTM), es considerado el mayor arco volcánico Neógeno de Norteamérica, extendiéndose por la parte central del país en una franja entre las latitudes 18° 30′ N y 21° 30′ N sobre el extremo sur de la placa norteamericana (Fig. 1). Está constituido por campos monogenéticos, estratovolcanes y mesetas volcánicas cuya distribución está controlada por la tectónica regional (Ferrari et al., 1999; Ferrari et al., 2012; Johnson & Harrison, 1989; Suter et al., 1999). La composición de las lavas muestra una gran variabilidad composicional, si bien predominan los términos calcoalcalinos (ver Ferrari et al., 2012; Gómez-Tuena et al., 2007 y referencias incluidas).
A) Situación del Cinturón Volcánico Mexicano (CVM) y del Campo Volcánico de Michoacán-Guanajuato (CVMG). B) Localización del volcán El Metate en el CVMG.
El Campo Volcánico de Michoacán-Guanajuato (CVMG) es uno de los dos campos monogenéticos más grandes del CVTM (Hasenaka, 1994; Hasenaka & Carmichael, 1985). Se sitúa en la parte centro-occidental del mismo y comprende cerca de 400 volcanes de tamaño medio (2–12 km de diámetro basal y 100–1000 m de altura; Hasenaka, 1994) y alrededor de 1000 afloramientos de menor entidad producidos por actividad monogenética, incluyendo maares, domos y conos de ceniza (Hasenaka & Carmichael, 1987), entre los cuales pueden encontrarse ~377 pequeños volcanes en escudo andesíticos (Hasenaka, 1994) y más raramente estratovolcanes (e.g. Tancítaro, Patamban), distribuidos en una superficie de unos 40000 km2. Se considera que el volcanismo monogenético en esta zona empezó hace ~3 Ma (Hasenaka & Carmichael, 1985), con una clara tendencia del frente volcánico a migrar hacia el suroeste. Predominan los términos calcoalcalinos, aunque también existen alcalinos y más raramente adakíticos (ver Ferrari et al., 2012; Gómez-Tuena et al., 2007 y referencias incluidas).
El Metate es un volcán en escudo que ocupa un área de algo más de 101 km2 en el sector sur del CVMG (Fig. 1). Produjo en torno a quince coladas de lavas de las que apenas existe información petrológica ni geoquímica (una única muestra de afinidad calcoalcalina descrita en Hasenaka & Carmichael, 1987). La edad de este volcán tampoco está bien definida, ya que solo existe una determinación radiométrica, basada en 14C, de 4700 ± 200 años (Hasenaka & Carmichael, 1985).
Métodos analíticos
Los análisis de elementos mayores en roca total se realizaron con un espectrómetro secuencial de rayos X Siemens SRS 3000 siguiendo el procedimiento descrito en Lozano-Santa Cruz & Bernal (2005) en el Laboratorio Universitario de Geoquímica Isotópica (LUGIS) del Instituto de Geología, Universidad Nacional Autónoma de México (UNAM).
La química mineral fue determinada mediante microsonda electrónica. Se utilizó una microsonda JEOL® “Superprobe JXA-8900 M” con cuatro analizadores WDS, del Centro Nacional de Microscopía Electrónica de la Universidad Complutense de Madrid. Los análisis se realizaron sobre láminas pulidas metalizadas con una película de grafito y las condiciones de medida fueron de 15 kV de diferencia de potencial y un haz de electrones de 20 nA de intensidad de corriente y 5 μm de diámetro. Los tiempos de medida fueron de 10 segundos en la posición del pico y 5 en la de cada fondo, con unos límites de detección de 0.5–6% en peso para óxidos con concentraciones >1.5% en peso y <10% para óxidos con concentraciones <1.5% en peso. Los patrones usados fueron: sillimanita, albita, almandino, kaersutita, microclina, ilmenita, fluorapatito, escapolita, aleación de Ni, cromita, gahnita, bentonita y estroncianita que proceden del Smithsonian Institute (Washington DC, EEUU) y de la Universidad de Harvard (Boston, EEUU). Se analizaron Si, Al, Fe, Mn, Mg, Ca, Na, K, Ti, Ni, Cr y P en olivino, plagioclasa, clinopiroxeno, ortopiroxeno y anfíbol. Excepto en el caso del anfíbol, los análisis se realizaron exclusivamente en aquellos cristales o zonas de los cristales que no mostraban evidencias de desequilibrio respecto a la matriz que los engloba.
Petrología y geoquímica
En la Tabla 1 se recoge la composición en elementos mayores de las rocas objeto de estudio y en las Tablas 2 a 6, se resume la composición de las fases consideradas en cada una de dichas muestras.
Análisis químicos de roca total (% en peso) de lavas del volcán El Metate
M20
M57
M60
M61
M63
M64
M65
M66
M70
M74
M77
SiO2
61.29
55.02
57.90
60.07
59.40
60.09
57.02
59.26
58.11
60.23
60.69
TiO2
0.67
0.92
0.80
0.69
0.70
0.68
0.90
0.68
0.73
0.69
0.60
Al2O3
17.69
17.48
18.27
18.00
17.36
17.42
17.29
18.08
17.53
17.83
18.60
Fe2O3t
5.38
7.63
6.11
5.66
6.21
5.89
7.30
5.71
6.67
5.81
5.45
MnO
0.09
0.13
0.11
0.09
0.11
0.10
0.13
0.10
0.10
0.09
0.08
MgO
3.09
5.78
3.36
3.19
4.01
3.40
5.05
3.17
4.65
3.39
3.13
CaO
6.59
7.56
7.10
6.58
6.43
6.11
6.84
6.75
6.94
6.06
5.94
Na2O
3.91
3.44
3.94
3.98
3.68
3.78
3.63
4.03
3.73
3.93
3.99
K2O
1.45
1.18
1.61
1.48
1.80
1.76
1.49
1.41
1.52
1.89
1.43
P2O5
0.21
0.05
0.11
0.00
0.02
0.01
0.05
0.04
0.21
0.27
0.17
LOI
0.34
0.50
0.26
0.00
0.24
0.25
0.00
0.30
0.22
0.25
0.52
Total
100.70
99.68
99.56
99.73
99.95
99.47
99.70
99.53
100.42
100.45
100.61
*LOI: pérdida por calcinación.
Análisis químicos, porcentaje de componente forsterita y temperatura calculada de fenocristales y microfenocristales de olivinos en lavas del volcán El Metate
M57
M65
1
2
3
4
5
6
1
Feno.
Microfeno.
Feno.
Microfeno.
Feno.
Feno.
Feno.
SiO2
38.96
38.92
39.07
39.94
36.84
37.99
37.42
Al2O3
—
—
—
0.01
0.22
0.04
—
FeOT
19.68
21.14
20.79
15.99
23.45
22.68
25.24
MnO
0.18
0.24
0.26
0.22
0.42
0.26
0.30
MgO
41.28
40.78
40.15
44.89
36.29
38.16
36.24
CaO
0.16
0.23
0.18
0.14
0.21
0.14
0.19
Na2O
0.01
—
—
0.01
0.01
0.03
0.03
K2O
—
0.00
0.01
0.01
0.01
0.01
0.02
TiO2
0.03
0.04
0.03
0.00
0.04
—
0.01
NiO
0.08
0.07
0.20
0.27
0.07
0.09
0.12
Cr2O3
0.02
—
—
—
—
—
0.02
P2O5
0.05
0.09
0.05
0.04
—
0.03
0.04
F
—
0.01
—
—
—
—
—
Cl
—
0.01
—
0.00
0.03
—
—
Total
100.44
101.52
100.73
101.51
97.55
99.42
99.62
% Fo
79
77
77
83
73
75
72
T(°C)
1212
1216
1218
1198
1232
1225
1225
Análisis químicos expresados en % en peso del óxido.
Todos los análisis en fenocristales han sido realizados en los bordes.
Todas las rocas consideradas presentan textura hipocristalina, inequigranular, porfídica y vacuolar (Fig. 2). La matriz es microcristalina y está compuesta por plagioclasa (Pl), ortopiroxeno (Opx), clinopiroxeno (Cpx), óxidos, olivino (Ol), y ocasionalmente vidrio. A excepción de los óxidos, las fases que componen la matriz aparecen asimismo como feno y microfenocristales, siendo Pl y piroxenos comunes a todas las muestras, mientras que el Ol aparece solo ocasionalmente. El anfíbol (Amp) también aparece como fenocristal en algunas de las muestras estudiadas. Por lo que respecta a su composición química (Tabla 1), las lavas de El Metate presentan un rango de variación en SiO2 relativamente reducido por lo que están representadas casi exclusivamente por términos andesíticos (Fig. 3), con bajas relaciones FeOT/MgO. Los relativamente altos contenidos de SiO2 (61–55%) y bajos de MgO (5.8−3%) y del índice #Mg (66–58) indican que en todos los casos se trata de rocas evolucionadas.
A) Fenocristales de olivino y clinopiroxeno en una matriz de plagioclasas de Tipo 4 (muestra M-57, nícoles cruzados). B–C) Ejemplos de plagioclasas de los tipos 1, 2 y 3 e inclusiones de apatito (muestra M-64, nícoles cruzados); D) Anfíbol con corona de opacita granular compuesta por Plg+Cpx+óxidos (muestra M-20, nícoles cruzados); E) pseudomorfo de anfíbol completamente reemplazado (muestra M-60, nícoles cruzados); F) Anfíboles con corona de opacita simplectítica (muestra M-57, nícoles paralelos).
Proyección de las lavas de El Metate (círculos rojos y verdes) y rocas del Campo Volcánico de Michoacán–Guanajuato (cruces, según datos en Gómez-Tuena et al., 2007 y referencias incluidas) en el diagrama TAS (Le Bas et al., 1992). Datos en porcentaje en peso.
En general las concentraciones de MgO, Fe2O3, y CaO, disminuyen respecto a SiO2 (Fig. 4), mientras que K2O, Na2O, P2O5 y en cierta medida Al2O3, tienden a aumentar. Sin embargo, la distribución de estas rocas en los diagramas SiO2-elemento revela la existencia de dos tendencias contrastadas. Como se puede observar sobre todo en el diagrama SiO2-MgO, un grupo de muestras (M-20, M-60, M-61 y M-66), se distribuyen a lo largo de una pauta de menor pendiente respecto a la trayectoria representada por el resto. Asimismo las pautas generales de enriquecimiento/empobrecimiento no son siempre coincidentes. Así el primer grupo de muestras exhibe un empobrecimiento en K2O, Na2O y Al2O3 al aumentar el contenido de SiO2, mientras que el resto muestra un enriquecimiento.
Diagramas SiO2-óxidos (en % en peso) para las lavas de El Metate. Se indican las dos pautas de diferenciación observadas.
Fases mineralesOlivino
Los olivinos son en su mayoría feno y microfenocristales (Fig. 2A), aunque también aparecen como microcristales formando la matriz. En general son subidiomorfos, de bordes rotos o con golfos de corrosión, así como con zonas internas huecas. Si bien al microscopio es difícil de distinguir, en la microsonda se ha podido detectar que los fenocristales pueden presentar bordes de recrecimiento y en algún caso, una fina corona de Opx. El olivino aparece exclusivamente en las muestras con mayor contenido en MgO (>4.5%) y muestran un rango composicional relativamente estrecho entre Fo83 y Fo72 (Tabla 2).
Plagioclasa
La plagioclasa se encuentra presente en todas las rocas estudiadas como fenocristales, microfenocristales, formando parte de la matriz y como inclusiones microcristalinas en fenocristales de otras fases. Es posible distinguir cuatro poblaciones principales de cristales de Pl (Fig. 2B, C y D):
Tipo 1:
Cristales idio-subidiomorfos que pueden ser tanto feno como microfenocristales o cristales de la matriz y que presentan macla de Carlsbad o polisintética. Esporádicamente algunos fenocristales pueden tener extinción ondulante o presentar un núcleo alterado.
Tipo 2:
Feno y microfenocristales con zonación oscilatoria muy marcada, en ocasiones con extinción ondulante sobreimpuesta, que pueden tener bordes corroídos y recrecimientos y/o núcleos alterados y con inclusiones. Con frecuencia en estas plagioclasas se encuentra gran cantidad de microcristales de apatito incluidos en los bordes recrecidos (que también se observan en la matriz en algunas muestras), siendo las inclusiones del núcleo indistinguibles.
Tipo 3:
Fenocristales totalmente alterados y de bordes corroídos, en los que se sigue apreciando el maclado, y con un borde de recrecimiento también corroído (textura de reabsorción generalizada tipo “spongy cellular, sieved” (ver Streck, 2008).
Tipo 4:
Microcristales que se encuentran en coronas que rodean a los cristales de anfíbol o como inclusiones en algunos piroxenos (Fig. 2D).
El rango composicional observado en las plagioclasas analizadas es relativamente amplio, entre An71 y An54, y corresponden en su mayor parte a labradoritas (Fig. 5), excepto dos ejemplares de bitownita identificados en la matriz microcristalina. A pesar de la existencia de cuatro poblaciones principales de plagioclasas, el rango composicional observado en los análisis de feno- y microfenocristales y microcristales (ver la sección de geotermometría) muestran un rango equivalente al del conjunto de plagioclasas, con cierta tendencia de los cristales de la matriz a proyectarse en el extremo más anortítico del espectro (labradorita-bitownita), mientras que los bordes de microfenocristales se proyectan en la región de tendencia albítica.
Diagramas de clasificación de plagioclasas (Ab-An-Or) y de piroxenos (Wo-En-Fe).
Las plagioclasas de Tipo 3 (de textura de reabsorción generalizada tipo “spongy cellular sieved”) no han podido ser caracterizadas geoquímicamente dado el alto grado de alteración que presentan. Tampoco se han caracterizado las plagioclasas presentes en las coronas de los anfíboles, ya que pueden representar productos de reacciones de desequilibrio.
Ortopiroxeno
El ortopiroxeno raramente se encuentra como fenocristal, pero es frecuente como microfenocristal, asi como en la matriz o formando parte del núcleo de un cristal con recrecimiento de clinopiroxeno.
Los cristales analizados se clasifican como enstatitas (Fig. 5) y tienen un rango composicional que oscila entre valores de #Mg (#Mg = Mg/(Mg+Fe+2) × 100) desde 81% a 73% (Tabla 4) para los microfenocristales, mientras que los cristales de la matriz se proyectan en un rango composicional algo más estrecho (79–74%).
Clinopiroxeno
El clinopiroxeno aparece como fenocristales, microfenocristales así como en la matriz. La mayor parte de los microfeno y microcristales tienen inclusiones de opacos, mientras que los fenocristales muestran zonas isótropas, principalmente en el núcleo, y se caracterizan por tener un estrecho borde de recrecimiento (Fig. 2A). Algunos cristales exhiben además características particulares: pueden encontrarse en continuidad óptica con Opx formando un zonado en parches, o bien hallarse en forma de agregados de cristales con zonación concéntrica, donde el núcleo está formado por Opx y las zonas externas por Cpx.
El valor #Mg de los Cpx analizados varía entre 97–80% (Tabla 5) y se proyectan como augitas en el diagrama de clasificación Mg-Ca-Fe (Fig. 5). Los cristales de la matriz y los microfenocristales presentan las composiciones más magnésicas, con rangos de #Mg entre 84-92%. En cambio, los bordes de los fenocristales analizados abarcan un amplio abanico de composiciones (#Mg = 80–97%).
Anfíbol
El Amp aparece como fenocristal y también como microfenocristal y muestra una compleja historia de alteración que incluye dos tipos de opacitización (sensu Plechov et al., 2008). Por un lado, se puede observar opacitización en coronas, que a su vez están compuestas por dos zonas: una parte interna que se desarrolla directamente en contacto con el anfíbol (opacita granular) y compuesta por Pl, Cpx y óxidos de entre 3 y 10 μm; y una parte externa de opacita simplectítica, que en ocasiones es la única que se desarrolla, compuesta por un agregado isótropo de minerales de grano fino cuya composición no pudo ser determinada. Por otro lado, también se observa opacitización en forma de áreas de descomposición volumétrica. Se denomina de este modo a zonas de reemplazamiento en el interior del cristal que han sido alteradas a través de planos de debilidad y que en este caso están compuestas mayoritariamente por cristales de clinopiroxeno que heredan la orientación estructural del anfibol al que reemplazan.
Atendiendo al tipo de corona de opacita que presentan, se pueden distinguir dos tipos de cristales: a) los casi euhedrales con tamaños que pueden ir desde fenocristales hasta microfenocristales, que no suelen presentar alteración en áreas de descomposición volumétrica y muestran únicamente coronas de opacita simplectítica (isótropa) (Fig. 2F), y b) los que presentan corona de opacita granular de Pl, Cpx y opacos muy desarrollada (Fig. 2D), que en ocasiones llega a reemplazar todo el mineral (pseudomorfos) (Fig. 2E). Estos cristales no suelen presentar opacita isótropa y son de tamaño inferior a los primeros. Todas las muestras con anfíbol incluyen ambos tipos, si bien los primeros son los más frecuentes.
En todos los casos se corresponden a anfíboles cálcicos y se clasifican como Mg-hastingsitas, a excepción de una pargasita (mayor contenido en Na respecto a Ca). No se han observado diferencias composicionales entre los diferentes tipos de cristales identificados petrográficamente.
Geotermobarometría
A partir de las diferentes formulaciones de geotermómetros y geobarómetros disponibles en la literatura (ver Putirka, 2008) se seleccionaron aquellos que mejor se adaptan en cada caso al rango composicional de las rocas objeto de estudio, así como a la paragénesis representada por Pl, Cpx, Opx, Ol y Amp.
La mayoría de estos geotermobarómetros estiman las condiciones P-T a partir de equilibrios químicos mineral-mineral o mineral-líquido. En consecuencia, para aplicarlos de forma estricta, solo se pueden considerar aquellos cristales o zonas de los mismos, que no presenten evidencias de desestabilización y se encuentren en contacto con la fase (sea mineral o líquido) con la que se encuentra en equilibrio. En consecuencia, los cálculos basados en equilibrios mineral-líquido se han limitado a las composiciones obtenidas en aquellas zonas de los minerales sin texturas de desequilibrio y que se encuentran en contacto con la matriz, asumiendo que la composición de la roca total representa de forma aproximada la composición del magma en equilibrio con las fases que contiene. Este es el caso de los termobarómetros para Ol, Cpx, Opx y Pl, que se basan en relaciones de intercambio entre el mineral y el fundido en equilibrio. Sin embargo, para el Amp las consideraciones anteriores no son aplicables debido a que el termobarómetro disponible utiliza únicamente la composición del mineral, por lo que se pueden utilizar tanto las composiciones obtenidas en el centro como en el borde del cristal.
Olivino
En el caso del olivino hemos adoptado la ecuación 22 de Putirka (2008), basada en el equilibrio olivino-líquido:T(°C)={15294.6+1318.8P(GPa)+2.4834[P(GPa)]2}/{8.048+2.8352lnDMgOl/liq+2.097 ln[1.5(CNML)]+2.575ln[3(CSiO2liq)]-1.41NM+0.222H2Oliq+0.5P(GPa)}
Esta ecuación es esencialmente una modificación del modelo de Beattie (1993), basado en el reparto de Mg entre olivino y líquido, pero que evita las sobreestimaciones en el cálculo de T en condiciones hidratadas.
Para testar la existencia de equilibrio entre olivino y el líquido, se utiliza el coeficiente de distribución Fe-Mg de la reacción:
Como demostraron Roeder & Emslie (1970), este parámetro varía muy poco con la T o la composición del par olivino-líquido y tiene un valor prácticamente constante de 0.30 ± 0.03 en sistemas basálticos a P<2–3 GPa.
Plagioclasa
En el caso de la plagioclasa, teniendo en cuenta la habitual presencia de agua en rocas calcoalcalinas (típicamente entre 1% y 4% en volcanes del CVMG; Johnson et al. 2010), hemos adoptado la ecuación 24a de Putirka (2008), que considera la posible presencia de agua en el líquido y reduce los errores respecto a calibraciones anteriores (p.ej. Ghiorso et al., 2002; Mathez, 1973; Putirka, 2005).104T(K)=6.4706+0.3128ln{XAnPlXCaOliq(XAlO1.5liq)2(XSiO2liq)2}-8.103(XSiO2liq)+4.87(XKO0.5liq)+1.5346(XAbPl)2+8.661(XSiO2liq)-3.341*10-2(P(kbar))+0.18047(H2Oliq)
Para los cálculos de presión hemos utilizado la ecuación. 25a de Putirka (2008), originaria de Putirka (2005), por coherencia con el termómetro.P(kbar)=-42.2+4.94×10-2T(K)+1.16×T(K)ln{XAbPlXAlO1.5liqXCaOliqXAnPlXNaO0.5liqXSiO2liq}-382.3(XSiO2liq)2+514.2(XSiO2liq)3-19.6ln(XAbPl)-139.8(XCaOliq)+287.2(XNaO0.5liq)+163.9(XKO0.5liq)
Como prueba de equilibrio del par Pl-liq se ha utilizado la constante de equilibrio KD(An-Ab)Pl-liq. Aunque este parámetro varía respecto a P, T y H2O, al ser dividido en dos intervalos de T se obtienen valores relativamente constantes de KD(An-Ab)Pl-liq = 0,10 ± 0,05 para T<1050 °C y KD(An-Ab)Pl-liq = 0,27 ± 0,11 para T >1050°C.
También se ha empleado la expresión 25b de (Putirka, 2008) ya que el equilibrio Pl-líquido proporciona un buen higrómetro cuando T es conocida:H2O(Wt)=28.95-0.0032T(°C)ln{XAnPlXCaOliq(XAlO1.5liq)2(XSiO2liq)2}-18.9(XKO0.5liq)+14.5(XMgOliq)-40.3(XCaOliq)+5.7(XAnpl)+0.108P(Kbar)
Puesto que estas expresiones requieren conocer previamente algunos parámetros de las otras, la resolución se ha efectuado mediante cálculo iterativo.
Ortopiroxeno
Para estimar las temperaturas de cristalización del ortopiroxeno hemos adoptado el termómetro calibrado por Putirka (2008) en su ecuación 28a. Esta expresión rectifica algunos problemas de sobreestimación de T asociados con el anterior termómetro de Beattie (1993), que hasta entonces era el único disponible, basado en el equilibrio Opx-líquido. La nueva calibración es aplicable en un amplio rango de P y entre 750–1600 °C, SiO2 = 33–77% y H2O = 0–14.2%, con un error de ± 39 °C (n = 793).104T(°C)=4.07-0.329[P(GPa)]+0.12[H2Oliq]+0.567ln[XFm2SiO2O6Opx(XSiO2liq)2(XFeOliq+XMnOliq+XMgOliq)2]-3.06[XMgOliq]-6.17[XKO0.5liq]+(Mg#liq)+2.57[XFeOpx]
donde Fm = Fe + Mn + Mg
Para los cálculos de P hemos seleccionado la ecuación 29a de Putirka (2008), que utiliza la correlación positiva de FmAl2SiO6 respecto a P (entre 0.001 y 0.05 kbar) y proporciona estimaciones con un error de ± 2.6 kbar.P(kbar)=-13.97+0.0129T(°C)+0.001416T(°C)ln[XNaAlSiO2O6OpxXNaO0.5liqXAlO1.5liq(XSiO2liq)2]-19.64(XSiO2liq)+47.49(XMgOliq)+6.99(XFeOpx)+37.37(XFmAl2SiO6Opx)+0.748(H2Oliq)+79.67(XNaO0.5liq+XKO0.5liq)
Como test de equilibrio del par Opx-líquido nos basamos en el diagrama de Rhodes (Rhodes et al., 1979) y el valor KD (Fe-Mg)Opx-liq = 0.29 ± 0.06 (para n = 785). Este valor no depende de P y T, pero disminuye ligeramente con el aumento en contenido en sílice de modo que KD (Fe-Mg)Opx-liq = 0,4805−0,3733*XSiliq.
Clinopiroxeno
Los barómetros para clinopiroxenos se basan en la reacción NaO0.5liq + AlO1.5liq + 2SiO2liq = NaAlSi2O6cpx, si bien los modelos más recientes han sido calibrados usando H2Oliq (en % en peso) como variable para describir mejor las muestras hidratadas (Grove & Juster, 1989; Kinzler & Grove, 1992; Patiño Douce, 2005; Putirka et al., 1996; Putirka, 2008; Scaillet & MacDonald, 2003; Sisson & Grove, 1993a; Sisson & Grove, 1993b; Walter & Presnall, 1994). De estas, hemos seleccionado para el cálculo de P la expresión 32c de Putirka (2008), por su buen ajuste respecto a valores experimentales (error estimado de ± 1.5 kbar para n = 99).P(kbar)=-57.9+0.0475T(K)-40.6(XFeoliq)-47.7(XCaTsCpx)+0.676(H2Oliq)-153(XCaO0.5liqXSiO2liq)+6.89[XAlCpxXAlO1.5lq]
Análogamente, para los cálculos de temperatura se ha utilizado la ecuación 33 de (Putirka, 2008), con errores estimados de ± 42 °C (para n = 320 hidratado) y ± 46 °C (para n = 854 anhidro).104T(K)=7.53-0.14ln(XldCpxXCaOliqXFmliqXDiHdCpxXNaliqXAlliq)+0.07(H2Oliq)-14.9(XCaOliqXSiO2liq)-0.08ln(XTiO2liq)-3.62(XNaO0.5liq+XKO0.5liq)-1.1(Mg#liq)-0.18ln(XEnFsCpx)-0.027P(Kbar)
Como test de equilibrio se ha calculado el valor de KD(Fe-Mg)Cpx-liq, que adopta un valor de 0,28 ± 0,08 mostrando poca variación respecto a la temperatura (Putirka, 2008).
Anfíbol
Para calcular los valores de P y T de cristalización de los anfíboles hemos adoptado el reciente termobarómetro de Ridolfi et al. (2010), que no se basa en el equilibrio con el líquido y representa el único modelo actualmente disponible para anfíboles en magmas calcoalcalinos. Esta expresión permite asimismo calcular el porcentaje de H2O del magma a partir del cual cristalizaría el anfíbol analizado.T=-151.487Si*+2,041
Como se ha indicado anteriormente, la aplicación de geotermómetros basados en equilibrios mineral-líquido plantea ciertas limitaciones ya que requiere contar con pares que no muestren signos de desequilibrio (p.ej. bordes corroídos o coronas de reacción). Además, si se acepta que la composición del líquido está representada por la roca total, es necesario asumir que el equilibrio se verifica únicamente con el borde de los minerales presentes en la roca. En consecuencia, las temperaturas obtenidas se limitan a la etapa final de cristalización de dichas fases. Para obtener las temperaturas a las que se produjo la cristalización de zonas internas de los minerales, sería necesario conocer la composición del líquido en equilibrio en el momento de la cristalización de dichas zonas.
De acuerdo con lo expuesto, en este trabajo los cálculos basados en geotermómetros mineral-líquido se han restringido a aquellos análisis obtenidos en minerales que no muestran evidencias texturales de desequilibrio respecto a la matriz y limitados bien a minerales que no muestran zonación (usualmente microcristales), o bien a bordes de fenocristales o microfenocristales. Asimismo, como se ha indicado más arriba, se han considerado únicamente aquellos pares mineral-líquido que proporcionan valores KD dentro de los rangos considerados de equilibrio en cada caso. En las Tablas 2 a 6 y en la Fig. 6 se recogen los resultados obtenidos de esta aproximación.
A) Diagrama de temperatura vs. SiO2 (% en peso) en roca total; y B) Diagrama de presión vs temperatura obtenidos para olivino, plagioclasa, ortopiroxeno, clinopiroxeno y anfíbol en las lavas de El Metate.
En el caso del Ol (Tabla 2), las temperaturas obtenidas oscilan entre 1232−1198°C, con una tendencia general a disminuir la temperatura de cristalización a medida que progresa la diferenciación (Fig. 6A).
Los cálculos para equilibrios Pl-liq proporcionan unas temperaturas y presiones de 1162–1126 °C, y 7.3−3.6 kbar respectivamente (Tabla 3). No se han observado diferencias significativas en los valores de T respecto al contenido en SiO2 entre los distintos tipos de plagioclasas ni entre fenocristales y microcristales. Por otro lado, la variación de P respecto a T sugiere un rango de presiones bastante más amplio que el que se observa en las temperaturas calculadas.
Análisis químicos, porcentaje de los componentes An y Ab, y valores de temperatura (T), presión (P) y contenido en agua (H2O%) calculados para cristales de plagioclasa en lavas del volcán El Metate
M20
M57
M60
M61
M63
M64
M65
M66
M77
1
1
1
1
1
1
2
1
1
2
3
1
Feno.
Microfeno.
Feno.
Feno.
Matriz
Feno.
Matriz
Matriz
Feno.
Feno.
Microfeno.
Matriz.
SiO2
53.03
53.58
52.51
52.90
49.91
54.71
52.18
50.54
51.87
52.87
52.45
54.14
Al2O3
28.22
29.05
29.01
28.35
30.59
27.84
28.42
30.72
28.96
28.02
28.62
28.98
FeOT
0.68
0.76
0.66
0.93
0.71
0.86
0.93
0.61
0.75
0.85
0.57
0.69
MnO
—
—
0.02
—
0.03
—
0.05
0.03
0.04
0.09
—
—
MgO
0.10
0.17
0.06
0.05
0.06
0.05
—
0.09
0.12
0.08
0.08
0.08
CaO
12.17
12.21
12.60
11.42
14.90
11.19
12.78
14.54
12.84
12.34
12.87
12.02
Na2O
4.49
3.88
4.70
4.93
3.43
5.17
4.53
3.31
4.37
4.81
4.65
4.56
K2O
0.15
0.17
0.18
0.07
0.15
0.36
0.30
0.15
0.13
0.11
0.11
0.13
TiO2
0.01
0.01
0.01
0.06
0.02
0.01
0.07
0.01
0.01
0.01
—
0.07
NiO
0.08
—
0.05
0.02
0.03
—
—
0.02
0.02
—
0.03
—
Cr2O3
0.01
—
0.00
0.00
—
0.03
—
—
0.00
—
—
—
P2O5
0.08
0.02
0.02
0.07
0.00
—
0.04
—
0.00
0.03
0.07
0.01
F
0.03
—
0.05
0.01
0.06
0.06
—
0.03
0.00
—
—
0.02
Cl
0.01
0.00
0.00
0.03
0.00
0.02
0.01
0.00
0.03
0.00
—
0.00
Total
99.05
99.84
99.87
98.83
99.87
100.29
99.30
100.06
99.12
99.21
99.46
100.70
% An
60.0
63.5
59.7
56.1
70.6
54.5
60.9
70.8
61.9
58.6
60.5
59.3
% Ab
40.0
36.5
40.3
43.9
29.4
45.5
39.1
29.2
38.1
41.4
39.5
40.7
T (°C)
1135
1148
1150
1148
1144
1126
1133
1162
1147
1143
1145
1128
H2O%
1.7
1.5
1.5
1.5
1.5
1.8
1.6
1.2
1.6
1.7
1.6
2.1
P (Kbar)
5.4
4.7
7.3
7.2
3.7
6.5
5.3
3.6
5.9
6.5
6.1
5.2
Análisis químicos expresados en % en peso del óxido. Todos los análisis en fenocristales han sido realizados en los bordes.
T calculada según ecuación 3; H2O calculada según ecuación 4; P calculada según ecuación 5.
Por su parte, la cristalización de los ortopiroxenos abarca un rango de temperaturas entre 1147 °C y 1027 °C (Tabla 4), con una tendencia general a la T al aumentar el contenido de SiO2 (Fig. 6A). En el gráfico P-T los resultados se disponen con pocas excepciones en una relación directa de disminución de P con respecto a T (Fig. 6B). Como en el caso de las plagioclasas, no se observan diferencias entre las dos trayectorias de evolución o los tipos de Opx identificados.
Análisis químicos, número de magnesio y valores de temperatura (T) y presión (P) calculados para cristales de ortopiroxeno en lavas del volcán El Metate
M60
M64
M65
M66
M70
M74
1
2
1
2
1
1
2
3
1
1
Microfeno.
Microfeno.
Feno.
Microfeno.
Matriz
Matriz
Microfeno.
Microfeno.
Microfeno.
Microfeno.
SiO2
55.30
55.11
52.95
52.65
54.22
53.69
53.13
54.49
53.91
55.74
Al2O3
1.22
0.93
0.80
1.96
0.24
1.58
2.19
0.92
0.46
0.83
FeOT
12.33
12.59
16.92
15.50
14.89
13.39
13.73
12.70
15.00
12.95
MnO
0.36
0.38
0.45
0.35
0.28
0.38
0.23
0.28
0.32
0.23
MgO
28.40
28.93
26.44
27.10
27.06
28.40
27.98
28.61
27.38
29.36
CaO
1.42
1.55
1.68
1.76
2.18
2.14
1.62
1.60
2.05
1.59
Na2O
0.05
0.05
0.02
0.03
0.05
0.05
0.02
0.04
0.04
0.04
K2O
0.01
0.00
0.00
0.01
0.01
0.03
—
—
0.00
0.00
TiO2
0.28
0.30
0.21
0.35
0.19
0.32
0.31
0.16
0.21
0.07
NiO
0.02
0.01
0.03
0.00
0.00
0.02
0.01
0.06
0.00
0.00
Cr2O3
0.00
0.00
0.00
0.03
0.01
0.01
0.02
—
0.00
0.02
P2O5
0.01
0.07
0.00
0.00
0.00
0.06
0.02
0.01
0.03
0.00
F
0.04
0.00
0.00
0.02
0.02
0.00
—
—
0.00
0.00
Cl
0.01
0.02
0.00
0.02
0.00
0.00
0.02
—
0.00
0.01
Total
99.42
99.92
99.50
99.77
99.15
100.08
99.24
98.86
99.39
100.85
#Mg
80
81
78
80
79
74
81
81
79
73
T (°C)
1078
1147
1027
1042
1092
1053
1104
1113
1071
1066
P (Kbar)
2.1
2.2
0.0
0.2
2.8
0.6
0.1
1.1
2.0
0.7
Análisis químicos expresados en % en peso del óxido. Todos los análisis en fenocristales han sido realizados en los bordes
T calculada según ecuación 6; P calculada según ecuación 7.
El clinopiroxeno presenta el rango más amplio de temperaturas, desde los 1147 °C a los 1018 °C (Tabla 5), con cierta tendencia a disminuir el valor de T con el aumento de SiO2 (Fig. 6A). La variación de P con respecto a la T sigue, como en el resto de fases analizadas, la relación directa habitual. Al igual que en el caso de Pl y Opx, tampoco se observan diferencias significativas en función de las diferentes tipologías de cristales analizados o de las pautas de evolución identificadas a partir de elementos mayores en roca total.
Análisis químicos, número de magnesio y valores de temperatura (T) y presión (P) calculados para cristales de clinopiroxeno en lavas del volcán El Metate
M60
M61
M63
M64
M66
M70
M74
M77
1
2
1
1
1
2
3
1
1
2
1
1
Matriz
Microfeno.
Feno.
Feno.
Feno.
Microfeno.
Microfeno.
Feno.
Feno.
Feno.
Feno.
Feno.
SiO2
51.36
49.18
51.12
49.80
50.48
50.313
51.941
50.64
50.20
50.79
49.91
50.90
Al2O3
2.85
3.89
3.32
3.50
3.01
2.706
2.302
2.81
3.66
2.98
3.45
2.69
FeOT
7.61
8.15
6.36
8.23
6.70
9.051
6.587
9.31
8.01
8.89
8.05
7.91
MnO
0.34
0.28
0.20
0.20
0.08
0.14
0.124
0.26
0.15
0.12
0.20
0.22
MgO
16.44
15.13
16.42
15.14
16.52
15.465
16.599
18.42
14.75
16.14
15.55
15.51
CaO
20.84
21.77
22.62
21.13
21.78
20.594
22.026
17.39
21.09
18.65
21.47
21.71
Na2O
0.36
0.45
0.36
0.35
0.31
0.337
0.261
0.31
0.36
0.33
0.40
0.27
K2O
0.06
—
0.01
0.03
0.00
—
0.012
0.00
0.00
0.03
—
0.00
TiO2
0.55
1.00
0.41
0.74
0.50
0.674
0.342
0.72
0.74
0.79
0.50
0.45
NiO
—
0.05
0.05
0.10
0.02
0.034
0.034
—
0.04
—
0.02
0.03
Cr2O3
—
0.07
0.06
0.07
0.29
—
0.094
0.03
0.07
—
0.03
—
P2O5
0.00
0.29
0.03
0.01
0.04
—
0.016
—
0.08
0.01
0.05
0.03
F
—
0.05
0.05
—
0.07
—
0.023
0.00
0.02
—
0.01
0.06
Cl
0.02
0.01
—
0.03
0.00
—
—
0.02
0.01
—
0.00
—
Total
100.43
100.31
101.00
99.33
99.79
99.31
100.36
99.90
99.17
98.71
99.62
99.76
#Mg
89
92
97
87
96
84
91
89
83
80
91
88
T (°C)
1111
1098
1069
1097
1045
1080
1069
1018
1128
1147
1099
1064
P (Kbar)
4.9
5.9
3.6
5.5
2.5
4.1
2.9
0.4
6.3
6.5
6.3
3.6
Análisis químicos expresados en % en peso del óxido. Todos los análisis en fenocristales han sido realizados en los bordes.
T calculada según ecuación 9; P calculada según ecuación 8.
Los anfíboles analizados indican un rango de variación de T para la cristalización de esta fase entre 995 y 922 °C (Tabla 6), disponiéndose en una pauta de pendiente negativa respecto al aumento del contenido en SiO2 de la roca correspondiente (Fig. 6A). Tanto los núcleos como los bordes analizados se proyectan en todo el rango de variación de T-SiO2, aunque en este caso se aprecia cierta disminución de T de centro a borde en los cristales zonados. Por otro lado, las muestras con menor SiO2 presentan los valores de P más altos, siguiendo una tendencia P-T análoga al resto de fases minerales (Fig. 6B).
Análisis químicos, clasificación y valores de temperatura (T), presión (P) y contenido en agua del fundico (H2O) calculados para cristales de anfíbol en lavas del volcán El Metate
M20
M60
M61
M66
M74
M77
1
2
1
1
1
2
1
1
Medio
Centro
a Centro
b Medio
c Borde
Centro
Centro
Centro
a Centro
b Borde
Centro
SiO2
45.04
45.13
42.96
41.72
42.57
43.69
42.58
43.68
43.54
44.23
44.839
Al2O3
11.38
11.19
12.51
12.90
12.66
12.72
11.97
10.74
12.19
11.25
11.001
FeOT
8.69
8.16
8.37
10.73
9.64
9.41
10.14
10.33
11.11
8.53
8.349
MnO
0.08
0.11
0.04
0.09
0.09
0.12
0.14
0.14
0.12
0.08
0.161
MgO
16.45
17.57
16.74
14.99
16.37
15.47
15.88
15.64
15.02
16.57
16.697
CaO
12.17
11.99
12.28
11.54
11.74
11.68
11.95
11.79
11.82
11.85
11.159
Na2O
2.18
2.09
2.40
2.36
2.55
2.35
2.38
2.39
2.25
2.23
2.152
K2O
0.38
0.37
0.54
0.51
0.51
0.49
0.39
0.37
0.52
0.61
0.382
TiO2
1.94
1.55
2.08
2.35
1.97
1.50
2.13
2.30
1.38
1.38
1.204
NiO
—
0.10
0.05
—
0.13
—
0.12
—
0.04
0.03
0.008
Cr2O3
0.07
0.34
0.29
0.01
0.02
—
0.01
0.03
0.04
0.13
0.095
P2O5
—
—
0.02
—
0.02
0.03
0.03
0.02
0.02
0.03
0.044
F
0.11
0.02
0.12
0.09
0.15
0.06
0.11
0.02
0.06
0.08
0.051
Cl
0.04
0.04
0.03
0.04
0.04
0.04
0.10
0.02
0.04
0.03
0.01
Total
98.52
98.64
98.41
97.32
98.44
97.54
97.93
97.47
98.14
97.03
96.15
Clasificación
magnesio-Hastingsita
magnesio-Hastingsita
magnesio-Hastingsita
magnesio-Hastingsita
magnesio-Hastingsita
Pargasita
magnesio-Hastingsita
magnesio-Hastingsita
magnesio-Hastingsita
magnesio-Hastingsita
magnesio-Hastingsita
T (°C)
944
943
995
990
990
964
978
922
954
946
943
P (Kbar)
2.9
2.7
3.9
4.5
4.0
4.2
3.5
2.8
3.7
3.0
2.7
H2O (%)
6.0
5.5
5.5
6.3
5.5
6.9
5.7
5.2
6.5
5.3
5.9
Análisis químicos expresados en % en peso del óxido.
T calculada según ecuación 10; P calculada según ecuación 11; H2O calculada según ecuación 12.
Discusión y conclusiones
Como se desprende de las variaciones composicionales basadas en elementos mayores, puestas de manifiesto en los diagramas SiO2-elemento (Fig. 4), las lavas de El Metate se han fraccionado a partir de dos pautas de diferenciación constrastadas. Sin embargo, por lo que respecta a sus características petrográficas, apenas se observan diferencias significativas entre ambas tendencias. Así, la presencia de olivino, que está limitada a una de las pautas de evolución, coincide con los términos de menor grado de diferenciación (es decir, con mayor MgO y #Mg, y menor SiO2), que están ausentes en la otra línea de diferenciación. Asimismo, si bien la presencia de anfíbol es característica de la serie con menor abundancia de MgO, dicha fase también se encuentra en las rocas con mayor grado de diferenciación del otro grupo (y por ende, con menor contenido relativo en MgO y K2O). A pesar de estas diferencias, en todos los casos Pl, Opx y Cpx, aparecen como fenocristales y por tanto deben haber participado en el proceso de diferenciación de ambas series, sin que tampoco se hayan detectado diferencias significativas en lo que se refiere a la composición o la textura de ninguna de las fases estudiadas entre ambos grupos de rocas.
Por lo que respecta a las características del proceso de diferenciación, la frecuente presencia de texturas de desequilibrio en muchas de las fases (p.ej. zonaciones complejas, bordes corroídos y recrecimientos en plagioclasas, bordes de opacita y coronas de reacción en anfíboles), sugieren una evolución relativamente compleja. La posibilidad más razonable, siendo además el escenario más habitualmente observado en el volcanismo monogenético de Michoacán (ver p.ej. Cebriá et al., 2011; Corona-Chávez et al., 2006; Chesley et al., 2002; Luhr & Carmichael, 1985; Urrutia-Fucugauchi & Uribe-Cifuentes, 1999), sugiere a una evolución en sistema abierto, que podría implicar tanto la asimilación de material xenolítico como inyecciones de nuevos aportes de magma en el sistema. Por este motivo resulta especialmente adecuada la aproximación adoptada aquí, considerando exclusivamente aquellos minerales que representan fases en equilibrio con el magma que los engloba y cuya composición se asume respresentada por la composición de la roca total. Obviamente, esto implica que los resultados termobarométricos obtenidos se refieren al momento final de la cristalización de cada fase y por tanto es posible que no aporten información sobre las etapas previas de cristalización, cuyo estudio requeriría establecer la composición de los respectivos líquidos en equilibrio.
En la fig. 6A se presentan las temperaturas obtenidas para cada una de las fases respecto al contenido en SiO2 de la roca, que se puede asumir como indicador del grado de diferenciación. Como se puede observar, la secuencia de cristalización se verifica en un relativamente estrecho margen de temperatura, comenzando con la aparición de Ol (1232−1198 °C), seguido de Pl (1162–1126 °C), Opx (1147−1027 °C) y Cpx (1147–1018 °C). En esta secuencia, el Ol deja de cristalizar con la aparición del Opx y los fenocristales de Pl con la del Cpx. Por su parte, el Amp cristalizaría en una última etapa, entre los 995°C y 922°C.
Los datos de P obtenidos (Fig. 6B) por los diferentes geobarómetros, indican que la cristalización se verifica entre los ~7 y 0 kbar, con una correspondencia directa entre P y T lo que sugiere que el proceso de diferenciación tuvo lugar desde ~25 km de profundidad hasta niveles superficiales.
Es de destacar que, como se ha descrito, los cristales de anfíbol presentan en todos los casos evidencias de desequilibrio, a pesar de lo cual se han podido calcular valores de T, P y contenido en agua gracias a la aplicación de un geotermobarómetro basado únicamente en la composición del mineral, por lo que estos datos solo aportan las condiciones de cristalización de los Amp, independientemente de que se encuentren o no en equilibrio con la roca huesped. En consecuencia, cabe la posibilidad de que los cristales de anfíbol representen xenocristales y que por tanto no se encuentren en equilibrio con el magma que los engloba. La procedencia de estos xenocristales podría resultar bien de la asimilación de las rocas atravesadas por el magma durante su ascenso a la superficie o bien de la mezcla con un pulso de magma más silíceo y de menor temperatura. Por otro lado, las condiciones geotermobarométricas calculadas y las variaciones en elementos mayores, también podrían explicarse si el Amp se genera en la última etapa en la secuencia de cristalización de estos magmas. Este supuesto se vería apoyado por la existencia de Amp en las dos tendencias de diferenciación identificadas, apareciendo siempre en los términos más evolucionados (con menor contenido en MgO) y en relación directa con un característico empobrecimiento en K2O (que sería consecuencia del fraccionamiento de esta fase). Asimismo, esta hipotesis se vería reforzada ya que la cristalización de Amp está limitada a presiones entre 5 y 2.5 kbar. Esto podría explicar su desequilibrio al alcanzar niveles superficiales, donde se produce la desgasificación del magma y una reducción drástica del contenido en H2O.
Al igual que en el caso del Amp, las evidencias de desequilibrio de otras fases pueden explicarse por desestabilización de dichas fases tras su cristalización durante etapas previas de la diferenciación, al variar las condiciones fisicoquímicas del magma en etapas posteriores de su evolución. Por ejemplo, en el caso de la Pl la desestabilización puede ser debida al incremento en el contenido de agua durante la cristalización, que como hemos visto se sitúa en torno al 1.6% durante su cristalización y alcanzaría hasta el 6.9% durante la cristalización de Amp. Esta es una situación relativamente frecuente en magmas calcoalcalinos y se ha observado también en volcanes cercanos como el Paricutin (Cebriá et al., 2011), donde su desaparición como fenocristal coincidiendo con la aparición de Opx, se interpretó como consecuencia del aumento en H2O, en este caso debido a la entrada de agua en el sistema por asimilación de un componente cortical.
En resumen, los datos preliminares obtenidos en este estudio permiten establecer las condiciones P-T-H2O bajo las que ha tenido lugar el proceso de diferenciación de las lavas de El Metate. Así, la secuencia de cristalización Ol ± Pl-Opx-Cpx( ± Amp?) parece iniciarse a unos 25 km de profundidad y temperaturas en torno a 1230 °C. El aumento del contenido en H2O que se produce en el sistema al progresar la cristalización pudo haber provocado en primer lugar que la Pl deje de ser una fase estable en el sistema y que cuando en torno a los 990°C se alcanza un porcentaje de agua de ~7% se inicie la cristalización del Amp, que dejará de ser estable a niveles superficiales. No obstante las texturas de desequilibrio que presentan estas dos fases minerales sugieren que algunos cristales de Pl y los de Amp pueden tener un origen xenolítico.
Estos datos preliminares constituyen la base necesaria para perfilar un modelo petrogenético, que permita definir con precisión los parámetros del proceso, pero requieren resolver previamente algunas de las incóginas planteadas y que se encuentran actualmente en estudio. Así por ejemplo, es necesario establecer precisiones sobre las variaciones composicionales y termobarométricas a nivel de cristal en aquellos casos donde se observan zonaciones, lo que requiere además establecer la composición de sus correspondientes líquidos en equilibrio. Esto permitiría confirmar o descartar si en algún caso se trata de xenocristales y si se ha producido o no la entrada de nuevos aportes de magma en el sistema. Esta situación también podrá verificarse a partir del estudio de los microcristales presentes en coronas de reacción, como las observadas en Amp. Finalmente, contar con un mayor número de datos permitirá confirmar que realmente no existen diferencias termobarométricas entre las dos pautas de diferenciación observadas a partir de los elementos mayores de las rocas, como parecen indicar los primeros datos aquí expuestos.
AGRADECIMIENTOS
Este trabajo ha sido financiado mediante el proyecto CGL2011-23422 del Ministerio de Economía y Competitividad de España. Asimismo se agradece la concesión a E. Losantos de una beca FPI (BES-2012-052390) asignada a dicho proyecto.
Agradecemos a Rufino Lozano-Santa Cruz y Laura Luna González, del Instituto de Geología (UNAM) por los análisis de FRX y por el apoyo en los trabajos de campo respectivamente, y a Alfredo Fernández Larios del Centro Nacional de Microscopía Electrónica (UCM) por su asistencia en los análisis por microsonda electrónica. Gracias especiales a Francisco Pérez Torrado (Universidad de Las Palmas de Gran Canaria) y a José Luis Arce (Instituto de Geología, UNAM) por la revisión del manuscrito original, que han mejorado el trabajo y han aportado nuevas ideas para futuras investigaciones.
ReferenciasBeattieP.Olivine-melt and orthopyroxene-melt equilibriaContributions to Mineralogy and Petrology19931151103111http://dx.doi.org/10.1007/BF00712982.CebriáJ.M.MartinyB.M.López-RuizJ.Morán-ZentenoD.J.The Parícutin calc-alkaline lavas: New geochemical and petrogenetic modelling constraints on the crustal assimilation processJournal of Volcanology and Geothermal Research20112011–4113125http://dx.doi.org/10.1016/j.jvolgeores.2010.11.011.Corona-ChávezP.Reyes-SalasM.Gardu-o-MonroyV.H.Israde-AlcántaraI.Lozano-Santa CruzR.Morton-BermeaO.Hernández-ÁlvarezE.Assimilation of granitic xenoliths in the Michoacán-Guanajuato volcanic field: The case of Arócutin, Michoacán, MexicoRevista Mexicana de Ciencias Geológicas2006232233245CostaF.AndreastutiS.Bouvet de MaisonneuveC.PallisterJ.S.Petrological insights into the storage conditions, and magmatic processes that yielded the centennial 2010 Merapi explosive eruptionJournal of Volcanology and Geothermal Research2013261209235http://dx.doi.org/10.1016/j.jvolgeores.2012.12.025.ChesleyJ.RuizJ.RighterK.FerrariL.Gómez-TuenaA.Source contamination versus assimilation: An example from the Trans-Mexican volcanic arcEarth and Planetary Science Letters20021953–4211221http://dx.doi.org/10.1016/S0012-821X(01)00580-5.DahrenB.TrollV.AnderssonU.ChadwickJ.GardnerM.JaxybulatovK.KoulakovI.Magma plumbing beneath Anak Krakatau volcano, Indonesia: evidence for multiple magma storage regionsContributions to Mineralogy and Petrology20121634631651http://dx.doi.org/10.1007/s00410-011-0690-8.FerrariL.López-MartínezM.Aguirre-DíazG.Carrasco-NúñezG.Space-time patterns of Cenozoic arc volcanism in central Mexico: From the Sierra Madre Occidental to the Mexican Volcanic BeltGeology1999274303306http://dx.doi.org/10.1130/0091-7613(1999)027<0303:STPOCA>2.3.CO;2.FerrariL.Orozco-EsquivelT.ManeaV.ManeaM.The dynamic history of the Trans-Mexican Volcanic Belt and the Mexico subduction zoneTectonophysics2012522–523122149http://dx.doi.org/10.1016/j.tecto.2011.09.018.GhiorsoM.S.HirschmannM.M.ReinersP.W.KressV.C.The pMELTS: A revision of MELTS for improved calculation of phase relations and major element partitioning related to partial melting of the mantle to 3 GpaGeochemistry, Geophysics, Geosystems200235135http://dx.doi.org/10.1029/2001GC000217.Gómez-TuenaA.Orozco-EsquivelM.T.FerrariL.Alaniz-ÁlvarezS.A.Nieto-SamaniegoA.F.>
Igneous petrogenesis of the Trans-Mexican Volcanic BeltGeological Society of America Special Paper 4222007Geological Society of America129181http://dx.doi.org/10.1130/2007.2422(05).GroveT.L.JusterT.C.Experimental investigations of low-Ca pyroxene stability and olivine-pyroxene-liquid equilibria at 1-atm in natural basaltic and andesitic liquidsContributions to Mineralogy and Petrology19891033287305http://dx.doi.org/10.1007/BF00402916.HasenakaT.Size, distribution, and magma output rate for shield volcanoes of the Michoacán-Guanajuato volcanic field, Central MexicoJournal of Volcanology and Geothermal Research1994631–21331http://dx.doi.org/10.1016/0377-0273(94)90016-7.HasenakaT.CarmichaelI.S.E.The cinder cones of Michoacán-Guanajuato, central Mexico: their age, volume and distribution, and magma discharge rateJournal of Volcanology and Geothermal Research1985251–2105124http://dx.doi.org/10.1016/0377-0273(85)90007-1.HasenakaT.CarmichaelI.S.E.The Cinder Cones of Michoacán-Guanajuato, Central Mexico: Petrology and ChemistryJournal of Petrology1987282241269http://dx.doi.org/10.1093/petrology/28.2.241.JohnsonC.A.HarrisonC.G.A.Tectonics and volcanism in Central Mexico: A landsat thematic mapper perspectiveRemote Sensing of Environment198928273286http://dx.doi.org/10.1016/0034-4257(89)90119-3.JohnsonE.R.WallaceP.J.CashmanK.V.Delgado GranadosH.Degassing of volatiles (H2O, CO2, S, Cl) during ascent, crystallization, and eruption at mafic monogenetic volcanoes in central MexicoJournal of Volcanology and Geothermal Research20101971–4225238http://dx.doi.org/10.1016/j.jvolgeores.2010.02.017.KeidingJ.K.SigmarssonO.Geothermobarometry of the 2010 Eyjafjallajökull eruption: New constraints on Icelandic magma plumbing systemsJournal of Geophysical Research: Solid Earth2012117B9B00C09KinzlerR.J.GroveT.L.Primary magmas of mid-ocean ridge basalts 1. Experiments and methodsJournal of Geophysical Research199297B568856906http://dx.doi.org/10.1029/91JB02840.Le BasM.J.Le MaitreR.W.WoolleyA.R.The construction of the Total Alkali-Silica chemical classification of volcanic rocksMineralogy and Petrology1992461122http://dx.doi.org/10.1007/BF01160698.Lozano-Santa CruzR.BernalJ.P.Characterization of a new set of eight geochemical reference materials for XRF major and trace element analysisRevista Mexicana de Ciencias Geológicas2005223329344LuhrJ.F.CarmichaelI.Jorullo Volcano, Michoacán, Mexico (1759–1774): The earliest stages of fractionation in calc-alkaline magmasContributions to Mineralogy and Petrology1985902142161http://dx.doi.org/10.1007/BF00378256.MathezE.Refinement of the Kudo-Weill plagioclase thermometer and its application to basaltic rocksContributions to Mineralogy and Petrology19734116172http://dx.doi.org/10.1007/BF00377654.Patiño DouceA.E.Vapor-absent melting of tonalite at 15–32 kbarJournal of Petrology2005462275290http://dx.doi.org/10.1093/petrology/egh071.PlechovP.Y.TsaiA.E.ShcherbakovV.D.DirksenO.V.Opacitization conditions of hornblende in Bezymyannyi volcano andesites (March 30, 1956 eruption)Petrology20081611935http://dx.doi.org/10.1134/S0869591108010025.PutirkaK.JohnsonM.KinzlerR.LonghiJ.WalkerD.Thermobarometry of mafic igneous rocks based on clinopyroxene-liquid equilibria, 0–30 kbarContributions to Mineralogy and Petrology1996123192108http://dx.doi.org/10.1007/s004100050145.PutirkaK.D.Igneous thermometers and barometers based on plagioclase + liquid equilibria: Tests of some existing models and new calibrationsAmerican Mineralogist2005902–3336346http://dx.doi.org/10.2138/am.2005.1449.PutirkaK.D.Thermometers and Barometers for Volcanic SystemsReviews in Mineralogy and Geochemistry200869161120http://dx.doi.org/10.2138/rmg.2008.69.3.RhodesJ.M.DunganM.A.BlanchardD.P.LongP.E.Magma mixing at mid-ocean ridges: Evidence from basalts drilled near 22°N on the mid-Atlantic RidgeTectonophysics1979551–23561http://dx.doi.org/10.1016/0040-1951(79)90334-2.RidolfiF.RenzulliA.PueriniM.Stability and chemical equilibrium of amphibole in calc-alkaline magmas: an overview, new thermobarometric formulations and application to subduction-related volcanoesContributions to Mineralogy and Petrology201016014566http://dx.doi.org/10.1007/s00410-009-0465-7.RoederP.L.EmslieR.F.Olivine-liquid equilibriumContributions to Mineralogy and Petrology197029275289http://dx.doi.org/10.1007/BF00371276.ScailletB.MacDonaldR.Experimental constraints on the relationships between peralkaline rhyolites of the Kenya Rift ValleyJournal of Petrology2003441018671894http://dx.doi.org/10.1093/petrology/egg062.SissonT.W.GroveT.L.Experimental investigations of the role of H2O in calc-alkaline differentiation and subduction zone magmatismContributions to Mineralogy and Petrology1993a1132143166http://dx.doi.org/10.1007/BF00283225.SissonT.W.GroveT.L.Temperatures and H2O contents of low-MgO high-alumina basaltsContributions to Mineralogy and Petrology1993b1132167184http://dx.doi.org/10.1007/BF00283226.StreckM.J.Mineral Textures and Zoning as Evidence for Open System ProcessesReviews in Mineralogy and Geochemistry2008691595622http://dx.doi.org/10.2138/rmg.2008.69.15.SuterM.ContrerasJ.Gómez-TuenaA.SiebeC.Quintero-LegorretaO.García-PalomoA.MacíasJ.L.Alaniz-ÁlvarezS.A.Nieto-SamaniegoA.F.FerrariL.Effect of strain rate in the distribution of monogenetic and polygenetic volcanism in the Transmexican volcanic belt: Comments and ReplyGeology1999276571575http://dx.doi.org/10.1130/0091-7613(1999)027<0571:EOSRIT>2.3.CO;2.Urrutia-FucugauchiJ.Uribe-CifuentesR.M.Lower-crustal xenoliths from the Valle de Santiago maar field, Michoacan-Guanajuato volcanic field, central MexicoInternational Geology Review1999411210671081http://dx.doi.org/10.1080/00206819909465192.WalterM.J.PresnallD.C.Melting behavior of simplified lherzolite in the system CaO-MgO-Al2O3-SiO2-Na2O from 7 to 35 kbarJournal of Petrology1994352329359http://dx.doi.org/10.1093/petrology/35.2.329.