Litio Meteorización de Silicatos en Marte
Litio Meteorización de Silicatos en Marte
Litio Meteorización de Silicatos en Marte
2018
sus isótopos como trazadores geoquímicos
TESIS DOCTORAL
2018
A mis padres,
hermanas
y a la pequeñita de la casa
Índice
Índice
Introducción………………………………………………………………………………. 1
Objetivos de investigación…………………………………………………………………... 6
Organización de la tesis…………………………………………………………………........ 7
Capítulo 1
Planteamiento del problema…………………………………………………… 11
1.1. La estabilidad del agua en Marte y en la Tierra: el punto de
divergencia…………………………………………………………………………………………… 13
1.1.1. La firma isotópica del litio como marcador de meteorización………. 15
1.1.2. El ciclo del litio en condiciones marcianas………………………….............. 16
1.1.3. La creación de sobresaturación en minerales primarios de Marte... 17
1.2. Los diversos mecanismos de formación de los minerales de arcilla en
Marte……………………………………………………………………………………………………. 20
1.2.1. La llegada del Curiosity a Marte: nuevas técnicas de análisis in
situ…………………………………………………………………………………………………….. 21
1.2.2. Problemas de caracterización de los minerales de arcilla…………….. 22
1.2.3. Minerales de arcilla formados a baja temperatura: estudio del
desorden estructural………………………………………………………………….............. 23
1.2.4. Estudio del estado de oxidación de los minerales de arcilla.………… 24
1.2.5. Relación entre la temperatura y el tiempo en el proceso de
cristalización…………………………………………………………………………………........ 25
1.3. Estabilidad temporal de las masas de agua en Marte………………............... 26
1.3.1. Modelos geoquímicos: Li y sus isótopos…...…………………………………. 26
Capítulo 2
Instrumentación y tratamiento de datos…………………………………. 27
2.1. Análogos de laboratorio del proceso de disolución/cristalización de
basalto a baja temperatura……………………………………………………………………. 31
2.1.1. Técnicas electroquímicas para el seguimiento in situ de
parámetros físico-químicos y adquisición de datos …………………………….. 32
2.2. Técnicas de difracción para el estudio de estructuras cristalinas:
muestras terrestres…..…………………………………………………………………………… 33
2.2.1. Difracción de Rayos X...………………………………………………………………. 35
2.2.2. Difracción de Neutrones de Alta Resolución………….…………………….. 35
2.2.3. Tratamiento de los datos difractométricos…………………………………. 36
2.3. Técnicas del estudio de muestras de Marte………………………………………. 38
2.3.1. Instrumento Chemistry & Mineralogy, CheMin…………………………….. 38
2.3.1.1. Toma de muestras………………………………………………………………... 39
i
Índice
ii
Índice
iii
Índice
iv
Índice
v
Capítulo 8
Introducción
Introducción
Por consiguiente, el viento y la niebla ácida son los principales agentes que
causan meteorización, y las reacciones entre los gases atmosféricos y los
minerales pueden llegar a estar controlando la tasa y los productos de la
meteorización química [John Carter et al., 2015]. Sin embargo, la superficie
de Marte presenta evidencias de alteración acuosa, como morfologías
fluviales y lacustres [Ansan et al., 2011; Caleb I. Fassett and Head, 2008;
Hynek et al., 2010] o presencia de minerales de arcillas [J Carter et al., 2010;
Ehlmann et al., 2011; Caleb I Fassett and Head, 2011; Morris et al., 2000;
Murchie et al., 2009; Poulet et al., 2005] cuyos estudios apoyan la hipótesis
3
Introducción
Uno de los ciclos geoquímicos más relevantes para este trabajo es el del
carbonato-silicato [Urey, 1952], de retroalimentación negativa entre los
factores climáticos y la velocidad de meteorización química de las rocas. Al
igual que la Tierra, Marte está compuesto de rocas ígneas, cuyos minerales
principales son los olivinos, piroxenos y plagioclasas [Mustard et al., 2008;
Rogers et al., 2007]. La meteorización de las rocas basálticas está controlada
por la mineralogía primaria, las propiedades microestructurales que
condicionará el movimiento del agua a través de la roca [Meunier et al., 2007]
y las características físico-químicas del medio (Tª, pH, composición
atmosférica, presencia de H2O, etc.).
4
Introducción
El registro isotópico del litio tiene una respuesta específica a los cambios en
la meteorización de las rocas silicatadas, otorgándose por ello el carácter de
un buen proxy [Fairén et al., 2015]. El Li es retenido por los minerales
silicatados pudiéndose usar como trazador para buscar una conexión entre
la ruptura de los minerales silicatados y la reducción del CO2. Investigaciones
geoquímicas de isótopos de litio proporcionan información relevante sobre
el equilibrio entre los minerales silicatados en las rocas primarias y los
procesos geoquímicos asociados con su meteorización en ambientes
acuosos [Chan et al., 1992; Decarreau et al., 2012; Hathorne and James,
2006].
5
Introducción
Objetivos de investigación
6
Introducción
Organización de la tesis
7
Introducción
Análisis, mediante modelos geoquímicos, de las posibles vías por las que los
procesos acuosos podrían haber dado lugar al fraccionamiento isotópico del
Li en Marte. Aportaciones al conocimiento sobre el grado y la extensión de
la meteorización en Marte, siendo útil para determinar las condiciones
ambientales sobre el planeta en el pasado.
Conclusiones
9
Capítulo 1
Planteamiento del problema
Planteamiento del problema
13
Capítulo 1
[Cabrol and Grin, 2010; J Carter et al., 2010], que son depresiones de gran
envergadura que pudieron haber albergado cuerpos de agua durante este
periodo. Esto implicaría la existencia en el pasado de un clima más húmedo
que el existente hoy en día, y una atmósfera de mayor espesor que impida la
evaporación inmediata del agua [John Carter et al., 2015; Ramirez et al.,
2014]. Si esto hubiese sido así, aumentarían las probabilidades de que Marte
fuese en su pasado un planeta habitable.
14
Planteamiento del problema
16
Planteamiento del problema
17
Capítulo 1
18
Planteamiento del problema
19
Capítulo 1
20
Planteamiento del problema
Las principales ventajas que aporta el rover a este trabajo son: (1) la
capacidad de detectar directamente el litio, por medio de su instrumento
Chemistry and Camera, ChemCam [Wiens et al., 2012]. (2) La disponibilidad
de un equipo de difracción de rayos-X, que permite obtener datos de
minerales de arcilla dentro de un contexto identificable y analizado
químicamente, gracias al conjunto de equipos que contiene a bordo.
21
Capítulo 1
22
Planteamiento del problema
Por otro lado, la baja resolución de los datos obtenidos por el CheMin, así
como la falta de toma de datos de los minerales de arcilla como fase única,
hacen realmente complicado el análisis cuantitativo de las mismas. En este
trabajo se van a presentar métodos que permiten una aproximación a la
información clave sobre la estructura y química del cristal, como por ejemplo
el Simulated Annealing o el FAULTS [Casas-Cabanas et al., 2006; Rodríguez-
Carvajal, 1993].
23
Capítulo 1
24
Planteamiento del problema
𝑑𝐺 𝜕𝑇 𝜕𝐺
= ( ) 𝑑𝑡 + ( ) 𝑑𝑇 [1.1]
𝑑𝑡 𝜕𝑡 𝜕𝑇
25
Capítulo 1
26
Capítulo 2
Instrumentación y tratamientos de datos
Instrumentación y tratamiento de datos
29
Instrumentación y tratamiento de datos
31
Capítulo 2
32
Instrumentación y tratamiento de datos
33
Capítulo 2
Es una técnica no destructiva que fue creada para obtener información del
estado físico-químico de las muestras, se fundamenta en la interferencia
producida entre el material sólido cristalino y la fuente de radiación –en este
caso, rayos-X o neutrones–. Cuando existe orden y periodicidad en la
muestra, se obtiene un dominio coherente de difracción. Los átomos
distribuidos en los planos de reflexión (hkl), separados entre sí por espacios
interplanares constantes, dhkl, actúan como centros dispersores. La longitud
de onda del haz monocromático ha de ser próxima al espacio interlaminar,
para que al someter la muestra a dicho haz se dé el patrón de interferencia
constructiva, en función de la dirección de reflexión. Como resultado, se
obtiene el difractograma de intensidad frente al ángulo de reflexión. Esto se
explica mediante la Ley de Bragg:
34
Instrumentación y tratamiento de datos
35
Capítulo 2
36
Instrumentación y tratamiento de datos
37
Capítulo 2
38
Instrumentación y tratamiento de datos
Las muestras son tamizadas con una malla de 150 μm. Cuando abre la
cubierta, se somete el embudo a vibración mientras va introduciendo 75
mm2 de la muestra, de los cuales a la celda se destinan tan sólo 15 mm2 para
realizar el análisis y el resto se emplea para homogenizar tanto el embudo
como la celda a usar. Dispone de dos tipos de celdas en la rueda
portamuestras: de mylar y de kapton. Las de kapton son más duraderas pero
interfieren con determinados tipos de minerales (p. ej. filosilicatos) por dar
una señal de difracción a ~6.7° 2θ. Placas piezoeléctricas agitan la muestra
para conseguir enfrentar caras aleatorias de los cristales al haz de rayos-X.
Tras la medida la muestra es desechada para reutilizar la celda.
Dispone de un ánodo de cobalto, Co, que produce rayos-X a 6.9 keV (Kα) y 7.6
keV (Kβ). El haz colimado atraviesa la muestra y los rayos-X difractados son
recogidos por el detector (CCD). Un criostato se encarga de mantener la
temperatura del detector entre -60°C y -50°C durante las medidas. Las
imágenes recogidas durante el análisis son de 600x582 píxeles. La
frecuencia normal de toma de datos es de una imagen cada 15 segundos. Los
39
Capítulo 2
2.3.1.3. Calibración
2.3.2.2. Calibración
41
Capítulo 2
2.3.3.3. Calibración
42
Instrumentación y tratamiento de datos
43
Capítulo 3
Estudio de los mecanismos de meteorización del
basalto mediante modelización geoquímica y
análogos de laboratorio
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
47
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
𝛺 = 𝐼𝐴𝑃/𝐾𝑆 [3.1]
50
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
Figura 3.1.
Variación de la
velocidad de
disolución de los
silicatos, en
función del pH y
de la temperatura
(K). Los silicatos
son más estables
en torno a un
valor neutro de
pH.
51
Capítulo 3
𝑑𝑚 𝐸𝑎 1 1
𝑛
𝐸𝑛 1 1
= −𝑆𝐴 (𝐾𝑎 𝑒 − 𝑅 (𝑇−298.15𝐾) 𝑎𝐻1+ + 𝐾𝑛 𝑒 − 𝑅 (𝑇−298.15𝐾)
𝑑𝑡
[3.3]
𝐸𝑏 1 1
𝑛
+ 𝐾𝑏 𝑒 − 𝑅 (𝑇−298.15𝐾) 𝑎𝐻2+ ) (1 − 𝑝 𝑞
𝛺 )
Si se supone que hay condiciones alejadas del equilibrio (es decir, cuando la
fase mineral está sobresaturada o subsaturada) el término de la
sobresaturación (1 − 𝛺), que contiene todas las dependencias
termodinámicas, pasa a ser cuantitativamente menos importante (en
comparación con la superficie reactiva). De modo que la cantidad de mineral
que precipita o se disuelve está principalmente controlada por las
constantes cinéticas y la superficie reactiva.
52
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
53
Capítulo 3
Figura 3.3. Evolución de una suspensión de basalto en agua hacia distintos estados
de pH en función de la superficie reactiva; (a) alcalinos con alta superficie reactiva
(b) ácidos en basalto sin fragmentar en presencia de volátiles (p. ej. CO 2, H2S, SO2).
55
Capítulo 3
56
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
Figura 3.4. Modelo geoquímico de los factores que inducen sobresaturación durante
al meteorización del basalto en Marte dando lugar a la precipitación de minerales
secundarios. Junto a las variables inductoras del proceso de sobresaturación:
evaporación y enfriamiento, actúan la disolución de volátiles en el agua
proporcionando los contra-iones necesarios. La extrusión de sales durante la
formación de hielo y la formación de gradientes de densidad van a dar al precipitación
en el fondo en función de las condiciones de reactividad del sustrato basáltico. Figura
modificada de [Fairén et al., 2017].
57
Capítulo 3
58
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
Tabla 3.2. Superficie reactiva de los minerales del basalto como función de la fracción
molar del mineral en la superficie del basalto, Xmin.
Fase SR (m 2/Kgagua) Comp. Modal (% vol.) X min mol0/kgagua
Anortita 103.40 50.00 0.50 0.10
Forsterita 14.75 15.00 0.20 0.07
Fayalita 9.90 15.00 0.10 0.06
Diopside 7.67 15.00 0.10 0.05
Enstantita 2.40 5.00 0.10 0.03
59
Capítulo 3
Uno de los objetivos de este estudio ha sido mostrar cómo las diferentes
velocidades de la generación de sobresaturación (y no solo un valor concreto
de este parámetro) pueden influir considerablemente en la secuencia y el
tipo de fases en el sistema, cuando se considera un enfoque cinético del
problema de la precipitación de fases secundarias.
Cabe destacar que uno de los principales problemas técnicos para analizar
la formación de minerales secundarios en el laboratorio, a partir de la
disolución del basalto, es cómo obtener suficientes cationes en solución para
alcanzar el estado de saturación para sulfatos o cloruros. Esto se debe a la
baja solubilidad de los minerales que contiene el basalto. Incluso, cuando los
volátiles atmosféricos llegan a acidificar suficientemente la solución –
promoviendo un aumento de la disolución– se alcanzan condiciones
próximas al equilibrio cuando se detiene la disolución del basalto. Esta baja
solubilidad conduce a valores bajos de fuerza iónica (es decir, la
conductividad medida se encuentra en el rango de 10–20 μS/cm).
60
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
𝑆
𝑄𝑊 = (∆𝑃 − ∆𝜋𝑜𝑠𝑚 ) ∙ 𝐾𝑊 ∙ [3.4]
𝑑
El reactor diseñado combina tres reservorios, de los cuales dos de ellos (T1
y T2) son utilizados para inducir la disolución del basalto, operando o bien
en modo “batch” o como reactor de flujo continuo. El tercer reservorio (T3)
es el lugar donde se produce la mezcla de la solución resultante de disolver
el basalto en cualquiera de las dos condiciones (próxima o alejada del
61
Capítulo 3
Figura 3.6. (a) Esquema del dispositivo diseñado para simular procesos de
evaporación a baja temperatura, combinando disolución de basalto en suspensión
dada en T1 y T2 y la solución salina resultante de la osmosis inversa, T3. (b) a)
Imagen del reactor y membranas de ósmosis inversa. b) Imagen del reservorio de
mezcla, T3 y las sondas de medida en continuo de conductividad y pH. c) Informática
para el registro continuo de datos.
62
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
63
Capítulo 3
𝑞
donde 𝑅 es la velocidad de disolución; υ𝑖 el coeficiente estequiométrico; ⁄𝐴
el volumen de solución a través del sistema por unidad de área; y el término
ΔC que viene dado por (𝐶𝑠𝑎𝑙𝑖𝑑𝑎 − 𝐶𝑒𝑛𝑡𝑟𝑎𝑑𝑎 ) representa la diferencia de
concentraciones entre la solución de salida del reactor y la solución de
entrada que, en estado estacionario, tiende a ser constante. En condiciones
alejadas del equilibrio, el valor determinado experimentalmente representa
la constante específica de disolución 𝐾𝑖 –para unas condiciones
determinadas de pH– tal y como está definida en la expresión [3.1], ya que el
término de sobresaturación (1 − 𝛺) al aproximarse a 1 no influye en el valor
medido.
64
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
Figura 3.7. Esquema del set-up para el funcionamiento del equipo como baño de flujo
continuo.
65
Capítulo 3
𝑑𝐶 𝑉
𝑅υ𝑖 = − ∙ [3.6]
𝑑𝑡 𝐴
66
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
67
Capítulo 3
Figura 3.8. Sistema de reflujo empleado para generar soluciones con concentraciones
cercanas a la concentración de saturación con el basalto.
68
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
punto. Los valores obtenidos fueron empleados para fijar unas condiciones
óptimas de la relación agua/roca en el reservorio T3 del reactor, donde se
produce el intercambio entre la solución procedente del equilibrio con el
basalto y la solución salina, procedente de las membranas de ósmosis.
69
Capítulo 3
70
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
71
Capítulo 3
Figura 3.10. Reactor de mezcla (T3) entre la solución procedente del basalto (T2) y el
agua salina resultante de la ósmosis inversa. El sistema permite modelizar diferentes
procesos de creación de sobresaturación, modelizando la evaporación a baja
temperatura, en condiciones similares a las que han podido tener lugar en Marte. El
sistema permite el seguimiento de la evolución de la conductividad, pH, ORP y oxígeno
disuelto en tiempo real.
72
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
también va a ser baja, con valores del orden de 20-40 μS/cm (para la masa
de basalto y el tamaño de grano empleados) por lo que el proceso de
generación de sobresaturación progresa muy lentamente. La ósmosis solo es
eficaz para aumentar la concentración cuando partimos de una
concentración inicial elevada. Sin embargo esta configuración es muy útil
para analizar las etapas iniciales del proceso de sobresaturación del basalto
así como los parámetros físico-químicos asociados. Un experimento de este
tipo, que combina la disolución en condiciones alejadas del equilibrio con la
creación de sobresaturación por ósmosis, se muestra en la figura 3.11. Las
gráficas de la figura 3.11a muestran el aumento de concentración en
ausencia de CO2 mientras que en el caso 3.11b el experimento de disolución
se realizó con una atmósfera de CO2, cuya PCO2 = 1 bar.
Figura 3.11. Etapas iniciales del proceso de creación de sobresaturación a partir del
basalto. Las dos figuras representan curvas de evolución de la conductividad, pH y ORP
en experimentos que combinan ciclos de ósmosis inversa y la concentración resultante
de la disolución del basalto en condiciones alejadas del equilibrio, en un reactor de
flujo continuo. (a) Sin CO2 (ac.), (b) con CO2 disuelto.
74
Estudio de los mecanismos de meteorización del basalto mediante modelización geoquímica y
análogos de laboratorio
75
Capítulo 3
este tipo de proceso va a corresponder a las sales más solubles, que van a
depender de la reacción con los volátiles atmosféricos que actúen como
contra-iones en el sistema.
3.5. Conclusiones
77
Capítulo 4
Minerales de arcilla detectados por Difracción de
Rayos-X sobre el cráter Gale
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
81
Capítulo 4
82
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
83
Capítulo 4
84
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
1) Grupo Bradbury
85
Capítulo 4
Las tres muestras (Tabla 4.1) poseen una consistencia similar: rocas
sedimentarias compactas de mudstones, que han sido muestreadas con la
rover’s drill. Las alícuotas fueron previamente tamizadas (<150 μm) y los
86
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
Tabla 4.1. Datos técnicos de las muestras analizadas en este estudio. El Sol de
descarga de datos, imágenes menores tomadas, tiempo de integración acumulado
durante varias noches de medida, tipo de ventana del portamuestra y rango de medida
del equipo.
Tiempo de
Muestra Sol Imágenes Menores Tipo de ventana Rango 2θ (paso)
integración
Marimba MB 1426 60 30.000 (8b) Mylar 2 – 52° (0.05)
Quela QL 1472 60 30.000 (5a) Kapton 2 – 52° (0.05)
Sebina SB 1499 51 25.500 (4b) Kapton 2 – 52° (0.05)
87
Capítulo 4
88
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
89
Capítulo 4
90
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
Las diferencias entre los patrones extraídos por el LBF, la ausencia de picos
abruptos –característicos de la estructura ideal de los minerales de arcilla–y
la ausencia de convergencia con el Simulated Annealing han sido
interpretadas como posibles desórdenes estructurales. La presencia de
bandas anchas, donde debería haber un pico bien definido, son indicativos
de cierto grado de desorden planar o de apilamiento.
91
Capítulo 4
92
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
93
Capítulo 4
94
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
95
Capítulo 4
97
Capítulo 4
98
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
Figura 4.8.
Intensidades
integradas del
mineral de arcilla,
extraídas de sus
correspondientes
patrones
generales,
empleando el
método de ajuste
Le Bail.
99
Capítulo 4
100
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
Tabla 4.4. Grados 2θ y dhkl de las principales reflexiones hkl de los patrones de
difracción ajustados por LBF en una primera aproximación a la información contenida
en el difractograma.
Sebina Quela Marimba
(hkl)
2 Theta dhkl (Å) I (cal.) 2 Theta dhkl (Å) I (cal.) 2 Theta dhkl (Å) I (cal.)
(0 0 1) 10.181 10.082 322.9 10.478 9.798 281.4 10.455 9.819 714.5
(0 0 2) 20.444 5.041 0 21.044 4.899 8.9 20.997 4.910 0.300
(0 2 0) 23.314 4.428 21.6 23.719 4.353 94.7 24.771 4.171 13.5
(-1 1 2) 28.164 3.677 87.6 28.046 3.692 82.3 28.581 3.624 141.2
(1 1 2) 33.598 3.095 4.5 34.636 3.005 0 33.939 3.065 4.2
(0 0 3) 30.876 3.361 165.5 31.795 3.266 49 31.724 3.273 114.9
101
Capítulo 4
Por otro lado, los picos (-112) y (112) pueden aumentar o disminuir en
función del orden. Si tenemos en cuenta que el refinamiento determinó
porcentajes en peso similares para la fase arcillosa en las tres muestras
(Tabla 4.3), las intensidades de la tabla 4.4 indican que la muestra más
ordenada podría ser la de Marimba, seguida por Sebina y por último la de
mayor desorden podría ser la de Quela .
102
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
Entre ellos los minerales de capa mezclada se pueden llegar a formar cuando
dos fases con estructura similar se ven favorecidas para su precipitación
bajo las mismas condiciones o cuando el proceso de difusión se ve
interrumpido por sedimentación, compactación o cementación.
103
Capítulo 4
104
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
Tabla 4.5. Átomos con sus coordenadas fraccionarias y ocupaciones del mineral
glauconítico en P1, correspondiendo a la capa sin rotar.
Átomo nº x y z Occ. Átomo nº x y z Occ.
K+ 1 0.000 0.000 0.000 0.80 Si4+ 38 0.008 0.670 0.228 0.90
K+ 2 0.500 0.500 0.000 0.80 Si4+ 39 0.508 0.830 0.228 0.90
Al3+ 3 0.332 0.333 0.500 0.26 Si4+ 40 0.492 0.830 0.772 0.90
3+ 4+
Al 4 0.668 0.333 0.500 0.26 Si 41 0.492 0.170 0.772 0.90
3+ 4+
Al 5 0.668 0.667 0.500 0.26 Si 42 0.508 0.170 0.228 0.90
Al3+ 6 0.332 0.667 0.500 0.26 Al3+ 43 0.008 0.330 0.228 0.10
Al3+ 7 0.832 0.833 0.500 0.26 Al3+ 44 0.992 0.330 0.772 0.10
Al3+ 8 0.168 0.833 0.500 0.26 Al3+ 45 0.992 0.670 0.772 0.10
3+ 3+
Al 9 0.168 0.167 0.500 0.26 Al 46 0.008 0.670 0.228 0.10
Al3+ 10 0.832 0.167 0.500 0.26 Al3+ 47 0.508 0.830 0.228 0.10
Fe3+ 11 0.332 0.333 0.500 0.44 Al3+ 48 0.492 0.830 0.772 0.10
Fe3+ 12 0.668 0.333 0.500 0.44 Al3+ 49 0.492 0.170 0.772 0.10
Fe3+ 13 0.668 0.667 0.500 0.44 Al3+ 50 0.508 0.170 0.228 0.10
3+ 2-
Fe 14 0.332 0.667 0.500 0.44 O 51 0.488 0.000 0.176 1.00
Fe3+ 15 0.832 0.833 0.500 0.44 O2- 52 0.512 0.000 0.824 1.00
Fe3+ 16 0.168 0.833 0.500 0.44 O2- 53 0.988 0.500 0.176 1.00
Fe3+ 17 0.168 0.167 0.500 0.44 O2- 54 0.012 0.500 0.824 1.00
3+ 2-
Fe 18 0.832 0.167 0.500 0.44 O 55 0.510 0.181 0.390 1.00
Fe2+ 19 0.332 0.333 0.500 0.20 O2- 56 0.490 0.181 0.610 1.00
Fe2+ 20 0.668 0.333 0.500 0.20 O2- 57 0.490 0.819 0.610 1.00
Fe2+ 21 0.668 0.667 0.500 0.20 O2- 58 0.510 0.819 0.390 1.00
Fe2+ 22 0.332 0.667 0.500 0.20 O2- 59 0.010 0.681 0.390 1.00
2+ 2-
Fe 23 0.832 0.833 0.500 0.20 O 60 0.990 0.681 0.610 1.00
Fe2+ 24 0.168 0.833 0.500 0.20 O2- 61 0.990 0.319 0.610 1.00
Fe2+ 25 0.168 0.167 0.500 0.20 O2- 62 0.010 0.319 0.390 1.00
Fe2+ 26 0.832 0.167 0.500 0.20 O2- 63 0.768 0.240 0.165 1.00
2+ 2-
Mg 27 0.332 0.333 0.500 0.10 O 64 0.232 0.240 0.835 1.00
2+ 2-
Mg 28 0.668 0.333 0.500 0.10 O 65 0.232 0.760 0.835 1.00
Mg2+ 29 0.668 0.667 0.500 0.10 O2- 66 0.768 0.760 0.165 1.00
Mg2+ 30 0.332 0.667 0.500 0.10 O2- 67 0.268 0.740 0.165 1.00
Mg2+ 31 0.832 0.833 0.500 0.10 O2- 68 0.732 0.740 0.835 1.00
2+ 2-
Mg 32 0.168 0.833 0.500 0.10 O 69 0.732 0.260 0.835 1.00
Mg2+ 33 0.168 0.167 0.500 0.10 O2- 70 0.268 0.260 0.165 1.00
Mg2+ 34 0.832 0.167 0.500 0.10 O2- 71 0.965 0.000 0.393 1.00
Si4+ 35 0.008 0.330 0.228 0.90 O2- 72 0.035 0.000 0.607 1.00
Si4+ 36 0.992 0.330 0.772 0.90 O2- 73 0.465 0.500 0.393 1.00
4+ 2-
Si 37 0.992 0.670 0.772 0.90 O 74 0.535 0.500 0.607 1.00
105
Capítulo 4
106
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
El modelo que mejor ajustó los datos observados (SB χ2 = 0.86, QL χ2 = 2.51
y MB χ2 = 5.04) incluye una secuencia de capas apiladas rotadas 120° entre
sí, con una probabilidad de existencia de 33%. La secuencia resultó ser la
misma para las tres muestras estudiadas, diferenciándose entre ellas en los
vectores de apilamiento de las capas.
Tabla 4.5. Datos de probabilidad de existencia de cada capa αij (siendo 1 la capa ideal,
2 la capa ideal rotada 120° y 3 la capa ideal rotada 240°), con sus vectores de
traslación asociados (Rxij, Ryij, Rzij) para cada una de las muestras.
Sebina Quela Marimba
i j αij Rxij Ryij Rzij αij Rxij Ryij Rzij αij Rxij Ryij Rzij
1 1 0 0.000 0.000 1.000 0 0.000 0.000 1.000 0 0.000 0.000 1.000
1 2 1 -0.711 0.000 1.000 1 -0.947 0.000 1.000 1 -1.079 0.000 1.000
1 3 0 0.000 0.000 1.000 0 0.000 0.000 1.000 0 0.000 0.000 1.000
107
Capítulo 4
Como se observa en la imagen del SAED (Fig. 4.12a) contiene una mezcla de
la fase en diferentes estados, cristalina representada por los spots brillantes
y amorfa simbolizada por la difuminación que presentan los anillos
concéntricos. Estos anillos también tienen que ver con los desórdenes de
apilamiento, debido a que la presencia de distintas orientaciones de las
108
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
capas dentro del mismo cristal sería semejante a una muestra policristalina
con orientación aleatoria.
109
Capítulo 4
Los minerales glauconíticos son del grupo de las micas alcalinas (esmectitas
de Fe) cuya fórmula química es la siguiente:
(K,Na)(Fe3+,Al,Mg,Fe2+)2(Si,Al)4O10(OH)2
111
Capítulo 4
y por otro, las celadonitas se puede dar por alteración hidrotermal del
basalto oceánico. La principal diferencia que genera la paragénesis entre
ambos minerales por ambientes de formación es la temperatura.
Entre estos minerales, puede que la opción que más se ajuste al ambiente
marciano presentado sea los minerales glauconíticos, su proceso de
cristalización se caracteriza por ser a bajas temperaturas y es facilitado por
la presencia de Mg2+ [Harder, 1980]. Por el contrario, la celadonita se
descartó porque es una fase de altas temperaturas, típica de alteración de
basalto y no de autigénesis en sedimentos. Siendo altamente viable la
presencia de interestratificados nontronita-glauconita.
112
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
113
Capítulo 4
114
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
El primer factor que ha de darse para que estos minerales cristalicen es que
se alcance la sobresaturación en el medio acuoso, en cuanto a los elementos
que lo forman (K, Fe, Mg, Al, Si y O). La principal fuente de los mismos va a
ser la meteorización de las rocas basálticas que condicionarán la
composición química del medio. El análisis del índice de alteración de las
rocas [Nesbitt and Young, 1989] de las localizaciones de las muestras indican
cierto grado de meteorización química, con valores superiores a los
reportados para la media de la corteza marciana, 37% [Taylor and
McLennan, 2009] (Fig. 4.14). Estos índices han sido calculados con los datos
obtenidos por el APXS a bordo del Curiosity (datos disponibles en la PDS). Un
estudio sobre este índice de alteración acuosa, junto a otros, se realizará en
el capítulo 7 por lo que aquí no se entrará en detalle.
Figura 4.14.
Representación del
Índice de Alteración
Química (IAQ) de la
zona de estudio.
115
Capítulo 4
116
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
117
Capítulo 4
Las condiciones óxicas reducen fuertemente la movilidad de los iones del Fe,
como indicativo estas localizaciones suelen presentar oxihidróxidos de Fe en
lugar de pirita. Al contrario que en la zona de estudio de este trabajo, en la
base del cráter -localizaciones John Klein y Cumberland- se han observado
bajas concentraciones de pirita y de akaganeita indicando posiblemente
condiciones anóxicas en el pasado. En John Klein, la presencia Magnetita
antigénica probablemente, es coherente con anoxia durante la deposición
[Schikorr, 1929; D T Vaniman et al., 2014], pudiendo haber sido la tasa de
deposición o sedimentación superior a la de oxidación por parte del O2 y
radiación UV de la superficie de la columna de agua.
118
Minerales de arcilla detectados por Difracción de Rayos-X sobre el cráter Gale
4.5. Conclusión
119
Capítulo 4
120
Capítulo 5
Estudio del fraccionamiento de isótopos de litio
con Difracción de Neutrones de Alta Resolución
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
123
Capítulo 5
124
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
125
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
127
Capítulo 5
Los minerales dioctaédricos, en los que dos de cada tres octaedros contienen
un catión en su interior (Al3+ en la estructura ideal), se caracterizan por
poseer una carga negativa significativa en sus capas silicatadas, debido
principalmente a sustituciones octaédricas. Dicha carga suele ser
contrarrestada por cationes de pequeño tamaño (p. ej. Li+), ya que son
susceptibles a cambios irreversibles tras sufrir un proceso de evaporación
[Greene-Kelly, 1952].
128
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
- Capa tetraédrica (T): formada por tetraedros de SiO4 unidos por los
bordes dando lugar a cavidades hexagonales.
- Capa octaédrica (O): formada por octaedros de AlO4(OH)2, que
compartiendo los bordes dan lugar, de igual modo que la anterior, a
formas hexagonales.
- Capa tetraédrica (T): dispuesta simétricamente a la primera con
respecto a la capa octaédrica.
- Espacio interlaminar entre cada conjunto T:O:T que puede albergar
diferentes cationes coordinados con moléculas de agua.
octaédricas [Decarreau et al., 2012; Jaynes and Bigham, 1987; Schultz, 1969;
Tettenhorst, 1962; Theng et al., 1997].
- Capa tetraédrica (T): formada por tetraedros de SiO4 unidos por los
bordes dando lugar a cavidades hexagonales.
- Capa octaédrica (O): formada por octaedros de AlO4(OH)2, que
compartiendo los bordes dan lugar, de igual modo que la anterior, a
formas hexagonales.
- Espacio interlaminar entre secuencias T:O donde se dan enlaces débiles
de tipo Van der Waals entre láminas T:O.
130
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
Figura 5.2.
Estructura de
la caolinita e
incorporación
de Li a sus
bordes
estructurales.
131
Capítulo 5
5.3. Metodología
La beidelita, (SBld-1, Idaho, EE.UU.) [Post et al., 1997], posee una capacidad
de intercambio catiónica de entre 70-100 meq/100g de arcilla. Su fórmula
estructural es la siguiente:
132
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
Es por ello que se empleó este procedimiento como una posible vía para
inducir la incorporación del litio a la estructura de las arcillas. Para ello, se
utilizó una disolución 3M de carbonato de litio (Li2CO3) en D2O/H2O. Se
mantuvo en agitación y se centrifugó repetidas veces, con el fin de
homogeneizar el catión de intercambio. Todas las disoluciones se realizaron
empleando una mezcla de agua deuterada y agua ligera des-ionizada
(D2O/H2O) con el fin de minimizar la dispersión incoherente asociado a los
hidrógenos de las moléculas de agua, y a los grupos hidroxilo de la
estructura. A continuación, se lavaron con D2O/H2O y se prepararon dos
133
Capítulo 5
En función del tamaño y carga del catión intercambiable, cada arcilla precisa
de una temperatura u otra para colapsar su estructura, siendo a partir de
200°C la adecuada para los isótopos de Li por su pequeño tamaño [Greene-
Kelly, 1952]. Por tanto, en las condiciones experimentales que se han
empleado es esperable que ocurra la migración total del Li al interior de las
estructuras cristalinas.
134
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
135
Capítulo 5
136
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
I(Q) = ∑ 𝑐𝛼 ̅̅̅
𝑏𝛼2 + F(Q) [5.1]
α
2
Δ𝐿𝑖 (𝑄) = 𝑐𝐿𝑖 (𝑏72𝐿𝑖 − 𝑏62𝐿𝑖 ) + 𝑐𝐿𝑖 (𝑏72𝐿𝑖 − 𝑏62𝐿𝑖 ) ∙ [𝑆𝐿𝑖𝐿𝑖 (𝑄) − 1]
[5.3]
+ Σ𝛼≠𝐿𝑖 2𝑐𝐿𝑖 𝑐𝛼 (𝑏7𝐿𝑖 − 𝑏6𝐿𝑖 ) ∙ [𝑆𝛼𝐿𝑖 (𝑄) − 1]
137
Capítulo 5
138
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
139
Capítulo 5
140
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
142
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
143
Capítulo 5
144
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
Como criterio de ajuste del modelo, se empleó el valor de χ2: una medida de
desviación entre el perfil calculado y el observado. Se considera aceptable si
es inferior a 5. En los ajustes realizados, se obtuvieron valores comprendidos
entre 2 y 4. Además, se comprobó la verosimilitud cristalográfica de la
estructura resultante. Para ello se representó la estructura resultante
empleando programas de representación estructural como el CrystalMaker
o el Vesta, y se analizó la geometría de los distintos poliedros que la
conforman.
145
Capítulo 5
146
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
Tabla 5.1. Fases presentes en cada muestra junto con sus porcentajes, volumen y
parámetros de celda unidad.
Vol. Celda
Muestra Fase a b c α β γ %
Unidad (Å3)
Caolinita 5.160 8.941 7.394 91.728 104.65 89.856 329.909 67.14
Caolinita_25
Cuarzo 5.526 5.526 5.401 90.000 90.000 120.00 142.808 32.86
Caolinita 5.162 8.942 7.396 91.659 104.66 89.868 330.121 54.43
Caolinita_300
Cuarzo 5.520 5.520 5.401 90.000 90.000 120.00 142.531 45.57
Beidelita_H 5.169 8.936 14.14 90.000 100.13 90.000 642.833 58.43
Beidelita_25 Beidelita_I 5.207 8.885 10.12 90.000 97.458 90.000 464.150 0.780
Cuarzo 4.898 4.898 5.393 90.000 90.000 120.00 112.048 40.79
Beidelita_H 5.218 8.812 11.93 90.000 100.29 90.000 546.789 11.76
Beidelita_300 Beidelita_I 5.140 8.968 10.05 90.000 96.121 90.000 460.445 4.710
Cuarzo 4.906 4.906 5.404 90.000 90.000 120.00 112.660 83.53
5.5.2.1. Beidelita
148
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
Figura 5.8. Refinamiento Rietveld. (a) Beidelita_25 (b) Beidelita_300. Las posiciones
de Bragg corresponden, en orden descendiente a la Beidelita_H, Beidelita_I y por
último al cuarzo.
149
Capítulo 5
Por otro lado, las fases más ordenadas de ambas muestras sí que muestran
un colapso estructural con espaciados basales de d001=10.03Å y d001=9.99Å
para la Beidelita_25 y la Beidelita_300, respectivamente.
Tabla 5.2. Coordenadas fraccionarias y ocupación de los átomos de las fases arcillosas
en la muestra Beidelita_25.
Beidelita_25
Posición Beidelita_H Beidelita_I
Átomo Wyckoff x y z occ. x y z occ.
Si1 8j 0.1480 0.1630 0.3090 0.9430 0.1790 0.1510 0.2830 0.9430
Al1 8j 0.1480 0.1630 0.3090 0.0570 0.1790 0.1510 0.2830 0.0570
6
Li 2c 0.0000 0.0000 0.5000 0.0000
7 0.0000 0.0000 0.5000 0.0094
Li 2c 0.0000 0.0000 0.5000 0.0000
Al2 4h 0.0000 0.3840 0.5000 0.2233 0.0000 0.3760 0.5000 0.2233
Fe2 4h 0.0000 0.3840 0.5000 0.0130 0.0000 0.3760 0.5000 0.0130
Mg2 4h 0.0000 0.3840 0.5000 0.0138 0.0000 0.3760 0.5000 0.0138
O1 4i 0.2180 0.0000 0.2610 0.2500 0.1070 0.0000 0.1920 0.2500
O2 8j 0.4135 0.1814 0.2443 1.0000 0.4420 0.2300 0.2360 1.0000
O3 8j 0.2010 0.2150 0.4370 1.0000 0.1950 0.2040 0.4440 1.0000
O-H 4i 0.6270 0.0000 0.3870 0.2500 0.7130 0.0000 0.4020 0.2500
6
Li 2b 0.0000 0.5000 0.0000 0.0016
7 0.0000 0.5000 0.0000 0.0000
Li 2b 0.0000 0.5000 0.0000 0.0184
151
Capítulo 5
Tabla 5.3. Coordenadas fraccionarias y ocupación de los átomos de las fases arcillosas
en la muestra Beidelita_300.
Beidelita_300
Posición Beidelita_H Beidelita_I
Átomo Wyckoff x y z occ. x y z occ.
Si1 8j 0.1480 0.1630 0.3090 0.9430 0.1790 0.1510 0.2830 0.9430
Al1 8j 0.1480 0.1630 0.3090 0.0570 0.1790 0.1510 0.2830 0.0570
6
Li 2c 0.0000 0.0000 0.5000 0.0090
7 0.0000 0.0000 0.5000 0.0094
Li 2c 0.0000 0.0000 0.5000 0.0007
Al2 4h 0.0000 0.3840 0.5000 0.2233 0.0000 0.3760 0.5000 0.2233
Fe2 4h 0.0000 0.3840 0.5000 0.0130 0.0000 0.3760 0.5000 0.0130
Mg2 4h 0.0000 0.3840 0.5000 0.0138 0.0000 0.3760 0.5000 0.0138
O1 4i 0.2180 0.0000 0.2610 0.2500 0.1070 0.0000 0.1920 0.2500
O2 8j 0.4060 0.2260 0.2620 1.0000 0.4420 0.2300 0.2360 1.0000
O3 8j 0.2010 0.2150 0.4370 1.0000 0.1950 0.2040 0.4440 1.0000
O-H 4i 0.6270 0.0000 0.3870 0.2500 0.7130 0.0000 0.4020 0.2500
6
Li 2b 0.0000 0.5000 0.0000 0.0001
7 0.0000 0.5000 0.0000 0.0000
Li 2b 0.0000 0.5000 0.0000 0.0068
5.5.2.1. Caolinita
Los resultados del refinamiento indican (Tabla 5.1) que no hay diferencias
significativas en los parámetros de celda entre la muestra Caolinita_25 y la
Caolinita_300, apoyando la hipótesis de que se trata de una estructura de
alta estabilidad. Con respecto a las coordenadas atómicas (Tabla 5.4), vemos
152
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
Figura 5.9. Refinamiento Rietveld. (a) Caolinita_25 (b) Caolinita_300. Las posiciones
de Bragg corresponden, en orden descendiente a la caolinita y al cuarzo.
153
Capítulo 5
5.6. Conclusiones
154
Estudio del fraccionamiento de isótopos de litio con difracción de neutrones de alta resolución
155
Capítulo 5
156
Capítulo 6
Los isótopos de litio como marcadores de
meteorización en Marte
Los isótopos de litio como marcadores de meteorización en Marte
160
Los isótopos de litio como marcadores de meteorización en Marte
La diferencia de masa relativa ( ̴15%) entre sus dos isótopos estables, 6Li y
7Li, conduce a un significante fraccionamiento entre el agua y las fases
sólidas, proporcionando información relevante sobre los procesos de
meteorización química [L Chan et al., 1992]. La disolución de basalto se da
de acuerdo a un mecanismo congruente, por sí sola no causa
fraccionamiento apreciable [Wimpenny et al., 2010a]; sin embargo, el
fraccionamiento isotópico tiene lugar como resultado de procesos de
intercambio iónico entre el agua y las fases solidas [Huh et al., 2004; Huh et
al., 1998; Kısakűrek et al., 2005; Taylor and Urey, 1938]. Los minerales
secundarios, en particular las arcillas, incorporan preferencialmente el
isótopo ligero en su red estructural, conduciendo al agua remanente a una
firma isotópica pesada [Misra and Froelich, 2012; Vigier et al., 2008; Zhang
et al., 1998]. El fraccionamiento isotópico también ocurre en oxihidróxidos
e hidróxidos, como la ferrihidrita o la gibbsita. Las variaciones en la firma
isotópica del litio, δ7Li, producidas por la alteración acuosa son apreciables
solo en los primeros milímetros de la superficie de roca basáltica. La δ7Li en
la superficie de la roca es más ligera que su interior, debido a la
incorporación preferencial del 6Li en arcillas u otros productos de alteración
ricos en óxidos sobre la superficie [Pistiner and Henderson, 2003].
161
Capítulo 6
162
Los isótopos de litio como marcadores de meteorización en Marte
163
Capítulo 6
164
Los isótopos de litio como marcadores de meteorización en Marte
166
Los isótopos de litio como marcadores de meteorización en Marte
167
Capítulo 6
lugares de fijación, efectos del tamaño del cristal y cinéticos, etc. [Verney-
Carron et al., 2011; Williams and Hervig, 2005]. Prácticamente, la única
información experimental disponible (y no siempre) sobre la captura de litio
en arcillas y su fraccionamiento isotópico es la variación de la capacidad de
intercambio (después de la incorporación permanente del Li a las capas
octaédricas) y el factor de fraccionamiento, α.
169
Capítulo 6
ELi+ [K + ]𝛽K+
K K/Li = [6.4]
EK+ [Li+ ]𝛽Li+
170
Los isótopos de litio como marcadores de meteorización en Marte
Los cálculos iniciales que envuelven la disolución del basalto bajo las
restricciones particulares del modelo (enfriamiento, evaporación, y pH), nos
permitió determinar la especiación iónica en disolución y el estado de
sobresaturación de fases sólidas y gases. La cantidad de sólidos que se
disuelven o precipitan en cada step es modelada cinéticamente. Las
expresiones de tasas usadas para calcular la disolución de silicatos fueron
basadas en la teoría del estado de transición [Lasaga, 2014], incluyendo
coeficientes de actividad elemental para modelar la dependencia del pH, y el
factor de activación de Arrhenius para la restricción de la influencia de la
temperatura.
dm −Ea 1
.( −
1
)
= SA [K H+ e R T 298.15 anH+
dt
−Ea 1
.( −
1
)
−Ea 1
.( −
1
) [6.6]
+ (K n e R T 298.15 K OH− e R T 298.15 am − ](Ω
OH
− 1)
dm
Donde dt
es la tasa de disolución/precipitación del mineral (mol/s-1); SA es
el área de la superficie reactiva (m2); KH+, Kn, KOH- son las constantes de tasa
específicas para la disolución en medios ácidos, neutros o básicos (mol·m-
2·s-1); R es la constante de gases; y T es la temperatura (K).
171
Capítulo 6
Ri Ci α−1
= ( ) [6.7]
R i−1 Ci−1
172
Los isótopos de litio como marcadores de meteorización en Marte
173
Capítulo 6
Óxidos e hidróxidos Ferrihidrita Carga generada por la ruptura de sus 100-700 0.998 [5]
Gibbsita bordes 100-500 0.986 [6]
[1] [Zhang et al , 1998]
[2] [Wimpenny et al , 2010]
[3] [Vigier et al., 2008]
[4] [Chan et al., 2002]
[5] [Pogge von Strandmann et al., 2010]
[6] [Pistiner and Henderson, 2003]
7
Li
(6 )
Li sample
δ7 Li = 7Li − 1 × 1000 [6.9]
(6 )
[ Li LSVEC ]
174
Los isótopos de litio como marcadores de meteorización en Marte
El efecto cooperativo entre las fases que precipitan y las que se disuelven
junto con su captura específica sobre la firma isotópica final de la disolución
restante se ilustra en la Figura 6.2. Para simplificar, se ha considerado un
sistema donde dos arcillas precipitan o se disuelven: beidelita y clorita
(ambas identificadas en Marte) [Le Deit et al., 2012]. En el primer caso (Fig.
6.2a) se asume que solo la beidelita captura isotopos de Li, mientras en el
segundo caso se permite que ambas arcillas lo hagan (Fig. 6.2b). La beidelita
sigue un comportamiento particular que involucra un proceso secuencial de
precipitación, disolución y precipitación de nuevo como función del pH del
sistema. Por el contrario, la clorita está continuamente precipitando a lo
largo de este rango de pH. Si la beidelita es la única fase que captura Li (Fig.
6.2a), la curva refleja la evolución de la tendencia de la δ7Li imitando la forma
de la curva que representa la cantidad de masa de la arcilla en el sistema,
aunque con un pequeño retraso. Esto es debido a la interacción entre el
proceso de captura y el de liberación. De hecho, incluso aunque la beidelita
175
Capítulo 6
177
Capítulo 6
178
Los isótopos de litio como marcadores de meteorización en Marte
179
Capítulo 6
180
Los isótopos de litio como marcadores de meteorización en Marte
donde la ratio entre las constantes de equilibrio fue restringida para que
fuese igual al factor de fraccionamiento dado por Zhang et al. [1998], y los
valores de log K fueron tomados de la base de datos PHREEQC [Merriam Jr
and Thomas, 1956].
181
Capítulo 6
182
Los isótopos de litio como marcadores de meteorización en Marte
183
Capítulo 6
Figura 6.3. Continuación c) y d) representan lo mismo que las Fig. 6.3a y 6.3b, pero
analizando el caso de la vermiculita. En este caso, el valor determinado
experimentalmente de la fracción restante es F=0.6636, con α=1.029±0.005.
184
Los isótopos de litio como marcadores de meteorización en Marte
185
Capítulo 6
186
Los isótopos de litio como marcadores de meteorización en Marte
a)
188
Los isótopos de litio como marcadores de meteorización en Marte
Figura 6.6. Firma isotópica del Li en agua resultante de la meteorización del basalto
en condiciones próximas al equilibrio y Tª constante.
189
Capítulo 6
porosidad [Sø et al., 2008] de 0.57 del basalto fracturado [McWorter and
Sunada, 1977]. La composición mineralógico, constantes de tasa, y
superficies reactivas empleadas en los cálculos se describen en la tabla 6.3.
El modelo se extiende a un periodo de tiempo equivalente a 4x105 años.
192
Los isótopos de litio como marcadores de meteorización en Marte
Figura 6.7. δ7Li para la caolinita, beidelita-k y clorita del modelo representado en la
figura 6.6. Dependencia de la cinética de precipitación/disolución y del pH.
193
Capítulo 6
194
Los isótopos de litio como marcadores de meteorización en Marte
195
Capítulo 6
196
Los isótopos de litio como marcadores de meteorización en Marte
197
Capítulo 6
Figura 6.9. δ7Li para la caolinita, beidelita-k y clorita del modelo representado en la
figura 6.8. Dependencia sobre la cinética de precipitación/disolución y del pH.
198
Los isótopos de litio como marcadores de meteorización en Marte
199
Capítulo 6
200
Los isótopos de litio como marcadores de meteorización en Marte
201
Capítulo 6
202
Los isótopos de litio como marcadores de meteorización en Marte
203
Capítulo 6
204
Los isótopos de litio como marcadores de meteorización en Marte
205
Capítulo 6
206
Los isótopos de litio como marcadores de meteorización en Marte
Ni el MSL ni ninguna otra misión previa a Marte han podido realizar hasta el
momento análisis de isótopos de Li. El MSL incluye un Laser-Induced
Breakdown Spectrometer (LIBS) como parte de su suite ChemCam. El LIBS es
muy sensible al Li; sin embargo, la separación isotópica de emisiones
atómicas para el Li es demasiado estrecha: de 0.015 nm a 670 nm frente a la
resolución de la ChemCam de 0.65 nm a la misma longitud de onda. Así
podría haber una pequeña posibilidad de ver un efecto isotópico si analizase
un oxido de Li puro (isótopos de B pueden ser detectados con el LIBS por
esta vía [Russo et al. 2011]), pero es imposible para el MSL. El MSL también
incluye un cromatógrafo de gas/espectrómetro de masas como parte de la
suite SAM, cuya función es detectar y caracterizar la composición isotópica
de volátiles, vaporizando muestras del suelo y de rocas mediante
calentamiento. Pero el Li no es volátil, SAM es incapaz de analizar sus
isótopos.
207
Capítulo 6
6.8. Conclusión
208
Capítulo 7
El ciclo geoquímico del litio y las “direcciones” de
meteorización en Marte
El ciclo geoquímico del litio y las “direcciones” de meteorización en Marte
211
Capítulo 7
212
El ciclo geoquímico del litio y las “direcciones” de meteorización en Marte
Figura 7.1. Evolución del océano primitivo en la Tierra según la hipótesis de Kempe y
Degens [1985]. La fijación del Ca+2 como carbonatos insolubles propiciaría una
primera etapa alcalina. Este primitivo carácter alcalino del océano sería balanceado
por la liberación de halógenos y sulfatos asociados a la desgasificación volcánica,
dando lugar al océano de Halita.
213
Capítulo 7
215
Capítulo 7
216
El ciclo geoquímico del litio y las “direcciones” de meteorización en Marte
217
Capítulo 7
sino que se usó la CREOSA específica del LIBS (Helio) usando el software C-
QuEST [Cousin et al., 2011].
Los espectros han sido normalizados al rango espectral, es decir, para cada
pixel la intensidad se ha dividido por la suma de las intensidades del rango
espectral del pico estudiado [Wiens et al., 2013], en este caso el rango es el
VNIR. Para poder cuantificar el Li se precisa realizar el ajuste, en cada
muestra, de su línea de emisión representativa, calcular el área de su pico y
a partir de ella determinar sus concentraciones. Para la deconvolución entre
las líneas de emisión tanto del elemento de interés como de las que
enmascaran su señal, se tuvo en cuenta los datos de la siguiente tabla:
Tabla 7.1. Describe las longitudes de onda de las líneas de emisión para el Li, el rango
de longitud de onda que se va a considerar y las longitudes de onda de las líneas de
emisión que van a interferir.
Elemento Longitud de onda (nm) Rango de longitud de onda (nm) Líneas de interferencia (nm)
Li 670.97 668.5 - 673.5 670.12 (Si); 671.95 (Ca)
218
El ciclo geoquímico del litio y las “direcciones” de meteorización en Marte
Coronation
% en peso
SiO2 45.80
TiO2 0.97
Al2O3 11.60
FeOT 17.90
MgO 4.00
CaO 8.40
Na2O 3.79
K2O 0.66
ppm ± σ
Li 16.20 ± 0.92
Pearson2
% en peso
SiO2 44.06
TiO2 1.65
Al2O3 8.53
FeOT 24.84
MgO 3.21
CaO 5.22
Na2O 3.19
K2O 1.21
ppm ± σ
Li 17.81 ± 1.54
John Klein
% en peso
SiO2 46.63
TiO2 0.94
Al2O3 9.87
FeOT 18.73
MgO 10.53
CaO 5.84
Na2O 2.74
K2O 0.37
ppm ± σ
Li 8.29 ± 0.54
219
Capítulo 7
Cumberland
% en peso
SiO2 46.82
TiO2 1.00
Al2O3 9.44
FeOT 19.72
MgO 8.79
CaO 5.40
Na2O 2.72
K2O 0.37
ppm ± σ
Li 10.62 ± 0.72
Windjana
% en peso
SiO2 45.43
TiO2 0.84
Al2O3 7.33
FeOT 17.97
MgO 9.48
CaO 4.79
Na2O 1.27
K2O 3.08
ppm ± σ
Li 8.51 ± 0.48
Confidence
% en peso
SiO2 49.55
TiO2 0.97
Al2O3 12.76
FeOT 19.93
MgO 3.57
CaO 3.98
Na2O 2.69
K2O 0.94
ppm ± σ
Li 14.12 ± 1.10
Mojave
% en peso
SiO2 53.80
TiO2 0.98
Al2O3 16.20
FeOT 14.80
MgO 4.53
CaO 4.50
Na2O 3.43
K2O 0.81
ppm ± σ
Li 19.33 ± 1.34
220
El ciclo geoquímico del litio y las “direcciones” de meteorización en Marte
Telegraph
% en peso
SiO2 43.75
TiO2 1.11
Al2O3 11.48
FeOT 18.52
MgO 3.68
CaO 8.15
Na2O 2.61
K2O 0.63
ppm ± σ
Li 9.79 ± 0.55
Gobabeb
% en peso
SiO2 43.47
TiO2 0.94
Al2O3 8.68
FeOT 19.66
MgO 6.85
CaO 7.64
Na2O 1.99
K2O 0.37
ppm ± σ
Li 11.65 ± 0.92
Oudam
% en peso
SiO2 50.33
TiO2 0.90
Al2O3 9.07
FeOT 19.14
MgO 5.47
CaO 3.06
Na2O 1.76
K2O 0.67
ppm ± σ
Li 5.91 ± 0.94
Marimba
% en peso
SiO2 48.91
TiO2 0.91
Al2O3 9.28
FeOT 19.34
MgO 4.69
CaO 4.84
Na2O 2.15
K2O 1.15
ppm ± σ
Li 13.19 ± 0.97
221
Capítulo 7
Quela
% en peso
SiO2 45.45
TiO2 0.85
Al2O3 9.23
FeOT 19.56
MgO 4.49
CaO 6.22
Na2O 2.40
K2O 0.67
ppm ± σ
Li 9.66 ± 0.81
Sebina
% en peso
SiO2 42.43
TiO2 0.83
Al2O3 8.47
FeOT 17.24
MgO 4.47
CaO 10.54
Na2O 1.62
K2O 0.62
ppm ± σ
Li 20.07 ± 1.15
222
El ciclo geoquímico del litio y las “direcciones” de meteorización en Marte
Las tendencias de los elementos mayoritarios del sedimento (Si, Fe, Mg, Ca,
Al, Na), su pérdida o ignición, nos proporciona una herramienta para
determinar el grado de meteorización que han sufrido las rocas basálticas a
lo largo de la trayectoria del rover. La interpretación se va a basar tanto en
el índice de alteración química (IAQ) como en el índice de alteración máfica
(IAM).
223
Capítulo 7
𝐴𝑙2 𝑂3
𝐼𝐴𝑄 (%) = [ ] × 100 [7.2]
(𝐴𝑙2 𝑂3 + 𝑁𝑎2 𝑂3 + 𝐶𝑎𝑂 + 𝐾2 𝑂)
224
El ciclo geoquímico del litio y las “direcciones” de meteorización en Marte
Figura 7.5. a) Representación gráfica de los valores IAQ de las distintas localizaciones
en un diagrama ternario A-CN-K. Tendencias de meteorización.
Por otro lado, las rocas según van sufriendo meteorización van presentando
un aumento en K [Fedo et al., 1995], lo que se reflejaría en una tendencia de
los datos hacia la esquina del K2O. En las muestras analizadas y
representadas en el gráfico A-CN-K se observa un enriquecimiento en K en
la muestra de Windjana, en comparación al resto de localizaciones, este no
se debe a que haya sufrido una mayor meteorización, sino a que posee un
mayor porcentaje de feldespato potásico (21.1% [Treiman et al., 2016]) en
comparación con las otras localizaciones (1-9%), determinados por
difracción de rayos-X, CheMin. La laterización caracteriza las últimas fases
de la meteorización. La meteorización química conduce a la formación de
arcillas (caolinita) y óxidos de hierro y aluminio. Como muestra la fig. 7.6, en
todas las localizaciones se observan índices de laterización similares y
225
Capítulo 7
226
El ciclo geoquímico del litio y las “direcciones” de meteorización en Marte
(𝐴𝑙2 𝑂3 )
𝐼𝐴𝑀𝑟𝑒𝑑𝑢𝑐𝑡𝑜𝑟 = 100 × [ ] [7.2]
𝐴𝑙2 𝑂3 + 𝐹𝑒𝑂(𝑇) + 𝑀𝑔𝑂 + 𝐶𝑎𝑂 + 𝑁𝑎2 𝑂 + 𝐾2 𝑂
227
Capítulo 7
(𝐴𝑙2 𝑂3 + 𝐹𝑒𝑂(𝑇) )
𝐼𝐴𝑀ó𝑥𝑖𝑐𝑜 = 100 × [ ] [7.3]
𝐴𝑙2 𝑂3 + 𝐹𝑒𝑂(𝑇) + 𝑀𝑔𝑂 + 𝐶𝑎𝑂 + 𝑁𝑎2 𝑂 + 𝐾2 𝑂
228
El ciclo geoquímico del litio y las “direcciones” de meteorización en Marte
Por último los diagramas A-L-F [Nesbitt and Wilson, 1992] representan la
tendencia simultánea de todos los elementos móviles de la muestra, Ca, Na,
K y Mg (Fig. 7.8). El vector de meteorización global sale de la esquina L, que
representa el contenido total de los elementos móviles, y según se aleja
mayor lixiviación y grado de meteorización. Los valores máximos del
IAMÓxico se representarían en el eje A-F.
229
Capítulo 7
elemente inmóvil en las muestras del cráter, pasando de formar parte del
mineral primario a uno secundario. La presencia de óxidos de Fe,
principalmente magnetita en las muestras de la base y hematite en las de
menor profundidad, junto con los minerales glauconíticos, puede ser la
causa de esta tendencia. La desviación de la muestra de Sebina hacia valores
de menor meteorización y perdida de Ca, Na y K, en comparación con
Marimba y Quela de composición similar, se debe a una gran fluctuación de
los datos tomados en cada spot, que puede ser debido a la coincidencia de
ellos con fases de sulfato de Ca (identificadas por DRX, capítulo 4).
231
Capítulo 7
Figura 7.10. a) Evolución del pH, δ7Li, y masa de agua evaporada a lo largo del
tiempo. b) Relación de la firma isotópica con respecto al sumatorio de los minerales
primarios y secundarios.
233
Capítulo 7
7.5. Conclusiones
235
Capítulo 7
236
Capítulo 8
Conclusiones
Conclusión
239
Conclusión
Se han combinado los análisis elementales del LIBS con datos de análisis
mineralógicos cuantitativos para trazar trayectorias de meteorización, a
partir de la desviación composicional/mineralógica. Los resultados sugieren
dos etapas en el proceso de meteorización: (i) una etapa inicial donde la
rápida disolución de silicatos primarios mantuvo el pH en condiciones
neutras a ligeramente alcalinas, dando como resultado la formación de
arcillas. (ii) Y, una evolución ácida, inducida por la existencia de volátiles
unida a un proceso de progresiva evaporación a baja temperatura.
241
Capítulo 8
Referencias bibliográficas
Referencias bibliográficas
Appelo, C., A. Vinsot, S. Mettler, and S. Wechner (2008), Obtaining the porewater
composition of a clay rock by modeling the in-and out-diffusion of anions and
cations from an in-situ experiment, Journal of Contaminant Hydrology, 101(1-4), 67-
76.
Beck, P., J. Barrat, M. Chaussidon, P. Gillet, and M. Bohn (2004), Li isotopic variations
in single pyroxenes from the Northwest Africa 480 shergottite (NWA 480): a record
of degassing of Martian magmas? 1, Geochimica et Cosmochimica Acta, 68(13),
2925-2933.
Bentor, Y. K., and M. Kastner (1965), Notes on the mineralogy and origin of
glauconite, Journal of Sedimentary Research, 35(1), 155-166.
Berner, R. A., and K. Caldeira (1997), The need for mass balance and feedback in the
geochemical carbon cycle, Geology, 25(10), 955-956.
245
Referencias bibliográficas
Bish, D., and R. Von Dreele (1989), Rietveld refinement of non-hydrogen atomic
positions in kaolinite, Clays and Clay Minerals, 37(4), 289-296.
Brant, C., L. Coogan, K. Gillis, W. Seyfried, N. Pester, and J. Spence (2012), Lithium
and Li-isotopes in young altered upper oceanic crust from the East Pacific Rise,
Geochimica et Cosmochimica Acta, 96, 272-293.
Bristow, T. F., et al. (2018), Clay mineral diversity and abundance in sedimentary
rocks of Gale crater, Mars, Science Advances, 4(6).
Brown, and Brindley (1980), Crystal structures of clay minerals and their X-ray
identification, Mineralogical Society London.
246
Referencias bibliográficas
Burton, K., and N. Vigier (2012), Lithium Isotopes as Tracers in Marine and
Terrestrial Environments, in Handbook of Environmental Isotope Geochemistry,
edited, pp. 41-59, Springer.
Carroll, D. (1959), Ion exchange in clays and other minerals, Geological Society of
America Bulletin, 70(6), 749-779.
Carter, J., D. Loizeau, N. Mangold, F. Poulet, and J.-P. Bibring (2015), Widespread
surface weathering on early Mars: A case for a warmer and wetter climate, Icarus,
248(0), 373-382.
Carter, J., F. Poulet, J.-P. Bibring, and S. Murchie (2010), Detection of hydrated
silicates in crustal outcrops in the northern plains of Mars, Science, 328(5986),
1682-1686.
Catling, D. C. (2007), Mars: Ancient fingerprints in the clay, Nature, 448(7149), 31-
32.
Chan, L. H., and F. A. Frey (2003), Lithium isotope geochemistry of the Hawaiian
plume: results from the Hawaii Scientific Drilling Project and Koolau volcano,
Geochemistry, Geophysics, Geosystems, 4(3).
247
Referencias bibliográficas
Chan, L.-H., J. C. Alt, and D. A. Teagle (2002), Lithium and lithium isotope profiles
through the upper oceanic crust: a study of seawater–basalt exchange at ODP Sites
504B and 896A, Earth and Planetary Science Letters, 201(1), 187-201.
Charlton, S. R., and D. L. Parkhurst (2011), Modules based on the geochemical model
PHREEQC for use in scripting and programming languages, Computers &
Geosciences, 37(10), 1653-1663.
Chevrier, V., F. Poulet, and J.-P. Bibring (2007), Early geochemical environment of
Mars as determined from thermodynamics of phyllosilicates, Nature, 448(7149),
60-63.
248
Referencias bibliográficas
Dauphas, N., F.-Z. Teng, and N. T. Arndt (2010), Magnesium and iron isotopes in 2.7
Ga Alexo komatiites: mantle signatures, no evidence for Soret diffusion, and
identification of diffusive transport in zoned olivine, Geochimica et Cosmochimica
Acta, 74(11), 3274-3291.
Deit, L., H. Ernst, F. Frank, P. Monica, R. A. Pio, and J. Ralf (2013), Sequence of infilling
events in Gale Crater, Mars: Results from morphology, stratigraphy, and mineralogy,
Journal of Geophysical Research: Planets, 118(12), 2439-2473.
249
Referencias bibliográficas
illite-smectites and other interstratified minerals, paper presented at Clays and Clay
Minerals, Citeseer.
Eberl, D. (1978), Reaction series for dioctahedral smectites, Clays and Clay Minerals,
26(5), 327-340.
Eberl, D. (1984), Clay mineral formation and transformation in rocks and soils, Phil.
Trans. R. Soc. Lond. A, 311(1517), 241-257.
250
Referencias bibliográficas
Fassett, C. I., and J. W. Head (2008), Valley network-fed, open-basin lakes on Mars:
Distribution and implications for Noachian surface and subsurface hydrology,
Icarus, 198(1), 37-56.
Fassett, C. I., and J. W. Head (2011), Sequence and timing of conditions on early Mars,
Icarus, 211(2), 1204-1214.
Fedo, C. M., H. Wayne Nesbitt, and G. M. Young (1995), Unraveling the effects of
potassium metasomatism in sedimentary rocks and paleosols, with implications for
paleoweathering conditions and provenance, Geology, 23(10), 921-924.
Flesch, G., A. Anderson Jr, and H. Svec (1973), A secondary isotopic standard for
6Li/7Li determinations, International Journal of Mass Spectrometry and Ion
Physics, 12(3), 265-272.
Gaines Jr, G. L., and H. C. Thomas (1953), Adsorption studies on clay minerals. II. A
formulation of the thermodynamics of exchange adsorption, The Journal of Chemical
Physics, 21(4), 714-718.
251
Referencias bibliográficas
Gellert, R., J. Campbell, P. King, L. Leshin, G. Lugmair, J. Spray, S. Squyres, and A. Yen
(2009), The alpha-particle-x-ray-spectrometer (apxs) for the mars science
laboratory (msl) rover mission, paper presented at The 40th Lunar and Planetary
Science Conference, LPI Contribution.
252
Referencias bibliográficas
Halevy, I., M. T. Zuber, and D. P. Schrag (2007), A sulfur dioxide climate feedback on
early Mars, Science, 318(5858), 1903-1907.
253
Referencias bibliográficas
Hill, R., and C. Howard (1987), Quantitative phase analysis from neutron powder
diffraction data using the Rietveld method, Journal of Applied Crystallography,
20(6), 467-474.
Hofmann, U., and R. Klemen (1950), Verlust der austauschfähigkeit von lithiumionen
an bentonit durch erhitzung, Zeitschrift für anorganische Chemie, 262(1‐5), 95-99.
Hower, J. (1961), Some factors concerning the nature and origin of glauconite,
American Mineralogist: Journal of Earth and Planetary Materials, 46(3-4_Part_1),
313-334.
Hoyer, M., N.-A. Kummer, and B. Merkel (2015), Sorption of Lithium on Bentonite,
Kaolin and Zeolite, Geosciences, 5(2), 127.
Huh, Y., L.-H. Chan, and O. A. Chadwick (2004), Behavior of lithium and its isotopes
during weathering of Hawaiian basalt, Geochemistry, Geophysics, Geosystems, 5(9),
Q09002.
Huh, Y., L.-H. Chan, L. Zhang, and J. M. Edmond (1998), Lithium and its isotopes in
major world rivers: implications for weathering and the oceanic budget, Geochimica
et Cosmochimica Acta, 62(12), 2039-2051.
Hynek, B. M., M. Beach, and M. R. Hoke (2010), Updated global map of Martian valley
networks and implications for climate and hydrologic processes, Journal of
Geophysical Research: Planets, 115(E9).
Jaynes, W., and J. Bigham (1987), Charge reduction, octahedral charge, and lithium
retention in heated, Li-saturated smectites, Clays and Clay Minerals, 35(6), 440-448.
254
Referencias bibliográficas
Kempe, S., and E. T. Degens (1985), An early soda ocean?, Chemical Geology, 53(1-
2), 95-108.
Kim, J., U. G. Nielsen, and C. P. Grey (2008), Local environments and lithium
adsorption on the iron oxyhydroxides lepidocrocite (γ-FeOOH) and goethite (α-
FeOOH): A 2H and 7Li solid-state MAS NMR study, Journal of the American Chemical
Society, 130(4), 1285-1295.
Kısakűrek, B., R. H. James, and N. B. Harris (2005), Li and δ7Li in Himalayan rivers:
Proxies for silicate weathering?, Earth and Planetary Science Letters, 237(3-4), 387-
401.
Lasaga, A. C. (2014), Kinetic theory in the earth sciences, Princeton university press.
Magna, T., U. Wiechert, and A. N. Halliday (2006), New constraints on the lithium
isotope compositions of the Moon and terrestrial planets, Earth and Planetary
Science Letters, 243(3–4), 336-353.
255
Referencias bibliográficas
Malin, M. C., and K. S. Edgett (2003), Evidence for Persistent Flow and Aqueous
Sedimentation on Early Mars, Science, 302(5652), 1931-1934.
Merriam Jr, C. N., and H. C. Thomas (1956), Adsorption studies on clay minerals. VI.
Alkali ions on attapulgite, The Journal of Chemical Physics, 24(5), 993-995.
Meunier, A., P. Sardini, J. Robinet, and D. Prêt (2007), The petrography of weathering
processes: facts and outlooks, Clay Minerals, 42(4), 415-435.
256
Referencias bibliográficas
of Mars: A potential landing site for the Mars Science Laboratory (MSL) mission,
Astrobiology, 10(7), 687-703.
Millot, R., N. Vigier, and J. Gaillardet (2010), Behaviour of lithium and its isotopes
during weathering in the Mackenzie Basin, Canada, Geochimica et Cosmochimica
Acta, 74(14), 3897-3912.
Misra, S., and P. N. Froelich (2012), Lithium isotope history of Cenozoic seawater:
changes in silicate weathering and reverse weathering, science, 335(6070), 818-
823.
Moore, D. M., and R. C. Reynolds (1989), X-ray Diffraction and the Identification and
Analysis of Clay Minerals, Oxford university press Oxford.
Morris, R. V., et al. (2016), Silicic volcanism on Mars evidenced by tridymite in high-
SiO<sub>2</sub> sedimentary rock at Gale crater, Proceedings of the National
Academy of Sciences, 113(26), 7071-7076.
257
Referencias bibliográficas
Murray, H. H. (2006), Chapter 2 Structure and Composition of the Clay Minerals and
their Physical and Chemical Properties, in Developments in Clay Science, edited by
H. H. Murray, pp. 7-31, Elsevier.
Nesbitt, H., and R. Wilson (1992), Recent chemical weathering of basalts, American
Journal of Science, 292(10), 740-777.
Odin, G. S., and P. Fullagar (1988), Chapter C4 geological significance of the glaucony
facies, in Developments in sedimentology, edited, pp. 295-332, Elsevier.
Parkhurst, D. L., and C. Appelo (1999), User's guide to PHREEQC (Version 2): A
computer program for speciation, batch-reaction, one-dimensional transport, and
inverse geochemical calculations.
258
Referencias bibliográficas
Parkhurst, D. L., and C. Appelo (2013), Description of input and examples for
PHREEQC version 3: a computer program for speciation, batch-reaction, one-
dimensional transport, and inverse geochemical calculationsRep. 2328-7055, US
Geological Survey.
Post, J., B. Cupp, and F. Madsen (1997), Beidellite and associated clays from the
DeLamar Mine and Florida Mountain area, Idaho, Clays and Clay Minerals, 45(2),
240-250.
259
Referencias bibliográficas
mudstone succession from the Murray formation, Gale crater, Mars, Earth and
Planetary Science Letters, 471, 172-185.
Rietveld (1969), A profile refinement method for nuclear and magnetic structures,
Journal of Applied Crystallography, 2(2), 65-71.
Ryu, J.-S., N. Vigier, S.-W. Lee, K.-S. Lee, and O. A. Chadwick (2014), Variation of
lithium isotope geochemistry during basalt weathering and secondary mineral
transformations in Hawaii, Geochimica et Cosmochimica Acta, 145, 103-115.
Ryu, J.-S., N. Vigier, S.-W. Lee, K.-S. Lee, and O. A. Chadwick (2014), Variation of
lithium isotope geochemistry during basalt weathering and secondary mineral
transformations in Hawaii, Geochimica et Cosmochimica Acta, 145, 103-115.
260
Referencias bibliográficas
Schikorr, G. (1929), On the reactions between iron, its hydroxides and water, Z.
Elektrochem, 35, 65-70.
Sears, V. F. (1992), Neutron scattering lengths and cross sections, Neutron news,
3(3), 26-37.
Seitz, H.-M., and A. B. Woodland (2000), The distribution of lithium in peridotitic and
pyroxenitic mantle lithologies—an indicator of magmatic and metasomatic
processes, Chemical Geology, 166(1-2), 47-64.
Seitz, H.-M., G. P. Brey, J. Zipfel, U. Ott, S. Weyer, S. Durali, and S. Weinbruch (2007),
Lithium isotope composition of ordinary and carbonaceous chondrites, and
differentiated planetary bodies: bulk solar system and solar reservoirs, Earth and
Planetary Science Letters, 260(3-4), 582-596.
Seitz, H.-M., G. P. Brey, S. Weyer, S. Durali, U. Ott, C. Münker, and K. Mezger (2006b),
Lithium isotope compositions of Martian and lunar reservoirs, Earth and Planetary
Science Letters, 245(1–2), 6-18.
Sø, H. U., D. Postma, R. Jakobsen, and F. Larsen (2008), Sorption and desorption of
arsenate and arsenite on calcite, Geochimica et Cosmochimica Acta, 72(24), 5871-
5884.
261
Referencias bibliográficas
Steefel, C. I., S. Carroll, P. Zhao, and S. Roberts (2003), Cesium migration in Hanford
sediment: a multisite cation exchange model based on laboratory transport
experiments, Journal of Contaminant Hydrology, 67(1–4), 219-246.
Stefá nsson, A., S. R. Gı́slason, and S. Arnó rsson (2001), Dissolution of primary
minerals in natural waters: II. Mineral saturation state, Chemical Geology, 172(3-4),
251-276.
Tang, Y.-J., H.-F. Zhang, and J.-F. Ying (2007a), Review of the lithium isotope system
as a geochemical tracer, International Geology Review, 49(4), 374-388.
Tang, Y.-J., H.-F. Zhang, E. Nakamura, T. Moriguti, K. Kobayashi, and J.-F. Ying
(2007b), Lithium isotopic systematics of peridotite xenoliths from Hannuoba, North
China Craton: Implications for melt–rock interaction in the considerably thinned
lithospheric mantle, Geochimica et Cosmochimica Acta, 71(17), 4327-4341.
Tardy, Y., J. Duplay, and B. Fritz (1987), Stability fields of smectites and illites as a
function of temperature and chemical compositionRep., Swedish Nuclear Fuel and
Waste Management Co.
Taylor, S. R., and S. McLennan (2009), Planetary crusts: their composition, origin and
evolution, Cambridge University Press.
262
Referencias bibliográficas
Taylor, T. I., and H. C. Urey (1938), Fractionation of the lithium and potassium
isotopes by chemical exchange with zeolites, The Journal of Chemical Physics, 6(8),
429-438.
Teng, F.-Z., W.-Y. Li, R. L. Rudnick, and L. R. Gardner (2010), Contrasting lithium and
magnesium isotope fractionation during continental weathering, Earth and
Planetary Science Letters, 300(1-2), 63-71.
Tertre, E., F. Hubert, S. Bruzac, M. Pacreau, E. Ferrage, and D. Prêt (2013), Ion-
exchange reactions on clay minerals coupled with advection/dispersion processes.
Application to Na+/Ca2+ exchange on vermiculite: Reactive-transport modeling,
batch and stirred flow-through reactor experiments, Geochimica et Cosmochimica
Acta, 112(0), 1-19.
Theng, B., S. Hayashi, M. Soma, and H. Seyama (1997), Nuclear magnetic resonance
and X-ray photoelectron spectroscopic investigation of lithium migration in
montmorillonite, Clays and Clay Minerals, 45(5), 718-723.
Thompson, G. R., and J. Hower (1975), Mineralogy of glauconite, Clays and Clay
Minerals, 23(4), 289-300.
Treacy, M., J. Newsam, and M. Deem (1991), A general recursion method for
calculating diffracted intensities from crystals containing planar faults, Proc. R. Soc.
Lond. A, 433(1889), 499-520.
263
Referencias bibliográficas
264
Referencias bibliográficas
Walker, P. Hays, and J. F. Kasting (1981), A negative feedback mechanism for the
long‐term stabilization of Earth's surface temperature, Journal of Geophysical
Research: Oceans, 86(C10), 9776-9782.
Warner, M. Sowe, S. Gupta, A. Dumke, and K. Goddard (2013), Fill and spill of giant
lakes in the eastern Valles Marineris region of Mars, Geology, 41(6), 675-678.
Williams, D., and K. Williams (1978), Electrophoresis and zeta potential of kaolinite,
Journal of Colloid and Interface Science, 65(1), 79-87.
Williams, L. B., and R. L. Hervig (2005), Lithium and boron isotopes in illite-smectite:
The importance of crystal size, Geochimica et Cosmochimica Acta, 69(24), 5705-
5716.
265
Referencias bibliográficas
Zhang, L., L.-H. Chan, and J. M. Gieskes (1998), Lithium isotope geochemistry of pore
waters from Ocean Drilling Program Sites 918 and 919, Irminger Basin, Geochimica
et Cosmochimica Acta, 62(14), 2437-2450.
Zhou, X., D. Liu, H. Bu, L. Deng, H. Liu, P. Yuan, P. Du, and H. Song (2018), XRD-based
quantitative analysis of clay minerals using reference intensity ratios, mineral
intensity factors, Rietveld, and full pattern summation methods: A critical review,
Solid Earth Sciences.
266
Capítulo 8
Anexos
Anexos
Índice
269
Anexos
ExCHANGE_MASTER_SPECIES
T T-
R R-
M M-
G G-
ExCHANGE_SPECIES
T- = T-
log_k 0.0
#.........................Definición de todas les especies de intercambio
PHASES
Fix_pH
H+ = H+
log_k 0
Fix_pe
e- = e-
log_k 0.0
# …………………Definición de fases con dos isotopos de Litio
Basalt (am)
SiAl0.36Fe0.19Mg0.28Ca0.26Na0.08K0.007Li0.001O3.31 + 2.62H+ + 0.69H2O =
0.36Al+3 +0.26Ca+2 +0.19Fe+2 + 2H2O+SiO2 +0.007K+ +0.001Li+ + 0.28Mg+2
+0.08Na+ + 0.008e-
log_k 2.23
-delta_H -51.039 kJ/mol # Calculated enthalpy of reaction Anorthite
Albite_Li
Na0.95Li0.05AlSi3O8 +4.0000 H+ = + 1.0000 Al+++ + 0.95Na+ + 0.05Li+ + 2.0000
H2O + 3.0000 SiO2
log_k 2.7645
-delta_H-51.8523 kJ/mol
# Calculated enthalpy of reaction Albite
# Enthalpy of formation: -939.68 kcal/mol
-analytic -1.1694e+001 1.4429e-002 1.3784e+004 -7.2866e+000 -
1.6136e+006
# -Range: 0-300
Anorthite_li
Ca0.5K0.9995Li0.0005Al2(SiO4)2 +8.0000 H+ = + 0.5Ca++ +0.9995K+ +
0.0005Li+ + 2.0000 Al+++ + 2.0000 SiO2 + 4.0000 H2O
log_k 26.5780
-delta_H -303.039 kJ/mol
# Calculated enthalpy of reaction Anorthite
270
Anexos
RATES
Albite_Li
-start
1 rem specific rate from Sverdrup, 1990, in kmol/m2/s
2 rem parm(1) = 10 * (A/V, 1/dm) (recalc's sp. rate to mol/kgw)
3 rem parm(2) = corrects for field rate relative to lab rate
4 rem temp corr: from p. 162. E (kJ/mol) / R / 2.303 = H in H*(1/T-1/298)
10 dif_temp = 1/TK - 1/298
20 pk_H = 12.5 + 3359 * dif_temp
30 pk_w = 14.8 + 2648 * dif_temp
40 pk_OH = 13.7 + 3359 * dif_temp
#41 rem ^12.9 in Sverdrup, but larger than for oligoclase...
50 pk_CO2 = 14.0 + 1677 * dif_temp
#60 pk_org = 12.5 + 1254 * dif_temp # ...rate increase for DOC
70 rate = 10^-pk_H * ACT("H+")^0.5 + 10^-pk_w + 10^-pk_OH * ACT("OH-")^0.3
71 rate = rate + 10^-pk_CO2 * (10^SI("CO2(g)"))^0.6
#72 rate = rate + 10^-pk_org * TOT("Doc")^0.4
80 moles = parm(1) * parm(2) * rate * (1 - SR("Albite_Li")) * time
81 rem decrease rate on precipitation
90 if SR("Albite") > 1 then moles = moles * 0.0000001
100 save moles
110 put(moles,30)
-end
Anorthite_Li
-start
10 dif_temp = 1/TK - 1/298
20 pk_H = 3.50 + 897 * dif_temp
30 pk_W = 9.12 + 930 * dif_temp
271
Anexos
272
Anexos
273
Anexos
274
Anexos
Al 6.381e-003
Ca 6.560e-005
Fe 2.666e-007
K 3.000e-005
Mg 3.145e-003
Na 3.000e-005
Si 8.721e-006
S 3.690e-05
Cl 1.2e-5
water 10
save solution 1
end
use solution 1
INCREMENTAL_REACTIONS True
REACTION 1
CO2(g) 0.3
# H2S(g) 0.00005
# Cl2(g) 0.0000005
# SO2(g) 0.000005
H2O -1
50 moles in 20000 steps
KINETICS 1
Beidellite-K
-formula Beidellite-K 1
-m 0
-m0 0
-tol 1e-011
Anorthite_Li
-formula Anorthite_Li 1
275
Anexos
-m 8.264
-m0 8.264
-parms 12680 0.67
-tol 1e-011
Enstatite
-formula Enstatite 1
-m 2.434
-m0 2.434
-parms 3562 0.67
-tol 1e-011
Fayalite
-formula Fayalite 1
-m 2.431
-m0 2.431
-parms 1225 0.67
-tol 1e-011
Forsterite
-formula Forsterite 1
-m 1.465
-m0 1.465
-parms 1644 0.67
-tol 1e-011
Chamosite-7A
-formula Chamosite-7A 1
-m 0
-m0 0
-tol 1e-011
Clinochlore-14A
-formula Clinochlore-14A 1
-m 0
-m0 0
-tol 1e-011
Beidellite-Ca
-formula Beidellite-Ca 1
-m 0
-m0 0
-tol 1e-011
Smectite-high-Fe-Mg
-formula Smectite-high-Fe-Mg 1
-m 0
-m0 0
-tol 1e-011
Kaolinite
276
Anexos
-formula Kaolinite 1
-m 0
-m0 0
-tol 1e-008
Ice
-formula Ice 1
-m 0
-m0 0
-tol 1e-008
Brucite
-formula Mg(OH)2 1
-m 0
-m0 0
-tol 1e-008
Nontronite-Mg
-formula Nontronite-Mg 1
-m 0
-m0 0
-tol 1e-008
Gibbsite
-formula Al(OH)3 1
-m 0
-m0 0
-tol 1e-008
Glauconite
-formula KFe3AlSi3O10(OH)2 1
-m 0
-m0 0
-tol 1e-008
Li_OH
-formula Li(OH):H2O 1
-m 0
-m0 0
-tol 1e-008
Zabuyelite
-formula Li2CO3 1
-m 0
-m0 0
-tol 1e-008
Li_sulfate
-formula Li2(SO4) 1
-m 0
-m0 0
277
Anexos
-tol 1e-008
Montmor-Ca
-formula Ca.165Mg.33Al1.67Si4O10(OH)2 1
-m 0
-m0 0
-tol 1e-008
CO2(g)
-formula CO2(g) 1
-m 0
-m0 0
-tol 1e-008
-steps 31536000000 in 20000 steps # seconds
-step_divide 100
-runge_kutta 6
-bad_step_max 500
-cvode true
-cvode_steps 50
-cvode_order 2
KNOBS
-iterations 100
-convergence_tolerance 1e-012
-tolerance 1e-012
-step_size 100
-pe_step_size 10
-diagonal_scale false
-delay_mass_water false
USER_PUNCH
-headings Li+ TOtLI ads1 ads2 ads3 ads4 ads5 adstot Li6w li7w DLi7W L6c1 l7c1
D7clay1 L6c2 l7c2 D7clay2 L6c3 l7c3 D7clay3 L6c4 l7c4 D7clay4 L6c5 l7c5 D7clay5
ttime Pevap SRCO2
-start
10 REM ------- ----fractionation constants------------------------
20 alfa1 = 1.011
30 alfa2 = 1.011
40 alfa3 = 1.011
50 alfa4 = 1.011
60 alfa5 = 1.011
70 clay1 = KIN("Beidellite-K")
80 clay2 = KIN("Kaolinite")
90 clay3 = KIN("Montmor-Ca")
100 clay4 = KIN("Clinochlore-14A")
278
Anexos
279
Anexos
280
Anexos
281
Anexos
282
Anexos
283
Anexos
# X 1e-005
KX Beidellite-K kinetic_reactant 0.00033
AlT3 Kaolinite kinetic_reactant 0.0002
CaR2 Montmor-Ca kinetic_reactant 0.000016
MgM2 Clinochlore-14A kinetic_reactant 0.0005
KG Glauconite kinetic_reactant 0.0001
-pitzer_exchange_gammas true
PRINT
SELECTED_OUTPUT
END
CALCULATE_VALUES
vtotmin
-start
3 DIM min$(5)
5 DiM Form$(5)
7 DIM Vol_min(5)
8 DIM Mol_Ini(5)
9 for i = 1 to 5 step 1
10 min$(1) = "Anorthite"
20 min$(2) = "Forsterite"
30 min$(3) = "Fayalite"
50 min$(4) = "Diopside"
60 min$(5) = "Enstatite"
61 mol_ini(1)= 4.264
63 mol_ini(2)= 1.4
65 mol_ini(3)= 1.4
67 mol_ini(4)= 1.4
69 mol_ini(5)= 1.4
320 form$(i) = PHASE_FORMULA(min$(i))
340 gram_min(i)= mol_ini(I)*gfw(form$(i))
360 next i
-end
USER_PRINT
-start
10 DIM min$(5)
284
Anexos
20 DIM form$(5)
30 DIM v_min(5)
40 DIM mol_ini(5)
45 DIM vperc_min(5)
47 DIM rho_min(5)
48 DIM gram_min(5)
49 DIM Xs_Min(5)
50 FOR r = 1 TO 1000000 STEP 100
60 gram_b = 0
70 vb = 0
71 rock$= "basalt"
90 REM --- MINERALS AND MODAL COMPOSITION OF THE ROCK (% vol) AS THE
STARTING VALUES -------
100 nmins = 5
110 vol_basalt = 2e-5 #(m3)
120 min$(1) = "Anorthite"
130 vperc_min(1) = 50
140 min$(2) = "Forsterite"
145 vperc_min(2) = 15
150 min$(3) = "Fayalite"
155 vperc_min(3) = 15
160 min$(4) = "Diopside"
165 vperc_min(4) = 15
170 min$(5) = "Enstatite"
175 vperc_min(5) = 5
190 rho_min(1) = 2.75e6 #gr/m3
200 rho_min(2) = 3.27e6
210 rho_min(3) = 4.39e6
220 rho_min(4) = 3.40e6
230 rho_min(5) = 3.19e6
240 Xs_Min(1) = 0.5
250 Xs_Min(2) = 0.2
260 Xs_Min(3) = 0.1
270 Xs_Min(4) = 0.1
280 Xs_Min(5) = 0.1
290 FOR i = 1 TO 5 STEP 1
295 form$(i) = PHASE_FORMULA$(min$(i))
297 v_min(i) = vperc_min(i)*vol_basalt/100
299 gram_min(i) = rho_min(i)*v_min(i)
305 mol_ini(i) = gram_min(i)/GFW(form$(i))
310 gram_b = gram_b+gram_min(i)
315 vB = vB + V_min(i)
#310 vb = vb+gram_min(i)/rho_min(i)
320 rho_bas = gram_b/vb
330 NEXT i
340 w = 1e-3/(1e-3+(VB))
285
Anexos
360 rs_v_b = r
370 ssp_bas2 = (w/(1-w))*(rs_v_b/rho_bas)
372 jmin=0
378 Print "---reactive surface of"rocks"as a function of ratio S/V379
Print" "
380 print " ", "S/V="rs_v_b"(1/m)," "size(Cubes)="6*(1/rs_v_b)*1000,
"mm"
381 print" "
382 print " " " Porosity = " w , "(dimensionlees)"
390 PRINT "Sspecific"rock$"="ssp_bas2, " ( m2/g)"
400 PRINT "Sreactive"rock$"(m2/kgw)=" gram_b*ssp_bas2
440 PRINT "grams"rock$ "="gram_b "in a Kgw"
450 PRINT "vol"rock$"="vb, "(m3)"
460 PRINT "rho"rock$ "="rho_bas "(g/m3)"
461 print " "
472 Print"-----Sreactive minerals in"rocks"as a function of wetting param Xmin--"
474 print" Phase name", " Reactive surface (m2/kgw)", " Modal comp (%
Vol)" " Xmin" " moles0/kgw"
475 print"-------------------------------------------------------------------------------------------
-"
476 FOR j = 1 TO 5 STEP 1
477 s_min = gram_min(j)*ssp_bas2*xs_min(j)
478 jmin = jmin+ s_min
#479 print jmin
480 PRINT min$(j), " ", s_min " " Vperc_min(j) " " xs_min(j) " "
Mol_ini(J)
485 next j
490 NEXT r
-end
SOLUTION 0
temp 0
pH 10.5
pe -8
redox pe
units mol/kgw
density 1
Al 0.0001738
Ca 0.0002285
286
Anexos
Na 0.0001285
K 0.0001285
Fe(2) 0.0003671
Mg 0.0004268
P 4.284e-007
S(-2) 5.447e-007
Si 0.0001501 as SiO2
-water 1
REACTION 1
H2O -0.01
CO2(g) 0.01
SO2(g) 0.001
Cl2(g) 0.001
SOLUTION 1-2
temp 0
pH 10
pe -8.5
redox pe
units mol/kgw
density 1
Al 0.0001738
Ca 0.0002285
Na 0.0001285
K 0.0001285
Fe(2) 0.0003671
Mg 0.0004268
P 4.284e-007
S(-2) 5.447e-007
Si 0.0001501 as SiO2
-water 1
SOLUTION 3-20
temp 0
pH 10.5
pe -8
redox pe
units mol/kgw
density 1
Al 0.0001738
Ca 0.0002285
Na 0.0001285
K 0.0001285
287
Anexos
Fe(2) 0.0003671
Mg 0.0004268
P 4.284e-007
S(-2) 5.447e-007
Si 0.0001501 as SiO2
-water 1
KINETICS 3-20
Anorthite
-formula CaAl2Si2O8 1
-m 0.098863
-m0 0.098863
-parms 103.4 0.67
-tol 1e-08
Diopside
-formula Diopside 1
-m 0.047109
-m0 0.047109
-parms 7.6704 0.67
-tol 1e-08
Enstatite
-formula Enstatite 1
-m 0.031782
-m0 0.031782
-parms 2.3989 0.67
-tol 1e-08
Fayalite
-formula Fayalite 1
-m 0.064636
-m0 0.064636
-parms 9.9038 0.67
-tol 1e-08
Forsterite
-formula Forsterite 1
-m 0.069737
-m0 0.069737
-parms 14.754 0.67
-tol 1e-08
Jarosite-K
-formula KFe3(SO4)2(OH)6 1
-m 0
-m0 0
-tol 1e-08
288
Anexos
Gypsum
-formula CaSO4:2H2O 1
-m 0
-m0 0
-tol 1e-08
Alunite
-formula KAl3(OH)6(SO4)2 1
-m 0
-m0 0
-tol 1e-08
Beidellite-Ca
-formula Ca.165Al2.33Si3.67O10(OH)2 1
-m 0
-m0 0
-tol 1e-08
Smectite-high-Fe-Mg
-formula Ca.025Na.1K.2Fe.5Fe.2Mg1.15Al1.25Si3.5H2O12 1
-m 0
-m0 0
-tol 1e-08
Clinochlore-14A
-formula Mg5Al2Si3O10(OH)8 1
-m 0
-m0 0
-tol 1e-08
Goethite
-formula FeOOH 1
-m 0
-m0 0
-tol 1e-08
Gibbsite
-formula Al(OH)3 1
-m 0
-m0 0
-tol 1e-08
Calcite
-formula CaCO3 1
-m 0
-m0 0
-parms 1640 0.67
-tol 1e-08
-steps 1 in 1 steps # seconds
289
Anexos
-step_divide 1
-runge_kutta 3
-bad_step_max 500
-cvode true
-cvode_steps 100
-cvode_order 5
USER_PRINT
-start
2 V1_ini=9.8863e-02 *GFW("CaAl2Si2O8")/2.76
3 V2_ini=4.7109e-02*GFW("CaMgSi2O6")/3.4
4 V3_ini=3.1782e-02*GFW("MgSiO3")/3.9
5 V4_ini=6.4636e-02*GFW("Mg2SiO4")/4.39
6 V5_ini=6.9737e-02 *GFW("Fe2SiO4")/3.24
7 VT_ini=V1_ini+V2_ini+V3_ini+v4_ini+v5_ini
20 V1_t=KIN("Anorthite")*GFW("CaAl2Si2O8")/2.76
30 V2_t=KIN("Diopside")*GFW("CaMgSi2O6")/3.4
40 v3_t=KIN("Enstatite")*GFW("MgSiO3")/3.9
50 V4_t=KIN("Fayalite")*GFW("Mg2SiO4")/4.39
60 v5_t=KIN("Forsterite")*GFW("Fe2SiO4")/3.24
70 VT_t=V1_t+V2_t+V3_t+v4_t+v5_t
80 weather_fract= VT_t/VT_ini
90 print "weathered fraction in Vol (cm3)=" weather_fract
-end
TRANSPORT
-cells 20
-shifts 1000
-time_step 3153600 # seconds
-lengths 20*0.01
-dispersivities 20*0.0015
-correct_disp true
-diffusion_coefficient 4e-10
-thermal_diffusion 1 4e-10
USER_PUNCH 1
-headings weather
-start
2 V1_ini=9.8863e-02 *GFW("CaAl2Si2O8")/2.76
3 V2_ini=4.7109e-02*GFW("CaMgSi2O6")/3.4
4 V3_ini=3.1782e-02*GFW("MgSiO3")/3.9
5 V4_ini=6.4636e-02*GFW("Mg2SiO4")/4.39
6 V5_ini=6.9737e-02 *GFW("Fe2SiO4")/3.24
290
Anexos
7 VT_ini=V1_ini+V2_ini+V3_ini+v4_ini+v5_ini
20 V1_t=KIN("Anorthite")*GFW("CaAl2Si2O8")/2.76
30 V2_t=KIN("Diopside")*GFW("CaMgSi2O6")/3.4
40 v3_t=KIN("Enstatite")*GFW("MgSiO3")/3.9
50 V4_t=KIN("Fayalite")*GFW("Mg2SiO4")/4.39
60 v5_t=KIN("Forsterite")*GFW("Fe2SiO4")/3.24
70 VT_t=V1_t+V2_t+V3_t+v4_t+v5_t
80 weather_fract= VT_t/VT_ini
90 punch "weathered fraction in Vol (cm3)=" weather_fract
-end
KNOBS
-iterations 100
-convergence_tolerance 1e-11
-tolerance 1e-15
-step_size 100
-pe_step_size 10
SELECTED_OUTPUT 1
291
Anexos
A.4. Ejemplo del script para realizar ajustes de estructuras con desorden de
apilamiento, mediante el programa FAULTS o DIFFAX.
TITLE
___________________________________
STRUCTURAL
FST_CMD SEQ 6 1 2 3 1 2 3
! a b c gamma
Cell 5.358788 8.507235 10.015397 90.00
1.00 1.00 1.00 0.00
!Laue symmetry
Symm -1
!number of layer types
NLAYERS 3
!layer width
Lwidth Infinite
LAYER 1
!Layer symmetry
LSYM None
292
Anexos
293
Anexos
294
Anexos
295
Anexos
LAYER 2
!Layer symmetry
LSYM None
296
Anexos
297
Anexos
298
Anexos
LAYER 3
!Layer symmetry
LSYM None
299
Anexos
300
Anexos
301
Anexos
STACKING
302
Anexos
!stacking type
RECURSIVE
!number of layers
INFINITE
TRANSITIONS
!layer 1 to layer 1
LT 0.000339 0.000000 0.000000 1.000000
-31.00 0.00 0.00 0.00
FW 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
0.00 0.00 0.00 0.00 0.00 0.00
!layer 1 to layer 2
LT 0.999661 -0.44800 0.000000 1.000000
31.00 0.00 0.00 0.00
FW 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
0.00 0.00 0.00 0.00 0.00 0.00
!layer 1 to layer 3
LT 0.000000 0.000000 0.000000 1.000000
0.00 0.00 0.00 0.00
FW 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
0.00 0.00 0.00 0.00 0.00 0.00
!layer 2 to layer 1
LT 0.000000 0.000000 0.000000 1.000000
0.00 0.00 0.00 0.00
FW 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
0.00 0.00 0.00 0.00 0.00 0.00
!layer 2 to layer 2
LT 0.000339 0.000000 0.000000 1.000000
-31.00 0.00 0.00 0.00
FW 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
0.00 0.00 0.00 0.00 0.00 0.00
!layer 2 to layer 3
LT 0.999661 -0.36475 -0.26011 1.000000
31.00 0.00 0.00 0.00
FW 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
0.00 0.00 0.00 0.00 0.00 0.00
!layer 3 to layer 1
LT 0.999661 -0.33464 0.358920 1.000000
31.00 0.00 0.00 0.00
FW 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
303
Anexos
CALCULATION
LMA
Corrmax 30
Maxfun 2500
Tol 0.100000E-04
Nprint 0
EXPERIMENTAL
!Filename Scale factor code
File _________.dat 1.000000000 11.00
Fformat XYSIGMA
!Linear interpolation
BgrInter ________.BGR
!Number of background patterns
BgrNum 0
304
Meteorización de silicatos en Marte: el litio y
2018
sus isótopos como trazadores geoquímicos