Sintesis Geologica y Cartografia Preliminar de Los Salares Utilizando Imagenes Aster
Sintesis Geologica y Cartografia Preliminar de Los Salares Utilizando Imagenes Aster
Sintesis Geologica y Cartografia Preliminar de Los Salares Utilizando Imagenes Aster
Fernando Cceres Jerez Eric Pirard Universit de Lige Dpartement des Goressources, Gotechnologies et Matriaux de Construction (GEOMAC) Sart Tilman B52/3 4000 Lige Belgique Tel: +32 43669525 Fax: +32 43669520 [email protected] Jorge Gutirrez Venegas Universidad Tcnica de Oruro Facultad Nacional de Ingeniera Carrera de Geologa Centro de Investigaciones en Teledeteccin (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 utilizacin de la tcnicas de teledeteccin para la cartografa de minerales evaporticos fue grandemente utilizado. En este trabajo se utilizan las imgenes VNIR y SWIR del captor ASTER con el objetivo de cartografiar los minerales secundarios presentes en el Salar de Coipasa. Mediante las tcnicas de desmezcla espectral aparte de identificar minerales presentes en la imagen, permite tambin cartografiarlos. En la imagen subset de Coipasa se pudo constatar la presencia de dos minerales evaporticos identificados en base a su firma espectral (yeso y mirabilita).La distribucin del yeso coincide con las descripciones de terreno realizadas durante un programa de exploracin.
Palabras claves: Salar de Coipasa, Sensores Remotos, Minerales Evaporiticos, Desmezcla Espectral.
1.- Introduccin
Este articulo es una contribucin ms, a la larga literatura que existe sobre el estudio de secuencias evaporticas utilizando sensores remotos. Los pticos (multi hiperespectrales) y radares demostraron su potencial en la cartografa y monitoreo de este tipo de ambiente. A esta experiencia podemos adjuntar la facilidad que tienen estas tcnicas de cubrir grandes regiones con un bajo costo (50 centavos de $us americanos/km2). Los resultados presentados en este trabajo tienen un carcter preliminar, los cuales estn siendo confirmados y extrapolados a todo el salar con la elaboracin de una tesis realizada por la egresada de la Carrera de Ingeniera Geolgica, Luca Juregui S.
contenido de boro en el salar de Uyuni, con este propsito los autores utilizan la relacin TM 4/TM 7. Bryant (1996) utiliza la tcnica de desmezcla espectral sobre imgenes Landsat TM de la regin Chott el Djerid (Tunes) con el objetivo de determinar la distribucin espacial de yeso, halita y carnalita. Bryant (2003) emplea imgenes AVHRR en el monitoreo de playas en el Norte y Sud del continente africano. En este tipo de ambiente los cambios hidrolgicos pueden ser detectados por el anlisis de imgenes satelitales, estas variaciones permiten determinar de manera indirecta los cambios de precipitacin y evaporacin del sistema. Camargo, et al. (2003) emplea la informacin de reflectancia adquirida in situ sobre el salar de Uyuni con el fin de calibrar las imgenes satelitales, los autores recomiendan la utilizacin de este tipo de ambientes por su elevada homogeneidad espectral y los dbiles efectos atmosfricos que existen en la zona. Otros autores emplean captores aerotransportados para el estudio de este tipo de ambiente, como es el caso de la regin de Death Valley (California). Con este propsito Crowley & Clark (1992) y Crowley (1993) utilizan imgenes que provienen del captor AVIRIS (Airborne Visible/Infrared Imaging Spectrometer) y Crowley, Mars & Hook (2000) emplea imgenes MASTER (MODIS/ ASTER Airborne Simulator).
La parte Norte del salar esta cubierto por un lago permanente alimentado por el ro Lauca, principal tributario del salar. En la poca de lluvias (noviembre a marzo) todo el salar es recubierto por una capa de salmuera de 15 cm de alto. Hacia el Sud del Salar de Coipasa se encuentra la costra de sal ms grande del planeta, el Salar de Uyuni (10.000 km2), su principal tributario el Ro Grande desemboca en la parte Sud del salar. El Sud del Altiplano, geogrficamente conocido como los Lipez, presenta numerosos lagos salados y salmueras (ms de 30) en el seno de pequeas cuencas de origen volcnico. El quimismo de las salmueras de esta regin es variado. En la figura 1 se puede observar la localizacin de las diferentes cuencas del Altiplano.
Figura 2: Seccin N-S del Altiplano Central para observar la relacin que existe entre los lagos y salares (modificado de Risacher, 1992)
El lago Titicaca alcanzo su nivel actual hace 10.000 aos atrs, en la fase rida del Holoceno (entre 7700 y 3650 aos atrs), su nivel bajo en 50 m. Es importante remarcar que esta fase fue puesta en evidencia en otras regiones de Amrica Latina. Risacher & Fritz (1992) igualmente reconocen esta fase rida en el Salar de Uyuni gracias a la presencia de polyhalita, mineral que precipita en salmueras concentradas cinco veces ms que las actuales. Despus de esta fase rida el lago Titicaca alcanzo su nivel actual de 3808 m. Dos niveles de endurecimientos alargados localizados entre 3700 y 3747 m tapizan toda la cuenca Pocoyu (denominacin proporcionada a la cuenca Poop - Coipasa - Uyuni por Wirrmann & Mourguiart, 1995 y Argollo & Mourguiart, 2000). En funcin a esta informacin pudieron ser puestos en evidencia dos palelagos en la parte Central del Altiplano, las fases fueron datadas 27.000 aos para el lago Minchin y 13.000 a 10.000 aos para el Tauca (figura 3). Fornari, Risacher & Fraud (2001) luego de la datacin de los 121 metros perforados en el salar de Uyuni confirman la existencia del evento Escara (191.000 aos atrs), una fase que no es tomada en cuenta por muchos autores. Las dataciones realizadas sobre las construcciones arreciferas son errticas a causa de la presencia de fenmenos de recristalizacin, sin embargo la correlacin de los mencionados palelagos con las costras algareas pudo claramente ser establecida. El estudio de las diatomeas como la presencia de Coprolitos de Artemia (pequeos 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 sequa hace 10.000 aos atrs 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 aos atrs, los ltimos trabajos realizados por Fornari, Risacher & Fraud (2001) determinan que el evento Coipasa es una oscilacin del gran lago Tauca.
El estudio de los fenmenos evaporticos esta en relacin a este ltimo punto; a grandes rasgos los mecanismos de saturacin que pueden conducir a las variaciones de concentracin de las aguas naturales son:
En este caso la constante del producto [CO3--][Ca++] [Ca++][HCO3-]2 ser la concentracin simultanea de calcio y de iones carbonatos. En el curso de la precipitacin de la calcita, el elemento inicial en exceso se concentrara y el elemento inicial en dficit se diluye. En la prctica las aguas que se encuentran en la naturaleza contienen a la vez carbonatos y bicarbonatos, por lo que es mejor trabajar con el ndice de alcalinidad (alc): alc = 2[CO3--] + [HCO3-] Si el alc es superior a 2[Ca++] en la solucin inicial, ella evolucionara hacia la zona alcalina, se concentrara progresivamente en lcalis y su pH se establecer entre 8 y 9. Si alc es inferior a 2[Ca++], entonces la solucin seguir el camino de la salinidad neutra, ella se enriquecer en calcio y su pH se estabilizara entre 7 8. Dentro el camino neutro el prximo mineral a precipitar ser el yeso, CaSO 4.2H 2O, la reaccin que controla esta precipitacin puede escribirse como: [Ca++][SO4--][H2O] = K3 Donde, K3 es el producto de solubilidad del yeso. La solucin se concentra en sulfatos o en calcio segn el origen de la precipitacin; si existe un exceso de uno u otro elemento, entonces la salmuera evolucionara hacia una fase de tipo Na-Ca-(Mg)-Cl o Na-(Mg)-Cl-SO4.
Figura 4: Diagrama terico seguido por las aguas en un fenmeno evaportico (Hardie & Eugster, 1970 y Risacher, 1992).
Dentro el camino alcalino el dficit de calcio limita la formacin de yeso y conduce a la formacin de silicatos magnsicos segn la reaccin: Mg++ + nH4SiO4 => Silicato Mg + H+ Los iones H+ liberados por la reaccin neutralizaran una parte de los iones carbonatos. Lo cual puede invertir la evolucin de la solucin hacia un camino salino neutro. La condicin para la existencia de esta inversin es de tener en el agua inicial: alc < 2[Ca++] + 2[Mg++]. Si esta reaccin no se verifica, entonces la solucin contina dentro del camino alcalino. Risacher & Fritz (1991) constatan que en un buen nmero de aguas del Altiplano no se verifica esta condicin, sin embargo evolucionan hacia una facie salina neutra. Esto quiere decir que la composicin qumica de partida no es una condicin suficiente para explicar el futuro qumico del agua. Otra hiptesis pone de manifiesto un segundo factor que influye en la evolucin de la salmuera, dicho componente esta en relacin con la interaccin del agua con los sedimentos. Esta teora esta basada en las observaciones efectuadas en los salares tipo carbonato sdico que tienen pocos sedimentos lacustres debajo la costra de sal.
Por el contrario los salares que bifurcan hacia condiciones de salindad neutra tienen gran estrato de sedimentos arcillosos, arcillos arenosos y salferos que encierran la napa subaflorante. Entonces, cuando las aguas de la fuente o del ri penetran dentro los sedimentos localizados sobre el borde del salar, algunos de los siguientes fenmenos pueden producirse: a) Disolucin del yeso: el yeso es siempre abundante en el borde del salar y tiene como origen la sequedad del lago precedente o la evaporacin por capilaridad de la napa ms profunda. Las aguas que entran en la cuenca adicionan Ca++ en la solucin, lo que hace disminuir la relacin: alc/ (2[Ca++] + 2[Mg++]) b) La transformacin del vidrio volcnico en arcilla: las aguas de aporte tienen en suspensin vidrio volcnico proveniente de las formaciones volcnicas de la cuenca. Este vidrio volcnico dentro del salar se transforma en arcilla con la liberacin de iones H+ que tienen la tendencia de neutralizar los carbonatos.
c) La composicin de las rocas de la cuenca: si estas rocas son cristalizadas, volcnicas, sin impurezas ni mineralizaciones, ellas no proporcionaran ningn in en cantidades apreciables a la cuenca. Entonces el anin mayor ser HCO3- de la atmsfera, el cual conduce sin remedio hacia el polo alcalino. Si por lo contrario en la cuenca existe mineralizacin o rocas sedimentarias, ellas liberan sulfatos y cloruros que empujan a las salmueras hacia una salinidad neutra, como es el caso de Coipasa.
En la poca de lluvia, el agua concentrada del lago inunda toda la costra de sal y los iones son dispersados en todas partes por difusin. En el Salar de Coipasa los ros tributarios tienen el papel de homogenizar las salmueras cuando la costra esta completamente inundada. Esta hiptesis explica mejor el enriquecimiento que existe en salmueras en la parte SE del salar y su empobrecimiento en la parte Oeste (zona sin grandes aportes de agua). Lebrun et al. (2002) pone de manifiesto una elevada diferencia de la relacin Mg/Li en el Salar de Coipasa con respecto al Salar de Uyuni. Esta diferencia tiene su origen en la composicin y proporcin de las rocas que se encuentran dentro la cuenca. El Li, Mg y B esencialmente provienen de las rocas volcnicas; por otro lado el Mg y K son abundantes en otro tipo de rocas. Si la proporcin de rocas volcnicas o su composicin es diferente de una cuenca a otra, esto inducir en la diferencia de la relacin Mg/Li en el agua de la cuenca. El relativo enriquecimiento de sulfatos en las salmueras de Coipasa esta directamente relacionado a las concentraciones de sus tributarios. El Ro Lauca tiene una mayor concentracin de sulfatos en comparacin al Ro Grande. El origen de este enriquecimiento de dichos sulfatos es difcil a verificar, estos pueden estar relacionados al viento, el cual transporta los sulfuros de la actividad volcnica de la zona y luego los deposita en el drenaje de la cuenca (Risacher & Fritz, 1991). Otra posible explicacin, es la presencia de sulfuros en la cuenca, como fue descrito en el salar de Hamburgo en Chile (Alpers & Whittemore, 1990).
Subs ys te ma
Banda N o.
Nive l de Quantificacin
1 2
VNIR
3N 3B 4 5 6 0.78- 0.86 0.78- 0.86 1.60- 1.70 2.145- 2.185 2 . 18 5 - 2 . 2 2 5
15
8 bits
60
SWIR
7 8 9 10 11 2.235- 2.285 2.295- 2.430 2.360- 2.430 8.125- 8.475 8.475- 8.825 8.925- 9.275 10.25- 10.95 10.95- 11.65
30
8 bits
60
TIR
12 13 14
90
16 bits
60
6.- Datos
El trabajo esta particularmente basado en las imgenes del captor ASTER (Advanced Spaceborne Thermal Emissin and Reflection Radiometer) que es uno de los cinco sensores embarcados sobre el satlite EOS-AM1 y denominado despus TERRA plataforma, puesto en orbita en 1998 y que remplaza al satlite EOS de la NASA. Su orbita a 700 km de altura, heliosincronica, cercana a los Polos y el tiempo que pasa a nivel del Ecuador permiten minimizar los efectos perjudiciales causados por las nubes. El satlite tiene una frecuencia de repeticin de 16 das, es decir que despus de 16 das toma una imagen de la misma regin. El captor ASTER fue construido por el Ministry of Internacional Trade and Industry (MITI) del Japn y esta compuesto de tres subsistemas: El primer subsistema VNIR (Visible and Near Infrared) cubre la regin visible e infrarrojo cercano del espectro electromagntico. El subsistema es un escner tipo pushbroom formado de dos telescopios uno a nadir que tiene tres bandas y otro localizado en la parte posterior con una simple banda; este canal es utilizado para las observaciones estereoscpicas. Los detectores CCD son lineales construidos en Si, cada captor posee 5000 detectores, los detalles de la resolucin espectral y espacial del subsistema son mostrados en la tabla 1. En segundo subsistema SWIR (Shortwave Infrared), con sus 6 bandas cubre el infrarrojo cercano y medio, los detectores son tambin de tipo pushbroom y opera con un simple telescopio que apunta a nadir, los datos son adquiridos por un espejo. Los detectores CCD son lineales y construidos en PtSi, cada banda espectral posee 2048 detectores. En la tabla 1 se observa la resolucin espectral y espacial del subsistema. El tercer subsistema TIR (Thermal Infrared) con sus 5 bandas mensura la parte termal del infrarrojo; es un escner mecnico compuesto de un simple telescopio a posicin fija observando a nadir y de un escner con un espejo localizado en la parte posterior. Cada banda utiliza 10 detectores PC construidos en HgCdTe, para preservar la dinmica de la imagen ella es codificada en 16 bits. La informacin espacial y espectral se puede observar en la tabla 1. Los datos que adquiere el satlite tienen una amplia cobertura espectral y una resolucin relativamente elevada, lo que permite discriminar una serie de materiales sobre la superficie (Yamaguchi et al. 1998). En este trabajo utilizamos una imagen adquirida el 19/03/ 2001 a 7 h 07 m 15 s (hora universal), el nivel de correccin es 1B y es una de las cuatro imgenes que cubre enteramente el salar.
En este trabajo solo se toma en cuenta las imgenes de los subsistemas VNIR y SWIR, no se considera el subsistema TIR ya que las mensuras dentro este subsistema y segn la ley de Boltzmann toman en cuenta tres propiedades fsicas del material (absorcin, reflexin y temperatura), por lo contrario los otros dos subsistemas solamente mensuran uno (reflexin). Con el fin de realizar un ensaye del mtodo presentado en este artculo, todo el trabajo es realizado sobre un subset de la imagen, el cual cubre un rea aproximada de 40 km2. La localizacin del subset con respecto al salar puede observarse en la figura 5. La zona fue elegida por la variedad mineralgica que se observa en superficie; en el punto H4 ; V-5 (19 25.347 latitud sud y 68 20.613 longitud W) el proyecto Coipasa aparte de la halita encontro tambin cristales de yeso.
La zona estudiada presenta una gran uniformidad espectral y mineralgica, nosotros pensamos que no existe una creacin de artefactos. El mtodo puede ser utilizado en regiones donde no existe informacin meteorolgica, como es el caso de Coipasa. Los modelos ms complejos proporcionan mejores resultados pero desgraciadamente requieren de mayor informacin; debemos mencionar tambin que la informacin espectroscpica adquirida sobre el terreno puede servir para eliminar el efecto atmosfrico.
En la literatura existen muchos mtodos para extraer los pxeles ms puros de una imagen, en este trabajo se utiliza un mtodo basado en la informacin espectral. Boardman et al., (1995) sugiere la tcnica PPI (Pxel Purity Index), la cual se basa en las propiedades geomtricas de un conjunto convexo y considera a cada pxel como un vector en un espacio a n-dimensiones (n-D). El algoritmo no es aplicado a la imagen original sino a la imagen resultante de una reduccin dimensional de datos basada en el mtodo MNF (Minimun Noise Fraction) desarrollada por Green, et al., (1988) El algoritmo comienza con la generacin de un vector aleatorio a n-D (denominado skewer), en el cual todos los pxeles de la imagen son proyectados, los pxeles que aparecen en la parte extrema del skewers son marcados y luego contados, este procedimiento es generalizado a un gran numero de skewer generados de igual manera. Al final de este proceso se obtiene una imagen con los pxeles ms contados, es decir los pxeles que fueron contados ms de 1000 veces son retenidos como posibles endmembers. Los candidatos a ser considerados como verdaderos endmembers son seleccionados en base de una comparacin con una librera espectral y explorando su posicin en el espacio n-D. Esta seleccin supervisada requiere un gran conocimiento de las tcnicas de teledeteccin y espectroscopia mineral. El mayor defecto de este mtodo es la pobre interactividad y los resultados no pueden ser empleados para analizar un gran nmero de imgenes. Muchos autores consideran que la naturaleza aleatoria de los skewers es otra gran limitacin del algoritmo (Plaza, et al., 2002). Todo el procedimiento seguido para la obtencin de los valores extremos se puede observar en la figura 6.
d) La mezcla molecular, acaece a nivel atmico, por ejemplo dentro de dos lquidos o entre lquido y slido, en este caso la mezcla es tambin considerada como no lineal. El tipo y la naturaleza de las mezclas son cruciales para el conocimiento de la seal y como se puede deducir del prrafo anterior la mayor parte de las mixturas existentes en la naturaleza son no lineales, lamentablemente los algoritmos que analizan este tipo de mezclas estn en un estado embrionario. Van der Meer (2001) considera que si la escala de la mezcla es microscpica y que si los fotones actan sobre un material en vez de actuar sobre muchos materiales, entonces la mezcla puede ser considerada como lineal. Las fuentes que contribuyen a las mezclas lineales son: a) Los sistemas de imgenes pticas que integran la reflexin de la luz sobre cada pxel. b) Todos los materiales presentes dentro el campo de visin del captor contribuyen a la reflexin de cada uno de los pxeles. c) La variacin de la luz en funcin de la topografa da como resultado la atenuacin de las seales. Considerando lo anteriormente enunciado una combinacin lineal de pxeles puros buscara a descomponer la mezcla de espectros de reflectancia presentes en cada pxel Ri en fracciones fi de sus endmembers Reij y puede ser escrito como:
Clas e Cloruros
M ine ral
Formula Qumica
Picos de Abs orcion (m) 1.93 1.48, 1.8, 1.9,1.3 1.46, 1.94, 2.42 0.99, 1.19, 1.48, 1.74, 1.96 1.43, 1.91 0.99, 1.22, 1.44, 1.94, 2.47 1.8, 2.1, 2.33 2.34
Sulfatos
Carbonatos
Dolomita
Boratos
Ulexita
RMS =
k =1
(R
n j =1
jk
R jk ') / n
2
Ri = f j Re ij + i
j =1
donde: n es el nmero de bandas espectrales y m es el numero de pxeles dentro la imagen (Van der Meer, de Jong & Bakker, 2001).
f
j =1
=1 y 0 f 1 i
donde Ri es la reflectancia de la mezcla espectral de cada pxel en la banda i, fi es la fraccin de cada pxel puro j calculado banda por banda, Reij es la reflectancia del endmember j en el canal i, i es el numero de banda, j es cada uno de los n endmembers de la imagen, i es el error residual o la diferencia entre la mensura y el modelo DN (Digital Number) en la banda i. La nica solucin para resolver esta ecuacin es minimizar el error residual mediante la tcnica de mnimos cuadrados. Este error debera ser en teora igual al ruido del instrumento, en el caso de que solamente un endmember este presente en la imagen (Van der Meer, de Jong & Bakker, 2001).
Figura 7: Espectros provenientes de la libreria espectral de la USGS comparados con los espectros remuestreados en funcin de las bandas del captor ASTER.
Algunos de los minerales evaporticos presentes en la tabla 2 son hidratados, por lo que es probable que las propiedades de absorcin estan causadas por las molculas de agua, que son un componente esencial de la estructura de estos minerales. La longitud de las bandas de absorcin de dichos minerales segn algunos autores esta ntimamente ligada a las O-H de las molculas de agua que constituyen estos minerales (Hunt et al., 1971). En las formas anhidras de estos minerales la absorcin es debido a la presencia de inclusiones fluidas y/o la presencia de agua adsorbida (Drake, 1995). Para el anlisis de los espectros presentes en la imagen se utilizo la librera espectral del USGS (http:// speclab.cr.usgs.gov), esta librera construida en laboratorio tiene una frecuencia de muestreo de 10 nm, en los cuales se puede observar en detalle los diferentes picos de absorcin (tabla 2) de los minerales evaporticos. Esta librera fue remuestreada siguiendo un modelo gausiano de respuesta espectral y de esta manera obtener las caractersticas espectrales del captor ASTER. Debemos mencionar tambin que no todos los minerales evaporticos estn presentes en las libreras espectrales, en dichas libreras solamente se encuentran los principales minerales evaporiticos, los ms raros y los que determinan en muchos casos el tipo de cuenca no estn presentes. En la figura 7 se puede observar los espectros de los diferentes minerales evaporticos susceptibles a estar presentes en el Salar de Coipasa como tambin la posicin de las diferentes bandas del captor ASTER. De la figura podemos deducir lo siguiente: La halita no presenta picos de absorcin, el espectro remuestreado es casi plano, esto hara que la desmezcla de la halite sea difcil de realizar, sin embargo se puede pensar que ella esta presente en todo el salar en una mezcla 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 absorcin (ms pronunciado para el yeso). Esta sutil diferencia puede ser utilizada en un anlisis espectral. La mirabilita y tenardita presentan diferentes signaturas espectrales, sin embargo el pico de absorcin de ambos minerales se encuentra en la misma regin 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 anlisis ambos minerales son tomados en cuenta. Los carbonatos como la calcita y la dolomita presentan picos de absorcin sobre las bandas 7 y 8, las cuales se deben esencialmente a la forma de vibracin de los aniones carbonatos. Lamentablemente en el salar no fueron descritos dichas facies minerales, segn la hiptesis sustentada por
Figura 8: Comparacin de la imagenen ASTER (3/2/1, RVB), con la imagen resultante de la transformacin 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 ros tributarios del salar, y as estos minerales pueden ser recubiertos por yeso o halita. Por ultimo la ulexita presenta un pico de absorcin en la banda 5, para que este mineral se forme tiene que existir un aporte de boro en la cuenca (generalmente una fuente hidrotermal). Este mineral no fue descrito en la superficie del salar, las nicas ocurrencias que existen en la zona se encuentran al Sud - Oeste del salar en la zona de Hizo. Si comparamos la deduccin con la informacin de la tabla 2, podemos observar que muchos de los picos de absorcin caractersticos de estos minerales estn en la zona de absorcin del agua lo que imposibilita ser detectados con el captor ASTER.
La figura 9 nos muestra los diferentes endmembers encontrados en la imagen, los nombres en la figura son el resultado de la comparacin con la librera espectral terica. La variacin entre la forma del espectro del yeso-1 y yeso2, puede tener su origen en muchos factores fsicos y qumicos, como el contenido de agua, tamao de las partculas, cambios en la composicin qumica del mineral evaportico, etc. Como el sensor capta la informacin a ms de 700 km de altitud, estos parmetros son difciles a mensurar, por lo que preferimos utilizar los dos espectros de referencia en el anlisis de la imagen. Una verificacin de terreno confirmara o negara la existencia de dichas poblaciones.
Figura 10: Imagenes resultantes de la desmezcla espectral. a) Distribucin de la abundancia de mirabilita, b) yeso-1, c) yeso-2 y d) error del mtodo.
9.- Conclusiones
Como dijimos en la introduccin de este trabajo y lo volvemos a recalcar en este acpite, el motivo del presente trabajo es mostrar a la comunidad cientfica boliviana la existencia de otros mtodos capaces de diferenciar minerales a un bajo costo y con una precisin aceptable.El captor ASTER con la informacin espectral que mensura es capaz de diferenciar los minerales secundarios presentes en Coipasa. El advenimiento de la tecnologa hyperespectral embarcada sobre los satlites, como ya es el caso de Hyperion, facilitara este trabajo y en un futuro muy cercano las cartas geoqumicas obtenidas a partir de los satlites sern muy comunes y de nuestra adaptacin a este nuevo tipo de tcnicas depender en cierta medida el futuro de la comunidad cientfica nacional.
10.- Agradecimientos
Los autores agradecen a las personas que participaron en el proyecto Coipasa: Sophie Fis, Andr Lannoye, CharlesEtienne 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. Ballivin, O. & Risacher, F. 1981. Los salares del Altiplano Boliviano. Mtodos de estudio y estimacin econmica. UMSA ORSTOM, 241p. Boardman, J. W. Kruse, F.A. & Green, R. O. 1995. Mapping 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. Bryant. R. G. 1996. Validated linear mixture modelling of Landsat TM data for mapping evaporite minerals on a playa surface: methods and applications. International Journal of Remote Sensing, 17 (2): 315-330. Bryant, R. G. 2003. Monitoring playa systems using remote sensing: methods and applications. AridLands. 51. (accessible via: http://ag.arizona.edu/OALS/ALN/aln51). Camargo, R. A. Ponzoni, F. J. Zullo, Jr. J. Queiroz, G. & Arnaud, Y. 2003. Characterization of the Salar de Uyuni for in-orbit satellite calibration. IEEE Transactions on Geoscience and Remote Sensing, 41 (6): 1461-1468. Chapman, J. E. Rothery, D. A. Francis, P. W. & Pontual, A. 1989. Remote Sensing of evaporite mineral zonation in salt flats (salars). International Journal of Remote Sensing, 10 (1): 245-255.
Figura 9: Espectros de los endmembers resultantes del anlisis, comparados con espectros de laboratorio (USGS) remuestreados en funcin de los canales de ASTER.
La zona de la miribilita se encuentra en la parte Norte de la imagen y en una poblacin menos intensa en la parte Este de la imagen. Este mineral en su forma cristalina fue encontrado en profundidad por el proyecto Coipasa y la ocurrecia de dicha poblacion necesita ser verificada sobre terreno. Cabe recalcar que la mayor parte de las descripciones superficiales fueron realizadas con una lupa de gelogo y diferenciar un agregado amorfo de mirabilita de un agregado amorfo de halita es un trabajo muy difcil de efectuar. La tercera zona a la que denominamos yeso-2 se encuentra entre las dos zonas anteriormente descritas. El error del mtodo se encuentra al Sud Oeste de la imagen, esta zona no tiene ninguna relacin con las dos primeras zonas.
Crowley, J. K. 1991. Visible and near-infrared (0.4 2.5 um) reflectance spectra of playa evaporite minerals. Journal of Geophysical Research, 96 (B10): 16231-16240. Crowley, J. K. 1993. Mapping playa evaporite minerals with AVIRIS data: a first report from Death Valley, California. Remote Sensing of Environment, 44: 337-356. Crowley, J.K. 1996. Mg- and K-Bearing borates and associated evaporites at Eagle Borax Spring, Death Valley, California: a spectroscopic exploration. Economic Geology, 91: 622-635. Crowley, J. K. & Clark, R. N. 1992. AVIRIS study of Death Valley evaporite deposits using least-squares band-fitting methods. In: Summaries of the 3th Annual JPL Airborne Geoscience Worshop, Jet Propulsion Laboraty, Pasadena, California. JPL Publication 92-14, 1: 29-31. Crowley, J. K. & Hook, S. J. 1996. Mapping playa evaporite minerals and associated sediments in Death Valley, California, with multispectral thermal infrared images. Journal of Geophysical Research, 101 (B1): 643-660. Crowley, J. K. & Mars, J. C. 2000. Mapping evaporite minerals in the Death Valley salt pan using MODIS/ASTER airborne simulator (MASTER) data. Fourteenth International Conference on Applied Geologic Remote Sensing, Las Vegas, Nevada,USA: 344-345. Drake, N. A. 1995. Reflectance spectra of evaporite minerals (400 2500 nm): applications for remote sensing. International Journal of Remote Sensing, 16 (14): 2555-2571. Ericksen, G. E. Vine, J. D. & Ballon, R. 1978. Chemical composition and distribution of lithium-rich brines in salar de Uyuni and nearby salars in Southwester Bolivia. Energy, 3: 355-363. Eugster, H. & Hardie, L. 1978. Saline Lakes. In Lakes, Chemistry, Geology, Physics. A. Lehnman ed. SpringerVerlag, 237-293. Fis, S. 1996. Contribution de la prospection gophysique a ltude du Salar de Coipasa (Oruro, Bolivie). Travail de fin dtudes dIngnieur Civil Gologue, Universit de Lige, Belgique, 151p. Fornari, M. Risacher, F. & Fraud. G. 2001. Dating of paleolakes in the central Altiplano of Bolivia. Palaeogeography, Palaeoclimatology, Palaeoecology, 172: 269-282. Fujisada, H. 1996. ASTER Level-1 data processing algorithm. IEEE Transactions on Geoscience and Remote Sensing, 36 (4): 1101-1112. Green, A. A. Berman, M. Switzer, M. & Craig, M. D. 1988. A transformation for ordering multispectral data in terms of image quality with implications for noise removal. IEEE Transactions on Geoscience and Remote Sensing, 26 (31): 65-74. Gutierrez, J. 1997. Aplicacin del G.P.S. y la gravimetra para un estudio preliminar del paleo-relieve interno del Salar de Coipasa. Tesis de Grado Carrera de Ingeniera Geolgica FNI UTO.
Hardie, L. & Eugster, H. 1970. The evolutions of closedbasin brines. Mineral Society American Special Paper, 3: 273-290. Houssa, C.-E. 1998. Analyse mineralogique des ressources du Salar de Coipasa, Bolivie. Travail de fin dtudes dIngnieur Civil Gologue, Universit de Lige, Belgique, 130p Hunt, G.R. Salisbury, J.W. & Lenhoff, C.J. 1971. Visible and near-infrared spectra of minerals and rocks, IV, Sulfides and Sulfates. Modern Geology, 3: 1-14. Iwasaki, A. Fujisada, H. Asao, H. Shindou, O. & Akagi, S. 2002. Enhancement of spectral separation Performance for ASTER/SWIR, Infrared Spaceborne Remote Sensing. IX, Proceedings SPIE Vol: 4486, Paper-No, 4486-06: 42-50. Keshave, N. & Mustard. J. F. 2002. Spectral unmixing. IEEE Signal Procesing Magazine,19 : 44-57. Lebrun, V. 1997. Etude gochimique du Salar de Coipasa, Bolivie. Travail de fin dtudes dIngnieur Civil Gologue, Universit de Lige, Belgique, 140p. Lebrun, V. Pacosillo, P. Gutierrez, J. Cceres, F. Pool, E. & Pirard, E. 2002. Geochemistry of bitter brines in the Salar de Coipasa (Bolivia). Aardkundige Mededelingen, 12: 21-24. Rettig, S. L. Davis, J.R. Smith, R. L. Smith, K. A. & Ericksen, G. E. 1979. Bolivia. U.S. Geological Survey Professional Paper, 321 p. Risacher, F. 1992. Gochimie des lacs sales et crotes de sel de lAltiplano bolivien. Sciences Gologiques Bulletin, 45 (3-4): 219p. Risacher, F. & Fritz, B. 1991. Quaternary geochemical evolution of the salars of Uyuni and Coipasa, Central Altiplano, Bolivia. Chemical Geology, 90: 211-231. Risacher, F. & Fritz, B. 1992. Mise en evidence dune phase climatique holocne extrmement aride dans lAltiplano Central par la prsence de la polyhalite dans le Salar de Uyuni (Bolivie). Comptes Rendus de lAcadmie des Sciences, Paris, 314 (t. II) : 1371-1377. Pacosillo, P. 1998 Estudio preliminar de la cuenca evaportica del Salar de Coipasa, en base al muestreo de ros, fuentes y salmueras. Tesis de Grado Carrera de Ingeniera Geolgica FNI UTO. Plaza, A. Martinez, P. Perez, R. M. & Plaza, J. 2002. Spatial/ spectral endmember extraction by multidimensional morphological operations. IEEE Transactions on Geoscience and Remote Sensing, 40 (9): 2025-2041. Pool, E. 1995. Etude prliminaire du Salar de Coipasa (OruroBolivia) sur base des donnes gologiques, gochimiques et minralogiques. Mmoire de Matrise en Gologie des Terrains Superficiels, Universit de Lige, Belgique, 136p. Sabins, F. F. 1997. Remote sensing, principles and interpretation. W. H. Freeman & Company, New York, 3rd. Ed. Sabins, F. F. 1999. Remote sensing for mineral exploration. Ore Geology Reviews, 14: 157-183.
USGS Geobol. 1992. Geology and Mineral Resources of the Altiplano and Cordillera Occidental, Bolivia. U.S. 195 209 Geological Survey Bulletin, 1975: Sabins F. F. & Miller , R. M. 1994. Resource assessment Salar de Uyuni and vinicity. Tenth Thematic Conference on Geologic Remote Sensing, Vol. I: I92-I103. Servant, M. & Fontes, J. C. 1978. Les lacs quaternaires des plateaux des Andes boliviennes premires interprtations paloclimatiques. Cahiers ORSTOM, Srie Gologie, 10 (1): 9-23. Sylvestre, F. Servant-Vildary, S. Servant. M. 1998. Le dernier maximum glaciaire (21 000-17000 14C ans B. P.) dans les Andes tropicales de Bolivie daprs ltude des diatomes. Comptes Rendus de lAcadmie des Sciences, Paris, 327 : 611-618. Sylvestre, F. Servant, M. Servant-Vildary, S. Causse, C. Fournier, M. & Ybert, J.-P. 1999. Lake-level chronology on the southern Bolivian Altiplano (18-23S) during lateglacial time and the Early Holocene. Quaternary Research, 51 (1): 54-66.
Van der Meer. F. 2001. Basic physics of spectrometry. In Imaging spectrometry: basic principles and prospective applications. F. Van der Meer & S. de Jong, ed. Kluwer Academic Publishers, 3-16. Van der Meer, F. de Jong, S. & Bakker, W. 2001. Imaging spectrometry: basic analytic techniques. In Imaging spectrometry: basic principles and prospective applications. F. Van der Meer & S. de Jong, ed. Kluwer Academic Publishers, 16-61. White, K. & Drake, N. 1993. Mapping the distribution and abundance of gypsum in South-Central Tunisia from Landsat Thematic Mapper data. Zeitschrift fur Geomorphologie, 37 (3): 309-325. Yamaguchi, Y. Kawkami, T. Kahle, A. B. Pniel, M. & Tsu, H. 1998. ASTER mission planning and operations concept. SpaceOps, International Symposium on Technical Interchange for Space Mission Operations and Ground Data Systems: 8p.