Boletín Geológico y Minero 135 (4)
ISSN-L: 0366-0176, eISSN: 2253-6167
https://doi.org/10.21701/bolgeomin/135.4/001

Determinación de origen geográfico en esmeraldas de Colombia, Zambia y Pakistán: análisis de fotoluminiscencia y correlación con elementos traza

Geographic origin determination in emeralds from Colombia, Zambia, and Pakistan: photoluminescence analysis and correlation with trace elements

 

1. Introducción

 

La esmeralda es un ciclosilicato que cristaliza en el sistema hexagonal, es la variedad verde de los berilos y tiene una fórmula general Be3Al2(SiO3)6. Esta debe su color a la sustitución del Al, que se encuentra en posición octaédrica en su estructura cristalina, por elementos cromóforos como Cr, V y Fe (Araújo et al., 2019Araújo Neto, J. F. de, Müller, A., Barreto, S. de B., Santos, L. C. M. de L., and Carrino, T. A. (2019). Mineralogical and gemological characterization of emerald crystals from Paraná deposit, NE Brazil: a study of mineral chemistry, absorption and reflectance spectroscopy and thermal analysis. Brazilian Journal of Geology, 49(3), e20190014. 10.1590/2317-4889201920190014.). Las esmeraldas se han formado en diferentes contextos geológicos y edades en la Tierra. El más conocido se da a partir de procesos magmáticos, como es el caso de Zambia, Brasil o Rusia, entre otros, (Giuliani et al., 1990Giuliani, G., Silva, L. J. H. D., & Couto, P. (1990). Origin of emerald deposits of Brazil. Mineralium Deposita, 25, 57-64. 10.1007/BF03326384.; Giuliani & Groat, 2019Giuliani, G., & Groat, L. A. (2019). Geology of corundum and emerald gem deposits: a review. Gems and Gemology, 55(4), 464-489. 10.5741/GEMS.55.4.464.; Alonso-Perez & Day, 2021Alonso-Perez, R., & Day, J. M. D. (2021). Rare Earth Element and Incompatible Trace Element Abundances in Emeralds Reveal Their Formation Environments. Minerals, 11(5), 513. 10.3390/min11050513.). Otras se forman por procesos metamórficos, como las esmeraldas provenientes de Pakistán o Austria, entre otros (Giuliani & Groat, 2019Giuliani, G., & Groat, L. A. (2019). Geology of corundum and emerald gem deposits: a review. Gems and Gemology, 55(4), 464-489. 10.5741/GEMS.55.4.464.; Alonso-Perez & Day, 2021Alonso-Perez, R., & Day, J. M. D. (2021). Rare Earth Element and Incompatible Trace Element Abundances in Emeralds Reveal Their Formation Environments. Minerals, 11(5), 513. 10.3390/min11050513.); y, finalmente, las esmeraldas que se encuentran hospedadas en rocas sedimentarias, como es el caso de las esmeraldas colombianas (Giuliani et al.,1993Giuliani, G., Cheilletz, A., Sheppard, S. M. F., & Arboleda, C. (1993). Geochemistry and origin of the emerald deposits of Colombia. En Fenoll Hach-Ali, Torres-Ruiz y Gervilla (Eds.), Current research in geology applied to ore deposits (pp. 105-108). Biennial SGA Meeting, 2, Granada (España), 1993/09/09-11. ISBN 84-338-1772-8., 1995Giuliani, G., Cheilletz, A., Arboleda, C., Carrillo, V., Rueda, F., & Baker, J. H. (1995). An evaporitic origin of the parent brines of Colombian emeralds: fluid inclusion and sulphur isotope evidence. European Journal of Mineralogy, 7(1), 151-165.; Banks et al.,2000Banks, D. A., Giuliani, G., Yardley, B. W. D., & Cheilletz, A. (2000). Emerald mineralisation in Colombia: Fluid chemistry and the role of brine mixing. Mineralium Deposita, 35, 699-713. 10.1007/s001260050273.; Giuliani & Groat, 2019Giuliani, G., & Groat, L. A. (2019). Geology of corundum and emerald gem deposits: a review. Gems and Gemology, 55(4), 464-489. 10.5741/GEMS.55.4.464.; Alonso-Perez & Day, 2021Alonso-Perez, R., & Day, J. M. D. (2021). Rare Earth Element and Incompatible Trace Element Abundances in Emeralds Reveal Their Formation Environments. Minerals, 11(5), 513. 10.3390/min11050513.). La determinación del origen de una esmeralda es un desafío para los laboratorios gemológicos en el mundo, debido al hallazgo de nuevos depósitos que pueden llegar a presentar características similares a esmeraldas de diferentes países, como es el caso de Musakashi, en Zambia, o Pansjhir, en Afganistán (Saeseaw et al., 2014SaeseawS., PardieuV., & SangsawongS. (2014). Three-Phase Inclusions in Emerald and Their Impact on Origin Determination. Gems & Gemology, 50(2), 114-132. 10.5741/GEMS.50.2.114.), las cuales presentan características fisicoquímicas similares a las esmeraldas colombianas. Por lo anterior, se han venido desarrollando diferentes métodos para determinar el origen de las esmeraldas a partir de su composición química. Trabajos como los de Abduriyim & Kitawaki (2006Abduriyim, A., & Kitawaki, H. (2006). Applications of laser ablation – inductively coupled plasma – mass spectrometry (LA-ICP-MS) to gemology. Gems and Gemology, 42(2), 98-118.), Bersani et al. (2014Bersani, D., Azzi, G., Lambruschi, E., Barone, G., Mazzoleni, P., Raneri, S., & Lottici, P. P. (2014). Characterization of emeralds by micro-Raman spectroscopy. Journal of Raman Spectroscopy, 45, 11-12. 10.1002/jrs.4524.), Kiefert & Karampelas (2011Kiefert, L., & Karampelas, S. (2011). Use of the Raman spectrometer in gemmological laboratories: Review. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 80, 119-124. 10.1016/j.saa.2011.03.004.), Moroz et al. (2000Moroz, I., Roth, M., Boudeulle, M., & Panczer, G. (2000). Raman Microspectroscopy and Flourescence of Emeralds from Various Deposits. Journal of Raman Spectroscopy, 31, 485-490.), Palke et al. (2019PalkeA., Saeseaw, S., RenfroN. D., SunZ., & McClure, F. (2019). Geographic origin determination of blue sapphire. Gems & Gemology, 55(4), 536-579. 10.5741/GEMS.55.4.536.), Karampelas et al. (2019Karampelas, S., Al-Alawi, A., Al-Shaybani, B., Mohamed, F., & Sangsawong, S. (2019). Emeralds from the most important occurrences: Chemical and spectroscopic data. Minerals, 9(9), 561. 10.3390/min9090561.) y Alonso-Perez et al. (2024Alonso-Perez, R., Day, J. M., Pearson, D. G., Luo, Y., Palacios, M. A.,Sudhakar, R., & Palke, A. (2024). Exploring emerald global geochemical provenance through fingerprinting and machine learning methods. Artificial Intelligence in Geosciences, 5, p. 100085. 10.1016/j.aiig.2024.100085.), entre otros, han logrado establecer algunos parámetros que permiten distinguir esmeraldas de Colombia de las de otros países, a partir de las concentraciones de Ga entre 5-50 partes por millón (ppm), Fe 150-3000 ppm, Na 2500-8000 ppm, Mg 3400-10000 ppm, Cs 5-30 ppm y también bajas concentraciones de litio y potasio. Cedeño et al. (2015Cedeño, C. J., Jiménez, J. F., Herreño, J. H., & Fortaleché, D. (2015). Progress on the study of parameters related to the origin of Colombian emeralds. Color, World Emer, (Special Issue), 88-97.) reportaron, en un estudio hecho en 530 esmeraldas de los principales distritos productores de Colombia, ubicados en ambos cinturones esmeraldíferos (occidental: CEOCC; y oriental: CEOR), una relación Cr/V de la siguiente manera: para el CEOR en el distrito de Chivor, se observaron valores de Cr/V entre 4 y 5 y en Gachalá de 2.5; mientras que, para el CEOCC en el distrito de Muzo, los valores de la relación son muy cercanos a 1.

Adicionalmente, Giuliani et al. (2019Giuliani, G., Groat, L. A., Marshall, D., Fallick, A. E., & Branquet, Y. (2019). Emerald deposits: A review and enhanced classification. Minerals, 9(2), 105. 10.3390/min9020105.) reportan datos con microsonda electrónica (EPMA) de diferentes localidades en el mundo, y resaltan que la mayoría de los lugares tienen una mayor concentración de Cr frente a V con las siguientes excepciones: Canadá, China, Pakistán, Noruega y Colombia (CEOCC). A grandes rasgos, una de las mayores características de la esmeralda colombiana, desde el punto de vista químico, y específicamente de la concentración de elementos cromóforos, es la baja concentración de Fe, además de similares concentraciones de Cr y V, hecho que comparte con otras zonas como Noruega, China y Afganistán.

Por otra parte, la fotoluminiscencia está tomando una mayor relevancia, debido a su capacidad de brindar información sobre la estructura cristalina de los materiales gemológicos: en diamantes, se utiliza para identificación y distinción de mejoramiento de color o simplemente para distinguir entre naturales y sintéticos (Wotherspoon et al., 2003Wotherspoon, A., Steeds, J. W., Catmull, B., & Butler, J. (2003). Photoluminescence and positron annihilation measurements of nitrogen doped CVD diamond. Diamond and Related Materials, 12, 652-657. 10.1016/S0925-9635(02)00229-7.; Breeding et al., 2009Breeding, C. M., & Shigley, J. E. (2009). The “type” classification system of diamonds and its importance in gemology. Gems & Gemology, 45(2), 96-111. 10.5741/GEMS.45.2.96.; Lim et al., 2010Lim, H., Park, S., Cheong, H., Choi, H.-M., & Kim, Y. C. (2010). Discrimination between natural and HPHT-treated type IIa diamonds using photoluminescence spectroscopy. Diamond and Related Materials, 19(10), 1254-1258.). Además, la fotoluminiscencia en las esmeraldas se usa para la identificación y caracterización de una banda generada por el Cr (Wood, 1965Wood, D. L. (1965). Absorption, fluorescence, and Zeeman effect in Emerald. The Journal of Chemical Physics, 42(10), 3404-3410. 10.1063/1.1695742.) que contiene dos picos (líneas de absorción), conocidos como R1 y R2, en aproximadamente 683 y 680 nm respectivamente. Estos picos han sido relacionados con el origen geográfico de las esmeraldas en trabajos previos como el de Schmetzer (2015Schmetzer, K. (2015). Photoluminescence of emeralds: Sample orientation procedure and correlation of the R1 peak position with SiO2 contents. Journal of Gemmology, 34, (5), 441-443.), Thompson et al. (2014Thompson, D. B., Kidd, J. D., Åström, M., Scarani, A., & Smith, C. P. (2014). A Comparison of R-line Photoluminescence of Emeralds from Different Origins. The Journal of Gemmology, 34(4), 334-343. 10.15506/jog.2014.34.4.334., 2017Thompson, D. B., Bayens, C. J., Morgan, M. B., Myrick, T. J., & Sims, N. E. (2017). Photoluminescence spectra of emeralds from Colombia, Afghanistan, and Zambia. Gems & Gemology, 53(3), 296-311. 10.5741/GEMS.53.3.296.) y Moroz et al. (2000Moroz, I., Roth, M., Boudeulle, M., & Panczer, G. (2000). Raman Microspectroscopy and Flourescence of Emeralds from Various Deposits. Journal of Raman Spectroscopy, 31, 485-490.).

Algunos autores han correlacionado la química de las esmeraldas con espectros de fotoluminiscencia. Moroz et al. (2000Moroz, I., Roth, M., Boudeulle, M., & Panczer, G. (2000). Raman Microspectroscopy and Flourescence of Emeralds from Various Deposits. Journal of Raman Spectroscopy, 31, 485-490.) sugirieron que la sustitución de Al por Fe y Mg, la cual es común para la esmeralda en depósitos de tipo esquisto, cambia la posición de los picos en los espectros y, en consecuencia, podría usarse para diferenciar esmeraldas de diferentes ambientes geológicos. Thompson et al. (2014Thompson, D. B., Kidd, J. D., Åström, M., Scarani, A., & Smith, C. P. (2014). A Comparison of R-line Photoluminescence of Emeralds from Different Origins. The Journal of Gemmology, 34(4), 334-343. 10.15506/jog.2014.34.4.334.) reportaron que los cambios en el pico R1 son generados por la concentración y posición del Si en el cristal, donde pueden ocurrir sustituciones de impurezas o defectos cristalinos. Sin embargo, Schmetzer (2015Schmetzer, K. (2015). Photoluminescence of emeralds: Sample orientation procedure and correlation of the R1 peak position with SiO2 contents. Journal of Gemmology, 34, (5), 441-443.) propone que el cambio en la concentración de Si podría surgir de algunas sustituciones en el lugar del Al, lo que produciría la incorporación de otro elemento, generando los cambios en el pico R1. Posteriormente, Thompson et al. (2017Thompson, D. B., Bayens, C. J., Morgan, M. B., Myrick, T. J., & Sims, N. E. (2017). Photoluminescence spectra of emeralds from Colombia, Afghanistan, and Zambia. Gems & Gemology, 53(3), 296-311. 10.5741/GEMS.53.3.296.) concluyeron que existe una fuerte correlación entre altos contenidos de Mg y Na con longitudes de onda más largas en el espectro (para el pico R1) de las esmeraldas pertenecientes al tipo esquisto, mientras que el caso contrario ocurre para las esmeraldas con baja concentración de estos elementos que pertenecen a otro tipo de ambientes, como el tipo colombiano, con valores de longitud de onda más cortos.

2. Materiales y métodos

 

Se analizaron un total de 82 muestras de la siguiente manera: 59 muestras procedentes de Colombia, de las cuales 29 muestras pertenecen al cinturón esmeraldífero occidental (CEOCC), y 30 muestras al cinturón esmeraldífero oriental (CEOR); además, se incluyeron 10 muestras de Pakistán (PK) y 13 muestras de Zambia (Zmb). A través de la observación de las esmeraldas de las diferentes localidades, se describió la reacción a luz UV de onda larga con lámpara ultravioleta de cada una de las muestras (Tabla 1). Para este análisis se usaron las lámparas marca Convoy con luz en 365 nm y la lámpara Varsun R 838 con luz en 365 nm.

