SINTESIS GEOLOGICA Y CARTOGRAFIA PRELIMINAR DE LOS

MINERALES EVAPORITICOS DEL SALAR DE COIPASA

UTILIZANDO IMÁGENES ASTER

Fernando Cáceres Jerez

Eric Pirard
Université de Liège
Département des Géoressources, Géotechnologies
et Matériaux de Construction (GEOMAC)
Sart Tilman B52/3
4000 Liège
Belgique
Tel: +32 43669525
Fax: +32 43669520
[email protected]
Jorge Gutiérrez Venegas
Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Carrera de Geología
Centro de Investigaciones en Teledetección (CEINTEL-UTO)
Tel: +591 225261375
Fax: +591 225272835
[email protected]

Abstract

The use of remote sensing techniques for mapping of evaporitic minerals has been widely promoted and tested in serveral
regions around the world. In this work a set of VNIR and SWIR images from the ASTER sensor have been used in order to map
secondary minerals known to occur in the Salar de Coipasa. Spectral unmixing techniques have been used to identify the
mineral components and map their abundancy. A subset of the Salar de Coipasa is presented wherein two evaporitic minerals
have been identified by their diffuse reflectance spectrum (mirabilite and gypsum). The spatial distribution of gypsum is in
good accordance with field descriptions made during the geological exploration programme.

Resumen

La utilización de la técnicas de teledetección para la cartografía de minerales evaporíticos fue grandemente utilizado. En
este trabajo se utilizan las imágenes VNIR y SWIR del captor ASTER con el objetivo de cartografiar los minerales secundarios
presentes en el Salar de Coipasa. Mediante las técnicas de desmezcla espectral aparte de identificar minerales presentes en
la imagen, permite también cartografiarlos. En la imagen subset de Coipasa se pudo constatar la presencia de dos minerales
evaporíticos identificados en base a su firma espectral (yeso y mirabilita).La distribución del yeso coincide con las
descripciones de terreno realizadas durante un programa de exploración.

Palabras claves: Salar de Coipasa, Sensores Remotos, Minerales Evaporiticos, Desmezcla Espectral.
1.- Introducción contenido de boro en el salar de Uyuni, con este propósito
Este articulo es una contribución más, a la larga literatura los autores utilizan la relación TM 4/TM 7.
que existe sobre el estudio de secuencias evaporíticas Bryant (1996) utiliza la técnica de desmezcla espectral sobre
utilizando sensores remotos. Los ópticos (multi – imágenes Landsat TM de la región Chott el Djerid (Tunes)
hiperespectrales) y radares demostraron su potencial en la con el objetivo de determinar la distribución espacial de
cartografía y monitoreo de este tipo de ambiente. A esta yeso, halita y carnalita. Bryant (2003) emplea imágenes
experiencia podemos adjuntar la facilidad que tienen estas AVHRR en el monitoreo de playas en el Norte y Sud del
técnicas de cubrir grandes regiones con un bajo costo (50 continente africano. En este tipo de ambiente los cambios
centavos de $us americanos/km2). hidrológicos pueden ser detectados por el análisis de
Los resultados presentados en este trabajo tienen un carácter imágenes satelitales, estas variaciones permiten determinar
preliminar, los cuales están siendo confirmados y de manera indirecta los cambios de precipitación y
extrapolados a todo el salar con la elaboración de una tesis evaporación del sistema.
realizada por la egresada de la Carrera de Ingeniería Camargo, et al. (2003) emplea la información de
Geológica, Lucía Jáuregui S. reflectancia adquirida in situ sobre el salar de Uyuni con el
fin de calibrar las imágenes satelitales, los autores
2.- Trabajos anteriores recomiendan la utilización de este tipo de ambientes por
La discusión científica del origen y geoquímica del salar de su elevada homogeneidad espectral y los débiles efectos
Coipasa comienza con Ericksen (1976) que publica los atmosféricos que existen en la zona.
primeros análisis de las salmueras superficiales; después Otros autores emplean captores aerotransportados para el
Retting, et al. (1979), Ballivián & Risacher (1981), Risacher & estudio de este tipo de ambiente, como es el caso de la
Fritz (1991) y USGS-Geobol (1992) estudiaron los salares y región de Death Valley (California). Con este propósito
sus salmueras. Crowley & Clark (1992) y Crowley (1993) utilizan
Entre 1996-1999 la Cooperación Universitaria belga junto a imágenes que provienen del captor AVIRIS (Airborne
la Universidad de Lieja (Bélgica) y a la Universidad Técnica Visible/Infrared Imaging Spectrometer) y Crowley, Mars
de Oruro (Bolivia) estudian el Salar de Coipasa con el fin de & Hook (2000) emplea imágenes MASTER (MODIS/
evaluar los recursos evaporíticos de la cuenca (Proyecto ASTER Airborne Simulator).
Coipasa). Con esta premisa el proyecto realizo dos estudios
principales: uno geofísico utilizando la gravimetría y el
estudio geoquímico de las salmueras; los dos trabajos fueron
realizados en una malla regular (2 km x 2 km) que cubre todo
el salar. A estos trabajos se incluyen 200m totales de
perforación, repartidos en 38 pozos en diferentes sectores
del salar; el pozo más profundo alcanza la profundidad de 15
m en el centro del salar. De este proyecto salieron los trabajos
de Pool (1995), Fis (1996), Gutierrez (1997), Lebrun (1997),
Pacosillo (1998), Houssa (1998) y Lebrun, et al. (2002).
Sylvestre, Servant-Vildary & Servant (1998) de la ORSTOM
en 1998 realizan una perforación de 4.87 m de profundidad
en el centro del salar con el propósito de estudiar el paleo-
clima de la región. En este resumen no incluimos los
innumerables trabajos realizados por YPFB y SERGEOMIN
que tuvieron como objetivo conocer la geología de la región.
Chapman, et al. (1989) es uno de los primeros en utilizar las
imágenes satelitales para el estudio de este tipo de ambientes,
a partir de imágenes Landsat Tematic Mapper (TM), dichos
autores emplean la técnica de decorrelation stretching con
el fin de obtener la información mineralogíca del Salar de
Atacama (Chile). White & Drake (1993) utilizan las imágenes
Landsat TM para cartografiar la abundancia y distribución
de yeso en la zona de Chott Fedjaj (Tunes)
Sabins & Millar (1994) y Sabins (1997 y 1999) emplean
imágenes Landsat TM y AVHRR (Advanced Very High
Figura 1: Localización de las principales cuencas del Altiplano
Resolution Radiometer) para la determinación del Boliviano.
3.- Los salares del Altiplano boliviano La parte Norte del salar esta cubierta por un lago permanente
El Altiplano es una cuenca endorreica de ±200.000 km a 2 alimentado por el río Lauca, principal tributario del salar. En
4000 msnm se encuentra entre la Cordillera Occidental y la época de lluvias (noviembre a marzo) todo el salar es
Oriental, por sus características meteorológicas que presenta recubierto de una capa de salmuera de 15 cm de alto.
la cuenca tiene las condiciones ideales para la formación Hacia el Sud del Salar de Coipasa se encuentra la costra de
de ambientes evaporíticos, por otra parte ella encierra un sal más grande del planeta, el Salar de Uyuni (10.000 km2),
sin numero de cuencas evaporíticas. su principal tributario el Río Grande desemboca en la parte
Sud del salar. El Sud del Altiplano, geográficamente
3.1.- Geología del Altiplano conocido como los Lipez, presenta numerosos lagos salados
Las Cordilleras que encierran la Cuenca Altiplánica son de y salmueras (más de 30) en el seno de pequeñas cuencas de
naturaleza diferente, así la Cordillera Oriental esta origen volcánico. El quimismo de las salmueras de esta
constituida por formaciones paleozoicas y por el contrario región es variado. En la figura 1 se puede observar la
la Cordillera Oriental que pertenece al Complejo Volcánico localización de las diferentes cuencas del Altiplano.
Neógeno – Cuaternario de los Andes Centrales.
El substrato Precámbrico y Paleozoico del Altiplano está 3.3.- Historia lacustre del Altiplano
recubierto por una espesa serie volcánica-sedimentaria (que Otro tema de controversia y de gran estudio en estos últimos
en algunas regiones tiene un espesor de 15.000 m) de edades años son los diferentes eventos lacustres que existieron en
que varían del Cretácico al presente. el Altiplano. En los salares de Uyuni y Coipasa y
alrededores se encuentra la memoria climática de los últimos
3.2.- Hidrología del Altiplano 200.000 años del Altiplano boliviano, tema que apasiona a
La mayor parte de las cuencas presentes dentro del Altiplano muchos científicos nacio-nales e internacionales. En la figura
se extienden según un eje Norte-Sud (figura 1). El lago 2 se observan las diferentes fases lacustres que se
Titicaca (3808 msnm) es la cuenca más septentrional y su desarrollan en el Altiplano norte (lago Titicaca) y centro. Se
salinidad media es de 1 g/l. El agua del lago es evacuada tiene demostrado que el nivel de los paleólagos decrece de
por el Río Desaguadero el cual desemboca en el lago Poopó, una fase lacustre a la siguiente (Servant & Fontes, 1978).
que es un lago salado con una salinidad promedio de 25 En la figura se puede observar que algunos de estos niveles
g/t. Al SW del lago Poopó se encuentra el Salar de Coipasa de agua no se extienden hacia el Sud y el umbral de salida
(3656 msnm), objeto de este estudio; este salar contiene del lago baja con el curso del tiempo; este proceso permite
una costra de halita de 2500 km2. que el nivel de los diferentes lagos descienda.

