Tasa de Transferencia Sedimentaria de La Serie de U238

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 92

UNIVERSIDAD DE CHILE

FACULTAD DE CIENCIAS FSICAS Y MATEMTICAS


DEPARTAMENTO DE GEOLOGA

TASAS DE TRANSFERENCIA SEDIMENTARIA


EN LA CUENCA DEL RO UBLE : ESTIMACIN
BASADA EN DESEQUILIBRIOS EN LA SERIE DEL URANIO-238

MEMORIA PARA OPTAR AL TTULO DE GELOGO

FRANCISCO DANIEL DEL CAMPO ROJAS

PROFESOR GUA:
MARTIN REICH MORALES

MIEMBROS DE LA COMISIN:
JACOBUS LE ROUX
SEBASTIEN CARRETIER

SANTIAGO DE CHILE
JULIO 2010
Resumen

Las series del U y del Th son tres cadenas de decaimiento radioactivo que terminan en
istopos estables de plomo, pasando por varios nucleidos radioactivos intermedios. Cualquier
sistema natural cerrado debera alcanzar despus de un cierto tiempo lo que se denomina el
equilibrio secular, que consiste en un estado estacionario donde las proporciones entre los
nucleidos de la serie son iguales a constantes conocidas. Una desviacin del equilibrio secular
(desequilibrio) se debe a algn proceso de fraccionamiento qumico o isotpico. El retorno al
equilibrio de un cierto nucleido de la serie (respecto del padre) debera ocurrir en un tiempo
proporcional a su vida media en un sistema cerrado, aunque en la naturaleza rara vez se
cumple esta condicin.
La cuantificacin de los desequilibrios en la serie del 238U en sistemas fluviales permite
acotar la escala de tiempo en que ocurren los fenmenos responsables del fraccionamiento y el
eventual retorno al equilibrio. Esta metodologa se ha aplicado en ros de larga trayectoria como
el Ganges (Granet et al., 2007) y el Amazonas (Dosseto et al., 2006a), pero no todava en ros
de trayectoria corta en un margen convergente como el chileno. Uno de los objetivos del
proyecto SESPEED, en el cual se enmarca esta memoria, es realizar esto ltimo. En esta
memoria se abord la cuantificacin de tiempos de transporte de sedimentos en la cuenca del
ro uble en Chile central. Para ello se midi las razones de actividad entre 238U, 234U, 230Th y
226
Ra (serie del 238U), en 8 muestras de sedimentos del lecho del ro estudiado.
Para caracterizar la evolucin de los sedimentos ro abajo y poder interpretar los
desequilibrios, se analiz en las muestras: (a) la composicin qumica de elementos mayores y
traza, (b) la mineraloga por difraccin de rayos X y (c) las razones isotpicas de Sr y Nd. De
forma complementaria se realiz anlisis de elementos mayores y traza en siete muestras de
agua del ro uble. Los datos qumicos sugieren que la principal fuente sedimentaria son las
rocas volcnicas andesticas a baslticas de las Formaciones Cola de Zorro, Trapa-Trapa y
Cura-Malln, habiendo tambin aportes menores del batolito grantico del curso alto del ro. Esta
interpretacin es apoyada por la composicin de tierras raras e isotpica de Sr (87Sr/86Sr) y Nd
(143Nd/144Nd) de los sedimentos, la cual es muy similar a la de rocas volcnicas de volcanes
activos geogrficamente cercanos. La composicin isotpica de Sr y Nd se mantiene dentro de
un rango pequeo si se compara con sedimentos del ro Amazonas, lo que sugiere que el
perodo de maduracin de los sedimentos es mucho ms corto en el ro uble que en el ro
Amazonas (hasta ~100 ka segn Dosseto et al., 2006a).
La evolucin qumica e isotpica de los sedimentos se interpreta como el resultado de los
fenmenos siguientes: aportes sedimentarios desde fuentes distintas a las rocas volcnicas
(batolito grantico en la cordillera, formaciones sedimentarias en la Depresin Central),
disolucin preferencial de minerales ferro-magnesianos, deposicin progresiva de partculas en
suspensin al entrar en el valle central.
El 234U est prcticamente en equilibrio secular respecto del 238U, aunque el 230Th presenta
ligeras desviaciones del equilibrio respecto del 234U, especialmente ro abajo. La razn U/Th en
los sedimentos se mantiene casi constante a lo largo del ro, y las concentraciones de U y Th en
el agua del ro son bajas (< 1 ppt y < 0,1 ppt respectivamente). Esto permite afirmar que tanto el
uranio como el torio son prcticamente inmviles en los sedimentos. La razn de actividad
(226Ra/230Th) no difiere de 1,00 en ms de 8%. Los desequilibrios son mucho menores que
aquellos medidos en los ros Ganges (Granet et al., 2007) o Amazonas (Dosseto et al., 2006a),
lo que apoya la hiptesis de tiempos de residencia sustancialmente menores. Bajo el supuesto
de aportes nulos (o prdida) de Ra entre las cuarta y la quinta muestras (tramo de 20 km), se
lleg a una tasa de transferencia sedimentaria mnima de ~14 m a-1 (20 km en ~1400 aos).