Tabla 1. Graduación de fluorescencia de las muestras de estudio. 
Table 1. Fluorescence grading of the studied samples. 
LOCALIDADReacción a luz UV (Fluorescencia 365 nm)
Rosa IntensoRosa medioRosa débilRojo IntensoRojo medioRojo débilInerte
Cinturón esmeraldífero occidental CEOCCEOMCATEDRAL_12X
EOMCATEDRAL_17XX
EOMCATEDRAL_18X
EOMCATEDRAL_19X
EOMCATEDRAL_20X
EOMCATEDRAL_21X
EOMCATEDRAL_22X
EOMCATEDRAL_23X
EOMCATEDRAL_24X
EOMCATEDRAL_25X
CUNAS_01_PL_600_PER_3X
CUNAS_02_PL_600_PERX
CUNAS_03_PL_600_PERX
CUNAS_04_PL_600_PERXX
CUNAS_06_PL_600_PERX
CUNAS_07_PL_600_PERX
CUNAS_18_PL_600_PERX
CUNAS_19_PL_600_PERXX
CUNAS_20_PL_600_PERX
CUNAS_21_PL_600_PERX
CUNAS_22_PL_600_PERX
EOPBCULATA_23X
EOPBCULATA_24X
EOPBCULATA_25X
EOPBCULATA_26X
EOPBCULATA_27X
EOPBCULATA_28X
EOPBCULATA_29X
EOPBCULATA_30X
Cinturón esmeraldífero oriental CEORTIGRES_2_PL_600X
TIGRES_7_PL_601X
TIGRES_8_PL_602X
TIGRES_9_PL_603X
TIGRES_11_PL_604X
TIGRES_CANU_01X
TIGRES_CANU_02X
TIGRES_CANU_03X
TIGRES_CANU_04X
TIGRES_CANU_05X
TIGRES_CANU_06X
ACHIOTE_T_GOCEN_01X
ACHIOTE_T_GOCEN_02X
ACHIOTE_T_GOCEN_03X
ACHIOTE_T_GOCEN_04X
ACHIOTE_T_GOCEN_05X
ACHIOTE_T_GOCEN_06X
ACHIOTE_T_GOCEN_07X
ACHIOTE_T_GOCEN_08X
ACHIOTE_T_GOCEN_09X
EOCHPALOARAÑADO_1X
EOCHPALOARAÑADO_2X
EOCHPALOARAÑADO_3X
EOCHPALOARAÑADO_4X
EOCHPALOARAÑADO_5X
EOCHPALOARAÑADO_6X
EOCHPALOARAÑADO_7X
EOCHPALOARAÑADO_8X
EOCHPALOARAÑADO_9X
EOCHPALOARAÑADO_10X
Pakistán (PK)PAKISTAN_SWAT_1X
PAKISTAN_SWAT_2X
PAKISTAN_SWAT_3X
PAKISTAN_SWAT_4X
PAKISTAN_SWAT_5X
PAKISTAN_SWAT_6X
PAKISTAN_SWAT_7X
PAKISTAN_SWAT_8X
PAKISTAN_SWAT_9X
PAKISTAN_MINGORA_18X
Zambia (ZMB)ZAMBIA_KAGEM_1_PERX
ZAMBIA_KAGEM_2_PERX
ZAMBIA_KAGEM_3_PERX
ZAMBIA_KAGEM_4_PERX
ZAMBIA_KAGEM_6_PERX
ZAMBIA_KAGEM_7_PERX
ZAMBIA_KAGEM_8_PERX
ZAMBIA_KAGEM_11_PERX
ZAMBIA_KAGEM_13_PERX
ZAMBIA_KAGEM_15_PERX
ZAMBIA_KAGEM_16_PERX
ZAMBIA_KAGEM_17_PERX
ZAMBIA_KAGEM_18_PERX

Los análisis de fotoluminiscencia (Tablas 2 y 3) fueron tomados con el equipo Horiba LabRAM HR Evolution equipment, con un tiempo de adquisición 5 s, acumulaciones de 10 scans, un rango de longitud de onda entre 660 y 700 nm, objetivo plan N 10x / 0.25 ∞ / - / FN22, grilla de difracción 1800 gr/mm, resolución 0.025 nm, láser 532 nm, poder del láser 250 mW, poder de excitación 0.5 mW y un diámetro del láser de 100 μm.

El procedimiento de toma de espectros por fotoluminiscencia se inició con la orientación de cada una de las muestras. En los cristales que tienen bien definida su forma de prisma hexagonal, resulta fácil identificar la dirección de toma del espectro, sea perpendicular o paralela al eje cristalográfico C (Fig. 1); sin embargo, las muestras en las que no se puede definir la orientación, presentan dificultad para determinar la dirección de sus ejes cristalográficos; por lo tanto, se procedió a usar como herramienta de apoyo, un polariscopio GIA modelo 415, instrumento que resulta útil para determinar la orientación cristalográfica de estas muestras.

media/001_001.jpg
Figura 1. Orientación cristalográfica de un cristal de esmeralda. Perpendicular al eje C (flecha roja), paralelo al eje C (flecha azul). Toma de espectro con láser verde (532 nm) orientado perpendicular al eje cristalográfico C (derecha). 
Figure 1. Crystallographic orientation of an emerald crystal. Perpendicular to the C axis (red arrow), parallel to the C axis (blue arrow). The laser (532 nm) is oriented perpendicular to the C axis (right). 

La esmeralda debe ser colocada con la superficie de la cara que coincida con la orientación del eje cristalográfico deseado (perpendicular o paralelo). El láser usado para este estudio corresponde a 532 nm, cuyo color es verde (Fig. 1, derecha) El tiempo de toma de cada espectro es de aproximadamente 3.20 min. En total se tomaron 160 espectros.

Posteriormente, a los espectros obtenidos por fotoluminiscencia se les realizó un pretratamiento en el programa Lab Spec6 HORIBA Scientific de la siguiente manera: corrección de la línea base con algoritmo usado de tipo polinomio de grado 7, máximo de puntos 203 y reducción de ruido de 4 puntos, tipo de suavizado y filtrado polinómico, parámetros del suavizado tamaño 31, grado 2. El análisis y detección de la posición de los picos (R), la intensidad (I), el área (A) y el ancho de altura medio (W) es llevado a cabo por medio de ecuaciones matemáticas usando distribución simétrica Gaussiana del programa Labspec 6 en puntos de mayor intensidad.

La determinación química de cada una de las muestras se realizó con el equipo de fluorescencia de rayos X de energía dispersiva (ED-XRF) PANalytical Epsilon 5, con las siguientes características: tubo de rayos X Sc/W, ánodo de W, potencia máxima 600 W, detector de compuestos de Ge, 30 mm2, temperatura del detector -196 °C, número de targets o filtros: 14; condiciones de trabajo de los filtros 25 a 100 kV/ 6 a 24 mA, ángulo de trabajo 45°, materiales de los filtros: Al2O3, CeO2, Ag, Mo, Zr, Ge, Fe, Ti, Al, B4C, Si, KBr, CaF2, W. Para el tratamiento estadístico se procedió de la siguiente manera: para los datos de ED-XRF se realizó un pretratamiento distribuyéndolos en elementos mayores, en unidades de porcentaje (Tabla 4) y los elementos traza en ppm (Tabla 5). Además, se compararon los datos de las esmeraldas de Pakistán y Zambia con estudios previos, con el fin de mejorar la representatividad de los datos presentados.

