Boletn geogrfico. Ao XXX N31 - 2008, pp 13 - 25. Departamento Geografa. Edicin especial: VII Jornadas Patagnicas de Geografa.

Universidad Nacional del Comahue. Neuqun. ISSN 0326-1735

VARIABILIDAD DE LA TEMPERATURA DE SUPERFICIE EN EL NW DE LA PROVINCIA DEL NEUQUEN A PARTIR DE INFORMACIN SATELITAL Marisa Gloria Cogliati1 Mariana Ins Surez2

Resumen La complejidad del terreno y la falta de datos meteorolgicos siguiendo el gradiente altitudinal obliga a la exploracin de nuevos mtodos o esquemas entre los que se encuentra el anlisis por medio de modelos o imgenes satelitales. El presente trabajo pretende describir el marco de informacin con el que se cuenta en el NW neuquino y la introduccin de anlisis espacial de distribucin de la temperatura de superficie a partir de imgenes NOAA- AVHRR con miras a la inclusin futura de nuevos esquemas de estudio. La temperatura del aire presenta gran variabilidad espacial, principalmente en los valores mnimos, una idea de esta distribucin puede asociarse a la temperatura de superficie que es obtenida a partir del anlisis de imgenes NOAA-AVHRR. Este trabajo presenta el anlisis de situaciones particulares de la temperatura de superficie de la zona NW de la Provincia de Neuqun utilizando imgenes de satlites NOAA en alta resolucin (AVHRR) e imgenes LANDSAT de la regin. La temperatura de superficie obtenida utilizando la informacin de los canales 4 y 5 del sensor AVHRR presenta la distribucin espacial de la temperatura del suelo que pone en evidencia la complejidad del comportamiento, asociado al relieve. Se observa la dependencia de las isotermas con la altura, alcanzando los menores valores de temperatura de superficie en los picos. Esta metodologa permitira incorporar informacin cuya escala considerara con mayor detalle que en escala sinptica de la temperatura de superficie de pxeles de aproximadamente 1 km2. El anlisis posterior de una mayor cantidad de situaciones, permitira realizar mapas estacionales de temperatura de superficie en una escala que no es posible obtener con el instrumental que se encuentra instalado en este momento, pudiendo asociarse posteriormente al comportamiento de la temperatura del aire.
1

Departamento de Geografa, Facultad de Humanidades, Universidad Nacional del Comahue. E-mail: [email protected] 2 Departamento de Geografa, Facultad de Humanidades, Universidad Nacional del Comahue. E-mail: [email protected]

13

Boletn Geogrfico

La mediciones de superficie de la temperatura media del aire en distintas estaciones de la regin en el perodo 1993 a 1997 presenta una gran amplitud anual (aprox 32C) en todas las estaciones. La temperatura media mxima es de 29 C en Puesto Vallejos y 18.9 C en Pampa de Chacaico mientras que la temperatura media mnima es -12.5 C en la misma estacin. Las estaciones ubicadas en el oeste de la regin (Cajn de los Chenques, Las Lagunas, Buta Malln, Estacin Chacaico) presentan menores valores de temperatura mnima media y mxima media, mientras que las estaciones ubicadas hacia el este y al sur de la regin presentan mayores valores. Palabras clave: topoclimatologa, temperatura de superficie, satlites