Figura 2: Sección N-S del Altiplano Central para observar la relación que existe entre los lagos y salares (modificado de Risacher,
1992)
El lago Titicaca alcanzo su nivel actual hace 10.000 años El estudio de los fenómenos evaporíticos esta en relación a
atrás, en la fase árida del Holoceno (entre 7700 y 3650 este último punto; a grandes rasgos los mecanismos de
años atrás) y su nivel bajo de 50 m. Es importante remarcar saturación que pueden conducir a las variaciones de
que esta fase fue puesta en evidencia en otras regiones de concentración de las aguas naturales son:
América Latina. Risacher & Fritz (1992) igualmente
reconocen esta fase árida en el Salar de Uyuni gracias a la 4.1.- Balance Hídrico
presencia de polyhalita, mineral que precipita en salmueras El balance hídrico cuantifica las aguas que entran y salen
concentradas cinco veces más que las actuales. Después de la cuenca y esta en paralelo al balance químico, el cual
de esta fase árida el lago Titicaca alcanzo su nivel actual de mensura los solutos que entran y salen de la cuenca. Del
3808 m. estudio de diferentes posibilidades hídricas se puede
Dos niveles de endurecimientos algareos localizados entre deducir:
3700 y 3747 m tapizan todo la cuenca Pocoyu (denominación a) La temperatura juega un rol importante en la existencia
proporcionada a la cuenca Poopó - Coipasa - Uyuni por de un medio evaporítico, condición que es más importante
Wirrmann & Mourguiart, 1995 y Argollo & Mourguiart, que la precipitación > evaporación; sin embargo esta
2000). En función a esta información pudieron ser puestos condición no es suficiente.
en evidencia dos paleólagos en la parte Central del Altiplano, b) El establecimiento de un medio evaporítico, es decir la
las fases fueron datadas 27.000 años para el lago Minchin y existencia de una cuenca geológicamente y
13.000 a 10.000 años para el Tauca (figura 3). Fornari, Risacher topográficamente cerrada e impermeable. La
& Féraud (2001) luego de la datación de los 121 metros impermeabilidad es un factor importante, ya que cuanto
perforados en el salar de Uyuni confirman la existencia del más una cuenca es taponada, más elevada será la salinidad
evento Escara (191.000 años atrás), una fase que no es de sus aguas.
tomada en cuenta por muchos autores. c) La salinidad del agua aumenta con la aridez del clima y
Las dataciones realizadas sobre las construcciones se puede constatar que: en climas no áridos las sales son
arreciferas son erráticas a causa de la presencia de eliminadas por los afluentes, en climas muy áridos el agua
fenómenos de recristalización, sin embargo la correlación no es suficiente y los aportes a la cuenca son débiles; las
de los mencionados paleólagos con las costras algareas aguas son concentradas en una cantidad muy reducida y la
pudo claramente ser establecida. ausencia total de precipitación conduce a una ausencia total
El estudio de las diatomeas como la presencia de Coprolitos del deposito salino.
de Artemia (pequeños crustaceos que viven en medios
salinos) dentro los sedimentos de estos paleolagos muestran
que ellos eran salados. Probablemente las costras de sal de
los salares de Uyuni y Coipasa provienen de la sequía hace
10.000 años atrás del lago Tauca (Risacher, 1992).
Sylvestre, et al. (1999) pone en evidencia otro evento lacustre
al cual denomina Coipasa que ocurrió hace 9500 y 8500
años atrás, los últimos trabajos realizados por Fornari,
Risacher & Féraud (2001) determinan que el evento
Coipasa es una oscilación del gran lago Tauca.

4.- Generación de un fenómeno evaporítico