Tabla 4. Composición de elementos mayores por XRF. El contenido de BeO (13,15 %) fue añadido manualmente como valor de referencia (Groatet al., 2008Groat, L. A., Giuliani, G., Marshall, D. D., & Turner, D. (2008). Emerald deposits and occurrences: A review. Ore Geology Reviews, 34(1-2), 87-112. 10.1016/j.oregeorev.2007.09.003.), para completar el cierre elemental. LOD= Límite de detección. 
Table 4. Major element composition by XRF. The BeO content (13.15 %) was manually added as a reference value (Groatet al., 2008Groat, L. A., Giuliani, G., Marshall, D. D., & Turner, D. (2008). Emerald deposits and occurrences: A review. Ore Geology Reviews, 34(1-2), 87-112. 10.1016/j.oregeorev.2007.09.003.), to ensure elemental closure. LOD= Limit of detection. 
ORIGENMUESTRADISTRITOBeO(%)Na2O(%)MgO(%)Al2O3(%)SiO2(%)TOTAL
CEOCCEOMCATEDRAL_12Muzo - Q13.150< LOD< LOD15.37066.82095.340
CEOCCEOMCATEDRAL_17Muzo - Q13.150< LOD4.64715.74465.04998.590
CEOCCEOMCATEDRAL_18Muzo - Q13.150< LOD< LOD15.78270.09199.023
CEOCCEOMCATEDRAL_19Muzo - Q13.150< LOD< LOD16.89367.46997.512
CEOCCEOMCATEDRAL_20Muzo - Q13.150< LOD< LOD15.91569.19098.255
CEOCCEOMCATEDRAL_21Muzo - Q13.150< LOD4.94016.19864.52198.809
CEOCCEOMCATEDRAL_22Muzo - Q13.150< LOD< LOD17.59968.18998.938
CEOCCEOMCATEDRAL_23Muzo - Q13.150< LOD< LOD17.51268.33598.997
CEOCCEOMCATEDRAL_24Muzo - Q13.150< LOD< LOD15.92569.82098.895
CEOCCEOMCATEDRAL_25Muzo - Q13.150< LOD< LOD18.23667.49698.882
CEOCCEOMACUNAS_01Maripí13.1501.6002.85012.92068.56099.080
CEOCCEOMACUNAS_02Maripí13.150< LOD2.69012.68067.57096.090
CEOCCEOMACUNAS_03Maripí13.1500.3903.15014.41067.09098.190
CEOCCEOMACUNAS_04Maripí13.1502.1102.47012.64066.93097.300
CEOCCEOMACUNAS_06Maripí13.1503.1203.68012.94064.48097.370
CEOCCEOMACUNAS_07Maripí13.1500.8303.01012.99067.78097.760
CEOCCEOMACUNAS_Q_18Maripí13.1504.2000.05013.93066.53097.860
CEOCCEOMACUNAS_19Maripí13.150< LOD1.71012.23070.54097.630
CEOCCEOMACUNAS_20Maripí13.1501.2902.48012.71067.56097.190
CEOCCEOMACUNAS_21Maripí13.1500.6800.11013.99069.00096.930
CEOCCEOMACUNAS_22Maripí13.1500.6300.37013.67070.03097.850
CEOCCEOPBCULATA_23_2Peñas Blancas13.150< LOD< LOD20.27562.06995.494
CEOCCEOPBCULATA_24Peñas Blancas13.150< LOD< LOD16.82067.82397.793
CEOCCEOPBCULATA_25Peñas Blancas13.150< LOD< LOD22.32263.75499.226
CEOCCEOPBCULATA_26Peñas Blancas13.150< LOD< LOD20.28665.19998.635
CEOCCEOPBCULATA_27_2Peñas Blancas13.150< LOD< LOD21.88363.63598.668
CEOCCEOPBCULATA_28Peñas Blancas13.150< LOD< LOD22.24663.39698.792
CEOCCEOPBCULATA_29Peñas Blancas13.150< LOD< LOD25.25460.72999.133
CEOCCEOPBCULATA_30Peñas Blancas13.150< LOD< LOD21.79564.26399.208
CEOCCPromedio13.1501.6502.47416.59266.68297.782
CEOCCDesviación Estándar0.0001.2871.5723.5992.5752.177
CEOREOCHTIGRES_02Chivor13.150< LOD< LOD15.62669.77098.546
CEOREOCHTIGRES_07Chivor13.150< LOD< LOD15.65065.60194.401
CEOREOCHTIGRES_08Chivor13.150< LOD< LOD18.39767.50199.048
CEOREOCHTIGRES_09Chivor13.150< LOD7.17416.29862.27498.896
CEOREOCHTIGRES_11Chivor13.150< LOD2.81119.29763.39798.655
CEORTIGRES_01Chivor13.150< LOD0.33213.78468.52295.788
CEORTIGRES_02Chivor13.150< LOD7.12115.63262.10898.011
CEORTIGRES_03Chivor13.150< LOD2.01514.21266.39695.773
CEORTIGRES_04Chivor13.150< LOD6.00013.95756.48989.596
CEORTIGRES_05Chivor13.150< LOD3.45014.59966.66197.860
CEORTIGRES_06Chivor13.150< LOD1.57512.61670.55997.900
CEORACHIOTE_1Somondoco13.150< LOD3.58021.73061.11099.570
CEORACHIOTE_2Somondoco13.150< LOD< LOD22.48063.89099.520
CEORACHIOTE_3Somondoco13.150< LOD< LOD16.67069.04098.860
CEORACHIOTE_4Somondoco13.150< LOD< LOD15.86069.03098.040
CEORACHIOTE_5Somondoco13.150< LOD< LOD21.02063.81097.980
CEORACHIOTE_6Somondoco13.1506.270< LOD15.28061.68096.380
CEORACHIOTE_7Somondoco13.150< LOD< LOD17.82067.39098.360
CEORACHIOTE_8Somondoco13.15014.060< LOD18.56053.01098.780
CEORACHIOTE_9Somondoco13.150< LOD< LOD18.52067.51099.180
CEOREOCHPALOARAÑA_01Macanal13.150< LOD0.18114.45771.49999.287
CEOREOCHPALOARAÑA_02Macanal13.150< LOD2.55415.39567.81798.916
CEOREOCHPALOARAÑA_03Macanal13.150< LOD< LOD14.63671.68599.471
CEOREOCHPALOARAÑA_04Macanal13.1502.0802.79314.73066.41599.168
CEOREOCHPALOARAÑA_05Macanal13.1503.2080.56415.68366.89699.501
CEOREOCHPALOARAÑA_06Macanal13.1501.3482.95816.10465.79499.354
CEOREOCHPALOARAÑA_07Macanal13.1500.321< LOD15.65169.74598.867
CEOREOCHPALOARAÑA_08Macanal13.1502.6242.86914.87765.55599.075
CEOREOCHPALOARAÑA_09Macanal13.1503.2262.39214.57565.80599.148
CEOREOCHPALOARAÑA_10Macanal13.150< LOD< LOD16.90969.03699.095
CEORPromedio13.1504.1423.02316.36865.86798.101
CEORDesviación Estándar0.0004.3682.1392.3964.1672.039
PKESM-PAKISTÁN_SWAT_01Swat13.150< LOD1.71113.61968.41796.897
PKESM-PAKISTÁN_SWAT_02Swat13.1501.4095.47212.04063.98196.052
PKESM-PAKISTÁN_SWAT_03Swat13.1500.0535.01211.90366.83496.952
PKESM-PAKISTÁN_SWAT_04Swat13.1501.1577.1859.86165.68697.039
PKESM-PAKISTÁN_SWAT_05Swat13.1507.3620.60712.80063.00496.923
PKESM-PAKISTÁN_SWAT_06Swat13.1502.7836.41912.95362.25597.560
PKESM-PAKISTÁN_SWAT_07Swat13.1506.8026.69011.59358.44796.682
PKESM-PAKISTÁN_SWAT_08Swat13.1507.1900.75411.20664.77597.075
PKESM-PAKISTÁN_SWAT_09Swat13.1507.3507.18511.88957.45997.033
PKESM-PAKISTÁN_SWAT_10Swat13.1503.5146.94812.22260.58596.419
PKPromedio13.1504.1804.79812.00963.14496.863
PKDesviación Estándar0.0003.0072.7121.0333.5430.407
ZMBZAMBIA_KAGEM_1Kafubu13.150< LOD< LOD13.79770.52797.474
ZMBZAMBIA_KAGEM_2Kafubu13.150< LOD< LOD36.08347.27696.509
ZMBZAMBIA_KAGEM_3Kafubu13.150< LOD< LOD14.97766.97595.102
ZMBZAMBIA_KAGEM_4Kafubu13.150< LOD< LOD22.67761.48897.315
ZMBZAMBIA_KAGEM_6Kafubu13.150< LOD< LOD24.59556.57594.320
ZMBZAMBIA_KAGEM_7Kafubu13.150< LOD< LOD17.61166.34497.105
ZMBZAMBIA_KAGEM_8Kafubu13.150< LOD< LOD18.45564.19695.801
ZMBZAMBIA_KAGEM_11Kafubu13.150< LOD8.34417.07649.96588.535
ZMBZAMBIA_KAGEM_13Kafubu13.150< LOD< LOD23.92257.78394.855
ZMBZAMBIA_KAGEM_15Kafubu13.150< LOD7.20512.17762.57695.108
ZMBZAMBIA_KAGEM_16Kafubu13.150< LOD< LOD18.60666.16897.924
ZMBZAMBIA_KAGEM_17Kafubu13.150< LOD< LOD14.40568.37295.927
ZMBZAMBIA_KAGEM_18Kafubu13.150< LOD< LOD13.26970.25596.674
ZMBPromedio13.1500.0007.77519.05062.19295.588
ZMBDesviación Estándar0.000.000.816.527.402.40
Tabla 5. Composición de elementos menores por XRF en ppm. LOD= Límite de detección. 
Table 5. Minor element composition by XRF in ppm. LOD= Limit of detection. 
ORIGENMUESTRADISTRITOSClKCaScTiVCrMnFeNiCuZnGaCr/V
CEOCCEOMCATEDRAL_12Muzo - Q659.931046.05773.1571941.80< LOD< LOD3280.003890.00163.031699.5420.1125.5955.8284.311.19
CEOCCEOMCATEDRAL_17Muzo - Q253.381485.081027.184695.59< LOD75.231112.181547.19< LOD960.1316.66< LOD< LOD58.581.39
CEOCCEOMCATEDRAL_18Muzo - Q186.911620.74936.81157.21355.2539.041856.412130.00< LOD1240.6917.0020.7020.5575.741.15
CEOCCEOMCATEDRAL_19Muzo - Q543.581673.501178.97512.351017.5059.6311610.004540.00< LOD1636.6018.9428.3232.0966.750.39
CEOCCEOMCATEDRAL_20Muzo - Q2358.701428.681002.16188.87333.92< LOD4950.002160.00< LOD781.65< LOD18.8429.1456.890.44
CEOCCEOMCATEDRAL_21Muzo - Q370.661136.82844.23243.85319.95< LOD4490.002220.001321.6319.25< LOD22.85< LOD68.790.49
CEOCCEOMCATEDRAL_22Muzo - Q0.002460.001536.48317.17288.84< LOD1719.381878.533.071285.6824.58< LOD33.2276.961.09
CEOCCEOMCATEDRAL_23Muzo - Q437.941198.171002.33252.50397.35< LOD1779.421937.85< LOD1366.2517.3024.99< LOD70.941.09
CEOCCEOMCATEDRAL_24Muzo - Q496.151741.641211.11472.26239.3343.232060.001931.10< LOD1074.3415.0323.8729.7073.410.94
CEOCCEOMCATEDRAL_25Muzo - Q0.00962.17660.67604.46245.76< LOD2950.003710.00< LOD1045.6529.9135.0142.5468.821.26
CEOCCEOMACUNAS_01Maripí210.24992.00593.56735.431047.00< LOD11000.008550.00< LOD756.755.0014.00127.0063.000.78
CEOCCEOMACUNAS_02Maripí43.25878.00448.2886.48915.00< LOD10820.008770.00< LOD616.872.0012.0016.0057.000.81
CEOCCEOMACUNAS_03Maripí380.841904.001133.15676.82858.00< LOD8100.002650.00< LOD934.408.0017.0028.0072.000.33
CEOCCEOMACUNAS_04Maripí259.102020.00909.84434.541037.0062.0011530.008890.00< LOD723.184.0014.0035.0058.000.77
CEOCCEOMACUNAS_06Maripí973.121836.001693.501622.371023.00282.008110.006610.00< LOD1112.055.0022.0037.0068.000.82
CEOCCEOMACUNAS_07Maripí217.45784.00544.58172.24951.00< LOD10000.008120.00< LOD690.317.0012.0053.0057.000.81
CEOCCEOMACUNAS_Q_18Maripí267.912290.001101.612744.45489.0027.006220.005070.0049.00876.356.0022.00178.0067.000.82
CEOCCEOMACUNAS_19Maripí384.84766.00522.16909.81792.00< LOD10430.007660.00< LOD812.002.0075.0018.0065.000.73
CEOCCEOMACUNAS_20Maripí881.012260.001354.801251.44686.0054.0010010.007180.00< LOD1291.092.00376.0026.0060.000.72
CEOCCEOMACUNAS_21Maripí730.842330.001493.431693.84713.0042.0011080.008110.00< LOD1306.484.00341.0025.0071.000.73
CEOCCEOMACUNAS_22Maripí267.912200.001079.193037.48811.00110.006150.004420.00< LOD893.8321.00105.0021.0078.000.72
CEOCCEOPBCULATA_23_2Peñas Blancas209.763100.002282.9022191.47< LOD63.382930.002380.00< LOD1489.7237.9833.5744.1029.220.81
CEOCCEOPBCULATA_24Peñas Blancas0.001895.021132.019133.88< LOD< LOD2300.002320.00< LOD932.8029.55< LOD< LOD36.821.01
CEOCCEOPBCULATA_25Peñas Blancas396.041688.56666.56306.51128.08< LOD983.271235.32< LOD931.0946.0045.3948.8220.571.26
CEOCCEOPBCULATA_26Peñas Blancas253.203320.001455.31386.84219.43< LOD1820.374200.00< LOD706.8544.7640.2441.8219.382.31
CEOCCEOPBCULATA_27_2Peñas Blancas0.002620.001884.43415.83< LOD< LOD2670.002240.0053.201937.3448.3834.50< LOD34.870.84
CEOCCEOPBCULATA_28Peñas Blancas676.821852.911185.11313.67202.35< LOD2240.002230.00158.401156.3259.5044.9258.3234.251.00
CEOCCEOPBCULATA_29Peñas Blancas0.001600.06963.48405.05162.50< LOD1805.541975.416.05800.6991.6738.6083.6733.161.09
CEOCCEOPBCULATA_30Peñas Blancas0.001725.66853.21847.28125.70< LOD1030.051625.0384.82674.0544.4048.9844.9017.831.58
CEOCC Promedio395.161752.241085.184370.74556.5877.965346.094144.15229.901025.9323.2557.5547.0356.660.78
CEOCC Desviación Estándar463.61655.48421.2613690.42374.5557.373920.672623.27245.98386.7221.6687.8637.8419.420.67
CEOREOCHTIGRES_02Chivor< LOD3500.001657.65240.49172.7233.851314.585610.00< LOD631.188.91256.3543.1733.704.27
CEOREOCHTIGRES_07Chivor< LOD2210.001000.07253.60233.4037.871221.764690.00< LOD731.1614.8022.0519.1837.483.84
CEOREOCHTIGRES_08Chivor< LOD2020.00784.86122.2675.170.00993.264370.00< LOD598.1524.0420.7225.7935.774.40
CEOREOCHTIGRES_09Chivor< LOD1487.00795.79419.39198.3821.651116.394650.00< LOD656.9322.4314.9123.2034.184.17
CEOREOCHTIGRES_11Chivor< LOD1538.00783.37605.87158.490.001084.824640.00< LOD1092.7127.6132.4928.3225.724.28
CEORTIGRES_01Chivor< LOD1914.001116.383795.06113.1254.58933.042750.00< LOD1204.962.6914.4614.1729.342.95
CEORTIGRES_02Chivor< LOD2740.003403.601306.74195.57111.08845.822770.00< LOD1398.809.9823.6932.4527.033.27
CEORTIGRES_03Chivor< LOD2140.001776.512908.83293.70177.111024.524090.00< LOD796.593.3621.0317.7027.673.99
CEORTIGRES_04Chivor< LOD3380.001842.931658.1188.1363.69590.751877.00< LOD953.4216.1322.3026.1016.363.18
CEORTIGRES_05Chivor< LOD2080.001967.455295.93216.5112.831296.035790.00< LOD705.714.7714.0617.2228.934.47
CEORTIGRES_06Chivor< LOD3450.001892.741550.90258.34166.101189.004850.00< LOD800.284.3316.1416.5231.624.08
CEORACHIOTE_1Somondoco800.921217.00415.070.00275.00< LOD843.00442.00< LOD654.64< LOD28.00< LOD< LOD0.52
CEORACHIOTE_2Somondoco< LOD< LOD332.060.00< LOD< LOD863.00469.00< LOD763.05< LOD< LOD< LOD41.000.54
CEORACHIOTE_3Somondoco1281.472530.001328.24857.64307.00118.00600.00229.00< LOD1087.5741.0037.0038.00< LOD0.38
CEORACHIOTE_4Somondoco4004.591587.001245.221572.34221.00213.00633.00370.00< LOD1888.38< LOD34.0030.0032.000.58
CEORACHIOTE_5Somondoco4485.141281.001079.191357.93391.00192.00704.00514.00< LOD1930.34< LOD< LOD< LOD34.000.73
CEORACHIOTE_6Somondoco8249.462440.001992.353573.51< LOD337.00735.00693.00< LOD2713.6720.0041.0030.0029.000.94
CEORACHIOTE_7Somondoco0.001237.001245.221429.40317.00< LOD626.00648.00< LOD1559.66< LOD31.00< LOD35.001.04
CEORACHIOTE_8Somondoco2242.571228.00830.151000.58311.00< LOD645.00412.00< LOD1174.99< LOD35.00< LOD< LOD0.64
CEORACHIOTE_9Somondoco1281.47482.00498.09357.35234.00N.D1066.00626.00N.D1025.32N.DN.DN.D36.000.59
CEOREOCHPALOARAÑA_01Macanal392.351288.40768.94437.85118.8376.84463.541649.78N.D663.658.4519.0612.7228.493.56
CEOREOCHPALOARAÑA_02Macanal836.962610.001557.14440.80143.8456.79696.602160.00N.D< LOD359.3310.3922.7516.1728.423.10
CEOREOCHPALOARAÑA_03Macanal190.631061.66546.66129.7252.797.75513.831696.5154.17368.478.8911.9612.8522.753.30
CEOREOCHPALOARAÑA_04Macanal506.901261.05957.101253.3461.9079.05434.831434.3527.79559.65< LOD16.7720.3119.723.30
CEOREOCHPALOARAÑA_05Macanal88.701393.05920.80242.2042.14< LOD298.73838.4712.14471.2610.3719.8915.1714.032.81
CEOREOCHPALOARAÑA_06Macanal176.031668.361185.90260.01126.5942.10568.851145.5815.93405.927.9021.5223.4822.312.01
CEOREOCHPALOARAÑA_07Macanal320.951685.73924.623423.4262.0151.06322.721259.3329.18772.879.8020.4430.4120.103.90
CEOREOCHPALOARAÑA_08Macanal275.842750.001884.43304.47105.8296.80431.471932.30< LOD326.516.4832.5127.1616.444.48
CEOREOCHPALOARAÑA_09Macanal541.191888.701529.71832.7412.6399.60337.61825.7211.90648.577.7519.3217.2311.032.45
CEOREOCHPALOARAÑA_10Macanal588.671563.51982.57619.87179.4576.67707.601916.89< LOD668.0322.8040.4435.1726.912.71
CEORPromedio1459.101918.331241.491208.35177.3488.56770.022178.3025.19920.3913.3132.9223.8527.592.83
CEORDesviación Estándar1782.81822.70633.111321.36105.2880.02298.711783.5112.23534.4610.0544.1612.2511.105.97
PKESM-PAKISTÁN_SWAT_01Swat1399.73596.63175.663050.0058.2420540.0040.6920.0849.607566.11< LOD1102.54< LOD< LOD0.49
PKESM-PAKISTÁN_SWAT_02Swat1123.751137.60259.873900.0085.0626430.00103.4428.6742.018549.47< LOD69.50< LOD< LOD0.28
PKESM-PAKISTÁN_SWAT_03Swat1603.31490.04180.483580.0044.4120470.0073.0923.1631.037329.0112.62572.98< LOD< LOD0.32
PKESM-PAKISTÁN_SWAT_04Swat1009.82602.33140.543110.0062.9620600.0048.9521.9924.8762785.14< LOD352.57< LOD< LOD0.45
PKESM-PAKISTÁN_SWAT_05Swat1519.47780.43< LOD2550.0072.9721410.0059.2720.5544.309967.8528.46784.27< LOD< LOD0.35
PKESM-PAKISTÁN_SWAT_06Swat859.15367.9163.732070.0075.9617210.0033.8227.4654.388178.08< LOD700.36< LOD< LOD0.81
PKESM-PAKISTÁN_SWAT_07Swat2270.004350.00195.592620.0090.3717650.0056.9330.5651.107328.31120.14868.72< LOD< LOD0.54
PKESM-PAKISTÁN_SWAT_08Swat1023.46230.72241.822170.00< LOD21060.0059.2830.2236.4566939.57< LOD174.41< LOD< LOD0.51
PKESM-PAKISTÁN_SWAT_09Swat1363.88345.41195.403690.0049.3519950.0096.4229.0260.737215.7126.10596.72< LOD< LOD0.30
PKESM-PAKISTÁN_SWAT_10Swat2430.001238.10281.344470.0059.8022670.00177.9217.4055.9860693.9311.841064.63< LOD< LOD0.10
PKPromedio1460.261013.92192.713121.0066.5720799.0074.9824.9145.0424655.3239.83628.670.000.000.33
PKDesviación Estándar527.181217.7086.93785.4425.692568.6842.474.8111.5326840.7836.91349.760.000.000.11
ZMBZAMBIA_KAGEM_1Kafubu< LOD623.741046.28633.20149.14< LOD210.46831.9462.4114085.9244.9021.2052.34< LOD3.95
ZMBZAMBIA_KAGEM_2Kafubu< LOD0.870.471299.10475.35< LOD89.093210.00< LOD8525.6979.4053.02105.69< LOD36.03
ZMBZAMBIA_KAGEM_3Kafubu< LOD1864.040.43947.43398.21< LOD276.161315.00< LOD16925.4894.2735.45122.8326.554.76
ZMBZAMBIA_KAGEM_4Kafubu< LOD0.240.20< LOD< LOD< LOD92.473160.000.0011253.3563.76< LOD147.91< LOD34.17
ZMBZAMBIA_KAGEM_6Kafubu< LOD0.450.372.21< LOD< LOD133.663770.00416.5413176.7095.3251.22140.87< LOD28.21
ZMBZAMBIA_KAGEM_7Kafubu< LOD932.290.21801.43997.71< LOD107.164670.00177.5112050.6660.43< LOD128.5322.5843.58
ZMBZAMBIA_KAGEM_8Kafubu< LOD0.560.610.27< LOD< LOD130.323860.00< LOD14673.4198.06< LOD131.00< LOD29.62
ZMBZAMBIA_KAGEM_11Kafubu< LOD1873.461281.68< LOD< LOD< LOD46.024040.00191.256280.6126.3654.92279.95< LOD87.79
ZMBZAMBIA_KAGEM_13Kafubu< LOD0.970.76< LOD< LOD< LOD209.311359.00< LOD17939.61< LOD< LOD71.76< LOD6.49
ZMBZAMBIA_KAGEM_15Kafubu< LOD1744.010.50< LOD386.65< LOD392.755260.00107.2420982.0087.9124.5439.3618.1613.39
ZMBZAMBIA_KAGEM_16Kafubu< LOD1353.371173.35< LOD< LOD< LOD24.714910.00< LOD7091.9252.7924.20402.8818.04198.72
ZMBZAMBIA_KAGEM_17Kafubu< LOD976.360.26< LOD< LOD< LOD418.312690.00< LOD17533.9631.4428.3149.4815.866.43
ZMBZAMBIA_KAGEM_18Kafubu< LOD0.280.27< LOD332.16< LOD400.762620.0014.5513890.0834.6529.5438.3816.406.54
ZMBPromedio#¡DIV/0!720.82269.65613.94456.540.00194.713207.38138.5013416.1164.1135.82131.6119.6016.47
ZMBDesviación Estándar0.00780.74513.89464.11300.670.00137.511412.52124.134384.2330.9620.46104.2910.5110.27