Introduccin El estudio de la temperatura y otras variables meteorolgicas en zonas de terreno complejo, tales como el noroeste de la Provincia del Neuqun, son comparativamente escasos en relacin a estudios realizados para zonas de valle o relieves llanos. Los principales factores se deben a la falta de datos meteorolgicos siguiendo el gradiente altitudinal y la complejidad del terreno. El rea de estudio seleccionada adems, por su gran extensin y problemas asociados a comunicaciones, presenta inaccesibilidad en pocas invernales. Por ello se hace necesario incorporar mtodos alternativos para el estudio de variables meteorolgicas. La bsqueda de nuevos mtodos se encuentra en un proceso de cambio permanente, adquiriendo gran importancia el anlisis por medio de modelos fsicos o de sensores remotos. Este trabajo presenta el anlisis de la temperatura de brillo y de superficie del satlite LANDSAT 7 ETM+ para una escena que comprende el NW de la Provincia de Neuqun, a partir de informacin del canal trmico, para describir la distribucin espacial de la temperatura de superficie en escala local o regional y en zonas con muy baja densidad de datos en superficie. El anlisis posterior comprendiendo una mayor cantidad de imgenes, permitir realizar mapas temticos de distinta escala temporal y su posterior correlacin con observaciones en superficie. En satlites como los LANDSAT TM y ETM+ es ms complejo obtener valores de temperatura de superficie que en el caso de los satlites NOAA- AVHRR dado que poseen una sola banda trmica y requieren para el clculo modelos de transferencia de radiacin y emisividad de superficie. Por esta razn la informacin de la banda trmica de los LANDSAT TM y ETM no ha sido extensivamente usada para mapeos de temperatura de superficie, aunque posea una resolucin espacial alta (Li y otros, 2000). Las imgenes satelitales LANDSAT TM han sido ampliamente estudiadas para diferentes aplicaciones. Las principales dificultades para el clculo de la temperatura de superficie residen en problemas con la calibracin, baja frecuencia temporal y dificultades en la seleccin de un algoritmo apropiado (Qin y otros, 2001). En 14

Variabilidad de la temperatura de superficie en el NW de la provincia del Neuqun ...

general la banda trmica del LANDSAT TM se usa para determinar la temperatura de brillo al nivel de satlite (Mansor, 1994; Saraf, 1995; Zhang, 1997) o simplemente se usa el valor del numero digital (DN) para las mismas aplicaciones (Ritchie, 1990; Oppenheimer, 1997). La temperatura del aire presenta una importante variabilidad espacial y temporal, asociada a variaciones de altura del terreno, y variacin diurna y estacional. Alturas mayores a 600 m son suficientes para causar diferencias verticales de elementos climticos y cobertura vegetal (Barry, 1981). Los factores geogrficos que influencian los elementos meteorolgicos en zonas de montaa son: la latitud, la altitud y la topografa. La radiacin solar, radiacin neta y la temperatura decrecen en mayores latitudes produciendo un descenso de las lneas de vegetacin y nieve en altura y, adems, influencian el ritmo de la variacin estacional y diurna. El efecto de la altitud sobre los elementos climticos se observa sobre la presin, la densidad del aire, la presin de vapor y la insolacin. La presencia de altos sistemas montaosos provoca diferencias notorias de temperatura entre las cumbres y el piso de los valles. Las oscilaciones diarias son tambin marcadas y se incrementan con la altura. El comportamiento anual de la temperatura es estacional, presentando los mximos valores en los meses de verano y los mnimos en los meses de invierno. La ubicacin en zona de latitudes medias y la influencia de las masas de aire determinan tal variabilidad. La pendiente y la orientacin del terreno marcan efectos fundamentales en el ingreso de radiacin y en las condiciones de temperatura (Barry, 1981) En zonas de altas montaas, especialmente sobre la lnea de vegetacin, la vida de plantas y animales es controlada por el clima cercano a la superficie (microclima). Adems, la distribucin espacial del microclima, forma un mosaico debido a los efectos de la topografa que causa topoclimas distintivos. La combinacin de efectos topo y microclimticos determina la variabilidad en la cobertura vegetal (Barry, 1981). La gran intensidad de radiacin solar en altas altitudes resulta en alta temperatura de superficie. Turner (1958) midi temperaturas de 80C en julio de 1957 en humus oscuro a 2.070 m de altura en tztal, Austria. El sitio tena una exposicin SW, con una pendiente de 35 y con una temperatura del aire a 2 m de 30 C. La diferencia entre la temperatura del suelo y del aire depende de las condiciones de tiempo meteorolgico. Mediciones detalladas a 2.072 m cerca de Obergrgl en Austria (Aulizky, 1962) mostr que existe una relacin lineal entre los valores medio y extremos de la temperatura del suelo y del aire cuando el suelo no est congelado. En este trabajo, se selecciona una metodologa para la obtencin de la temperatura de superficie, a partir de la banda 6 de los satlites LANDSAT 7 ETM+ en la zona del NW de la provincia del Neuqun y de la informacin disponible en superficie.