Los conocimientos actuales de los procesos evaporíticos
fueron grandemente descritos en los trabajos de Hardie &
Eugster (1970 y 1978) y Risacher (1992). En este acápite
presentamos un resumen de los mencionados artículos.
La alteración de las rocas por la lluvia cargan naturalmente
las aguas superficiales de especies químicas disueltas, en
el trayecto que estas realizan de las nacientes hacia el centro
de la cuenca. Con este proceso los principales solutos que
son trasportados por los ríos son: Na+, K+, Ca++, Fe++, Mg++,
HCO3-, CO3--, SO4-- y H4SiO4. Los solutos acumulados dentro
del agua pueden quedarse en solución o precipitar en razón
de las modificaciones termodinámicas del medio:
modificación del Eh, pH, P, T y el aumento de la
concentración hasta la saturación por evaporación. Figura 3: Extensión del paleólago Tauca.
4.2.- Geoquímica de las evaporitas continen- En este caso la constante del producto [CO3--][Ca++] ó
[Ca++][HCO3-]2 será la concentración simultanea de calcio
tales
y de iones carbonatos. En el curso de la precipitación de la
Los trabajos de Hardie & Eugster (1970 y 1978) y Risacher
calcita, el elemento inicial en exceso se concentrara y el
(1992) ponen en manifiesto la ligazón que existe entre las
elemento inicial en déficit se diluye. En la práctica las aguas
fases químicas que se depositan y la historia sedimentaria
que se encuentran en la naturaleza contienen a la vez
de la cuenca, a continuación proporcionamos un resumen
carbonatos y bicarbonatos, por lo que es mejor trabajar con
de un proceso evaporitico (figura 4):
el índice de alcalinidad (alc):
El primer mineral que precipita es la calcita (CaCO3), que es
alc = 2[CO3--] + [HCO3-]
un mineral insoluble. Dos reacciones pueden guiar a esta
Si el alc es superior a 2[Ca++] en la solución inicial, ella
precipitación, los iones carbonatos (CO3--) ó los iones
evolucionara hacia la zona alcalina, se concentrara
bicarbonatos (HCO3-).
progresivamente en álcalis y su pH se establecerá entre 8 y
Ca++ + CO3-- => CaCO3
9. Si alc es inferior a 2[Ca++], entonces la solución seguirá
por la ley de acción de masas
el camino de la salinidad neutra, ella se enriquecerá en calcio
[CO3--][Ca++] = K1
y su pH se estabilizara entre 7 – 8.
y
Dentro el camino neutro el próximo mineral a precipitar será
Ca++ + 2HCO3- => CaCO3 + H2O + CO2
el yeso, CaSO 4.2H 2O, la reacción que controla esta
o
precipitación puede escribirse como:
[Ca++][HCO3-]2 = K2[H2O]PCO2
[Ca++][SO4--][H2O] = K3
Donde [H2O] es la actividad del agua que tiene un valor ±1,
Donde, K3 es el producto de solubilidad del yeso. La
PCO2 es la presión parcial de CO2 dentro la atmósfera y es un
solución se concentra en sulfatos o en calcio según el origen
valor constante. K1 y K2 son el producto de solubilidad
de la precipitación; si existe un exceso de uno u otro
dentro del agua y ellos dependen solamente de la
elemento, entonces la salmuera evolucionara hacia una fase
temperatura y de la presión. Cuando los iones carbonatos o
de tipo Na-Ca-(Mg)-Cl o Na-(Mg)-Cl-SO4.
bicarbonatos son predominantes dentro la solución una de
las dos ecuaciones es suficiente.

Figura 4: Diagrama teórico seguido por las aguas en un fenómeno evaporítico (Hardie & Eugster, 1970 y Risacher, 1992).
Dentro el camino alcalino el déficit de calcio limita la Por lo contrario los salares que bifurcan hacia condiciones
formación de yeso y conduce a la formación de silicatos de salinad neutra tienen un gran estrato de sedimentos
magnésicos según la reacción: arcillosos, arcillos – arenosos y salíferos que encierran la
Mg++ + nH4SiO4 => Silicato – Mg + H+ napa subaflorante. Entonces, cuando las aguas de la fuente
Los iones H+ liberados por la reacción neutralizaran una o del rió penetran dentro los sedimentos localizados sobre
parte de los iones carbonatos. Lo cual puede invertir la el borde del salar, algunos de los siguientes fenómenos
evolución de la solución hacia un camino salino neutro. La pueden producirse:
condición para la existencia de esta inversión es de tener en a) Disolución del yeso: el yeso es siempre abundante en el
el agua inicial: alc < 2[Ca++] + 2[Mg++]. Si esta reacción no borde del salar y tiene como origen la sequedad del lago
se verifica, entonces la solución continúa dentro del camino precedente o la evaporación por capilaridad de la napa
alcalino. más profunda. Las aguas que entran en la cuenca adicionan
Risacher & Fritz (1991) constatan que en un buen número Ca++ en la solución, lo que hace disminuir la relación: alc/
de aguas del Altiplano no se verifica esta condición, sin (2[Ca++] + 2[Mg++])
embargo evolucionan hacia una facie salina neutra. Esto b) La transformación del vidrio volcánico en arcilla: las
quiere decir que la composición química de partida no es aguas de aporte tienen en suspensión vidrio volcánico
una condición suficiente para explicar el futuro químico proveniente de las formaciones volcánicas de la cuenca.
del agua. Otra hipótesis pone de manifiesto un segundo Este vidrio volcánico dentro del salar se transforma en
factor que influye en la evolución de la salmuera, dicho arcilla con la liberación de iones H+ que tienen la tendencia
componente esta en relación con la interacción del agua de neutralizar los carbonatos.
con los sedimentos. Esta teoría esta basada en las observa-
ciones efectuadas en los salares tipo carbonato sódico que
tienen pocos sedimentos lacustres debajo la costra de sal.

Figura 5: Localización de la región estudiada.

c) La composición de las rocas de la cuenca: si estas rocas En la época de lluvia, el agua concentrada del lago inunda
son cristalizadas, volcánicas, sin impurezas ni toda la costra de sal y los iones son dispersados en todas
mineralizaciones, ellas no proporcionaran ningún ión en partes por difusión. En el Salar de Coipasa los ríos tributarios
cantidades apreciables a la cuenca. Entonces el anión mayor tienen el papel de homogenizar las salmueras cuando la costra
será HCO3- de la atmósfera, el cual conduce sin remedio esta completamente inundada. Esta hipótesis explica mejor
hacia el polo alcalino. Si por lo contrario en la cuenca existe el enriquecimiento que existe en salmueras en la parte SE del
mineralización o rocas sedimentarias, ellas liberan sulfatos salar y su empobrecimiento en la parte Oeste (zona sin
y cloruros que empujan a las salmueras hacia una salinidad grandes aportes de agua).
neutra, como es el caso de Coipasa. Lebrun et al. (2002) pone de manifiesto una elevada diferencia
de la relación Mg/Li en el Salar de Coipasa con respecto al
5.- El Salar de Coipasa Salar de Uyuni. Esta diferencia tiene su origen en la
El Salar de Coipasa a diferencia del Salar de Uyuni, tiene en composición y proporción de las rocas que se encuentran
su parte Norte, en la desembocadura del río Lauca, un lago dentro la cuenca.
permanente que juega un rol importante en la geoquímica El Li, Mg y B esencialmente provienen de las rocas
del salar. volcánicas; por otro lado el Mg y K son abundantes en otro
Así mismo dentro el salar existe una isla (volcán Pucarani) tipo de rocas. Si la proporción de rocas volcánicas o su
que ocupa la parte NW del salar, sobre la cual se encuentra composición es diferente de una cuenca a otra, esto inducirá
el pueblo civil de Coipasa (figura 5). en la diferencia de la relación Mg/Li en el agua de la cuenca.
Los principales tributarios del salar aparte del Río Lauca El relativo enriquecimiento de sulfatos en las salmueras de
son el Río Sabaya al Norte, el Río Laca Jahuira al Este y el Coipasa esta directamente relacionado a las concentraciones
Río Moscoma al Sud. de sus tributarios. El Río Lauca tiene una mayor concen-
tración de sulfatos en comparación al Río Grande. El origen
5.1.- Geoquímica del salar de este enriquecimiento de dichos sulfatos es difícil a
verificar, estos pueden estar relacionados al viento, el cual
Las salmueras de Coipasa según la clasificación de Eugster
transporta los sulfuros de la actividad volcánica de la zona y
& Hardie (1978) son de tipo Na-Cl-(SO4), la presencia entre
luego los deposita en el drenaje de la cuenca (Risacher &
paréntesis de SO4 significa que los sulfatos son un menor
Fritz, 1991). Otra posible explicación, es la presencia de
constituyente del salar.
sulfuros en la cuenca, como fue descrito en el salar de
Risacher (1992) describe 2.5 m de espesor de costra de sal
Hamburgo en Chile (Alpers & Whittemore, 1990).
constituida por halita en su parte superior y por mirabilita
(Na2SO4.10H2O) en la parte inferior. En las 38 perforaciones
realizadas por el proyecto Coipasa se encontró solamente
10 a 60 cm de costra de halita. En la parte inferior, entre 1 à Subs ys te ma
Banda Rango Re s olución Nive l de Tamaño de la
N o. Es pe ctral (um) Es pacial (m) Quantificación image n (km)
1.5 m el salar esta generalmente constituido por sedimentos
lacustres, estos sedimentos están compuestos de arcillas 1 0.52- 0.60