Finalmente, se realizaron análisis de componentes principales (PCA) en RStudio de los datos obtenidos por los diferentes métodos: se desarrolló un PCA para los espectros de fotoluminiscencia, otro para la concentración química de las esmeraldas y dos más para la combinación de la química y la fotoluminiscencia (orientados de manera perpendicular y paralelo al eje C), con la finalidad de determinar diferencias que permitan agrupar esmeraldas de diferentes proveniencias. Para hacer la combinación de datos de diferentes orígenes (fotoluminiscencia y química), se escalan estos datos con el fin asegurar que todas las variables influyan de manera similar en el PCA. Esto previene que algunas de las variables dominen sobre el análisis y además facilita la interpretación de las relaciones entre las observaciones. El proceso de escalado se realizó en RStudio mediante el método de Z-score, en donde cada variable se transforma para que tenga una media de 0 y una desviación estándar de 1. Esto permite que cada variable tenga una contribución equitativa al análisis, independientemente de su escala original. Posteriormente, se calcula la matriz de covarianza o de correlación de los datos; a continuación, se procede con la descomposición de valores propios utilizando la ecuación característica de una matriz para obtener los componentes principales, ordenados de mayor a menor según la cantidad de varianza que explican. Posteriormente, se seleccionan los primeros componentes principales (PC1 y PC2) que capturan la mayor parte de la varianza total, permitiendo así la reducción de la dimensionalidad del conjunto de datos. Finalmente, se transforman los datos originales en el espacio de los componentes principales para facilitar la interpretación y el análisis de patrones significativos presentes en los datos. Este análisis se hace debido a que las esmeraldas de distintas regiones geológicas tienen composiciones químicas únicas, como consecuencia de las variaciones en los minerales y elementos presentes en sus lugares de origen. El PCA permite separar las muestras basándose en esas diferencias en sus composiciones químicas; además, esas variaciones en la composición, también se ven reflejadas en la respuesta de fotoluminiscencia, debido a que la presencia de elementos como Mg y Fe puede influir en la intensidad y la forma de los picos de emisión en los espectros (Moroz et al., 2000Moroz, I., Roth, M., Boudeulle, M., & Panczer, G. (2000). Raman Microspectroscopy and Flourescence of Emeralds from Various Deposits. Journal of Raman Spectroscopy, 31, 485-490.; Thompson et al., 2014Thompson, D. B., Kidd, J. D., Åström, M., Scarani, A., & Smith, C. P. (2014). A Comparison of R-line Photoluminescence of Emeralds from Different Origins. The Journal of Gemmology, 34(4), 334-343. 10.15506/jog.2014.34.4.334.)

3. Resultados

 

3.1. Fluorescencia con lámpara UV

 

Las esmeraldas que muestran mayor reacción a la luz UV son las del CEOR de Colombia y las que menos muestran reacción son las de Zambia y Pakistán. Los colores correspondientes a la reacción generada con la lámpara UV en las esmeraldas de Colombia, están en el rango de rosa intenso a rosa débil o rojo de medio a débil. Las esmeraldas de Zambia y Pakistán están en el rango de inerte a rojo débil (Fig. 2, Tabla 1).

media/001_002.jpg
Figura 2 Respuesta a excitación con lámpara UV (365 nm) de las muestras de estudio. 
Figure 2 The reaction of the samples to a UV lamp (365 nm). 

3.2. Fotoluminiscencia

 
3.2.1. PERPENDICULARMENTE AL EJE C
 

Los espectros tomados con una orientación perpendicular al eje C (Fig. 3), presentan una intensidad notablemente mayor en el pico R1 respecto al R2. Sin embargo, hay que destacar que, para las esmeraldas de Colombia, los picos más sobresalientes son R1 y R2 ubicados a 683 nm y 680 nm respectivamente, mientras que para las esmeraldas de Pakistán y Zambia el pico R1 se encuentra ubicado a 684 nm aproximadamente; además, en estas esmeraldas (Pakistán y Zambia) se presenta una banda con centro alrededor de 675 nm (denominada en este trabajo R3) la cual tiene mayor intensidad que el pico R2. Vale la pena aclarar que la banda R3 se genera por la resta de la línea base con el espectro original, y está altamente relacionada con una caída de la intensidad en una longitud de onda cercana a 681.80 nm (entre los picos R2 y R1). Entre menos caída de intensidad en esa zona, la banda R3 se vuelve más alta, como lo que ocurre con las esmeraldas de Zambia y Pakistán, mientras que, en las esmeraldas colombianas, esta banda no se hace tan evidente. Además, la posición de los picos cambia entre las diferentes áreas geográficas: para las esmeraldas colombianas el pico R1 se encuentran en una longitud de onda menor a 684 nm, mientras que, para las esmeraldas de los otros países, esta longitud de onda es mayor a 684 nm.

media/001_003.jpeg
Figura 3 Espectros orientados perpendicularmente al eje cristalográfico C. 
Figure 3 Spectra obtained perpendicular to the C axis. 

El promedio de la posición de los picos R1 y R2 muestra rangos que pueden ayudar a distinguir entre las diferentes áreas de estudio (Tabla 2): las esmeraldas de CEOR muestran un promedio del pico R1 de 683.603 nm, las de CEOCC 683.735 nm, las de Pakistán están en 684.187 nm y las de Zambia en 684.132 nm. Es decir que: 1) para las esmeraldas colombianas los picos de R1 se encuentran en una longitud de onda por debajo de 684 nm mientras que las esmeraldas de los otros dos países están por encima de 684 nm; 2) entre las esmeraldas colombianas, también hay distinción, debido a que las muestras de CEOR presentan longitudes de onda en general menores a las de CEOCC. Lo contrario ocurre con el pico R2, en donde las esmeraldas colombianas presentan mayores longitudes de onda, en promedio, 680.60 nm en comparación con los otros dos países los cuales tienen un promedio de 680.25 nm. La distancia entre los picos (R1-R2) es un factor que también presenta diferencias entre las áreas de estudio: las muestras colombianas se encuentran muy cercanas a 3.1 nm mientras que las de otros países están más cercanas a 3.9 nm (Fig. 3, Tabla 2). El ancho de altura medio del pico, FWHM por sus siglas en inglés (Full Width at Half Maximum), está representado por (W) y es un parámetro que también muestra diferencias. Las esmeraldas del CEOR tienen menor ancho que las de todas las demás localidades: 1.7 nm en R1 y 1.2 nm en R2, mientras que el mayor ancho lo presentan Pakistán y Zambia, con un ancho alrededor de R1 2.8 nm y R2 de 1.8 nm (Fig. 3, Tabla 2).