15

Boletn Geogrfico

rea de Estudio El rea de estudio comprende el noroeste de la Provincia del Neuqun, comprendiendo los Departamentos Minas y Chos Malal desde los 36 10 S a 37 30 S y los 69 30 W a 71 15 W. La conformacin fsica del rea presenta rasgos particulares. La topografa est dominada por tres cordones orogrficos paralelos alineados longitudinalmente y separados entre s por depresiones que conforman valles fluviales. Estos cordones son la Cordillera de los Andes al oeste, la Cordillera del Viento en el centro y el Macizo Tromen en el este. Las alturas mximas se ubican en los macizos del centro y el este de la regin y corresponden al volcn Domuyo (4.709 m s.n.m.) y al volcn Tromen (3978 m s.n.m.). Las lneas de nivel del terreno y la zona de estudio se detallan en la Figura 1(U.S. Geological Survey, Earth Resources Observation and Science (EROS)).. Datos utilizados En este estudio se analiz una imagen LANDSAT 7 con el sensor ETM+ (Enhanced Thematic Mapper Plus) adquirida el 7 de febrero de 2001 (path 232/row 086). Comprendiendo el sector del NW de la provincia de Neuqun representado en la Figura 1. La imagen fue creada por U.S. Geological Survey (USGS) y contiene los datos en formato GeoTIFF (Geographic Tagged Image-File Format ), producto L1G.

Figura 1: Curvas de nivel del sector NW de la Provincia de Neuqun (las lneas negras representan el lmite provincial e internacional de la provincia). Datos: USGS, Shuttle Radar Misin.

16

Variabilidad de la temperatura de superficie en el NW de la provincia del Neuqun ...

Metodologa Un xito limitado se ha alcanzado para determinar la temperatura de superficie utilizando datos trmicos de los satlites LANDSAT Thematic Mapper (TM) y ETM+. Los satlites TM y ETM+ tienen una banda infrarroja (10.44 - 12.42 m), por lo que el uso de un algoritmo para correccin del tipo split-window es imposible, pero la alta resolucin espacial (120 m para TM y 60 m para ETM+) hace el estudio trmico infrarrojo muy atractivo desde el punto de vista local y regional. En este trabajo se analiza la temperatura de brillo del canal trmico del satlite LANDSAT 7 con el sensor ETM + y se estima la temperatura de superficie a partir de datos de emisividad calculados por medio del ndice de vegetacin (NDVI) y la estimacin de la trasmitancia atmosfrica (Zhang y otros, 2006). La escena de satlite se proces utilizando el software para procesamiento de imgenes raster, Idrisi versin Kilinmanjaro. Tratamiento de Imgenes Para el clculo de la temperatura de brillo del canal 6 se efectu el clculo de la radianza espectral. A partir de los niveles digitales de la banda trmica, se calcul la radianza espectral (Lsensor) utilizando la Ec. 1: Lsensor = gain * DN + bias (1)

Donde Lsensor es la radianza espectral en la banda trmica en W/(m2 ster mm); gain es la pendiente de la funcin de conversin radianza/DN; DN es el nmero digital en un determinado pxel a partir de un producto L1G; y bias es la ordenada al origen de la ecuacin e conversin radianza/DN (Landsat Project Science Office, 2002). Los valores de gain y bias pueden encontrarse en el archivo encabezado de las imgenes ETM+. Los valores de radianza de la banda trmica de los satlites ETM+ fueron transformados a temperatura de brillo en el tope de la atmsfera, usando constantes de calibracin proporcionadas por la oficina del Proyecto Landsat (2002) de acuerdo a la Ec (2):

Tsensor

K2 ln( K 1 / Lsensor 1)

(2)

Donde Tsensor es la temperatura efectiva del satlite (temperatura de brillo) en K y K1 y K2 son constantes de calibracin de pre lanzamiento, para Landsat 7 ETM+, K1 = 666.09 W/(m2 sr mm) y K2 =1282.71 K (Landsat Project Science Office, 2002). Los valores de sesgo (bias) y ganancia (gain) se obtienen del archivo de encabezado (header).