más o menos arenosas y alternancias de arenas. Dichos 2 0.63- 0.69

VNIR 15 8 bits 60
sedimentos pueden o no contener cristales de mirabilita en 3N 0.78- 0.86
diferentes proporciones. A 1.8 m de profundidad en el pozo 3B 0.78- 0.86
7 y a 1.3 m en los pozos 4a, 4b y 4c (figura 5), las perfora- 4 1.60- 1.70
ciones alcanzaron cristales de mirabilita; este mismo manto
5 2.145- 2.185
tiene un metro de espesor en el pozo 4c. En todos los pozos
6 2 . 18 5 - 2 . 2 2 5
realizados para muestrear la salmuera superficial la costra SWIR 30 8 bits 60
de halita no excede de los 80 cm; bajo la costra se encuentra 7 2.235- 2.285

una capa de arcilla lodosa negra (Lebrun, et al. 2002). 8 2.295- 2.430

El contenido de Mg, K y Li en el salar proviene del lavado 9 2.360- 2.430

de las rocas presentes en la cuenca, ellos son transportadas 10 8.125- 8.475
por los ríos tributarios hacia el salar; dichos elementos no 11 8.475- 8.825
son precipitados como parte constituyente de ningún
TIR 12 8.925- 9.275 90 16 bits 60
mineral evaporítico. Si comparamos con el Salar de Uyuni,
13 10.25- 10.95
los ríos tributarios juegan otro papel en la geoquímica del
salar, en Coipasa, el lago permanente opera como un 14 10.95- 11.65