Tabla 2. Datos de tratamiento de los espectros tomados perpendicularmente al eje C. R1= pico de emisión localizado entre 683 y 685 nm. R2= pico de emisión localizado entre 679 a 681 nm. R3= pico localizado cercano a 674 nm. I = intensidad de los picos, W= ancho de pico (en altura media), CEOCC= Cinturón esmeraldífero occidental, CEOR= Cinturón esmeraldífero oriental, PK=Pakistán, ZMB= Zambia. 
Table 2. Data treatment of the spectra taken perpendicularly to the C axis. R1= emission peak located between 683 and 685 nm. R2= emission peak located between 679 and 681 nm. R3= peak located near 674 nm. I = peak intensity, W= peak width (at half height), CEOCC= Western Emerald Belt, CEOR= Eastern Emerald Belt, PK= Pakistan, ZMB= Zambia. 
ORIGENMUESTRAR1R2I1I2W1W2ÁREA 1ÁREA R2R1-R2R1/R2I1/I2W1/W2A1/A2
CEOCCEOMCATEDRAL_12683.568680.5484.590.8501.651.5201.7061.1728.337.7602.059.9303.0201.0042.7801.4564.048
CEOCCEOMCATEDRAL_17683.723680.5742.010.410834.8962.0161.4064.314.5001.249.6503.1491.0052.4081.4343.453
CEOCCEOMCATEDRAL_18683.658680.6002.523.350713.6301.8611.2764.998.910969.1103.0581.0043.5361.4595.158
CEOCCEOMCATEDRAL_19683.813680.470883.254284.1672.2231.4852.089.730449.1133.3431.0053.1081.4974.653
CEOCCEOMCATEDRAL_20683.555680.6135.673.7901.661.1901.6291.1989.836.4102.117.5302.9421.0043.4151.3604.645
CEOCCEOMCATEDRAL_21683.710680.6527.208.7002.019.0301.9641.40615.073.3003.021.5103.0581.0043.5701.3974.989
CEOCCEOMCATEDRAL_22683.710680.5482.338.0201.038.2502.0161.3545.017.5901.496.5503.1621.0052.2521.4893.353
CEOCCEOMCATEDRAL_23683.684680.6656.142.9701.829.0201.9391.45812.676.5002.838.6403.0191.0043.3591.3304.466
CEOCCEOMCATEDRAL_24683.568680.5483.783.3501.554.5201.7061.1986.871.2001.981.9703.0201.0042.4341.4243.467
CEOCCEOMCATEDRAL_25683.684680.6651.789.590533.3141.8871.4063.594.350798.1143.0191.0043.3561.3424.504
CEOCCCUNAS_01_PL_600_PER_3683.853680.588166.90074.5032.4301.667431.636132.1833.2651.0052.2401.4583.265
CEOCCCUNAS_02_PL_600_PER683.930680.536103.48345.5242.5581.589281.82576.9893.3941.0052.2731.6103.661
CEOCCCUNAS_03_PL_600_PER683.840680.627246.67183.6172.3001.562603.962139.0533.2131.0052.9501.4724.343
CEOCCCUNAS_04_PL_600_PER683.827680.679321.04994.0502.2741.666777.271166.8263.1481.0053.4141.3654.659
CEOCCCUNAS_06_PL_600_PER683.801680.653641.205172.6012.0931.4061.428.820258.3013.1481.0053.7151.4895.532
CEOCCCUNAS_07_PL_600_PER683.853680.653295.41284.2082.2231.614698.889144.6973.2001.0053.5081.3774.830
CEOCCCUNAS_18_PL_600_PER683.711680.8892.395.460720.4622.0161.2245.140.860923.0502.8221.0043.3251.6475.569
CEOCCCUNAS_19_PL_600_PER683.802680.627871.589273.7642.2751.5622.110.260455.2973.1751.0053.1841.4564.635
CEOCCCUNAS_20_PL_600_PER683.814680.6661.038.960290.0362.3261.6142.572.610498.3763.1481.0053.5821.4415.162
CEOCCCUNAS_21_PL_600_PER683.763680.640488.330164.4652.2231.6411.155.490287.2003.1231.0052.9691.3554.023
CEOCCCUNAS_22_PL_600_PER683.736680.626430.961142.4082.0681.615948.567244.7393.1101.0053.0261.2813.876
CEOCCEOPBCULATA_23683.774680.7174.672.2201.289.2802.1451.58810.669.3002.179.3903.0571.0043.6241.3514.896
CEOCCEOPBCULATA_24683.658680.5483.239.8301.340.6201.9131.3286.597.2301.895.1303.1101.0052.4171.4403.481
CEOCCEOPBCULATA_25683.632680.5612.847.6301.186.6801.9131.3285.798.9701.677.5103.0711.0052.4001.4413.457
CEOCCEOPBCULATA_26683.607680.5485.521.9902.337.4801.8101.22410.637.3003.044.9203.0591.0042.3621.4793.493
CEOCCEOPBCULATA_27683.826680.6654.812.0201.290.1702.1971.61411.252.4002.216.9303.1611.0053.7301.3615.076
CEOCCEOPBCULATA_28683.723680.65213.721.7003.965.4002.0421.43229.825.9006.044.0303.0711.0053.4601.4264.935
CEOCCEOPBCULATA_29683.723680.6657.859.6902.148.8601.9641.38016.434.5003.156.0603.0581.0043.6581.4245.207
CEOCCEOPBCULATA_30683.762680.6134.901.0401.721.5902.0681.48410.787.4002.720.0103.1491.0052.8471.3933.966
CEOCCPromedio683.735680.6223.155.8771.018.8022.0621.4456.584.9461.491.1313.1131.0053.0661.4294.372
CEOCCDesviación Estándar0.0950.0763.077.449925.0840.2210.1576.492.4841.359.8930.1130.0000.5150.0780.712
CEORTIGRES_2_PL_600683.659680.6012.844.950854.6841.8351.2245.558.0301.113.2403.0581.0043.3291.5004.993
CEORTIGRES_7_PL_601683.582680.6272.841.630813.0181.6291.1714.926.4101.013.7502.9551.0043.4951.3904.860
CEORTIGRES_8_PL_602683.672680.6663.028.490910.6781.8351.2765.916.1901.236.4603.0061.0043.3261.4394.785
CEORTIGRES_9_PL_603683.659680.6141.029.560401.6451.8091.2761.983.030545.4343.0451.0042.5631.4183.636
CEORTIGRES_11_PL_604683.621680.6142.669.200910.0271.7841.2765.067.9501.235.8203.0071.0042.9331.3984.101
CEORTIGRES_CANU_01683.634680.6141.665.450599.3671.6801.1982.978.640764.0193.0201.0042.7791.4033.899
CEORTIGRES_CANU_02683.647680.6142.050.300623.1701.8091.2503.949.080828.9373.0331.0043.2901.4484.764
CEORTIGRES_CANU_03683.569680.6141.495.520525.8591.6291.2242.592.720684.9402.9551.0042.8441.3313.785
CEORTIGRES_CANU_04683.621680.5751.611.340558.7181.7061.1462.926.250681.3493.0461.0042.8841.4894.295
CEORTIGRES_CANU_05683.634680.588860.303326.3461.7581.2241.609.740425.1133.0461.0042.6361.4363.787
CEORTIGRES_CANU_06683.686680.6547.859.8302.489.2401.7831.30214.920.2003.448.9603.0321.0043.1581.3704.326
CEORACHIOTE_T_GOCEN_01683.620680.5744.848.4102.184.7401.7841.2769.205.5902.967.1603.0461.0042.2191.3983.102
CEORACHIOTE_T_GOCEN_02683.620680.6268.393.4002.770.2701.7321.27615.474.0003.761.6702.9941.0043.0301.3584.114
CEORACHIOTE_T_GOCEN_03683.529680.6006.841.9802.221.0001.5251.11911.108.6002.646.4302.9291.0043.0811.3634.198
CEORACHIOTE_T_GOCEN_04683.620680.6135.542.9501.752.0701.7321.25010.218.9002.330.6003.0071.0043.1641.3864.385
CEORACHIOTE_T_GOCEN_05683.632680.6135.996.8801.905.0601.7581.25011.220.2002.534.0903.0191.0043.1481.4074.428
CEORACHIOTE_T_GOCEN_06683.632680.6136.294.9102.001.3701.7581.22411.777.8002.606.5503.0191.0043.1451.4374.519
CEORACHIOTE_T_GOCEN_07683.620680.6267.969.0702.532.2901.7581.25014.911.2003.368.1102.9941.0043.1471.4074.427
CEORACHIOTE_T_GOCEN_08683.632680.62610.596.5003.541.7601.8091.32820.409.9005.005.7703.0061.0042.9921.3634.077
CEORACHIOTE_T_GOCEN_09683.620680.6267.176.9002.371.0501.7581.27613.428.9003.219.5702.9941.0043.0271.3784.171
CEOREOCHPALOARAÑADO_1683.595680.60114.580.2004.104.1801.7581.25027.283.5005.459.8902.9941.0043.5531.4074.997
CEOREOCHPALOARAÑADO_2683.531680.5369.422.6303.665.5201.6811.19816.856.0004.673.4202.9951.0042.5711.4033.607
CEOREOCHPALOARAÑADO_3683.569680.60121.219.7006.209.7401.7321.22439.125.9008.088.2902.9681.0043.4171.4164.837
CEOREOCHPALOARAÑADO_4683.582680.60115.564.4004.462.7401.7581.27629.127.0006.060.4102.9811.0043.4881.3784.806
CEOREOCHPALOARAÑADO_5683.595680.61412.627.1003.582.8401.7321.22423.280.9004.666.7202.9811.0043.5241.4164.989
CEOREOCHPALOARAÑADO_6683.543680.60114.501.0004.194.1701.6811.22425.940.6005.462.9702.9421.0043.4571.3734.748
CEOREOCHPALOARAÑADO_7683.569680.60120.270.8005.817.5701.7321.22437.376.3007.577.4802.9681.0043.4841.4164.933
CEOREOCHPALOARAÑADO_8683.543680.58813.676.1003.957.8901.7061.22424.841.6005.155.7302.9551.0043.4551.3944.818
CEOREOCHPALOARAÑADO_9683.543680.5758.490.2902.686.4101.6811.22415.188.1003.499.4402.9681.0043.1601.3734.340
CEOREOCHPALOARAÑADO_10683.505680.58815.373.8004.612.2001.6291.22426.656.8006.008.0602.9171.0043.3331.3314.437
CEORPromedio683.603680.6077.911.4532.452.8541.7321.23714.528.6683.235.6792.9961.0043.1211.4014.372
CEORDesviación Estándar0.0460.0255.796.8421.662.9960.0690.04410.612.5182.188.9230.0370.0000.3330.0390.483
PKPAKISTAN_SWAT_1684.176680.2631.794.740419.7232.8672.0075.478.070896.8353.9131.0064.2761.4286.108
PKPAKISTAN_SWAT_2684.240680.3153.646.270686.3342.8151.98110.926.7001.447.2503.9251.0065.3131.4217.550
PKPAKISTAN_SWAT_3684.214680.2762.111.310465.0492.8671.9556.443.510967.8203.9381.0064.5401.4666.658
PKPAKISTAN_SWAT_4684.098680.1852.205.420584.4782.7641.8516.489.7001.151.6603.9131.0063.7731.4935.635
PKPAKISTAN_SWAT_5684.098680.1591.580.470413.8862.7391.7994.607.200792.5593.9391.0063.8191.5225.813
PKPAKISTAN_SWAT_6684.176680.1982.706.100452.8042.7641.8517.961.610892.1493.9781.0065.9761.4938.924
PKPAKISTAN_SWAT_7684.176680.2112.458.640363.0682.7641.7997.233.570695.1523.9651.0066.7721.53710.406
PKPAKISTAN_SWAT_8684.227680.1591.982.080279.1902.8411.8775.994.520557.9184.0681.0067.0991.51310.744
PKPAKISTAN_SWAT_9684.227680.2763.360.100577.5562.9451.92910.532.4001.185.9003.9511.0065.8181.5278.881
PKPAKISTAN_MINGORA_18684.240680.262862.823161.3872.9191.6682.680.560286.5643.9781.0065.3461.7509.354
PKPromedio684.187680.2302.270.795440.3482.8281.8726.834.784887.3813.9571.0065.2731.5158.007
PKDesviación Estándar0.0530.055823.402153.3920.0710.1022.514.039332.4090.0460.0001.1680.0911.907
ZMBZAMBIA_KAGEM_1_PER684.188680.223925.723143.2952.2711.8252.697.950278.3273.9651.0066.4601.2459.693
ZMBZAMBIA_KAGEM_2_PER684.278680.2621.638.030245.4602.2712.0335.268.970531.2764.0161.0066.6731.1179.918
ZMBZAMBIA_KAGEM_3_PER684.265680.3141.020.970181.9332.2711.9033.227.990368.4593.9511.0065.6121.1948.761
ZMBZAMBIA_KAGEM_4_PER684.226680.3272.442.480372.0442.2711.8777.252.040743.1823.8991.0066.5651.2109.758
ZMBZAMBIA_KAGEM_6_PER684.071680.5483.316.580631.3402.5321.9548.937.5801.312.9903.5231.0055.2531.2966.807
ZMBZAMBIA_KAGEM_7_PER684.149680.4835.690.340991.6162.6351.77115.958.5001.869.8803.6661.0055.7381.4878.535
ZMBZAMBIA_KAGEM_8_PER684.239680.3272.083.040376.7222.9192.0596.471.730825.7573.9121.0065.5291.4177.837
ZMBZAMBIA_KAGEM_11_PER683.801680.7304.783.7101.342.7302.1711.61411.055.1002.306.9003.0711.0053.5631.3454.792
ZMBZAMBIA_KAGEM_13_PER684.278680.285496.07079.1862.9441.9751.554.690166.4733.9931.0066.2651.4919.339
ZMBZAMBIA_KAGEM_15_PER684.188679.7911.086.850118.9102.5831.1372.987.900145.9884.3971.0069.1402.27120.467
ZMBZAMBIA_KAGEM_16_PER683.801680.7304.125.0101.224.7302.1191.6149.305.5502.104.1703.0711.0053.3681.3134.422
ZMBZAMBIA_KAGEM_17_PER684.239680.210573.565144.4823.0482.0081.860.990308.7584.0291.0063.9701.5186.027
ZMBZAMBIA_KAGEM_18_PER683.994680.418487.935289.4624.1792.4762.170.450762.9733.5761.0051.6861.6882.845
ZMBPromedio684.132680.3582.205.408472.4552.6321.8656.057.649901.9333.7751.0065.3711.4308.400
ZMBDesviación Estándar0.1680.2441.754.309438.0160.5590.3094.337.637755.5670.3830.0011.8790.2984.288

Por otra parte, las intensidades (I) en las esmeraldas del CEOR son mayores que en las muestras de las demás áreas de estudio, con valores por encima de 7000 cuentas llegando hasta aproximadamente 25000 cuentas. En el CEOCC, los valores que se presentan son intermedios, y oscilan entre 2000 hasta los 13000 aproximadamente, con un promedio de 3155 cuentas. Mientras que las más bajas las presenta Pakistán y Zambia con un promedio de 2270 y 2200 cuentas respectivamente. Es necesario indicar que los valores encontrados para el área y la intensidad de los picos presentan rangos bastantes amplios, lo que denota poco agrupamiento en los datos (desviaciones estándar altas). Esto puede ser debido, entre otros aspectos, a la saturación de los cristales. Sin embargo, la relación de estos valores I1/I2 y A1/A2 presenta una menor desviación (Fig. 3, Tabla 2).

3.2.2. PARALELO AL EJE C
 

Con respecto a los espectros tomados paralelamente al eje C (Fig. 4), la posición del pico R1 en las muestras de Zambia y Pakistán siguen presentándose en longitudes mayores a las muestras de Colombia. También muestran diferencias con respecto al pico R2, el cual se encuentra ubicado a una longitud de onda menor que las muestras de Colombia, es decir, las muestras de Pakistán y Zambia se encuentran en promedio a 680.254 y 680.386 nm respectivamente, mientras que las colombianas se encuentran en promedio a 680.480 nm. Las intensidades, aunque disminuyen bastante con respecto a la orientación perpendicular, también presentan una alta desviación que puede ser originada, como se mencionó con anterioridad, por la saturación de la muestra. Sin embargo, la relación entre intensidades y áreas es un mejor parámetro para este estudio. El ancho de altura media de los picos R1 y R2 presenta una diferencia entre las esmeraldas de Colombia, con promedio entre 2.2 y 1.4 nm respectivamente, y las de Pakistán y Zambia, con promedio de los picos R1 y R2 de 3.7 nm - 2.3 nm y 3.5 nm - 2.2 nm respectivamente (Fig. 4, Tabla 3).