17

Boletn Geogrfico

La temperatura calculada en la Ec. (2) no es la temperatura de superficie, si no que es la temperatura de brillo al tope de la atmsfera (Zhang y otros, 2006). Para obtener la temperatura de superficie a partir de la temperatura de brillo, podemos utilizar varios modelos distintos, sin embargo, dichos modelos requieren de extensa informacin de superficie que no siempre est disponible en nuestra zona. Qin y otros (2001) desarrollaron el algoritmo de mono ventana, que se expresa:

Tsup

a 6 (1 C 6 D6 ) [b6 (1 C 6 D6 ) C 6 D6 ]Tsensor D6Ta , (3) C6

Donde Tsup es la temperatura de superficie en K y Tsensor es la temperatura de brillo en K calculado a partir de la banda 6 del satlite Landsat ETM+, Ta es la temperatura efectiva media atmosfrica (K), a6 y b6 son constantes: a6 es -67.355351 y b6 es 0.458606 para la temperatura de superficie entre 273.5 y 343.5 K. C6 and D6 deben calcularse de acuerdo a Ec. (4) y Ec.(5):

C 6 6 D6 (1 6 )[1 (1 ) 6 ]

(4) (5)

Donde es la emisividad de la superficie y 6 es la transmitancia atmosfrica. 6 puede ser estimado a partir del contenido de vapor de agua y Ta puede calcularse a partir de la ecuacin lineal para la atmsfera estndar en el verano de latitudes medias (6): (6) Donde Ta es la temperatura media efectiva de la atmsfera y To es la temperatura del aire cercana a la superficie. Clculo de la emisividad, a partir del ndice verde: De acuerdo con Van de Griend and Owes (1993) la relacin entre los valores de NDVI y los de emisividad pueden expresarse segn (7) para rangos de NDVI entre 0.157 y 0.727. 1.0094 0.047 ln( NDVI ) (7) donde es la emisividad y NDVI es el ndice de Vegetacin en Diferencias Normalizadas calculado a partir de (8).

Ta 16.0110 0.92621To

NDVI

(banda 4 banda3) (banda 4 banda3)

(8)

18

Variabilidad de la temperatura de superficie en el NW de la provincia del Neuqun ...

Si el valor de NDVI se encuentra fuera del rango (0.1570.727), se divide en 5 intervalos y se puede incluir los valores estimados de emisividad de acuerdo a Zhang y Wang (2001). Estos valores emisividad corresponden a 0.99 cuando el ndice de NDVI est entre 0.727 y 1.0 ; a 0.955 cuando el NDVI es entre 0.157 y -0.18; a 0.965 cuando el NDVI es de -0.18 a -1.0 y a 0.985 cuando el ndice de NDVI es de 1.0. Al efectuar el clculo del ndice de vegetacin es necesario adecuar la resolucin espacial de las bandas 3 y 4 a una resolucin similar a la de la banda 6, que se reduce a la mitad. Muchos factores afectan la recuperacin de la temperatura de superficie a partir de los datos del infrarrojo trmico del satlite, pero algunos de ellos como la transmitancia, la humedad del aire, la radianza entrante y saliente son difciles de obtener. Discusin y Resultados

Figura 2: Radianza espectral calculada mediante Ec. (1) en la imagen satelital Landsat ETM+, del 7 de febrero de 2001 (path 232/row 086).

19

Boletn Geogrfico

Figura 3: Temperatura de Brillo en el tope de la atmsfera calculada mediante Ec. (2) en la imagen satelital Landsat ETM+ del 7 Febrero de 2001 (path 232/row 086). La imagen de temperatura de brillo del canal trmico del LANDSAT 7 ETM+ (Figura 3) correspondiente al NW de Neuqun indica que la temperatura presenta una variabilidad espacial, que se correlacionara bien con las diferencias de altura de la zona que se presentan en la Figura 1. Las zonas ms bajas del terreno, que coincidiran con el piso de los valles fluviales, presentan valores altos de temperatura, de alrededor de los 30 C a 37 C (ver Figura 3). En las zonas de mayor altura, tales como las cumbres de los cordones montaosos de la Cordillera del Viento y el Macizo Tromen, se exhiben los valores ms bajos de temperatura de brillo, que en general estn en el orden de los 1 C a los 8 C. En estos casos deber analizarse desde el punto de vista de la cobertura de nieve, punto que ser analizado en futuros trabajos. Los mnimos valores considerados corresponderan a geoformas de altura ubicadas principalmente en el sector sudoeste de la imagen. Esto se asociara a los factores de altitud y latitud.