concentrador durante la estación de estío, pero con mayor Tabla 1: Información espectral y espacial del captor de ASTER.
eficacia que los sedimentos deltaicos del salar de Uyuni
(Hardie & Eugster, 1970).
6.- Datos En este trabajo solo se toma en cuenta las imágenes de los
El trabajo esta particularmente basado en las imágenes del subsistemas VNIR y SWIR, no se considera el subsistema
captor ASTER (Advanced Spaceborne Thermal Emissión TIR ya que las mensuras dentro este subsistema y según la
and Reflection Radiometer) que es uno de los cinco sensores ley de Boltzmann toman en cuenta tres propiedades físicas
embarcados sobre el satélite EOS-AM1 y denominado del material (absorción, reflexión y temperatura), por lo
después TERRA plataforma, puesto en orbita en 1998 y que contrario los otros dos subsistemas solamente mensuran
remplaza al satélite EOS de la NASA. uno (reflexión).
Su orbita a 700 km de altura, heliosincronica, cercana a los Con el fin de realizar un ensaye del método presentado en
Polos y el tiempo que pasa a nivel del Ecuador permiten este artículo, todo el trabajo es realizado sobre un subset
minimizar los efectos perjudiciales causados por las nubes. de la imagen, el cual cubre un área aproximada de 40 km2. La
El satélite tiene una frecuencia de repetición de 16 días, es localización del subset con respecto al salar puede
decir que después de 16 días toma una imagen de la misma observarse en la figura 5. La zona fue elegida por la variedad
región. mineralógica que se observa en superficie; en el punto H4 ;
El captor ASTER fue construido por el Ministry of V-5 (19° 25.347’ latitud sud y 68° 20.613’ longitud W) el
Internacional Trade and Industry (MITI) del Japón y esta proyecto Coipasa aparte de la halita encontro también
compuesto de tres subsistemas: cristales de yeso.
El primer subsistema VNIR (Visible and Near Infrared) cubre
la región visible e infrarrojo cercano del espectro 7.- Método
electromagnético. El subsistema es un escáner tipo 7.1.- Pretratamiento de las Imágenes
“pushbroom” formado de dos telescopios uno a nadir que Sobre las imágenes brutas existen varios errores técnicos
tiene tres bandas y otro localizado en la parte posterior con puestos en evidencia después que el satélite se encontraba
una simple banda; este canal es utilizado para las en orbita (errores de paralax y crosstalk), con este objetivo
observaciones estereoscópicas. Los detectores CCD son se tuvo que crear todo un procedimiento para la calibración
lineales construidos en Si, cada captor posee 5000 de las mismas. Primeramente se corrigió el error de crosstalk
detectores, los detalles de la resolución espectral y espacial de la imagen; este problema es creado por el pésimo
del subsistema son mostrados en la tabla 1. aislamiento de algunos componentes electrónicos los cuales
En segundo subsistema SWIR (Shortwave Infrared), con reflectan la luz de sol hacia los captores y generan
sus 6 bandas cubre el infrarrojo cercano y medio, los fantasmas en las imágenes (Iwasaki, et al., 2002). Este
detectores son también de tipo “pushbroom” y opera con problema es observado en las bandas 4 y 9 del subsistema
un simple telescopio que apunta a nadir, los datos son SWIR. El método de corrección empleado en este trabajo
adquiridos por un espejo. Los detectores CCD son lineales fue desarrollado por los mismos autores.
y construidos en PtSi, cada banda espectral posee 2048 El error de paralax se encuentra sobre todas las bandas del
detectores. En la tabla 1 se observa la resolución espectral subsistema SWIR y es generado por la gran distancia que
y espacial del subsistema. existe entre captores, el error es de 0.5 píxeles para cada 900
El tercer subsistema TIR (Thermal Infrared) con sus 5 bandas m de elevación (Fujisada, 1998). Este problema fue corregido
mensura la parte termal del infrarrojo; es un escáner mediante la ortorectificación de todas las imágenes, para
mecánico compuesto de un simple telescopio a posición ortorectificar las imágenes se utilizo un modelo digital del
fija observando a nadir y de un escáner con un espejo terreno construido en base a la información adquirida por el
localizado en la parte posterior. Cada banda utiliza 10 mismo satélite
detectores PC construidos en HgCdTe, para preservar la
dinámica de la imagen ella es codificada en 16 bits. La
información espacial y espectral se puede observar en la
7.2.- Corrección Atmosférica
Con el fin de eliminar el efecto atmosférico que se encuentra
tabla 1.
en la imagen se utilizo el método IARR (Internal Average
Los datos que adquiere el satélite tienen una amplia
Relative Reflectance), este método calcula el espectro medio
cobertura espectral y una resolución relativamente elevada,
de toda la imagen y luego divide este espectro por cada
lo que permite discriminar una serie de materiales sobre la
uno de los píxeles de la imagen. De esta manera la imagen
superficie (Yamaguchi et al. 1998).
de radianza es transformada en información de relativa
En este trabajo utilizamos una imagen adquirida el 19/03/
reflectancia.
2001 a 7 h 07 m 15 s (hora universal), el nivel de corrección
Este método genera artefactos en las zonas de fuerte
es 1B y es una de las cuatro imágenes que cubre enteramente
absorción y puede llevar después a falsas interpretaciones
el salar.
de la imagen (Van der Meer, de Jong & Bakker, 2001).
La zona estudiada presenta una gran uniformidad espectral y En la literatura existen muchos métodos para extraer los
mineralógica, nosotros pensamos que no existe una creación píxeles más puros de una imagen, en este trabajo se utiliza
de artefactos. El método puede ser utilizado en regiones un método basado en la información espectral. Boardman et
donde no existe información meteorológica, como es el caso al., (1995) sugiere la técnica PPI (Píxel Purity Index), la cual
de Coipasa. Los modelos más complejos proporcionan se basa en las propiedades geométricas de un conjunto
mejores resultados pero desgraciadamente requieren de convexo y considera a cada píxel como un vector en un
mayor información; debemos mencionar también que la espacio a n-dimensiones (n-D). El algoritmo no es aplicado a
información espectroscópica adquirida sobre el terreno la imagen original sino a la imagen resultante de una
puede servir para eliminar el efecto atmosférico. reducción dimensional de datos basada en el método MNF
(Minimun Noise Fraction) desarrollada por Green, et al., (1988)
7.3.- Desmezcla Espectral El algoritmo comienza con la generación de un vector
La desmezcal espectral es el método por el cual un espectro aleatorio a n-D (denominado skewer), en el cual todos los
mensurado de un píxel es descompuesto en una colección píxeles de la imagen son proyectados, los píxeles que
de espectros constituyentes llamados “endmembers” y un aparecen en la parte extrema del skewers son marcados y
conjunto de las fracciones correspondientes denominada luego contados, este procedimiento es generalizado a un
“abundancia”; este ultimo parámetro indica la proporción gran numero de skewer generados de igual manera. Al final
de cada endmember presente en un píxel (Keshva & Mustard de este proceso se obtiene una imagen con los píxeles más
2002). El método se puede subdividir en tres etapas, las contados, es decir los píxeles que fueron contados más de
cuales son enteramente independientes la una de la otra: 1000 veces son retenidos como posibles endmembers. Los
candidatos a ser considerados como verdaderos
endmembers son seleccionados en base de una comparación
7.3.1.- Determinación de los valores extremos con una librería espectral y explorando su posición en el
Un endmember puede traducirse literalmente como un valor espacio n-D.
extremo y técnicamente se traduce como un píxel puro, es Esta selección supervisada requiere un gran conocimiento
decir un píxel que contenga solamente un material. de las técnicas de teledetección y espectroscopia mineral. El
Generalmente dicho píxel corresponde a un objeto mayor defecto de este método es la pobre interactividad y
macroscópico y/o microscópico presente en la escena, como los resultados no pueden ser empleados para analizar un
por ejemplo el agua, suelo, vegetación, etc. El término también gran número de imágenes. Muchos autores consideran que
engloba la utilización de información espectral proveniente la naturaleza aleatoria de los skewers es otra gran limitación
de librerías espectrales obtenidas en condiciones de del algoritmo (Plaza, et al., 2002). Todo el procedimiento
laboratorio. seguido para la obtención de los valores extremos se puede
observar en la figura 6.

7.3.2.- Modelo de la Desmezcla

