Acid-Gangue Interactions in Heap Leach Operations
Acid-Gangue Interactions in Heap Leach Operations
Acid-Gangue Interactions in Heap Leach Operations
Resumen: La lixiviación en pilas representa una quinta parte de la producción mundial de cobre, proveniente
principalmente de minerales pórfidos, sin embargo, las recuperaciones de metales a menudo no son óptimas.
La ganga, y su interacción con el ácido, juega un papel importante en tales procesos. Por lo tanto, es necesario
comprender adecuadamente los minerales de ganga presentes en el mineral, sus relaciones texturales
relativas a la distribución del tamaño de partícula, reactividad con ácido bajo diferentes condiciones y
relación con litotipos y alteración geológica en el yacimiento, para predecir el comportamiento del mineral en
la trituración, operaciones de aglomeración, curado y lixiviación en pilas. Las herramientas mineralógicas
disponibles para la caracterización son la difracción de rayos X de rutina, la microscopía óptica, la
microscopía electrónica de barrido automatizada y el microanálisis de la sonda de electrones, junto con los
avances más recientes en imágenes infrarrojas hiperespectrales y tomografía computarizada de rayos X. El
uso integrado de estas técnicas permite abordar la abundancia de minerales, las relaciones de textura y la
química mineral en todo el rango de tamaños de partículas y aglomerados. Además, los resultados de
lixiviación de diagnóstico se pueden interpretar mejor cuando se calibran con datos mineralógicos sólidos. La
vinculación de los atributos del mineral, el comportamiento metalúrgico y su distribución en el yacimiento
forma una parte integral de un enfoque geometalúrgico para predecir y abordar los cambios durante el
proceso de lixiviación en pilas. Investigaciones adicionales deberían abordar los fundamentos de la reacción
de ganga con ácido fuerte y la degradación estructural concomitante durante los procesos de curado y
aglomeración, y cómo esto difiere de la reactividad de ácido de ganga en condiciones ácidas más débiles,
combinado con la temperatura y los efectos del flujo de fluido de la lixiviación en pilas. La caracterización
previa y posterior es necesaria para comprender y cuantificar los efectos de las variables para la reactividad
del ácido de ganga en estas diversas operaciones. Los resultados de la caracterización deberían conducir a un
refinamiento de la jerarquía de la reactividad mineral de la ganga en diferentes condiciones de conminución y
lixiviación, y cómo estos a su vez se relacionan con los tipos de mineral identificados que se extraen.
Keywords: reactividad al ácido de ganga; mineralogía; lixiviación en pilas; curación; aglomeración; textura.
1. Introducción
Dada la disminución en las leyes del mineral con un incremento proporcional en la complejidad para el
procesamiento, la lixiviación en pilas ha demostrado, en las últimas décadas, que es un medio
económicamente viable para extraer materias primas como cobre, zinc, níquel y oro [1]. En la actualidad, los
minerales pórfidos de cobre representan la mayor parte de la producción de cobre, y la lixiviación en pilas
representa ~ 21% de la producción mundial de cobre [2]. A pesar de esto, las recuperaciones a menudo no son
óptimas [1]. Para comprender los procesos clave que tienen lugar durante las operaciones unitarias
involucradas en la tecnología de lixiviación en pilas, que incluyen trituración, aglomeración y curado, y la
lixiviación en pilas, la mineralogía es fundamental y ayuda a predecir el comportamiento del mineral durante
estas operaciones unitarias [3-5].
A este respecto, no solo los minerales del mineral deberían generar interés para la recuperación del metal,
sino que es importante entender la mineralogía de la ganga, particularmente desde la perspectiva de
Tabla 1. Importantes minerales de silicato presentes en los depósitos de pórfidos de cobre [8] .
Figura 1. Minerales formadores de rocas comunes ordenados por cristalización y cinética de disolución en
soluciones de lixiviación de sulfato ácido [8].
Las estructuras de silicato se pueden dividir en seis tipos (Tabla 2). Los orto-, ciclo- y sorosilicatos son en
gran parte solubles, dando como resultado una disolución congruente, es decir, disolución estequiométrica
para liberar cationes y sílice simultáneamente. Dichas reacciones pueden dar como resultado la formación de
gel de sílice, debido a la polimerización rápida de la sílice en solución. La formación de gel puede causar
ceguera y evitar la percolación, por lo tanto, es perjudicial en las operaciones de lixiviación en pilas [14,15]. Los
ino-, filo y tectosilicatos son en gran parte insolubles, pero con el tiempo tiene lugar una disolución
incongruente que da como resultado la liberación del catión metálico y un residuo silíceo restante, es decir, la
disolución no es estequiométrica [16].
Ti 12,058 Mn 3745
4+ 2+
Zr 11,033 Cu 3598
4+ 2+
Ca 3510
2+
Al 7201–7858 Pb 3469
3+ 2+
Y3 6552 Ba 3213
+ 2+
Be 4774 Li 1469
+
2+
Zn 3937 N 1347
2+ a+
Ni 3887 K 1251
2+ +
Fe 3845
2+
Otro ejemplo de importancia práctica en la lixiviación en pilas es la medida en que Al forma parte de la
estructura en los tectosilicatos. Los tectosilicatos que contienen Al en la relación de Al: Si de al menos 2: 3 serán
más reactivos y susceptibles a la gelatinización que aquellos con relaciones Al: Si más bajas, y explican una
mayor reactividad de anortita (rica en calcio) sobre la albita (rica en sodio) plagioclasa. De manera similar, los
filosilicatos que contienen hierro férrico apreciable en las láminas (por ejemplo, nontronita) serán susceptibles al
ataque ácido y la gelatinización [15,19].
Además de esto, el trabajo de Free [20] mostró que la biotita consume ácido más rápidamente a un pH más bajo
en comparación con el clorito, y los estudios exhaustivos de Jansen y Taylor [21] mostraron que el clorito es un
consumidor significativo de ácido a largo plazo. Dado que la biotita es el consumidor ácido más reactivo en
ácidos fuertes, clorita fue clasificada como la segunda en experimentos con titanita, ilmenita, epidota,
hornblenda, moscovita, plagioclasa, feldespato potásico, rutilo y anatasa [18], por lo que la disolución de clorita
puede contribuir al ácido consumo en lixiviación en pilas Lo mismo parece ser cierto para el epidoto en la
lixiviación a largo plazo [21]. En base a las observaciones de los diversos trabajadores anteriores, la Tabla 4
enumera el orden anticipado de reactividad para los minerales puros seleccionados, así como las reacciones
involucradas. Vale la pena señalar que los minerales de arcilla, particularmente los del grupo de esmectita, ya
sea como productos de alteración del mineral o productos de la reacción de álabe en el montón, tienen el
potencial de mantener el ácido en sus estructuras [22,23]. Si bien no es estrictamente el consumo de ácido, la
capacidad de adsorción de ácido hace que estos minerales sean importantes ya que el ácido no está disponible
para reaccionar con los minerales de Cu. La disolución de la ganga durante la interacción ácida conducirá a una
sucesión de residuos y productos precipitados, con la liberación concomitante de elementos tales como Ca, Mg,
Na, K y sílice soluble. Se espera que los productos finales sean óxidos hidratados y cuarzo secundario [24,25].
Por lo tanto, el conocimiento del ensamblaje mineral inicial, la reactividad con el ácido y los productos
potenciales que se formarán, son los primeros parámetros importantes para predecir la ganga.
comportamiento mineral en lixiviación.
Tabla 4. Reactividad relativa propuesta de los silicatos de ganga importantes encontrados comúnmente en
los depósitos de pórfidos de cobre, en condiciones de lixiviación en pilas .
Mineral Solubility
Readily
Soluble:
#biotite: (H,K)2(Mg,Fe)2Al2(SiO4)3 + 10H+ = 2K+ + 2Al+++ + 3H4SiO4 + 2(Fe,Mg)++
Moderately soluble:
chlorite: (Fe,Mg)3[(Si,Al)2O5](OH)4 + 6H+ = 3(Mg,Fe)++ + H4SiO4 + Al+++ + H2O
epidote: Ca2Al2O.(Al,Fe)OH[Si2O7][SiO4] + 13H+ = 2Ca++ + 2Al+++ + (Al,Fe)+++ + 3H4SiO4 + H2O
#pyroxe Ca(Mg,Fe)(SiO3)2 + 4H+ + 2H2O = Ca++ + 2H4SiO4 + Mg++
ne:
amphibo Ca2(Mg,Fe)4Al[Si7AlO22](OH)2 + 18H+ + 4H2O = 2Ca++ + 4Mg++ + 2Al+++ + 7H4SiO4
le:
#anorthi CaAl2Si2O8 + 10H+ = 2Ca++ + 2Al+++ + 2H4SiO4
te:
Poorly soluble:
Muscovite, albite, K-feldspar; quartz, clays *, anhydrite
# Reactions from [21]; * note potential to hold acid through absorption; + mono-, ++ di- and +++ trivalent
cations.
Distribución del tamaño del grano mineral y relación con la distribución del tamaño de partícula y el
modo de trituración
Cristalidad mineral
Partículas aglomeradas vs. sueltas
Resistencia ácida
Porosidad
Temperatura
Para el tamaño del grano mineral, como es intuitivo, cuanto menor es el tamaño del grano del mineral,
mayor es la reactividad esperada. El trabajo de Snäll y Liljefors [18] mostró diferencias en la lixiviabilidad de
minerales para tamaños de grano de 2 mm frente a polvos en minerales separados, y concluyó que para
minerales menos solubles en ácido, el tamaño de grano desempeñó un papel importante para la reactividad.
Los minerales lixiviados en el orden de más fácil a más difícil fueron: biotita> clorito de Mg> esfeno> ilmenita>
epidota> hornblenda> moscovita> plagioclasa (Ca)> albita / oligoclasa> microclina> rutilo> anatasa, para un
tamaño de grano dado. Conclusiones similares sobre el tamaño del grano mineral y la reactividad con el ácido
fueron alcanzadas por Free [20]. Además, Terry [14] mostró que el clorito de grano fino reaccionó más rápido
en HCl que el clorito bien cristalizado, lo que resaltó la importancia de la cristalinidad para determinar la
reactividad, siendo las fases poco cristalinas más reactivas que las bien cristalizadas de la misma composición.
Sin embargo, se debe considerar el tamaño del grano mineral con respecto a la accesibilidad de la superficie del
mineral al ácido; esto es una función de la exposición mineral en la partícula, que a su vez está relacionada con
la distribución del tamaño de partícula. Los estudios suelen considerar esto desde la perspectiva de los
minerales minerales que contienen cobre [26,27], pero no desde la perspectiva mineral de la ganga. Dada la
heterogeneidad que se puede esperar en la mineralogía de las partículas grandes (tamaño superior 25 mm),
junto con la distribución de la porosidad, un enfoque central cada vez menor de la reactividad de las partículas
[28] probablemente sea una simplificación excesiva. Por lo tanto, la distribución del tamaño del grano mineral,
en relación con la distribución del tamaño de partícula (PSD) del mineral, podría ser un parámetro probado
durante las pruebas en columna, así como también en el curado con ácido.
En la aglomeración, la distribución del tamaño de partícula del mineral triturado se estrecha para
proporcionar una porosidad óptima en el montón, junto con una introducción del ácido (ácido de curado)
antes de la lixiviación en pilas [1]. La aglomeración también reduce los efectos de segregación que podrían
afectar negativamente la permeabilidad en el montón. Bouffard [29] mostró que la adición de humedad, el PSD
del mineral y el tiempo de residencia afectan la calidad de los aglomerados producidos, todos los cuales tienen
implicaciones para la lixiviación en pilas. La humectación, el crecimiento, la consolidación y la rotura son los
pasos clave que deben darse para ayudar a
corregir las características de aglomeración. El consumo de ácido parece aumentar significativamente en
presencia de agua, por lo que el contenido de humedad es un factor que afecta la reactividad del ácido de
ganga durante el curado. La cantidad de finos y minerales arcillosos blandos requiere monitoreo, ya que estos
pueden afectar la humectabilidad y la retención de humedad en los aglomerados [29]. Por lo tanto, es
importante saber la mineralogía de ganga y las relaciones de textura en el mineral, ya que determinan la rotura
y, por lo tanto, el PSD alimentado al proceso de aglomeración.
Se ha demostrado que los tipos de trituración utilizados para la trituración tienen un efecto sobre la PSD
y la generación de grietas, la última de las cuales es deseable en partículas grandes para facilitar el acceso a los
minerales del mineral para la lixiviabilidad [26,27,30]. Por supuesto, la presencia de grietas y poros aumenta la
reacción por difusión, y también es aplicable en el contexto de la reactividad del ácido de ganga. Una
investigación exhaustiva de las relaciones texturales, incluida la distribución del tamaño del grano mineral en
el mineral, junto con las pruebas de conminución y consumo de ácido, debería ayudar a determinar las
relaciones entre la distribución del tamaño de grano, PSD y la reactividad, teniendo en cuenta los tipos de
minerales presentes.
El curado ácido se realiza antes de la lixiviación en pilas para acelerar la reacción con minerales de cobre,
a través de la reacción con H2SO4 concentrado. Sin embargo, una consecuencia clave para los minerales de
ganga es que el ácido fuerte deshidrata sílice polimérica en la ganga, inhibiendo así la reacción posterior de la
ganga y, por lo tanto, la disolución de sílice [19,31] bajo condiciones de lixiviación en pilas. Dado que la
aglomeración / curado tiene lugar con un contacto mucho mayor de mineral y ácido, y con una mayor fuerza
de ácido que los fluidos de percolación en el montón, se espera que las tasas de reactividad sean más altas
durante el curado. El estudio de Free [20] mostró diferentes tasas de reactividad para el clorito, la biotita y la
montmorillonita pulverizados puros en función de la fuerza del ácido, en presencia de altos índices de agua:
sólidos; la relación entre la fuerza del ácido y la reactividad de ganga no está bien documentada en la
literatura, sin embargo, particularmente considerando diferentes tamaños de partícula de partida en la
aglomeración (determinada por la capacidad de trituración del mineral), en comparación con los montones
(determinados por el tamaño del aglomerado). La cantidad de ácido que asegura un curado efectivo tampoco
se entiende bien desde una perspectiva fundamental; las dosificaciones de ácido óptimas para el curado se
determinan más empíricamente para una mineralogía de ganga determinada, el pH de la solución de
lixiviación y la PSD [9,32].
Trabajo reciente de Lu et al. [33] señaló la escasez en investigación fundamental para informar la práctica
de curado y aglomeración con respecto a las reacciones ácidas, y mostró que el curado con ácido debe tener
lugar con la adición de humedad para obtener una buena aglomeración y contacto superficial efectivo con
H2SO4 para reaccionar con partículas. La adición de humedad para una reactividad efectiva parece encontrar
apoyo en el trabajo de Free [20], para el cual se usaron proporciones de agua alta: sólidos. El trabajo de Lu et al.
[33] confirmaron además las observaciones de Snäll y Liljefors [18] de que el clorito, particularmente cuando el
Fe sustituye al Mg, es un filosilicato relativamente reactivo, que se degrada durante el curado con ácido. La
lixiviación de columna posterior eliminó de manera efectiva los productos formados a través de la interacción
con el ácido. Esta información sugiere que la reactividad ácida a largo plazo de clorita descrita por Jansen y
Taylor [21] puede acortarse considerablemente en condiciones de curado, de modo que el clorito, aunque
moderadamente reactivo en condiciones de lixiviación en pilas (Tabla 4), puede considerarse altamente
reactivo bajo ácido regímenes de curación. Por lo tanto, la caracterización cuidadosa del mineral antes y
después del curado debe establecer qué minerales reaccionan durante la duración del curado y en qué medida.
En particular, la difracción de rayos X puede verificar la descomposición estructural de los minerales,
particularmente en las fracciones más finas del tamaño de partícula [34]. Además, se necesitaría llevar a cabo
tal trabajo en los minerales de mineralogía de ganga diferente para establecer diferencias en la reactividad de
la ganga en condiciones de curado y acumulación. Esto también abordaría la cuestión de la cantidad de ácido
que se debe agregar para un curado exitoso.
Dentro del montón, la dinámica estructural y el flujo de fluidos deberían tenerse en cuenta, dada la
reactividad por flujo de ácido canalizado frente a difusivo en el montón [1]. Tales flujos se han modelado
utilizando la tomografía de emisión de positrones y la tomografía computarizada de rayos X [1,35]. La
reactividad con la ganga conllevaría una combinación de reactividad bajo estos regímenes de flujo de fluidos
que no se tiene en cuenta adecuadamente en las pruebas a escala de laboratorio. Además, la reactividad del
ácido de ganga da como resultado la lixiviación para producir la precipitación de fases tales como sulfatos de
Al y jarositas, que pueden conducir a obstrucción del espacio poroso y, por lo tanto, a una menor
permeabilidad del montón [36]. También se deben considerar los efectos de la temperatura cuando se emplea
biolixiviación en el montón, dado el clima más frío bajo el cual las operaciones
puede necesitar ser llevado a cabo. Los efectos de temperatura no parecen estar bien establecidos para las
interacciones entre ganga y ácido en el nivel fundamental. De nuevo, se realizan pruebas empíricas para
evaluar las diferencias cinéticas en la lixiviabilidad del mineral en condiciones operativas de verano e invierno
[9].
Un estudio más reciente [38] relacionó los tipos de roca y alteración en un depósito de pórfido de cobre
con las recuperaciones de cobre y el consumo de ácido, a través de una caracterización mineralógica detallada,
pruebas cinéticas de minicolumnas y pruebas de neutralización automatizadas. Los aspectos operativos se
superpusieron a las características fundamentales de las rocas, considerando la estacionalidad (precipitación y
temperatura sobre el período operativo) y la duración y frecuencia del riego. Los hallazgos clave del estudio
fueron:
- La diorita alterada potásica (roca ígnea básica rica en biotita) mostró el mayor consumo de ácido y la
menor recuperación de cobre de ~ 60%
- Las mejores recuperaciones de cobre (~ 80%) se obtuvieron de monzonita de cuarzo alterada potásica
(rocas ácidas ricas en cuarzo con alteración potásica del feldespato K, en lugar de biotita )
- Teniendo en cuenta los estilos de alteración del pórfido de cobre, la alteración fílica de la diorita mostró
un bajo consumo de ácido (presencia de cuarzo / moscovita / sericita)
- La alteración argílica no mostró un consumo significativo de ácido, debido a la presencia de caolinita, que
está en equilibrio con las soluciones de lixiviación. La recuperación de cobre fue eficiente desde tales
rocas.
- La alteración argílica avanzada, por otro lado, generó problemas de percolación y devolvió una pobre
recuperación de cobre, que los autores atribuyeron a la presencia de esmectita, la retención de ácido en
las estructuras de lámina internas y la captura de cobre a través del intercambio iónico [22]. El ácido
también puede mantenerse en las estructuras laminares externas de las arcillas [23].
- Las rocas brechadas en el depósito mostraron una mejor porosidad que las rocas plutónicas, lo que
conduce a un ataque más penetrante de fluidos, más formación de minerales arcillosos y posteriores
problemas de percolación. Además, se descubrió que el cuarzo secundario encapsula los minerales de
cobre, lo que provoca una recuperación pobre del cobre.
- Se encontró que la reactividad mineral de Gangue era inicialmente rápida, seguida de una disminución
en la tasa de reactividad. La diferencia entre las recuperaciones esperadas y reales se debió a factores
operacionales. Sin embargo, las pruebas de minicolumna permitieron una buena previsibilidad para el
comportamiento del mineral y permitieron el pronóstico de producción de las operaciones mineras .
El estudio [38] relaciona muy bien los efectos mineralógicos fundamentales discutidos en la sección 2, y
muestra que el tipo de roca y el tipo de alteración en un depósito de pórfido de cobre son factores valiosos para
predecir el comportamiento del mineral en la plataforma de lixiviación. En este caso, un enfoque
geometalúrgico exige un modelo geológico sólido que sustenta las características fundamentales del
comportamiento del mineral, a partir de las cuales se puede predecir mejor el comportamiento en la lixiviación
en pilas, lo que permite tomar decisiones informadas sobre planificación de la producción, almacenamiento y
manejo durante las operaciones. Como la caracterización de los minerales de pórfidos de cobre en términos de
roca y tipos de alteración es muy conocida en la geología económica [9,10,12,13], la aplicación de este
conocimiento para la lixiviación en pilas solo puede ser beneficiosa.
El enfoque geometalúrgico puede comprender los niveles estratégicos y tácticos, según lo definido por
MacKay et al. [39]. El enfoque estratégico comprende la capacidad de predicción a partir de datos en
profundidad derivados de pruebas de laboratorio (por ejemplo, pruebas de columna), mineralogía,
geoquímica y otras pruebas (la
base de datos geometalúrgica) para modelar la variabilidad espacial del mineral y así identificar el valor sobre
la vida de la mina. El enfoque táctico, por otro lado, tiene un enfoque de corto a mediano plazo que también
hace uso de la base de datos geometalúrgica, pero incorpora nuevos datos derivados de las operaciones
cotidianas, muy similar a lo previsto por Baum et al. [4], para mejorar la eficiencia del proceso (que abarca
todas las operaciones de la unidad) en términos de valor realizado .
5. Conclusiones
El estudio de las interacciones entre ganga y ácido es intrínseco a las operaciones de lixiviación en pilas,
dada la composición de ganga> 95% de los minerales tratados para la recuperación de cobre, particularmente
en minerales de pórfido. La reactividad relativa de los silicatos comunes proporciona una plataforma desde la
cual predecir el comportamiento del mineral al nivel mineral-químico fundamental. Este conocimiento debe
ampliarse con la caracterización sistemática completa de los minerales antes y después de la reacción durante
las diferentes operaciones de la unidad, y relacionándose con los tipos de minerales presentes en el yacimiento.
Para fines de curado, el conocimiento de la composición mineral inicial del mineral triturado y la relación
de la distribución del tamaño del grano mineral con la distribución del tamaño de partícula se debe probar
mediante una minuciosa investigación mineralógica. Los tipos de minerales presentes y sus relaciones de
textura, deben estar más relacionados con el método de trituración usado, la fuerza del ácido aplicada y el
ácido consumido durante la duración del curado. El consumo de ácido parece aumentar con la adición de
agua, por lo que el contenido de humedad es otro factor a considerar para la reactividad del ácido de ganga en
estas operaciones. Por lo tanto, la caracterización previa y posterior al curado mediante técnicas mineralógicas
es necesaria para explicar las reacciones y la posible degradación estructural de los minerales bajo diferentes
condiciones de curado. Los resultados de estas investigaciones determinarán la adición óptima de ácido y agua
para el curado efectivo de diferentes tipos de mineral.
Una comparación de la reactividad del ácido de ganga en condiciones de curado frente a lixiviación en
pilas sería importante para refinar la jerarquía de la reactividad de la ganga bajo estos diferentes regímenes.
Las condiciones del montón para la evaluación incluyen parámetros tales como la adición de ácidos de menor
fuerza (en comparación con las condiciones de curado), la temperatura y la dinámica del flujo de fluidos que
deberían estar relacionados con el consumo de ácido de ganga y la textura del mineral.
Un componente vital de las investigaciones comparativas anteriores debe ser el uso integrado de la
mineralogía para caracterizar minerales y productos en diversas operaciones unitarias. La caracterización
detallada del mineral inicial constituye la base para identificar problemas clave que pueden seguirse con
técnicas específicas. En particular, las pruebas de diagnóstico de lixiviación deben calibrarse contra la
información mineralógica para proporcionar resultados realistas. Se prevén herramientas de rutina como XRD,
microscopía electrónica de barrido óptica y (automatizada) y microanálisis con sonda de electrones, así como
capacidades más avanzadas de generación de imágenes por infrarrojos y XCT para lograr estos objetivos.
La recuperabilidad del cobre por lixiviación en pilas se ha entendido tradicionalmente a través de la
experiencia en lugar de la aplicación científica fundamental. Sin embargo, el potencial emocionante ofrecido
por un enfoque geometalúrgico debería permitir una caracterización más robusta en la comprensión de las
relaciones fundamentales entre minerales, partículas, minerales, aglomerados, fuerza y cantidad de ácido,
dinámica de montón espacial y temporal para el proceso de lixiviación en pilas, y la relación final de estos
parámetros con la distribución de litología y alteración en el cuerpo mineral. El modelado y la simulación de
estos aspectos, junto con otros criterios de diseño, deberían conducir a un mejor control operativo y una
recuperación de metal optimizada de las operaciones diarias.
Agradecimientos: se deben divulgar todas las fuentes de financiación del estudio. Indique claramente las
subvenciones que ha recibido en apoyo de su trabajo de investigación. Indique claramente si recibió fondos para
cubrir los costos de publicación en acceso abierto.
© 2018 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC BY) license (http://creativecommons.org/licenses/by/4.0/).