media/001_004.jpg
Figura 4. Espectros orientados paralelo al eje cristalográfico C. 
Figure 4. Spectra obtained parallel to the C axis. 
Tabla 3. Datos de tratamiento de los espectros tomados paralelamente al eje C (las abreviaturas son las mismas de la Tabla 2). 
Table 3. Data of spectra taken parallel to the C-axis (abbreviations are the same as in Table 2). 
ORIGENMUESTRAR1R2I1I2W1W2ÁREA 1ÁREA R2R1-R2R1/R2I1/I2W1/W2A1/A2
CEOCCEOMCATEDRAL_12683.594680.444844.2531.034.9302.4751.3542.224.4301.492.0403.1501.0050.8161.8281.491
CEOCCEOMCATEDRAL_17683.684680.496588.419599.3572.3271.4591.457.460930.5143.1881.0050.9821.5951.566
CEOCCEOMCATEDRAL_18683.671680.509442.763503.2952.3461.4851.105.660795.3053.1621.0050.8801.5801.390
CEOCCEOMCATEDRAL_19683.723680.353284.205317.5572.7321.642826.373554.9093.3701.0050.8951.6641.489
CEOCCEOMCATEDRAL_20683.607680.496810.2161.014.4102.3981.4062.068.0001.518.6103.1111.0050.7991.7051.362
CEOCCEOMCATEDRAL_22683.710680.496429.807541.0352.2171.4591.014.330839.9693.2141.0050.7941.5201.208
CEOCCEOMCATEDRAL_23683.658680.509838.099855.6782.3271.4852.076.1201.352.1403.1491.0050.9791.5681.535
CEOCCEOMCATEDRAL_24683.568680.457755.148943.9172.5271.3802.031.1601.386.9503.1111.0050.8001.8311.464
CEOCCEOMCATEDRAL_25683.452680.509354.353453.8061.8621.250702.364603.8272.9431.0040.7811.4901.163
CEOCCCUNAS_01_PL_600_PER_3683.724680.575122.244125.9382.5861.693336.464226.9403.1491.0050.9711.5271.483
CEOCCCUNAS_02_PL_600_PER683.685680.536154.942117.7882.0421.432336.827179.5783.1491.0051.3151.4261.876
CEOCCCUNAS_03_PL_600_PER683.647680.497136.969148.7422.3011.563335.523247.4503.1501.0050.9211.4721.356
CEOCCCUNAS_06_PL_600_PER683.517680.41980.91091.5591.9911.224171.505119.3173.0981.0050.8841.6271.437
CEOCCCUNAS_07_PL_600_PER683.595680.56256.49866.3072.5871.745155.555123.1733.0331.0040.8521.4821.263
CEOCCCUNAS_18_PL_600_PER683.685680.445287.008325.2192.5781.459787.457504.9943.2401.0050.8831.7671.559
CEOCCCUNAS_19_PL_600_PER683.711680.406231.620265.5352.6291.615648.161456.5863.3051.0050.8721.6271.420
CEOCCCUNAS_20_PL_600_PER683.647680.406201.666257.2622.1921.250470.475342.3763.2411.0050.7841.7531.374
CEOCCCUNAS_21_PL_600_PER683.711680.445166.170167.3992.8091.824496.855324.9943.2661.0050.9931.5401.529
CEOCCEOPBCULATA_23683.594680.457273.839296.7372.1981.354640.699427.8013.1371.0050.9231.6231.498
CEOCCEOPBCULATA_24683.658680.470711.354814.3612.4231.3801.834.8901.196.5903.1881.0050.8741.7551.533
CEOCCEOPBCULATA_25683.620680.496909.729964.8422.2501.4322.2501.471.1103.1241.0050.9431.5710.002
CEOCCEOPBCULATA_26683.607680.4961.326.6701.515.5902.2501.3803.177.1002.226.7603.1111.0050.8751.6301.427
CEOCCEOPBCULATA_27683.581680.457845.1621.062.6902.3281.3022.093.9401.473.1003.1241.0050.7951.7871.421
CEOCCEOPBCULATA_28683.555680.483537.773612.3291.9651.3021.124.900848.7293.0721.0050.8781.5091.325
CEOCCEOPBCULATA_29683.529680.496818.1171.027.0702.1721.3021.891.9301.423.5903.0331.0040.7971.6681.329
Promedio683.629680.477488.317564.9342.3401.4471.120.417842.6943.1531.0050.8641.6171.330
Desviación Estándar0.0710.050337.090399.4770.2400.158837.429574.5490.0900.0000.1110.1160.319
CEORTIGRES_2_PL_600683.543680.52365.12268.2021.8101.380125.465100.2053.0201.0040.9551.3111.252
CEORTIGRES_7_PL_601683.608680.510161.289128.7401.8871.276324.029174.8633.0981.0051.2531.4791.853
CEORTIGRES_8_PL_602683.582680.52365.40276.0792.0691.328144.013107.5573.0591.0040.8601.5581.339
CEORTIGRES_11_PL_604683.543680.53690.89985.3261.9131.432130.087130.0873.0071.0041.0651.3361.000
CEORTIGRES_CANU_01683.517680.523417.814522.8162.0951.276931.682710.1212.9941.0040.7991.6421.312
CEORTIGRES_CANU_02683.582680.49792.624102.4882.2241.459219.276159.1293.0851.0050.9041.5251.378
CEORTIGRES_CANU_03683.530680.484164.371204.3672.6821.563469.179339.9893.0461.0040.8041.7161.380
CEORTIGRES_CANU_04683.595680.536101.689126.5502.3211.380251.201185.9143.0591.0040.8041.6811.351
CEORTIGRES_CANU_05683.466680.51092.303117.6751.8621.250182.941156.5762.9561.0040.7841.4901.168
CEORTIGRES_CANU_06683.505680.497125.249140.9171.9651.250262.042187.5013.0081.0040.8891.5721.398
CEORACHIOTE_T_GOCEN_01683.542680.4831.106.0101.383.4702.0431.1982.405.0801.764.2403.0591.0040.7991.7051.363
CEORACHIOTE_T_GOCEN_02683.568680.4961.496.4201.742.9302.0171.2503.212.6202.319.1203.0721.0050.8591.6131.385
CEORACHIOTE_T_GOCEN_03683.542680.509635.089754.3091.9651.2501.328.5401.003.6703.0331.0040.8421.5721.324
CEORACHIOTE_T_GOCEN_04683.607680.509785.974829.7102.1201.3541.773.8801.196.0703.0981.0050.9471.5661.483
CEORACHIOTE_T_GOCEN_05683.516680.496623.620811.1101.9911.1981.321.8101.034.2503.0201.0040.7691.6621.278
CEORACHIOTE_T_GOCEN_06683.568680.470845.210967.2122.3281.3802.094.1701.421.3103.0981.0050.8741.6861.473
CEORACHIOTE_T_GOCEN_07683.645680.457865.644990.7852.4231.3802.232.9901.455.8203.1881.0050.8741.7561.534
CEORACHIOTE_T_GOCEN_08683.632680.4701.152.9701.163.2802.3271.3542.856.2701.677.0803.1621.0050.9911.7181.703
CEORACHIOTE_T_GOCEN_09683.581680.509581.732733.2382.1981.3021.361.1501.016.3203.0721.0050.7931.6881.339
CEOREOCHPALOARAÑADO_1683.466680.4321.513.0101.732.3301.9661.2503.165.6602.305.4503.0341.0040.8731.5721.373
CEOREOCHPALOARAÑADO_2683.466680.4451.227.8901.367.2601.9401.2762.535.2801.857.4603.0211.0040.8981.5201.365
CEOREOCHPALOARAÑADO_3683.492680.4581.090.8101.165.2701.9921.3022.313.0601.614.9903.0341.0040.9361.5301.432
CEOREOCHPALOARAÑADO_4683.479680.4451.471.6201.932.3502.2761.3023.565.9102.678.8703.0341.0040.7621.7481.331
CEOREOCHPALOARAÑADO_5683.453680.4971.146.600495.4401.4741.0941.798.890576.8222.9561.0042.3141.3483.119
CEOREOCHPALOARAÑADO_6683.531680.4711.302.5401.217.4101.9391.3022.688.9101.687.4103.0601.0041.0701.4891.594
CEOREOCHPALOARAÑADO_7683.505680.4452.197.4302.061.6701.8361.2504.294.4202.743.4903.0601.0041.0661.4691.565
CEOREOCHPALOARAÑADO_8683.531680.4321.822.1802.014.5002.0691.2764.012.8502.737.0003.0991.0050.9051.6211.466
CEOREOCHPALOARAÑADO_9683.505680.4451.313.8901.618.1402.1211.3022.966.0402.243.2703.0601.0040.8121.6281.322
CEOREOCHPALOARAÑADO_10683.427680.4321.762.2202.247.4601.9141.1983.590.4502.866.3102.9951.0040.7841.5981.253
Promedio683.535680.484838.539924.1742.0611.3041.812.3411.256.9273.0511.0040.9071.5801.442
Desviación Estándar0.0550.033628.057706.6980.2290.0921.330.534944.3760.0520.0000.2870.1190.358
PKPAKISTAN_SWAT_2683.711680.3281.322.3001.323.9603.6832.3205.183.9503.269.8703.3831.0050.9991.5871.585
PKPAKISTAN_SWAT_3683.621680.172113.239115.9893.7612.269453.346280.1223.4491.0050.9761.6581.618
PKPAKISTAN_SWAT_4683.763680.263529.699509.3634.3562.2422.456.0301.215.7803.5001.0051.0401.9432.020
PKPAKISTAN_SWAT_5683.776680.289416.365399.2353.8602.4511.710.9901.041.5803.4871.0051.0431.5751.643
PKPAKISTAN_SWAT_6683.750680.237521.774521.0443.6312.3992.016.8301.330.5903.5131.0051.0011.5141.516
PKPAKISTAN_SWAT_7683.737680.081462.937434.7533.7342.2431.840.1101.038.1503.6561.0051.0651.6651.772
PKPAKISTAN_SWAT_8683.801680.289507.196480.7043.5542.4251.918.6201.240.7603.5121.0051.0551.4661.546
PKPAKISTAN_SWAT_9683.763680.2891.351.1801.359.3703.5032.2945.037.6803.319.8903.4741.0050.9941.5271.517
Promedio683.741680.254638.093625.1393.7632.3612.522.8211.563.7783.4871.0051.0211.5941.613
Desviación Estándar0.0520.082416.646424.2050.2490.1211.566.3961.032.1430.0780.0000.0360.1490.162
ZMBZAMBIA_KAGEM_2_PER683.749680.340518.145481.8303.7852.6072.087.8301.337.2603.4091.0051.0751.4521.561
ZMBZAMBIA_KAGEM_3_PER683.749680.379680.223683.0403.5542.5032.573.5101.819.6403.3701.0050.9961.4201.414
ZMBZAMBIA_KAGEM_4_PER683.749680.327255.038252.2153.2972.033895.101545.8313.4221.0051.0111.6221.640
ZMBZAMBIA_KAGEM_6_PER683.749680.314598.185594.0653.3742.1642.148.5601.368.1603.4351.0051.0071.5601.570
ZMBZAMBIA_KAGEM_7_PER683.749680.4961.328.8601.185.7703.6312.4245.136.4803.059.0903.2531.0051.1211.4981.679
ZMBZAMBIA_KAGEM_8_PER683.749680.379847.645692.0793.7852.3723.415.5301.747.5803.3701.0051.2251.5961.954
ZMBZAMBIA_KAGEM_11_PER683.762680.483429.105415.8312.8091.7191.282.880761.0003.2791.0051.0321.6341.686
ZMBZAMBIA_KAGEM_13_PER683.697680.275809.476799.5573.6572.1383.151.5001.819.5003.4221.0051.0121.7111.732
ZMBZAMBIA_KAGEM_15_PER683.749680.431286.208251.2523.9652.8681.207.930767.0233.3181.0051.1391.3821.575
Promedio683.744680.386654.343609.2263.5092.2782.476.4361.485.9783.3591.0051.0741.5411.667
Desviación Estándar0.0200.080351.716310.3400.3540.3451.417.476819.9240.0690.0000.0840.1120.155
3.2.3. ANÁLISIS DE PCA
 

Con ayuda de un PCA, se identifican principalmente 2 grupos en función de la posición del pico R1 y la banda R3: el primero está conformado por las muestras colombianas CEOCC y CEOR y el segundo está conformado por las muestras de Pakistán y Zambia (Fig. 5). El componente principal de mayor contribución está relacionado con la posición del pico R1, aunque no es significativamente distinguible debido a que hay muestras de Colombia, Pakistán y Zambia que comparten una proporción en el eje de PC1. El segundo componente de mayor contribución está relacionado con la banda R3, en donde las muestras de Pakistán y Zambia presentan mayor intensidad.

media/001_005.jpg
Figura 5. Biplot de puntos (rangos en longitud de onda para R1 y R3) de los espectros tomados perpendicularmente. 
Figure 5. Point’s biplot of the spectra taken perpendicular to the C axis (ranges in wavelength for R1 and R3). 

Para el área (A), intensidad (I), posición del pico (R), distancia entre los picos (R1-R2) y ancho de altura medio (W) se realizó un PCA (Fig. 6). El PC1 tiene 61.8 % de contribución, además, permite distinguir significativamente las esmeraldas de Colombia frente a las de Pakistán y Zambia. Para las muestras de Colombia, las variables que más aportan son la intensidad y las áreas de los picos, mientras que para Pakistán y Zambia las variables que más aportan son: posición del pico R1, la altura media del pico, la relación entre intensidades, la distancia entre la posición de los picos y la relación de las áreas. El PC2 no es significativamente distinguible, ya que muestras de Colombia, Pakistán y Zambia comparten un rango en este eje.

media/001_006.jpg
Figura 6. Biplot con los datos obtenidos del tratamiento de los espectros de fotoluminiscencia orientados perpendicularmente. 
Figure 6. Biplot with the data obtained from the photoluminescence spectra (oriented perpendicularly). 

En el PCA con los datos de tratamiento estadístico orientado paralelamente (Fig. 7), se observa que el PC1 permite distinguir las muestras de Colombia de las de Zambia y Pakistán. Hay variables en los espectros de Zambia y Pakistán que tienen mucho aporte en esta distinción, estas son: el ancho de altura media de los dos picos (W), la relación de posición R1/R2, la distancia entre R1-R2 y la posición de R1. Para las esmeraldas colombianas, la variable que más aporta en la separación es la posición de R2. Hay otras variables que no aportan significativamente a la separación entre las muestras como son: I1, I2, Área R1, Área R2, I1/I2, W1/W2.

media/001_007.jpg
Figura 7. Biplot con los datos obtenidos del tratamiento de los espectros de fotoluminiscencia orientados paralelamente. 
Figure 7. Biplot with the data obtained from the photoluminescence spectra (oriented parallelly). 

3.3. Fluorescencia de rayos X (FRX)

 

Los resultados de FRX de los elementos mayores se recogen en la Tabla 4, mientras que los elementos menores se presentan en la Tabla 5. En promedio, las concentraciones de SiO₂ y Al₂O₃ constituyen entre 75 y 83 % de la composición total en todas las muestras, sin embargo, las esmeraldas de Colombia presentan valores en promedio más altos, cercanos a 81 %. El contenido de BeO (13,15 %) fue añadido manualmente como valor de referencia (Groat et al., 2008Groat, L. A., Giuliani, G., Marshall, D. D., & Turner, D. (2008). Emerald deposits and occurrences: A review. Ore Geology Reviews, 34(1-2), 87-112. 10.1016/j.oregeorev.2007.09.003.), con el fin de completar el cierre elemental. El 90 % de las muestras de Pakistán presentan concentraciones de Na2O y MgO con un rango que oscila entre 5 y 14 %. Sin embargo, la lectura de Na2O se debe tomar con precaución ya que el equipo de XRF, aunque lo detecta, puede presentar dificultades para cuantificar.