En la naturaleza existen cuatro tipos de mezcla espectral (Van
der Meer, 2001):
a) Las mezclas de superficie, en las cuales los componentes
dispersos y los objetos que están dentro del campo de vista
del captor pueden ser separados ópticamente, este tipo de
mezclas pueden ser consideradas como lineales.
b) Las mezclas intimas, se producen cuando los materiales
están en contacto con la superficie de dispersión, por ejemplo
los granos de mineral en los suelos o dentro la rocas. En
función de las propiedades ópticas de cada uno de los
componentes, el resultado de la señal es una combinación
no lineal de los espectros.
c) Coatings, este tipo de mezcla ocurre cuando el material en
estudio esta enterrado por otros materiales en forma de
estratificación, en este caso la dispersión y la transmisión
varían en función del espesor del material y la longitud de
onda. En este tipo de desmezcla las propiedades espectrales
de los materiales son consideradas como no lineales.
Figura 6: Procedimiento seguido para la extracción de los endmembers
y la abundancia de la imagen.
d) La mezcla molecular, acaece a nivel atómico, por ejemplo Formula
Clas e M ine ral Picos de Abs orcion (µm)
dentro de dos líquidos o entre líquido y sólido, en este caso Química
la mezcla es también considerada como no lineal. Cloruros Halita NaCl 1.93
El tipo y la naturaleza de las mezclas son cruciales para el
conocimiento de la señal y como se puede deducir del párrafo Anhidrita CaSO4 1.48, 1.8, 1.9,1.3
anterior la mayor parte de las mixturas existentes en la Basanita CaSO4.1/2H2O 1.46, 1.94, 2.42
naturaleza son no lineales, lamentablemente los algoritmos
Sulfatos Yeso CaSO4.2H2O 0.99, 1.19, 1.48, 1.74, 1.96
que analizan este tipo de mezclas están en un estado
embrionario. Van der Meer (2001) considera que si la escala Tenardita Na2SO4 1.43, 1.91
de la mezcla es microscópica y que si los fotones actúan
Mirabilita Na2SO4.10H2O 0.99, 1.22, 1.44, 1.94, 2.47
sobre un material en vez de actuar sobre muchos materiales,
entonces la mezcla puede ser considerada como lineal. Calcita CaCO3 1.8, 2.1, 2.33
Carbonatos
Las fuentes que contribuyen a las mezclas lineales son: Dolomita CaMg(CO3)2 2.34
a) Los sistemas de imágenes ópticas que integran la reflexión
Boratos Ulexita NaCaB5O9.8H2O 0.99, 1.46, 1.96, 2.23
de la luz sobre cada píxel.
b) Todos los materiales presentes dentro el campo de visión Tabla 2: Caracteristicas de absorción de los principlales minerales
del captor contribuyen a la reflexión de cada uno de los evaporíticos.
píxeles.
c) La variación de la luz en función de la topografía da como
resultado la atenuación de las señales.
∑ (R − R jk ') / n
n
Considerando lo anteriormente enunciado una combinación 2
jk
lineal de píxeles puros buscara a descomponer la mezcla de m
RMS = ∑
j =1
espectros de reflectancia presentes en cada píxel Ri en
fracciones fi de sus endmembers Reij y puede ser escrito k =1 m
como: donde: n es el número de bandas espectrales y m es el
n numero de píxeles dentro la imagen (Van der Meer, de Jong
Ri = ∑ f j Re ij + ε i & Bakker, 2001).
j =1

este proceso puede existir bajo las siguientes condiciones: 8.- Discusión y resultados
n
8.1.- Análisis Espectral
∑f
j =1
j =1 y 0 ≤ f ≤1
i
Antes de comenzar con el análisis de la imagen es necesario
comprender las propiedades espectrales de los minerales a
analizar. En la literatura existen numerosos trabajos que
donde Ri es la reflectancia de la mezcla espectral de cada
utilizan la espectrometría para analizar minerales
píxel en la banda i, fi es la fracción de cada píxel puro j
evaporíticos. Tal es el caso de Crowley (1991) y Drake (1995),
calculado banda por banda, Reij es la reflectancia del
Crowley (1996) que trabajan entre 0.4 a 2.5 µm y Crowley y
endmember j en el canal i, i es el numero de banda, j es cada
Hook (1996) que trabajan por encima de 2.5 µm del espectro
uno de los n endmembers de la imagen, åi es el error residual
electromagnético.
o la diferencia entre la mensura y el modelo DN (Digital
Tomando en cuenta la descripción realizada en el párrafo
Number) en la banda i. La única solución para resolver esta
geoquímica del Salar de Coipasa, los minerales que podemos
ecuación es minimizar el error residual mediante la técnica
encontrar en esta cuenca son resumidos en la tabla 2.
de mínimos cuadrados. Este error debería ser en teoría igual
La mayoría de estos minerales presenta espectros de
al ruido del instrumento, en el caso de que solamente un
absorción entre 1 a 2.5 µm, esta absorción se presenta bajo
endmember este presente en la imagen (Van der Meer, de
dos procesos distintos: la transición electrónica que toma
Jong & Bakker, 2001).
en cuenta la valencia del elemento (electrón en el ultimo
nivel orbital) y la vibración molecular (energía interatómica)
7.3.3.- Error de la Dezmezcla que implica las ligazones moleculares (Crowley, 1991). La
El error residual para cada píxel de la imagen y sobre todas mayor parte de las características espectrales observadas
las bandas puede ser dado por el error cuadrático medio o en los minerales evaporíticos tienen un origen vibratorio y
RMS (Root Mean Square), representado sobre la forma de son atribuidos generalmente a la forma de vibración interna
una imagen, el cual es calculado de la diferencia entre el de ciertos minerales, como es el caso de los boratos,
modelo (Rjk) y la mensura (Rjk’) y puede ser escrito como: carbonatos y de los aniones hidroxilos (Crowley, 1991).
 Algunos de los minerales evaporíticos presentes en la tabla
2 son hidratados, por lo que es probable que las propiedades
de absorción estan causadas por las moléculas de agua,
que son un componente esencial de la estructura de estos
minerales. La longitud de las bandas de absorción de dichos
minerales según algunos autores esta íntimamente ligada a
las O-H de las moléculas de agua que constituyen estos
minerales (Hunt et al., 1971). En las formas anhidras de estos
minerales la absorción es debido a la presencia de
inclusiones fluidas y/o la presencia de agua adsorbida
(Drake, 1995).
Para el análisis de los espectros presentes en la imagen se
utilizo la librería espectral del USGS (http://
speclab.cr.usgs.gov), esta librería construida en laboratorio
tiene una frecuencia de muestreo de ±10 nm, en los cuales
se puede observar en detalle los diferentes picos de
absorción (tabla 2) de los minerales evaporíticos. Esta
librería fue remuestreada siguiendo un modelo gausiano de
respuesta espectral y de esta manera obtener las
características espectrales del captor ASTER.
Debemos mencionar también que no todos los minerales
evaporíticos están presentes en las librerías espectrales, en
dichas librerías solamente se encuentran los principales
minerales evaporiticos, los más raros y los que determinan
en muchos casos el tipo de cuenca no están presentes.
En la figura 7 se puede observar los espectros de los
diferentes minerales evaporíticos susceptibles a estar
presentes en el Salar de Coipasa como también la posición
de las diferentes bandas del captor ASTER. De la figura
podemos deducir lo siguiente:
La halita no presenta picos de absorción, el espectro
remuestreado es casi plano, Que volverá la desmezcla de la
halite muy difícil a realizar, sin embargo se puede pensar
que ella esta presente en todo el salar en una mescla intima
con el yeso y o la mirabilita.
Los espectros de los sulfatos como es el caso de la anhidrita,
basanita y yeso tienen igual similitud, la única diferencia es
la profundidad del pico de absorción (más pronunciado para
el yeso). Esta sutil diferencia puede ser utilizada en un
análisis espectral.
La mirabilita y tenardita presentan diferentes signaturas
espectrales, sin embargo el pico de absorción de ambos
minerales se encuentra en la misma región espectral. Estos
minerales no fueron descritos en la superficie del salar y a
este remarque debemos agregar que la tenardita cristaliza
cuando la mirabilita comienza a deshidratarse. De todas
maneras en el análisis ambos minerales son tomados en
cuenta.
Los carbonatos como la calcita y la dolomita presentan picos
Figura 7: Espectros provenientes de la libreria espectral de la USGS de absorción sobre las bandas 7 y 8, las cuales se deben
comparados con los espectros remuestreados en función de las bandas
esencialmente a la forma de vibración de los aniones
del captor ASTER.
carbonatos. Lamentablemente en el salar no fueron descritos
dichas facies minerales, según la hipótesis sustentada por
Figura 8: Comparación de la imagenen ASTER (3/2/1, RVB), con la imagen resultante de la transformación MNF (primeros tres componetes). En
esta ultima imagen se puede observar las diferentes poblaciones de minerales presentes en la zona.

Lebrun, et al (2002) dichas facies cristalizan en los ríos La figura 9 nos muestra los diferentes endmembers
tributarios del salar y que dichos minerales pueden ser encontrados en la imagen, los nombres en la figura son el
recubiertos por yeso o halita. resultado de la comparación con la librería espectral teórica.
Por ultimo la ulexita presenta un pico de absorción en la La variación entre la forma del espectro del yeso-1 y yeso-
banda 5, para que este mineral se forme tiene que existir un 2, puede tener su origen en muchos factores físicos y
aporte de boro en la cuenca (generalmente una fuente químicos, como el contenido de agua, tamaño de las
hidrotermal). Este mineral no fue descrito en la superficie partículas, cambios en la composición química del mineral
del salar, las únicas ocurrencias que existen en la zona se evaporítico, etc. Como el sensor capta la información a más
encuentran al Sud - Oeste del salar en la zona de Hizo. de 700 km de altitud, estos parámetros son difíciles a
Si comparamos nuestra deducción a la información de la mensurar, por lo que preferimos utilizar los dos espectros
tabla 2, podemos observar que muchos de los picos de de referencia en el análisis de la imagen. Una verificación de
absorción característicos de estos minerales están en la terreno confirmara o negara la existencia de dichas
zona de absorción del agua lo que imposibilita ser detectados poblaciones.
con el captor ASTER.
8.3.- Desmezcla Espectral
8.2.- Valores extremos En la figura 10 se puede observar los resultados de la
Siguiendo todo el proceso descrito en el párrafo 7.3.1, un desmezcla espectral como los errores proporcionados por
primer parámetro de cuantificación es la imagen MNF en el método. El paso siguiente es la verificación de terreno,
falso color. Esta imagen representa más del 80% de la como este es un trabajo preliminar por el momento no existe
información contenida en toda la imagen. En esta escena un trabajo de verificación.
podemos observar también las diferentes poblaciones que De todas maneras en los resultados se puede observar la
existen en la zona. proporción de cada mineral en la imagen y claramente se
Una comparación de la imagen MNF con la imagen ASTER pueden diferenciar tres zonas:
(bandas 3/2/1 en RVB), puede verse en la figura 8. El yeso-1 tiene dos poblaciones, una al Este y otra al Oeste
Después del análisis PPI fueron retenidas solamente 3 de la imagen, la población del este coincide con las
endmembers que tienen algunas características espectrales descripciones del proyecto Coipasa (Lebrun, 1997).
similares a los minerales de la librería espectral de la USGS.
Figura 10: Imagenes resultantes de la desmezcla espectral. a) Distribución de la abundancia de mirabilita, b) yeso-1, c) yeso-2 y d) error del método.
 9.- Conclusiones
