Introducción
Colombia sobresale en el mercado de las gemas por la alta calidad de sus esmeraldas; su finura se deriva del carácter particular de sus yacimientos, originados en un contexto geológico poco común. Los depósitos colombianos representan hasta el momento, la única fuente de esmeraldas encontradas en una cuenca sedimentaria afectada por hidrotermalismo no asociado a magmatismo (Cheilletz y Giuliani, 1996; Giuliani et al., 2019; Groat et al., 2008). Aunque existen reportes de otras localidades con ocurrencias en escenarios geológicos similares (Giuliani et al., 2019; Hewton et al., 2013), solo los depósitos de Colombia mantienen una importancia económica e histórica. Sin embargo, al observar las cifras de producción de las últimas décadas, existe una constante tendencia a la baja. En 2005 se produjeron ~6.746.000 quilates y en 2019 ~654.500 quilates (UPME, 2020). Ante este panorama, se hace necesario intensificar la exploración de tal manera que nuevos descubrimientos puedan contribuir al aumento de la producción. Para ello es imprescindible disponer de un conocimiento geológico detallado de las minas actuales; partiendo del principio de que no se puede buscar lo que no se conoce. A pesar de que la comunidad científica ha contribuido con abundantes publicaciones sobre el origen y formación de los yacimientos esmeraldíferos, los estudios son generales y existen pocos reportes sobre los rasgos geológicos de una mina en particular. El objetivo de este trabajo es proveer información geoquímica, mineralógica y estructural sobre la mineralización de esmeraldas de la mina La Pava, Muzo-Quípama; para promover la definición de parámetros geológicos que conlleven a realizar una exploración satisfactoria.
Contexto geológico
Los depósitos esmeraldíferos en Colombia se distribuyen a lo largo de dos corredores en la Cordillera Oriental de los Andes (Moreno et al., 2009; Reyes et al., 2006; Terraza, 2019). La actual configuración tectónica de la Cordillera Oriental está ligada a diferentes procesos que han tenido lugar especialmente durante los últimos 150 Ma (Horton et al., 2020); se destaca la extensa acumulación de sedimentos principalmente marinos durante el Cretácico inferior, como parte del relleno de una cuenca de retro-arco (Teixell et al., 2015). En el Cretácico superior la dinámica extensional se vio interrumpida y la sedimentación empezó a ser controlada por el inicio del acortamiento cortical de la cuenca. Durante el Cenozoico, la compresión continuó y dio origen a secuencias predominantemente continentales. Finalmente, se dio paso a la exhumación y exposición de la actual cadena montañosa desde el Neógeno (Bayona, 2018; Horton et al., 2020; Kammer et al., 2020; Mora, 2007). Como resultado, la Cordillera Oriental está dominada por la presencia de sedimentitas, con exposiciones menores del basamento metamórfico y registros de actividad magmática muy restringidos (Bernet et al., 2016; Vásquez et al., 2010).
El distrito Muzo-Quípama es el principal productor de esmeraldas en la actualidad y se localiza a pocos kilómetros al NW del municipio de Muzo (Figura 1); hace parte del Cinturón Esmeraldífero Occidental junto con los distritos de Maripí, Coscuez, Peñas Blancas y Yacopí-La Palma (Reyes et al., 2006; Terraza, 2019). La mineralización de esmeraldas se asocia con pelitas y calizas del Valanginiano al Barremiano (Cheilletz y Giuliani, 1996; Giuliani et al., 1995; Mantilla et al., 2007; Reyes et al., 2006; Terraza, 2019). En la región de Muzo-Quípama la secuencia estratigráfica, en orden cronológico (de más antiguo a más joven), está compuesta por: la Formación Rosablanca que consta de una sucesión mayoritariamente calcárea de aproximadamente 1000 m de espesor; sobre ella descansan capas delgadas de lodolitas, limolitas y arcillolitas de la Formación Furatena, cuyo espesor oscila en 1000 m. Posteriormente se encuentra la Formación Muzo, la cual corresponde a un paquete de 200-300 m de espesor conformado por lodolitas carbonosas, lodolitas calcáreas y niveles delgados y esporádicos de calizas. Finalmente, sobre la Formación Muzo se encuentran capas gruesas de limolitas y lodolitas silíceas de la Formación Capotes con 700 m de espesor (Reyes et al., 2006; Terraza, 2019). De acuerdo con diversos autores, en la zona de estudio la ocurrencia de esmeraldas se restringe exclusivamente a niveles de lodolitas calcáreas de la Formación Muzo (Mantilla et al., 2007; Moreno et al., 2009; Pignatelli et al., 2015; Reyes et al., 2006; Terraza, 2019). Desde el punto de vista tectónico, el distrito Muzo-Quípama se encuentra bajo la influencia de estructuras desarrolladas en un régimen transpresivo (Branquet et al., 1999; Laumonier et al., 1999; Reyes et al., 2006); la falla del Río Minero, la falla de Itoco y la falla de Peñas Blancas, son las estructuras mayores que configuran el área (Figura 1). El Anticlinal Alto La Chapa-Borbúr constituye un pliegue regional que ejerce influencia en la localización de zonas mineralizadas (Reyes et al., 2006). En La Pava se ha descrito una estructura de tipo flor positiva que permite la exposición de rocas de la Formación Muzo (Branquet et al., 1999); a escala más detallada, las labores mineras están ubicadas en un dominio estructural caracterizado por pliegues apretados, particularmente el Anticlinal y Sinclinal de La Pava (Figura 2), los cuales se desarrollan en el bloque colgante de la falla El Águila, una estructura de cabalgamiento de alto ángulo con vergencia al occidente (González-Durán et al., 2019).
Metodología
Se realizó un levantamiento estratigráfico a escala 1:100, para lo cual se utilizó un bastón de Jacob con acople de rayo láser para medir los espesores de las capas. Se estudiaron dos secciones estratigráficas, una a lo largo del cauce del río Itoco y la otra sobre la quebrada La Pava. En ambas columnas se tomaron muestras de roca en intervalos de dos metros. Se seleccionaron 101 muestras de rocas, provenientes de las tres poblaciones -río Itoco, quebrada La Pava y túneles mineros- para determinar la composición química. Los análisis se realizaron en los laboratorios del Servicio Geológico Colombiano. Para la cuantificación de los óxidos de elementos mayores -SiO₂, Al₂O₃, TiO₂, Fe₂O₃, MgO, CaO y Na₂O-, se emplearon muestras fundidas con una mezcla de meta y tetraborato de litio, mientras que para los elementos menores -Zr, Ge, Nb, Mo, Sn, Sb y Hf- se utilizaron muestras prensadas. La cuantificación se llevó a cabo en un equipo de fluorescencia de rayos X, FRX Panalytical Axios Mineral. Los análisis de elementos trazas -Ni, Cu, As, Li, Cs, Be, V, Cr, Ba, Mn, Sc, Co, Zn, Ga, Rb, Sr, Cd, In, Y, Tl, Pb, Bi, Th y U- se realizaron en un espectrómetro de masas con plasma acoplado inductivamente, ICP-MS Perkin Elmer NEXION, luego de una digestión multiácida convencional (4 ácidos: HF, HNO3, HCl y HClO4).
Adicionalmente, se realizaron análisis petrográficos de 35 secciones delgadas pulidas en un microscopio Olympus BX-51. Al tratarse de rocas pelíticas, se acudió a los diagramas de clasificación de Folk (1980) para rocas con componentes terrígenos. Así mismo, se seleccionaron 55 cristales de esmeraldas para análisis químicos mediante fluorescencia de rayos X. La cuantificación se realizó en un equipo de energía dispersiva PANalytical Epsilon 5. Ambos procedimientos se efectuaron en las instalaciones del Centro de Desarrollo Tecnológico de la Esmeralda Colombiana.
Resultados
Mineralización
La unidad predominante en el área es la Formación Muzo, presentándose como una secuencia dominada por limolitas negras monótonas. Se observan abundantes afloramientos a lo largo del cauce del río Itoco, la quebrada La Pava y el antiguo tajo de explotación a cielo abierto. La deformación tectónica es intensa, evidenciada por un fuerte clivaje -que generalmente enmascara la dirección de estratificación-, pliegues apretados y abundantes zonas de brechamiento. La mineralización de esmeraldas está asociada con zonas de brechas y venas con tendencia N10E, las cuales se disponen de manera paralela al plano axial del Anticlinal de La Pava; esto como resultado de un emplazamiento por deslizamiento flexural intracapas. La mayor concentración de venas productivas se distribuye en el núcleo del anticlinal y se proyecta hacia la charnela (Figura 2).
En el río Itoco aflora una capa de limolitas de color gris claro a beige, que exhiben una fuerte alteración hidrotermal (albitización y carbonatización) y que han sido afectadas por brechamiento tectónico; ese nivel se denomina “albititas” (Figura 3A). A escala de afloramiento estas rocas se presentan como un nivel de brechas estratiformes de seis metros de espesor; son de carácter polimíctico puesto que contienen clastos de rocas blancas albitizadas y limolitas negras levemente alteradas (Figura 3B); así mismo, algunos fragmentos de albititas presentan plegamientos tipo chevron (Figura 3C). El diámetro de los clastos varía entre 2 a 120 cm, con textura caótica y matriz compuesta por material triturado. Las albititas tienen forma lenticular y su espesor puede disminuir hasta alcanzar unos pocos centímetros -como es el caso en la quebrada La Pava-. Al interior de los túneles las albititas se encuentran altamente tectonizadas, son comunes los replegamientos, brechamientos y saltos por fallas locales (Figura 3D). Por el alto contraste de color que representan, las albititas son usadas como niveles guía por los mineros, debido a que las venas mineralizadas generalmente ocurren en zonas aledañas.
De manera local, se observan brechas hidrotermales con matriz de albita, cuyos cristales alcanzan los 2 cm de diámetro -exhiben textura “grano de arroz”-, y en menor proporción carbonatos, pirita y barita; estas brechas afectan las limolitas negras debajo de las capas de albititas. Aunque las brechas son locales, las zonas de albitización son de amplia extensión. Gran parte de las zonas albitizadas se encuentran sobreimpuestas por vetillas y venas de carbonatos -calcita y dolomita- con ocurrencia esporádica de esmeraldas (Figura 3E). Estas venas son a menudo drusiformes y el tamaño de los cristales es mayor a 1 cm de diámetro; las esmeraldas suelen ocurrir en espacios abiertos -en compañía de calcita, dolomita y albita- y típicamente se presentan como prismas hexagonales y di-hexagonales, en ocasiones con terminaciones de múltiples facetas con pirámides hexagonales y bi-hexagonales (Figura 4). Las esmeraldas de La Pava se caracterizan por su alta claridad y brillo, mientras que su color varía desde un verde levemente saturado a moderado, siendo el verde intenso muy poco común.
Las venas portadoras de esmeraldas no suelen superar los 25 cm de espesor, y con frecuencia se observan fenómenos de reapertura de antiguas venas de calcita fibrosa o albita. En algunos sectores hay ocurrencia de venas drusiformes de cuarzo y dolomita que pueden alcanzar los 120 cm de espesor. Ocasionalmente se encuentran sulfuros como pirita y calcopirita, aunque en bajas cantidades; también suelen encontrarse minerales tardíos que rellenan espacios tales como calcita -con habito “granizo”- y cristales euhedrales de barita dorada.
A nivel petrográfico las rocas hospedantes de la mineralización varían entre limolitas de grano medio a fino, arcillolitas y lodolitas con abundante contenido de materia orgánica. Se observan algunos niveles fosilíferos con presencia de foraminíferos planctónicos; también son comunes las intercalaciones entre niveles con distintos tamaños de grano y la afectación de estructuras sedimentarias por microdeformación (Figura 5A). A nivel composicional los principales constituyentes son las plagioclasas, micas y granos de cuarzo; como minerales traza se identificaron monacitas, zircones, rutilos, agregados de pirita framboidal y carbonatos.
Las rocas mineralizadas, especialmente las muestras colectadas en los túneles exhiben dos estilos de alteración hidrotermal muy definidos. El primer estilo corresponde a la albitización y se expresa mediante un intenso vetilleo y brechamiento hidrotermal, así como el reemplazamiento de los minerales formadores de roca (plagioclasas y micas) por albita. Las albitas se presentan comúnmente como cristales tabulares y maclados con crecimiento drusiforme (Figura 5B). Es claro distinguir múltiples generaciones de albita que se sobreimponen entre sí; en algunos sectores se observan vetillas que previamente contenían moscovita que han sido reabiertas por inyecciones de albita (Figura 5C). Si bien, la albita es el mineral predominante (más del 99%), también suelen encontrarse esporádicamente cristales de xenotimo, pirita dodecaédrica, turmalina, rutilo y fluorapatito. Este último puede exhibir cristales de gran tamaño (1-3 cm de diámetro), de color rosado característico de las zonas esmeraldíferas. En las albititas, la albitización es intensa y consiste en un reemplazamiento casi total de la roca original. El reemplazamiento genera blanqueamiento, en el sentido que el color de la nueva roca es generalmente gris claro o beige, un hecho que contrasta con el color negro original de las limolitas (Figura 3).
El segundo estilo de alteración hidrotermal corresponde a la carbonatización. Este evento es posterior a la albitización y se caracteriza por la abundante ocurrencia de carbonatos: dolomita y calcita. Ambos minerales se encuentran como cristales romboédricos al interior de venas, conformando la matriz de brechas hidrotermales o como cristales diseminados en las rocas albitizadas (Figura 5D, 5E). Es común observar cristales de albita con bordes corroídos por el reemplazamiento parcial por carbonatos; también ocurren otras texturas indicativas de reemplazamiento como pseudomorfos e inclusiones relícticas. A diferencia de la albitización que consta de vetillas delgadas, los carbonatos suelen encontrarse en venas con espesores que oscilan entre los 2 y los 25 cm. Las rocas circundantes a las venas generalmente reaccionan ante el ácido clorhídrico por el alto contenido de calcita.
Otros minerales que también hacen parte de este ensamblaje, pero que ocurren en mucho menor proporción, son cuarzo, moscovita, albita, pirita dodecaédrica, calcopirita, esfalerita y esmeralda (Figura 5F, 5G). La ocurrencia de esmeraldas en vetillas compuestas únicamente por albita es casi nula. En las albititas, la carbonatización generalmente está acompañada por la ocurrencia de azufre nativo, junto con cristales de cuarzo, los cuales particularmente contienen abundantes inclusiones de fragmentos de rocas, cristales de moscovita y fragmentos de anhidrita que posiblemente sean relictos de la etapa de albitización.
Finalmente, se presenta una etapa hidrotermal tardía caracterizada por la ocurrencia de barita, siderita y yeso (Figura 5H). Estos minerales se encuentran, en su mayoría, como relleno de espacios abiertos y por lo tanto exhiben formas euhedrales. Por otro lado, la oxidación supérgena ha dado origen a la ocurrencia de goethita, malaquita y covelina a partir de la alteración de pirita y calcopirita. Los carbonatos que contienen hierro (dolomita y siderita) también se ven afectados por la oxidación y suelen mostrar pátinas, reemplazamientos y textura boxwork con óxidos e hidróxidos de hierro. La Figura 6 muestra una secuencia paragenética esquemática y generalizada de la mineralización de La Pava.
Geoquímica de rocas
Para el análisis geoquímico de roca total, las muestras se dividieron en tres grupos con base en el criterio geológico de las observaciones en campo. El primer grupo abarca 68 muestras de limolitas negras recolectadas en la ribera del río Itoco, las cuales corresponden a una zona con poca manifestación de alteración hidrotermal (a excepción del nivel de albititas) y por lo tanto se asume que los resultados representan valores cercanos a la composición original de las rocas antes de la mineralización. El segundo grupo corresponde a 19 muestras de limolitas negras recolectadas en la quebrada La Pava; allí se observaron rasgos claros de alteración hidrotermal, entre ellos zonas de vetilleo y brechas hidrotermales. En tercer lugar, se encuentran 14 limolitas negras fuertemente alteradas hidrotermalmente, adquiridas al interior de los túneles en cercanías a venas portadoras de esmeraldas; por lo tanto, constituyen testigos de primera mano del metasomatismo responsable de formar las esmeraldas. Las Tablas 1 y 2 muestran un resumen estadístico de los resultados de los análisis geoquímicos para los tres grupos de muestras; los datos pueden ser consultados en los Anexos 1 y 2.
En general, las muestras están dominadas por las concentraciones de SiO₂ (40,6-57,1%) y Al₂O₃ (11,3- 22,6%) con cantidades menores de otros elementos como Fe₂O₃, Na₂O, K₂O, CaO, TiO₂ y MgO, mientras que la concentración promedio de P₂O₅ no sobrepasa el 1%. La variación de la mayoría de los datos es alta en cada grupo de muestras, en especial en elementos de fácil movilidad como Ca, Fe y Mg donde la desviación estándar supera en magnitud en algunos casos al valor medio. El promedio de pérdida de componentes por ignición (LOI) oscila entre 8,5% y 18,8%, lo cual sugiere un aporte considerable de componentes orgánicos y elementos volatilizables. En cuanto a elementos traza, la variación sigue siendo alta para la mayoría, principalmente en elementos litófilos como Mn y V, y aún más acentuada en elementos calcófilos como Mo, Zn, Cd, Ni, y As. La Figura 7 muestra dos columnas litoestratigráficas levantadas en el río Itoco y la quebrada La Pava; estas incluyen una representación gráfica de los porcentajes relativos de elementos mayores y registros de algunos elementos trazas (Be, Cr/V, Mn e Y). La sección del río Itoco abarca la parte superior de la Formación Muzo y la base de la Formación Capotes; el contacto se encuentra aproximadamente en el metro 150. El atributo más sobresaliente está relacionado con el nivel de albititas comprendido entre los 130-140 m. El principal rasgo que diferencia este nivel de toda la sección es el aumento en la concentración de Ca, Mg, Mn e Y. De manera recíproca, las concentraciones de Ti, K y Be disminuyen considerablemente. La relación Cr/V se mantiene relativamente constante en la parte superior de la sección, mientras que en la periferia de las albititas y en la parte inferior exhibe fluctuaciones. En adición a los elementos expuestos en la gráfica, casi la totalidad de los elementos analizados muestran cambios notables en las albititas. Por otra parte, el Fe y el Ca presentan un aumento gradual hacía la parte alta de la sección, el cual contrasta con una disminución en el contenido de Na.
En la sección de la quebrada La Pava, se observa una mayor abundancia en las concentraciones de Ca, Na e Y en comparación con el río Itoco; de igual forma, las oscilaciones de la mayor parte de los elementos también son más recurrentes. En el intervalo comprendido entre 20-110 m (el bloque más afectado por la actividad hidrotermal), se observa un incremento sustancial en los valores de Ca, localmente de Mg, y una disminución drástica en el contenido de Na; sin embargo, hacia la parte superior, la concentración de Na aumenta de tal manera que alcanza proporciones cercanas al 95%, específicamente en el metro 100. El Be exhibe empobrecimientos en varios sectores sin mostrar zonas de incrementos sustanciales (un comportamiento similar al del Cr y V), mientras que el Y tiende a mostrar incrementos con valores cercanos a los alcanzados en el nivel de albititas del río Itoco. El Mn muestra un enriquecimiento en el intervalo de 40-50 m, el cual se correlaciona de manera directa con el Mg. La relación Cr/V presenta fluctuaciones a lo largo de toda la sección. Los contenidos de Fe y Ti no presentan variaciones considerables ni se observa una tendencia.
La Figura 8 muestra una comparación gráfica entre la composición promedio de las rocas del río Itoco, normalizada con respecto a la composición de la corteza continental superior provista por Rudnick y Gao (2014). Gran parte de los elementos exhiben ligeras variaciones; sin embargo, las anomalías más notables corresponden al enriquecimiento de Mo - más de ≈50 veces- junto con aumentos menores en los contenidos de As, Cd, V y Bi. Por otro lado, los empobrecimientos más notables corresponden a los valores de Mg, Mn, Co, Sr, Y y Pb. Cabe resaltar que el Be, elemento indispensable para la formación de las esmeraldas, exhibe tan solo un ligero aumento -alrededor de 2 veces- con respecto al valor medio de la corteza.
En adición, se construyó un diagrama multielementos para visualizar gráficamente las principales anomalías geoquímicas. Se tomó como línea base o valor de fondo los valores promedio del río Itoco. El eje “y” representa la diferencia entre el promedio de los valores del río Itoco, La Pava y los túneles, sobre la desviación estándar del río Itoco (Figura 9). Las líneas grises punteadas representan dos desviaciones estándar, de modo que los valores que son mayores a 2 son vistos anomalías positivas y los valores menores a -2 como anomalías negativas (Gandhi y Sarkar, 2016).
SiO2 | Al2O3 | TiO2 | Fe2O3T | MgO | CaO | Na2O | K2O | P2O5 | LOI | Total | Zr | Ge | Nb | Mo | Sn | Sb | Hf | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Río Itoco | ||||||||||||||||||
n | 68 | 68 | 68 | 68 | 68 | 42 | 68 | 68 | 51 | 68 | 68 | 68 | 36 | 67 | 57 | 26 | 21 | 23 |
Min | 23,7 | 5,2 | 0,2 | 0,2 | 0,1 | 0,1 | 0,2 | 0,5 | 0,03 | 3,9 | 93,5 | 55 | 2 | 11 | 2 | 8 | 12 | 8 |
Max | 70,4 | 28,5 | 1,5 | 33,8 | 6,0 | 32,8 | 9,3 | 4 | 0,46 | 23,6 | 99,9 | 310 | 4 | 36 | 368 | 16 | 35 | 19 |
Prom | 57,1 | 22,6 | 1,1 | 4,0 | 0,5 | 1,8 | 2,3 | 2,1 | 0,14 | 8,5 | 99,5 | 168 | 3 | 24 | 46 | 9 | 19 | 10 |
Med | 58,6 | 23,7 | 1,1 | 2,4 | 0,3 | 0,2 | 2,0 | 2,2 | 0,12 | 8,0 | 99,7 | 158 | 3 | 25 | 20 | 9 | 15 | 9 |
De | 7,7 | 4,4 | 0,3 | 5,6 | 0,8 | 5,9 | 1,3 | 0,6 | 0,09 | 3,0 | 0,9 | 58 | 1 | 7 | 65 | 2 | 7 | 3 |
Quebrada La Pava | ||||||||||||||||||
n | 19 | 19 | 19 | 19 | 18 | 18 | 19 | 16 | 17 | 19 | 19 | 19 | 10 | 19 | 19 | 2 | 17 | 4 |
Min | 9,7 | 2,9 | 0,1 | 0,2 | 0,1 | 0,1 | 0,1 | 0,1 | 0,03 | 2,6 | 96,8 | 25 | 2 | 5 | 3 | 8 | 14 | 8 |
Max | 70,5 | 24,2 | 1,1 | 13,7 | 14,4 | 27,2 | 10,6 | 4,9 | 0,58 | 38,4 | 99,9 | 178 | 5 | 25 | 815 | 10 | 31 | 11 |
Prom | 56,1 | 16,2 | 0,6 | 3,6 | 1,1 | 3,6 | 3,2 | 2,4 | 0,27 | 12,9 | 99,3 | 106 | 3 | 15 | 170 | 9 | 19 | 9 |
Med | 58,1 | 16,4 | 0,6 | 3,8 | 0,3 | 0,5 | 1,2 | 2,3 | 0,28 | 12,0 | 99,6 | 105 | 2 | 15 | 110 | 9 | 17 | 9 |
De | 13,0 | 5,2 | 0,3 | 2,9 | 3,3 | 7,2 | 3,4 | 1,2 | 0,20 | 7,0 | 0,8 | 37 | 1 | 5 | 206 | 1 | 5 | 1 |
Túneles | ||||||||||||||||||
n | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | - | 14 | 14 | - | 4 | 10 |
Min | 21,9 | 5,3 | 0,1 | 1,2 | 0,3 | 2,9 | 0,4 | 0,3 | 0,07 | 7,8 | 96,2 | 32 | - | 5 | 2 | - | 12 | 8 |
Max | 59,6 | 18,3 | 1,0 | 8,6 | 9,2 | 28,3 | 7,4 | 4 | 0,62 | 28,8 | 98,7 | 168 | - | 19 | 179 | - | 22 | 125 |
Prom | 40,6 | 11,3 | 0,5 | 4,1 | 3,7 | 13,6 | 3,8 | 1,2 | 0,30 | 18,8 | 97,8 | 85 | - | 12,0 | 31 | - | 17 | 26 |
Med | 40,8 | 11,6 | 0,4 | 4,1 | 2,3 | 15,3 | 3,5 | 0,8 | 0,30 | 19,5 | 97,8 | 90 | - | 12 | 4 | - | 17 | 13 |
De | 11,1 | 3,9 | 0,3 | 1,9 | 3,2 | 8,3 | 2,3 | 1,1 | 0,16 | 6,5 | 0,8 | 40 | - | 5 | 58 | - | 4 | 35 |
Lc | 0,36 | 0,18 | 0,04 | 0,1 | 0,1 | 0,1 | 0,11 | 0,08 | 0,024 | 0,1 | - | 2 | 2 | 2 | 2 | 8 | 12 | 8 |
Elementos mayores y menores reportados en % en peso; elementos trazas reportados en ppm en peso (Zr, Ge, Nb, Mo, Sn, Sb y Hf). LOI: % pérdida de peso por ignición; Fe2O3T hace referencia al hierro total (Fe2+ + Fe3+). n: número de muestras, Min: valor mínimo, Max: valor máximo, Prom: promedio, Med: mediana, De: desviación estándar y Lc: límite de cuantificación.
Ni | Cu | As | Li | Cs | Be | V | Cr | Ba | Mn | Sc | Co | Zn | Ga | Rb | Sr | Cd | In | Tl | Pb | Bi | Th | U | ||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Río Itoco | ||||||||||||||||||||||||
n | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 | 68 |
Min | 3 | 6 | 2 | 1 | 0,4 | 0,6 | 44 | 17 | 37 | 4 | 6,1 | 0,2 | 6 | 7,1 | 22 | 23 | 0,1 | 0,01 | 1,2 | 0,2 | 2 | 0,1 | 1,4 | 0,9 |
Max | 390 | 78 | 166 | 179 | 15,0 | 7,3 | 2051 | 290 | 714 | 1713 | 41,7 | 64 | 627 | 36,8 | 191 | 173 | 9,0 | 0,10 | 14,2 | 4,3 | 25 | 1,4 | 22,4 | 11,3 |
Prom | 32 | 18 | 10 | 55 | 8,1 | 3,5 | 409 | 111 | 357 | 144 | 16,4 | 7 | 46 | 27,7 | 122 | 108 | 0,6 | 0,05 | 3,6 | 1,4 | 7 | 0,5 | 12,6 | 2,2 |
Med | 22 | 13 | 7 | 56 | 8,2 | 3,4 | 227 | 100 | 365 | 43 | 16,0 | 2 | 19 | 29,1 | 121 | 117 | 0,4 | 0,04 | 2,8 | 1,4 | 6 | 0,5 | 13,3 | 1,7 |
De | 49 | 15 | 20 | 29 | 3,4 | 1,1 | 458 | 50 | 124 | 266 | 4,7 | 10 | 96 | 5,7 | 38 | 37 | 1,1 | 0,02 | 2,5 | 0,7 | 5 | 0,2 | 4,3 | 1,5 |
Quebrada La Pava | ||||||||||||||||||||||||
n | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 18 | 19 | 19 |
Min | 4 | 8 | 2 | 1 | 0,1 | 0,3 | 12 | 5 | 10 | 2 | 0,3 | 0,4 | 6 | 1,6 | 1 | 16 | 0,1 | 0,005 | 0,8 | 0,1 | 2 | 0,1 | 0,8 | 0,3 |
Max | 289 | 99 | 100 | 72 | 10,2 | 4,6 | 2122 | 216 | 587 | 689 | 24,2 | 12 | 273 | 35,4 | 211 | 290 | 9,6 | 0,13 | 20,4 | 6,5 | 66 | 0,9 | 34,1 | 21,2 |
Prom | 110 | 38 | 41 | 19 | 3,2 | 2,2 | 811 | 127 | 286 | 92 | 12,5 | 6 | 72 | 21,3 | 95 | 59 | 1,6 | 0,04 | 7,5 | 2,6 | 15 | 0,4 | 12,7 | 6,5 |
Med | 98 | 41 | 35 | 11 | 2,7 | 2,1 | 717 | 145 | 315 | 23 | 13,0 | 7 | 56 | 22,1 | 102 | 44 | 1,0 | 0,04 | 7,1 | 2,6 | 9 | 0,3 | 10,7 | 4,5 |
De | 74 | 24 | 27 | 21 | 2,7 | 1,2 | 620 | 64 | 169 | 172 | 6,4 | 4 | 68 | 9,1 | 60 | 61 | 2,2 | 0,04 | 5,4 | 1,6 | 145 | 0,3 | 9,9 | 5,6 |
Túneles | ||||||||||||||||||||||||
n | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 14 | 13 | 14 | 14 |
Min | 23 | 10 | 3 | 2 | 0,3 | 0,5 | 35 | 13 | 17 | 100 | 5,9 | 4 | 62 | 8,0 | 7 | 26 | 1,0 | 0,01 | 5,5 | 0,04 | 3 | 0,1 | 2,3 | 0,8 |
Max | 371 | 35 | 106 | 14 | 5,1 | 3,3 | 1897 | 251 | 388 | 1740 | 17,8 | 21 | 2232 | 33,1 | 190 | 111 | 28,4 | 0,10 | 17,3 | 4,7 | 39 | 0,4 | 19,8 | 11,1 |
Prom | 78 | 23 | 21 | 6 | 1,6 | 1,6 | 257 | 59 | 106 | 798 | 10,2 | 10 | 440 | 15,0 | 53 | 69 | 6,9 | 0,04 | 11,5 | 0,7 | 11 | 0,2 | 6,1 | 4,1 |
Med | 42 | 23 | 11 | 6 | 1,0 | 1,3 | 110 | 39 | 84 | 635 | 9,4 | 9 | 323 | 11,9 | 32 | 69 | 3,4 | 0,04 | 11,4 | 0,2 | 8 | 0,2 | 5,1 | 2,8 |
De | 95 | 7 | 27 | 3 | 1,5 | 0,9 | 489 | 62 | 94 | 591 | 3,5 | 4 | 545 | 7,6 | 55 | 26 | 7,5 | 0,03 | 3,0 | 1,2 | 10 | 0,1 | 4,3 | 3,1 |
Lc | 0,5 | 1 | 0,4 | 0,5 | 0,05 | 0,1 | 0,2 | 0,5 | 1 | 0,8 | 0,05 | 0,8 | 4 | 0,2 | 0,2 | 1 | 0,08 | 0,005 | 0,05 | 0,01 | 2 | 0,05 | 0,1 | 0,01 |
n: número de muestras, Min: valor mínimo, Max: valor máximo, Prom: promedio, Med: mediana, De: desviación estándar, Lc: límite de cuantificación.
Al graficar los valores promedio de la quebrada La Pava y los túneles, se aprecia un comportamiento similar en la mayoría de los elementos. Sin embargo, las fluctuaciones se acentúan en un mayor grado en las muestras de los túneles. Las variaciones más notables observadas en las muestras de los túneles corresponden con enriquecimientos en los valores de MgO, CaO, P₂O₅, Zn, Cd, Mn e Y; en contraste con empobrecimientos en SiO₂, Al₂O₃, TiO₂, K₂O, Nb, Li, Cs, Be, Ba y Ga. Pocos elementos permanecen inalterados, tal es el caso del Fe₂O₃ y el Co, cuyo comportamiento no registra fluctuaciones apreciables en ninguno de los dos grupos; de la misma manera, algunos elementos presentan tendencias no uniformes como por ejemplo el K₂O, Mo, V, Cr, Mn, Tl y Th, los cuales se enriquecen en un grupo y empobrecen en otro. Elementos como el Na₂O, Ni, Cu, As, Pb y U presentan ligeros enriquecimientos en ambos grupos, mientras que Zr, Sc, Rb, Sr y Bi exhiben un empobrecimiento sustancial.
Geoquímica de esmeraldas
La Tabla 3 muestra los datos de la composición química de 55 esmeraldas de La Pava. En promedio, las concentraciones de SiO₂ y Al₂O₃ constituyen el 86%; mientras que el contenido de BeO (13,15%) fue añadido manualmente como valor de referencia, tomado de Groat et al. (2008), con el fin de completar el cierre elemental. De esta manera, las concentraciones medias de otros como elementos como K₂O (0,08%), CaO (0,01%) y SO₃ (0,07%), no sobrepasan el 1% del total. El Na₂O solo fue detectado en 5 muestras con valores particularmente altos (3-6%), los cuales probablemente se deban a la presencia de inclusiones de albita. Con respecto a los elementos traza, cinco elementos fueron cuantificados en todas las muestras y sus respectivos promedios son: Cl (1203 ppm), Fe (608 ppm), V (396 ppm), Cr (288 ppm) y Ga (27 ppm). Otros elementos presentes en menores proporciones, con menor frecuencia, incluyen: Sc, Ti, Mn, Ni, Cu y Zn. Los datos de los elementos cromóforos (V, Cr y Fe) reflejan que el Fe es el elemento con la concentración más alta, casi superando por el doble a los demás elementos (Cr y V); a su vez, la relación Cr/V es en promedio 0,71 e indica la predominancia del V sobre el Cr. Con el objetivo de visualizar las proporciones relativas de los principales elementos cromóforos de las muestras de La Pava y compararlas con esmeraldas de otras minas, se elaboraron dos diagramas ternarios (Figura 10). Los datos de los otros distritos fueron tomados de la base de datos del Centro de Desarrollo Tecnológico de la Esmeralda Colombiana publicada en Cedeño et al. (2015). La Figura 10A muestra datos de esmeraldas de los principales distritos en Colombia, y se observa que las muestras de La Pava se agrupan en el campo de esmeraldas ricas en Fe, un comportamiento que difiere sustancialmente de la tendencia general para esmeraldas colombianas. La Figura 10B muestra una comparación con datos de las principales minas del distrito de Muzo-Quípama; aunque la dispersión es amplia, el rango ocupado por las muestras de La Pava muestra una gran similitud con datos de las minas Futuro y Túnel 5.
Muestra | % Óxidos (peso) | Trazas (ppm) | |||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
BeO* | SiO₂ | Al₂O₃ | Na₂O | K₂O | CaO | SO₃ | Total | Cl | Sc | Ti | V | Cr | Fe | Mn | Ni | Cu | Zn | Ga | |||
EP-01 | 13,15 | 69,49 | 16,70 | ND | 0,12 | 0,28 | ND | 99,74 | 1349 | ND | 64 | 201 | 221 | 525 | ND | ND | ND | 32 | 16 | ||
EP-02 | 13,15 | 70,42 | 16,03 | ND | 0,07 | 0,03 | 0,07 | 99,77 | 932 | ND | ND | 297 | 145 | 589 | ND | ND | 20 | ND | 26 | ||
EP-03 | 13,15 | 69,71 | 16,08 | ND | 0,20 | 0,21 | 0,11 | 99,45 | 3030 | ND | 114 | 496 | 327 | 968 | ND | ND | 37 | 54 | 31 | ||
EP-04 | 13,15 | 71,53 | 14,94 | ND | 0,08 | 0,02 | ND | 99,73 | 1086 | 123 | ND | 407 | 292 | 517 | 2 | ND | 15 | ND | 28 | ||
EP-05 | 13,15 | 70,14 | 16,25 | ND | 0,09 | 0,06 | 0,03 | 99,73 | 1169 | ND | ND | 429 | 291 | 563 | ND | 8 | 17 | ND | 24 | ||
EP-06 | 13,15 | 70,73 | 15,79 | ND | 0,06 | 0,03 | ND | 99,76 | 971 | 92 | ND | 343 | 210 | 514 | ND | ND | ND | 21 | 26 | ||
EP-07 | 13,15 | 70,59 | 15,72 | ND | 0,07 | 0,10 | 0,07 | 99,70 | 1115 | 41 | ND | 439 | 397 | 646 | ND | 12 | 19 | ND | 28 | ||
EP-08 | 13,15 | 70,48 | 15,89 | ND | 0,07 | 0,16 | ND | 99,75 | 1249 | 113 | 32 | 172 | 158 | 410 | ND | ND | ND | 146 | 18 | ||
EP-09 | 13,15 | 69,32 | 17,12 | ND | 0,07 | 0,06 | ND | 99,71 | 1054 | ND | ND | 455 | 305 | 705 | 13 | ND | 20 | 41 | 26 | ||
EP-10 | 13,15 | 65,10 | 16,15 | 5,24 | 0,05 | 0,03 | 0,05 | 99,78 | 744 | 108 | ND | 322 | 170 | 562 | ND | ND | ND | 25 | 25 | ||
EP-11 | 13,15 | 70,56 | 15,98 | ND | 0,06 | 0,02 | ND | 99,76 | 785 | 94 | ND | 414 | 303 | 510 | ND | ND | 18 | 22 | 25 | ||
EP-12 | 13,15 | 70,20 | 16,16 | ND | 0,08 | 0,07 | 0,06 | 99,72 | 1161 | ND | ND | 396 | 273 | 623 | 1,2 | ND | 21 | ND | 25 | ||
EP-13 | 13,15 | 64,89 | 15,59 | 6,07 | 0,07 | ND | ND | 99,77 | 1011 | ND | ND | 337 | 218 | 489 | ND | ND | 17 | ND | 23 | ||
EP-14 | 13,15 | 70,65 | 15,74 | ND | 0,10 | 0,04 | ND | 99,69 | 1387 | ND | ND | 433 | 303 | 649 | 4,6 | ND | ND | 32 | 29 | ||
EP-15 | 13,15 | 69,85 | 15,44 | ND | 0,12 | 0,89 | 0,10 | 99,54 | 1534 | ND | 49 | 439 | 291 | 1269 | 54 | ND | 19 | 368 | 32 | ||
EP-16 | 13,15 | 70,23 | 16,20 | ND | 0,08 | 0,05 | ND | 99,71 | 1096 | ND | ND | 447 | 325 | 678 | ND | ND | ND | ND | 35 | ||
EP-17 | 13,15 | 70,06 | 16,14 | ND | 0,13 | 0,10 | 0,05 | 99,64 | 1646 | 70 | 50 | 452 | 302 | 662 | 17 | ND | 18 | 116 | 31 | ||
EP-18 | 13,15 | 70,45 | 15,93 | ND | 0,11 | 0,04 | 0,05 | 99,72 | 1275 | 94 | ND | 298 | 193 | 525 | ND | 26 | 147 | ND | 29 | ||
EP-19 | 13,15 | 66,62 | 15,41 | 4,49 | 0,07 | 0,04 | ND | 99,79 | 862 | ND | ND | 296 | 157 | 537 | ND | ND | 17 | ND | 27 | ||
EP-20 | 13,15 | 71,16 | 15,32 | ND | 0,08 | 0,03 | ND | 99,75 | 1082 | ND | ND | 342 | 171 | 609 | ND | ND | 23 | 33 | 32 | ||
EP-21 | 13,15 | 70,23 | 16,09 | ND | 0,11 | 0,09 | ND | 99,67 | 1484 | 107 | ND | 471 | 307 | 641 | 8 | ND | 17 | ND | 29 | ||
EP-22 | 13,15 | 71,17 | 15,22 | ND | 0,08 | 0,05 | 0,05 | 99,72 | 1038 | 57 | ND | 366 | 315 | 695 | ND | ND | 16 | ND | 26 | ||
EP-23 | 13,15 | 70,58 | 15,75 | ND | 0,09 | 0,07 | 0,07 | 99,70 | 1472 | ND | ND | 411 | 286 | 499 | ND | 11 | ND | 28 | 26 | ||
EP-24 | 13,15 | 71,03 | 15,36 | ND | 0,07 | 0,12 | ND | 99,73 | 911 | 142 | ND | 447 | 306 | 554 | ND | ND | 21 | 49 | 30 | ||
EP-25 | 13,15 | 70,22 | 15,83 | ND | 0,09 | 0,34 | 0,06 | 99,69 | 1105 | ND | ND | 470 | 347 | 786 | ND | ND | ND | 63 | 31 | ||
EP-26 | 13,15 | 69,99 | 16,49 | ND | 0,07 | 0,06 | ND | 99,75 | 1032 | ND | ND | 315 | 196 | 585 | ND | 15 | ND | 35 | 23 | ||
EP-27 | 13,15 | 70,59 | 15,80 | ND | 0,06 | 0,11 | ND | 99,71 | 1189 | ND | ND | 455 | 318 | 610 | ND | ND | 15 | 32 | 29 | ||
EP-28 | 13,15 | 70,63 | 15,89 | ND | 0,08 | 0,02 | ND | 99,78 | 906 | 82 | ND | 299 | 163 | 535 | ND | ND | ND | ND | 28 | ||
EP-29 | 13,15 | 66,96 | 15,82 | 3,71 | 0,08 | 0,04 | ND | 99,75 | 936 | ND | ND | 395 | 257 | 599 | ND | ND | 13 | ND | 24 | ||
EP-30 | 13,15 | 70,08 | 16,36 | ND | 0,06 | 0,04 | 0,07 | 99,75 | 766 | 86 | ND | 361 | 289 | 652 | ND | ND | ND | ND | 25 | ||
EP-31 | 13,15 | 69,64 | 16,85 | ND | 0,08 | ND | ND | 99,72 | 1042 | 111 | ND | 466 | 266 | 545 | 7 | 16 | 18 | 93 | 33 | ||
EP-32 | 13,15 | 71,08 | 15,31 | ND | 0,09 | 0,09 | ND | 99,72 | 1164 | ND | ND | 431 | 290 | 586 | ND | ND | ND | 45 | 29 | ||
EP-33 | 13,15 | 68,89 | 17,59 | ND | 0,06 | 0,13 | ND | 99,81 | 831 | ND | ND | 191 | 223 | 426 | ND | ND | 23 | ND | 16 | ||
EP-34 | 13,15 | 69,27 | 17,04 | ND | 0,06 | 0,24 | ND | 99,75 | 792 | 204 | ND | 369 | 285 | 585 | 13 | ND | ND | ND | 24 | ||
EP-35 | 13,15 | 70,54 | 16,02 | ND | 0,05 | ND | ND | 99,76 | 597 | 115 | ND | 478 | 327 | 624 | ND | ND | ND | ND | 30 | ||
EP-36 | 13,15 | 70,06 | 16,50 | ND | 0,05 | 0,03 | ND | 99,79 | 675 | 112 | ND | 255 | 148 | 578 | ND | 13 | ND | 62 | 28 | ||
EP-37 | 13,15 | 69,42 | 16,95 | ND | 0,11 | 0,04 | ND | 99,66 | 1776 | 79 | ND | 439 | 295 | 490 | ND | 13 | 30 | ND | 30 | ||
EP-38 | 13,15 | 70,06 | 16,43 | ND | 0,09 | ND | ND | 99,72 | 1093 | 81 | ND | 459 | 327 | 580 | ND | ND | ND | ND | 32 | ||
EP-39 | 13,15 | 69,60 | 16,42 | ND | 0,13 | 0,32 | 0,07 | 99,69 | 1515 | ND | ND | 435 | 310 | 585 | ND | ND | 21 | ND | 28 | ||
EP-40 | 13,15 | 69,33 | 17,10 | ND | 0,08 | 0,03 | ND | 99,69 | 1312 | 115 | ND | 417 | 295 | 596 | ND | ND | 29 | 35 | 24 | ||
EP-41 | 13,15 | 70,09 | 16,39 | ND | 0,05 | 0,07 | ND | 99,74 | 621 | 216 | ND | 452 | 301 | 645 | ND | ND | 19 | 41 | 24 | ||
EP-42 | 13,15 | 68,67 | 17,91 | ND | 0,05 | ND | ND | 99,78 | 878 | ND | ND | 384 | 264 | 483 | ND | ND | ND | ND | 21 | ||
EP-43 | 13,15 | 70,86 | 15,66 | ND | 0,06 | ND | ND | 99,74 | 1100 | ND | ND | 372 | 186 | 643 | ND | ND | ND | ND | 28 | ||
EP-44 | 13,15 | 70,93 | 15,62 | ND | 0,06 | ND | ND | 99,75 | 1090 | ND | ND | 366 | 209 | 542 | ND | ND | ND | ND | 31 | ||
EP-45 | 13,15 | 70,56 | 15,85 | ND | 0,08 | 0,07 | ND | 99,70 | 1388 | 108 | ND | 343 | 176 | 633 | ND | ND | ND | ND | 33 | ||
EP-46 | 13,15 | 71,04 | 15,34 | ND | 0,09 | 0,05 | ND | 99,67 | 1502 | 220 | ND | 370 | 226 | 666 | ND | ND | ND | ND | 31 | ||
EP-47 | 13,15 | 70,70 | 15,66 | ND | 0,09 | 0,09 | ND | 99,68 | 1354 | 234 | ND | 431 | 294 | 569 | ND | ND | 24 | ND | 27 | ||
EP-48 | 13,15 | 69,80 | 16,66 | ND | 0,08 | 0,03 | ND | 99,72 | 1197 | ND | ND | 417 | 302 | 557 | ND | 13 | 19 | 38 | 29 | ||
EP-49 | 13,15 | 69,76 | 16,46 | ND | 0,12 | 0,13 | ND | 99,63 | 2020 | ND | ND | 485 | 364 | 595 | ND | ND | ND | ND | 27 | ||
EP-50 | 13,15 | 70,47 | 15,73 | ND | 0,07 | 0,06 | 0,08 | 99,55 | 827 | 197 | ND | 813 | 1507 | 723 | ND | ND | 23 | 31 | 54 | ||
EP-51 | 13,15 | 67,10 | 16,32 | 3,09 | 0,07 | 0,03 | ND | 99,75 | 994 | ND | ND | 394 | 305 | 507 | ND | ND | 15 | ND | 27 | ||
EP-52 | 13,15 | 70,29 | 15,86 | ND | 0,09 | 0,28 | ND | 99,66 | 1373 | ND | ND | 413 | 304 | 788 | ND | 21 | 21 | 59 | 31 | ||
EP-53 | 13,15 | 70,91 | 15,51 | ND | 0,11 | 0,03 | ND | 99,70 | 1546 | 97 | ND | 373 | 287 | 432 | ND | ND | 18 | ND | 24 | ||
EP-54 | 13,15 | 69,07 | 17,27 | ND | 0,12 | 0,04 | ND | 99,65 | 1695 | ND | ND | 421 | 285 | 761 | ND | ND | ND | ND | 31 | ||
EP-55 | 13,15 | 69,02 | 17,23 | ND | 0,16 | 0,03 | ND | 99,59 | 2430 | ND | ND | 437 | 272 | 624 | ND | 15 | 23 | ND | 30 |
ND: no detectado. *Valor de referencia tomado de Groat et al. (2008).
Discusión
Control estructural
La mineralización de esmeraldas en La Pava muestra una asociación íntima con venas y brechas hidrotermales emplazadas durante un evento deformativo compresional, que ha generado estructuras como el Anticlinal de La Pava, un pliegue por propagación de falla descrito inicialmente por González-Durán et al. (2019). El papel de los anticlinales -particularmente aquellos vinculados a la propagación de fallas inversas- en la formación de yacimientos esmeraldíferos, ha sido documentado por diversos autores en áreas contiguas como Muzo y Coscuez (Branquet et al., 1999; Laumonier et al., 1999; Pignatelli et al., 2015). Existen distintas posturas sobre si existe una relación genética entre las fallas y la mineralización. Branquet et al. (1999) proponen que las fallas de cabalgamiento actúan como zonas de despegue donde se promueve la migración y focalización de fluidos mineralizantes. Bajo este modelo, la sobrepresión de fluidos, en zonas de pliegues por propagación de falla, puede conllevar al eventual brechamiento y mineralización; Cheilletz y Giuliani (1996) describen niveles de albititas brechadas y estratiformes similares a los observados en el río Itoco, y los interpretan como el resultado de la acción de brechamiento tectónico que ha favorecido el metasomatismo. Por otro lado, diversos autores plantean que no existe tal relación y que la mineralización es controlada por un arreglo de estructuras mesoscópicas derivadas de un evento de deformación frágil y penetrativo (Mora et al., 2005; Terraza, 2019).
No obstante, las características descritas en La Pava favorecen notablemente la intervención de una falla de cabalgamiento en la génesis de la mineralización. Dado que la falla del Águila contribuyó al desarrollo del Anticlinal de La Pava, una estructura que originó deslizamientos flexurales y generación de espacios que eventualmente fueron ocupados por la mineralización en el bloque colgante; esta falla también fue determinante para la posterior exhumación de la Formación Muzo en el área de La Pava. Así mismo, el Anticlinal de la Pava exhibe rasgos de un pliegue disarmónico o asimétrico; una característica descrita por Cheilletz y Giuliani (1996) en estructuras del distrito de Muzo-Quípama. Ante esta configuración, las zonas productivas tenderían a estar confinadas hacia la charnela de los pliegues, y es allí donde las labores de exploración deberían enfocarse.
Implicaciones estratigráficas
Desde el punto de vista estratigráfico, el horizonte de albititas brechadas constituye un marcador litogeoquímico indiscutible, que por su naturaleza puede ser fácilmente identificado en campañas de perforación, no solamente por su color sino por sus contrastantes características geoquímicas. Además, el nivel de albititas brechado y estratiforme es una característica litoestratigráfica indicativa de la Formación Muzo (Mantilla et al., 2007; Pignatelli et al., 2015; Reyes et al., 2006).
Al comparar los datos geoquímicos de la sección del río Itoco (Formación Muzo y Capotes) con valores disponibles en la literatura de unidades correlacionables, específicamente la Formación Paja, no se observan diferencias significativas en los datos generales. Particularmente, Campos y Roser (2007) reportan enriquecimientos en V, Cr y Ni, y empobrecimientos notables en CaO y Na₂O - tendencias que no sobresalen en el río Itoco-; los valores promedio de elementos involucrados en la mineralización (Be: 2,42 ppm (n=18), Y: 33,5 ppm (n=24), P₂O₅: 0,25% (n=24) y la relación Cr/V de 0,35) tampoco exhiben diferencias sustanciales.
Por otra parte, Rivera et al. (2018) reportan valores promedio de algunos elementos traza, entre los cuales se destacan enriquecimientos en Zn, Ni, Cd y una relación Cr/V de 0,43; además, usaron elementos sensibles a cambios redox para inferir condiciones paleoambientales. Al aplicar estos mismos índices en datos del río Itoco, se obtienen rangos similares, lo cual indicaría que no existieron variaciones composicionales ni de quimio-facies significativas durante la diagénesis de estas rocas. Cabe resaltar que la relación Cr/V es anómalamente baja en las rocas de la Formación Paja, en comparación con otras unidades circundantes (Rivera et al., 2018).
Asociaciones minerales
Por otra parte, las asociaciones paragenéticas muestran grandes similitudes con las descripciones publicadas sobre otras minas de esmeraldas en la región (Banks et al., 2000; Cheilletz et al., 1994; Giuliani et al., 2000; Mantilla et al., 2007; Ottaway et al., 1994; Romero, 1998). Una diferencia notable es la presencia de xenotimo en el ensamblaje de albitización; este mineral también ha sido reportado por Romero et al. (manuscrito en revisión), en la mina Cunas y allí también se encuentra asociado a una etapa de albitización. Distintos minerales portadores de elementos de tierras raras (REE) han sido descritos en zonas esmeraldíferas, entre ellos parisita-(Ce), synchysita-(Ce), bastnasita y monacita (Banks et al., 2000; Cheilletz y Giuliani, 1996; Pignatelli et al., 2015; Romero y Hernández, 1999; Romero, 1998; Toloza et al., 2019; Verbel, 2017); así como de minerales portadores de flúor -fluorita, fluorapatito y parisita-(Ce)- (Banks et al., 2000; Mantilla et al., 2007; Romero, 1998). El papel del flúor es crítico dado que actúa como ligando del berilio (Banks et al., 2000; Verbel, 2017), y promueve la lixiviación, transporte y precipitación de la esmeralda. Además del Be, también se ha planteado que los complejos fluorocarbonados tienen la capacidad de transportar REE y dar origen a la precipitación de fluorocarbonatos como la parisita- (Ce) (Verbel, 2017; Williams-Jones et al., 2012). En La Pava, la asociación entre albita, fluorapatito y xenotimo podría sugerir el involucramiento del flúor como ligando durante la albitización. De igual manera, durante la etapa de carbonatización, la presencia de minerales ricos en flúor como fluorapatito, sumado a la alta concentración de carbonatos; podrían indicar una actividad considerable de fluorucarbonatos que actuaron como ligandos, y habrían permitido el enriquecimiento necesario de Be para precipitar esmeraldas. Aunque en este estudio no se identificaron los minerales parisita-(Ce) y fluorita, no se descarta su ocurrencia en La Pava.
Por otro lado, la presencia de monacita autigénica podría indicar que la fuente de REE son las rocas encajantes y que no se requiere de aportes de componentes externos para dar origen a la mineralización; una situación similar se podría plantear para el Be, con base en los valores encontrados en las rocas del río Itoco.
Metasomatismo
El balance de masas realizado con los tres grupos de muestras sugiere que los procesos hidrotermales tuvieron la capacidad de movilizar la mayoría de los elementos analizados (Figura 8 y 9). Las variaciones más notorias, que se aprecian particularmente en las muestras de los túneles, están relacionadas con enriquecimientos en Mg, Ca, Zn, Cd, Mn, Y, U y empobrecimientos en Si, Al, Ti, Ba y Ga. Los altos valores de Mg, Ca y Mn podrían corresponder a la introducción de carbonatos (especialmente dolomita) durante la etapa de carbonatización. El Zn y Cd posiblemente estén alojados en esfaleritas y el Y en xenotimo. El U posiblemente este alojado en el xenotimo o en un mineral no identificado. Por otra parte, el empobrecimiento de Al, Ga y Si podría corresponder a la lixiviación que ocurre durante la carbonatización, específicamente cuando los carbonatos disuelven y reemplazan la albita. Del mismo modo, la deficiencia en Ti y Ba podría explicarse por la lixiviación, removilización y posterior precipitación de rutilo y barita (Figura 5H). Cabe resaltar que en las muestras del río Itoco (Figura 8), los empobrecimientos de Mg, Mn e Y, podrían sugerir un leve grado de metasomatismo, y que esas rocas no representan la composición exacta y original como se presume.
En la sección del río Itoco, donde el metasomatismo fue pobremente desarrollado, el contenido promedio de Be es de 3,5 ppm; mientras que en la quebrada La Pava es de 2,2 ppm y en los túneles de 1,6 ppm -valores muy similares a los reportados por Giuliani et al. (1999) y Pignatelli et al. (2017) para rocas de diferentes distritos esmeraldíferos colombianos-. Esto indica un porcentaje de lixiviación que oscila entre 37%- 55% en peso para el Be en zonas de intensa alteración hidrotermal. Giuliani et al. (1999) reportan para la zona productiva de Chivor, un porcentaje de lixiviación de Be del ≈18% en peso; mientras que Cheilletz y Royant (2001) reportan un valor de hasta 80%; lo cual indica que existen diferencias sustanciales en el grado de lixiviación en los distritos esmeraldíferos, un hecho posiblemente derivado de distintos regímenes de interacción fluido-roca. Además, este parámetro puede resultar útil al momento de realizar estimación de recursos, pues al integrarse con información volumétrica, podría aportar luces sobre la cantidad de Be que precipitaría como esmeralda. Un factor que debe esclarecerse es el porcentaje de Be efectivo que puede llegar a ser precipitado; no todo el Be lixiviado formará esmeraldas, ya que pueden surgir pérdidas en el transporte, así como por la captación por otros minerales como filosilicatos (Cheilletz y Royant, 2001).
Maya et al. (2004) realizan un compendio sobre el estado del arte de la exploración geoquímica en las regiones esmeraldíferas. Se contemplan diferentes escalas y métodos de prospección; para el caso de muestreo de rocas semidetallado, se contemplan los trabajos de Beus y Mineev (1972), PNUD (1975), Beus (1979), Escovar (1979), Mendoza-Parada (1996) y Giuliani et al. (1999); en donde se realizaron balances de masa comparando zonas con alteración hidrotermal vs zonas no alteradas. Los principales indicadores de zonas de lixiviación (considerados como anomalías negativas) son el Li, Cs, Be, K, Si, Al, Ti, P, REE, Rb, Sc, V y Cr. Por otra parte, las zonas mineralizadas tienden a presentar concentraciones altas de Na, Ca, Mg, Mn, S, y localmente Y, F y P. En los datos de La Pava, el comportamiento del Li, Cs, Be, Si, Al y Ti es consistente con un patrón de lixiviación muy definido; además el Zr, Nb y Ga también exhiben un comportamiento similar, aunque no hayan sido descritos en otros sectores. Por lo tanto, estos elementos podrían ser empleados como anomalías negativas en la prospección local. Elementos que han mostrado buenos resultados en otras áreas como el K, Ba, Rb y Sc, no son los mejores indicadores en los datos obtenidos de La Pava. En cuanto al enriquecimiento, los valores de P₂O₅, Y y U, son los indicadores más consistentes tanto en La Pava como en los túneles; a pesar de que elementos como el Ca, Mg, Mn, Cd y Zn exhiben la misma tendencia e incluso alcanzan picos elevados, solamente sobresalen a escala muy detallada y por lo tanto deberían tratarse con cautela. Para el caso del Y, debido a su estrecha asociación con el ensamblaje de albitización (a manera de xenotimo); podría constituir un marcador geoquímico fácilmente reconocible mediante análisis petrográficos. Una situación similar ocurre con el fluorapatito; el cual exhibe una relación cercana con las zonas productivas y es un mineral de fácil identificación incluso en campo. Los altos valores de P₂O₅ pueden corresponder con la presencia de fluorapatito o xenotimo; además, debido a la similitud entre el comportamiento geoquímico del P y el F (sumado a que el apatito alberga F en su estructura), se infiere que este último elemento podría también convertirse en un indicador geoquímico de zonas mineralizadas a escala local, algo que ha sido propuesto por algunos autores (Beus y Mineev, 1972; Verbel, 2017).
Otros indicadores geoquímicos como la relación K/Na, empobrecimientos en Mo, Pb, Cr y V; altas concentraciones de Na y Ca, no representan parámetros claros para la prospección local en La Pava. A nivel regional algunos elementos (Na, Sc y Be) exhiben enriquecimientos anómalos que permiten demarcar los distritos esmeraldíferos (Guerrero et al., 2020). Del mismo modo, la relación K/Na es un indicador útil para escalas más regionales y en otro tipo de muestras como suelos y aguas (Beus, 1979; Mendoza-Parada, 1996); entre tanto que los valores de Cr y V exhiben una distribución muy aleatoria debido a su fuerte lixiviación y posible incorporación en micas (Cheilletz y Royant, 2001; Giuliani et al., 2000).
Según Giuliani et al. (1999) el empobrecimiento en elementos como Be, Cr, V, P₂O₅, Ti, y REE en rocas aledañas a la mineralización, sería el resultado de la distribución parcial dentro del relleno de las venas; donde se encontrarían conformando fases minerales (esmeralda, fluorapatito, rutilo y parisita-(Ce)) o incorporándose como trazas. Este hecho respalda el uso de anomalías negativas y asociaciones minerales en las venas, como criterios de exploración de esmeraldas.
Composición química de las esmeraldas
Al comparar las concentraciones promedio de Fe de las esmeraldas de La Pava -608 ppm-, con el rango de los distritos colombianos -772 ppm a 3086 ppm- (Cedeño et al., 2015), se evidencia que las muestras de la Pava están por debajo del promedio, es decir, son pobres en Fe. De igual manera sucede con el Cr y el V; para La Pava los valores promedio que se calcularon son 288 ppm de Cr y 389 ppm de V; mientras que los promedios de esmeraldas colombianas oscilan entre 716 ppm a 3922 ppm para V y 1244 ppm a 5012 ppm para Cr. Al observar a primera vista las proporciones relativas de la Figura 10 puede parecer que las esmeraldas de la Pava sean ricas en Fe; sin embargo, las comparaciones realizadas permiten esclarecer que estas muestras son pobres en elementos cromóforos, en particular Cr y V. Como consecuencia, su distribución en el diagrama ternario se localiza en la esquina dominada por Fe; un comportamiento insólito para esmeraldas colombianas, pero recurrente en algunas minas del distrito Muzo-Quípama.
El bajo contenido de elementos cromóforos en las esmeraldas de La Pava puede ser el resultado de: i) bajo contenido de estos elementos en las rocas fuente; ii) cambios significativos en las variables fisicoquímicas (temperatura, presión y composición del fluido mineralizante) que no permitieron su movilización; iii) incorporación en otros minerales. Al comparar las concentraciones de Cr y V del río Itoco con datos publicados en otras zonas esmeraldíferas (Beus y Mineev, 1972; Cheilletz y Giuliani, 1996; Pignatelli et al., 2017), no se observan diferencias significativas; una situación similar se observa al revisar los datos de presión y temperatura reportados para esmeraldas de La Pava (Cheilletz y Giuliani, 1996; Giuliani et al., 2000; Jiménez, 2020; Ramírez, 2020; Toloza, 2017). Además, las fluctuaciones en las concentraciones de estos elementos son altas en las muestras de La Pava y los túneles, indicando un alto grado de lixiviación, lo cual excluye la posibilidad de una baja oferta o limitada movilidad durante el metasomatismo. Por otra parte, las micas -particularmente la moscovita- son minerales capaces de captar cantidades sustanciales de Cr y V, alcanzando valores entre 0,1% a 6% e incluso mayores (Imai et al., 1980; Morata et al., 2001; Randive et al., 2015), los cuales son claramente mayores a las concentraciones que pueden albergar las esmeraldas (< 2%) (Groat et al., 2008). En la secuencia paragenética de La Pava, se observan múltiples episodios de cristalización de moscovita, previos a la depositación de la esmeralda (Figura 6), por lo tanto, podría plantearse que los elementos cromóforos fueron en gran manera sustraídos por la moscovita, de tal manera que, al momento de la cristalización de la esmeralda, su disponibilidad era limitada. Cheilletz et al. (1994) soporta este hecho dado que reporta concentraciones de Cr: 208 ppm y V: 892 ppm en micas del distrito Muzo-Quípama.
La relación promedio de Cr/V en las esmeraldas (0,71) difiere notablemente con respecto a la de las rocas del río Itoco (0,51) y las zonas mineralizadas (0,22-0,36), hecho que sugiere que hay una mayor captación de Cr en las esmeraldas en comparación con la disponibilidad en las rocas. Posiblemente, esta relación esté influenciada por una mayor captación de V en las micas y por lo tanto un enriquecimiento relativo en el fluido mineralizante de Cr. La relación Cr/V en los análisis reportados por Cheilletz et al. (1994) es de 0,23, un valor muy cercano al que se halla en las zonas mineralizadas de La Pava. Este comportamiento también se ha observado en el Cinturón Esmeraldífero Oriental, donde las relaciones Cr/V son normalmente mayores a 1 en las esmeraldas, pero menores en las rocas reservorio (Cedeño et al., 2015; Pignatelli et al., 2017).
Conclusiones
En la mina La Pava existen una serie de parámetros indicadores de la mineralización de esmeraldas a diferentes escalas. A escala semidetallada: pliegues tipo anticlinal por propagación de fallas inversas; niveles estratigráficos dominados por pelitas, que contengan horizontes albitizados con anomalías geoquímicas que indiquen un alto grado de lixiviación y removilización de una amplia gama de elementos. El metasomatismo de Na y Ca (albitización y carbonatización) es un rasgo típico de bloques mineralizados. A escala local, las anomalías geoquímicas positivas de Mg, P, Cd, Zn, Mn, e Y; negativas de Si, Al, Ti, K, Zr, Nb, Li, Cs, Be, Ba y Ga; en conjunto con la localización de zonas de charnela con presencia de venas, cuyas asociaciones paragenéticas incluyan albita, fluorapatito, xenotimo, calcita, dolomita, rutilo y otras fases ricas en F y REE como parisita-(Ce). Otros minerales como barita, sulfuros y filosilicatos, pueden ocurrir como accesorios y deben tratarse con precaución antes de considerarse como indicadores.
Las esmeraldas de La Pava son pobres en V y Cr, con concentraciones de Fe cercanas a las de otras esmeraldas colombianas. La relación Cr/V no refleja fielmente la proporción original de estos elementos en las rocas encajantes, lo cual puede ser el resultado de una incorporación diferencial en las micas que modifica la disponibilidad de Cr y V al momento de la precipitación de las esmeraldas.
Los parámetros analizados (geoquímica, mineralogía y control estructural) en La Pava exhiben diferencias sustanciales con datos reportados en otras minas y distritos, lo cual demuestra la variabilidad y complejidad de los yacimientos esmeraldíferos. Esto no implica que no se pueda realizar extrapolación de variables, sino más bien que se debe estudiar cada caso en particular para poder establecer los parámetros comunes y así lograr a generar un modelo.