20

Variabilidad de la temperatura de superficie en el NW de la provincia del Neuqun ...

Figura 4: NDVI calculado mediante Ec. (8) para la imagen satelital LANDSAT ETM del 7 Febrero de 2001 (path 232/row 086) para el NW de la provincia de Neuqun.

Figura 5: Emisividad calculada mediante Ec. (7) para la imagen satelital LANDSAT7 ETM+ del 7 Febrero de 2001 (path 232/row 086) para el NW de la provincia de Neuqun. 21

Boletn Geogrfico

Figura 6: Temperatura de superficie (C) calculada mediante Ec. (3) para la imagen satelital LANDSAT7 ETM+ del 7 Febrero de 2001 (path 232/row 086) para el NW de la provincia de Neuqun. El clculo de NDVI en el NW de Neuqun expresado en la Figura 4 indica una variabilidad espacial notoria de la cobertura vegetal, que disminuye de oeste a este en forma progresiva. Los valores mximos del orden de los 0.4 a 0.5 se encuentran principalmente en el sector oeste de la escena y tambin en los valles fluviales que siguen los cursos de agua tributarios del ro Cur Leuv. Los valores mnimos entre 0.1 y 0.4- se observan en gran parte de la imagen: sector centro y este. La variabilidad espacial que presenta la cobertura vegetal se asociara en parte a las diferencias de precipitacin que se presenta a sotavento de la cadena montaosa. Los ndices de cobertura en el resto de la imagen presentan un valor ms homogneo, con valores entre 0.1 y 0.4. Esta distribucin presenta una variacin del ndice de cobertura vegetal en esa zona, que estara asociado a la influencia de los valores de precipitacin y humedad. La Figura 6, presenta la temperatura de superficie en C calculada aplicando la mono ventana de Qin (2001), el procedimiento requiere de una estimacin de la trasmitancia atmosfrica en funcin de la humedad del aire y la temperatura media del aire simultnea con el momento de la obtencin de la imagen. Los datos medios se estimaron a partir de datos de los reanlisis del NCEP (Nacional Centres of Atmospheric Research). Este procedimiento corrige los valores en el tope de la atmsfera a partir de la consideracin de las condiciones atmosfricas.

22

Variabilidad de la temperatura de superficie en el NW de la provincia del Neuqun ...

Conclusiones La determinacin de la temperatura de superficie utilizando datos trmicos de los satlites LANDSAT ETM+, reviste una mayor dificultad que en el caso de los satlites NOAA. Los satlites LANDSAT ETM+ tienen una banda infrarroja (10.44 - 12.42 m), por lo que el uso de un algoritmo para correccin del tipo split-window es imposible. Sin embargo, la alta resolucin espacial (60 m para el sensor ETM+) hace el estudio trmico infrarrojo de estos satlites, muy atractivo desde el punto de vista local y regional. En este trabajo se analiza la temperatura de brillo y temperatura de superficie del canal trmico del satlite LANDSAT 7 con sensor ETM + en una escena que comprende el NW de la provincia de Neuqun. La temperatura de superficie en C fue calculada aplicando la mono ventana de Qin (2001), el procedimiento requiere de una estimacin de la trasmitancia atmosfrica en funcin de la humedad del aire y la temperatura media del aire simultnea con la obtencin de la imagen. Este anlisis brinda informacin distribuida espacialmente de las variaciones de temperatura de superficie en pxeles con 60 m de resolucin, mostrando concordancia con variabilidad espacial asociada a la topografa. Este hecho permitira realizar anlisis de zonas con baja densidad de informacin meteorolgica o difcil acceso. Sin embargo, debe contarse para su adecuado tratamiento de datos meteorolgicos de superficie tales como contenido de humedad atmosfrica y temperatura media del aire cercano a la superficie, simultneos al pasaje del satlite.