Como dijimos en la introducción de este trabajo y lo
volvemos a recalcar en este acápite, el motivo del presente
trabajo es mostrar a la comunidad científica boliviana la
existencia de otros métodos capaces de diferenciar minerales
a un bajo costo y con una precisión aceptable.El captor
ASTER con la información espectral que mensura es capaz
de diferenciar los minerales secundarios presentes en
Coipasa.
El advenimiento de la tecnología hyperespectral embarcada
sobre los satélites, como ya es el caso de Hyperion, facilitara
este trabajo y en un futuro muy cercano las cartas
“geoquímicas” obtenidas a partir de los satélites serán muy
comunes y de nuestra adaptación a este nuevo tipo de
técnicas dependerá en cierta medida el futuro de la
comunidad científica nacional.

10.- Agradecimientos
Los autores agradecen a las personas que participaron en
el proyecto Coipasa: Sophie Fis, André Lannoye, Charles-
Etienne Housse, Pierre Francus, Guy Wansard y Pedro
Pacosillo.

11.- Referencias
Alpers, N. C. Whittemore, D. O. 1990. Hydrogeochemistry
and stable isotopes of ground and surface waters from two
adjacent closed basins, Atacama Desert, Nothern, Chile.
Applied Geochemistry, 5: 719-734.
Ballivián, O. & Risacher, F. 1981. Los salares del Altiplano
Boliviano. Métodos de estudio y estimación económica.
UMSA – ORSTOM, 241p.
Figura 9: Espectros de los endmembers resultantes del análisis,
comparados con espectros de laboratorio (USGS) remuestreados en Boardman, J. W. Kruse, F.A. & Green, R. O. 1995. Mapping
función de los canales de ASTER. target signatures via partial unmixing of AVIRIS data. In
Summaries of 5th JPL Airborne Earth Science Workshop,
Jet Propulsion Laboratory, Pasadena, California. Publication
95-1: 23-26.
La zona de la miribilita se encuentra en la parte Norte de la Bryant. R. G. 1996. Validated linear mixture modelling of
imagen y en una población menos intensa en la parte Este Landsat TM data for mapping evaporite minerals on a playa
de la imagen. Este mineral en su forma cristalina fue surface: methods and applications. International Journal of
encontrado en profundidad por el proyecto Coipasa y la Remote Sensing, 17 (2): 315-330.
ocurrecia de dicha poblacion necesita ser verificada sobre Bryant, R. G. 2003. Monitoring playa systems using remote
terreno. sensing: methods and applications. AridLands. 51.
Cabe recalcar que la mayor parte de las descripciones (accessible via: http://ag.arizona.edu/OALS/ALN/aln51).
superficiales fueron realizadas con una lupa de geólogo y Camargo, R. A. Ponzoni, F. J. Zullo, Jr. J. Queiroz, G. &
diferenciar un agregado amorfo de mirabilita de un agregado Arnaud, Y. 2003. Characterization of the Salar de Uyuni for
amorfo de halita es un trabajo muy difícil de efectuar. in-orbit satellite calibration. IEEE Transactions on
La tercera zona a la que denominamos yeso-2 se encuentra Geoscience and Remote Sensing, 41 (6): 1461-1468.
entre las dos zonas anteriormente descritas. Chapman, J. E. Rothery, D. A. Francis, P. W. & Pontual, A.
El error del método se encuentra al Sud Oeste de la imagen, 1989. Remote Sensing of evaporite mineral zonation in salt
esta zona no tiene ninguna relación con las dos primeras flats (salars). International Journal of Remote Sensing, 10
zonas. (1): 245-255.
Crowley, J. K. 1991. Visible and near-infrared (0.4 – 2.5 um) Hardie, L. & Eugster, H. 1970. The evolutions of closed-
reflectance spectra of playa evaporite minerals. Journal of basin brines. Mineral Society American Special Paper, 3:
Geophysical Research, 96 (B10): 16231-16240. 273-290.
Crowley, J. K. 1993. Mapping playa evaporite minerals with Houssa, C.-E. 1998. Analyse mineralogique des ressources
AVIRIS data: a first report from Death Valley, California. du Salar de Coipasa, Bolivie. Travail de fin d’études
Remote Sensing of Environment, 44: 337-356. d’Ingénieur Civil Géologue, Université de Liège, Belgique,
Crowley, J.K. 1996. Mg- and K-Bearing borates and 130p
associated evaporites at Eagle Borax Spring, Death Valley, Hunt, G.R. Salisbury, J.W. & Lenhoff, C.J. 1971. Visible and
California: a spectroscopic exploration. Economic Geology, near-infrared spectra of minerals and rocks, IV, Sulfides and
91: 622-635. Sulfates. Modern Geology, 3: 1-14.
Crowley, J. K. & Clark, R. N. 1992. AVIRIS study of Death Iwasaki, A. Fujisada, H. Asao, H. Shindou, O. & Akagi, S.
Valley evaporite deposits using least-squares band-fitting 2002. Enhancement of spectral separation Performance for
methods. In: Summaries of the 3th Annual JPL Airborne ASTER/SWIR, Infrared Spaceborne Remote Sensing. IX,
Geoscience Worshop, Jet Propulsion Laboraty, Pasadena, Proceedings SPIE Vol: 4486, Paper-No, 4486-06: 42-50.
California. JPL Publication 92-14, 1: 29-31. Keshave, N. & Mustard. J. F. 2002. Spectral unmixing. IEEE
Crowley, J. K. & Hook, S. J. 1996. Mapping playa evaporite Signal Procesing Magazine,19 : 44-57.
minerals and associated sediments in Death Valley, Lebrun, V. 1997. Etude géochimique du Salar de Coipasa,
California, with multispectral thermal infrared images. Journal Bolivie. Travail de fin d’études d’Ingénieur Civil Géologue,
of Geophysical Research, 101 (B1): 643-660. Université de Liège, Belgique, 140p.
Crowley, J. K. & Mars, J. C. 2000. Mapping evaporite Lebrun, V. Pacosillo, P. Gutierrez, J. Cáceres, F. Pool, E. &
minerals in the Death Valley salt pan using MODIS/ASTER Pirard, E. 2002. Geochemistry of bitter brines in the Salar de
airborne simulator (MASTER) data. Fourteenth International Coipasa (Bolivia). Aardkundige Mededelingen, 12: 21-24.
Conference on Applied Geologic Remote Sensing, Las Vegas, Rettig, S. L. Davis, J.R. Smith, R. L. Smith, K. A. & Ericksen,
Nevada,USA: 344-345. G. E. 1979. Bolivia. U.S. Geological Survey Professional
Drake, N. A. 1995. Reflectance spectra of evaporite minerals Paper, 321 p.
(400 – 2500 nm): applications for remote sensing. Risacher, F. 1992. Géochimie des lacs sales et croûtes de sel
International Journal of Remote Sensing, 16 (14): 2555-2571. de l’Altiplano bolivien. Sciences Géologiques Bulletin, 45
Ericksen, G. E. Vine, J. D. & Ballon, R. 1978. Chemical (3-4): 219p.
composition and distribution of lithium-rich brines in salar Risacher, F. & Fritz, B. 1991. Quaternary geochemical
de Uyuni and nearby salars in Southwester Bolivia. Energy, evolution of the salars of Uyuni and Coipasa, Central
3: 355-363. Altiplano, Bolivia. Chemical Geology, 90: 211-231.
Eugster, H. & Hardie, L. 1978. Saline Lakes. In Lakes, Risacher, F. & Fritz, B. 1992. Mise en evidence d’une phase
Chemistry, Geology, Physics. A. Lehnman ed. Springer- climatique holocène extrêmement aride dans l’Altiplano
Verlag, 237-293. Central par la présence de la polyhalite dans le Salar de
Fis, S. 1996. Contribution de la prospection géophysique a Uyuni (Bolivie). Comptes Rendus de l’Académie des
l’étude du Salar de Coipasa (Oruro, Bolivie). Travail de fin Sciences, Paris, 314 (t. II) : 1371-1377.
d’études d’Ingénieur Civil Géologue, Université de Liège, Pacosillo, P. 1998 Estudio preliminar de la cuenca evaporítica
Belgique, 151p. del Salar de Coipasa, en base al muestreo de ríos, fuentes y
Fornari, M. Risacher, F. & Féraud. G. 2001. Dating of salmueras. Tesis de Grado Carrera de Ingeniería Geológica
paleolakes in the central Altiplano of Bolivia. FNI – UTO.
Palaeogeography, Palaeoclimatology, Palaeoecology, 172: Plaza, A. Martinez, P. Perez, R. M. & Plaza, J. 2002. Spatial/
269-282. spectral endmember extraction by multidimensional
Fujisada, H. 1996. ASTER Level-1 data processing algorithm. morphological operations. IEEE Transactions on Geoscience
IEEE Transactions on Geoscience and Remote Sensing, 36 and Remote Sensing, 40 (9): 2025-2041.
(4): 1101-1112. Pool, E. 1995. Etude préliminaire du Salar de Coipasa (Oruro-
Green, A. A. Berman, M. Switzer, M. & Craig, M. D. 1988. A Bolivia) sur base des données géologiques, géochimiques
transformation for ordering multispectral data in terms of et minéralogiques. Mémoire de Maîtrise en Géologie des
image quality with implications for noise removal. IEEE Terrains Superficiels, Université de Liège, Belgique, 136p.
Transactions on Geoscience and Remote Sensing, 26 (31): Sabins, F. F. 1997. Remote sensing, principles and
65-74. interpretation. W. H. Freeman & Company, New York, 3rd.
Gutierrez, J. 1997. Aplicación del G.P.S. y la gravimetría para Ed.
un estudio preliminar del paleo-relieve interno del Salar de Sabins, F. F. 1999. Remote sensing for mineral exploration.
Coipasa. Tesis de Grado Carrera de Ingeniería Geológica Ore Geology Reviews, 14: 157-183.
FNI – UTO.
USGS – Geobol. 1992. Geology and Mineral Resources of Van der Meer. F. 2001. Basic physics of spectrometry. In
the Altiplano and Cordillera Occidental, Bolivia. U.S. 195 - Imaging spectrometry: basic principles and prospective
209 Geological Survey Bulletin, 1975: applications. F. Van der Meer & S. de Jong, ed. Kluwer
Sabins F. F. & Miller , R. M. 1994. Resource assessment – Academic Publishers, 3-16.
Salar de Uyuni and vinicity. Tenth Thematic Conference on Van der Meer, F. de Jong, S. & Bakker, W. 2001. Imaging
Geologic Remote Sensing, Vol. I: I92-I103. spectrometry: basic analytic techniques. In Imaging
Servant, M. & Fontes, J. C. 1978. Les lacs quaternaires des spectrometry: basic principles and prospective applications.
plateaux des Andes boliviennes premières interprétations F. Van der Meer & S. de Jong, ed. Kluwer Academic
paléoclimatiques. Cahiers ORSTOM, Série Géologie, 10 (1): Publishers, 16-61.
9-23. White, K. & Drake, N. 1993. Mapping the distribution and
Sylvestre, F. Servant-Vildary, S. Servant. M. 1998. Le dernier abundance of gypsum in South-Central Tunisia from
maximum glaciaire (21 000-17000 14C ans B. P.) dans les Landsat Thematic Mapper data. Zeitschrift fur
Andes tropicales de Bolivie d’après l’étude des diatomées. Geomorphologie, 37 (3): 309-325.
Comptes Rendus de l’Académie des Sciences, Paris, 327 : Yamaguchi, Y. Kawkami, T. Kahle, A. B. Pniel, M. & Tsu, H.
611-618. 1998. ASTER mission planning and operations concept.
Sylvestre, F. Servant, M. Servant-Vildary, S. Causse, C. SpaceOps, International Symposium on Technical
Fournier, M. & Ybert, J.-P. 1999. Lake-level chronology on Interchange for Space Mission Operations and Ground Data
the southern Bolivian Altiplano (18°-23°S) during late- Systems: 8p.
glacial time and the Early Holocene. Quaternary Research,
51 (1): 54-66.