i
Rsum
Les sries de l'U et le Th sont des chanes de dsintgration radioactive qui se composent
de plusieurs nuclides radioactifs, se terminant en isotopes stables de plomb. Tout systme
ferm devrait atteindre aprs un certain temps ce qu'on appelle l'quilibre sculaire, un tat
stationnaire o les ratios des nuclides de la srie sont gales des constantes connues. Un
cart par rapport l'quilibre sculaire est d un processus de fractionnement chimique ou
isotopique. Le retour l'quilibre d'un certain nuclide dans la srie par rapport au pre devrait
se produire dans un temps proportionnel sa demi-vie dans un systme ferm, bien que dans la
nature cette dernire condition est rarement remplie.
La quantification des dsquilibres dans la srie de l'238U dans des systmes naturels
permet d'estimer le temps, dans lequel se produisent les phnomnes responsables du
fractionnement. Cette mthodologie a t applique quelques fleuves ayant un long trajet, par
exemple le Gange (Granet et al., 2007) et l'Amazone (Dosseto et al., 2006a), mais pas encore
dans des rivires courtes d'une marge convergente. Le stage dont ce rapport donne un compte
rendu, a t fait dans le cadre d'un projet de recherche sur la quantification des temps de
transport sdimentaire dans le versant ouest des Andes Centrales. Spcifiquement, les
chantillons analyss sont des sdiments provenant de huit points le long de la rivire uble au
centre du Chili. Le travail de laboratoite a t fait au Laboratoire d'Hydrologie et de Gochimie
de Strasbourg (LhyGeS), dans le cadre du projet SESPEED financ par le CNRS. Les isotopes
mesurs appartiennent la srie de l'238U (238U-234U-230Th-226Ra).
Pour caractriser l'volution des sdiments en aval et interprter les dsquilibres dans les
chantillons les analyses suivantes ont t faites: (a) les concentrations des lments majeurs
(par ICP-AES) et traces (par ICP-MS), (b) la minralogie par diffraction des rayons X et (c) les
rapports isotopiques de Sr et Nd (par MC-ICP-MS). Une analyse complmentaire des lments
majeurs et traces a t ralise dans sept chantillons d'eau de la rivire uble. La composition
chimique et la signature isotopique de Sr-Nd des sdiments suggrent que la principale source
sdimentaire sont les roches volcaniques andsitiques basaltiques des Formations Trapa-
Trapa et Cura-Mallin, ayant galement des contributions subordonnes du batholite granitique
de la chane. En normalisant les donnes des sdiments par rapport l'oxyde d'aluminium
(Al2O3) et le Thorium (Th), il y a une tendance la baisse pour le Fe, Mg, Ca, Mn, Sc, V et Co,
et l'augmentation de K, Rb et Ba, ce qui est interprt comme le rsultat de trois phnomnes:
la mtorisation des minraux mafiques, les apports provenant du batholite granitique de la
chane, et le dpt progressif des particules en suspension entrant dans la valle.
La composition isotopique du Sr (87Sr/86Sr) et Nd (143Nd/144Nd) est trs semblable celle
des roches volcaniques des volcans actifs proches gographiquement, et reste dans une petite
gamme par rapport l'volution des sdiments dans le bassin du fleuve Amazone. Ceci suggre
que la priode de maturation des sdiments est beaucoup plus courte dans la rivire uble que
dans le fleuve Amazone (jusqu' ~ 100 ka selon Dosseto et al., 2006a).
L'234U est pratiquement l'quilibre sculaire par rapport l'238U, mais le 230Th montre des
lgres dviations de l'quilibre, surtout en aval. En outre, le ratio U/Th dans les sdiments reste
peu prs constant le long de la rivire, la concentration d'U et Th dans l'eau de la rivire est
faible (<1 ppt et <0,1 ppt, respectivement). Tout cela indique que l'uranium et le thorium sont
pratiquement immobiles dans les sdiments tudis. Le rapport d'activit (226Ra/230Th) dans les
sdiments commence partir de 1,00 (quilibre sculaire) et augmente jusqu' une valeur
maximale de 1,082, en forte baisse 1,045 et reste ce niveau-ci. Aprs l'interprtation des
rsultats un taux de transfert de sdiments minimum de 14 m par an a t estim (20 km dans
environ 1400 ans).

ii
Agradecimientos

Agradezco profundamente el apoyo incondicional de toda mi familia, a lo largo de toda mi


vida, y particularmente durante mi carrera, gratitud especialmente dedicada a mis padres Cecilia
y Guillermo, as como a mi hermano Felipe y a mis abuelitos Marta, Rodolfo, Eileen y Mario.

Les doy mil gracias a mis compaeros, profesores auxiliares y amigos de Geologa, de la
facultad y de otros lados, por compartir todas esos momentos que ya suman aos, en la sala de
clases o en el patio de los naranjos, en terreno, de vacaciones y en la ciudad, aprendiendo,
rindo y creciendo juntos (la lista es larga, as que rmese de paciencia o pase al siguiente
prrafo): lvaro Contreras (don Carter), Andrs Quezada (don Ramn), Alejandra Serey, Aranzaz
Bulnes (tut), Carlos Rozas, Carmen Bustamante (penita), Carolina Guzmn (chica), Carolina
Miralles (caro), Carolina Valenzuela (caro), Constanza Urresty (cony), Claudio Vergara (ciego),
Daniel Bertin, Daniel Daz, Daniel Gonzlez, Daniela Jerez, David Bustos (guatn), Diego Garca
(carcass), Diego Salas (larry), Estban Leiva (chuck), Estban Salazar (chala), Felipe Astudillo,
Felipe Coloma (waren), Felipe Orellana (felipn), Felipe Tapia, Fernanda lvarez (las dos feas),
Fernando Lpez (fea), Fernando Poblete (pjaro), Francisco Delgado, Gabriel Orozco, Gonzalo
Huerta (chalo), Hernn Bobadilla, Hugo Ulloa (pollo), Ismael Murillo (isma), Jacqueline Puratich
(jacky), Javier lvarez, Javier Corts (buho), Javier Valenzuela (nacho), Javiera Morand, Javiera
Vidal, Joaqun Cceres, Jose Escribano, Jos Ignacio Snchez (nachiwis), Joseline Tapia, Leonardo
Baeza (leo), Luis Lpez (mito), M. Carolina Soto, M. Eliana Lorca, M. Pa Rodriguez, Mara Nquira,
Mara Vargas, Marisol Lara, Martn Jacques (mahoma), Martn Lemus, Matas Pea (ataxia), Manuel
Olivares, Mauricio Muoz, Miguel Bravo, Natalia Garrido (naty), Nicols Anguita, Oscar Benavente,
Pablo Gonzlez, Pablo Snchez, Paola Martnez, Patricia Galvez, Paula Castillo, Paula Muoz,
Paulina Cecci, Rodolfo Lagos, Rodolfo Ferrando, Rodrigo Luca, Rodrigo Otero (rotero), Sebastin
Nuez, Sergio Caldern (checho), Susana Henrquez (susy), Tania Sauma, Toms Opazo (topazo),
Valeria Ramrez (vale), Valentina Flores (vale), Violeta Tolorza (la violle) ...

Expreso mi gratitud hacia mis profesores y maestros por su paciencia, motivacin por
ensearnos y calidad humana, particularmente a mis "profes" de Geologa y Geofsica (en
orden cronolgico respecto a mi carrera): Francisco "Pancho" Herv, Jaime Campos, Francisco
Munizaga, Alfonso Rubilar, Emilio Vera, Gabriel Vargas, Csar Arriagada, Martin Reich, Diego
Morata, Luisa Pinto, Denis Legrand, Andrs Pavez, Mauricio Caldern, Jacobus Le Roux,
Miguel ngel Parada, Jos Cembrano, Brian Townley, Sbastien Carretier, Luis Lara, Reinaldo
Charrier, Sofa Rebolledo, Sergio Seplveda, Carlos Palacios y Victor Maksaev.

Agradezco la valiosa y sacrificada labor de los funcionarios del departamento de Geologa


(y de la facultad), sin ellos nada de esto sera posible: Mara Rosa, Arnoldo Quilodrn (el quilo),
Carlos Gmez (don Carlos), Rosa Flores (Rosita),

Les doy las gracias a mis profesores gua Martin Reich, Sbastien Carretier y Jacobus Le
Roux, por su amabilidad, disponibilidad y paciencia.

Agradezco a todo el equipo del Laboratorio de Hidrologa y Geoqumica de Estrasburgo


(LHyGeS: Laboratoire d'Hydrologie et de Gochimie de Strasbourg), especialmente a Franois
Chabaux, Eric Pelt y Mathieu Granet, por haberme acogido amablemente y guiarme durante mi
estada en esa institucin entre octubre de 2009 y marzo de 2010. Al IRD y el CNRS por el
financiamiento de mi pasanta en el marco del projecto SESPEED.

Muchsimas gracias a todos ustedes de corazn!!

iii
Contenidos

CAPTULO 1: Introduccin 9
1.1 Formulacin........................................................................................................................... 9
1.2 Objetivos .............................................................................................................................. 11

CAPTULO 2: Series del uranio y del torio 12


2.1 Generalidades..................................................................................................................... 12
2.2 Causas de desequilibrio en la serie del U-238 .............................................................. 15
2.3 Desequilibrios en ambientes acuosos ............................................................................ 16

CAPTULO 3: Hidrologa, morfologa y clima 19


3.1 Hidrologa ............................................................................................................................ 19
3.2 Relieve ................................................................................................................................. 19
3.3 Clima .................................................................................................................................... 21

CAPTULO 4: Marco geolgico 22


4.1 Paleozoico ........................................................................................................................... 22
4.2 Mesozoico............................................................................................................................ 22
4.3 Cenozoico............................................................................................................................ 22
4.4 Rocas intrusivas ................................................................................................................. 25

CAPTULO 5: Mtodos analticos 27


5.1 Molienda de muestras ....................................................................................................... 27
5.2 Mineraloga por difraccin de rayos X ............................................................................ 27
5.3 Manipulaciones qumicas.................................................................................................. 28
5.4 Preparacin de filamentos ................................................................................................ 30
5.5 Espectrometra de masas .................................................................................................. 32

CAPTULO 6: Resultados 39
6.1 Muestras .............................................................................................................................. 39
6.2 Elementos mayores y traza (ICP-AES, ICP-MS) .......................................................... 39
6.3 Mineraloga por difraccin de rayos X ............................................................................ 42

iv
6.4 U, Th y Ra: Razones de actividad y concentraciones (TIMS) .................................... 43
6.5 Razones isotpicas de Sr y Nd (TIMS y MC-ICPMS) .................................................. 44

CAPTULO 7: Discusin 46
7.1 Elementos mayores y traza e istopos de Sr y Nd....................................................... 46
7.2 Desequilibrios en la serie del U-238 ............................................................................... 53

CAPTULO 8: Interpretaciones y conclusiones 61


8.1 Qumica e istopos de Sr y Nd ........................................................................................ 61
8.2 Desequilibrios en la serie del U-238 ............................................................................... 62

Referencias 64

Anexo A: Grficos de elementos mayores y traza 69

Anexo B: Resinas cromatogrficas 83


B.1 Resinas intercambiadoras de iones................................................................................ 83
B.2 Resinas de extraccin....................................................................................................... 87

Anexo C: Dilucin isotpica 88

Anexo D: Clculo de razones isotpicas 91

v
Indice de Figuras

Figura 1: Vidas medias de los nucleidos de las series del uranio y del torio (Bourdon et al.,
2003)................................................................................................................................. 13

Figura 2: Relieve y red de drenaje del ro Itata. Los puntos muestreados se indican con crculos
rojo.................................................................................................................................... 20

Figura 3: Precipitacin y temperatura medias mensuales en la Estacin Meteorolgica de


Chilln............................................................................................................................... 21

Figura 4: Unidades geolgicas de la cuenca del ro Itata. Los puntos rojos son los puntos de
muestreo........................................................................................................................... 23

Figura 5: Esquema del funcionamiento de un espectrmetro de masa Thermo Scientific Triton


(TIMS)............................................................................................................................... 33

Figura 6: Esquema de las trayectorias de los iones bajo la influencia del electroimn. ............. 34

Figura 7: Elementos mayores y traza normalizados respecto a Th y Al2O3. .............................. 46

Figura 8: Perfiles de tierras raras de sedimentos (semilog) normalizados respecto a rocas


volcnicas ......................................................................................................................... 50

Figura 9: Perfiles de tierras raras de muestras del ro uble normalizados resp. a basaltos
ocenicos.......................................................................................................................... 50

Figura 10: Nd vs. 87Sr/86Sr: comparacin con sedimentos del ro Amazonas. Tomado de Viers et
al. (2008)........................................................................................................................... 51

Figura 11: Grfico Nd vs. 87Sr/86Sr. Los tringulos corresponden a rocas de volcanes activos
cercanos ........................................................................................................................... 52

Figura 12: Grfico (234U/238U) vs. distancia a lo largo del thalweg del ro uble (Itata para el
punto en 155 km). ............................................................................................................. 53

Figura 13: Grfico (230Th/234U) vs. distancia a lo largo del thalweg del ro uble (Itata para el
punto en 155 km). ............................................................................................................. 54

Figura 14: Grfico U/Th y (230Th/232Th) vs. distancia a lo largo el ro uble. .............................. 54

Figura 15: Grfico (230Th/232Th) vs. (238U/232Th) en muestras de sedimentos del ro uble........ 55

Figura 16: Grfico (226Ra/230Th) vs. distancia a lo largo del thalweg del ro uble (Itata para el
punto en 155 km) .............................................................................................................. 56

Figura 17: Grfico Ba/Th vs. distancia a lo largo del thalweg del ro uble en sedimentos. ...... 57

vi
Figura 18: Grfico Ba/Z vs. distancia a lo largo del thalweg del ro uble en muestras de agua.
.......................................................................................................................................... 59

Figura 19: Grfico (226Ra/232Th) vs. (238U/232Th) en sedimentos del ro uble............................ 60

Figura 20: Grfico (226Ra/238U) vs. (230Th/238U) en sedimentos del ro uble. ............................ 60

Figura 21: Esquema de las estructura de la resina (fuente: pgina web de EICHROM)............87

vii
Indice de Tablas

Tabla 1: Pesaje de muestras y spike, utilizando una balanza analtica con una precisin de 5 g.
.......................................................................................................................................... 28

Tabla 2: Protocolo de ataques cidos ....................................................................................... 29

Tabla 3: Procedimiento de impregnacin de filamentos. ........................................................... 31

Tabla 4: Comparacin del instrumental de espectrometra de masas ....................................... 37

Tabla 5: Concentraciones de elementos mayores en las muestras de sedimentos del ro uble


(% peso). .......................................................................................................................... 39

Tabla 6: Concentraciones de elementos mayores en las muestras de agua del Ro uble (ppm).
.......................................................................................................................................... 40

Tabla 7: Concentraciones de elementos traza en las muestras de sedimentos. ....................... 41

Tabla 8: Concentraciones de elementos traza y algunos mayores en muestras de agua del ro


uble ................................................................................................................................ 42

Tabla 9: Interpretacin de los difractogramas de rayos X en roca total de sedimentos del ro


uble. ............................................................................................................................... 42

Tabla 10: Interpretacin de difractogramas de las fracciones arcillosas (<2 m): muestras del
ro uble. .......................................................................................................................... 43

Tabla 11: Razones de actividad entre nucleidos de la serie del U-238 (238U, 234U, 230Th y 226Ra)
y 232Th ............................................................................................................................... 44

Tabla 12: Concentraciones de U, Th y 226Ra en muestras de sedimento del ro uble obtenidas


por TIMS. .......................................................................................................................... 44

Tabla 13: Razones isotpicas de Sr y Nd en las muestras de sedimentos obtenidas por TIMS e
ICP-MS. ............................................................................................................................ 45

Tabla 14: Proporciones entre La y tierras raras pesadas (Tb, Ho, Er y Tm) en sedimentos del
ro uble. .......................................................................................................................... 51

Tabla 15: Resinas usadas para la separacin qumica ............................................................. 83

Tabla 16: Clasificacin de resinas segn carga elctrica del grupo funcional ........................... 84

Tabla 17: Funcin de alguna de las resinas de extraccin de EICHROM...................................87

viii
CAPTULO 1: Introduccin

1.1 Formulacin

A medida que un macizo de roca se meteoriza progresivamente desde la superficie, la zona


expuesta a la atmsfera es fracturada y alterada, y el material es erosionado dando origen a
sedimentos. El lapso entre el comienzo de la meteorizacin de un volumen de roca, y la erosin
del mismo se denomina tiempo de meteorizacin. Considerando esta lnea temporal de base,
el tiempo de transferencia sedimentaria corresponder a la edad de la erosin.

La cuantificacin de esta escala temporal es de gran importancia no solo en reas como


la sedimentologa, geomorfologa fluvial e hidrologa, sino tambin en tectnica y evolucin de
macizos cordilleranos (Bourdon et al., 2003). Por ejemplo, en los modelos de dinmica
morfolgica de cuencas de drenaje, es uno de los parmetros esenciales (ej. Crave y Davy,
2001; Verhaar et al., 2008). Los resultados de estos modelos muestran que la evolucin
temporal de una red de drenaje depende fuertemente de las capacidades de transporte y
erosin de los ros de una cuenca. Dicha capacidad es descrita, entre otras cantidades, por un
tiempo caracterstico de transporte.

Recientemente, algunos autores han estudiado la meteorizacin y la pedognesis


(formacin de suelo) a travs de metodologas basadas en los desequilibrios de la serie del
uranio (ej. Chabaux et al., 2003; Dosseto et al., 2006a). Tanto la meteorizacin in situ, es decir,
la descomposicin qumica y fsica de rocas frescas en su lugar de afloramiento (por efectos
oxidativos y de disolucin acuosa), como la maduracin de sedimentos en presencia de agua
producen un fraccionamiento (o separacin qumica) de elementos, en particular del uranio (U)
y del torio (Th) contenidos inicialmente en los minerales constituyentes de la roca. Dichos
fraccionamientos qumicos, junto con un eventual fraccionamiento isotpico de U, generan
desequilibrios en las series del U y del Th, por ejemplo, la del uranio-238 (Bourdon et al., 2003).

Como se explica ms adelante (ver 2.1), el tiempo de regreso al equilibrio en un sistema


cerrado es proporcional a la vida media del segundo istopo de la cadena de decaimiento (234U
en este caso). La duracin de los procesos que producen fraccionamiento en diversos
ambientes, en particular el fluvial, es del mismo orden de magnitud que las vidas medias de tres
nucleidos de la serie del uranio-238 (238U): ~245 ka para el uranio-234 (234U), ~75 ka para el
torio-230 (230Th) y ~1,6 ka para el radio-226 (226Ra). Gracias a ello, es posible medir
desequilibrios significativos en muestras de sedimentos fluviales. Luego, como explican
Chabaux et al. (2003), a partir de su adecuada interpretacin y modelamiento se puede estimar
(o al menos acotar) tiempos de maduracin y transporte sedimentarios (ej. Granet et al., 2007;
Dosseto et al., 2006a,b,c; Vigier et al., 2001; Plater et al., 1988, 1994). Otros estudios han
logrado caracterizar cuantitativamente la dinmica de la meteorizacin y de la formacin de
suelos (ej. Braun et al., 1993; Gueniot et al., 1988a,b).

Es importante tambin mencionar, que existen otras aproximaciones para abordar la


pregunta del tiempo de transferencia sedimentaria. Por ejemplo, los nucleidos cosmognicos
(TCN: Terrestrial Cosmogenic Nucleids) han sido medidos en sedimentos para cuantificar tasas
de erosin promedio a la escala de cuenca (ej. Granger et al.,1996; Schaller et al., 2001; Riebe
et al., 2001a,b; Schaller et al., 2002; Riebe et al., 2004; Schaller et al., 2004; von Blanckenburg
et al., 2004; von Blanckenburg, 2005; Kirchner et al., 2006, Carretier et al., 2009), como tambin

9
edades de exposicin de depsitos aluviales (ej. Ritz et al., 1995; Anderson et al., 1996; Repka
et al., 1997; Brown et al., 1998; Hancock et al., 1999; van-der Woerd et al., 2002; Siame et al.,
2004; Ritz et al., 2006; Regard et al., 2006). Sin embargo, no existe hasta la fecha un mtodo
establecido para la cuantificacin a escala geolgica de los tiempos de transporte de
sedimentos, siendo sta una lnea de investigacin relativamente nueva.

El enfoque de los desequilibrios en las series del U ha sido usado en estudios de transporte
de sedimentos en ros de larga trayectoria y con extensas hoyas hidrogrficas, por ejemplo, en
tributarios del ro Amazonas (Dosseto et al., 2006a,b,c) y en tributarios del ro Ganges (Granet
et al., 2007). Sin embargo, este enfoque no ha sido empleado an en ros de trayectoria corta
en un margen continental activo, por ejemplo en la vertiente occidental de los Andes.

A lo largo de la vertiente occidental de los Andes Centrales, donde se enmarca la regin de


inters de este trabajo (cuenca del ro uble), las principales caractersticas tectnicas relativas
a la subduccin de la placa de Nazca bajo la placa Sudamericana (espesor litosfrico, ngulo
de subduccin, ngulo, tasa de convergencia, historia de la deformacin, ...)1, as como el clima
(precipitaciones, temperaturas, historia climtica, ...) presentan variaciones desde importantes a
enormes. Por otro lado, la morfologa de un ro y su cuenca dependen a escala geolgica de
varios factores: del clima (Vandenberghe, 2002; Huisink, 1999; Bogaart et al., 2002; Bogaart y
van Balen, 2000), de la tectnica regional (Riebe et al., 2001a,b), de las litologas presentes y
de la actividad volcnica. Por ello, no es extrao encontrar a lo largo de Chile, cuencas con una
amplia gama de morfologas. Es intuitivo plantear entonces que la capacidad de transporte de
un ro y la tasa de erosin promedio en su cuenca, dependen a escala geolgica, de forma
directa, y en mayor o menor medida, de todos los factores climticos y tectnicos mencionados,
as como de las litologas y la actividad volcnica.

En esta Memoria de Ttulo se aborda la cuantificacin de tiempos de transporte de


sedimentos a lo largo del ro uble utilizando los desequilibrios de la serie del uranio-238 en
muestras de sedimentos. Este estudio constituye un primer paso dentro de un proyecto de
investigacin ms amplio, el proyecto SESPEED, que pretende probar dos mtodos para medir
velocidades de transporte de sedimentos gruesos fluviales (carga de fondo) en los Andes
Centrales entre los 30S y los 38S (ros Bo-Bo, Tigiririca, Cachapoal, Limar, Choapa, etc.):
desequilibrios en la serie del U-238 y nucleidos cosmognicos. Este proyecto es financiado por
el IRD de Francia (Institut de Recherche pour le Dveloppement) y dirigido por el Dr. Carretier
del LMTG de Toulouse (Laboratoire des Mchanismes et Transferts en Gologie). Participan del
proyecto investigadores del LMGT y del LHyGeS de Estrasburgo (Laboratoire d'Hydrologie et de
Gochimie de Strasbourg), con la colaboracin de dos acadmicos del Departamento de
Geologa (FCFM, Universidad de Chile).

Para esta memoria se escogi el ro uble, por dos motivos: (a) se contaba con suficientes
muestras para estudiar la evolucin de sus sedimentos a lo largo de su curso, y (b) al presetar
su cuenca una escorrenta superficial grande y un clima hmedo respecto de los ros
muestreados para el proyecto, se esperaba a priori encontrar en el uble desequilibrios
mximos dentro de este conjunto de muestras.

1
A lo largo de la vertiente occidental de los Andes Centrales el espesor elstico litosfrico aumenta de sur a norte
(Tassara et al., 2007). Tambin existe en esa zona un segmento de flat-slab (subduccin plana) entre los 30S y los
33S con una gran variacin del ngulo de subduccin, que coincide espacialmente con la zona de los valles
transversales. Respecto de la historia de la deformacin, adems del alzamiento de los Andes, un buen ejemplo es el
Sistema de Fallas de Atacama, que jug un rol fundamental durante el Negeno entre los 20S y los 27S, actuando
como una barrera para los sedimentos provenientes del este (Riquelme et al., 2003).

10
1.2 Objetivos

Objetivo general
Estimar tiempos de transporte fluvial en la cuenca del ro uble, con el propsito de
caracterizar cuantitativamente su dinmica sedimentaria.

Objetivos especficos
Se emple un enfoque similar al usado por Granet et al. (2007) y Dosseto et al. (2006c), es
decir, medir, interpretar y modelar los desequilibrios de la serie del 238U en sedimentos fluviales
del ro uble. Se formularon los siguientes objetivos especficos.

1. Estudiar la qumica de elementos mayores y traza, la mineraloga y la isotopa de Sr-Nd


en las muestras, con el propsito de caracterizar la evolucin de los sedimentos.

2. Obtener, por espectrometra de masas, las razones de actividad entre nucleidos de la


serie del 238U (238U-234U-230Th-226Ra), as como las proporciones isotpicas de Sr y Nd,
en muestras de sedimento tomadas en distintos puntos del cauce principal del ro uble,
y otros puntos de la cuenca del ro Itata.

3. Usar los antecedentes geoqumicos y geolgicos, para constreir el comportamiento de


los istopos de U, Th y Ra, y a partir de ello interpretar los desequilibrios medidos y
estimar tiempos de residencia sedimentaria.

11
CAPTULO 2: Series del uranio y del torio

2.1 Generalidades

Las denominadas series del uranio y del torio corresponden a tres cadenas de
desintegracin radioactiva, que comienzan con un nucleido de la serie de los actnidos con una
vida media (t) larga, especficamente el uranio-238 (238U, t 4,47 Ga), el uranio-235 (235U,
t0,707 Ga) y el torio-232 (232Th, t14,1 Ga), y terminan con un istopo estable de plomo
(206Pb, 207Pb y 204Pb, respectivamente; ver Figura 1). Entre el padre (actnido) y el hijo (plomo)
existen pasos intermedios de decaimiento, los cuales tienen vidas medias mucho menores,
desde microsegundos a pocos cientos de miles de aos. En un sistema cerrado, es decir, sin
intercambio de materia con su medio, las concentraciones Nk de los nucleidos de cualquier serie
de desintegracin radioactiva, satisfacen el siguiente sistema de ecuaciones diferenciales
(Bourdon et al., 2003).
dN 1
= 1 N 1
dt
dN k
= k N k + k 1 N k 1 k 2
dt
Las constantes k, llamadas constantes de desintegracin, definen las velocidades de
decaimiento radioactivo de cada reaccin en la cadena, siendo inversamente proporcionales a
las vidas medias2. Aplicando la transformada de Laplace a estas ecuaciones, Bourdon et al.
(2003) llegaron a la solucin explcita de las cantidades Nk como funciones del tiempo.
m
t
N m (t ) = Aim e
i
m {1,..., n} t 0
i=1
i m
i k
Aim N 0j
j=1 m k = j;k i k i

Naturalmente, se cumple la siguiente condicin de consistencia.

N m (0 ) = N m0 i.e.
m
N m0 = Aim m {1,..., n}
i=1

Equilibrio secular y desequilibrio:

Se define primero el concepto de equilibrio radioactivo, y enseguida el de desequilibrio. Por


simplicidad, consideremos los dos primeros nucleidos en la serie, es decir el nucleido padre A1 y
su hijo A2, cuyas concentraciones en un sistema cerrado satisfacen las ecuaciones de arriba,
i.e.

2
= ln(2)/ t

12
Figura 1: Vidas medias de los nucleidos de las series del uranio y del torio (Bourdon et al., 2003).

13
t
1
N 1 = N 10 e
1 2t 1 t
N 2 = N 20 N 10 e + N 10 e 1
2 1 2 1
De la Figura 1, es claro que 2 >>1 , luego T1 >> T2 . Consideremos el tiempo t que ha
transcurrido desde que el sistema estuvo abierto, es decir desde que realiz intercambios de
elementos con el medio por ltima vez. Si t es suficientemente grande para que exp(-2t) sea
despreciable, pero no as exp(-1t), se puede hacer la siguiente aproximacin.
1 N 10 2t N
2 N 2 = 2 N 2
0
e + 1 1
1 1 / 2 1 1 / 2
t
( )
... 2 N 20 1 N 10 e 2 + 1 N 1 1 N 1
Si se considera la ecuacin diferencial para la cantidad N2, se notar que la condicin
anterior implica directamente dN2/dt 0. Esta condicin se denomina equilibrio radioactivo o
equilibrio secular (Bourdon et al., 2003), y consiste para un sistema cerrado, en un equilibrio
entre la tasa de formacin de nucleidos A2 a partir de A1 (N11), y la tasa de formacin de
nucleidos A3 a partir de A2 (N22). En este contexto, la actividad de un nucleido A1 se define
como N11, y siendo A1 el padre de A2, (A2/A1) denota la razn de actividad ("activity ratio") del
sistema A2-A1, dada por (A2/A1) (N22)/(N11) (Bourdon et al., 2003). La condicin de equilibrio
radioactivo se reduce a la igualdad entre las actividades del nucleido padre y la del hijo. El
equilibrio secular equivale a una razn de actividad (A2/A1) igual a 1.

Un desequilibrio en la cadena o serie de desintegracin se define como cualquier


desviacin de la condicin de equilibrio secular. Tomando en cuenta el caso del sistema A1-A2,
la siguiente aproximacin es vlida si 1<< 2.
t
(
2 N 2 2 N 20 1 N 10 e ) 2
+ 1 N 1
A A2 t
luego 2 1 1 e 2
A1 A1 0

Es directo constatar que el proceso de retorno al equilibrio radioactivo depende


esencialmente de 2, sean cuales sean las condiciones iniciales de desequilibrio (2N02 - 1N01
0). Y el tiempo de retorno al equilibrio ser TEq 6(t)2 = 6ln(2)/2 despus del cual, si el
sistema ha permanecido cerrado, se habr restablecido el equilibrio radioactivo. As, TEq es de
unas seis veces la mxima vida media de los nucleidos intermedios de la cadena (todos menos
el primero, el padre). Visto de otro modo, la velocidad de reaccin (regreso al equilibrio) est
siempre determinada por la constante de velocidad ms lenta. Por ende, para la serie del
uranio-238, despus de seis veces la vida media del uranio-234 (aprox.1,5 Ma) un sistema
cerrado estar en equilibrio secular. La clave de este mtodo est en que las diferencias de
comportamiento geoqumico de los elementos involucrados: U, Th, Ra, etc. (ej. comportamiento
en solucin, adsorcin, etc), producen fraccionamiento y ruptura del equilibrio.

14
2.2 Causas de desequilibrio en la serie del U-238

Como se dijo en la seccin 1.1, la serie del 238U es de particular inters para su aplicacin
al estudio de fenmenos que ocurren en escalas de tiempo del orden de decenas a cientos de
miles de aos, hasta 1 Ma (similar a las vidas medias del 234U, 230Th y 226Ra). Tal es el caso de
la maduracin y el transporte de sedimentos fluviales, y es por eso que se concentra la atencin
en esta serie. Los desequilibrios en la serie del 238U son producidos mediante dos categoras de
mecanismos: el fraccionamiento qumico y los efectos del alpha-recoil.

2.2.1 Fraccionamiento qumico

Algunos procesos naturales, tales como los cambios de fase, la meteorizacin o la fusin
parcial de una roca, producen difusin qumica entre el sistema (rocas, granos minerales,
magma, sedimentos, agua, etc.) y su medio, de forma tal que el intercambio de elementos es
diferenciado. Es decir, hay elementos que entran o salen del sistema en forma preferencial
respecto de otros, modificando las proporciones entre los elementos que afectan (sus
abundancias relativas). Este fenmeno se llama fraccionamiento qumico, y depende de las
caractersticas geoqumicas de los elementos en cuestin (en el caso de la serie del 238U stos
son U, Th, Ra, Pa, ...), as como de las condiciones qumicas y fsicas del sistema y su entorno.
Si se trata de aquellos elementos presentes en una serie de desintegracin radioactiva, es
evidente que el fraccionamiento provoca desequilibrio radioactivo3. Para los ambientes
superficiales, en particular el fluvial, las variables que controlan el fraccionamiento son: la
temperatura, la composicin del agua, la actividad de microorganismos, la presencia de
coloides, presencia de minerales capaces de adsorber elementos, etc. Los mecanismos de
fraccionamiento qumico en ambientes fluviales son detallados ms adelante (ver 2.3).

2.2.2 Alpha-recoil

Este fenmeno ocurre con aquellos nucleidos radioactivos con modo de decaimiento alfa,
es decir, cuyo mecanismo implica la desintegracin de su ncleo en uno ms liviano (en 4
u.m.a.) y un ncleo de He o partcula alfa (e.g. 238U, 234U o 230Th en la serie del 238U). Cuando
ocurre la desintegracin radioactiva del nucleido padre, debido a la fuerza de Coulomb (y a las
fuerzas nucleares) entre el nucleido hijo y la partcula alfa, stos se alejan uno del otro
violentamente. La partcula alfa es mucho ms liviana que el nucleido hijo, y es eyectada a una
gran velocidad, mucho mayor que la velocidad de retroceso (recoil) del hijo, por conservacin
de momentum y energa. El ncleo de He causa una estela recta de daos en la estructura
cristalina, lo que puede facilitar localmente la difusin de elementos dentro del cristal, en
particular del nucleido hijo y/o padre. Por otro lado, el nucleido hijo puede ser dislocado de su
posicin cristalina, propiciando su movilizacin. Si el padre estaba cerca del borde del cristal,
esa dislocacin puede significar la salida directa del hijo del cristal hacia una fase vecina (ej.

3
Tambin existe el fraccionamiento isotpico, i.e. el intercambio entre el sistema y su medio de un istopo con
respecto a otro en forma diferenciada. La naturaleza diferenciada de dicho intercambio, se debe a las diferencias
relativas de masa entre los istopos de uno o ms elementos. Por ejemplo, el oxgeno tiene tres istopos estables
16 17 18
( O, O y O), siendo el ms liviano tambin el ms abundante. Cuando el agua de mar se evapora, aquellas
16
molculas de H2O que contienen slo O cambian de fase con mayor probabilidad que el resto de las molculas
de agua. Sin embargo, es sabido que este fenmeno afecta principalmente a los elementos de baja masa
atmica, y por ende bajo nmero atmico. Los istopos de elementos ms pesados, como el uranio o el torio,
tienen diferencias relativas de masa tan pequeas, que el efecto es en general despreciable.

15
otro cristal, una solucin acuosa o la atmsfera).

En la serie del 238U, este mecansmo es el responsable del fraccionamiento isotpico entre
238
el U y el 234U, pues propicia la movilidad del nucleido hijo respecto de la del padre.

2.3 Desequilibrios en ambientes acuosos

2.3.1 Movilidad del Th

Chabaux et al. (2003) afirman que el Th es altamente insoluble y casi inmvil en un perfil de
meteorizacin con pH>3-4, aunque puede llegar a ser mvil en presencia de cidos orgnicos.
Segn los mismos autores existen dos explicaciones a este comportamiento: (a) el Th se
encuentra principalmente en minerales residuales primarios y secundarios, incluyendo
minerales accesorios relictos, y posiblemente en minerales de arcilla (Middleburgh et al., 1988;
Braun et al., 1993); y (b) el Th se encuentra fuertemente ligado o adsorbido a minerales
neoformados como partculas de arcilla y/o hidrxidos de Fe (Middleburgh et al., 1988; Gueniot
et al., 1988a,b). As, por un lado el Th se da en minerales resistentes a la meteorizacin y la
erosin, y por otro, si llega a ser disuelto, es muy probable que no migre significativamente al
ser adsorbido en forma muy estable a ciertos minerales. Sin embargo, por su alta afinidad a la
materia orgnica, puede ser movilizado con sta, as como mediante algunos complejos
inorgnicos, como sulfato, fluoruro, fosfato e hidrxido. El Th slo se encuentra en la naturaleza
como un catin tetravalente (Chabaux et al., 2003).

2.3.2 Movilidad del U

A diferencia del Th, el U es altamente mvil. Los datos termodinmicos indican que en
condiciones superficiales, las especies de U estables en presencia de H2O ms abundantes son
los estados de oxidacin UIV y UVI (Chabaux et al., 2003). En condiciones oxidantes, el U ocurre
como UVI, formando el ion uranilo (uranyl) UO22+, y forma fcilmente complejos con carbonato e
hidrxido, y tambin con fosfato y fluoruro (Chabaux et al., 2003). La formacin de estos
compuestos de uranilo aumenta sustancialmente la solubilidad del U, por ejemplo, la solubilidad
de la uraninita aumenta debido a la presencia de CO2 y la consecuente formacin de
carbonatos de uranilo (Chabaux et al., 2003). Ligandos orgnicos, tales como los cidos
hmico, flvico y ctrico, tambin pueden unirse con fuerza al U y modificar sustancialmente su
movilidad (Gascoyne, 1992; Lenhart et al., 2000; Montavon et al., 2000).

El fenmeno de adsorcin o retencin en la superficie de minerales es de gran relevancia


en el comportamiento del U en ambientes acuosos, pudiendo limitar su movilidad. La capacidad
de adsorcin de U disminuye desde los hidrxidos de Fe y el gel de slice, a las arcillas y micas,
y a los palos (Chabaux et al., 2003). Por ejemplo, en ambientes reductores, la bentonita es un
buen fijador de uranio, y parece ser un catalizador de la reduccin del UVI adsorbido en UIV
insoluble por materia orgnica (Giaquinta et al., 1997). El pH de la solucin, entre otras varias
variables fsico-qumicas, influencia fuertemente la capacidad de retencin de U por los
minerales (Chabaux et al., 2003). Los mecanismos involucrados en la inmovilizacin parcial del

16
U en suelos y perfiles de meteorizacin son: (1) su concentracin en minerales accesorios
relictos, y (2) adsorcin a o coprecipitacin con hidrxidos de Fe amorfos, minerales de arcilla
y/o materia orgnica (Moreira-Nordemann y Sieffermann, 1979; Tieh et al., 1980; Gueniot et al.,
1988a,b; Middleburgh et al., 1988).

La actividad de microorganismos es otro fenmeno que puede afectar la movilidad del U en


la naturaleza. Se ha mostrado que bacterias reductoras de hierro y azufre son capaces de
reducir UVI en UIV (Fredrikson et al., 2000; Lovley et al., 1991; Lovley y Phillips 1992; Spear et
al., 2000). Segn Suzuki y Banfield (1999), los numerosos mecanismos de acumulacin de U
por microbios, se subdividen en dos grandes categoras: dependientes del metabolismo (va la
accin de enzimas) e independientes del metabolismo (ej. interacciones de algn complejo de U
con zonas de la clula cargadas elctricamente). La fijacin de radio-nucleidos por coloides
(macromolculas, virus y fragmentos de microbios: de 1 a 1000 nm de dimetro) produce de
forma variable, y a veces importante, fraccionamiento en el sistema U-Th (Chabaux et al.,
2003).

2.3.3 Movilidad del Ra

La solubilidad del Ra es generalmente baja excepto en aguas salinas y reducidas (no es el


caso de ros) donde puede formar complejos con cloruros y sulfatos (Dickson, 1985; Langmuir y
Melchior, 1985; Sturchio et al., 1993) o con carbonatos en solucin acuosa (Rihs y Condomines,
2002). El Ra tambin puede ser adsorbido bajo determinadas condiciones, lo que puede limitar
su movilidad (Chabaux et al., 2008).

2.3.4 Alpha-recoil

Finalmente, pero no de forma menos relevante, el fenmeno del alpha-recoil puede afectar
la movilidad del tomo hijo (Chabaux et al., 2003).

Los daos en la estructura cristalina (y/o la oxidacin de 238UIV en 234UVI) debido a las
reacciones de desintegracin radioactiva, aumentan la movilidad del 234U respecto al 238U,
hacindolo ms propenso a ser incorporado a una solucin acuosa (Chabaux et al., 2003). Los
resultados de varios estudios, tanto en silicatos como en minerales radioactivos, indican que el
nucleido hijo puede ser movilizado ms fcilmente debido al "recoil" y que la eficiencia de su
incorporacin en solucin es funcin del tiempo de "annealing" o tiempo de recuperacin de la
estructura cristalina. Este mecanismo de fraccionamiento es ms efectivo, mientras ms corta
sea la escala temporal sobre la cual ocurre la disolucin del nucleido hijo (a su vez proporcional
a la vida media del padre) en comparacin con el tiempo de "annealing"4. En consecuencia,
este mecanismo es ms eficiente para nucleidos padres con vidas medias relativamente cortas
(Chabaux et al., 2008).

La eyeccin directa de tomos desde fases minerales hacia una fase fluida es uno de los
mecanismos de liberacin de 234U en aguas naturales (Kigoshi, 1971) y de Th, Ra y Rn en
aguas subterrneas (Sun y Semkow, 1998), y da cuenta del desequilibrio 238U-234U (Chabaux et
al., 2003). El 226Ra puede ser directamente eyectado a la fase acuosa por "alpha-recoil" de

4
Para la muscovita el tiempo de "annealing" puede llegar a los 100 ka, mientras que es de 15 ka para la monacita
(fosfato rico en tierras raras, entre ellas Th).

17
nucleos de 230Th. Para una distribucin homognea de 230Th en el volumen de un grano mineral,
se predice una proporcin entre tomos de 226Ra eyectados y producidos entre 1,310-4 y 310-2
en funcin del tamao y la forma del grano (Chabaux et al., 2008). Sin embargo, se han
observado elevadas concentraciones de 222Rn (un gas, hijo directo del 226Ra mediante
decaimiento alfa) en ciertos materiales naturales, lo que sugiere que en estos casos el nucleido
padre se concentra en los bordes de los granos (Krishnaswami y Seidemann, 1988; Swarzenski
et al., 2003). En este contexto, la adsorcin de nucleidos de la serie del 238U a granos minerales
adquiere gran relevancia.

18
CAPTULO 3: Hidrologa, morfologa y clima

El ro uble, de donde vienen casi todas las muestras, es el principal tributario del ro Itata,
que fluye de este a oeste en la parte norte de la VIII Regin del Bo-Bo, en Chile Central. La
cuenca hidrogrfica del ro Itata se sita entre las latitudes 3610'S y 3720'S y las longitudes
7100'W y 7255'W aproximadamente, y abarca una superficie de 11.294 km2 (ver Figura 2).
A continuacin se describen las principales caractersticas hidrolgicas, morfolgicas y
climticas de la cuenca de este ro.

3.1 Hidrologa

El ro principal se origina en la alta cordillera, la confluencia de los ros Cholgun, que viene
desde el norte, y Hupil, desde el sur, a pocos kilmetros al norte del ro Laja. En un recorrido
de 85 km con rumbo al NNW hasta su junta con el uble, capta la aguas de sus principales
tributarios que constituyen una red de drenaje paralela de direccin al poniente y cuyas
cabeceras se encuentran en la zona de La Montaa. De ella forman parte los ros Danicalqui,
Diguilln y Larqui.

Poco ms arriba de su junta con el Danilcalqui, el Itata presenta un salto de 20-25 m de


altura. A partir de l, ro fluye en la Depresin Intermedia entre riberas bajas y cultivables. En el
curso inferior, despus de trasponer la Cordillera de la Costa, los bancos de arena obligan a la
corriente a extenderse considerablemente en vegas y su profundidad disminuye. El ancho en la
boca del Itata vara de 100 a 200 m en verano, y llega a 300 m en invierno.

El principal afluente del Itata es el uble, cuya cuenca cubre 5.097 km2. Nace al pie del
paso de Buraleo, al oriente de los nevados de Chilln (3212 m s.n.m). Desarrolla su curso
superior en direccin al NNW y a 40 km de su origen recibe desde el norte su principal afluente,
el ro Los Sauces. Tras un recorrido de 155 km, el uble se junta en el borde oriental de la
Cordillera de la Costa al Itata en una localidad llamada Confluencia.

Despus de la confluencia con el uble, el Itata desarrolla su curso inferior en el corazn


de la Cordillera de la Costa, donde slo recibe esteros menores, siendo el ro Lonqun el de
mayor caudal. Drena un rea de cierta importancia (1.075 km2) de los cerros costeros al sur y al
este de la ciudad de Quirihue.

3.2 Relieve

Se reconocen los siguientes rasgos geomorfolgicos relevantes, de orientacin NNE-SSW


(Gajardo, 1981; Muoz y Niemeyer, 1984).

o Llanuras litorales: constituidas por sedimentos marinos terciarios y fluvio-marinos del


cuaternario, se elevan hasta alturas de 300 m s.n.m. y presentan dunas, cordones
litorales, marismas y paleolechos.
o Cordillera de la Costa: muestra los rasgos de un macizo cordillerano en estado de
madurez, atravesado por ros que desembocan en el mar. De relieve aplanado, contiene
algunos cerros que dominan el paisaje con alturas de 550 m a 770 m s.n.m.

19
Figura 2: Relieve y red de drenaje del ro Itata. Los puntos muestreados se indican con crculos rojo.
La red de drenaje y los lmites de cuenca fueron tomados del mapa digital del IGM. El relieve se traz a
partir de un archivo "raster" de la topografa de Sudamrica de 90 m de resolucin.

20
o Depresin Central: amplia llanura de unos 150 m de altitud media y de superficie
ondulada hacia el sur. Cerca de Chilln alcanza su mximo ancho este-oeste con 40 km.
o Precordillera: es la unidad morfolgica transicional entre la alta montaa por el este y la
Depresin Central por el oeste. Su altura vara entre los 400 m y 500 m s.n.m., y est
profundamente disectada por los causes que la atraviesan. Al norte del ro Cato
(~364'S), el quiebre de pendiente entre esta unidad y la Depresin Central es
fcilmente identificable, mientras que al sur de este ro el relieve pasa gradualmente de
la montaa a la Depresin Central.
o Cordillera Principal: sus contrafuertes ms occidentales, cuyas cumbres ms altas
alcanzan los 1800 a 2000 m s.n.m., constituyen un relieve abrupto. Ms al este el relieve
es ms suave y la mayora de las cumbres no sobrepasan los 2500 m s.n.m, alcanzando
excepcionalmente los 3000 m s.n.m.

3.3 Clima

El clima de la cuenca del ro Itata es mediterrneo y con al menos dos meses consecutivos
con dficit hdrico en verano. El patrn de variacin ombroclimtica es el caracterstico para la
zona centro sur de Chile: las precipitaciones tienden a ser mayores al oeste de ambas
cordilleras, pues ellas constituyen barreras para los frentes hmedos que llegan al continente
desde el Ocano Pacfico, diferencindose sectores de ombroclima hmedo y subhmedo, con
montos de precipitacin que varan en torno a los 1.000 mm anuales (Fuenzalida, 1971).

Sobre la base del rgimen pluviomtrico y trmico, se distinguen dos tipos climticos:
Mediterrneo pluviestacional - ocenico, en casi toda el rea de la cuenca, en que la amplitud
trmica anual es inferior a 20 C.
Mediterrneo pluviestacional continental, en los pocos sectores en que la altitud se eleva por
sobre los 2.000 m s.n.m, donde - por efecto de la altitud, la amplitud trmica anual suele ser
superior a 20 C.
La Estacin Meteorolgica Chilln registra una temperatura promedio anual de 14,08 C,
con una mnima de 7,6 C y una mxima de 20,56 C (promedios mensuales). La precipitacin
promedio es de 85,43 mm y el total anual de agua cada es de 1.025 mm (ver Figura 3).

Figura 3: Precipitacin y temperatura medias mensuales en la Estacin Meteorolgica de Chilln


(Informe DGA, 2004).

21
CAPTULO 4: Marco geolgico

En la zona geogrfica de la cuenca del ro Itata afloran rocas metamrficas asignadas al


Paleozoico, rocas sedimentarias y volcnicas, marinas y/o continentales del Mesozoico-
Cenozoico y sedimentos no consolidados del Cenozoico (Gajardo, 1981).

4.1 Paleozoico

Serie Oriental (Basamento Metamrfico)


Esquistos, filitas, gneises y metarenitas producto de metamorfismo de baja presin-alta
temperatura. Corresponden a las facies de esquistos verdes y granulita, de oeste a este
aparecen sucesivamente la biotita, la andalucita y la sillimanita (Aguirre et al., 1972). La base de
esta unidad es desconocida. La sobreyacen parcialmente y en discordancia angular y de
erosin unidades sedimentarias del Cretcico y Cenozoico. Edades radiomtricas del
metamorfismo sugieren una edad mnima paleozoica (Gajardo, 1981).

4.2 Mesozoico

Estratos de Pocillas - Coronel del Maule - Quirihue


Trisico

De arriba hacia abajo esta secuencia est constituda por: alternancia de andesitas
porfricas, areniscas cuarcferas, brechas volcnicas y riolitas, chert y alternancia de areniscas
cuarcferas y lutitas con intercalaciones de conglomerados de cuarzo (Lagno, 1979). De
acuerdo a la flora fsil encontrada por Moreno-Roa et al. (1976), se le asigna una edad rtica a
las rocas sedimentarias de los Estratos de Pocillas.

Formacin Rio Blanco


Cretcico

Lavas andesiticas e intercalaciones de brechas piroclsticas y epiclsticas. Gardeweg


(1980) le asigna una edad cretcica tarda en base a una datacin K-Ar de granitoides que la
intruyen (60,5 0,3 Ma). La base de esta unidad no se encuentra en la zona y est cubierta
discordantemente por la Formacin Cola de Zorro (Gajardo, 1981).

4.3 Cenozoico

Formacin Cura-Malln
Eoceno-Mioceno Inferior

Muoz y Niemeyer (1984) definen dos miembros: Ro Queuco (volcanoclstico) y Malla-Malla


(epiclstico). El Miembro Ro Queuco contiene casi la totalidad de los afloramientos de la
Formacin Cura-Malln en el rea de estudio, y est compuesto por tobas (de cristales, y en

22
Figura 4: Unidades geolgicas de la cuenca del ro Itata. Los puntos rojos son los puntos de muestreo
(Gajardo, 1981; Ferraris, 1981; Gajardo y Niemeyer, 1984).

23
menor cantidad vtreas y de lapilli), brechas (lticas, con fragmentos de andesitas y dacitas
porfricas y pmez) y, de forma subordinada, areniscas, conglomerados, lutitas y coladas de
lava. El Miembro Malla-Malla est formado por areniscas, conglomerados y lutitas, y en menor
proporcin por tobas, brechas, calizas y niveles carbonosos. La formacin es cubierta
discordantemente por la Formacin Cola de Zorro (Gajardo, 1981).

Formacin Cola de Zorro


Plioceno Superior-Pleistoceno

Andesitas, andesitas baslticas, de naturaleza calcoalcalina, y de forma ms escasa


basaltos, dacitas y rocas piroclsticas (tobas lticas y de cristales, y brechas).

En las andesitas el mineral mfico ms abuntante es el clinopiroxeno, con o sin presencia


de ortopiroxeno y/u olivino. La textura ms comn es porfrica con masa fundamental intersertal
a pilotaxtica, siendo en ocasiones traqutica. Los minerales de alteracin son iddingsita en
ortopiroxenos y arcilllas en plagioclasas. En las andesitas baslticas, el olivino es menos
abundante que el clino- y el ortopiroxeno. La textura es porfdica, por lo general seriada con
masa fundamental de diversos tipos (intergranular, intersertal, hialopiltica o pilotaxtica). Los
minerales de alteracin corresponden a iddingsita, bowlingita y montmorillonita. En los basaltos,
los fenocristales de plagioclasa predominan sobre el olivino y los piroxenos, que estn
esencialmente en la masa fundamental.

Formacin Trapa-Trapa
Mioceno Medio-Superior

Incluye andesitas, andesitas baslticas y dacitas calcoalcalinas, adems de tobas y


brechas, localmente arenosas o conglomerdicas. Las andesitas tienen textura
predominantemente porfrica, con vesculas y amgdalas. Los fenocristales de plagioclasa
(andesina-labradorita) son abundantes y alcanzan tamaos de hasta 1 cm, siendo los de
piroxeno ms pequeos y escasos. La masa fundamental es normalmente hialopiltica,
apareciendo tambin felstica intersertal, intergranular y pilotaxtica, donde se observan
microlitos de plagioclasa, clinopiroxeno, minerales opacos y cantidades variables de vdrio
(Muoz y Niemeyer, 1984). En algunas muestras, Muoz y Niemeyer (1984) observaron
anfbola (hornblenda), ortopiroxeno y olivino, en cantidades subordinadas.

Formacin Mininco
Plioceno

Sedimentitas lagunares y fluviales, representadas por areniscas tobceas, limolitas y


conglomerados, con intercalaciones de arcillolitas y tobas. En las areniscas, el componente
volcnico es principalente andestico, encontrndose tambin ceniza en la matriz. (Gajardo,
1981).

Formacin La Montaa
Pleistoceno

Corresponde principalmente a depsitos morrnicos compuestos por sedimentos no-


consolidados de escasa compactacin, arcillosos, de color rojo a pardo-rojizo, que contienen de
forma catica fragmentos de rocas volcnicas, baslticas y granticas, de distinto grado de
redondeamiento y en diferente estado de alteracin. El tamao de los clastos vara entre 7 y 40
cm de dimetro promedio, habiendo en algunos sectores bloques de hasta 5 m. Esta unidad

24
presenta intercalaciones de arcillas, arenas y cenizas volcnicas bien estratificadas que se
interpretan como glaciolacustres. Se dispone en discordancia de erosin, sobre rocas
estratificadas e intrusivas, mesozoicas y cenozoicas, y sobre los sedimentos de la Formacin
Mininco (Gajardo, 1981).

Sedimentos de terrazas fluviales


Holoceno

Depsitos sedimentarios distribuidos en los valles principales, formando amplias superficies


aterrazadas, y con espesores mximos de 20 m. Estn compuestos por gravas y ripios, con
clastos redondeados a subredondeados de andesitas y granitos de dimetro promedio entre 8 y
30 cm, y proporciones variables de arena y limo.

Sedimentos de valles actuales


Holoceno

Materiales no consolidados, incluyen gravas, ripios y arenas y forman el relleno actual de


los cauces fluviales. Estos depsitos alcanzan su mxima extensin areal en la Depresin
Central y la Cordillera de la Costa, donde los ros pierden parte de su capacidad de transporte.

Sedimentos del cono del rio Laja


Holoceno

Depsitos de arenas baslticas y andesticas, compuestos por materiales piroclsticos


provenientes del volcn Antuco y centros adyacentes, cuyo transporte fue fluvial y elico
(Gajardo, 1981; Ferraris, 1981).

4.4 Rocas intrusivas

Batolito de la Costa (Pzg)


Paleozoico Superior-Trisico Inferior

Este complejo plutnico est compuesto por rocas de composicin petrogrfica variable
entre tonalita y granodiorita, con sectores menos extensos de diorita cuarcfera y granito
(Gonzlez-Bonorino, 1970; Herv et al., 1976). Estas rocas tienen grandes cristales de
microclina, y han estado sometidas a intensos fenmenos cataclsticos, evidencia de lo cual es
el fracturamiento de los cristales de cuarzo y feldespato (Escobar et al., 1977). El Batolito de la
Costa intruye a las rocas del Basamento Metamrfico e infrayace con discordancia de erosin a
las rocas sedimentarias y sedimentario-volcnicas del Trisico, lo que permite asignarle una
edad paleozoica tarda a trisica temprana (Gajardo, 1981).

Intrusivos de la Cordillera de la Costa (Jg)


Jursico

Cuerpos intrusivos de composicin similar al Batolito de la Costa, pero de textura ms fina.


Dataciones radiomtricas por el mtodo K-Ar indican edades entre 207 Ma y 172 Ma (Herv y
Munizaga, 1978).

25
Intrusivos de la Cordillera Andina (Kg)
Cretcico

Cuerpos intrusivos granticos, tonalitas y granodioritas. Dataciones radiomtricas por el


mtodo K-Ar indican edades mnimas de 83,9 3,8 Ma y 85,4 5,2 Ma (Gajardo, 1981).

26
CAPTULO 5: Mtodos analticos

5.1 Molienda de muestras

Se realizaron dos etapas de molienda. La primera permiti obtener para cada muestra
aprox. 20 a 30 g de polvo < 100 m, parte del cual (5 a 7 g aprox.) se volvi a moler hasta
obtener un polvo < 50 m. Para corroborar la calidad de la molienda, se hizo pasar con xito
uno de los frascos de polvo < 50 m por un tamiz de 20 m.

Primera etapa: se us una mquina agitadora, la que funciona moviendo en crculos un


recipiente de gata con una cavidad cilndrica y un disco del mismo material. La muestra se
introduce en pequeas cantidades cada vez (5 g) en el recipiente a un costado del disco, se
cierra el recipiente y se lo fija dentro de la maquina que se hace funcionar durante 5 min. a 700
rpm. Las partculas de la muestra son golpeadas y aplastadas entre el disco y las paredes del
recipiente. Al final de cada ciclo de molienda, se abre el recipiente y se hace pasar el polvo por
un tamiz plstico de 100 m con la ayuda de un pincel. Cada vez, al lograr hacer pasar toda la
muestra a travs del tamiz, se guard el polvo resultante en frascos plsticos nuevos de
polietileno previamente limpiados con HCl y agua qumicamente muy pura o "agua MilliQ"5
(MQ). Entre la molienda de cada muestra se realiz un minucioso protocolo de limpieza del
material. Primero se moli 2,5 g de cuarzo puro durante 5 min, para eliminar las partculas
incrustadas en el gata, descartando el polvo resultante. Enseguida se limpi el material (pincel,
tamiz y recipiente y disco de gata) con una escobilla, una esponja y detergente, enjuagando
con abundante agua de la llave y luego con agua MQ. Finalmente, para secar ms rpidamente
el material, se us etanol y un secador de pelo.

Segunda etapa: aqu se utiliz un recipiente muy similar y algo ms pequeo, pero esta vez
usando esferas de gata (de 5 y 20 mm de dimetro) en vez del disco. Se moli de 5 a 7 g del
polvo obtenido en la primera etapa (<100 m) y se hizo pasar el polvo resultante por un tamiz
plstico de 50 m. Se us una mquina agitadora, moliendo cada vez 1 a 2 g de muestra a 700
rpm, siendo todo el resto del protocolo idntico al de la etapa anterior.

5.2 Mineraloga por difraccin de rayos X

Se enviaron las muestras molidas (a < 100 m, ver 6.1) al laboratorio de difraccin de rayos
X de la Universidad de Estrasburgo, para determinar las especies minerales presentes en las
muestras.

Primero se realizaron anlisis de roca total, con la ayuda de un difractmetro de RX


Brker (modelo D5000), en modo / 2 (barrido de 3 a 65, por pasos de 0,02 de una duracin
de 2 s, antictodo de Cu, 40 kV 30 mA).

Para la determinacin de arcillas se agreg agua desmineralizada a cada muestra, y

5
El agua MilliQ (MQ) se prepara a partir de agua bidestilada, hacindola pasar por un sistema de ultrafiltrado
microscpico que elimina casi la totalidad de los iones en solucin. Al final, esta agua tiene una resistividad
elctrica de entre 18 y 20 M/cm. Esta agua fue usada en todas las manipulaciones qumicas (ver 6.3).

27
mediante centrifugacin, se recuper y decant la fraccin arcillosa (< 2 m). Eliminando el
agua, qued una pasta hmeda que fue esparcida con una esptula sobre una lmina de vidrio,
con el propsito de orientar los granos de arcilla paralelos a la lmina, de modo de reforzar la
reflexin de la familia de planos {0 0 1}. Se dej secar las lminas al aire. Enseguida se
realizaron 4 caracterizaciones (4 difractogramas) de cada muestra:
o sin tratamiento
o despus de tratamiento con etil-glicol
o despus de tratamiento con hidracina
o despus de calentamiento a 490C durante 4 horas.

El propsito de realizar mediciones con y sin estos tratamientos, es de diferenciar ciertos


minerales arcillosos que no son disociables con slo un difractograma (por ejemplo la refleccin
(001) de la kaolinita a ~7,15 coincide con la refleccin (002) de la clorita).

Las muestras fueron caracterizadas con la ayuda de un un difractmetro de RX Brker (modelo


D5000), en modo / 2 (barrido de 3 a 15 o 3 a 30, por pasos de 0,02 de una duracin de 1
s, antictodo de Cu, 40 kV 30 mA).

5.3 Manipulaciones qumicas

Enseguida se describe las principales etapas de preparacin de muestras y separacin


qumica de los elementos de inters.

5.3.1.- Ataques cidos


Se introdujo y pes aprox. 200 mg de polvo de muestra tamizado a < 50 m, dentro de
frascos de tefln especial para isotopa marca Savillex, luego se introdujo y luego pes algunas
gotas de "spike"6 (6 gotas del spike de U y Th, y 3 del de Ra; ver Tabla 1).

Ataques U-Th Ataques Ra


Muestra muestra spike U-Th muestra spike Ra
[mg] [mg] [mg] [mg]
CH-08-01 217,17 187,92 191,76 88,32
CH-08-02 181,65 182,93 198,45 64,00
CH-08-03 192,98 174,90 224,47 64,27
CH-08-04 195,16 196,30 200,54 53,68
CH-08-05 174,14 190,36 213,82 62,70
CH-08-06 204,48 179,29 205,80 60,62
CH-08-08 197,55 190,05 227,33 63,93
CH-08-09 205,83 187,00 293,92 88,97
CH-08-09 dup 225,27 214,20 202,83 64,80
CH-08-Bl 0,00 62,26 0,00 60,66

Tabla 1: Pesaje de muestras y spike, utilizando una balanza analtica con una precisin de 5 g.

6
Un spike es una solucin artificial que contiene istopos naturales y artificiales en concentraciones conocidas, y
proporciones muy distintas a las muestras naturales. Se usa en la tcnica de dilucin isotpica que permite medir
concentraciones de istopos (ver Anexo B).

28
Da Pasos Propsito
Introduccin y pesaje de ~200 mg de muestra (polvo < 50 m) en
frascos de tefln. Hacer mezcla
para dilucin
Introduccin y pesaje de spike (6 gotas para U-Th y 3 gotas para Ra)
isotpica.
1 +1 ml de HNO3 6N, procurando humectar bien el polvo.
+ 4 ml de HF; se cerr los frascos y se los someti a ultrasonido Romper los
durante al menos 10 minutos (para separar las partculas). enlaces Si-O, es
decir, disolver los
A 70 C durante toda la noche y la maana del da siguiente. silicatos.
2 +1 ml de HClO4
Favorecer las
Evaporar a 70 C durante un da y medio. reacciones de
3 oxidacin para
Cerrar y calentar a 150 C durante una noche
eliminar la
4 materia orgnica
Dejar a 120 C durante 2 das.

29
Verificar que no quede ninguna partcula negra (restos de materia
5 orgnica); de haber alguna, agregar 0,5 ml de HNO3 bidestilado, 0,5 Disolver los gels
ml de HCl bidestilado y 1 ml de HF, cerrar y calentar a 120 C, luego a de slice, fluoruros
150 C hasta que desaparezcan las partculas negras. de Fe, minerales
Agregar 1 ml de H3BO3, 2 ml de HCl 6N y 3 ml de H20, cerrar y pesa-dos y xidos

Tabla 2: Protocolo de ataques cidos


ellos tenan un aspecto difano y no contenan partculas visibles.

6 mantener a 100 C durante medio da. Evaporar a 110 C en el da y de Al.


150 C en la noche.
Verificar que no reste ninguna partcula con la linterna, si queda
7 algo, repetir los pasos del da 5 y 6.
Disolver la costra en 1 ml de HNO3 7,5N, cerrar y someter a Eliminar restos de
ultrasonido durante 10 min. Evaporar a 120 C, luego 150 C. Disolver cidos de la
en 1 ml de HNO3 7,5N, cerrar y someter a ultrasonido durante 10 min., solucin (salvo
dejar a 70 C durante una noche. cido ntrico).
Enseguida se realiz los ataques cidos; siguiendo el protocolo mostrado en la Tabla 2.

8
Trasvasijar a frascos ms pequeos. Centrifugar a 4000 rpm durante Decantar los cris-
10 min. tales precipitados.
original. Algunos minerales son especialmente resistentes a los cidos, por ejemplo el circn,
Es muy importante haber logrado un molido homogneo y suficientemente fino para obtener

que adems tiene un alto contenido de uranio y torio, de ah la importancia del molido. Los
ataques cidos aqu realizados fueron exitosos, puesto que las soluciones obtenidas mediante
una digestin completa de los minerales: al final no debe restar ninguna partcula del polvo
Tanto los ataques cidos como las separaciones qumicas (detalladas a continuacin), se
realizaron con guantes y delantal bajo una campana en una sala limpia. Una sala limpia es
una pieza donde el aire contiene muy pocas partculas en suspensin, gracias a varios filtros de
aire, a dos antesalas con sus respectivas puertas y a un riguroso protocolo de limpieza y
mantencin. Las campanas contienen an menos partculas en suspensin.

5.3.2.- Separacin y purificacin qumica

Se utilizaron dos tcnicas de separacin y purificacin de elementos qumicos (en este caso
U, Th, Ra, Sr y Nd). Para el U y el Th se us cromatografa de intercambio inico, basada en las
propiedades fsico-qumicas de las resinas intercambiadoras de iones. Y para el Sr, el Nd y el
Ra, se usaron resinas de extraccin. Para ms detalles ver el Anexo B: Resinas
cromatogrficas. Al final del proceso se obtuvo para cada muestra y cada uno de estos
elementos, una solucin acuosa del elemento (en HCl o en HNO3) purificada (casi sin otros
elementos). A estas soluciones se les llama fracciones de Ra, Sr, Nd, U y Th. Enseguida se
evapor completamente las fracciones. De tratarse de una solucin de HCl, se volvi a disolver
el residuo slido en HNO3 y se lo evaporar nuevamente, para eliminar los grupos cloruro (los
anlisis de espectrometra de masa se realizan sobre una solucin en cido ntrico puro).

5.4 Preparacin de filamentos

Fijacin a los pies de filamento


Para sujetar los filamentos de forma firme y segura se usa una pieza llamada pie de
filamento. Se trata de una placa metlica plana de aproximadamente 20x40x1 mm, con dos
bornes o electrodos aislados elctricamente de la placa y entre ellos. Existen pies simples y
dobles, y se los utiliza en dos tcnicas diferentes de ionizacin trmica (ver 5.5.3 ms abajo).

Los filamentos utilizados son de Re, de 40 m de espesor por 2 mm de ancho (vienen en


tiras de ~15 cm de largo). Su preparacin comienza por cortarlos en segmentos de 2 cm de
largo. Luego, con una pinza, se pone cada segmento en una herramienta especial y se lo dobla
en forma de C en ngulos rectos. Luego se fija el pie de filamento a la herramienta mediante
un tornillo, quedando las puntas del filamento sobre las puntas de los bornes del pie. Enseguida
se solda el filamento a los bornes con un arco elctrico. Antes de soldar, el pie de filamento
debe haber sido limado, para eliminar restos de Re (el Re no suelda sobre si mismo), lavado
(baos de agua con detergente y de agua destilada) y secado, para eliminar las molculas
orgnicas de la superficie de los electrodos (como se explica ms adelante en 5.5.3, estas
molculas pueden provocar interferencias isobricas).

Degasificacin
Consiste en calentar sucesivas veces los filamentos (haciendo pasar por ellos una corriente
de hasta 4,5 A) en una cmara de vacio (P < 210-7 bar). Ello vaporiza y elimina la mayor parte
de los elementos voltiles. Para lograr la presin requerida a partir de la presin atmosfrica, se
bombea el aire durante unas 6 a 7 horas.

30
U Th

Depositar una pequea gota de H3PO4 (cido fosfrico) al 2% Depositar una pequea gota de H3PO4 al 2%

Subir la intensidad de corriente hasta 0,7 A. Esperar 2 minutos Subir la intensidad de corriente hasta 0,7 A. Esperar 2 minutos

Depositar tres pequeas gota de grafito una a una a 0,5 A Depositar tres pequeas gota de grafito una a una a 0,5 A

Redisolver la muestra a 70C en 2,5 l de HNO3 1N Redisolver la muestra a 70C en 2,5 l de HNO3 1N
Depositar la muestra, la primera gota a 0,5 A, el resto a 0,8 A. Depositar la muestra, la primera gota a 0,5 A, el resto a 0,7 A.
Esperar 10 segundos. Esperar 10 segundos.
Subir progresivamente la intensidad de corriente hasta 1 A. Subir progresivamente la intensidad de corriente hasta 1,5 A.
Esperar 10 s. Esperar 10 s.
Subir progresivamente la intensidad de corriente hasta 1,5 A. Subir progresivamente la intensidad de corriente hasta 1,9 A.
Esperar 20 s. Esperar 10 s.

Bajar progresivamente la intensidad de corriente hasta 0 A. Bajar progresivamente la intensidad de corriente hasta 0 A.

31
Sr Ra

Depositar una pequea gota de Ta2O5 (xido de tantalio) al 2% Depositar una pequea gota de Ta2O5 al 2%

Subir la intensidad de corriente hasta 0,3 A. Esperar 5 minutos Subir la intensidad de corriente hasta 0,3 A. Esperar 2 minutos

Redisolver la muestra a 70C en HNO3 1N Redisolver la muestra a 70C en 2,5 l HNO3 1N

Tabla 3: Procedimiento de impregnacin de filamentos.


Depositar 1 l de la muestra a 0,8 A. Esperar 10 segundos. Depositar la muestra a 0,5 A. Esperar 2 minutos.
Subir progresivamente la intensidad de corriente, hasta ver un Subir progresivamente la intensidad de corriente hasta 1,5 A.
color rojo oscuro (~ 2,2 A). Esperar 10 s. Esperar 10 s.
Subir progresivamente la intensidad de corriente hasta 2,0 A.
Bajar progresivamente la intensidad de corriente hasta 0 A.
Esperar 5 s.
Depositar una pequea gota de Ta2O5 al 2% a 0,5 A, esperar
1 minuto.
Impregnacin con la muestra

Para analizar los isotopos de U, Th, Ra y Sr, se deposit las fracciones respectivas sobre
filamentos de Re previamente preparados, siguiendo un protocolo que depende del elemento a
analizar (ver Tabla 3). La tcnica consiste en depositar la fraccin en pequeas gotas sobre el
filamento, evaporando cada gota antes de depositar la siguiente (para ello se hace pasar una
corriente elctrica baja por el filamento, calentndolo). Sobre el centro del filamento se deposita
previamente un impregnante (xido de tantalio o cido fosfrico), para contener las gotas de la
fraccin sobre esa zona, y eventualmente un material poroso que acta como una matriz, en
cuyo volumen se esparcen los tomos del elemento a analizar. En el caso de uranio es
necesario usar una suspensin de grafito como matriz, para ionizar los tomos de uranio de
forma progresiva.

5.5 Espectrometra de masas

Las mediciones de proporciones y concentraciones isotpicas fueron obtenidas por el autor


de esta memoria, en el Laboratorio de Hidrologa y Geoqumica de Estrasburgo (LhyGeS:
Laboratoire d'Hydrologie et de Gochimie de Strasbourg), en la Universidad de Estrasburgo,
Francia. Se usaron dos tcnicas diferentes de espectrometra de masas, siendo la nica
diferencia entre ambas el mtodo de ionizacin.

Para los istopos de U, Th, Ra y Sr se us un espectrmetro Thermo Scientific modelo


Triton, que emplea la ionizacin trmica (TIMS = Thermal Ionization Mass Spectrometer),
midiendo seis cantidades diferentes en cada muestra7:
o 234U/235U, que permite determinar 234U/238U.
o 235U/233U, que sirve para determinar [U].
o 232Th/229Th, que sirve para determinar [Th].
o 230Th/232Th, que junto a [Th] y [U] permite obtener (230Th/238U).
o 226Ra/228Ra, que permite calcular [226Ra], y junto a [Th], determinar (226Ra/230Th).
o 87Sr/86Sr, que es directamente la cantidad de inters.
Las frmulas que usadas para calcular los valores de concentracin y razones de actividad
estn en el Anexo D: Clculo de razones isotpicas. Las frmulas de concentracin se basan en
la tcnica de dilucin isotpica, explicada en el Anexo C.

Por otra parte, los istopos de Nd fueron analizados con un espectrmetro de multi-
coleccin (MC-ICPMS) Thermo Scientific modelo Neptune, que emplea la ionizacin por plasma
acoplado inductivamente (ICPMS: Inductively Coupled Plasma Mass Spectrometer). El
instrumento entrega directamente la razn isotpica de inters 143Nd/144Nd.

Ambos espectrmetros utilizados fueron construidos por el mismo fabricante y exceptuando


la fuente de iones, tienen exactamente las mismas partes y la misma configuracin de
funcionamiento (ver Figura 5). Por ello, se presenta primero un esquema del funcionamiento
general del espectrmetro desde la fuente inica en adelante (los mtodos de ionizacin se
presentan en 5.5.3).

7
Aqu "muestra" se refiere a muestra original con "spike". Para ms detalles, referirse al Anexo C: Dilucin isotpica.

32
Figura 5: Esquema del funcionamiento de un espectrmetro de masa Thermo Scientific Triton (TIMS).

5.5.1 Trayectoria de los iones

Una vez generados, los iones son acelerados horizontalmente a grandes velocidades8
mediante una diferencia de potencial V0 de 10 kV aprox., hacia una cmara de vaco9 a travs
de un orificio de ~1mm de dimetro. Durante este primer tramo, los iones son focalizados por
lentes inicos (campos elctricos calibrados perpendiculares al campo acelerador), es la ptica
inica de la Figura 5.

Enseguida, los iones entran en una zona donde el campo elctrico externo es nulo y el
campo magntico es constante y vertical hacia arriba (inducido por un electroimn). El campo
magntico ejerce sobre los iones una fuerza proporcional a su intensidad B0, horizontal y
siempre perpendicular a la velocidad de stos (fuerza de Lorentz). Por lo tanto, los iones
describen un arco de circunferencia en este segundo tramo.

8
Aqu, los istopos analizados estn ionizados con una carga 1+ y tienen masas entre 86 y 238 u.m.a. aprox.. En
-1 5 -1
consecuencia, la diferencia de potencial de 10 kV los acelera a velocidades entre 9104 ms y 1,510 ms .
-1
Expresado de otro modo, entre 265 y 440 veces la velocidad del sonido en el aire seco a 25C (340 ms ).
9 -5
En esta parte, el instrumento debe mantener una presin de 10 mbar durante todo el procedimiento. En el caso
del MC-ICPMS la ionizacin es llevada a cabo a presin atmosfrica, y en consecuencia la cmara de vaco debe
ser bombeada durante todo el procedimiento. De lo contrario, los iones interactuaran con las molculas del aire
en su trayectoria.

33
10
Figura 6: Esquema de las trayectorias de los iones bajo la influencia del electroimn .

El radio de curvatura de dicho arco depender del valor de la velocidad del ion (a su vez
funcin de V0 y de la razn masa-carga), de B0 y de la razn masa-carga del ion. La magnitud
de la fuerza de Lorentz es qB0v, que es igual a mv2/R (fuerza centrpeta) por la segunda ley de
Newton. Si adems consideramos que en promedio los iones parten del reposo, la velocidad
promedio a la cual se aceleran los iones est dada, por igualar la energa cintica (mv2/2) con la
energa potencial (qV0). As, se obtienen dos ecuaciones (para R y v):

m v q m 2V0
R= , v 2V0 luego R
q B0 m q B0

En la Figura 6, los iones describen arcos de circunferencia debido a un campo magntico,


el cual es distinto de cero slo en el cuadrante superior derecho del plano (para x>0 e y>0), por
lo que las trayectorias son rectas del eje "y" hacia la izquierda y del eje "x" hacia abajo. El tramo
circular de cada trayectoria varia entorno de una trayectoria media radio de "R" y una apertura
angular de 90. Tomando en cuenta las relaciones geomtricas expuestas en la Figura 6, la
desviacin angular "" de un ion respecto del eje y est dada por

R m
arcsen 1 arcsen 1 donde =
R' ' q

De esta manera, suponiendo que la energa cintica de los iones tiene poca variacin, los
iones de un mismo istopo tendrn todos una desviacin angular similar respecto del eje y
(Figura 6), y llegarn al mismo detector como un haz de partculas (ver Figura 5). Garantizar la

10
El origen de coordenadas coincide con el centro del cuarto de circunferencia rojo (radio R), que representa la
trayectoria media de los iones. Dependiendo de su relacin masa-carga, cada ion es desviado de la trayectoria
roja. El arco de circunferencia negro (radio R' ) es una trayectoria deflectada hacia dentro (ion ms liviano).

34
focalizacin de ese haz de partculas, y as maximizar el nivel de la seal, es precisamente la
funcin de la ptica inica.

5.5.2 Deteccin

Los equipos aqu empleados cuentan con dos tipos de instrumentos de deteccin: La copa
de Faraday y el contador de iones.

Copas de Faraday
Tanto el espectrmetro Triton, como el Neptune cuentan con diez copas de Faraday
dispuestas como se muestra en la Figura 5. Son instrumentos mviles (tienen motores muy
precisos) y su posicin es calibrada cada vez para optimizar la intensidad de la seal. Se trata
de un electrmetro conectado a una jaula de Faraday (material conductor que envuelve un
volumen) hacia cuyo interior se dirigen los iones. Cada ion (o grupo de iones) que entra(n) en la
jaula puede verse como un impulso de corriente localizado en el tiempo, el que induce(n) un
campo electromagntico variable (por las ecuaciones de Maxwell). El campo elctrico provoca
una corriente de electrones sobre la superficie de la jaula proporcional a la corriente de iones en
su interior. Es la tensin (o fuerza electromotriz) debido a esa corriente elctrica la que mide el
electrmetro.

Contador de Iones
Ms all de la copa central, en el sentido de avance de los iones, se ubica el contador de
iones. Entre las copas y el contador, hay un filtro de masa de cuadrupolo (quadrupole mass
filter) que selecciona los iones en funcin de su relacin masa-carga. Esta seleccin permite
obtener peaks de seal mejor definidos en torno a la relacin masa carga fsica, y menos ruido
de fondo. El campo de cuadrupolo es generado por cuatro barras cilndricas de metal cargadas
elctricamente, dispuestas paralelamente en un arreglo cuadrado. Slo los iones con
determinadas razones masa-carga tendrn trayectorias estables entre las cuatro barras. Dichos
iones no impactan directamente al detector, sino a un multiplicador de electrones (SEM:
Secondary Electron Multipliyer)11. El resto de los iones son atrados a alguna barra o desviados,
lejos del instrumento de deteccin. La deteccin se realiza finalmente en forma indirecta,
mediante el conteo de los pulsos de electrones secundarios (un pulso elctrico = un ion). Si bien
el contador de iones es capaz de percibir seales mucho ms dbiles que un copa de Faraday,
no est exento de limitaciones. Por ejemplo, la resolucin temporal de conteo del SEM es finita,
e igual al "tiempo muerto" instrumental (tiempo mnimo necesario para "resolver" o diferenciar
dos iones que se descargan sobre el SEM muy cerca en el tiempo).

Los instrumentos fueron usados de dos maneras diferentes para medir razones isotpicas:

A. Modo multicoleccin: consiste en focalizar dos o ms haces de iones (cada uno


corresponde a un istopo) en la misma cantidad de copas de Faraday, y medir la razn
isotpica a partir de las razones entre las respuestas instrumentales (en Volts). Como se

11
Un multiplicador de electrones es un tubo curvo metlico, con forma de cono o cuerno, con un alto voltaje opuesto
en carga al ion que se mide. Cada descarga de un ion en el metal cargado, produce la emisin de varios
electrones secundarios dentro del tubo, que a su vez impactan el interior del multiplicador y provocan la emisin
de an ms electrones secundarios, y as sucesivamente. Segn Jarvis et al. (1992) al final del tubo saldran del
8
orden de 10 electrones por cada ion que entra al tubo.

35
explica en 5.5.1, el espectrmetro de masa separa los iones segn su masa de una
manera que recuerda a un prisma que descompone un haz de luz blanca en sus
componentes segn longitud de onda.

B. Modo SEM/Faraday: se focaliza el haz de iones de un istopo al SEM, y se mide (en


nmero de tomos), y luego se dirige el haz del otro istopo hacia la copa central, y se
mide (en Volt). Antes de la medida, debe calibrarse la respuesta instrumental del
contador de iones respecto de la copa central, para que el computador que controla en
instrumento pueda calcular la razn isotpica.

5.5.3 Tcnicas de Ionizacin

Ionizacion trmica (TI = Thermal Ionization)

Como se explica en el 5.2.4 (Preparacin de filamentos), la muestra est impregnada en la


superficie de un filamento de un metal muy refractario (Re). Para calentar el filamento (hasta
~2000-2500 C), se hace pasar una corriente elctrica a travs de l. Debido a la alta
temperatura, parte de los tomos de la muestra se ionizan, el resto ebulle antes de recibir
suficiente energa para su ionizacin.

Una de las desventajas de esta tcnica es su lmite de deteccin, por ejemplo de ~20 ppm
para el 238U, siendo ampliamente menos sensible que el ICP-MS con ~10 ppt para el mismo
istopo. Tambin cabe mencionar, que no todos lo elementos son ionizados fcilmente con esta
tcnica, restringiendo su aplicabilidad a aquellos elementos y compuestos cuyo potencial de
ionizacin sea relativamente bajo. Adems, la ionizacin no es igualmente eficiente para todos
los elementos, siendo generalmente menor al 1%. Si bien es importante corregir por
fraccionamiento isotpico experimental (particularmente en la fuente inica), este efecto es en
promedio menor y ms consistente que en otras tcnicas.

Plasma Acoplado Inductivamente (ICP = Inductively Coupled Plasma)

Primero se nebuliza la muestra (en forma de solucin en HNO3 0,75N), separndola en


pequeas gotas, gracias a su aspiracin (un flujo de argn de alta velocidad produce la
diferencia de presin necesaria). Enseguida el aerosol pasa a una cmara donde las gotas ms
gruesas son drenadas (Jarvis et al., 1992), saliendo de ella tpicamente slo el ~2% ms fino
(Olesik, 1996). Esta nube de finas gotas es incorporada a un flujo de argn mantenido entre tres
tubos concntricos (antorcha), que termina en una llama de plasma12 muy caliente (hasta
~10000 K en su interior). El plasma de argn es afectado por un campo electromagntico
variable (usualmente ~27 MHz), que es generado por una bobina de induccin y es la fuente de
energa del plasma. Se dice que el plasma est acoplado al campo electromagntico. El plasma
caliente separa los compuestos presentes en la solucin original en tomos individuales y luego
los ioniza en una proporcin significativa.

12
Un plasma es un gas elctricamente conductor.

36
Tcnicas Ventajas Desventajas
1. Lmite de deteccin muy bajo 1. Fraccionamiento de masa
(~10 ppt). instrumental mayor.
Por
2. Capacidad de ionizar todo y 2. Ocurrencia de interferencia
plasma
mayor eficiencia de ionizacin. isobrica.
(ICP)
3. Medicin rpida y semi- 3. La contaminacin del Ar impide
automatizada. analizar ciertos elementos (p.e. Sr)
Ionizacin

1. Lmite de deteccin de ~20 ppm


1. Los iones producidos tienen un para el 238U.
rango de energa acotado, y no 2. No todos lo elementos pueden
hay necesidad de usar un filtro ser ionizados.
Termal (TI) de energa. 3. La eficiencia de ionizacin es
2. El uso de soluciones mono- ms baja (generalmente <1%)
elementales evita la interferencia 4. La preparacin de los filamentos
isobrica. y su anlisis demandan mucho
trabajo y tiempo.
Multiplica- 1. Medicin directa y absoluta de No se puede usar para medir
dor de los iones muestras concentradas (se satura
Deteccin

electrones 2. Lmite de deteccin muy bajo. y se daa).

La estabilidad de las condiciones Necesita de un gran nmero de


Copas de
de operacin permite mayor iones incidentes para detectar una
Faraday
precisin (mtodo ms robusto). seal.

Tabla 4: Comparacin del instrumental de espectrometra de masas


Una de las grandes ventajas de esta tcnica es su lmite de deteccin extremadamente
bajo, pudiendo llegar a menos de 10 ppt. Sin embargo, uno de los inconveniente de usar
plasma para la ionizacin es que produce fraccionamiento de masa instrumental", debiendo
corregir los resultados por este efecto. Por otro lado, los istopos de Sr no pueden medirse con
esta tcnica, dado que en el Ar del plasma hay grandes cantidades de Rb, y el 87Rb interfiere
con las mediciones de 87Sr por tener casi la misma masa que ste ltimo.

5.5.4 Procedimiento de medicin

TIMS Triton

1. Aumentar progresivamente la intensidad de corriente que pasa por el filamento para


calentarlo, hasta que ste comienza a emitir iones13. En ese momento abrir la compuerta
que separa la fuente de iones del resto del instrumento.
2. Calibrar la posicin de las copas de Faraday para los istopos a medir14.
3. Focalizar: optimizar el nivel de seal en funcin de los parmetros de la ptica inica
(tensiones).
4. Calibrar la relacin entre razn masa-carga y magnitud del campo magntico (B0)15.

13
En la prctica se espera que el filamento alcance una determinada temperatura (medida con un piranmetro).
14
Para cada grupo de istopos (U, U-Th, Ra, Sr) las posiciones ptimas de las copas han sido pre-calibradas y
guardadas en archivos en el computador que controla en instrumento.

37
5. En modo SEM/Faraday debe verificarse que el nivel de ruido electrnico del SEM sea bajo
(menos de 7-10 cuentas por minuto), y que no haya masas parsitas (peaks de seal en
masas atpicas), que son indicio de eventuales interferencias isobricas.
6. En modo SEM/Faraday, se necesita calibrar el contador de iones respecto de las copas.
7. Lanzar la medicin, tpicamente 200 ciclos de medicin, de lo cual se hace una estadstica
(promedio y error estndar interno). La seal debe mantenerse estable durante todo el
procedimiento, o al menos nunca bajar de un cierto nivel mnimo (100 cuentas por segundo
para el SEM).

ICP-MS Neptune

La medicin de 143Nd/144Nd, es mucho ms simple, porque no hay filamento y la seal es se


estabiliza por si sola. El procedimiento es casi totalmente automatizado y dura ~10 minutos por
muestra en vez de ~2 horas por muestra como en el Triton. Concretamente se realizan los
pasos 2-7 recin detallados, sin necesidad de monitorear la temperatura de un filamento ni el
nivel de la seal.

15
Al variar B0 barriendo un cierto intervalo de valores, el haz de iones se deflecta pasando de un lado a otro del
detector. El mximo de seal debera ocurrir justo en valor de B0 correpondiente a la razn masa-carga fsica del
ion (hay una equivalencia entre B0 y m/q). Esta condicin permite calibrar la relacin B0 vs. m/q.

38
CAPTULO 6: Resultados

6.1 Muestras
Las muestras analizadas corresponden a sedimentos fluviales de granulometra media a
fina (< 4mm) tomados del lecho del ro uble (cdigos CH08-01 a 06 excepto 03, y CH08-09) e
Itata (cdigo CH08-08), y a muestras de agua tomadas en los mismos lugares (cdigos CH08-
01 a 09, salvo 02 y 06); ver Figura 2. El muestreo se realiz en distintos puntos del cauce del
ro, en la vertiente oriental de la Cordillera de la Costa y en la Cordillera Principal (Figura 2). La
muestra CH08-03 fue tomada de un pequeo tributario del ro uble sobre roca grantica. Junto
con la muestra CH08-08 son las nicas que no fueron tomadas del lecho del ro uble.

6.2 Elementos mayores y traza (ICP-AES, ICP-MS)


Estos anlisis se realizaron sobre muestras molidas y fundidas las que fueron sometidas a
espectrometra de emisin atmica (para los elementos mayores) y espectrometra de masas
(para los elementos traza).

En las Tablas 5 y 6 se muestran las concentraciones de los elementos mayores en las


muestras de sedimentos (en % de peso) y de agua (en ppm) respectivamente. El indice LOI
(Loss On Ignition) corresponde a la prdida de los componentes voltiles de las muestras.

LOI SiO2 Al2O3 MgO CaO Fe2O3 MnO TiO2 Na2O K2O P2O5
110C 1000C % peso

LD - - 0,3 0,02 0,008 0,01 0,02 0,0001 0,001 0,05 0,004 0,06
CH08-01 0,90 2,56 58,28 15,14 3,32 4,77 9,15 0,15 1,17 3,07 1,31 0,18
CH08-02 0,90 2,65 58,94 15,65 3,02 4,79 7,96 0,15 1,00 3,36 1,41 0,16
CH08-03 0,81 2,61 64,17 15,25 1,85 3,26 5,76 0,10 0,75 3,08 2,24 0,13
CH08-04 0,83 2,33 57,75 15,18 3,23 4,87 9,67 0,15 1,25 3,19 1,38 0,18
CH08-05 0,83 2,38 59,03 15,68 3,03 4,86 8,13 0,14 1,05 3,29 1,39 0,17
CH08-06 0,56 1,68 62,69 15,14 2,31 4,43 6,56 0,12 0,87 3,63 1,80 0,21
CH08-09 0,66 1,70 63,76 14,93 2,20 4,26 6,05 0,10 0,82 3,43 1,93 0,15
CH08-08 0,40 0,73 63,56 14,69 2,95 5,47 6,38 0,11 0,83 3,12 1,61 0,14

LD = Lmite de deteccin
Tabla 5: Concentraciones de elementos mayores en las muestras de sedimentos del ro uble (% peso).

Las muestras de sedimentos contienen principalmente slice (58-64%), Al2O3 (15%) y Fe2O3
(6-10%) junto a menores contenidos de CaO (3-5%), MgO (2-3%) y Na2O (3%). El orden
relativo de abundancia de estos elementos se mantiene en las distintas muestras. Esto
corresponde aproximadamente a la composicin de rocas corticales de composicin intermedia
(andestica) en concordancia con la composicin de las rocas a partir de las cuales se formaron
los sedimentos del ro uble, an cuando no se cuente con datos geoqumicos de esas rocas.

39
Se trata de rocas volcnicas andesticas a baslticas y en menor medida granitos de la
Cordillera Principal, como tambin aquellas formaciones sedimentarias de la Depresin
Intermedia producto de la erosin de las primeras (ver 4.1).

COD K Na Ca Mg Fe Zn Cu F- Cl- SO42- NO3-


ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm
CH08-01 0,47 0,32 2,18 4,56 0,91 <LQ <LQ <LQ 0,03 0,51 2,45 0,52
CH08-03 0,42 0,12 1,16 2,70 0,43 <LQ <LQ <LQ 0,01 0,28 2,56 0,34
CH08-04 0,57 0,34 2,21 4,36 0,99 <LQ <LQ <LQ 0,02 0,54 2,74 0,33
CH08-05 0,51 0,37 2,31 4,44 1,04 <LQ <LQ <LQ 0,02 0,58 2,83 0,55
CH08-07 5,63 1,75 6,71 7,77 3,64 <LQ <LQ <LQ 0,05 2,42 4,10 -
CH08-09 1,50 0,90 3,76 5,02 1,89 <LQ <LQ <LQ 0,03 1,49 3,44 0,84
CH08-08 3,58 1,85 12,57 5,44 2,43 <LQ <LQ <LQ 0,03 5,41 12,97 -

LQ = Lmite de cuantificacin

Tabla 6: Concentraciones de elementos mayores en las muestras de agua del Ro uble (ppm).

La concentracin de elementos mayores en las muestras de agua presenta mayor


variabilidad respecto de las muestras de sedimento. Los elementos ms abundantes en el agua
son el Ca, Na y Mg, los que de acuerdo al anlisis de aniones se encontraran en forma de
sulfatos, nitratos o cloruros. Esto contrasta con la composicin de los sedimentos, que
contienen principalmente Si, Al y Fe, y en menor medida Ca, Mg, Na y K. que podran
encontrarse en forma de xidos u otros compuestos insolubles o poco solubles.

En la Tabla 7 se muestra el contenido de elementos traza en las muestras de sedimentos.

Sc V Cr Co Ni Cu Zn As Rb Sr Y Zr
ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm
LD* 0,40 0,30 0,40 0,07 4,00 2,00 20,00 0,20 0,70 1,00 2,00 0,40
CH08-01 25,87 173,50 49,96 20,28 35,69 68,86 140,74 1,62 35,49 315,64 24,88 153,06
CH08-02 23,21 149,98 30,91 19,48 30,53 62,67 129,27 1,50 38,89 322,56 28,36 150,82
CH08-03 18,69 115,28 18,73 12,62 19,26 48,90 93,78 1,22 72,48 243,99 20,74 165,06
CH08-04 25,65 204,35 50,81 20,43 28,79 42,96 111,12 1,84 39,14 313,58 24,29 155,28
CH08-05 23,13 164,06 47,47 18,72 24,78 38,90 91,06 1,54 38,95 328,23 23,77 151,90
CH08-06 20,77 128,91 28,40 13,73 49,73 30,67 79,57 1,93 51,92 324,77 26,21 171,48
CH08-09 17,20 115,01 21,13 13,20 20,07 26,19 70,25 1,71 53,84 323,15 21,83 163,64
CH08-08 19,11 145,68 44,84 16,86 78,34 48,26 102,77 1,61 46,37 376,60 18,39 116,90

Nb Mo Sb Cs Ba Hf Ta W Pb Th U
ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm
LD* 0,10 0,09 0,04 0,60 0,50 0,06 0,00 0,10 0,20 0,02 0,00
CH08-01 4,91 0,63 0,68 1,66 329,34 3,86 0,37 5,06 17,05 3,58 1,04
CH08-02 4,62 0,72 0,58 2,06 349,03 3,84 0,35 4,88 14,66 3,89 1,11
CH08-03 4,02 0,94 0,92 2,98 465,20 4,28 0,33 4,56 15,52 5,44 1,56
CH08-04 5,14 0,76 0,70 2,16 332,48 3,91 0,37 4,54 12,58 3,85 1,16
CH08-05 4,50 0,69 0,66 2,17 342,55 3,79 0,34 4,35 10,28 3,92 1,08
CH08-06 5,34 0,96 0,66 2,63 428,17 4,34 0,41 6,66 10,49 4,83 1,29
CH08-09 4,87 0,69 0,54 2,67 487,08 4,04 0,38 4,83 10,63 4,54 1,24
CH08-08 3,58 0,60 0,18 1,68 420,98 2,97 0,27 4,72 9,19 3,14 0,82

40
La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu
ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm
LD* 0,04 0,03 0,00 0,02 0,01 0,01 0,30 0,10 0,01 0,002 0,01 0,003 0,001 0,002
CH08-01 16,42 33,22 4,52 19,01 4,47 1,38 3,80 0,61 4,26 0,93 2,38 0,39 2,54 0,41
CH08-02 15,27 33,51 4,53 19,55 4,82 1,49 3,95 0,70 4,69 1,03 2,60 0,43 2,65 0,41
CH08-03 15,45 32,22 4,16 16,45 3,60 1,06 3,24 0,51 3,54 0,77 1,92 0,34 2,23 0,36
CH08-04 16,18 34,44 4,56 18,87 4,41 1,33 3,74 0,61 4,06 0,90 2,26 0,39 2,40 0,38
CH08-05 14,82 31,42 4,29 18,12 4,09 1,27 3,46 0,61 3,90 0,84 2,10 0,36 2,19 0,37
CH08-06 18,55 38,52 5,15 21,20 4,83 1,39 4,20 0,70 4,49 0,99 2,43 0,42 2,60 0,42
CH08-09 15,45 32,99 4,41 17,90 4,09 1,21 3,44 0,60 3,83 0,83 2,03 0,35 2,18 0,35
CH08-08 12,71 26,51 3,57 15,03 3,37 1,10 3,03 0,49 3,09 0,66 1,66 0,29 1,78 0,27

* Lmite de deteccin

Tabla 7: Concentraciones de elementos traza en las muestras de sedimentos.

La siguiente tabla presenta la concentracin en elementos traza (y algunos mayores como


el Na, Mg, Si, Ca) de las muestras de agua obtenidos mediante espectrometra de masas (ICP-
MS).
CH08-01 CH08-03 CH08-04 CH08-05 CH08-07 CH08-08 CH08-09
m m m m m m m
ppb %
ppb %
ppb %
ppb %
ppb %
ppb %
ppb %
Li 0,94 2,9 0,18 8,6 0,91 1,3 0,90 2,5 0,34 0,8 1,02 5,1 0,65 6,3
B 9,60 1,5 1,36 3,5 9,66 1,5 10,33 1,3 18,32 0,9 23,44 4,3 14,11 1,2
Na 2094 2,2 1115 1,0 2117 2,0 2244 1,4 6753 0,7 10460 1,7 3460 1,7
Mg 461,8 1,3 241,4 2,7 506,5 1,0 532,7 0,9 2142 1,7 1268 0,2 941,6 1,8
Al 9,410 1,0 5,393 1,9 9,87 1,6 10,40 0,4 14,40 4,3 11,09 0,6 6,559 0,2
Si 5978 1,2 3633 2,4 5960 2,6 6271 0,2 8701 1,1 8261 0,8 7210 0,6
K 372,4 1,2 191,2 3,7 386,3 2,3 414,3 1,2 1600 0,9 1319 1,6 801,1 0,2
Ca 5030 0,8 2997 1,1 4885 1,1 4995 1,0 9163 0,9 6226 0,2 5709 0,5
Sc 1,185 1,5 0,711 2,2 1,10 1,8 1,163 1,7 1,561 0,9 1,476 1,4 1,317 0,4
V 1,640 1,8 0,315 1,3 1,44 2,3 1,595 0,3 4,117 0,7 2,477 1,4 1,803 1,9
Cr 0,061 15 0,053 7,6 0,049 4,6 0,064 9,5 0,097 6,0 0,058 3,0 0,049 17
Mn 0,891 0,9 0,176 2,6 0,781 1,1 0,822 1,4 3,885 1,2 9,396 1,9 5,419 0,2
Fe 5,336 4,1 5,205 4,7 5,98 4,1 5,658 3,5 28,50 12 5,531 4,9 7,004 0,7
Co 0,010 19 0,008 25 0,008 21 0,008 49 0,211 5,4 0,033 3,1 0,030 9,5
Ni 0,035 35 0,030 43 0,021 3,8 0,023 33 0,383 2,6 0,099 3,5 0,047 19
Cu 0,204 2,1 0,161 5,8 0,197 5,1 0,193 2,3 2,840 0,9 0,645 2,2 0,531 2,6
Zn 2,586 2,4 3,639 1,6 4,65 5,2 4,316 1,7 4,016 0,8 3,850 1,0 2,785 3,2
As 1,689 3,5 0,115 7,9 1,57 4,8 1,601 5,6 1,596 1,8 1,292 2,4 1,167 2,7
Rb 0,935 2,3 0,561 1,3 0,933 0,7 0,979 0,3 6,680 1,5 3,675 1,8 1,941 1,0
Sr 18,89 0,9 13,05 0,8 19,05 1,9 19,89 0,8 77,00 1,0 38,81 2,7 30,42 1,0
Yb 0,0157 8,0 0,0204 12 0,0142 6,1 0,0145 8,0 0,0902 3,2 0,0145 9,0 0,014 6,0
Zr 5,56E-3 47 6,43E-3 14 2,86E-3 34 2,3E-3 23 0,0906 4,6 0,0041 9,7 5,8E-3 26
Mo 0,2802 2,0 0,3029 4,8 0,262 1,8 0,275 4,0 0,339 5,3 0,317 4,8 0,280 4,9

41
Cs 0,0335 2,6 0,0340 9,3 0,0258 2,4 0,0230 10 0,0139 4,1 0,0314 6,3 0,020 5,5
Ba 1,641 0,7 1,373 3,0 1,477 2,0 1,521 1,6 20,06 0,9 4,84 0,9 3,310 3,6
La 7,53E-3 12 8,89E-3 12 5,99E-3 20 7,16E-3 4,1 0,0285 6,5 0,0054 11 0,006 16
Ce 8,99E-3 4,5 6,68E-3 26 8,40E-3 9,2 9,68E-3 14 0,0597 6,7 0,0077 17 0,0095 6,6
W 0,0405 2,8 0,0140 6,9 0,0343 13 0,0379 10 0,0064 29 0,0197 12 0,0208 0,3
Pb 0,0106 15 0,0135 12 0,0097 13 8,24E-3 7,5 0,0154 6,6 0,0071 8,2 0,0045 31
Th 9,97E-6 >100 6,20E-4 26 3,38E-4 85 -6,8E-5 >100 1,67E-3 1,1 -3,4E-4 47 7,2E-5 >100

U 4,86E-3 24 4,21E-3 6,3 4,70E-3 22 4,32E-3 9,6 0,0419 3,0 0,0099 5,7 0,0074 20

Tabla 8: Concentraciones de elementos traza y algunos mayores en muestras de agua del ro uble
con sus respectivos errores estndar de medicin en % (Std), obtenidas a partir de anlisis ICP-MS. Los
errores estndar por sobre el 15% estn destacados en rojo. Se descartaron los elementos cuyos errores
estndar son mayores que 15% para todas las muestras, por ejemplo Nb y Eu.

Para extraer mayor informacin de los datos de elementos traza (Tablas 7 y 8), es
necesario normalizarlos dividiendo por un elemento inmvil o poco mvil. En el Anexo A, se
muestran grficos de elementos normalizados respecto de Al y Th, tomando en cuenta que en
general son elementos inmviles o menos mviles que el resto.

6.3 Mineraloga por difraccin de rayos X

La mineralogia de las muestras en estudio (principalmente cuarzo y plagioclasas) es


consistente con la hiptesis del origen preponderantemente volcnico de los sedimentos. La
ausencia de feldespato-K (salvo por la muestra CH08-03, tomada sobre un batolito grantico)
sugiere que el aporte de los granitos es minoritario.
Feldespato K
Plagioclasas

Magnetita

Magemita
Hematita
Piroxeno

Pirofano
Anfbola

Ceolitas
Cuarzo

Arcillas
Muestra

CH08-01 XX Xx X X X X X? X
CH08-02 XX Xx X X X X X? X
CH08-03 XX Xx X X x? X X
CH08-04 XX Xx x? X X?
CH08-05 XX Xx X? x? X X X? X? X
CH08-06 XX Xx x? X X X? X
CH08-08 XX Xx x? X X? X
CH08-09 XX Xx x? X X X

XX: Fase mayoritaria


X: Fase minoritaria

Tabla 9: Interpretacin de los difractogramas de rayos X en roca total de sedimentos del ro uble.
Las muestras fueron molidas (< 100 m) previamente. Los signos de interrogacin "?" indican la
sospecha de la presencia del mineral respectivo.

42
En la fraccin arcillosa se encuentra clorita, kaolinita, illita y smectita. Las arcillas provienen
probablemente de la meteorizacin de los feldespatos (plagioclasas esencialmente) y otros
silicatos (piroxenos, anfbolas).

Illita/Smectita

Plagioclasas
Smectita
Kaolinita

Anfbola

Ceolitas
Cuarzo
Clorita
Muestra

Illita
CH08-01 + + + + + + + +
CH08-02 + + + + + + + +
CH08-03 + + + + + + +
CH08-04 + + + + + +
CH08-05 + +? + + + + + +
CH08-06 + + + + +
CH08-08 + + + + +
CH08-09 + +? + + + +

Tabla 10: Interpretacin de difractogramas de las fracciones arcillosas (<2 m): muestras del ro uble.
Las cruces indican la presencia de cada mineral. Los signos de interrogacin "?" indican la sospecha de
la presencia del mineral respectivo.

6.4 U, Th y Ra: Razones de actividad y concentraciones (TIMS)

Enseguida se presentan las razones de actividad entre istopos de la serie del 238U
calculadas a partir de los datos de espectrometra de masas por ionizacin trmica16. La
posicin geogrfica relativa de las muestras, de este a oeste, se indica por las letras a-h. Los
desequilibrios de 226Ra son mayores que los del 230Th o el 234U, que se mantienen relativamente
cerca de la condicin de equilibrio secular.

La Tabla 12 muestra las concentraciones de U, Th y Ra en las muestras de arena


determinadas por dilucin isotpica (ver anexo C). CH08-Bl corresponde a un control de
contaminacin debido a las manipulaciones qumicas (blanco). Los casilleros vacos
corresponden a anlisis fallidos.

16 235 238 230 232


Las cantidades medidas directamente fueron las siguientes razones isotpicas U/ U, Th/ Th .... y
226
usando la tcnica de dilucin isotpica se pudo calcular [U], [Th], [ Ra] y las razones de actividad de inters.

43
234 238 230 232 230 234 226 230
( U/ U) ( Th/ Th) ( Th/ U) ( Ra/ Th)
Muestra
Media m Media m Media m Media m

CH-08-01 a 0,9948 0,0029 0,8301 0,0025 1,0116 0,0047 1,00042 0,0076


CH-08-02 b 0,9981 0,0030 0,8047 0,0020 1,0054 0,0045 1,02875 0,0075
CH-08-03 c 1,0115 0,0030 0,6851 0,0017 1,0191 0,0045 1,08686 0,0078
CH-08-04 d 1,0058 0,0028 0,8186 0,0020 0,9861 0,0042 1,08201 0,0084
CH-08-05 e 0,9980 0,0035 0,7803 0,0022 0,9897 0,0049 1,04525 0,0076
CH-08-06 f 1,0083 0,0028 - - - - - -
CH-08-08 h 0,9960 0,0027 0,7915 0,0023 1,0222 0,0044 1,06534 0,0078
CH-08-09 g - - - - - - 1,04637 0,0084
CH-08-09 dup g 1,0073 0,0029 0,7655 0,0022 1,0202 0,0047 1,04637 0,0084
238 234 230 226 232
Tabla 11: Razones de actividad entre nucleidos de la serie del U-238 ( U, U, Th y Ra) y Th
en muestras de sedimentso del ro uble. Las razones de actividad destacadas en negrilla cuantifican
directamente los desequilibrios radioactivos en la serie del U-238 (razn de actividad igual a 1,0 equivale
a equilibrio secular). A la derecha del cdigo de cada muestra se indica mediante una letra su posicin
relativa en el curso del ro de este a oeste (a=ro arriba, h = ro abajo).

226
[U] ppm [Th] ppm [ Ra] fg/g
Muestra
Media m Media m Media m
CH-08-01 a 0,9616 0,0025 3,5375 0,0032 326,6 2,2
CH-08-02 b 0,9742 0,0025 3,6862 0,0044 339,3 2,2
CH-08-03 c 1,2136 0,0034 5,5405 0,0049 458,7 2,9
CH-08-04 d 1,0054 0,0025 3,6962 0,0034 364,0 2,6
CH-08-05 e 0,9849 0,0025 3,7830 0,0042 343,0 2,2
CH-08-06 f 1,1497 0,0028 - - 417,8 2,7
CH-08-08 h 0,7187 0,0016 2,8052 0,0021 263,0 1,7
CH-08-09 g - - - - 376,2 2,4
CH-08-09 dup g 1,0470 0,0026 4,2652 0,0056 379,8 2,7
226
[U] pg [Th] pg [ Ra] pg
CH-08-Bl Media m Media m Media m
16,7 0,1 87 2 - -
226
Tabla 12: Concentraciones de U, Th y Ra en muestras de sedimento del ro uble obtenidas por TIMS.
A la derecha del cdigo de cada muestra se indica mediante una letra su posicin relativa en el curso del
ro de este a oeste (a=ro arriba, h = ro abajo).

6.5 Razones isotpicas de Sr y Nd (TIMS y MC-ICPMS)

Las razones isotpicas de Sr (87Sr/86Sr) y Nd (143Nd/144Nd) obtenidas para los sedimentos


(ver Tabla 13) presentan valores muy homogneos, con una ligera tendencia a valores ms
radiognicos ro abajo (descartando el puntos anmalo "c"). Por tratarse de istopos estables o
con vidas medias muy largas (miles de millones de aos), estos datos permiten identificar el
orgen de los sedimentos.

44
Muestra
87 86
Sr/ Sr 2m
143
Nd/
144
Nd Nd 2m
-6
10
CH08-01 a 0,703924 5 0,512826 3,667 0,1
CH08-02 b 0,703919 4 0,512827 3,687 0,1
CH08-03 c 0,704079 4 0,512826 3,667 0,1
CH08-04 d 0,703927 3 0,512817 3,492 0,1
CH08-05 e 0,703934 5 0,512835 3,843 0,1
CH08-06 f 0,703979 4 0,512817 3,492 0,1
CH08-08 h 0,704510 4 0,512780 2,770 0,1
CH08-09 g 0,704284 5 0,512819 3,531 0,1
CH08-09dup g 0,704288 5 0,512827 3,687 0,1

Tabla 13: Razones isotpicas de Sr y Nd en las muestras de sedimentos obtenidas por TIMS e ICP-MS.
CH08-09 y CH08-09dup son dos mediciones de la misma muestra (prueba de reproductibilidad). A la
derecha del cdigo de cada muestra se indica mediante una letra su posicin relativa en el curso del ro
de este a oeste (a=ro arriba, h = ro abajo).

45
CAPTULO 7: Discusin

7.1 Elementos mayores y traza e istopos de Sr y Nd

A continuacin se describe la evolucin qumica e isotpica de los sedimentos del ro uble


desde los Andes hacia el mar, y enseguida se discute e interpreta.

En el Anexo A se presentan grficos de concentraciones de elementos qumicos


normalizados respecto a elementos inmviles (Al y Th)17, en funcin de la distancia a lo largo
del thalweg (camino seguido por el curso del ro). La Figura 7 contiene una seleccin de
algunos de dichos grficos. Los puntos verdes representan muestras tomadas del lecho del ro
uble. Estos puntos muestran que ro abajo hay una disminucin del contenido relativo de Mg,
Fe, Ca, Mn, Co, Sc, V, Cr y Ti, un aumento de Ba, K y Rb (ver Figura 7). Tambin hay una
variacin leve a nula (dentro de las barras de error) del contenido relativo de U, Th, Si, Na, Sr,
Nd y Mo. La muestra CH08-03 (punto rojo a la izquierda en los grficos del Anexo A) fue tomada
en un pequeo tributario del uble sobre roca grantica (anotacin de terreno), y es ms pobre
en Fe, Mg, Ca, Mn, Co, y otros metales, y ms rica en K, Ba y Mo que otras muestras
geogrficamente cercanas. La muestra CH08-08, tambin marcada en rojo, fue tomada ro
abajo de la confluencia con el ro Itata y tambin presenta valores anmalos.

La composicin de elementos mayores en los sedimentos es cercana a la composicin de


las andesitas, lo cual es consistente con la geologa de la cuenca del ro uble: las posibles
fuentes sedimentarias son andesitas, andesitas baslticas y basaltos, y en menor medida
granitos. El perfil de tierras raras (Lantnidos) apoya esta hiptesis (ver Figuras 8 y 9), puesto
que ste es muy similar, incluso en sentido absoluto (concentracin), al perfil de tierras raras de
rocas de volcanes activos geogrficamente cercanos (Antuco y Laguna del Maule). La anomala
de Eu que se da en algunos granitos y pegmatitas no se observa en la Figura 9.

Figura 7: Elementos mayores y traza normalizados respecto a Th y Al2O3.


(contina en las 3 pginas siguientes...)

17
Debido al espacio que ocupan los grficos (14 pginas), se decidi inclurlos fuera de texto, en un anexo.

46
Figura 7: (...continuacin: 1/3) Elementos mayores y traza normalizados respecto a Th y Al2O3.

47
Figura 7: (...continuacin 2/3) Elementos mayores y traza normalizados respecto a Th y Al2O3.

48
Figura 7: (...continuacin: 3/3) Elementos mayores y traza normalizados respecto a Th y Al2O3.

49
Figura 8: Perfiles de tierras raras de sedimentos (semilog) normalizados respecto a rocas volcnicas
de la vertiente occidental de los Andes entre los 36 y 37,5 de latitud S (promedio sobre datos de los
volcanes Antuco y Laguna del Maule, tomados de Hickey et al. (1986)). En en eje "X", 1=La, ..., 14=Lu.

Figura 9: Perfiles de tierras raras de muestras del ro uble normalizados resp. a basaltos ocenicos
(promedio global). Los datos de normalizacin se tomaron de Arvalo y McDonough (2010). En en eje
"X", 1=La, ..., 14=Lu, el eje "Y" est en escala logartmica.

50
Muestra La/Tb La/Ho La/Er La/Tm
CH08-01 4,24 4,09 4,68 4,05
CH08-02 3,45 3,43 3,98 3,44
CH08-03 4,77 4,65 5,45 4,47
CH08-04 4,18 4,17 4,87 4,05
CH08-05 3,83 4,09 4,79 4,01
CH08-06 4,19 4,35 5,18 4,29
CH08-09 4,07 4,29 5,17 4,24
CH08-08 4,09 4,48 5,20 4,30

Tabla 14: Proporciones entre La y tierras raras pesadas (Tb, Ho, Er y Tm) en sedimentos del ro uble.
Notar que la muestra destacada en negrilla (CH08-03) presenta valores mximos sistemticamente,
sugiriendo mayor aporte de granitos para esa muestra.

La contribucin de los granitos como fuente sedimentaria sera relativamente baja, ya que
los granitos contienen significativamente menos Mg, Fe y Ca y ms slice que las andesitas y
presentan un perfil de tierras raras ligeramente distinto, marcadamente ms enriquecido en
tierras raras livianas. En ese mismo sentido, la Tabla 14 indica que la muestra CH08-03 recibi
aportes de una fuente enriquecida en tierras raras livianas (La, Ce, Pr, Nd, Pm y Sm) respecto
de las pesadas (Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb y Lu), o bien empobrecida en tierras raras
pesadas, posiblemente el Batolito Grantico de la Cordillera.

87 86
Figura 10: Nd vs. Sr/ Sr: comparacin con sedimentos del ro Amazonas. Tomado de Viers et al. (2008)
quienes a su vez lo adaptaron de Roddaz et al. (2005). Los puntos rojos encerrados en el valo naranja
representan las muestras de sedimentos del ro Itata.

51
La composicin de istopos de Sr (87Sr/86Sr) y Nd (143Nd/144Nd) se mantiene dentro de un
rango pequeo si se compara con la evolucin de sedimentos en la cuenca del ro Amazonas
(ver Figura 10), lo cual es coherente con que los sedimentos de ese ro atraviesan distancias
mucho mayores que los sedimentos del ro Itata, pasando probablemente por un perodo de
maduracin mucho ms largo, de una intensidad similar o mayor. En la Figura 10 la
composicin de las muestras del ro Itata est contenida en una pequea zona del campo ms
"primitivo" (ms cercano a la roca madre: volcanitas andesticas a baslticas). En la Figura 11,
los sedimentos analizados, especialmente antes de entrar en la Depresin Central (muestras
CH08-01, -02, -04, -05 y -06, en la Figura 11 a, b, d, e y f), presentan signaturas de Sr-Nd muy
similares a aquellas de rocas volcnicas de volcanes activos geogrficamente cercanos a la
zona de estudio (36S y 37,5S, especialmente este ltimo). Conforme se avanza ro abajo
desde la Cordillera Principal (a, b, d, e) hacia la Depresin Central (f) y la Cordillera de la Costa
(g), la signatura de los sedimentos se vuelve ligeramente ms radiognica. Lo mismo se
observa en la muestra CH08-03 o "c" pero esta vez como anomala pues sta se tom a slo 7
km ro abajo de "b" (entre "e" y "f" hay aprox. 60 km de distancia). La muestra CH08-08 (h) fue
tomada ro abajo de la confluencia con el ro Itata, en el batolito grantico paleozoico.

87 86
Figura 11: Grfico Nd vs. Sr/ Sr. Los tringulos corresponden a rocas de volcanes activos cercanos
geogrficamente tomados de Futa y Stern (1988) y Hickey et al. (1986). En la simbologa se indica la
latitud a la que corresponde cada dato. Los puntos rojos son las muestras de sedimentos del ro uble
(indexados de la "a" a la "h" ordenados de este a oeste). Se han includo barras de error de dos veces el
error estndar externo.

Los resultados sugieren que la principal fuente sedimentaria son las rocas volcnicas de la
Cordillera Principal, con aportes subordinados desde el Batolito Grantico de la Cordillera.

52
7.2 Desequilibrios en la serie del U-238
Los valores de la Tabla 11 se usaron para calcular las razones de actividad (234U/238U),
( Th/234U) y (226Ra/230Th), que fueron graficadas en funcin de la distancia a lo largo del cauce
230

del ro uble desde su nacimiento (Figuras 12, 13 y 15).

(234U/238U)
La razn de actividad (234U/238U) se mantiene cerca de la unidad (condicin de equilibrio),
con algunas fluctuaciones. El 234U y el 238U tienen un comportamiento qumico idntico, porque
la diferencia relativa de masas es depreciable y casi no hay fraccionamiento isotpico. Sin
embargo, el fenmeno de "alpha-recoil" propicia la movilidad del 234U respecto al 238U desde los
minerales que contienen uranio hacia fases secundarias o el agua del ro. Aparentemente este
mecanismo no es eficiente en este caso.

234 238
Figura 12: Grfico ( U/ U) vs. distancia a lo largo del thalweg del ro uble (Itata para el punto en 155
km).

(230Th/238U)
La razn de actividad (230Th/238U) flucta en torno de uno (ver Figura 13), prxima al
equilibrio secular. Junto a esto, la razn U/Th presenta una ligera tendencia a disminuir a lo
largo del ro, exceptuando el punto anmalo CH08-03 (ver Figura 14). La razn U/Th presenta
un patrn casi idntico al de (230Th/232Th) (ver Figura 14), lo cual no es sorprendente puesto que
la proporcin entre estas dos cantidades es (230Th/238U) salvo constantes18, y esta ltima es una
cantidad que vara muy poco. La naturaleza radiognica del 230Th (viene solamente del 234U), y

18 238 232
En la naturaleza, el 99,8% del uranio es U, y prcticamente el 100% del torio es Th.

53
230 234
Figura 13: Grfico ( Th/ U) vs. distancia a lo largo del thalweg del ro uble (Itata para el punto en 155
km).

230 232
Figura 14: Grfico U/Th y ( Th/ Th) vs. distancia a lo largo el ro uble.
Los valores han sido divididos por su promedio. El tamao de los puntos azules
en el eje "y" equivale aproximadamente a cuatro veces el error estndar externo.

54
la cuasi-estabilidad del 232Th (tiene una vida media de ~14 Ga, similar a la edad del universo),
son la razn de que los cambios en 230Th/232Th estn fuertemente correlacionados con U/Th19.

En ambientes superficiales el torio es casi siempre insoluble e inmvil (no necesariamente


el uranio). La ligera pendiente negativa indica sin embargo, que el uranio se empobrece lenta-
pero sistemticamente respecto del torio en los sedimentos ro abajo. De manera consistente,
en las muestras de agua del ro uble la concentracin de U se mantiene por debajo de 1 ppt
(Tabla 8), y es del orden de 10 o ms veces la concentracin de Th, y entonces el uranio es
disuelto aunque en bajas concentraciones. La movilidad del uranio puede verse limitada, por
ejemplo a travs de la adsorcin sobre minerales como hidrxidos de hierro, gel de slice,
arcillas y micas. Esto ltimo parece ser el caso, dado que la razn U/Th es casi constante en los
sedimentos y la concentracin de U es muy baja en el agua.

Por completitud se aade el grfico (230Th/232Th) vs. (238U/232Th) (Figura 15).

230 232 238 232


Figura 15: Grfico ( Th/ Th) vs. ( U/ Th) en muestras de sedimentos del ro uble.
La barra de colores indica el cdigo de cada muestra (01=CH08-01, ...).

(226Ra/230Th)
La razn de actividad 226Ra-230Th en los sedimentos parte desde 1,00 (equilibrio secular) y
aumenta hasta alcanzar un valor mximo de 1,082 27 km ro abajo (ver Figura 16; obviar la
muestra anmala CH08-03 en d=33km). Unos 22 km ms adelante, baja hasta 1,045 y se
mantiene en ese nivel hasta la confluencia con el ro Itata, despus de la cual aumenta
ligeramente. Si bien estos desequilibrios son mayores que en el caso del 230Th-234U o el 234U-
238
U, cabe mencionar que siguen siendo pequeos en comparacin a aquellos medidos en
partculas en suspensin del ro Amazonas por Dossetto et al. (2006a) ((226Ra/230Th) = 0,86-
1,00), o en sedimentos del lecho del ro Ganges por Granet et al. (2007) ((230Th/238U) = 0,89-
1,45).

19
En el caso del punto CH08-03 la anomala se debe probablemente a distintas fuentes sedimentarias (granito en
vez de andesitas).

55
Por otro lado, es interesante ver la evolucin del Ba, dado que tiene un comportamiento
qumico muy similar al Ra, elemento del cual no se tiene datos de concentracin. La razn
Ba/Th en los sedimentos se mantiene aproximadamente constante a lo largo del ro uble,
dentro de las barras de error (ver Figura 17), con un aumento leve (a nulo) entre d=120 km y
d=150 km (el ltimo punto est ro abajo de la confluencia con el ro Itata y no entra en el
anlisis). En las muestras de agua del ro uble se observa un patrn similar (ver Figura 18,
todas las muestras salvo la segunda -CH08-03- y la ltima -CH08-08- son del ro uble), pero
en el primer tramo del ro el Ba dismuye y luego se estabiliza y aumenta entre d=70 km y d=120
km.

226 230
Figura 16: Grfico ( Ra/ Th) vs. distancia a lo largo del thalweg del ro uble (Itata para el punto en
155 km)

Tanto el 230Th como el 226Ra son istopos puramente radiognicos, es decir, son producto
exclusivo de la desintegracin radioactiva del 238U. Entonces, suponiendo que las fuentes
sedimentarias (rocas frescas) cumplen la condicin de equilibrio secular20, las variaciones en la
razn de actividad 226Ra-230Th en los sedimentos se deben slo a mecanismos que fraccionan
Ra respecto de Th. Los mecanismos ms probables, dado el ambiente geolgico, son los
siguientes.
1. Aportes de sedimentos con una composicin isotpica diferente.
2. Transferencia diferencial de Ra y Th desde los minerales que contienen U (fuente de 226Ra
y 230Th) hacia el agua. Esto puede ocurrir por disolucin o por "alpha-recoil".
3. Fijacin de Ra y/o Th sobre la superficie de minerales (por adsorcin) desde solucin.
Como se explica en 2.3.4, el fenmeno de "alpha-recoil" puede fraccionar 226Ra respecto de
su padre de forma eficiente en este caso.

20
Recordar que el equilibrio secular se alcanza si el sistema ha estado cerrado durante al menos 5 a 6 vidas medias
230
(en este caso del Th, es decir, 450 ka).

56
4. "Fraccionamiento" de fases minerales que contienen Th y/o Ra21.

Ms all del motivo de los cambios en (226Ra/230Th) de los sedimentos a lo largo del cauce,
interesa acotar el tiempo de residencia de stos entre la entrada a la Depresin Central (CH08-
04 a -05) y la vertiente oriental de la Cordillera de la Costa (CH08-09) unos 70 km ro abajo.
Obviando la muestra anmala (la tercera), solamente entre entre los puntos cuarto y quinto hay
una disminucin de (226Ra/230Th)22. Esto se puede explicar por el efecto conjunto de la
desintegracin radioactiva del 226Ra, la deposicin de partculas en suspensin al entrar al valle
central y los aportes de sedimentos desde las Formaciones Mininco y La Montaa23. Si se
debiera a la desintegracin radioactiva de 226Ra exclusivamente (en el caso extremo de sistema
cerrado), el tiempo de residencia de los sedimentos en este tramo de poco ms de 20 km
tendra que ser de 1370 aos aprox. Este resultado es solamente hipottico. Para el clculo se
us la frmula aproximada dada al final de la seccin 2.1:

A2 A2 t 226
Ra 226
Ra t
1 1 e 2 1 1 e 226
A1 A1 0
230
Th 230
Th 0
R 1 226
Ra
t ln 0 / 226 donde R
R 1
230
Th

Figura 17: Grfico Ba/Th vs. distancia a lo largo del thalweg del ro uble en sedimentos.

21
Suponiendo que hay al menos dos minerales portadores de estos elementos (por ejemplo un mineral magmtico
primario y uno secundario que adsorbe), stos probablemente tendrn proporciones de Ra/Th distintas entre s. Si
se modifica la proporcin en la que se encuentran estos minerales en el sedimento de forma significativa, ello
fraccionar Ra respecto de Th.
22
Slo en tramos dondo hay una disminucin (cambio negativo) de esta cantidad, puede atriburse ese cambio a
desintegracin radioactiva. En el caso opuesto, los aportes de Ra predominan sobre la desintegracin.
23
Tal como se interpret la evolucin qumica e isotpica de Sr y Nd.

57
Si se asumiera que no hay aportes (slo prdida) de Ra respecto de Th en este tramo, el
tiempo que le habra tomado a los sedimentos en atravesar estos 20 km es necesariamente
menor a 1370 aos (tiempo para sistema cerrado). El contenido de Ba de los sedimentos se
mantiene constante en este tramo (ver Figura 17), sugiriendo que el Ra tampoco debera variar.
Sin embargo, se necesitaran datos qumicos ms precisos para testear esta hiptesis al nivel
de precisin de los desequilibrios.

Para explicar que (226Ra/230Th) no est ms cerca de 1,0 al llegar a la Cordillera de la Costa
que a la entrada del Valle Central, es necesario que haya aportes qumicos sistemticos a lo
largo del valle con Ra/Th > 1. Por otro lado, en la muestra CH08-08, tomada ro abajo de la
confluencia con el ro Itata y en medio del Batolito Grantico de la Costa, la razn (226Ra/230Th)
es ligeramente superior que en CH08-09, probablemente producto de los aportes desde la
mitad sur de la cuenca del Itata y/o desde el Batolito de la Costa. Una variante sencilla de la
ecuacin de desintegracin radioactiva para el 226Ra, sera incluir un trmino constante a la
derecha que representa una tasa de aporte de 226Ra.
d 226Ra
= 226 226 Ra + 230 230 Th + f 226 si 230
Th cte= 230 Th0 entonces
dt
226
( )
Ra (t ) 226 Ra 0 230 230 Th0 f 226 e 226t + 230 230 Th0 + f 226
226
Ra 226 226 Ra (t ) 226 226 Ra (t ) 226
Ra f 226 226t
e f 226
(t ) = 1 +1+

230
Th 230 230Th(t ) 230 230Th0

230
Th 0 230 230Th0
230 230Th0

El valor de (226Ra/230Th)=1,045 en CH08-05 y CH08-09, puede ser interpretado como la


constante que se suma a la exponencial por la derecha, i.e. f226/(230230Th0) ~0,045. En este
caso, f226/230Th0 ~ 4,1210-7 aos-1. Como (226Ra/230Th)0 prcticamente coincide con el valor final
en este modelo, se trata de una especie de estado estacionario.

Si se usa el valor de f226/230Th0 recin calculado para el tramo entre d=50 km (CH08-04) y
d=140 km (CH08-09), resolviendo la ecuacin, resulta
226
Ra 226
Ra t
(t ) F e 226 + F

230
Th

230
Th 0

f 226 226
Ra
con F = 1+ = 1,045 y R =
230 230 Th 0
230
Th
Pero esta vez R0 = 1,082, R(t) = 1,046 ~ F. En este caso el tiempo es muy grande, que en
la prctica significa que el tiempo de residencia calculado sera al menos 5 o 6 veces la vida
media del 226Ra, es decir, ~ 8000-10000 aos (8-10 ka) en un tramo de 90 km. El problema con
esta estimacin es que a slo 20 km del punto inicial hay un punto intermedio donde R = 1,045
= F, y para el cual se acaba de dar otra interpretacin en la que F=1. Luego, el modelo de flujo
constante no es satisfactorio. Sin embargo, es claro que debe haber un flujo neto entrante (no
uniforme) de 226Ra (f226) para explicar los datos de (226Ra/230Th) en el ltimo tramo, a menos que
la transferencia sedimentaria sea tan rpida que el sistema no alcance a equilibrarse en 90 km
(tiempo despreciable respecto de 8-10 ka). En realidad siempre existe una distribucin de f226 a
lo largo del cauce (en funcin del tiempo) que permite mantener constante el valor de
(226Ra/230Th) en todo momento. Para obtenerla, slo hace falta igualar d(226Ra/230Th)/dt a cero
en la ecuacin de arriba y resolver, aunque sera muy raro que ocurra sto. Consideremos un
escenario intermedio, donde el flujo entrante de 226Ra es lineal en el tiempo, F vale 1,0 en t=0

58
(flujo nulo al principio) y F* al final.

226 Ra
d 230

226
230 Th0 226 Ra + F (t ) si F (t ) = 1 + ( F * 1)
t
entonces
dt 230 230 Th0 T
Ra
226
226 226 Ra 226 Ra 1 a 226t 1 a
(t ) e + a t + donde
230
Th 230 Th0 Th 0 226
230 230
226
F * 1
a= viene de la e.d .
T226

Dado un valor de F*, se puede buscar el tiempo de residencia t=T imponiendo


226
Ra
(T )= 1,045

230
Th

Incluso en este modelo simple de dos parmetros, hay infinitas combinaciones de F* y T


que cumplen la condicin de arriba.

Por completitud, se aaden los siguientes grficos: Ba normalizado respecto de varios


elementos en muestras de agua vs. distancia a lo largo del ro uble (Figura 18), (226Ra/232Th)
vs. (238U/232Th) (Figura 19) y (226Ra/238U) vs. (230Th/238U) (Figura 20).

Figura 18: Grfico Ba/Z vs. distancia a lo largo del thalweg del ro uble en muestras de agua.
Z = K, Si, Na, Ca, Mg y Fe. El Ba tiene un comportamiento geoqumico similar al Ra.

59
226 232 238 232
Figura 19: Grfico ( Ra/ Th) vs. ( U/ Th) en sedimentos del ro uble
La barra de colores indica el cdigo de cada muestra (01=CH08-01, ...).

226 238 230 238


Figura 20: Grfico ( Ra/ U) vs. ( Th/ U) en sedimentos del ro uble.
La barra de colores indica el cdigo de cada muestra (01=CH08-01, ...).

60
CAPTULO 8: Interpretaciones y conclusiones

8.1 Qumica e istopos de Sr y Nd

8.1.1 Orgen de los sedimentos

La composicin qumica e isotpica de Sr y Nd de los sedimentos, sugiere que la principal


fuente sedimentaria son las rocas volcnicas de las Formaciones Cola de Zorro, Trapa-Trapa y
Cura-Malln, habiendo tambin aportes subordinados del batolito grantico.

8.1.2 Evolucin de los sedimentos

Las variaciones de este a oeste (ro abajo) en la composicin qumica e isotpica de Sr y


Nd del sedimento muestreado (< 2mm), pueden ser explicadas como una combinacin de los
siguientes procesos/fenmenos:
A. Variabilidad en la composicin de las fuentes del sedimento. Las distintas fuentes se
mezclan y dan origen a composiciones intermedias.
B. Seleccin de minerales durante el transporte (deposicin y transporte selectivos).
C. La disolucin de elementos desde los minerales del sedimento (asociada a la alteracin
de minerales por ejemplo).
D. La eventual precipitacin de nuevas fases minerales desde el agua.

Enseguida se interpreta en dos tramos del ro: Cordillera Principal y Valle Central

Cordillera Principal
Las muestras de sedimentos muestran una gran similitud con las rocas volcnicas que
conforman el macizo montaoso, respecto de su composicin de elementos mayores y de
istopos de Sr y Nd. El carcter anmalo de las muestra CH08-03 y CH08-08, y la evolucin
qumica e isotpica del resto de las muestras en este tramo, pueden ser explicadas en gran
medida por el aporte de sedimento proveniente de la erosin del batolito y la disolucin de
elementos desde el sedimento.

Depresin Central y la Cordillera de la Costa


Entre los puntos CH08-05 y CH08-06 se producen cambios qumicos significativos, por
ejemplo la disminucin de Fe, Mg, Ca, Mn, Sc, V y Co, y el aumento de K, Rb y Ba.
Geogrficamente esto corresponde a la entrada a la Depresin Central, justamente donde el
cambio de pendiente es ms violento. Mientras la pendiente del ro disminuye, tambin lo hace
la velocidad promedio del torrente, depositndose progresivamente las partculas en
suspensin, eventualmente aumentando la proporcin de arcillas. Esto podra modificar la
composicin qumica del sedimento. No puede descartarse que los aportes sedimentarios
desde las Formaciones Sedimentarias Mininco y La Montaa (en su mayora depsitos de
antiguos glaciares; son sedimentos de orgen volcnico), sean tambin responsables del
cambio qumico de los sedimentos del lecho del ro. Esos aportes constaran de materiales ms
antiguos y posiblemente ms maduros (en el sentido sedimentario), que los sedimentos que

61
vienen directamente de la cordillera, por lo tanto ms ricos en cuarzo y arcillas, y ms pobres en
los elementos propios de los minerales ferromagnesianos. Por ltimo, la disolucin de
elementos influye de forma directa en la evolcin qumica de los sedimentos.
La muestra CH08-09 fue tomada ro abajo de la confluencia con el ro Itata recibiendo
aportes de sedimento desde la mitad sur de la hoya hidrogrfica del ro Itata y desde el Batolito
de la Costa, ello explica su naturaleza anmala respecto de los otros datos.

En sntesis, la evolucin de los sedimentos se interpreta como el resultado de los siguientes


fenmenos: aportes sedimentarios minoritarios desde el batolito grantico en la cordillera,
deposicin progresiva de partculas en suspensin al entrar en el valle central, disolucin
diferenciada de elementos desde los sedimentos (especialmente en el Valle Central) y aportes
desde las formaciones sedimentarias en la Depresin Central.

El carcter acotado de la evolucin de la composicin qumica, mineralgica e isotpica de


los sedimentos, as como el ambiente geolgico y la longitud del ro uble (160 km), son
argumentos a favor de la hiptesis de tiempos de residencia relativamente cortos24 para los
sedimentos de este ro.

8.2 Desequilibrios en la serie del U-238

(234U/238U)
Los anlisis de sedimentos indican que no hubo fraccionamiento isotpico de U.

(230Th/238U)
Los desequilibrios son pequeos (casi dentro de las barras de error). Tanto el Th como el U
parecen ser poco mviles (especialmente el primero). No se pudo extraer informacin temporal
de los datos.

(226Ra/230Th)
La desviacin del equilibrio secular ((226Ra/230Th) = 1) nunca supera el valor de 0,085, que
es mucho menor que las desviaciones del equilibrio medidas en los ros Amazonas y Ganges, lo
cual apoya la hiptesis de tiempos de residencia sustancialmente ms cortos.

Slo el tramo entre los puntos cuarto (CH08-04) y quinto (CH08-05), de 20 km de longitud,
permiti extraer alguna informacin temporal, estimndose un tiempo de residencia mximo25
de 1400 aos. Cabe mencionar que esta estimacin no es en ningn caso una cuantificacin
certera de un tiempo de residencia sedimentario, sino slo una primera aproximacin de
aquello. Dada la complejidad del fenmeno en estudio, sera necesario contar con ms datos y
acotar cuantitativamente varias variables (flujos qumicos por ejemplo) para poder modelar la
serie del 238U en los sedimentos de forma ms satisfactoria.

Sin embargo, estos resultados preliminares adquieren su verdadera importancia en su


contexto global: el proyecto de investigacin sobre transferencia sedimentaria en varias
cuencas de los Andes Centrales. Es de particular inters la comparacin con los resultados a
obtenerse para el ro Choapa26 (y otros ros) cuyo clima es marcadamente ms seco respecto

24
Por ejemplo respecto de los ros Ganges y Amazonas, con tiempos de residencia del orden de 100 ka.
25
Bajo el supuesto de que no hay aportes (slo prdida) de Ra en ese tramo. Se necesitara de datos qumicos ms
precisos para testear este supuesto.
26
Muestras en anlisis por el grupo de Estrasburgo (LHyGeS).

62
de aquel del uble, con dinmicas de transporte sedimentario a priori diferentes. Intuitivamente,
menos lluvias y menos escorrenta superficial podran significar una tasa de transporte ms baja
y tiempos de residencia mayores en promedio.

Esta metodologa mostr ser limitada en el caso del ro uble, estando sujetos los
resultados finales (tasas de transferencia sedimentaria) a ms informacin. Su validez y
aplicabilidad para ros cortos podra depender a priori de la intensidad de la meteorizacin en la
cuenca, y por ende de las condiciones climticas. Los anlisis de sedimentos de cuencas
andinas ms ridas contemplados en el proyecto SESPEED permitirn eventualmente conclur
al respecto.

63
Referencias

Aguirre, L., Herv, F., Godoy, E., 1972. Distibution of metamorphic facies in Chile. An outline.
Kristalinikum, vol. 9, 7-19.

Anderson, R.S., Repka, J.L., Dick, G.S., 1996. Explicit treatment of inheritance in dating depositional
surfaces using in situ 10Be and 26Al. Geology, vol. 24, 4751.

Arvalo, R., McDonough, W.F., 2010. Chemical variations and regional diversity observed in MORB.
Chem. Geol., vol. 271, 70-85.

Bogaart, P.W., Van Balen, R.T., 2000. Numerical modelling of the response of alluvial rivers to Quaternary
climate change. Global and Planetary Change, vol. 27, 124141.

Bogaart, P.W., Tucker, G.E., de Vries, J.J., 2002. Channel network morphology and sediment dynamics
under alternating periglacial and temperate regimes: a numerical simulation study. Geomorphology, vol.
54, 257-277.

Bourdon, B., Turner, S., Henderson, G.M., Lundstrom, C.C., 2003. Introduction to U-series Geochemistry.
In: Bourdon, B., Henderson, G.M., Lundstrom, C.C., Turner, S.P. (Eds.), Uranium-series Geochemistry,
Reviews in Mineralogy and Geochemistry, vol. 52. Geochemical Society-Mineralogical Society of America,
Washington, pp. 1-21.

Braun, J.-J., Pagel, M., Herbillon, A., Rosin, C., 1993. Mobilization and redistribution of REEs and thorium
in a syenitic lateritic profile: a mass balance study. Geochim Cosmochim Acta, vol. 57, 4419-4434

Brown, E., Bourles, D., Burchfiel, B., Deng, Q., Li, J., Molnar, P., Raisbeck, G., Yiou, F., 1998. Estimation
of slip rates in the southern Tien Shan using cosmic ray exposure dates of abandoned alluvial fans. Geol.
Soc. Am. Bull., vol. 110, 377386.

Carretier, S., Regard, V., Soual, C., 2009. Theoretical cosmogenic nuclide concentration in river bed load
clasts: Does it depend on clast size?. Quaternary Geochronology, vol 4, 108-123

Cecioni, G., Westermann, G., 1968. The Triassic/Jurassic marine transition of Coastal Central Chile.
Pacific Geology, No. 1, 41-75.

Chabaux, F., Bourdon, B., Riotte, J., 2008. U-Series geochemistry in weathering profiles, river waters and
lakes. In: Krishnaswami, S. y Kirk Cochran, J., U-Th Series nuclides in Aquatic Systems. Radioactivity in
the Environment, vol. 13. Elsevier, pp. 49-104.

Chabaux, F., Riotte, J., Dequincey, O., 2003. U-Th-Ra fractionation during weathering and river transport.
In: Bourdon, B., Henderson, G.M., Lundstrom, C.C., Turner, S.P. (Eds.), Uranium-series Geochemistry,
Reviews in Mineralogy and Geochemistry, vol. 52. Geochemical Society-Mineralogical Society of America,
Washington, pp. 553-576.

Crave, A., Davy P., 2001. A stochastic precipiton model for simulating erosion/sedimentation dynamics.
Computers & Geosciences, vol. 27, 815827.

Dickson, B. L., 1985. Radium isotope in saline seepages, southwestern Yilgarn, western Australia.
Geochim.Cosmochim. Acta, vol. 49, 361366.

Dosseto, A., Bourdon, B., Gaillardet, J., Allgre, C.J., Filizola, N., 2006a. Time scale and conditions of
weathering under tropical climate: study of the Amazon basin with U-series. Geochim.Cosmochim. Acta,
vol. 70, 71-89.

64
Dosseto, A., Turner, S.P., Douglas, G.B., 2006b. Uranium-series isotopes in colloids and suspended
sediments: timescale for sediment production and transport in the Murray-Darling River system. Earth
Planet. Sci. Lett., vol. 246, 418-431.

Dosseto, A., Bourdon, B., Gaillardet, J., Maurice-Bourgoin, L., Allgre, C.J., 2006c. Weathering and
transport of sediments in the Bolivian Andes: Time constraints from uranium-series isotopes. Earth Planet.
Sci. Lett., vol. 248, 759-771.

Escobar, F., Guzmn, R., Vieira, C., 1977. Avance geolgico de las Hojas Rancagua-Curic, Talca-
Linares, Chanco, Concepcin-Chilln. Inst. Invest. Geol. (indito), 56 pp. Santiago.

Ferraris, F., 1981. Hoja Los Angeles-Angol. 1 mapa escala 1:250.000. Carta Geolgica de Chile N5,
Preliminar. Instituto de Investigaciones Geolgicas.

Fredrikson, J.K., Zachara, J.M., Kennedy, D.W., Duff, M.C., Gorby, Y.A., Li, S.M.W., Krupka, K.M., 2000.
Reduction of U(VI) in goethite (-FeOOH) suspensions by a dissimilatory metal-reducing bacterium.
Geochim. Cosmochim. Acta, vol. 64, 3085-3098.

Fuenzalida, P.H., 1971. Climatologa de Chile. Universidad de Chile, Departamento de Geofsica,


Publicacin interna de la Seccin de Meteorologa, 73 pp.

Futa, K., Stern, C., R., 1988, Sr and Nd isotopic and trace element compositions of Quaternary volcanic
centers of the southern Andes, Earth Planet. Sci. Lett., vol. 88, 253-262.

Gajardo, A., 1981. Hoja Concepcin-Chilln. 1 mapa escala 1:250.000. Carta Geolgica de Chile N4,
Preliminar. Instituto de Investigaciones Geolgicas.

Gardeweg, M., 1980. Geologa del rea del Nevado de Longav, VII Regin del Maule, Chile. Memoria de
Ttulo. Univ. Chile, Depto. Geol., Santiago, 247pp.

Gascoyne, M., 1992. Geochemistry of the actinides and their daughters. In: Uranium-series disequilibrium:
Application to Earth, Marine, and Environmental Sciences. Ivanovich M, Harmon RS (eds), Oxford
Sciences Publications, Oxford, pp. 34-61.

Gonzlez-Bonorino, F., 1970. Series metamrficas del basamento cristalino de la Cordillera de la Costa,
Chile Central. Univ. Chile, Depto. Geol., Public., No. 37, 67 pp.

Granet, M., Chabaux, F., Stille, P., France-Lanord, C., Pelt, E., 2007. Time-scales of sedimentary transfer
and weathering processes from U-series nuclides: Clues from the Himalayan rivers. Earth Planet. Sci.
Lett., vol. 261, 389406.

Granger, D., Kircher, J., Finkel, R., 1996. Spatially averaged long-term erosion rates measured from in
situ-produced cosmogenic nuclides in alluvial sediment. J. Geol., vol. 104, 249257.

Gueniot, B., Munier-Lamy, C., Berthelin, J., 1988a. Geochemical behavior of uranium in soils, part I.
Influence of pedogenetic processes on the distribution of uranium in aerated soils. J. Geochem. Explor.,
vol. 31, 21-37.

Gueniot, B., Munier-Lamy, C., Berthelin, J., 1988b. Geochemical behavior of uranium in soils, part II.
Distribution of uranium in hydromorphic soils and soil sequences. Application for surficial prospecting. J.
Geochem. Explor., vol. 31, 39-55.

Hancock, G., Anderson, R., Chadwick, O., Finkel, R., 1999. Dating fluvial terraces with 10Be and 26Al
profiles: application to the Wind River, Wyoming. Geomorphology, vol. 27, 4160.

65
Herv, F., Munizaga, F., 1978. Evidencias geocronolgicas de un magmatismo intrusivo trisico superior-
jursico en la Cordillera de la Costa de Chile entre los 3530' y 3630' Lat. S. In Congr. Geol. Argent.,
Neuqun, 1978, No. 7, Actas, vol. 2, pp. 43-52.

Herv, F., Munizaga, F., Mantovani, M., et al., 1976. Edades Rb-Sr neopaleozoicas del basamento
cristalino de la Cordillera de Nahuelbuta. In Congr. Geol. Chileno, Santiago, No. 1, Actas, vol. 2, F19-F26.
Hickey R., Frey, F., Gerlach, D., Lopez-Escobar, L., 1986, Multiple sources for basaltic arc rocks from the
Southern Volcanic Zone of the Andes (34-41S): Trace element and isotopic evidence for contributions
from subducted oceanic crust, mantle and continental crust, Journal of Geophysical Research, vol. 91, n
B6, 5963-5983.

Huisink, M., 1999. Changing fluvial styles in response to climate change. Examples from the Maas and
Vecht during the Weichselian PleniFand Lateglacial. Thesis Vrije Universiteit, Amsterdam, 127pp.

Informe DGA (Direccin General de Aguas), 2004. Diagnostico y clasificacin de los cursos y cuerpos de
agua segn objetivos de calidad, cuenca del ro Itata.

Jarvis, K. E.; A. L. Gray; and R. S. Houk, 1992. Handbook of Inductively Coupled Plasma Mass
Spectrometry. Chapman and Hall: New York.

Kigoshi, K., 1971. Alpha-recoil thorium-234: dissolution into water and the uranium-234/uranium-238
disequilibrium in nature. Science, vol. 173, 47-48.

Kirchner, J., Riebe, C., Ferrier, K., Finkel, R., 2006. Cosmogenic nuclide methods for measuring long-term
rates of physical erosion and chemical weathering. J. Geochem. Explor., vol. 88, 296299.
222 40 39 37
Krishnaswami, S., and D. E. Seidemann, 1988. Comparative study of Rn, Ar, Ar and Ar leakage
from rocks and minerals: Implications for the role of nanopores in gas transport through natural silicates.
Geochim.Cosmochim. Acta, vol. 52, 655658.

Lagno, L., 1979. Geologa del cuadrngulo Pocillas, VII Regin. Taller de Ttulo, Univ. Chile, Depto. Geol.
(indito). Santiago.

Langmuir, D., and D. Melchior., 1985. The geochemistry of Ca, Sr, Ba and Ra sulfates in some brines from
the Palo Duro Basin Texas. Geochim Cosmochim Acta, vol. 49, 24232432.

Lenhart, J.J., Cabanis, S.E., McCarthy, P., Honeymann, B.D., 2000. Uranium (VI) complexation with
citric,humic and fulvic acids. Radiochim. Acta, vol. 58, 5455-5463.

Lovley, D., Phillips, E.J.P., 1992. Reduction of uranium by Desulfovibrio desulfuricans. Appl. Environ.
Microbio., vol. 58, 850-856.

Lovley, D., Phillips, E.J.P., Gorby Y.A., Landa E.R., 1991. Microbial reduction of uranium. Nature, vol. 350,
413-416.

Middleburg, J.J., van der Weijden, C.H., Woittez, J.R.W., 1988. Chemical processes affecting the mobility
of major, minor and trace elements during weathering of granite rocks. Chem. Geol., vol. 68, 253-273.

Montavon, G., Mansel, A., Seibert, A., Keller, H., Kratz, J.V., Trautmann, N., 2000. Complexation studies of
2+
UO2 with humic acid at low metal ion concentrations by indirect speciation methods. Radiochim. Acta,
vol. 88, 17-24.

Moreira-Nordemann, L.M., 1977. Etude de la vitesse d'altration des roches au moyen de l'uranium utilis
comme traceur naturel. Application deux bassins du Nord Est du Brsil. Thse, Universit Paris VI,
France.

66
234 238
Moreira-Nordemann, L.M., 1980. Use of U/ U disequilibrium in measuring chemical weathering rate of
rocks. Geochim. Cosmochim. Acta, vol. 44, 103-108.

Moreira-Nordemann, L.M., Sieffermann, G., 1979. Distribution of uranium in soil profiles of Bahia state,
Brazil. Soil Science, vol. 127, 275-280.

Moreno-Roa, H., Herv, F., Godoy, E., et al. 1976. Reconocimiento geolgico del rea de Pocillas
(Cauquenes, VII Regin). Univ. Chile, Depto. Geol. (indito), Santiago, 9pp.

Muoz, J., Niemeyer, H., 1983. Hoja Laguna del Maule. 1 mapa escala 1:250.000. Carta Geolgica de
Chile N64. Servicio Nacional de Geologa y Minera.

Olesik, John W. 1996. Fundamental Research in ICP-OES and ICPMS. Analytical Chemistry, vol. 68,
469A-474A.

Plater, A.J., Dugdale, R.E., Ivanovich, M., 1988. The application of uranium series disequilibrium concepts
to sediment yield determination. Earth Surf. Process. Landf., vol. 13, 171-182.

Plater, A.J., Ivanovich, M., Dugdale, R.E., 1992. Uranium series disequilibrium in river sediments and
waters: the significance of anomalous activity ratios. Appl. Geochem., vol. 7, 101-110.

Plater, A.J., Dugdale, R.E., Ivanovich, M., 1994. Sediment yield determination using uranium-series
radionuclides: the case of the Wash and Finland drainage basin, eastern England. Geomorphology, vol.
11, 41-56.

Regard, V., Bellier, O., Braucher, R., Gasse, F., Bourles, D., Mercier, J., Thomas, J., Abbassi, M.,
Shabanian, E., Soleymani, S., 2006. 10Be dating of alluvial deposits from Southeastern Iran (the Hormoz
Strait area). Palaeogeogr. Palaeoclimatol. Palaeoecol., vol. 242, 3653.

Repka, J.L., Anderson, R.S., Finkel, R.C., 1997. Cosmogenic dating of fluvial terraces, Fremont River,
Utah. Earth Planet. Sci. Lett., vol. 152, 5973.

Riebe, C., Kirchner, J., Granger, D., Finkel, R., 2001a. Minimal climatic control on erosion rates in the
Sierra Nevada, California. Geology, vol. 29, 447450.

Riebe, C., Kirchner, J., Granger, D., Finkel, R., 2001b. Strong tectonic and weak climatic control of long-
term chemical weathering rates. Geology, vol. 29, 511514.
234 238
Riotte, J., Chabaux, F., 1999. ( U/ U) activity ratios in freshwaters as tracers of hydrological processes:
the Strengbach watershed (Vosges, France). Geochim. Cosmochim. Acta, vol. 63, 1263-1275.

Riquelme, R., Martinod, J., Hrail, G., Darrozes, J., Charrier, R., 2003. A geomorphological approach to
determining the Neogene to Recent tectonic deformation in the Coastal Cordillera of northern Chile
(Atacama). Tectonophysics, vol. 361, 255 275.
226 228 224
Rihs, S., and M. Condomines, 2002. An improved method for Ra isotope ( Ra, Ra, Ra)
measurements by gamma spectrometry in natural waters: Application to CO2-rich thermal waters from the
French Massif Central. Chem. Geol., vol. 182, 409421.

Ritz, J.-F., Brown, E.T., Bourle` s, D.L., Philip, H., Schlupp, A., Raisbeck, G.M., Yiou, F., Enkhtuvshin, B.,
1995. Slip rates along active faults estimated with cosmic-rayexposure dates: application to the Bogd fault,
Gobi-Altai, Mongolia. Geology, vol. 23, 10191022.

Ritz, J.-F., Vassallo, R., Braucher, R., Brown, E., Carretier, S., Bourles, D., 2006. Using in situ produced
10
Be to quantify active tectonics in the Gurvan Bogd mountain range (Gobi-Altay, Mongolia). In: Siame, L.,

67
Bourles, D., Brown, E. (Eds.), In Situ Produced Cosmogenic Nuclides and Quantification of Geological
Processes. Geological Society of America Special Paper 415, 87110.

Roddaz, M., Baby, P., Brusset, S., Hermoza, W., Darrozes, J., 2005. Forebulge dynamics and
environmental control in western Amazonia: the case study of the arch of Iquitos (Peru). Tectonophysics,
vol. 399, 87108.

Schaller, M., von Blanckenburg, F., Hovius, N., Kubik, P., 2001. Large-scale erosion rates from in situ-
produced cosmogenic nuclides in European river sediments. Earth Planet. Sci. Lett., vol. 188, 441458.

Schaller, M., von Blanckenburg, F., Veldkamp, A., Tebbens, L., Hovius, N., Kubik, P., 2002. A 30000 yr
10
record of erosion rates from cosmogenic Be in Middle European river terraces. Earth Planet. Sci. Lett.,
vol. 204, 307320.

Siame, L., Bellier, O., Braucher, R., Sebrier, M., Cushing, M., Bourles, D., Hamelin, B., Baroux, E., de
Voogd, B., Raisbeck, G., Yiou, F., 2004. Local erosion rates versus active tectonics: cosmic ray exposure
modelling in Provence (south-east France). Earth Planet. Sci. Lett., vol. 220, 345-364.

Spear, J.R., Figueroa, L.A., Honeyman, B.D., 2000. Modeling the removal of uranium (VI) under variable
sulfate concentrations by sulfate-reducing bacteria. Appl Envir Microbiology, vol. 66, 3711-3721.

Sturchio, N. C., J. K. Bohlke, and F. J. Markun, 1993. Radium isotope geochemistry of thermal waters,
Yellowstone National Park, Wyoming, USA. Geochim. Cosmochim. Acta, vol. 57, 12031214.

Sun, H., Semkow, T.M., 1998. Mobilization of thorium, radium and radon nuclides in ground water by
successive alpha-recoils. J Hydrol., vol. 205, 126-136.

Suzuki, Y., Banfield, J.F., 1999. Geomicrobiology of uranium. Rev. Min. Geochem., vol. 38, 393-432.

Swarzenski, P.W., D. Porcelli, P. S. Andersson, and J. M. Smoak, 2003. The behaviour of U- and Th series
nuclides in groundwater. In: U-series Geochemistry (Eds B. Bourdon, G. M. Henderson, C. C. Lundstrom,
and S. P. Turner). Reviews in Mineralogy and Geochemistry, vol. 52, 317361.

Tassara, A., Swain, C., Hackney, R., Kirby, F., 2007. Elastic thickness structure of South America
estimated using wavelets and satellite-derived gravity data. Earth Planet. Sci. Lett., vol. 253, 17-36.

Tieh, T., Ledger, E., Rowe, M., 1980. Release of uranium from granitic rocks during in situ weathering and
initial erosion (Central Texas). Chem. Geol., vol. 29, 227-248.

van-der Woerd, J., Tapponnier, P., Ryerson, F., Meriaux, A., Meyer, B., Gaudemer, Y., Finkel, R., Caffee,
M., Zhao, G., Xu, Z., 2002. Uniform postglacial slip-rate along the central 600 km of the Kunlun Fault
(Tibet), from Al-26, Be-10, and C-14 dating of riser offsets, and climatic origin of the regional morphology.
Geophys. J. Int., vol. 148, 356388.

Verhaar, P., Biron, P., Ferguson, R., Hoey, T., 2008. A modified morphodynamic model for investigating the
response of rivers to short-term climate change. Geomorphology, vol. 101, 674682.

Vigier, N., Bourdon, B., Turner, S., Allgre, C.J., 2001. Erosion timescales derived from U-decay series
measurements in rivers. Earth Planet. Sci. Lett., vol. 193, 549563.

von Blanckenburg, F., 2005. The control mechanisms of erosion and weathering at basin scale from
cosmogenic nuclides in river sediment. Earth Planet. Sci. Lett., vol. 237, 462479.

von Blanckenburg, F., Hewawasam, T., Kubik, P., 2004. Cosmogenic nuclide evidence for low weathering
and denudation in the wet, tropical highlands of Sri Lanka. J. Geophys. Res., vol. 109, F03008, 1-22.

68
Anexo A: Grficos de elementos mayores y traza

Considerando que el Al y el Th son en general menos mviles que el resto de los


elementos, los datos qumicos de las muestras de sedimento de lecho fueron normalizados
respecto de estos dos elementos. Enseguida se dividi los valores normalizados por su
promedio para tener variaciones en torno de la unidad. Los puntos verdes corresponden a
muestras tomadas en el lecho del ro uble, y los rojos en otro lugar de la cuenca del ro Itata
(CH08-03 se tom en un pequeo tributario sobre roca grantica, y CH08-09 se tom ro abajo
de la confluencia con el ro Itata). Si bien el la proporcin Al2O3/Th se mantiene en un valor
estable a ligeramente decreciente ro abajo (dentro de las barras de error) para los puntos
verdes, los puntos rojos muestran valores fuera de la tendencia. Es importante tener en cuenta
esta observacin, al comparar los valores normalizados respecto del Al con aquellos
normalizados respecto del Th.

69
70
71
72
73
74
75
76
77
78
79
80
81
82
Anexo B: Resinas cromatogrficas

En este apartado se explica el funcionamiento de las resinas cromatogrficas utilizadas.


stas pertenecen a dos categoras: (1) resinas intercambiadoras de iones y (2) resinas de
extraccin (ver Tabla 15).

Categora Nombre Fabricante granulometria grupo funcional Uso


comercial (micrones)
amonio cuaternario Separacion y
intercambiadoras

AG1x8 Biorad 200-400 + purificacion de U y


(-NR3 )
Th
de iones
Resinas

TEMEX Separacion de Ra
Touzart - y de Ba (eliminar
Biorad 200-400 sulfonato (-SO3 )
Matignon elementos
50Wx8 mayores)

Separacion Ra de
Sr - Spec EICHROM 50-100 No tiene Ba. Separar Sr de
Resinas de extraccin

otros elementos
Separacion de
27
LREE (entre
TRU - Spec EICHROM - No tiene ellas Nd) de
elementos
mayores.
Purificacion de
Ln - Spec EICHROM - No tiene LREE y
separacion de Nd

Tabla 15: Resinas usadas para la separacin qumica

B.1 Resinas intercambiadoras de iones


Un intercambiador de iones es un slido capaz de intercambiar parte de sus iones con los
iones provenientes de una solucin acuosa. Para acelerar los intercambios, el slido debe estar
finamente dividido, a fin de presentar la mxima superficie de contacto con la solucin (granos
finos). En la naturaleza existen minerales intercambiadores de iones, por ejemplo las zeolitas.
Sin embargo, en la practica seria imposible utilizarlos para cromatografa, debido a numerosos
inconvenientes: complejidad de los fenmenos de intercambio, frecuente falta de
reproductibilidad, inestabilidad frente a los cidos y bases fuertes. Por eso se usan resinas
sintticas.

Las resinas son polmeros orgnicos sintticos muy resistentes a los cidos y las bases, los
oxidantes y los reductores. En general se trata de poliestireno, sobre el cual estn ligados
grupos funcionales ionizados o ionizables que le confieren la propiedad de intercambiador de
iones. Segn la carga del grupo funcional ionizado ligado a la resina, se diferencian dos clases
de resinas (ver Tabla 16).

27
LREE = Light Rare Earth Element s (tierras raras livianas).

83
Grupos funcionales
- -
sulfonato (SO3 ), carboxilato (CO2 ),
Resina intercambiadora de -
aminodiacetato (N(CH2CO2 )2 ), fosfonato
cationes 2- -
(PO3 ), fosfinato (PO2 ), etc..
+
amonio cuaternarios (NRH3 ), amonio
Resina intercambiadora de + +
terciarios (NR2H ), fosfonio (PR ), sulfonio
aniones +
( S ), etc...

Tabla 16: Clasificacin de resinas segn carga elctrica del grupo funcional

Cuando una resina determinada es sumergida en una solucin cargada de iones, en


funcin de la composicin de la solucin y de su pH, sta aceptar selectivamente
determinados iones ligndolos a su estructura qumica. Qu iones acepta y qu iones
depender no slo del (de los) grupo(s) funcional(es) de la resina, sino que tambin de la
composicin de la solucin. El principio de funcionamiento de una resina intercambiadora de
iones, consiste en explotar su especificidad para aceptar determinados elementos (iones) y no
otros, actuando como un verdadero filtro que retiene los elementos de inters. Es por eso que la
cromatograda se realiza en columnas, por las cuales se hace pasar la muestra en forma de
solucin, a un pH tal que la resina tomar los iones que se quiere separar. Enseguida se agrega
cido en cantidades y molaridades precisas y calibradas experimentalmente, eluyendo primero
los elementos que no interesan, para finalmente recuperar el o los elementos de inters
(comnmente la recuperacin se logra agregando cido cuya molaridad es mucho mayor a la
del cido usado para eluir).

Volviendo a la descripcin de las resinas como materiales, stas se componen de una red
de poliestireno a la cual se adosan grupos funcionales. El poliestireno es un copolmero de
estireno y de divinilbenzeno. El estireno forma cadenas polimerizadas alargadas y el
divinilbenceno forma ramificaciones o puentes entre dichas cadenas. Mientras ms contenido
de div. (divinilbenzeno) tenga la resina, mas ramificada y mas cerrada ser su estructura
tridimensional. El contenido de div. se expresa comnmente por la letra X seguida del
porcentaje en peso respectivo. Una resina X10 por ejemplo tiene 10% de div. En la practica, las
resinas intercambiadoras de iones estn limitadas a las resinas X1 a X30, y comnmente X8 a
X16.

La estructura macromolecular de una resina contiene cavidades, o poros, a las cuales


penetran molculas de agua y iones disueltos al sumergirla en una solucin acuosa. La solucin
al interior de la resina se encuentra confinada, sometida a las fuerzas elsticas de la red
macromolecular. Esta presin de confinamiento es igual a la presin osmtica de la solucin
interior. A su vez y a modo de fuerza de reaccin, la solucin interior ejerce una presin sobre la
red molecular de la resina, lo cual aumenta su volumen (la infla). Mientras mayor es el
contenido de divinilbenzeno de una resina con un grupo funcional dado, menor es su porosidad
y menor es su cambio de volumen. Y para una resina con un contenido de div. fijo, la presin
interna depende del grupo funcional y ms an del ion libre contenido en la resina. Son las
resinas con menor contenido en div. aquellas donde este fenmeno ocurre de forma mas
pronunciada. A grosso modo, para resinas >X10 las variaciones de volumen de una forma
inica a otra son muy dbiles.

Se define la capacidad de intercambio de una resina como el numero de grupos funcionales


por cantidad (masa, volumen) de resina seca. En general se expresa en moles por kg de resina
seca. De manera similar se define la concentracin de cualquier soluto al interior de la resina,
como la cantidad de estos por kg de resina seca. Los iones no estn libres dentro de una

84
resina, sino complexados o ligados a un grupo funcional determinado, a un radical
intercambiador.

Propiedades de las resinas en presencia de soluciones: intercambio de iones.

Al introducir una resina intercambiadora de iones en una solucin de un electrolito, pueden


producirse simultneamente dos fenmenos.

(1) Si el ion libre o liberable al interior de la resina es distinto del ion en solucin, pero tiene la
misma carga que este ultimo, se producirn intercambios entre ambos iones (iones
intercambiables). Dicho intercambio inico depende esencialmente de los grupos funcionales de
la resina.

(2) Una parte del electrolito (aniones y cationes) penetra al interior de la resina, siendo
compartido por dos fases. El mismo fenmeno tiene lugar con los no-electrolitos (molculas). El
fenmeno de penetracin de solutos en la resina se denomina generalmente como equilibrio de
Donnan.

Afinidad relativa de una resina por diferentes iones

Cuando una resina, por ejemplo que contiene iones de la forma P+, es puesta en contacto
con una solucin que contiene iones Q+, tendrn lugar intercambios de P+ y Q+ en las dos fases
hasta llegar a un estado de equilibrio (ausencia de intercambios). El equilibrio que se alcanza et
reproductible e independiente del sentido en el cual los iones se han intercambiado. Una vez
que se ha establecido el equilibrio, se constata que uno de los dos iones, digamos Q+, se
encuentra en mayor proporcin que el otro dentro de la resina. Se dice que la resina tiene una
mayor afinidad por el ion Q+ que por el ion P+.
Para cuantificar la afinidad de una resina por un ion Aa+/- respecto de otro ion B b+/- de carga del
mismo signo, debe considerarse la reaccin de intercambio siguiente.
ARa + BSb = ASa + BRb
donde y son nmeros primos relativos tales que : = a:b (coeficientes estequiomtricos
para balancear las cargas). La constante de equilibrio de intercambio inico (cantidad
adimensional y termodinmica) es:

K 0
=
[B ]R [A]S

[B ]
= R


[ A]S



[B ]S [ A]R
A/ B
[B ]S [A]R
donde [X]Y se refiere a la actividad de X en la fase Y (R: resina, S: solucin). Esta constante
depende de la temperatura a la cual ocurren los intercambios.

Es ms fcil utilizar concentraciones en lugar de actividades, y definir una constante


aparente de intercambio.

(C B )R (C A )S (C ) (C A )S
K A/ B = = B R


(C )


(C B )S (C A )R (C B )S A R
donde (C X)Y se refiere a la concentracin de X en la fase Y.
La constante KA/B tiene una relacin directa con los coeficientes de particin de los iones A y
B entre la resina y la solucin (PA y PB respectivamente):
(C ) (C ) (PB )
PA = A R y PB = B R K A/ B =
(C A )S (C B )S (PA )

85
La distancia entre las afinidades de una resina por un par de iones Aa+/- y B b+/- no esta dada
directamente por la constante de equilibrio respectiva, sino por su logaritmo:
[B ] [ A]
( )
log K A0 / B = log R log R
[B ]S [ A]S
Si los dos iones tienen la misma carga, todos los coeficientes estequiomtricos valdrn 1 y
en consecuencia el logaritmo de la constante de intercambio (sea termodinmica o aparente)
ser igual a la diferencia entre dos trminos

Variacin de KA/B (cte. aparente) con la composicin ionica de la resina

Consideremos un par de iones A y B que, al equilibrio, son compartidos entre una resina
y una solucin. Si CR es la concentracin de ambos iones dentro de la resina, CRx y CR(1-x)
sern las concentraciones de A y B en la resina respectivamente, siendo x la fraccin de
resina bajo la forma A. Experimentalmente se constata que el valor de la constante aparente
de intercambio KA/B varia frecuente y sustancialmente con la composicin inica de la resina, es
decir con x.

Constante termodinmica de intercambio

La constante de intercambio de dos iones en una resina determinada depende de dos


factores: la concentracin total y la proporcin de los dos iones en la solucin en equilibrio con
la resina. Sin embargo, existe una constante intrnseca cuyo valor es realmente independiente
de los factores y que tiene un significado termodinmico. Sin embargo, es ms prctico utilizar
la constante aparente de intercambio porque relaciona cantidades directamente medibles
(concentraciones en vez de actividades).

86
B.2 Resinas de extraccin

La cromatografa de extraccin es una tcnica que permite separar selectivamente


elementos desde distintos tipos de muestras. Es un tipo de cromatografa lquido-lquido en la
cual la fase estacionaria es un extractante o una solucin de un extractante, incluidos sobre un
soporte inerte. Combina la selectividad de la extraccin lquido-lquido con las ventajas de la
cromatografa en columna (proceso en etapas mltiples y facilidad de manejo). La unin de un
ion metlico a una resina de este tipo involucra la formacin de un complejo organometlico a
partir del ion hidratado respectivo, lo que permite el diseo de distintos esquemas de separacin
eficientes y selectivos. Las resinas de extraccin se preparan cubriendo o impregnando un
soporte slido con el extractante.

Los 3 principales componentes de un sistema de cromatografa de extraccin son el


soporte inerte, la fase estacionaria y la fase mvil (Figura 21). El soporte inerte puede ser slica
porosa o un polmero orgnico de 50-150 m de dimetro. La fase estacionaria se compone de
extractantes lquidos, solos o en mezclas, en tanto que la fase mvil puede ser una solucin
cida (p.e. cido clorhdrico o ntrico) o bien un complejante que permite aumentar la
selectividad del proceso de extraccin. La inmovilizacin del extractante en el soporte es un
proceso fsico de inclusin sin que ocurra una reaccin qumica entre ambos. Esto hace que las
propiedades de complejacin y selectividad del extractante en la columna sean muy parecidas a
las que presenta el extractante en solucin.

Resina Selectividad

Actinidos(III, IV, VI),


TRU-Spec
Ln(III)

Sr-Spec Sr, Pb

Ln-Spec Ln(III)

UTEVA U(VI)

Th(IV), Np(IV), Pu(IV),


TEVA
Tc(VII))

Res. para
Actinidos
Actinidos

Figura 21: Esquema de las estructura de la resina Tabla 17: Funcin de algunas resinas de
(fuente: pgina web de EICHROM) extraccin de EICHROM

87
Anexo C: Dilucin isotpica

La dilucin isotpica es una tcnica que permite medir de forma muy precisa la
concentracin de un istopo en una muestra. Para explicar su principio de funcionamiento, es
necesario definir primero el trmino "spike".
Spike: Solucin que contiene el istopo de inters (el que desea medirse en un muestra),
digamos "A", as como un istopo artificial (o natural) "B", en concentraciones conocidas
(medidas y calibradas) y cuya proporcin sea muy distinta de la que espera encontrarse en la
muestra (composicin isotpica natural).

Pasos a seguir:
o Mezclar el spike con la muestra, midiendo con una balanza analtica la masa de spike y
la masa de muestra introducidos (precisin relativa del orden de 10-5).
o Realizar las manipulaciones qumicas necesarias para obtener a partir de la mezcla, la
fraccin del (de los) elemento(s) cuyos istopos quiere analizarse.
o Medir mediante espectrometra de masa la proporcin entre los tomos de "A" (AM) y los
de "B" (BM) presentes en la mezcla.
o Realizar los clculos que se explican enseguida.

Cantidades relevantes:

Sean AM , Am y AS los nmeros de tomos de "A" en la mezcla (AM), en la muestra


(Am) y en el spike (AS) respectivamente.
Sean B M , Bm y BS los nmeros de tomos de "B" de la mezcla (BM), de la muestra
(Bm) y del spike (BS) respectivamente.
AM A AS
Sean RM = , R= m y RS = las proporciones atmicas de la mezcla (RM), de
BM Bm BS
la muestra (R) y del spike (RS) respectivamente.
Sean mm y mS las masas de la muestra (mm) y del spike (mS) respectivamente.
Sean M A y M B las masas atmicas de los istopos "A" (MA) y "B" (MB)
respectivamente.
Sean [ AM ], [ Am ] y [ AS ] las concentraciones de "A" de la mezcla ([AM]), de la muestra
([Am]) y del spike ([AS]) respectivamente.
Sean [ BM ], [ Bm ] y [ BS ] las concentraciones de "B" de la mezcla ([BM]), de la muestra
([Bm]) y del spike ([BS]) respectivamente.

Desarrollo:

Necesariamente (por consistencia con la definicin de concentracin) se cumplen las siguientes


ecuaciones.

88
Am M A = [ Am ] mm AS M A = [ AS ] mS y AM = Am + AS
B m M B = [ Bm ] m m BS M B = [ B S ] m S y BM = Bm + BS
Luego

mm m m
+ [ AS ] S
[ Am ] [ Am ] + [ AS ] S
A + AS MA MA mm
RM = m = =
Bm + B S m m [ Bm ] M A m MA
[ Bm ] m + [ BS ] S + [ BS ] S
MB MB MB mm M B

[ Am ] M B [ AS ] M B
y como R = RS = , sigue que
[ Bm ] M A [ BS ] M A

mS M A
[ Am ] + [ BS ] RS
mm M B
RM =
[ Am ] m MA
+ [ BS ] S
R mm M B

Reordenando los trminos se llega a la frmula de dilucin isotpica:


1
m MA
[ Am ] = [ BS ] S (RM RS ) 1 RM
mm M B R

A continuacin se mostrar, bajo ciertos supuestos, que el factor (1-RM/R)-1 es muy cercano
a la unidad y por ende slo es necesario estimar de forma gruesa el valor de "R" (sea
reemplazndolo por "RN" o haciendo otro supuesto).

Supuestos:
Bm RS
Bm << BS << 1 y RS << R << 1
BS R

Si bien estas dos condiciones no pueden verificarse a priori, s puede esperarse que
cumplan si se satisfacen las dos condiciones siguientes.

A. El valor de "R" y "Bm" no se alejan sustancialmente de los valores promedio en el


entorno natural muestreado: "RN" y "BN" respectivamente. De forma explcita:
R = RN y Bm = B N
con [101 ,10] y [0,10]
RS B
B. 10 << 1 y 10 N << 1
RN BS
Estas dos condiciones implican necesariamente,
Bm B B RS 1 RS R
= N 10 N << 1 y = 10 S << 1
BS BS BS R RN RN

89
Bm RS R
pdq. si << 1 y << 1 entonces M << 1
BS R R
Es claro que "RM" ser una cantidad intermedia entre "RS" y "R", en efecto:
1 Am 1 AS
+
A + AS B B Bm B S
RM = m = S m = (1 ) R + RS
Bm + B S 1 1
+
B S Bm
1
Bm 1
donde = =
1 1 B
+ 1+ m
B S Bm BS

i.e., "RM" es una combinacin convexa entre "RS" (=1) y "R" (=0). Sigue que
RM R
= 1 + S 1 donde 1
R R
RM RS
Entonces << 1. q.e.d .
R R

90
Anexo D: Clculo de razones isotpicas

En la naturaleza, la razn 238U/235U es prcticamente constante y vale 137,88. Este valor se


usa para calcular 234U/238U a partir de 234U/235U:
234
U 234
U 1
= 234

238
U 235
U 137,88 238
La frmula de dilucin isotpica, aplicada a los istopos 238U y 233U (va el 235U), resulta en
mSpike M 238U 235
238 233
[U muestra ] [ U muestra ] =[ U spike ] (1 + 137,88) 233U mezcla
mmuestra M 233U U mezcla
232 229
Y para el Th - Th, la frmula de dilucin isotpica queda
mSpike M 232Th 232
Thmezcla
[Thmuestra ] [ 232Thmuestra ] =[ 229Thspike ] 229
mmuestra M 229Th Thmezcla
226
La concentracin de Ra se calcul con la frmula siguiente.
m Spike M 226 Ra
[ 226 Ramuestra ] =[ 226 Ra spike ] (Rmezxla R spike ) / (1 Rmezcla / Rmuestra )
m muestra M 228 Ra
226
Ra
donde R = 228
Ra
Para las razones isotpicas entre Th y U se usaron las siguientes frmulas.

1
230
Th 230 232
Th
=

232
Th 232
230
Th
U [U ] 137,88 M 232Th 238
238

Th [Th] 1 + 137,88 M 238U 232
232

1
230
Th 230
Th 238
U
=

238
U 232
Th 232
Th
1
230
Th 230
Th 234
U
=

234
U 238
U 238
U
1
226
Ra [ 226Ra] M 232Th 226 230
Th
=

230
Th [Th] M 226 Ra 232
232
Th

91

También podría gustarte