El PCA, con los datos de concentración química de las muestras estudiadas (Fig. 8), muestra que el PC1 es el de mayor contribución, con un 22 %, y permite diferenciar dos grupos de esmeraldas. El primero está formado por las muestras colombianas CEOR y CEOCC, en el cual los elementos de mayor aporte son los cromóforos Cr y V, además de Ga, Sc y K. El segundo grupo está conformado por las muestras de Pakistán y Zambia, en los cuales los elementos de mayor aporte son Ni, Ti, Fe y Zn. En cuanto al PC2, con una contribución de 15.6 %, permite distinguir entre los dos grupos anteriormente descritos, pero, además, entre ellos existe una subdivisión adicional. En el primer grupo es posible diferenciar entre las esmeraldas de Colombia CEOCC y CEOR, de tal forma que para el CEOCC los elementos de mayor aporte son el Cr, V, Sc y Ga, mientras que para CEOR son el Si, Cl y K. En el otro grupo el PC 2 permite diferenciar las muestras de Zambia, en las que los elementos de mayor aporte son Zn, Ni y la relación Cr/V que es mucho más alta que las otras áreas de estudio, mientras que en Pakistán el Ti, Fe y el Cu son los elementos de mayor aporte.

media/001_008.jpeg
Figura 8. Biplot de concentraciones de elementos químicos que componen las muestras de estudio (elementos mayores, menores y traza). 
Figure 8. Chemical concentration biplot of elements that compound the samples (major, minor, and trace elements). 

Por otra parte, utilizando los resultados de fotoluminiscencia, se presentan los datos de tratamiento de los espectros y la concentración de elementos químicos con respecto a su orientación. Con la orientación perpendicular (Fig. 9), se observa que el PC1 permite distinguir las muestras de origen colombiano de las no colombianas, utilizando parámetros como R1, la relación I1/I2 y la presencia de elementos como Ti, Fe y K, entre otros. No obstante, existen variables que contribuyen a ambos componentes principales y que también ayudan a diferenciar poblaciones en el PC2. En este contexto, elementos químicos como Cr, V, Sc y Ga son útiles para distinguir entre el CEOCC y el CEOR. Además, al utilizar el área de R2 e indicadores como I1 e I2 se facilita esta diferenciación. Asimismo, la separación entre las muestras de Zambia y Pakistán se ve facilitada por la presencia de elementos como Ti, Ni y Zn y la relación I1/I2.

media/001_009.jpg
Figura 9. Biplot de concentración química y datos de tratamiento de espectros orientados perpendicularmente. 
Figure 9. Biplot compiling chemistry and photoluminescence data (oriented perpendicularly). 

Con orientación paralela al eje C y la composición química (Fig. 10), se pueden distinguir las esmeraldas de Colombia de las de Pakistán y Zambia, es decir, se forman dos grupos: el primero, con las esmeraldas de Colombia y el segundo con las muestras no colombianas. Para el grupo de Colombia, las variables de mayor aporte en el PC1 son la posición de R2, la concentración de K, SiO2 y Cl. Para las esmeraldas de Zambia y Pakistán las variables que más aportan en el PC1 son el ancho de altura media de los dos picos, la relación de R1/R2, la distancia entre R1 y R2 y la posición de R1.

media/001_010.jpeg
Figura 10. Biplot con datos de tratamiento de espectros tomados paralelamente y concentración química. 
Figure 10. Biplot compiling chemistry and photoluminescence data (oriented parallelly). 

4. Discusión

 

4.1. Espectroscopía de fotoluminiscencia

 

La posición de los picos R1 y R2 son parámetros que permiten no solo identificar la meta-estabilidad del Cr (Wood, 1965Wood, D. L. (1965). Absorption, fluorescence, and Zeeman effect in Emerald. The Journal of Chemical Physics, 42(10), 3404-3410. 10.1063/1.1695742.) sino observar diferencias entre muestras. Los resultados por fotoluminiscencia muestran que la posición R1 diferencia las muestras colombianas de Zambia y Pakistán por encontrarse en longitudes de onda menores, con promedio de 683.650 nm, rangos que contrastan con Zambia y Pakistán que se encuentran en longitudes de onda mayores, por encima de 684.000 nm (Moroz et al., 2000Moroz, I., Roth, M., Boudeulle, M., & Panczer, G. (2000). Raman Microspectroscopy and Flourescence of Emeralds from Various Deposits. Journal of Raman Spectroscopy, 31, 485-490.; Thompson et al., 2014Thompson, D. B., Kidd, J. D., Åström, M., Scarani, A., & Smith, C. P. (2014). A Comparison of R-line Photoluminescence of Emeralds from Different Origins. The Journal of Gemmology, 34(4), 334-343. 10.15506/jog.2014.34.4.334., 2017Thompson, D. B., Bayens, C. J., Morgan, M. B., Myrick, T. J., & Sims, N. E. (2017). Photoluminescence spectra of emeralds from Colombia, Afghanistan, and Zambia. Gems & Gemology, 53(3), 296-311. 10.5741/GEMS.53.3.296.; Schmetzer, 2015Schmetzer, K. (2015). Photoluminescence of emeralds: Sample orientation procedure and correlation of the R1 peak position with SiO2 contents. Journal of Gemmology, 34, (5), 441-443.; Karampelas et al., 2019Karampelas, S., Al-Alawi, A., Al-Shaybani, B., Mohamed, F., & Sangsawong, S. (2019). Emeralds from the most important occurrences: Chemical and spectroscopic data. Minerals, 9(9), 561. 10.3390/min9090561.). Además, la posición de R1 puede llegar a ser un parámetro para la distinción entre CEOR y CEOCC en Colombia, debido a que el CEOCC presenta valores de longitud de onda en general mayores, y valores de intensidad menores que los de CEOR. Además, la presencia de la banda R3 se manifiesta como una disminución en la intensidad entre las longitudes de onda entre los picos R1 y R2, de modo que este parámetro podría convertirse en el futuro en otro criterio útil para diferenciar entre las esmeraldas colombianas y las extranjeras. Cabe destacar que esta caída de energía aún no ha sido documentada en investigaciones previas (Moroz et al., 2000Moroz, I., Roth, M., Boudeulle, M., & Panczer, G. (2000). Raman Microspectroscopy and Flourescence of Emeralds from Various Deposits. Journal of Raman Spectroscopy, 31, 485-490.; Thompson et al., 2014Thompson, D. B., Kidd, J. D., Åström, M., Scarani, A., & Smith, C. P. (2014). A Comparison of R-line Photoluminescence of Emeralds from Different Origins. The Journal of Gemmology, 34(4), 334-343. 10.15506/jog.2014.34.4.334., 2017Thompson, D. B., Bayens, C. J., Morgan, M. B., Myrick, T. J., & Sims, N. E. (2017). Photoluminescence spectra of emeralds from Colombia, Afghanistan, and Zambia. Gems & Gemology, 53(3), 296-311. 10.5741/GEMS.53.3.296.; Schmetzer, 2015Schmetzer, K. (2015). Photoluminescence of emeralds: Sample orientation procedure and correlation of the R1 peak position with SiO2 contents. Journal of Gemmology, 34, (5), 441-443.; Karampelas et al., 2019Karampelas, S., Al-Alawi, A., Al-Shaybani, B., Mohamed, F., & Sangsawong, S. (2019). Emeralds from the most important occurrences: Chemical and spectroscopic data. Minerals, 9(9), 561. 10.3390/min9090561.) sin embargo, no se abordará en este documento.

Los resultados geoquímicos de este estudio fueron comparados con algunos estudios previos en las regiones de Pakistán y Zambia. Para las esmeraldas de Swat Valley y Mohmand (Pakistán), Kazmi & Snee (1989Kazmi, A. H., & Snee, L. W. (1989). Emeralds of Pakistan: Geology, Gemology, and Genesis. Geological Survey of Pakistan, Professional Papers.) validan la presencia de minerales y elementos traza similares a los de este estudio, destacando altas concentraciones de Fe y Mg como indicadores clave para diferenciar entre esmeraldas colombianas y pakistaníes. Además, los autores también reportaron unas altas concentraciones de elementos cromóforos (Cr y V) para las esmeraldas de esta región. Es importante destacar que la presencia de Cr y V en la estructura del berilo causa una típica fluorescencia roja (Fig. 2, Tabla 1), pero si el Fe3+ está presente en el cristal, este efecto de fluorescencia es suprimido debido a un fenómeno llamado quenching (Nassau, 1978Nassau, K. (1978). The origins of color in minerals. American Mineralogist, 63(3-4), 219-229.; Kane & Liddicoat, 1985Kane, R. E., & Liddicoat, T. (1985). The Biron hydrothermal synthetic emerald. Gems & Gemology, 21(3), 156-70.; Groat et al., 2008Groat, L. A., Giuliani, G., Marshall, D. D., & Turner, D. (2008). Emerald deposits and occurrences: A review. Ore Geology Reviews, 34(1-2), 87-112. 10.1016/j.oregeorev.2007.09.003.).

Para las esmeraldas de Zambia, los estudios de Zwaan et al. (2005Zwaan, J. C., Seifert, A. V., Vrána, S., Laurs, B. M., Anckar, B., Simmons, W. B., Falster, A. U., Lustenhouwer, W. J., Muhlmeister, S., Koivula, J. I., & Garcia-Guillerminet, H. (2005). Emeralds from the Kafubu area, Zambia. Gems & Gemology, 41(2), 116-148. 10.5741/GEMS.41.2.116.) y Giuliani et al. (2019) proporcionaron una gran cantidad de datos sobre la composición geoquímica y las condiciones geológicas del área de Kafubu. Adicionalmente, Zwaan et al. (2005Zwaan, J. C., Seifert, A. V., Vrána, S., Laurs, B. M., Anckar, B., Simmons, W. B., Falster, A. U., Lustenhouwer, W. J., Muhlmeister, S., Koivula, J. I., & Garcia-Guillerminet, H. (2005). Emeralds from the Kafubu area, Zambia. Gems & Gemology, 41(2), 116-148. 10.5741/GEMS.41.2.116.) documentaron las características geoquímicas específicas, como la relación entre el Cr y el V, que concuerdan con las presentadas para este estudio.

El área y HWFM de los picos en los espectros adquiridos en ambas orientaciones pueden ser usados como parámetros adicionales, con el propósito de distinguir las muestras de Zambia y Pakistán que se caracterizan por tener mayores HWFM que las muestras colombianas.

La relación entre las intensidades I1/I2 en los picos representa otro parámetro diferenciador de las esmeraldas colombianas, con promedio de 3.2, frente a las muestras de Zambia y Pakistán, que presentan promedios cercanos a 5.2. Esto coincide con el trabajo de Thompson et al. (2014Thompson, D. B., Kidd, J. D., Åström, M., Scarani, A., & Smith, C. P. (2014). A Comparison of R-line Photoluminescence of Emeralds from Different Origins. The Journal of Gemmology, 34(4), 334-343. 10.15506/jog.2014.34.4.334.), quienes observaron que esta relación es mayor en esmeraldas tipo esquisto, como las muestras de Zambia, con promedio de 3.4 aproximadamente, mientras que las colombianas tienen una relación de I1/I2 de 2.85, ambos valores diferentes al presente estudio. Sin embargo, estas diferencias se deben a que el estudio de Thompson et al. (2014Thompson, D. B., Kidd, J. D., Åström, M., Scarani, A., & Smith, C. P. (2014). A Comparison of R-line Photoluminescence of Emeralds from Different Origins. The Journal of Gemmology, 34(4), 334-343. 10.15506/jog.2014.34.4.334.) solo emplea esmeraldas del CEOOC. En esmeraldas de este cinturón (CEOCC) la relación I1/I2 disminuye hacia un promedio de 2.9 a 3 (muy cercano al obtenido por estos autores). Por otra parte, Thompson et al. (2017Thompson, D. B., Bayens, C. J., Morgan, M. B., Myrick, T. J., & Sims, N. E. (2017). Photoluminescence spectra of emeralds from Colombia, Afghanistan, and Zambia. Gems & Gemology, 53(3), 296-311. 10.5741/GEMS.53.3.296.) demuestra que el ángulo de incidencia del láser con respecto al eje C del cristal, repercute en un aumento en el valor de la relación I1/I2, llegando a obtener valores por encima de 4.5, muy cercanos a los presentados en este artículo (5.2) para las esmeraldas de Pakistán y Zambia. Las muestras de Colombia en general presentan una mayor intensidad frente a las muestras de Zambia (Thompson et al., 2017Thompson, D. B., Bayens, C. J., Morgan, M. B., Myrick, T. J., & Sims, N. E. (2017). Photoluminescence spectra of emeralds from Colombia, Afghanistan, and Zambia. Gems & Gemology, 53(3), 296-311. 10.5741/GEMS.53.3.296.) y Pakistán, llegando a tener el doble o hasta 5 veces más intensidad. Finalmente, la relación entre las áreas de los picos en orientación perpendicular, donde las muestras de Zambia y Pakistán, con promedio de 8.3, son mayores que las de Colombia, que tienen un promedio de 4.3, permite identificar otro parámetro diferenciador.

4.2. Composición química de las esmeraldas

 

Al comparar las concentraciones promedio de las esmeraldas de Colombia (CEOR y CEOCC), Zambia y Pakistán, se evidencia una clara diferencia entre las concentraciones de los elementos cromóforos; por ejemplo, la baja concentración de Fe de las esmeraldas colombianas (Cedeño et al., 2015Cedeño, C. J., Jiménez, J. F., Herreño, J. H., & Fortaleché, D. (2015). Progress on the study of parameters related to the origin of Colombian emeralds. Color, World Emer, (Special Issue), 88-97.) frente a las esmeraldas de Zambia y Pakistán, que llegan a tener hasta 10 veces más concentración de Fe que las colombianas. Además, las muestras de Pakistán son las de mayor concentración de todas las áreas de estudio, con un promedio de 23000 ppm de Fe.

La relación Cr/V difiere entre las áreas de estudio. Las esmeraldas de Zambia tienen una relación promedio de 16, debido a las concentraciones muy bajas de V, cercanas a 250 ppm, mientras que los contenidos de Cr están alrededor de 3000 ppm. Por su parte, en Pakistán, aunque la relación es 0.33, lo cual podría indicar un enriquecimiento en V, las muestras presentan una concentración baja de ambos elementos (Cr y V), siendo ligeramente más alta en promedio la de V. Las esmeraldas de Colombia tienen unas concentraciones altas de los dos elementos Cr y V, en un rango de 500 a 9000 ppm, y presentan una relación en el CEOCC cercana a 1 y en el CEOR cercana a 2.5-3, que está acorde con los valores presentados por Cedeño et al. (2015Cedeño, C. J., Jiménez, J. F., Herreño, J. H., & Fortaleché, D. (2015). Progress on the study of parameters related to the origin of Colombian emeralds. Color, World Emer, (Special Issue), 88-97.).