23

Boletn Geogrfico

Referencias AULITZKY, H. (1962) Die Bodentemperaturverhltnisse einer zentralalpinen Hochgebirgs-Hangstation II. Archiv. Met. Geophys. Biokl., B, vol. 11, 301-62. BARRY, R.G. (1981) Mountain weather and climate, Methuen. London and New York. 309 pp. CRESSWELL, M. P.; MORSE, A. P.; THOMSON, M. C.AND CONNOR. S. J. (1999). Estimating surface air temperature, form Meteosat land surface temperatures, using an empirical solar zenith angle model. Int. J. Remote Sensing, Vol. 20, No. 6, 1125-1132. HUMES, K. S.; KUSTAS, W. P.; MORAN, M. S.; NICHOLS, W. D. AND WELTZ, M. A. (1994) Variability of emissivity and surface temperature over a sparsely vegetated surface. Water Resources Research, vol. 30, 12991310. LI, F.; JACKSON, T.; KUSTAS, W.; SCHMUGGE, T.; FRENCH, A.; COSH, M. and Bindlish, R. (2004) Deriving land surface temperature from Landsat 5 and 7 during SMEX02/SMACEX. Remote Sensing of Environment, vol. 92, issue 4, 30 September 2004, 521 534. MANSOR, S. B; CRACKNELL, A. P; SHILIN, B. V. and GORNYI, V. I. (1994) Monitoring of underground coal fires using thermal infrared data. International Journal of Remote Sensing, vol. 15, 16751685. OPPENHEIMER, C. (1997) Remote sensing of the color and temperature of volcanic lakes. International Journal of Remote Sensing, Vol. 18, 537. PRICE, J. C. (1984) Land surface temperature measurements from the split window channels of the NOAA 7 Advanced Very High Resolution Radiometer. Journal of Geophysical Research, vol. 89, 72317237. QIN, Z.; KARNIELI, A. and BERLINER, P. (2001) A mono-window algorithm for retrieving land surface temperature from Landsat TM data an its application to the Israel Egypt border region. International Journal of Remote Sensing, vol. 22, n 18, 3719 3746. RITCHIE, J. C.; COOPER, C. M. and SCHIEBE, F. R. (1990) The relationship of MSS and TM digital data with suspended sediments, chlorophyll, and temperature in Moon Lake,Mississippi. Remote Sensing of Environment, vol. 33, 137148. SARAF, A. K.; PRAKASH, A.; SENGUPTA, S. and GUPTA, R. P. (1995) Landsat TM data forestimating ground temperature and depth of subsurface coal fire in the Jharia coafield, India. International Journal of Remote Sensing, vol.16, 21112124. TURNER, H. (1958) Maximal tempemperaturen oberflchennaher Bodenschichten au der alpinen Waldgrenze. Wetter u. Leben, vol. 10, 1-12. VAN DE GRIEND, A. A. and OWE, M. (1993) On the relationship between thermal emissivity and the normalized different vegetation index for natural surfaces. International Journal of Remote Sensing, vol. 14, No 6, 11191131. YANG, J. and WANG, Y.Q. (2001) Estimation of land surface temperature using landsat-7 etm+ thermal infrared and weather station data . Proceedings of Huangshan International Thermal Infrared Remote Sensing Workshop, July 14-17, Huangshan, Anhui, P.R.China. 24

Variabilidad de la temperatura de superficie en el NW de la provincia del Neuqun ...

ZHANG, J, WANG, Y. and LI, YAN (2006) C++ program for retrieving land surface temperature from the data of Landsat TM/ETM+ band6. Computers & Geosciences, vol. 32 17961805. ZHANG, X.; VAN GENDEREN, J. L. and KROONENBERG, S. B. (1997) A method to evaluate the capability of Landsat 5 TM band 6 data for sub pixel coal fire detection. International Journal of Remote Sensing, vol. 18, 32793288.

25