Además, las muestras de CEOCC presentan una fuerte correlación positiva entre los elementos Cr, V, Ga y Sc, y esta agrupación presenta una correlación negativa con muestras de Pakistán, que tienen una fuerte correlación de los elementos Fe, Cu y Ti (Fig. 8). El Fe ha sido identificado como un parámetro diferenciador por su alta concentración en esmeraldas tipo esquisto (Krzemnicki et al., 2021Krzemnicki, M. S., Wang, H. A. O., and Büche, S. (2021). A new type of emerald from Afghanistan’s Panjshir Valley. The Journal of Gemmology, 37(5), 474-495. 10.15506/JoG.2021.37.5.474.). Además, contrasta con la alta concentración de Cr y V para esmeraldas de Colombia, hecho que ya ha sido identificado en estudios previos (Groat et al., 2008Groat, L. A., Giuliani, G., Marshall, D. D., & Turner, D. (2008). Emerald deposits and occurrences: A review. Ore Geology Reviews, 34(1-2), 87-112. 10.1016/j.oregeorev.2007.09.003.; Cedeño et al., 2015Cedeño, C. J., Jiménez, J. F., Herreño, J. H., & Fortaleché, D. (2015). Progress on the study of parameters related to the origin of Colombian emeralds. Color, World Emer, (Special Issue), 88-97.; Karampelas et al., 2019Karampelas, S., Al-Alawi, A., Al-Shaybani, B., Mohamed, F., & Sangsawong, S. (2019). Emeralds from the most important occurrences: Chemical and spectroscopic data. Minerals, 9(9), 561. 10.3390/min9090561.; Alonso-Perez & Day, 2021Alonso-Perez, R., & Day, J. M. D. (2021). Rare Earth Element and Incompatible Trace Element Abundances in Emeralds Reveal Their Formation Environments. Minerals, 11(5), 513. 10.3390/min11050513.). Las esmeraldas de CEOR presentan una fuerte correlación positiva entre los elementos K, Cl y Si. Además, esta agrupación presenta una correlación negativa con muestras de Zambia, cuyas variables de mayor aporte son Zn y Ni y la relación Cr/V. En cuanto a la relación Cr/V, de acuerdo con Karampelas et al. (2019Karampelas, S., Al-Alawi, A., Al-Shaybani, B., Mohamed, F., & Sangsawong, S. (2019). Emeralds from the most important occurrences: Chemical and spectroscopic data. Minerals, 9(9), 561. 10.3390/min9090561.), hay una relación más alta para Zambia que se puede deber a la baja concentración de V. Por otro lado, el contenido de K representa una de las variables que más aportan al CEOR; sin embargo, según lo reportado por Abduriyim & Kitawaki, (2006Abduriyim, A., & Kitawaki, H. (2006). Applications of laser ablation – inductively coupled plasma – mass spectrometry (LA-ICP-MS) to gemology. Gems and Gemology, 42(2), 98-118.) y Saeseaw et al. (2019Saeseaw, S., Renfro, N. D., Palke, A. C., Sun, Z., & McClure, S. F. (2019). Geographic Origin Determination of Emerald. Gems & Gemology, 55(4), 614-646. 10.5741/GEMS.55.4.614.), el contenido de K es mayor en esmeraldas de Zambia, contrario a los datos obtenidos en este estudio, aunque es necesario tener en cuenta que el K puede presentar problemas para su cuantificación con el uso de XRF, debido a que tiene una baja eficiencia de excitación en comparación con otros elementos.

5. Conclusiones

 

Las esmeraldas colombianas reaccionan a la luz ultravioleta en mayor proporción que las esmeraldas extranjeras debido a la cantidad de Cr presente en las muestras y menor cantidad de Fe, que genera en ocasiones el efecto quenching (Nassau, 1978Nassau, K. (1978). The origins of color in minerals. American Mineralogist, 63(3-4), 219-229.; Kane & Liddicoat, 1985Kane, R. E., & Liddicoat, T. (1985). The Biron hydrothermal synthetic emerald. Gems & Gemology, 21(3), 156-70.), razón por la cual las esmeraldas del CEOR presentan colores más intensos en presencia de esta luz. El análisis de dicha intensidad debería ser objeto de investigación en estudios posteriores, mediante espectroscopia UV-Vis-NIR. Las esmeraldas del CEOCC presentan una reacción intermedia, y finalmente las muestras de Zambia y Pakistán no reaccionan o reaccionan muy poco cuando son sometidas a la luz ultravioleta.

El uso de la fotoluminiscencia permite distinguir entre esmeraldas colombianas y esmeraldas de Pakistán y Zambia debido a las diferencias en la longitud de onda presentes en la posición de los picos R1 y R2, además de diferentes parámetros como son la intensidad de los picos, área de los picos, ancho medio de los picos, intensidad de R3, la distancia entre R1 y R2 y algunas relaciones entre estas variables como R1-R2, I1/I2, etc. Además, permite pensar que la fotoluminiscencia se presenta como una posible herramienta que permitirá diferenciar entre ambos cinturones esmeraldíferos colombianos, CEOR y CEOCC.

Debido a las diferencias presentes en los ambientes de formación que imprimen una huella geoquímica distintiva para cada uno de los depósitos de esmeraldas, es posible realizar distinciones a partir de sus concentraciones químicas. Un caso son los bajos valores de Fe, Ti, Cu, Zn y Ni, en Colombia, mientras que existe una mayor concentración de estos elementos para esmeraldas provenientes de Pakistán y Zambia. Entre las esmeraldas colombianas, la presencia de K, Cl y Sc, además de la relación entre los elementos cromóforos (Cr/V), permite diferenciar muestras del CEOR y el CEOCC.

Es importante resaltar que, en algunas ocasiones, las esmeraldas pueden traslaparse en algunos de los parámetros utilizados para diferenciarlas. Por ello, se recomienda implementar todas las herramientas disponibles para su correcta diferenciación, incluyendo análisis de fluorescencia, fotoluminiscencia y composición química.

Agradecimientos

 

Agradecemos a la Federación Nacional de Esmeraldas, FEDESMERALDAS, junto con las asociaciones ACODES, APRECOL, ASOCOESMERALDA y al Ministerio de Minas y Energía, MINENERGIA, por permitir y apoyar la investigación de las esmeraldas colombianas de las que se derivan trabajos de investigación como el presente estudio.

Referencias

 

1 

Abduriyim, A., & Kitawaki, H. (2006). Applications of laser ablation – inductively coupled plasma – mass spectrometry (LA-ICP-MS) to gemology. Gems and Gemology, 42(2), 98-118.

2 

Alonso-Perez, R., & Day, J. M. D. (2021). Rare Earth Element and Incompatible Trace Element Abundances in Emeralds Reveal Their Formation Environments. Minerals, 11(5), 513. https://doi.org/10.3390/min11050513.

3 

Alonso-Perez, R., Day, J. M., Pearson, D. G., Luo, Y., Palacios, M. A.,Sudhakar, R., & Palke, A. (2024). Exploring emerald global geochemical provenance through fingerprinting and machine learning methods. Artificial Intelligence in Geosciences, 5, p. 100085. https://doi.org/10.1016/j.aiig.2024.100085.

4 

Araújo Neto, J. F. de, Müller, A., Barreto, S. de B., Santos, L. C. M. de L., and Carrino, T. A. (2019). Mineralogical and gemological characterization of emerald crystals from Paraná deposit, NE Brazil: a study of mineral chemistry, absorption and reflectance spectroscopy and thermal analysis. Brazilian Journal of Geology, 49(3), e20190014. https://doi.org/10.1590/2317-4889201920190014.

5 

Banks, D. A., Giuliani, G., Yardley, B. W. D., & Cheilletz, A. (2000). Emerald mineralisation in Colombia: Fluid chemistry and the role of brine mixing. Mineralium Deposita, 35, 699-713. https://doi.org/ 10.1007/s001260050273.

6 

Bersani, D., Azzi, G., Lambruschi, E., Barone, G., Mazzoleni, P., Raneri, S., & Lottici, P. P. (2014). Characterization of emeralds by micro-Raman spectroscopy. Journal of Raman Spectroscopy, 45, 11-12. https://doi.org/10.1002/jrs.4524.

7 

Breeding, C. M., & Shigley, J. E. (2009). The “type” classification system of diamonds and its importance in gemology. Gems & Gemology, 45(2), 96-111. https://doi.org/10.5741/GEMS.45.2.96.

8 

Cedeño, C. J., Jiménez, J. F., Herreño, J. H., & Fortaleché, D. (2015). Progress on the study of parameters related to the origin of Colombian emeralds. Color, World Emer, (Special Issue), 88-97.

9 

Giuliani, G., & Groat, L. A. (2019). Geology of corundum and emerald gem deposits: a review. Gems and Gemology, 55(4), 464-489. https://doi.org/10.5741/GEMS.55.4.464.

10 

Giuliani, G., Silva, L. J. H. D., & Couto, P. (1990). Origin of emerald deposits of Brazil. Mineralium Deposita, 25, 57-64. https://doi.org/10.1007/BF03326384.

11 

Giuliani, G., Cheilletz, A., Sheppard, S. M. F., & Arboleda, C. (1993). Geochemistry and origin of the emerald deposits of Colombia. En Fenoll Hach-Ali, Torres-Ruiz y Gervilla (Eds.), Current research in geology applied to ore deposits (pp. 105-108). Biennial SGA Meeting, 2, Granada (España), 1993/09/09-11. ISBN 84-338-1772-8.

12 

Giuliani, G., Cheilletz, A., Arboleda, C., Carrillo, V., Rueda, F., & Baker, J. H. (1995). An evaporitic origin of the parent brines of Colombian emeralds: fluid inclusion and sulphur isotope evidence. European Journal of Mineralogy, 7(1), 151-165.

13 

Giuliani, G., Groat, L. A., Marshall, D., Fallick, A. E., & Branquet, Y. (2019). Emerald deposits: A review and enhanced classification. Minerals, 9(2), 105. https://doi.org/10.3390/min9020105.

14 

Groat, L. A., Giuliani, G., Marshall, D. D., & Turner, D. (2008). Emerald deposits and occurrences: A review. Ore Geology Reviews, 34(1-2), 87-112. https://doi.org/10.1016/j.oregeorev.2007.09.003.

15 

Kane, R. E., & Liddicoat, T. (1985). The Biron hydrothermal synthetic emerald. Gems & Gemology, 21(3), 156-70.

16 

Karampelas, S., Al-Alawi, A., Al-Shaybani, B., Mohamed, F., & Sangsawong, S. (2019). Emeralds from the most important occurrences: Chemical and spectroscopic data. Minerals, 9(9), 561. https://doi.org/10.3390/min9090561.

17 

Kazmi, A. H., & Snee, L. W. (1989). Emeralds of Pakistan: Geology, Gemology, and Genesis. Geological Survey of Pakistan, Professional Papers.

18 

Kiefert, L., & Karampelas, S. (2011). Use of the Raman spectrometer in gemmological laboratories: Review. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 80, 119-124. https://doi.org/10.1016/j.saa.2011.03.004.

19 

Krzemnicki, M. S., Wang, H. A. O., and Büche, S. (2021). A new type of emerald from Afghanistan’s Panjshir Valley. The Journal of Gemmology, 37(5), 474-495. https://doi.org/10.15506/JoG.2021.37.5.474.

20 

Lim, H., Park, S., Cheong, H., Choi, H.-M., & Kim, Y. C. (2010). Discrimination between natural and HPHT-treated type IIa diamonds using photoluminescence spectroscopy. Diamond and Related Materials, 19(10), 1254-1258.

21 

Moroz, I., Roth, M., Boudeulle, M., & Panczer, G. (2000). Raman Microspectroscopy and Flourescence of Emeralds from Various Deposits. Journal of Raman Spectroscopy, 31, 485-490.

22 

Nassau, K. (1978). The origins of color in minerals. American Mineralogist, 63(3-4), 219-229.

23 

Palke A., Saeseaw, S., Renfro N. D., Sun Z., & McClure, F. (2019). Geographic origin determination of blue sapphire. Gems & Gemology, 55(4), 536-579. https://doi.org/10.5741/GEMS.55.4.536.

24 

Saeseaw S., Pardieu V., & Sangsawong S. (2014). Three-Phase Inclusions in Emerald and Their Impact on Origin Determination. Gems & Gemology, 50(2), 114-132. https://doi.org/10.5741/GEMS.50.2.114.

25 

Saeseaw, S., Renfro, N. D., Palke, A. C., Sun, Z., & McClure, S. F. (2019). Geographic Origin Determination of Emerald. Gems & Gemology, 55(4), 614-646. https://doi.org/10.5741/GEMS.55.4.614.

26 

Schmetzer, K. (2015). Photoluminescence of emeralds: Sample orientation procedure and correlation of the R1 peak position with SiO2 contents. Journal of Gemmology, 34, (5), 441-443.

27 

Thompson, D. B., Kidd, J. D., Åström, M., Scarani, A., & Smith, C. P. (2014). A Comparison of R-line Photoluminescence of Emeralds from Different Origins. The Journal of Gemmology, 34(4), 334-343. https://doi.org/10.15506/jog.2014.34.4.334.

28 

Thompson, D. B., Bayens, C. J., Morgan, M. B., Myrick, T. J., & Sims, N. E. (2017). Photoluminescence spectra of emeralds from Colombia, Afghanistan, and Zambia. Gems & Gemology, 53(3), 296-311. https://doi.org/10.5741/GEMS.53.3.296.

29 

Wood, D. L. (1965). Absorption, fluorescence, and Zeeman effect in Emerald. The Journal of Chemical Physics, 42(10), 3404-3410. https://doi.org/10.1063/1.1695742.

30 

Wotherspoon, A., Steeds, J. W., Catmull, B., & Butler, J. (2003). Photoluminescence and positron annihilation measurements of nitrogen doped CVD diamond. Diamond and Related Materials, 12, 652-657. https://doi.org/10.1016/S0925-9635(02)00229-7.

31 

Zwaan, J. C., Seifert, A. V., Vrána, S., Laurs, B. M., Anckar, B., Simmons, W. B., Falster, A. U., Lustenhouwer, W. J., Muhlmeister, S., Koivula, J. I., & Garcia-Guillerminet, H. (2005). Emeralds from the Kafubu area, Zambia. Gems & Gemology, 41(2), 116-148. https://doi.org/10.5741/GEMS.41.2.116.