Source: http://www.slideshare.net/senarap/vol28-junio09
Timestamp: 2015-03-06 09:33:30+00:00

Document:
Evaluacion del impacto ambiental de...
by Sonny Swokowsky B...
Manual EIA Evaluación impactos ambi...
by smpacheco1
by fhbernal
128221 views
Impactos de cambios de uso de suelo...
Soza diego.2012.Metodologia de mode...
by Diego Soza
by Julio Cortes Leyama
Inventario hídrico Río san Francisc...
Articulocientificotesishidalgo 2012
by Victor Hidalgo Co...
VOL. 28, Nº 1 JUNIO 2009 ISSN 0717-2915SOCIEDAD LATINOAMERICANA DE PERCEPCIÓN REMOTA Y SISTEMAS DE INFORMACIÓN ESPACIALSOCIEDADE LATINO-AMERICANA EM SENSORIAMENTO REMOTO E SISTEMAS DE INFORMACAO ESPACIAILATINAMERICAN SOCIETY FOR REMOTE SENSING AND SPACE INFORMATION SYSTEMS Ciudad de Puerto Vallarta - México Imágen de satélite QuickBird de DigitalGlobe 2.
SOCIEDAD LATINOAMERICANA DE PERCEPCION REMOTA Y SISTEMAS DEINFORMACION ESPACIAL SOCIEDADE LATINO-AMERICANA EM SENSORIAMENTO REMOTO E SISTEMAS DE INFORMACAO ESPACIAI LATINAMERICAN SOCIETY ON REMOTE SENSING AND SPACE INFORMATION SYSTEM SELPER VOL. 28 Nº 1 JUNIO 2009 ISSN 0717-2915 FOTO PORTADA Puerto Vallarta - México Imagen de satélite QuickBird de DigitalGlobe, producto Pansharpened resolución espacial 0.60 metros, en composición Color Natural (3,2,1) de la Ciudad de Puerto Vallarta tomada en el mes de Octubre de 2008. Los sectores que se destacan son el Puerto y el Aeropuerto. Esta ciudad fue sede del VIII Simposio Selper.EJEMPLAR DE DISTRIBUCION GRATUITA. PROHIBIDA SU VENTA 3.
REVISTA02 VOL. 28 Nº 1 JUNIO 2009 DIRECCIONES DE CONTACTO SELPER DIRECTORIO SELPER, SEDE CUBA 2008 - 2010 UCT GEOCUBA IC Calle 4 No. 304, e/ 3ra y 3ra A Miramar, Playa, La Habana, Cuba CP 11 300 E-mail investigaciones@geocuba.co.cu PRESIDENTE Eloy Pérez Cuba VICE-PRESIDENTE VICE-PRESIDENTE VICE-PRESIDENTE Silvia Lucero Casas Myriam Ardila Torres Pedro Luis García Mexico Colombia Cuba SECRETARIO GENERAL RESPONSABLE DE DIVULGACIÓN ELECTRÓNICA TESOREROMaría Antonia García Cisneros Walter F. Sione Gastón Hernández Cuba Argentina Cuba COMITÉ DE RELACIONES INTERNACIONALES COMITÉ EDITORIAL COMITÉ DE EDUCACIÓN Paulo Roberto Martini (Brasil) María Cristina Serafini (Argentina) Tania María Sausen (Brasil) Presidente Presidente Presidente Aurelie Sand – Francia Miriam Esther Antes – Argentina Joao Avila - Brasil Maurizio Fea – Italia Leila María Fonseca – Brasil Maria Antonia García Cisnero - Cuba Alfredo Cuello – Argentina Francisca Celia González – Argentina Andres Cárdenas - Colombia Victor Barrena - Perú Fabián Lozano – México Ethel Rubín de Celís Llanos - Perú Ramiro Salcedo – Venezuela Jorge Luis Martin – Cuba Myriam Ardila Torres / Colombia Freddy Flores. – Venezuela María Elena Pezoa - Chile Olga Piedad Rudas - Colombia Anyul del Pilar Mora - Colombia Luis Geraldo Ferreira - Brasil Pedro Martínez Fernández - Cuba CAPÍTULOS CONSTITUIDOS GUATEMALA CAPÍTULOS ESPECIALES Carlos Alberto Duarte ARGENTINA Ingeniería Virtual ALEMANIA Francisca Celia González Ruta 4, 6-49 Zona 4, Oficina 14 Klaus Reiniger Universidad Nacional del Sur Ciudad de Guatemala (01004), Guatemala DLR Departamento de Geología Tel: 502 - 334-1039/4038 D-8031 Oberpfaffenohfen San Juan 670 (8000) Fax: 502 - 331-9390 Alemania Bahía Blanca, Argentina E-mail: chduarte@hotmail.com Tel: 49- 8153- 281.189 Tel: 54 - 291 - 459 5102 - int. 4360 Fax: 49- 8153- 281.443 Fax: 54 - 291 - 459 5127 GUYANA FRANCESA E-mail: ghgonzal@criba.edu.ar Laurent PolidoriI CANADÁ Directeur de Recherche IRD / US ESPACE 140 Fritz P. Dubois BOLIVIA Institut de Recherche pour le Développement 25 Nidland Rodolfo Ayala (ex-ORSTOM) Crs Nepean Ontario Kh2-8n2 Carrera de Topografía y Geodesia Route de Montabo - BP 165 - 97323 Cayenne Ontario, Canadá Facultad Técnica cedex Tel: 613- 596-4164 Universidad Mayor de San Andrés Tel. (+594) 594 29 92 81 Fax: 613- 723-9626 Av. Arce 2299 1º Piso Fax (+594) 594 31 98 55 La Paz, Bolivia E-mail:polidori@cayenne.ird.fr ESPAÑA Tel.: 591-2-2441401 José L. Labrandero MÉXICO Consejo Superior de Investigaciones Científicas (CSIC) BRASIL Silvia Lucero Casas Pinar 25- Madrid 28006, España Leila María García Fonseca UAM Agronomía y Ciencias Tel: 34- 411.10.98 INPE Universidad Autónoma de Tamaulipas Fax: 34- 562.55.67 Av. Dos Astronautas 1758, Sao José dos Campos Centro Universitario Victoria San Pablo, Brasil Cd. Victoria Tamaulipas, México FRANCIA Tel: 55 - 12-39456000 E-mail: scasas@uat.edu.mx Aurelie Sand E-mail: leila@dpi.inpe.br CNES URUGUAY 31401- Toulouse Cedex 9 CHILE Antonio Alarcón France Gina Ghio Servicio de Sensores Remotos Aeroespaciales Tel: 33-1- 447.67500 AYS Ltda Fuerza Aérea Uruguaya Fax: 33-1-447.67849 Calle Teatinos 258 of 811 Ruta 101 s/n Km. 19500 E-mail: aurelie.sand@cnes.fr Santiago, Chile. Carrasco, Canelones, Uruguay Tel: 562 - 6903768 Tel.: 598 -2 601 4083 HOLANDA Fax: 562 - 6903790 Fax: 598 -2 601 4090 Carlos Valenzuela E-mail: jnovikov@codelco.cl E-mail: director@ssra.gub.uy ITC 350 Boulevard 1945, P.O.X. 6. 7500 AA COLOMBIA VENEZUELA Enschede, Holanda Myriam Ardila Ramiro Salcedo Tel.: 31 53 874-444 PROSIS S.A. Centro de Procesamiento Digital del Instituto de Fax: 31 53 874-400 Calle 94 A n° 11 A-70 Ingenieria en Caracas Santa Fe de Bogota, Colombia Apdo. Postal 40200 / Caracas, Venezuela ITALIA Tel: 57-1-6357272 Tel/fax: 58 - 212 - 903 -4682 Maurizio Fea Fax: 57-1-6357190 E-mail: ramiros@fii.org ESA/ESRIN E-mail: mardila@prosis.com Via Galileo Galilei, 00044 PERÚ Frascati, Italia CUBA Victor Barrena Arroyo Casilla Postal 64 Eloy Pérez García Universidad Nacional Agraria La Molina Tel: 39 - 694180940 UCT GEOCUBA IC Av. La Universidad s/n Fax: 39 - 694180942 Calle 4 No. 304, e/ 3ra y 3ra A La Molina, Lima, Perú E-mail: Maurizio.Fea@esa.int Miramar, Playa, La Habana, Cuba CP 11 300 Tel / Fax: 51-1-349-5647 anexo 232/349-2041 investigaciones@geocuba.co.cu E-mail: vbarrena@lamolina.edu.pe USA Patricia M. Ravelo ECUADOR CAPÍTULOS EN FORMACIÓN SPOT José R. Aguiar V. Estados Unidos CLIRSEN Tel: 1-800-ask-spot ext. 137 PARAGUAY Edif. Instituto Geográfico Militar, Piso 4 Fax: 703-648.1813 Sergio M. Burgos Sosa Seniergues s/n y Paz y Miño E-mail: ravelo@spot.com IPPA Apartado Postal 17-08-8216 Dr. César Sánchez 431 Quito, Ecuador San Lorenzo, Paraguay Tel.: 593- 2- 254-3193 Tel/Fax: 595- 21-574909 Fax: 593- 2- 255-5454 Email: smburgos@highway.com.py E-mail: clirsen@clirsen.com 4.
REVISTA VOL. 28 Nº 1SELPER INDICE TEMÁTICO Y COMITÉ EDITORIAL JUNIO 2009 03 COMITÉ EDITORIAL ÍNDICE TEMÁTICO María Cristina Serafini (Argentina) PRODITEL Universidad Nacional de Luján Cruce rutas 5 y ex 7 (6700) Luján, Buenos Aires, Argentina Tel: 54-2323-423171 int 248 SEGUIMIENTO DE LA REGENERACIÓN DE ÁREAS Fax: 54-2323-425795 QUEMADAS A PARTIR DE IMÁGENES E-mail: proditel@mail.unlu.edu.ar HIPERESPECTRALES AHS Miriam Esther Antes (Argentina) PRODITEL Federico González-Alonso Universidad Nacional de Luján Cruce rutas 5 y ex 7 Margarita Huesca Luján, Buenos Aires, Argentina Tel: 54-2323-423171 int 248 Sergio Martinez Sanchez-Palencia .........................................................5 Fax: 54-2323-425795 E-mail: proditel@mail.unlu.edu.arx Leila María Fonseca (Brasil) EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO INPE DE SUELO EN LA HIDROLOGÍA DE LA CUENCA DEL Av. Dos Astronautas 1758, Sao José dos Campos, Sao Paulo, Brasil RÍO SAN JUAN EN NUEVO LEÓN, MÉXICO Tel: 55 - 12-39456000 E-mail: leila@dpi.inpe.br Axayácatl Maqueda Estrada Fabián Lozano (Mexico) Jianhong Ren Instituto Tecnológico y de Estudios Superiores de Monterrey , Diego Fabián Lozano García ............................................................. 13 Av. Euganio Garza Sada # 2501 sur, Col. Tecnológico, Monterrey, Nuevo León, México ANÁLISE DOS FATORES QUE CONTROLAM A Tel: 52 - 81 - 8358 - 1400 ext 5275 Fax: 52 - 81 - 8358 - 6280 DISTRIBUIÇAO ESPACIAL DA TURBIDEZ NO E-mail: dflozano@itesm.mx RESERVATÓRIO DE ITUPARARANGA- SP Jorge Martín Chiroles (Cuba) E-mail: ciencia@geocuba.com.cu Leticia Sabo Boschi Maria de Lourdes Bueno Trindade Galo Jorge Luis Martín (Cuba) E-mail: ciencia@geocuba.com.cu Evlyn Márcia Leão de Moraes Novo..................................................... 24 Francisca Celia González (Argentina) Universidad Nacional del Sur Departamento de Geología ADECUACIÓN DEL DEM DEL SRTM DE 90 METROS San Juan 670 (8000) Bahía Blanca, Argentina PARA USOS HIDROLÓGICOS EN LLANURAS Tel: 54 - 291 - 459 5102 - int. 4360 Sergio Masuelli Fax: 54 - 291 - 459 5127 E-mail: ghgonzal@criba.edu.ar Pablo M. Vazquez Freddy Flores (Venezuela) Omar Faure .......................................................................................... 35 Fundación Instituto de Ingeniería Carretera Vieja de Baruta, Sector Sartenejas, Urb. Monte Elena II Caracas, Venezuela Tel: 58 2-903 4661-4610 Fax: 58 2- 903 4780 E-mail: freddyf@fii.org COMITÉ DE EVALUADORES Héctor del Valle (CENPAT - Argentina) Francisca González (UNS - Argentina) Mirta Aída Raed (CSR, UNLu - Argentina) Editado por: SELPER Internacional Walter Sione (UNLu - Argentina) Universidad Nacional de Luján Rutas 5 y ex 7, (6700) Luján - Bs. As. - ARGENTINA EJEMPLAR DE DISTRIBUCIÓN GRATUITA PROHIBIDA SU VENTA 5.
REVISTA VOL. 28 Nº 1 PRESENTACIÓN DEL PRESIDENTE DE SELPER04 JUNIO 2009 SELPER PLAN EDITORIAL SELPER 2008 - 2010 PLAN EDITORIAL SELPER A partir de las decisiones adoptadas en el marco del XIII Simpo- sio Latinoamericano de Percepción Remota y Sistemas de Infor- mación Espacial, llevado a cabo en La Habana, Cuba, en setiem- bre de 2008, la edición de la Revista SELPER está disponible en la página de nuestra Sociedad: http://www.selper.org. En esta oportunidad hacemos llegar la publicación del volumen 28 Número 1, donde se incluyen trabajos que han sido enviados por especialistas en el área de Teledetección y SIG que desarro- llan sus actividades en diferentes instituciones de América Latina y nos brindan la oportunidad de conocer los alcances de estas herramientas aplicadas al estudio del espacio geográfico y los recursos naturales. Estos trabajos han sido remitidos a expertos evaluadores para su revisión y aprobación. En este sentido es necesario destacar el invalorable apoyo brindado por el grupo de evaluadores, a quie- nes agradecemos muy especialmente por su labor. NORMAS PARA LOS AUTORES Los artículos recibidos serán enviados a tres (3) f) Resultados obtenidosexpertos en la temática para su revisión. Los g) Conclusionestrabajos aprobados serán publicados en estricto h) Bibliografía: sólo se incluirá la citada en elorden, de acuerdo a las fechas de llegada de las texto. Se indicarán los autores, por ordencontribuciones. alfabético, año, título, revista o medio donde fue Los idiomas oficiales SELPER son: Español, publicado, incluyendo volumen y páginas, cuandoPortugués e Inglés. corresponda. Los trabajos deberán estructurarsecontemplando las siguientes secciones: Los títulos y subtítulos de cada sección deberán e s ta r c l a r a m e n t e i n d i c a d o s ( y a s e a c o n a) Título del trabajo. Nombre de los autores y numeración o tamaño de letras). Las tablas, fotosdirecciones completas y figuras deberán ser suficientemente nítidas, b) Resumen (no más de 150 palabras) llevar un título y estar numeradas en formaindicando al final las palabras claves. Deberá consecutiva.incluirse en Español o Portugués, además de Se deberá enviar una copia del trabajo enInglés formato Word y una copia papel. La extensión c) Introducción total del trabajo no deberá superar las 12 páginas, d) Objetivos (DIN-A4). e) Metodología empleada y materiales Los trabajos se enviarán a: editorial@selper.org 6.
VOL. 28 Nº 1 REVISTA05 JUNIO 2009 SEGUIMIENTO DE LA REGENERACIÓN DE ÁREAS QUEMADAS.. SELPER Seguimiento de la regeneración de áreas quemadas a partir de imágenes hiperespectrales AHS Federico González-Alonso (alonso@inia.es) (1) Margarita Huesca (alonso@inia.es) (1) Sergio Martinez Sanchez-Palencia (2) (1) Centro de Investigación Forestal. Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria INIA, Crta. De la Coruña, km. 7,5, 28040 Madrid, Spain (e-mail: alonso@inia.es ) (2) Dirección General del Medio Natural. Junta de Comunidades de Castilla La Mancha. Toledo.RESUMEN INTA. This image is the most complete information from the affected area in Europe and likely in the world.En el presente trabajo se lleva a cabo un estudio en The main goal of this research is to burn areasprofundidad de la extensión de terreno afectada por el mapping and analyze forest fire severity as accurateincendio forestal que se produ jo en la provincia de as possible using hyperspectral data. Burned forestedGuadalajara (España) entre los días 16 y 21 de julio areas have patterns of varying burn severity as ade 2005, estudiando el grado de afección del mismo consequence of various topographic, vegetation, andcon el fin de extraer conclusiones para el futuro. El meteorological factors. These patterns are detectedfuego se inicio en el paraje la Cueva de Casares, en el and mapped using satellite data.municipio de la Riba de Saelice y tuvo muy graves An accurate burn severity map can have significantconsecuencias ecológicas y personales. La implications for defining management strategies. Theteledetección supone una herramienta potente y results of this study show that AHS images wereversátil, que suministra información muy valiosa y que accurate in order to analyze burnedpuede facilitar enormemente la toma de decisiones areas and to describe forest fire severity, whichpara la gestión forestal post -incendio. El objetivo proves the usefulness of AHS data for forest fireprincipal de este estudio es realizar una cartografía de area monitoring.superficie quemada y grados de severidad, delincendio de Guadalajara, a partir de datos INTRODUCCIÓNhierespectrales analizando la capacidad de estossensores para realizar un seguimiento de áreas Cada año se estima una pérdida de 10 a 15 millonesquemadas, para ello el Laboratorio de Teledetección de hectáreas de bosque de regiones boreales ydel CIFOR-INIA dispone en propiedad de una imagen templadas a causa de los incendios forestales, siendode la zona quemada captada por el sensor en España la principal causa de destrucción dehiperespectral aerotransportado AHS del INTA. La recursos naturales, solamente en el año 2006, másmencionada imagen hiperspectral constituye la de ciento cincuenta mil hectáreas fueron afectadasinformación más completa de la zona afectada de que por este fenómeno en nuestro país (MIMAM, 2006).se pude disponer actualmente en Europa y muy Los incendios forestales favorecen el cambioprobablemente en el mundo. Los resultados climático al emitir gases de efecto invernadero,obtenidos muestran la utilidad de los sensores concretamente dióxido y monóxido de carbono,hiperespectrales como el AHS-INTA para realizar metano y óxidos de nitrógeno y, a su vez, el cambioeste tipo de estudios. climático está incrementando el riesgo de incendios forestales, tendencia que seguirá aumentando en lasABSTRACT próximas décadas (Levine, 1991). El área Mediterránea es una de las zonas donde se espera unBurned areas monitoring using hyperespectral data efecto más significativo del cambio climático con unThe present research studies the potentiality of incremento de la aridez, lo que conllevaría unhyperespectral data to monitor burned areas. The aumento del número de incendios. De igual forma, elstudy area belongs to a coniferous forest area, fuego afecta notablemente al suelo a través delsituated in Guadalajara province (Spain) where a incremento de la temperatura y la deposición delarge forest fire occurred in July 2005. Several studies cenizas, la desaparición de la vegetación y de losshow remote sensing techniques provided an useful horizontes orgánicos total o parcialmentetool for forest fire monitoring which help decision modificando la acción de la lluvia sobre la superficiemaking management. An AHS image, which cover the del suelo, lo que podría ocasionar dinámicas deaffected area, has been acquired by the Remote degradación del mismo de una importante magnitud,Sensing Laboratory of CIFOR-INIA delivered from modificar el ciclo hidrológico y favorecer la 7.
VOL. 28 Nº 1 REVISTA06 JUNIO 2009 SEGUIMIENTO DE LA REGENERACIÓN DE ÁREAS QUEMADAS.. SELPERescorrentía. Estimar y evaluar la superficie afectada El presente estudio se centra en el incendio forestalpor un incendio de una forma rápida y precisa es de ocurrido en la provincia de Guadalajara (España)vital importancia a la hora de identificar áreas de entre los días 16 y 21 de julio de 2005, dicho incendiointervención prioritaria para evitar o disminuir en la tuvo muy graves consecuencias ecológicas ymedida de lo posible ese proceso degradativo del personales. La magnitud y gravedad de este incendiosuelo (González-Alonso et al., 2007). hace que su estudio presente un gran interés, ya queLa magnitud del impacto de un incendio forestal está podrían extraerse conclusiones que ayudarían adirectamente relacionada con la severidad del fuego comprender la interacción entre factores que provocódefiniéndose dicho término como la medida del nivel un desenlace trágico desde el punto de vista humanode daño ocasionado a la vegetación, materia y ecológico. Dichas conclusiones ayudarían a preverorgánica y suelo. La severidad del fuego es un y tratar de evitar situaciones de alto riesgo detérmino descriptivo que integra los cambios físicos, incendios en el futuro. El objetivo principal de estequímicos y biológicos que implica un incendio en el trabajo es estudiar la potencialidad de las imágenesecosistema (White et al., 1996). Tradicionalmente la captadas por el sensor hiperespectral aero-superficie afectada por un incendio así como la transportado AHS del INTA para evaluar la extensiónseveridad ocasionada por el mismo han sido del terreno afectado, así como el grado de severidadcartografiados con trabajo de campo o desde un ocasionado por el incendio, realizándose tambiénavión, el contorno exterior de la superficie afectada es una primera estimación de la regeneración de lahabitualmente delimitando mediante la utilización de vegetación sobre el área afectada, y con toda estaun GPS. Con la llegada de las imágenes de satélite la información extraer conclusiones para el futuro quesuperficie afectada por el incendio así como la pueda facilitar enormemente la toma de decisionesseveridad del mismo puede ser cartografiada de una para la gestión forestal post-incendio.forma rápida, económica y efectiva.El uso de la teledetección proporciona una cobertura ÁREA DE ESTUDIOespacial y temporal adecuada para este tipo detrabajos reduciendo el costoso trabajo de campo El presente estudio se centra en el incendio forestal(Chuvieco et al., 2004). La integración de datos que se produjo en la provincia de Guadalajara,espaciales en un Sistema de Información Geográfica situada en la zona central de la Península Ibérica,(SIG) permite almacenar y procesar una gran entre los días 16 y 21 de julio de 2005. El fuego secantidad de datos, así como llevar a cabo análisis inició en el paraje conocido como la Cueva deespaciales con información procedente de varias Casares, en el municipio de la Riba de Saelice, afuentes y con distintos formatos (Sunar y Özkan, causa de una barbacoa. . El área afectada cuenta con2001). Las imágenes de satélites han sido una orografía complicada, caracterizada por unampliamente utilizadas para la estimación de entramado de valles encajados con pendientessuperficies afectadas por este fenómeno como lo importantes, a esto hay que sumarle que durante losmuestran numerosos estudios (Merinode- Miguel et días del incendio las condiciones climáticas fueronal., 2005; González-Alonso et al., 2005, Barbosa et al., muy desfavorables, con una humedad relativa que1999; Martín et al., 2002), sin embargo la llegó a bajar del 20 por ciento, los vientos llegaron adiscriminación de distintos grados de severidad ser de 30 a 50 kilómetros por hora y la temperaturadentro de esa área estimada no ha sido tan estudiada. superior a 30 grados centígrados favoreciendoEn los estudios de severidad del fuego a partir de enormemente la propagación del incendio,datos satelitales se han empleado principalmente dificultando así notablemente la lucha contar elimágenes de media a alta resolución espacial, como mismo, según estudios realizados por expertos,las proporcionadas por los satélites, Landsat, Spot parece que dicho incendio se produjo bajo algunas deentre otros (Millar et al., 2007; De-Santis et al., 2007). las condiciones más adversas que se han podido darRecientes avances en la tecnología d e la en la zona desde que se tienen registrosteledetección han conducido al desarrollo de climatológicos. También es importante destacar lainstrumentos hiperespectrales, capaces de medir alta continuidad espacial de la vegetación presentesimultáneamente la radiación reflejada en una amplia en la zona, dominada por Pinus pin asater y Cistusgama de longitudes de onda. En estos sistemas, se laurifolius entre otras. Siendo la primera especie unaadquieren un gran número de bandas espectrales de las que presenta una mayor capacidad depara la misma área de la superficie terrestre, inflamabilidad. La superficie afectada alcanza lasfacilitando así discriminar rasgos de la cubierta que 12.875 hectáreas, de las que 3.147 pertenecen alno serían perceptibles en determinadas regiones del Parque Natural del Alto Tajo lo que representa un 2.9espectro. Estos sensores de alta resolución espectral por ciento de la superficie total del parque. La principalincluyen bandas muy estrechas, discriminando vegetación afectada corresponde a un pinar de pinorasgos de ab sorción que no están accesibles a resinero con un soto bosque de reb ollo y quejigo y porradiómetros más groseros. Hasta ahora no se han un matorral acompañante dominado por jara, esteparealizados muchos estudios que utilicen datos y brezinal. También es importante destacar la muertehiperespectrales en la cartografía de superficies de 11 personas pertenecientes a un retén al llevar aquemadas y grados de severidad (Parra-Uzcátegui, cabo las tareas de extinción del mismo.2005). 8.
VOL. 28 Nº 1 REVISTA07 JUNIO 2009 SEGUIMIENTO DE LA REGENERACIÓN DE ÁREAS QUEMADAS.. SELPERMATERIALEl laboratorio de Teledetección del CIFOR -INIAdispone en propiedad de dos imágenes de la zonaquemada captada por el sensor hiperespectralaerotransportado AHS del INTA (Instituto Nacional deTécnica Aeroespacial). El sensor AHS es un sensorhiperespectral aeroportado de barrido mecánico(whiskbroom scanner). El AHS está instalado a bordod e l a a e r o n a v e C A S A C - 2 1 2 d e l I N TA ,adecuadamente equipada para realizar campañas deteledetección. El AHS posee 80 bandas, cubriendolas regiones espectrales VIS/NIR, SWIR, MIR y TIR ycuyos anchos de banda varían según la región.El primer vuelo se realizó el 6 de octubre de 2005,antes de las primeras lluvias, por lo que refleja lasituación del terreno tras el paso del fuego. Lamencionada imagen hiperespectral constituye la Figura 2. Mosaico, en falso color, formado por las sieteinformación más completa de la zona afectada de la pasadas del vuelo AHS-INTA realizado elque se puede disponer actualmente en Europa y muy 6 de Octubre de 2005probablemente en el mundo. La campaña fuecompuesta por 7 pasadas paralelas, con rumbos que MÉTODOScoinciden con el plano solar de forma que se eviten oen todo caso minimicen posibles problemas de Trabajo de campoiluminación diferencial. Además los ejes de vuelo de Se han realizado cuatro visitas al área afectada por elpasadas se trazan para que las imágenes posean incendio del Rodenal, los días 6, 17 y 25 de octubre decierto solape que permitirá componer a posteriori un 2005 y el 1 de junio de 2006, cada una de ellas de unomosaico de imágenes (figura 2). El segundo vuelo se o más días de duración En la primera de estas visitasllevó a cabo el 6 de octubre de 2006, exactamente un se realizó un reconocimiento general del áreaaño después de la campaña anterior. afectada por el incendio y se tomaron datosPara mantener estables las condiciones de destinados a la elaboración de un muestreoadquisición de las imágenes se conservó la estadístico del área. En las siguientes visitas seconfiguración de vuelo de la campaña anterior, tomaron los datos de campo en las parcelas elegidasmanteniendo los rumbos de las pasadas y la hora de siguiendo el protocolo de toma de datos de campoadquisición de las imágenes de forma que se elaborado previamente. Debido a las limitaciones deminimicen los efectos de las diferencias en la tiempo y personal fue necesario plantear unailuminación que pueden producirse sobre las campaña de muestreo que permitiera evaluar lasimágenes. variables necesarias de una forma rápida, para poderEn la figura 1 puede verse la orientación y cobertura cubrir el máximo de superficie. Asimismo, fuede las pasadas y en la figura 2 se muestra el mosaico necesario estratificar la superficie de estudio segúnde la imagen del año 2005 utilizando una composición los criterios que se consideraron más convenientesen falso color. (tipo de vegetación, grado de afección, cercanía a carreteras) para muestrear las zonas de mayor interés. Se trató de cubrir con puntos de muestreo todas las combinaciones posibles de grado de severidad y tipo de vegetación previa al incendio. Se tomaron dos tipos de puntos de muestreo: Los Puntos A definidos como los localizados en zonas quemadas para estimar el grado de severidad causado por el incendio. Y los Puntos B, puntos tomados en zonas no quemadas, utilizados para definir zonas de entrenamiento (endmember) de vegetación sana y suelo desnudo (quemado y no quemado). Estos puntos se tomaron en los alrededores de la zona afectada por el incendio, considerando puntos que comprendían todos los tipos de formaciones vegetales presentes en la zona. Los datos tomados Figura 1. En esta figura están representados en azul las de cada parcela del grupo B fueron: una fotografía delpasadas realizadas, y en rojo la trayectoria del vuelo de toma aspecto de la zona en la que se localiza el punto,de imágenes hiperespectrales AHS realizado el 6 de octubre coordenadas GPS del punto de muestreo, de 2005 sobre la zona afectada por el incendio del Rodenal composición florística previa al incendio, tipo (Guadalajara). En amarillo aparece el contorno exterior estructural de vegetación previa al incendio, según la aproximado del incendio. Sobre cartografía E 1: 250 000 9.
VOL. 28 Nº 1 REVISTA08 JUNIO 2009 SEGUIMIENTO DE LA REGENERACIÓN DE ÁREAS QUEMADAS.. SELPERclave de Marchetti y la altura de hojarasca y mantillo Estimación de la severidad del incendio forestal(hasta llegar al suelo mineral). En los puntos del grupo La severidad del incendio ha sido estimada medianteA además de los datos anteriores se estimó el grado la utilización del índice “Normalized Burn Ratio”de severidad. basado en la combinación de las bandas del infrarrojoLa información de campo será utiliza para validar la cercano y el infrarrojo medio utilizando el siguientecartografía obtenida de la estimación de los grados de algoritmo:severidad del incendio. NBR = (? icr - ? im) / (?icr + ?im);Corrección radiométrica y geométrica de lasimágenes AHS Donde ?icr y ?im representan las reflectancias delLas imágenes AHS de los años 2005 y 2006 han sido infrarrojo cercano e infrarrojo medio respectivamente.corregidas geométricamente utilizando un método de La banda del infrarrojo cercano es sensible a lacorrelación con ventanas móviles a lo largo de las actividad fotosintética de la vegetación mientras queimágenes para evaluar así su ajuste geométrico y la del infrarrojo medio proporciona información delpermitir una corrección geométrica con una precisión contenido de agua esto hace que el índice NBR seamayor que la conseguida con los métodos especialmente sensible a la detección de cambiostradicionales (Prado-Ortega, 2007). Reducir los entre la vegetación viva y el contenido de humedaderrores en la geolocalización de los píxeles es de gran así como de algunos cambios producidos en el sueloimportancia en los estudios multitemporales. El paso como consecuencia de la actividad del fuego (Key yde radiancias a reflectancias fue realizado con la Benson, 2005). El índice NBR fluctúa entre -1 y +1.ayuda del software ATCOR4 para cada una de las Valores bajos de índice indican severidad altapasadas que comprende el vuelo del 2005. El cálculo mientras que valores altos severidad baja ode reflectancias para las pasadas del vuelo del año inexistente.2006 se ha realizado mediante un análisis estadístico Para una mejor comprensión de los resultadosde las radiancias y reflectancias del 2005 y su obtenidos, la imagen será clasificada utilizando unaposterior aplicación a las pasadas del vuelo del 2006. clasificación no supervisada con el algoritmo IsoData, en tres grados de severidad, severidad alta, media yEstimación de la superficie afectada por el baja. La validación de dicha cartografía se realizaráincendio forestal mediante un análisis por matrices de confusión entrePara la delimitación del perímetro y la estimación de los resultados obtenidos y las parcelas del trabajo dela superficie afectada se realizó un análisis de campo anteriormente descritas.“Matched Filtering” sobre la imagen para enfatizar lazona quemada, este método asigna a cada píxel un Primera estimación de la regeneración de lavalor de probabilidad de pertenecer a una clase, dicha vegetaciónclase se ha definido con un conjunto de píxeles de El estudio de la regeneración se ha realizadoentrenamientos o “endmember”, dicho endmember mediante el uso del índice de vegetación dese definió meditante la selección 100 píxeles de la diferencia normalizada NDVI (Díaz-Delgado et al.,imagen que presentaran un grado de quemado 1998), dicho índice se basa en la combinación de lasconsiderable. A continuación se realizaron unas bandas roja e infrarroja cercana del espectroclasificaciones no supervisadas utilizando el electromagnético. El NDVI detecta la actividadalgoritmo Isodata, con dos clases quemado y no- fotosintética o grado de verdor existente en un píxelquemado, modificándose los umbrales inferiores del permitiendo así el estudio del desarrollo de lahistograma hasta alcanzar la máxima precisión vegetación. NDVI varía entre -1 y 1, cuanto másposible. Este método ha sido seleccionado por los elevado sea el índice mayor será la actividadbuenos resultados obtenidos en trabajos anteriores fotosintética registrada en el píxel y de forma contraria(Merino-de-Miguel et al., 2005). cuanto más bajo sean los valores de NDVI mayorLos resultados obtenidos con las imágenes AHS presencia de suelo desnudo encontraremos. El NDVIserán posteriormente comparados con los resultados esta formulado de la siguiente forma:de aplicar la misma metodología a imágenesmultiespectrales con distintas resolucionesespaciales. NDVI = (?irc – ?r) / (? irc + ? r); Donde ?icr y ?r representan las reflectancias del infrarrojo cercano y rojo respectivamente. Dichas bandas están centradas a las longitudes de onda de 0.8 ìm la del infrarrojo cercano y a 0.6 ìm la del rojo. El grado de regeneración ha sido calculado utilizando como algoritmo la diferencia de NDVI entre los años 2006 y 2005. Figura 3. Mosaico AHS del 6 de octubre de 2005, en falsocolor, y superficie clasificada como superficie afectada por el incendio en rojo 10.
VOL. 28 Nº 1 REVISTA09 JUNIO 2009 SEGUIMIENTO DE LA REGENERACIÓN DE ÁREAS QUEMADAS.. SELPERLa siguiente tabla muestra un resumen de la AWiFS y AHS y la relación con los datos oficialessuperficie afectada por dicho incendio estimada a proporcionados por el Ministerio de Medio Ambiente.partir de las imágenes MODIS, MERIS, LANDSAT, Tabla 1. Superficie afectada (ha) por el incendio del Rodenal en GuadalajaraSe puede observar como la superficie estimada conlas imágenes AHS se acerca con mayor precisión alas cifras proporcionadas por el Ministerio de MedioAmbiente, quedando así justificada 9 la utilización delos datos hiperespectrales para la estimación delperímetro así como del área afectada tras un incendioforestal.Estimación de la severidad del incendio forestalLa estimación de los grados de severidad se harealizado a partir del índice NBR, explicadoanteriormente, como consecuencia de la altaresolución espectral de las imágenes AHS, antes delcálculo de dicho índice se realizó un análisis para ladeterminación de las bandas que mejor Figura 4. Cartografía de los grados de severidad delrepresentaban las zonas del espectro infrarrojo incendio del Rodenal en Guadalajara, a partir del índicecercano e infrarrojo medio que intervienen en el NBR1cálculo de dicho índice. El criterio elegido para dicha En la figura 4 se puede observar la distribuciónselección de bandas, esta basado en aquellas que espacial de la severidad del incendio del Rodenalpresenten un menor ruido y una máxima variabilidad. ocurrido en Guadalajara en el 2005. La severidad seLos resultados obtenidos muestran que la banda que ha definido como el impacto del fuego sobre lamejor representa la zona del espectro del infrarrojo vegetación distinguiéndose tres grados de severidadcercano es la banda número 13 centrada a una de incendio. Severidad baja representada por loslongitud de onda de 0.8 ìm y para el infrarrojo medio colores blancos, severidad media con tonos amarillasla s bandas 31, 38 y 41 centradas a unas longitudes y finalmente los tonos anaranjados representan lasde onda de 2.13, 2.22 y 2.25 ìm respectivamente. De zonas de severidad alta. Se puede observar que laesta forma se calcularon tres índices de severidad zona más afectada por dicho incendio se sitúa en unadefinidos por NBR1, NBR2 y NBR3. Se realizaron de franja central del área afectada con una orientaciónesta forma tres clasificaciones de grados de suroeste – noreste. La severidad predominante es laseveridad distintas correspondientes a los tres clasificada como severidad media representada poríndices de severidad calculados. Utilizando el trabajo 7.286,0841 ha seguida por las zonas de severidadde campo como validación de la cartografía obtenida alta con 2.544,8490 ha, siendo la superficiese comprobó que la clasificación más precisa se clasificada como severidad baja la que comprendealcanzaba con el índice de severidad NBR1 calculado una menor superficie con 3.152,5740 ha. Analizandoa partir de la banda 31 del infrarrojo medio centra da a el modelo digital de terreno se puede observar la gran2.13 ìm., obteniéndose una fiabilidad del 70%. dependencia entre grados de severidad y relieve. LasLa siguente figura muestra la cartografía de los zonas de mayores pendientes por lo generalgrados de severidad obtenida utilizando el índice muestras un grado de severidad mayor que las zonasNBR1. más llanas. La información obtenida de la cartografía presentada, puede ayudar en gran medida al gestor en los trabajos de restauración de la zona afectada, así como en la planificación de futuros planes de prevención de incendios. 11.
VOL. 28 Nº 1 REVISTA10 JUNIO 2009 SEGUIMIENTO DE LA REGENERACIÓN DE ÁREAS QUEMADAS.. SELPERPrimera estimación de la regeneraciónLa siguiente figura representa la diferencia de NDVIentre los años 2006 y 2005 Figura 5. Cartografía de la primera estimación de la Figura 6. Cruce de los grados de severidad del incendio con regeneración en el periodo 2005-2006 en el incendio del los grados de regeneración de la vegetación Rodenal de Guadalajara Las zonas rojas son aquellas donde coincide una severidad alta con una regeneración baja, luego son En la figura 5 se puede observar la distribución las zonas de prioridad a la hora de realizarespacial del grado de regeneración observado en el actuaciones de mejora. Dichas zonas de prioridadperiodo 2005 -2006, calculado a partir del índice de están representadas en la figura 6 mediante unavegetación NDVI. Para una mejor interpretación de elipse negra, el área 1 ocupa una superficie de 1.531los resultados se han definido tres grados de ha, el área 2, 376 ha y el área 3 ocupa una superficieregeneración representados por tonos verde oscuro de 124 ha. Esta cartografía facilita en gran medida lalas zonas con una regeneración alta, los verdes claro toma de decisiones para la gestión forestalrepresentan zonas donde se ha producido una postincendio.regeneración media y las zonas donde no ha habidoregeneración o esta ha sido muy baja están CONCLUSIONESrepresentadas con tonos amarillos. Las zonas con ungrado de regeneración mayor se sitúan en la parte Tras producirse un incendio de las características delnorte del área afectada mientras que las menos incendio de Guadalajara surge la necesidad deregeneradas se presentan en la zona sur definiendo conocer la extensión de terreno afectado, el grado deuna franja con orientación noroeste – sureste. Se afección, y de estudiar en profundidad lo ocurrido, conpuede observar que la mayor parte de la zona el fin de extraer conclusiones para el futuro. Esafectada por el incendio presenta un grado de asimismo prioritario en los años siguientes analizarregeneración bajo (5.219,9105 ha) o medio cómo avanza la regeneración, con el fin de gestionar(4.850,0069 ha), mientras que las zonas con una la masa de la forma más adecuada. Para todas estasregeneración alta representan una superficie de cuestiones que se les plantean a los gestores2.913,5826 ha. forestales en las zonas afectadas por incendios, elLa opinión de los expertos ha permitido confirmar el presente estudio muestra que la teledeteccióngrado de precisión obtenido en la cartografía supone una herramienta potente y versátil, queanteriormente expuesta. El trabajo de campo que suministra información muy valiosa ya que puedeserá proporcionado por la Junta De Castilla La facilitar enormemente la toma de decisiones para laMancha en breve, permitirá una validación con un gestión forestal post-incendio.mayor rigor científico. Respecto a la capacidad de las imágenesFinalmente se realizó un cruce entre la severidad del hiperespectrales AHS para discriminar la superficieincendio y los grados de regeneración de la afectada por el incendio de Guadalajara podemosvegetación para localizar las zonas de principal señalar la gran precisión de los resultadosactuación, siendo estas representadas por aquellas comparados con los obtenidos con otros sensoreszonas donde la severidad del incendio ha sido multiespectrales, lo que confirma nuestra hipótesis declasificada como alta y la regeneración del mismo partida de la utilidad de este tipo de imágenes para labaja o prácticamente nula. La siguiente figura obtención de cartografías de áreas quemadas..representa el cruce de ambas cartografías Los resultados obtenidos en el estudio de laseñalándose las tres zonas de principal actuación severidad, nos lleva de igual forma a confirmar la(figura 6). capacidad de este sensor para discriminar entre diversos grados de severidad provocados por un incendio forestal. La clasificación de una escena 12.
VOL. 28 Nº 1 REVISTA11 JUNIO 2009 SEGUIMIENTO DE LA REGENERACIÓN DE ÁREAS QUEMADAS.. SELPERposterior al incendio, utilizando el índice de severidad conjunto de todos estos factores permitiríaNBR1, ofrece unos resultados bastante fiables comprender mejor el comportamiento del fuego y susdemostrando la utilidad de la teledetección para la efectos, y extraer conclusiones aplicables a la gestiónobtención de este tipo de cartografías aportando una de futuros incendios. Se considera también de granherramienta de apoyo al gestor de una forma fiable, interés realizar un seguimiento riguroso de larápida y económica comparada con los métodos regeneración que se produzca en la zona durante lostradicionales de trabajo de campo y foto- próximos años, empleando imágenes AHS y datosinterpretación. tomados de campo. Esto permitirá analizar las pautasConclusiones que pueden ser extrapoladas al estudio de regeneración de las distintas especies, yque se ha realizado de estimación de los grados de relacionarlas con los grados de severidad y laregeneración de la vegetación. temperatura del suelo tras el paso del fuego.Finalmente, se considera esencial explotar lainformación contenida en la imagen AHS, con el fin de AGRADECIMIENTOSobtener cartografía detallada de los distintos gradosde severidad causados por el incendio, de la El presente trabajo se ha realizado en el marco deltemperatura del suelo (muy influyente en la Convenio CC06-021 firmado entre el INIA y la Juntaregeneración natural) y conocer la relación entre la de Comunidades de Castilla La Mancha para elseveridad producida y la biomasa acumulada en la estudio del incendio del rodenal de Guadalajara.zona entre los Inventarios Nacionales Forestales Los autores agradecen al Instituto Nacional deSegundo y Tercero, así como la relación entre Técnica Aeroespacial la ayuda prestada en laseveridad y especies presente antes del incendio. adquisición y corrección de las imágenesEstos datos podrán combinarse también con aquellos hiperespectrales AHS empleadas en la realización derelativos a la evolución y características del incendio este trabajo.que se encuentran disponibles en los ServiciosForestales de Castilla – La Mancha. El estudio 13.
VOL. 28 Nº 1 REVISTA12 JUNIO 2009 SEGUIMIENTO DE LA REGENERACIÓN DE ÁREAS QUEMADAS.. SELPERREFERENCIAS Implications edited by J. S. Levine (Cambridge, USA: MIT press), pp. xxv – xxx.BARBOSA. P.M., GREGOIRE, J.L. y PEREIRA, J.M.1999. An algorithm for extracting burned areas from MARTÍN, M.P., DÍAZ DELGADO, R., CHUVIECO y G.time series of AVHRR GAC data applied at a VENTURA, G. 2002. Burned land mapping usingcontinental scale. Remote Sensing of Environment. NOAA-AVHRR and TERRA-MODIS. Forest Fire69: 253-263. Research & Wildland Fire Safety (Viegas. Ed). Millpress, Rotterdam.CHUVIECO, E. y MARTÍN, M. P. 2004. NuevasTecnologías para la estimación del riesgo de MERINO-DE-MIGUEL, S., GONZÁLEZ-ALONSO, F.,incendios forestales. Consejo Superior de GARCÍA-GIGORRO, A., ROLDÁN-ZAMARRÓN, A. yInvestigaciones Científicas. Instituto de Economía y CUEVAS, J.M. 2005. Is it possible to timely andGeografía. Colección de estudios ambientales y accurately estimate wildfire burnt areas using remotesocioeconómicos. sensing techniques? En: New Strategies for European Remote Sensing (Oluic, Ed). Millpress,DE SANTIS, A. y CHUVIECO, E. 2007. Burn severity Rotterdam.estimation from remotely sensed data : Performanceof simulation versus empirical models. Remote MILLER, J. D. y THODE, A. E. 2007. Quantifying burnSensing of Environment, 108: 422-435. severity in a hete rogeneous landscape with a relative version of the delta normalized burn ratio (DNBR).DÍAZ-DELGADO, R., SALVADOR, R., y PONS, X. Remotes Sensing of Environment. 109: 66-80.1998. Monitoring of plant community regenerationalter fire by remote sensing. In: L. Traboud (Ed.), Fire MIMAM. 2006. Ministerio de Medio Ambientemanagement and landscape ecology. Fairfield, WA: (España), Estadística incendios forestales.International Association of Wildland Fire. pp 315-324 http://www.mma.es/portal/secciones/biodiversidad/d efensa_incendios/estadisticas_incendios/in dex.htmGONZÁLEZ-ALONSO, F., MERINO-DE-MIGUEL, S., (accessed on March 2006).ROLDÁN-ZAMARRÓN, A., GARCÍA-GIGORRO, S. yCUVAS, J. M. 2007. MERIS Full Resolution data for PARRA UZCÁTEGUI, A. 2005. Estimación de nivelesmapping level-of-damage by forest fires: the Valencia de severidad del fuego a partir de datosde Alcántara event in August 2003. International hiperspectrales E0-1 Hyperion. Tesis doctoral,Journal of Remote Sensing. 28: 797-809. Departamento de Geografía, Universidad de Alcalá de Henares, Madrid, España.GONZÁLEZ-ALONSO, F., MERINO-DE-MIGUEL, S.,ROLDÁN-ZAMARRÓN, A., GARCÍA-GIGORRO, S. y PRADO ORTEGA, E. 2007. Mejoras en la evaluaciónCUEVAS, J.M. 2005. Algoritmos de estimación de y corrección de errores de corregistro en imágenessuperficies quemadas a partir de compuestos MODIS. hiperspectrales AHS para estudios multitemporales.XI Congreso Nacional de Teledetección, 21-23 Trabajo de investigación tutelada, Departamento deseptiembre 2005. Puerto de la Cruz. Tenerife. Geografía, Universidad de Alcalá de Henares, Madrid, España.KEY y BENSON. 2005. C. H. Key and N. C. Benson,Landscape assessment: Remote sensing of severity, SUNAR, F. y ÖZKAN, C. 2001. Forest fire analysisthe Normalized Burn Ratio. In: D. C. Lutes et al., with remote sensing data. International Journal ofEditors, FIREMON: Fire effects monitoring and Remote Sensing, 12: 2265 – 2277.inventory system, General Technical Resport, RMRS-GTR-164-CD:LA1-LA51, USDA Forest Service, WHITE, J.D., RYAN, K. C., KEY, C. C. y RUNNING,Rocky Mountain Research Station, Ogden, UT (2005). S.W. 1996. Remote sensing of forest fire severity and vegetation recovery. International Journal of WildlandLEVINE, J. S. 1999. Introduction In Global Biomasa Fire, 6: 125-136Burning: Atmospheric, Climatic and Biospheric 14.
VOL. 28 Nº 1 REVISTA13 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPER Evaluación del impacto del cambio de uso de suelo en la hidrología de la cuenca del Río San Juan en Nuevo León, México. Axayácatl Maqueda Estrada 1, Jianhong Ren 2, Diego Fabián Lozano García 3 1 Investigador del Instituto Tecnológico y de Estudios Superiores de Monterrey, Av. Eugenio G. Sada 2501, Monterrey, México, CP 64849 axa.maqueda@itesm.mx Tel. (52) 8183582000 ext 5288 2 Investigadora de la Universidad Texas A&M Kingsville, Environmental Engineering Department, 700 University Blvd. MSC 367, Kingsville, TX, USA. CP 78363 jianhong.ren@tamuk.edu Tel (361) 5932798 3 Coordinador del Laboratorio de Sistemas de Información Georreferenciada del ITESM Monterrey, Av. Eugenio G. Sada 2501, Monterrey, México CP 64849 dflozano@itesm.mx Tel 8183582000 ext 5288.RESUMEN ABSTRACTLa relación entre el tipo de cobertura del suelo de una This project evaluates the effect of potential landcovercuenca y su comportamiento hidrológico ha sido change on surface runoff generation in the San Juandescrita por muchos autores de diversas River Watershed caused by population growth.publicaciones. El tipo de cobertura del suelo, los A GIS based computer model, the Automatedbosques en particular, son uno de los factores de Geospatial Watershed AssessmentTool (AGWA),mayor influencia para la hidrología de una cuenca. was used. Historical land cover type and river flowLos bosques incrementan la infiltración del agua en el data available for the San Juan River watershed weresuelo, este proceso reduce el escurrimiento used as model inputs and for model calibration,superficial por lo tanto reduce el potencial de respectively. The computer model discretizes theinundación. La conversión de áreas boscosas a áreas watershed into six small units and predicts theurbanas o para la agricultura incrementará el watershed hydrology for future land coverescurrimiento de los ríos. Este artículo presenta un conditions. Simulation results show minor changes incaso de estudio conducido a nivel cuenca. El water flows at the watershed scale through the studypropósito de este proyecto es evaluar el efecto del period because subwatersheds that have land coverpotencial cambio de cobertura del suelo en la cuenca change make up only 23.5% of the study area. All thedel Río San Juan causado por el crecimiento de las impact is concentrated in two subdivided regions (theáreas urbanas. El modelo utilizado fue AGWA urban areas and the Huajuco Canyon regions).(Automated Geospatial Watershed Assessment Tool) In the urban areas, there is an increase up to 63.5% inbasado en SIG. Datos históricos de cobertura de surface flow when the area covered by the impervioussuelo y escurrimientos del Río San Juan fueron surfaces grows 18.7% by 2020. For theutilizados como información para alimentar el modelo. Huajuco Canyon, there is a similar tendency. FloodEl modelo permitió dividir la cuenca en seis regiones e hazard may increase due to increases in surfaceintentar predecir los escurrimientos de agua en base runoff during the wet season. There may be a locala condiciones futuras de cobertura del suelo. Los scarcity of surface water due to the dramatic decreaseresultados del modelo muestran que sólo hay in percolation and baseflow during the dry season incambios pequeños en los escurrimientos a nivel the Huajuco Canyon and urban areas. Thecuenca durante el periodo estudiado ya que sólo se land cover change effect during an extreme rainfallhan transformado subcuencas que cubren el 23.5 % event was evaluated using the rainfall event thatdel área de estudio. Sin embargo, todo el impacto occurred during Hurricane Gilbert in 1998. Simulationestá concentrado en dos subcuencas (Áreas Urbanas results show that land cover has a minimum influencey Cañón del Huajuco). on the volume of surface flow produced by this rainfallEn la subcuenca de Áreas Urbanas se predice un event.incremento de hasta 63.5% en el escurrimientosuperficial cuando el área cubierta por superficies INTRODUCCIÓNimpermeables aumente en 18.7% para el año 2020.Para el Cañón del Huajuco hay una tendencia similar. La degradación del medio ambiente en México se haPor lo tanto, el peligro de inundación aumenta debido vuelto una preocupación a nivel nacional. Esteal incremento del escurrimiento superficial durante la problema puede afectar la gobernabilidad de latemporada de lluvias y puede haber una escasez sociedad. Problemas como la erosión acelerada dellocal de agua debido a la caída de la infiltración y flujo suelo, la deforestación, la escasez de agua y labase del río durante la época de estiaje en las pérdida de biodiversidad son considerados ahorasubcuencas de Áreas Urbanas y Cañón del Huajuco. como el origen de muchos conflictos sociales en México. Debido a estos conflictos, los temas 15.
VOL. 28 Nº 1 REVISTA14 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPERrelacionados con el manejo de bosques y recursos demostrada por Descheemaeker et al. (2006). Estehídricos son considerados como asuntos de artículo presenta un caso de estudio a nivel cuenca.seguridad nacional (Cotler et al., 2004). La cuenca del El propósito de este proyecto fue evaluar el efecto delRío San Juan comprende un área de 32,789 cambio de cobertura de suelo en la hidrología de lakilómetros cuadrados y está localizada en los estados cuenca (infiltración, escurrimiento superficial, y flujode Nuevo León y Coahuila como se muestra en la base del río). El modelo basado en SIG (sistemas defigura 1. El principal usuario del agua de la cuenca es información geográfica), AGWA fue utilizado. Datosla ciudad de Monterrey. El área metropolitana de la históricos del tipo de cobertura de suelo yciudad tiene ahora cerca 4 millones de habitantes y se escurrimientos del Río San Juan fueron utilizadosespera que haya 5.1 millones en el año 2020. Para como entradas y calibración para el modelo.cumplir con las necesidades de la población crecientees necesario desarrollar al menos 25,000 hectáreas ENFOQUE Y METODOLOGÍAde nuevas áreas urbanas adicionales a las actuales57,500 hectáreas (Guajardo, 2003). El crecimiento de Descripción del área de estudiolas áreas urbanas implica la transformación de El área de estudio es una fracción de la cuenca delmatorral, bosque o zonas agrícolas en materiales Río San Juan en México como se muestra en laimpermeables como concreto y asfalto. Una Figura 1. El Río San Juan representa el 6.35% delimportante correlación entre el tipo de cobertura de área total de drenaje de la cuenca del Río Bravo. Lasuelo en una cuenca y su hidrología ha sido descrita extensión del área de estudio es de 3,488 kilómetrosen varias publicaciones (e.g. Taniguchi, 1997). La cuadrados. La lluvia promedio anual es de 414 mmcobertura del suelo, particularmente los bosques, es (CNA, 2005). 57% del area de la cuenca tiene ununo de los factores que más influyentes en la clima semicálido con una temperatura promediohidrología de una cuenca. anual entre 18 y 22 grados Celsius. La otra parte de laLos bosques permiten una alta infiltración de agua en cuenca (43%) tiene un clima semiárido con unael suelo, que minimiza el escurrimiento superficial y temperatura promedio anual superior a los 22 gradostransforma el potencial escurrimiento superficial en Celsius. La escasez de agua en la cuenca es causadaagua subterránea evitando inundaciones. El agua princi palmente por factores climáticos como altossubterránea fluye más lentamente hacia manantiales niveles de evaporación debido a altas temperaturasy corrientes principales, y evita que las corrientes promedio anuales. Además, la mayor parte de la lluviasuperficiales se sequen durante el estiaje ocurre entre Julio y Octubre cuando se presentan las(Kostadinov et al., 1998). Un incremento significante temperaturas más altas del año. Por lo tanto el aguaen escurrimientos de ríos causado por la conversión es escasa a pesar de que hay lluvias acumuladasde bosques a áreas agrícolas ha sido observado por anualmente relativamente altas, de 600 a 1200 mmvarios investigadores previamente (e.g. Heil et al. (INEGI, 1990). Vientos dominantes del Golfo de2003; Williamson et al., 1987). Una reducción México chocan con las montañas de la cuenca ysignificativa del escurrimiento superficial después del concentran la mayor parte de la lluvia en la cara esteincremento de cobertura vegetal del suelo fue de la zona montañosa (Figura 2). Figura 1. Área de estudio Figura 2. Distribución de la lluvia en el área de estudio y de los índices utilizados.El área de estudio fue dividida en seis subcuencas cubre 527.8 km2 o 15% del área de estudio, toda estaque se nombraron Zona Agrícola, Bosque, Área subcuenca está en una zona montañosa con alturasSuburbana, Cañón del Huajuco, Cañón de la superiores a los 3,000 m.s.n.m. A pesar de existirHuasteca y Áreas Urbanas. La Figura 3 muestra el pequeñas áreas agrícolas la mayor parte de laárea de estudio por subcuenca y de fondo una imagen subcuenca está cubierta por bosques. El Cañón delLandSAT de 1998. La Zona Agrícola que está cubierta Huajuco está al sur de Monterrey y abarca 370 km 2 opor pastizales y huertas citrícolas cubre 694 km 2 o 10.5% del área de estudio y se encuentra en un19.7% del área de estudio. La mayor parte de esta proceso de acelerada urbanización. La cobertura delsubcuenca son áreas planas a excepción de algunas suelo es bosques en las montañas y matorral en losmontañas en su extremo oeste. La región de Bosque valles. Las nuevas áreas urbanas se encuentran en 16.
VOL. 28 Nº 1 REVISTA15 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPERlos valles. El Cañón de la Huasteca es una subcuenca en pendientes pronunciadas. La subcuenca de Áreasaislada por las altas montañas, por lo tanto la lluvia es Urbanas es pequeña comparada con el resto del áreaconsiderablemente menor que en las otras de estudio, comprende solamente 445 km 2 o 13% delsubcuencas. Abarca 1,145 km 2 o 33% del área de total y está cubierta principalmente por la ciudad deestudio. Esta subcuenca está cubierta en su totalidad Monterrey. La cobertura del suelo en esta área espor montañas. La cobertura del suelo es escasa, principalmente asfalto y concreto. La subcuenca dematorral desértico en pendientes pronunciadas. Esta Área Suburbana cubre sólo 9% del total y estásubcuenca está cubierta en su totalidad por montañas. formada por desarrollos urbanos de baja densidad yLa cobertura del suelo es escasa, matorral desértico zonas agrícolas. Figura 3. Imagen LandSAT del área de estudio y su división en 6 subcuencasEl modelo de lluvia-escurrimiento de lluvia corresponden a 10 estacionesLa herramienta AGWA 1.5 fue utilizada en este climatológicas instaladas en el área de estudioestudio (Miller et al., 2002). AGWA contiene a su vez (Figura 4). La lluvia se distribuyó en la cuenca deel modelo SWAT (Soil & Water Assessment Tool, the acuerdo al criterio de Thiessen (Miller et al., 2002).SWAT model) (Arnold, 2000). SWAT es un modelo Los registros del caudal del Río San Juan de 1976 adiseñado para predecir los efectos de prácticas de 1996 fueron extraídos del sistema BANDAS (IMTA a,manejo de la tierra (agricultura, pastoreo, desarrollo 1999).urbano) en la hidrología de cuencas grandes condiferentes tipos de suelo y por periodos más grandesa un año. Es un modelo que utiliza la lluvia acumuladadiaria, el tipo de suelo, el tipo de cobertura de suelo yla topografía como entrada. El modelo SWAT produceresultados anuales y mensuales para infiltración,escurrimiento superficial, flujo base y flujo total de lacuenca (Miller et al., 2002). Utiliza el método delNúmero de Curva desarrollado por el SoilConservation Service del Departamento deAgricultura de E.U.A. para estimar el escurrimientosuperficial en función de la lluvia, el suelo y el tipo decobertura del suelo (USDA, 1986).Datos de lluvia y caudal del ríoLa lluvia acumulada diaria del año 1970 al año 2000fue utilizada como entrada el modelo. Estos datos se Figura 4. Localización de las estaciones climatológicasobtuvieron del sistema ERIC (IMTA, 1999). Los datos usadas en el estudio, cuadrícula 30x30 km 17.
VOL. 28 Nº 1 REVISTA16 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPERDatos de suelo y 5 fueron realizadas para evaluar el efecto delLos mapas de suelo producidos por el INEGI (INEGI, cambio de la cobertura del suelo en la hidrología de la2000) escala 1:250,000 fueron desarrollados cuenca. La lluvia promedio diaria del año 1970 al añoutilizando la clasificación de la FAO (FAO, 1974). 182 2000 y la cobertura del suelo de los años 1974, 1998 ydiferentes tipos de suelo se encuentran en el área de 2020 fueron utilizadas como entradas para la corridasestudio de acuerdo a la clasificación de la FAO. Estos 3, 4 y 5 del modelo SWAT.tipos de suelo fueron reclasificados en cuatro gruposhidrológicos (A, B, C, D) de acuerdo al trabajo de RESULTADOSLevick et al. (2002) e INEGI (2000). Debido a la escaladel mapa de suelos puede haber un error en el Grupos hidrológicos del suelomodelo, la magnitud de este error no es cuantificable. La Figura 5 muestra la reclasificación de los tipos de suelo. 54.4% del área de estudio está formada porCobertura de suelo suelos tipo A, que están formados por gravas y arenasExiste información de cobertura del suelo para los que tienen una alta tasa de infiltración. La mayor parteaños 1974 y 1998. La clasificación de la cobertura del de los suelos tipo A se encuentran en las zonassuelo en México fue desarrollada considerando la montañosas. Los suelos tipo A permiten la infiltraciónapariencia de la vegetación (Moreno, 2002). La del agua en el suelo de la cuenca. Por lo tanto soninformación de la cobertura del suelo fue reclasificada muy importantes para el flujo base de los ríos. Losbajo los criterios desarrollados por el SCS (USDA, suelos tipo D son compuestos por arcillas de textura1986). Ya que esta clasificación realizada por INEGI fina tienen la menor capacidad de absorber agua yen México provee solamente información sobre el tipo cubren el 24.5% del área de estudio. Se encuentrande vegetación y no sobre el uso actual del suelo en zonas con poca pendiente. Los suelos tipo B y Cpuede haber errores que no se pueden estimar en formados por limos cubren el 11.2 y el 9.7% del áreaesta reclasificación. La cobertura del suelo para el de estudio respectivamente.año 2020 fue estimada utilizando como base losresultados del Análisis Estratégico del Área Cambios en la cobertura del suelo durante elMetropolitana de Monterrey (Guajardo, 2003). En periodo de estudioeste estudio se estima que la población crecerá 30% y La figura 6 muestra los datos de cobertura del sueloel área urbana crecerá 40% del año 2000 al 2020. de los años 1974, 1998 y 2020. El porcentaje de contribución de cada tipo de cobertura al área total seSimulaciones utlilizando el Modelo SWAT muestra en la figura 7. Del año 1974 a 1998 hubo unAntes de calibrar el modelo SWAT primero se decremento del 5.1% en el matorral desértico en eldeterminó el porcentaje de escurrimiento superficial y área de estudio. En el año 1974 sólo 2.2% de laflujo base del río San Juan utilizando registros de una cuenca estaba cubierta por áreas urbanas. Para elestación hidrométrica a la salida de la cuenca. Se año 1998 6.8% de la cuenca eran áreas urbanas yutilizó un módulo del modelo SWAT para hacer esta para el 2020 10.9% serán áreas impermeables.separación (Arnold, 1995). Los flujos diarios de 1970 Entonces, los cambios más notables en la coberturaa 1994 fueron utilizados para calcular el flujo base. es la reducción del matorral desértico y el incrementoDespués, se realizaron seis corridas del modelo de las áreas urbanas. El crecimiento del área urbanaSWAT. La primera fue para calibrar el modelo, la lluvia es mostrado en la Figura 8, donde las áreas amarillasdiaria, el flujo base del río, y la cobertura del suelo fueron pastizales, matorral desértico o zonaspara el año 1974 fueron utilizados como entradas del agrícolas en 1998 y para el 2020 muy probablementemodelo. La calibración se realizó de acuerdo a lo sean áreas urbanas. Las áreas verdes fueronrecomendado por los desarrolladores del modelo matorral desértico en 1998 y se espera que seanSWAT, se modificaron algunos parámetros de bosques en el año 2020. El área convertida en áreasentrada (Neitsch et al., 2000). La segunda corrida del urbanas en el año 2020 es el 18.7 % del área de lamodelo fue para verificar la calibración del modelo y subcuenca Áreas Urbanas. Para la subcuenca Cañónlas coberturas del suelo de los años 1974 y 1998 del Huajuco el 19.7% de su área será convertida defueron utilizadas para simular los flujos del río para los pastizal o matorral a áreas urbanas para el año 2020.años 1976 y 1994 respectivamente. Las corridas 3, 4 18.
VOL. 28 Nº 1 REVISTA17 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPER Figura 5. Tipos hidrológicos de suelo en el área de estudio, cuadrícula de 20x20 km Figura 6. Tipos de cobertura de suelo para los años A) 1974, B) 1998 y C) 2020. Cuadrícula de 25x25 km 19.
VOL. 28 Nº 1 REVISTA18 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPER Figura 7. Contribuciones porcentuales de cada tipo de cobertura de suelo al área total de la cuenca para los años A) 1974, B) 1998, y C) 2020 Figura 8. Mayores Cambios en la cobertura del suelo del año 1998 al añ0 2020Variación de la lluvia para 1975, 489 mm para 1998 y 632 mm para 1999. ALos datos de lluvia están disponibles del año 1970 al causa de las variaciones entre los años, unaaño 2000 para todas las estaciones climatológicas. simulación utilizando los datos de lluvia de 1974 noLa lluvia acumulada mensual de los años cercanos a puede ser comparada con una simulación utilizandola producción de los datos de cobertura del suelo los datos de 1998. Para prevenir la interferencia deusados en el estudio se muestran en la Figura 9. El las variaciones de la lluvia en el efecto del cambio demes con la precipitación más alta ocurrieron en la cobertura vegetal, la precipitación promedio diariaSeptiembre para el año 1974, Julio para 1975, Agosto fue utilizada como entrada para el modelo para laspara 1998 y Octubre para 1999. La precipitación corridas, 3, 4 y 5.media anual es de 463 mm para el año 1974, 683 mm 20.
VOL. 28 Nº 1 REVISTA19 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPER Figura 9. Datos de precipitación mensual colectados en las estaciones climatológicas en el área de estudio para los años A) 1974, B) 1975, C) 1998 y D) 1999Calibración y verificación del modelo SWAT para simulados se muestran en la figura 10A. La diferenciael escurrimiento superficial porcentual se encuentra entre -3.9% y 39.3% duranteEl flujo base y el escurrimiento superficial estimado los meses de lluvia (Septiembre, Octubre,utilizando el modelo SWAT fue de 39.5% y 60.5% del Noviembre) y en un rango de -148% a 11.8% durantetotal del flujo del río respectivamente. La Figura 10A el resto del año. La diferencia para el acumuladomuestra los resultados de la calibración para el anual es de 1.5%. Los resultados de la verificación delescurrimiento superficial para el año 1974. Los modelo se muestran en la Figura 10B y 10Cparámetros de entrada que fueron usados en la respectivamente. La diferencia porcentual para elcalibración final del modelo fueron available soil water acumulado anual es de 2.1% para el año 1986 y 6.8%capacity = 0.28, soil evaporation compensation factor para el año 1994. La diferencia porcentual para los= 0.65, treshold depth of water in shallow aquifer for acumulados mensuales es 27.8% para el año 1986 yevaporation to occur = 1020 mm y baseflow recession 41.7% para el año 1994 durante los meses de lluvias.constant = 0.0313. Las diferencias porcentuales entrelos escurrimientos superficiales medidos y los 21.
VOL. 28 Nº 1 REVISTA20 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPER Figura 10. A) Resultados de la calibración del modelo SWAT para los escurrimientos superficiales mensuales para 1974, B) resultados de la verificación del modelo para el escurrimiento mensual superficial para 1986, C) resultados de la verificación del modelo para el escurrimiento mensual superficial para 1994. Las barras azules y los valores indican la diferencia porcentual entre los escurrimientos simulados y medidos.Escurrimiento simulado promedio anual incremento en el escurrimiento superficial entre elLa figura 11A muestra el escurrimiento simulado año 1974 y el año 2020 se encuentra en laspromedio anual para toda el área de estudio. Los subcuencas Áreas Urbanas y Cañón del Huajuco.cambios en la cobertura del suelo inducen un Hay que recordar que las superficies impermeablesincremento del 0.7% y del 4.2% en el escurrimiento en la subcuenca Áreas Urbanas aumentaron 18.7%superficial y un decremento del 5.0% y 6.7% en el flujo de 1974 a 2020. Sin embargo este cambio en labase para el año 1998 y el año 2020 comparado a subcuenca Áreas Urbanas causa un incremento del1974. Sin embargo, el escurrimiento total ( la suma 63.5% en el escurrimiento superficial para el añodel escurrimiento y el flujo base) a la salida de la 2020 en esta subcuenca. Para el Cañón del Huajuco,cuenca disminuye 2.4% y 1.8% en los años 1998 y un incremento del 197% en áreas impermeables2020 respectivamente a pesar de los cambios en la causa un incremento del escurrimiento superficial delcobertura vegetal en las subcuencas de Áreas 17.8% para el año 2020.Urbanas y Cañón del Huajuco respectivamente. Lafigura 11B muestra los cambios en el escurrimientosimulado para las seis subcuencas. El mayor 22.
VOL. 28 Nº 1 REVISTA21 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPERUna comparación de los cambios en la infiltración anual. Las mayores diferencias en el flujo base entreentre las seis subcuencas es mostrada en la Figura el año 1974 y el 2020 se encuentran en las11C. Las mayores caídas se encuentran en el Cañón subcuencas Áreas Urbanas (-54.7%).Esto esdel Huajuco (10.4%) y Áreas Urbanas (66.4%) para el causado por los pastizales y el matorral desértico queaño 2020. En el caso de la subcuenca Bosque, la se espera sean convertidos en superficiesinfiltración que aumenta 5.2% produce una reducción impermeables en la subcuenca Áreas Urbanas parael escurrimiento superficial del 15.4% para el año el 2020 tal y como se mostró en la Figura 8.2020. La figura 11D muestra el flujo base simulado Figura 11. A) Escurrimiento simulado anual y cambios porcentuales en los años 1998 y 2020 comparados con el año 1974 para toda el área de estudio, B) escurrimiento superficial simulado promedio anual y cambios porcentuales para los años 1998 y 2020 comparados con el año 1974 para las seis subcuencas, C) infiltración simulada promedio anual y cambios porcentuales para los años 1998 y 2020 comparados con el año 1974 para las seis subcuencas, D) flujo base simulado promedio anual y cambios porcentuales para los años 1998 y 2020 comparados con el año 1974 para las seis subcuencas 23.
VOL. 28 Nº 1 REVISTA22 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPERDISCUSIÓN Y CONCLUSIONES para el año 2020. Por lo tanto un incremento en el escurrimiento superficial muestra la posibilidad delLos datos del tipo de suelo y cobertura de suelo incremento en la ocurrencia de inundaciones en ladisponibles en México fueron adaptados para que estación de lluvias. Los resultados muestran unafuncionaran como entradas para el modelo SWAT. tendencia hacia la reducción del flujo base en elLos datos de lluvia fueron adecuados a pesar de que Cañón del Huajuco y las Áreas Urbanas. Estefaltaban unos cuantos días en los registros del decremento significa una reducción en el caudal deperiodo de estudio de 30 años. Los resultados de la las corrientes superficiales perennes que puedecalibración del modelo SWAT y su verificación fueron derivar en escasez local de agua superficial. A pesarlos esperados. Esto permitió que el modelo SWAT de todos los cambios en la cobertura del suelo elfuera utilizado para evaluar el efecto del cambio en la escurrimiento total del área de estudio permanececobertura del suelo en la hidrología del área de casi sin variación, un decremento del 2.4% para elestudio. Los resultados de la simulación muestran año 1998 y 1.8% para el año 2020. Estas diferenciasque el escurrimiento superficial se incrementó tan pequeñas se deben a que los cambios en laconsiderablemente en las subcuencas Áreas cobertura del suelo han ocurrido en dos subcuencasUrbanas y Cañón del Huajuco, como resultado de la que sólo cubren 23.5% del área total de estudio. Latransformación esperada de matorral y pastizales en cobertura de suelo sin alterar en las subcuencassuperficies impermeables para el año 2020. El Cañón de la Huasteca y Bosque sirven de “buffer” queincremento en el escurrimiento superficial en la absorben los impactos de la urbanización en otrassubcuenca Áreas Urbanas es de 63.5% cuando el subcuencas. Particularmente, los tipos de suelo A enárea de superficies impermeables aumenta 18.7%. el Cañón de la Huasteca y Bosque permiten una fácilPara la subcuenca Cañón del Huajuco un incremento infiltración del agua en el suelo independiente dedel 19.7% en las superficies impermeables causa un cambios en la cobertura vegetal del suelo.aumento del 17.8% en el escurrimiento superficial 24.
VOL. 28 Nº 1 REVISTA23 JUNIO 2009 EVALUACIÓN DEL IMPACTO DEL CAMBIO DE USO DE SUELO SELPERBIBLIOGRAFÍA · Kostadinov Stadimir, S.Vucicevic, S. Markovic, M Nikolic. 1998. Runoff regime in small watersheds with· Arnold Jeff. 2000. SWAT, Soil and Water Assessment different degrees of forest cover. Fourth InternationalTool Software. Temple, TX. USDA Agricultural Conference on Headwater Control, Merano, Italy. A. A.Research Service. Balkema Editor.· Comisión Nacional del Agua. 2005 Síntesis de · Levick, L. R., D. J. Semmens, D. P. Guertin, I. S.Estadísticas del Agua en México Mexico D.F. Burns. S. N. Scott. C. L. Unkrich, D. C. Goodrich. 2002. Adding global· Cotler Helena, Bocco Gerardo. 2004. Manejo soils data to the Automated Geospatial WatershedIntegral de Cuencas en Mexico. Mexico D.F., Instituto Assessment Tool. Temple, AZ. USDA -ARSNacional de Ecologia, pp. 11, 41. Southwest Watershed Research Center.· Descheemaeker, K., J. Nyssen, J. Poesen, D. Raes, · Miller, S.N., D.J. Semmens, R.C. Miller, M.M. Haile, B. Muys, S. Deckers. 2006. Runoff on slopes Hernandez, D.C. Goodrich, W.P. Miller. W.G. Kepner,with restoring vegetation: A case study from the Tigray D.W. Ebert. 2002. GIS-Based Hydrologic Modeling:highlands, Etiopia. Journal of Hydrology. Volume 331, The Automated Geospatial Water Assessment Tool,Issues 1.2, pp. 219-241. Documentation and user manual. Tucson, AZ.· FAO, United Nations Food and Agriculture · Moreno Jose Alberto. 2002. Mapa de uso de suelo yOrganization. 1974. Legend of the soil map of the vegetación para los años 1974 y 1998, escalaworld. Paris, France. UNESCO. 1,250:000, Mexico D.F. INEGI / SEMARNAT.· Guajardo, A. A. 2003. Crecimiento urbano enM o n t e r r e y, A n á l i s i s e s t r a t é g i c o d e l a r e a · Neitsch S.L., J.G. Arnold, J.R. Kiniry, R. Srinivasan,metropolitana de Monterrey. Monterrey, Mexico. J.R. Williams. 2000. Soil and Water Assessment ToolITESM. User´s Manual. Chapter 13, pp. 125 SWAT input data: .WGN. Temple, TX. Blackland Research Center.· Heil Costa, M., A. Botta, J. A. Cardillo. 2003. Effectsof large scale changes in land cover on the discharge · Singh Vijay P. 1992. Elementary hydrology.of Tocantins River, Southwestern Amazonia. Journal Englewood Cliffs, NJ. Prentice Hallof Hydrology, Volume 283, Issues1-4, pp. 206-217. · Taniguchi, M. 1997. Subsurface water responses to· IMTA, Instituto Mexicano de Tecnología del Agua. land cover/use changes, and overview. Boston, MA.1999. Extractor Rápido de Información Climatologica, Kluwer Academic Publishers.ERIC II, Cuernavaca, Mexico. IMTA. · USDA, United States Department of Agriculture.· IMTA a, Instituto Mexicano de Tecnología del Agua 1986. Technical Release 55 – Urban hydrology for1999. Banco nacional de datos de aguas superficiales, small watersheds. Washington, DC. NaturalBANDAS, Cuernavaca, Mexico. IMTA. Resources Conservation Service – Conservation Engineering Division.· INEGI , Instituto Nacional de Estadistica Geografia eInformatica. 1990. Guia para la interpretacion · Williamson D.R., R.A. Stokes, J.K. Ruprecht. 1987.cartografica, climatología Response of input and output of water and chloride to clearing for agriculture. Journal of Hydrology. Volume· INEGI , 2000. Mapa de edafología, escala 1:250,000 94. Issues 1-2. pp. 1-28 25.
VOL. 28 Nº 1 REVISTA24 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPER Análise dos fatores que controlam a distribuição espacial da turbidez no reservatório de Itupararanga – SP Letícia Sabo Boschi1; Maria de Lourdes Bueno Trindade Galo1; Evlyn Márcia Leão de Moraes Novo2 1 Programa de Pós-Graduação em Ciências Cartográficas – FCT/Unesp leticiasabo@yahoo.com.br; mlourdes@fct.unesp.br 2 Programa de Pós-Graduação em Sensoriamento Remoto – INPE - evlyn@ltid.inpe.brRESUMO the presence of suspended solids, organic and inorganic matter, microscopic organisms and algae inA turbidez mede a atenuação da luz causada pela the water column. These materials may originate frompresença de matérias sólidas em suspensão, matéria the soil – when there is no riparian vegetation –,orgânica e inorgânica, organismos microscópicos e mining, industries, domestic sewage withoutalgas na coluna d’água. Esses materiais podem treatment (UFRRJ, 2007), and the resuspension oforiginar-se do solo (quando não há mata ciliar); da bottom sediments caused by wind action. Takingmineração; das indústrias ou do esgoto doméstico these controlling factors into account, we tested thesem tratamento (UFRRJ, 2007) e da re-suspensão de hypothesis that the spatial distribution of turbidity insedimentos de fundo pela ação de ventos. Levando- the reservoir Ituparaganga-SP is determined by these em conta esses fatores controladores, testou-se a distance between a given class of turbidity and thehipótese de que a distribuição espacial da turbidez no source of the material that gives rise to it. To do so, wereservatório de Ituparaganga-SP é determinada pela draw up a map of potential soil erosion – by crossingdistância entre uma dada classe de turbidez e a fonte the maps of slope, soil types, and land cover – ado material que lhe dá origem. Para isso, turbidity map and a map of distances between classesconfeccionou-se um mapa de potencial do solo à of turbidity and classes of potential soil erosion. Theerosão, cruzando os mapas de declividade, tipos de analysis of the results showed that it is not possible tosolos e cobertura do solo; um mapa de turbidez e um infer the spatial distribution of turbidity based only onmapa de distâncias entre as classes de turbidez e as the distance between the classes of turbidity and theclasses de potencial de erosão do solo. A análise dos classes of potential erosion, for the maps show thatresultados mostrou que não é possível inferir a the classes of higher turbidity do not come close to thedistribuição espacial da turbidez com base somente classes of high soil erosion potential, but to thena distância das classes de turbidez às classes de classes of forest of the land cover map despite havingpotencial à erosão, pois os mapas mostram que as low to medium erosion potential. This means thatclasses de maior turbidez não se apresentaram other hydrodynamic factors – not just the proximity –próximas às classes de alto potencial do solo à must control the turbidity, such as reservoir depth,erosão, mas sim, próximas às classes de floresta do size and shape of the water body in relation to themapa de cobertura do solo, mesmo possuindo dominant directions of wind, among other factors.potencial de erosão de baixo a médio. Isso significa While erosion is a process that occurs over time inque outros fatores hidrodinâmicos e não apenas a certain seasons, the turbidity may be a temporaryproximidade devem controlar a turbidez, tais como phenomenon resulting from a pulse of sediment orprofundidade do reservatório, dimensão e forma do from the sporadic action of the wind. These aspectscorpo d’água em relação às direções dominantes de should be analyzed in future work.vento, entre outros fatores. Keywords: Factors controlling the spatial distributionEnquanto, a erosão é um processo que ocorre ao of turbidity, potential soil erosion, geoprocessing.longo do tempo em determinadas estações do ano, aturbidez pode ser um fenômeno passageiro INTRODUÇÃOdecorrente de um pulso de sedimentos ou da açãoesporádica do vento. Esses aspectos deverão ser O barramento de um rio, com a finalidade de geraranalisados em trabalhos futuros. energia elétrica, irrigação, abastecimento d’ água ePalavras-chave: Fatores controladores da outros fins, é denominado reservatório. Contudo, nãodistribuição espacial da turbidez, potencial de erosão somente aspectos benéficos devem serdo solo, geoprocessamento. considerados, uma vez que a construção de barragens afeta o ambiente previamente existente,ABSTRACT alterando o funcionamento natural dos sistemas hidrológico, biológico, geológico e climático (VIEIRATurbidity measures the attenuation of light caused by et. al, 2002). 26.
VOL. 28 Nº 1 REVISTA25 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPERA qualidade ambiental dos reservatórios está a distribuição espacial da turbidez no reservatório éintrinsecamente associada às condições de determinada pela distância das classes de turbidezqualidade de água de seus rios contribuintes, à às classes de tendência de carreamento deocupação e uso do solo no entorno do reservatório, sedimentos. Para testar essa hipótese foramaos tipos de solos, à declividade, às atividades formulados os seguintes objetivos específicos:econômicas implantadas e às condições sanitárias produzir um mapa de potencial do solo à erosão, ado entorno e do tempo de residência de suas águas. partir do cruzamento dos mapas de declividade, tiposDe acordo com Mota (1995), os programas mais de solos e cobertura do solo; produzir um mapa derecentes de gestão dos recursos hídricos têm dado classes de turbidez a partir de medidas “in situ”maior destaque às medidas preventivas de controle, coletadas em 200 estações amostrais do reservatóriopois são mais eficientes e menos onerosas, evitando- e produzir um mapa de distâncias entre as classes dese prejuízos econômicos e sociais devido à poluição. turbidez e as classes de potencial do solo à erosão.Logo, o planejamento adequado do uso e ocupaçãodo solo tem sido uma melhor forma de prevenir a TURBIDEZ COMO INDICADOR DA QUALIDADEpoluição ambiental. Nesse contexto, o planejamento DA ÁGUAe gestão dos recursos hídricos dependemfundamentalmente de informações confiáveis que As águas pluviais ao escoarem pelo solo podemsão indispensáveis para um adequado carregar impurezas para os corpos d`água que vãogerenciamento dentro de bases sustentáveis. depender de vários fatores: uso do solo na área,O Brasil é um país com poucos sistemas de fatores hidrológicos, cobertura do solo, entre outros.monitoramento de qualidade da água e que, dispõe O processo erosivo no solo pode resultar node uma quantidade reduzida de informações sobre o transporte de partículas para os corpos d’água,estado de seus corpos da água, principalmente em contribuindo para modificações na qualidade, erelação às suas dimensões territoriais, diferenças principalmente a elevação da turbidez (MOTA, 1995).geográficas regionais e magnitude dos problemas de A turbidez representa o grau de interferência com apoluição. Este fato, mostra a necessidade do país passagem da luz através da água, conferindo umadesenvolver sua tecnologia de redes de aparência turva à mesma, causada pela presença demonitoramento de qualidade da água e reconhecer a matérias sólidas em suspensão (silte, argila, sílica,importância de investir no melhoramento e ampliação colóides), matéria orgânica e inorgânica finamentede seus bancos de informações no setor de recursos divididas, organismos microscópicos e algas. Ahídricos (REBOUÇAS, 2002). origem desses materiais pode ser o solo (quando nãoAs atividades de monitoramento e avaliação da há mata ciliar); a mineração (como a retirada de areiaqualidade de superfícies aquáticas podem ser ou a exploração de argila); as indústrias; ou o esgotorealizadas utilizando-se tecnologias ou técnicas, doméstico, lançado no manancial sem tratamentocomo o uso de sensoriamento remoto e análise de (UFRRJ, 2007).imagens de satélites. Se por um lado, existe uma Há distinção entre a matéria suspensa ou sedimento,limitação técnica dos dados multiespectrais que precipita rapidamente, daquela que se precipitaadquiridos sistematicamente, a partir de plataformas vagarosamente e que provoca a turbidez. A turbidez éorbitais para aplicações com água, a dimensão do encontrada em quase todas as águas de superfície,espelho d’água de muitos reservatórios brasileiros mas está normalmente ausente nas águaspode inviabilizar programas de monitoramento da subterrâneas.qualidade da água, face ao alto custo da As águas de lagos, lagoas, açudes e represasimplementação de um esquema de amostragem apresentam, em geral, baixa turbidez, porém variáveladequado. Devido ao seu caráter sinóptico e em função dos ventos e das ondas que, nas rasas,repetitivo, as imagens de satélite podem se constituir podem revolver os sedimentos do fundo. Via de regra,em um recurso auxiliar valioso no processo de alocar após uma chuva forte, as águas dos mananciais deadequadamente pontos de amostragem no superfície ficam turvas, graças ao carreamento dosreservatório e mesmo de mapear a distribuição sedimentos das margens pela enxurrada. Assim, osespacial de alguns componentes presentes no corpo solos argilosos e as águas em movimentação,d’água, correlacionados com sua variabilidade ocasionam turbidez.espectral, como é o caso da turbidez. A turbidez, além de reduzir a penetração da luz solarO aumento de turbidez nos períodos de maior na coluna d´água, prejudicando a fotossíntese dasprecipitação pode ocorrer em função do grande algas e plantas aquáticas submersas, pode recobriraporte de material que é carreado pelas chuvas para os ovos dos peixes e os invertebrados bênticos. Oso corpo d’água em questão, e pela mistura na água sedimentos em suspensão podem carrear nutrientescausada pelo vento, ocasionando uma re-suspensão e pesticidas, obstruindo as guelras dos peixes e atéde nutrientes dos sedimentos, em regiões menos interferir na habilidade do peixe em se alimentar e seprofundas (ROCHA, 2003). defender dos seus predadores. As partículas emNesse contexto, o objetivo geral desse trabalho é suspensão localizadas próximo à superfície podemanalisar os fatores que influenciam a istribuição absorver calor adicional da luz solar, aumentando aespacial da turbidez da água do reservatório de temperatura da camada superficial da água (UFRRJ,Itupararanga e a hipótese a ser investigada é a de que 2007). 27.
VOL. 28 Nº 1 REVISTA26 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPEROs padrões de qualidade da água referem-se a um Itupararanga, produzido por Queiroz (2006) ecerto número de parâmetros capazes de refletir a realizado a partir de imagens obtidas pelo sensorpresença efetiva ou em potencial de algumas CCD do satélite Cbers (bandas 2, 3 e 4),substâncias ou microorganismos que possam disponibilizadas gratuitamente pelo INPE – Institutocomprometer a sua qualidade estética ou sanitária de Pesquisas Espaciais (INPE, 2006), referente à(REBOUÇAS, 2002). coordenada órbita/ponto 155/126, com resoluçãoPara uma melhor representação das classes de espacial de 20 metros, obtidas em 21 de janeiro deturbidez é melhor determiná-las baseadas em algum 2006;padrão estabelecido. Por isso, de acordo com a • Mapa esquemático de solos produzido pelo ComitêResolução no 20, de 18 de junho de 1986 do da Bacia Hidrográfica do Rio Sorocaba e Médio TietêCONAMA (Conselho Nacional do Meio Ambiente) o – CBH, na escala 1:250000, em formato .pdf;Reservatório de Itupararanga pertence à classe 2 • Curvas de nível vetorizadas de cartas topográficas(BRASIL, 2006). O CONAMA estabelece ainda que a na escala 1:50000 de Sorocaba e São Roque, queturbidez permitida na classificação das águas vai até abrangem a bacia do reservatório de Itupararanga;100 UNT (Unidades Nefelométricas de Turbidez). • Dados de Turbidez, coletados em 21 de julho de 2006, pelos pesquisadores da FCA (Faculdade deMETODOLOGIA Ciências Agrárias)/ Unesp de Botucatu, engenheiros agrônomos Drs. Eduardo Negrisoli e Marcelo RochaA área de estudos do presente trabalho refere-se à Corrêa. Nesse levantamento no reservatório deárea de influência do reservatório de Itupararanga, Itupararanga, foram coletados 200 elementosEstado de São Paulo, cuja delimitação foi feita com amostrais posicionados com GPS de navegação.base na identificação das áreas que drenam em Estes pontos de coleta foram distribuídos por todo odireção ao reservatório (bacia hidrográfica da represa reservatório e caracterizados quanto à Turbidez dade Itupararanga). Desta maneira, não foi considerada água.toda bacia do rio Sorocaba, principal formador doreservatório de Itupararanga (QUEIROZ, 2006). Definição das classes de turbidez A Figura 1 a seguir, apresenta a distribuição dos 200Dados de entrada pontos amostrais referentes à turbidez, coletados por• Mapa de cobertura do solo do reservatório de todo o reservatório, em 21 de julho de 2006. Figura 1 – Distribuição dos pontos amostrais ao longo do reservatório de ItupararangaA partir dos pontos de turbidez, foi gerada uma grade Os intervalos utilizados para definir as classes detriangular (TIN) e na seqüência, uma grade retangular. turbidez foram definidos por Saito (2007) e sãoA partir dessa foi realizado o fatiamento, produzindo apresentados na Tabela 1 a seguir.assim, o mapa de classes de turbidez (Figura 2), oqual foi recortado a partir do limite do reservatório. Tabela 1 – Classes de turbidez usando intervalos definidos por Saito (2007) 28.
VOL. 28 Nº 1 REVISTA27 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPER Figura 2 – Mapa de classes de turbidez usando intervalos definidos por Saito (2007)Justifica-se o uso dos intervalos definidos por Saito desses não terem sido representativos, pois o(2007) ao invés dos intervalos baseados no padrão reservatório de Itupararanga não apresenta valoresestabelecido pelo CONAMA (Tabela 2), pelo fato de alta turbidez (acima de 100 FTU). Tabela 2 – Classes de Turbidez usando intervalos do CONAMAGeração do mapa de Declividade seqüência para gerar a declividade. Realizando oA partir da importação das curvas de nível foi possível fatiamento da declividade produziu-se o mapa degerar a grade triangular (TIN), utilizando como linhas classes de declividade (Figura 3), o qual foi recortadode quebra os afluentes do reservatório e em seguida a partir do limite da bacia do reservatório degerar a grade retangular, a qual foi utilizada na Itupararanga. 29.
VOL. 28 Nº 1 REVISTA28 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPER Figura 3 – Mapa de classes de declividade da bacia do reservatório de ItupararangaAnalisando a Figura 3 verifica-se que a bacia do Geração do mapa de solosreservatório de Itupararanga não apresenta muitas Vetorizando o mapa esquemático de solos (escalaáreas com alta declividade. Portanto, pode-se inferir 1:250000) gerou-se uma mapa de tipos de solosque a área de estudos é uma região predo- referente somente à bacia do reservatório deminantemente plana. Itupararanga (Figura 4). Figura 4 – Mapa de tipos de solos da bacia do reservatório de Itupararanga 30.
VOL. 28 Nº 1 REVISTA29 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPERPelo fato do mapa de tipos de solo da bacia do de potencial do solo à erosão foi realizada umareservatório de Itupararanga ser oriundo de um mapa operação booleana utilizando a linguagem LEGALproduzido na escala 1:25000, foram encontrados (Software Spring), na qual foram cruzados os mapassomente dois tipos de solo na bacia do reservatório produzidos (mapa de tipos de solos da bacia dode Itupararanga, são eles: predomínio de Argissolos reservatório e mapa de declividade da bacia doVermelho-Amarelos Distróficos (solo não ácido), reservatório) e o mapa de cobertura do solo (Figura 5),texturas média/argilosa e argilosa, A moderado (PVA- produzido por Queiroz (2006).1) e Predomínio de Latossolos Vermelho-Amarelos Gerando assim, o mapa de potencial do solo à erosãoDistróficos, textura argilosa A moderado (LVA-5). (Figura 6).Mapa de potencial do solo à erosãoA fim de produzir um mapa que indicasse as classes Figura 5 – Mapa de cobertura do solo da bacia do reservatório de Itupararanga 31.
VOL. 28 Nº 1 REVISTA30 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPER Figura 6 – Mapa de potencial do solo do reservatório de Itupararanga à erosãoComparando o mapa de potencial do solo à erosão predominantemente dos tipos:(Figura 6) com o mapa de classes de declividade floresta, pastagem e culturas. O que faz com que não(Figura 3) é possível verificar que no mapa de haja grandes variações, inferindo a essa variável umapotencial à erosão, a variável que teve maior participação não muito significativa. A Figura 7 aparticipação foi a declividade. Isso é justificável pelo seguir apresenta uma comparação entre os doisfato do mapa de solos apresentar somente dois tipos mapas.de solo e pelo fato da cobertura do solo ser Figura 7 – Comparação entre o mapa de classes de declividade (a) e o mapa de potencial do solo à erosão (b) 32.
VOL. 28 Nº 1 REVISTA31 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPERMapa de distâncias foram selecionadas as classes individualmente ePara testar a hipótese de que a distribuição espacial então gerados os mapas de distância entre a classeda turbidez no reservatório é resultado do aporte de de alto potencial à erosão, às classes de turbidez esedimentos e de que essa distribuição é determinada assim sucessivamente.pela distância das classes de turbidez às classes de Como resultado desses mapas de distâncias,potencial de transporte de sedimentos, foram obtiveram-se grades numéricas, as quais em seguidagerados os mapas de distâncias referentes a cada foram fatiadas nas classes denominadas: próxima,uma das classes de potencial à erosão. intermediária e distante das áreas de potencial doUm mapa de distância é um tipo de análise de solo à erosão. Os intervalos referentes a cada umaproximidade que apresenta zonas com larguras das classes definidas é apresentado na Tabela 3 aespecificadas (distâncias) em torno de um ou mais seguir e os mapas temáticos representando oselementos do mapa. Para gerar tal mapa, utilizou-se o mapas de distância são mostrados nas Figuras 8, 9 emapa temático de potencial do solo à erosão, do qual 10. Tabela 3 – Intervalos das classes de distância às classes de turbidez Figura 8 – Mapa de distância entre as classes de turbidez e a classe de alto potencial à erosão 33.
VOL. 28 Nº 1 REVISTA32 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPER Figura 9– Mapa de distância entre as classes de turbidez e a classe de médio potencial à erosão Figura 10 – Mapa de distância entre as classes de turbidez e a classe de baixo potencial à erosãoAnalisando as Figuras 8, 9 e 10 nota-se que a classe distribuição das classes de turbidez não depende tãodenominada próxima (verde escuro) às áreas de alto somente da distância a que essas se apresentam daspotencial do solo à erosão são quase inexistentes, classes (alto, médio e baixo) de potencial do solo àpara os três cruzamentos com o mapa de classes de erosão. Uma ilustração clara desta situação pode serturbidez. Isso se justifica pelo fato das classes serem visualizada na Figura 11 (a) a seguir, onde as classesgerais e pelo fato do fenômeno ter uma escala de de turbidez definidas como maior não se apresentamvariabilidade mais fina, sendo assim, os intervalos de próximas às áreas de alto potencial à erosão. A Figuraclasse identificados acabam não tendo sensibilidade. 11 (b) confirma este fato, pois mostra que a classeA partir dos mapas, pode-se inferir então, que a denominada próxima não existe nesta mesma área. 34.
VOL. 28 Nº 1 REVISTA33 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPER Figura 11 – (a) mapa de potencial à erosão associado ao mapa de classes de turbidez e (b) mapa de potencial à erosão associado ao mapa de distância da classe denominada alto potencial às classes de turbidezA Figura 12 (a) a seguir mostra que as classes de floresta e próximas às áreas de baixo e médiomaior turbidez apresentam-se próximas às áreas de potencial do solo à erosão. Figura 12 – (a) associação das classes de turbidez com a cobertura do solo (floresta) e (b) associação das classes de turbidez com o mapa de potencial à erosãoCONCLUSÕES mas também ao fato de que na área estudada, a turbidez não é um problema, assim como tambémLevando em consideração as variáveis abordadas, não existem áreas críticas de erosão.as operações realizadas, os mapas produzidos e o No presente trabalho, a declividade e a cobertura dofato de que reservatórios sofrem influência de quatro solo mostraram-se bastante responsáveis peladimensões (lateral, longitudinal, vertical e distribuição espacial da turbidez, como mostra atemperatura) pode-se concluir que a hipótese Figura 12. A Figura 12 (a) mostra que classes deestudada, na qual foi analisada se a distribuição maior turbidez apresentam-se próximas à classe deespacial da turbidez no reservatório pode ser floresta do mapa de cobertura do solo, mesmo essasdeterminada pela distância das classes de turbidez áreas apresentando baixo e médio potencial do solo ààs classes de potencial de transporte de sedimentos, erosão, o que significa que a floresta é responsávelnão é verdadeira, como mostra a Figura 11. pelo aporte de material orgânico, que provoca aNo desenvolvimento do trabalho verificou-se que a turbidez. Analisando a Figura 12 (b) verifica-se que osdistribuição espacial da turbidez também é resultado afluentes apresentam-se também comoda atuação de outras variáveis, como declividade, responsáveis pelo carreamento de material.tipos de solo, cobertura do solo, morfologia da área de Como recomendação para trabalhos futuros sugere-estudo (profundidade e atuação do vento), porém, o se que se possua uma maior quantidade de variáveis,fato da relação entre as classes de turbidez e as para que se possa inferir de forma mais adequadaclasses de potencial do solo à erosão ter sido tênue, quais os fatores que mais contribuem para anão se deve tão somente aos fatores hidrodinâmicos, distribuição espacial da turbidez. 35.
VOL. 28 Nº 1 REVISTA34 JUNIO 2009 ANÁLISE DOS FATORES QUE CONTROLAM A DISTRIBUIÇÃO ... SELPERBIBLIOGRAFIAS MOTA, S. Preservação e Conservação de RecursosBarbosa, C. C. F. Álgebra de mapas e suas Hídricos. 2oEdição. Rio de Janeiro: ABES, 1995.aplicações em sensoriamento remoto egeoprocessamento. Dissertação de Mestrado em MOTA, S. Introdução à engenharia ambiental. Rio deSensoriamento Remoto. INPE, São José dos Janeiro: ABES, 2003. 419p.Campos 1997. Disponível em: <http://www.dpi.inpe.br/teses/claudio/CAP1.html.>. Acesso em: 2007. REBOUÇAS, A. C. BRAGA, B.; REBOUÇAS, A. C.; TUNDISI, J. G. (Org.). Águas doces no Brasil: capitalBRASIL. Ministério do Meio Ambiente-Conselho ecológico, uso e conservação. 2º Edição. São Paulo:Nacional do Meio Ambiente. Resolução no20 de 18 Escrituras Editora, 2002.de junho de1986, Resoluções do CONAMA, Brasília, DF, 1aEd, ROCHA, R.R.A. & THOMAZ, S. M. Variação temporal2006. Disponível em : de fatores limnológicos em ambientes da planície de <http://www.mma.gov.br/port/conama/processos/61 inundação do alto rio Paraná (PR/MS – Brasil). 2004.AA3835/LivroConama.pdf>. Acesso em: fev. 2007. Acta Scientiarum, p. 261- 271.CAMARA, G. et al. Álgebra de Mapas. Cap.8. Notas SOS ITUPARARANGA. Caracterização Geoam-de aula. 2007. biental da Bacia da Represa de Itupararanga, Bacia Hidrográfica do Rio Sorocaba – SP. Disponível em:Câmara, G. Modelos, linguagens e arquiteturas para <http://www.sositupararanga.com.br/biblioteca/index.banco de dados geográficos. São José dos Campos. asp>. Acesso em: fev. 2007.Tese (Doutorado em Computação Aplicada) - InstitutoNacional de Pesquisas Espaciais - INPE, 1995. SAITO, E. A..Mapeamento dos níveis de turbidez na água a partir da classificação de imagensCARACTERIZAÇÃO geral da UGRH 10, 2000. multiespectrais e multisensor, utilizando redesDisponível em: neurais artificiais. Relatório de atividades parciais<http://www.rededasaguas.org.br/comite/relsmtseg. PIBIC 2006/2007.pdf>. Acesso em: nov. 2006. TUNDISI, J. G. & STRASKRABA, M., GerenciamentoESTEVES, F.A. Fundamentos de Limnologia. 1998. da Qualidade da Água de Represas. 2000. RIMA.Interciência/FINEP, Rio de Janeiro: 575 p. 251p.FELGUEIRAS, C. A. e CÂMARA, G.. Modelo UFRRJ – Universidade Federal Rural do Rio deNumérico do Terreno. Cap.7. Notas de aula. 2007. Janeiro. Pesquisa em Turbidez. Disponível em: <http://www.ufrrj.br/institutos/it/de/acidentes/turb.htQUEIROZ, R. P. Mapeamento da cobertura do solo m>. Acesso em: 2007.na bacia do reservatório de Itupararanga – SP. 78p.Trabalho de Graduação – Universidade Estadual VIEIRA, M. S.; FERREIRA, J. R.; CASTRO, P. M. G.;Paulista, SP. 2006. ROCHA, A. A. Aspectos da química da água e do sedimento do reservatório de Ibitinga (São Paulo –INPE – Instituto Nacional de Pesquisas Espaciais. Brasil – 21° 45’ S e 48° 50’ W). São Paulo: Boletim doCatálogo de Imagens CBERS. Disponível em: Instituto de Pesca, 2002. 15p.<www.dpi.inpe.br/catalogo>. Acesso em: 2006. 36.
VOL. 28 Nº 1 REVISTA35 JUNIO 2009 ADECUACIÓN DEL DEM DEL SRTM DE 90 METROS ... SELPER Adecuación del DEM del SRTM de 90 metros para usos hidrológicos en llanuras Sergio Masuelli;1 Pablo M. Vazquez2 y Omar Faure3 1. Instituto Mario Gulich, CONAE-UNC, Ruta C45 km 8 Córdoba (5187) Argentina, Prof. Asociado, sergio.masuelli@conae.gov.ar 2. INTA Rauch, EEA Cuenca del Salado. Av. Belgrano 416 (7203) Rauch, Buenos Aires Argentina. Investigador, pmvazquez@correo.inta.gov.ar 3. UTN Fac. Reg. Concepción del Uruguay, Ing. Pereyra 676 C. P. E3264BTD Concepción del Uruguay Entre Ríos República Argentina, Prof. Asociado, ofaure@frcu.utn.edu.ar Palabras claves: DEM, SRTM, hidrología de llanuras, filtrosRESUMEN centímetros. En tal sentido, en el presente trabajo mostramos elEl tener a disposición el DEM del SRTM de 90m ha DEM resultante luego de diversas mejoras para sucambiado significativamente la calidad de la utilización como base topográfica para un modeloinformación topográfica disponible, dado que es una hidrológico en llanura. Estas incluyen el filtrado delmejora sustantiva respecto del GTOPO 30 (900m de monte (usando imágenes Landsat), el filtrado delresolución espacial) y de las cartas del IGM en ruido de baja frecuencia y amplitud mencionado,muchos casos también, puesto que su alta resolución ruido de alta frecuencia en general y correccioneslo hace sensible a cambios en la topografía debido a para cuerpos de agua.obras humanas.Sin embargo, el DEM del SRTM que se baja por ABSTRACTinternet desde la página de la Universidad deMaryland tiene esencialmente la información en un The accessibility to the SRTM 90m DEM has greatlyestado crudo por cuanto la altura que muestra para changed the quality of the available topographiccada píxel tiene que ver con la primera superficie information. It is a significant improvement if it isreflectante, es decir por ejemplo que está dando compared to the GTOPO 30 (900m spatial resolution)cuenta del tope de la foresta y no del suelo en si and also, to the IGM charts, since its high resolutionmismo. Además es reconocido por el USGS que makes it sensitive to the any modification in theestos datos tienen algunas fallas, como ser datos topography of the observed areas as a result of man-extremadamente altos o extremadamente bajos, made changes.especialmente en zonas que involucran cuerpos de However, the SRTM DEM of that can be downloadedagua pues su nivel de backscatter es muy bajo y por lo from the webpage of the University of Marylandtanto su procesamiento es generalmente muy ruidoso. basically provides information wich has not beenPor ello se recomienda al usuario que se interiorice completely processed, as the height per pixel it showscon las características del DEM (resolución, precisión, has to do with the reflecting surface that appears first,generación de datos espurios, etc.) para juzgar la that is, it accounts for the top of the forest and not ofutilidad de esta información para la aplicación the ground itself. In addition, USGS has recognizedespecífica que se desea. Es claro que los datos there are some faults in these data: they may beespurios que tiene pueden determinar la total extremely high or low, especially in areas where thereimposibilidad para un uso específico mientras que are masses of water involved, as their backscatterpara otra aplicación resultar éstos irrelevantes. level is too low and consequently, their processing isLa resolución vertical (relativa) declarada es mejor generally very noisy. For this reason, users areque 10m, aunque algunos estudios la sitúan en recommended to be thoroughly familiar with thealrededor de 2m para la versión de 30m, y los datos characteristics of the DEM (resolution, accuracy,están digitalizados al metro. Por encima de esta generation of spurious data, etc) to be able to evaluateresolución observamos que para zonas muy planas the usefulness of the information for any specificestá presente un ruido de baja frecuencia de 1 metro application. It is clear that the spurious data mayde amplitud que es notable en zonas de llanura en determine the impossibility to use it for certaindirección NO-SE aproximadamente. applications while they are irrelevant for others.Por todo lo antes mencionado, el uso directo de este The vertical resolution (relative) reported is betterDEM para modelos hidrológicos en zonas de llanuras than 10 m, although some studies set it around 2 m forpuede ser muy problemático. Por cuestiones de the 30m version and data are digitalized to the meter.estabilidad numérica para pendientes inferiores a Beyond that resolution, we observed a low frequency1:1000 y para celdas de 100m se necesitan noise of 1 m of amplitude in very flat areas. This isprecisiones verticales del orden de algunos remarkable in the plains located in a NW - SE direction. 37.
VOL. 28 Nº 1 REVISTA36 JUNIO 2009 ADECUACIÓN DEL DEM DEL SRTM DE 90 METROS ... SELPERTo sum up, the direct use of this DEM for hydrologic importante cualidad puede, sin embargo, ser la causamodels in flat areas may be really complex. For de un gran número de inconvenientes al tratar de usarnumeric stability reasons in slopes below 1:1000 and esta información directamente. La altura que muestrafor 100 m cells, vertical precisions in the order of some para cada píxel corresponde a la primera superficiecentimeters are required. In this sense, in the present reflectante, es decir por ejemplo que está dandowork we show the resulting DEM after including many cuenta del tope de la foresta y no del suelo en siimprovements for its use as topographical basis for a mismo. Para nuestra aplicación en hidrología, dentrohydrologic model in plains. These consist of the forest del ejemplo anterior, puede darse la paradoja de quefiltering (using Landsat images), the filtering of the low los cursos naturales de drenaje, al poseerfrequency and amplitude noise and the high generalmente monte en sus inmediaciones, tenganfrequency noise in general and corrections for water en el DEM una altura mayor que su entorno y por lomasses. tanto, para el DEM no procesado, significar un obstáculo en vez de una vía de escurrimiento.1. INTRODUCCIÓN Por esta y otras causas es que, en la documentación complementaria al DEM, el USGS (ServicioLa posibilidad de disponer libremente con el Modelo Geológico de los EEUU) recomienda al usuario quede Elevación Digital (DEM) de 90m , obtenido a partir se interiorice con las características del DEMde la Misión del Shuttle de Topografía por medio de (resolución, precisión, generación de datos espurios,Radar (SRTM), a través de la página web de la etc.) para juzgar la utilidad de esta información para laUniversidad de Maryland, ha cambiado aplicación específica que se deseasignificativamente la calidad de la información (http://hydrosheds.cr.usgs.gov/datasets). Es clarotopográfica disponible para diversos usos. Además, que los datos espurios que tiene pueden determinarpermite una excelente base de trabajo para el uso la total imposibilidad para un uso específico mientrascombinado de sensores, pues se halla agrupada en el que para otra aplicación resultar éstos irrelevantes.mismo sistema de grilla que las imágenes Landsat También se advierte que estos datos tienen algunasortorectificadas del mismo sitio web, y con un muy alto fallas, como ser datos extremadamente altos ogrado de corregistro entre ellas. En principio, sus extremadamente bajos, especialmente en zonas quecaracterísticas permiten el uso de esta información involucran cuerpos de agua pues su nivel detopográfica sin necesidad de contar con especialistas backscatter es muy bajo para las imágenes SLCen interferometría SAR (INSAR). (Rodríguez et al. (Single Look Complex) originales y por lo tanto el2005) DEM obtenido es generalmente muy ruidoso (Lehner,La resolución vertical (relativa) declarada es mejor et al., 2006)que 10m, aunque algunos estudios la sitúan en El objetivo del presente trabajo es el uso de estaalrededor de 2m para la versión de 30m (Lee et al. información como base topográfica para modelos2005), y los datos están digitalizados al metro. Para hidrológicos de la subcuenca B4 del río Salado en lazonas de llanura los errores mencionados son provincia de Buenos Aires.sensiblemente más bajos pero siempre del orden delmetro. 2. METODOLOGÍASu relativamente elevada resolución espacialcomparada con su antecesor GTOPO 30 (DEM de 2.1 Área de estudio900m de resolución espacial) y de las cartas del IGM En la Figura 1 se muestra el área de estudio.(Instituto Geográfico Militar de Argentina) permite la (Vázquez y otros, 2006)visualización de detalles topográficos como sermodificaciones antrópicas del terreno. Esta Figura 1: Ubicación geográfica de la Subcuenca B4 del río Salado 38.
VOL. 28 Nº 1 REVISTA37 JUNIO 2009 ADECUACIÓN DEL DEM DEL SRTM DE 90 METROS ... SELPER2-2. Problemas que presenta el DEM de SRTM formando un cuadriculado. El origen de estas franjasEn la figura 2a se muestra un sector, debidamente se debe a que las imágenes usadas para obtener esteecualizado pero sin ningún procesamiento, del DEM DEM son de modo ScanSAR (Thomas Farr,bajado del sitio de la Universidad de Maryland vicedirector del SRTM comunicación personal).correspondiente a la escena Landsat 225-86. En ella También es claro observar líneas largas quepodemos observar claramente franjas paralelas en corresponden a caminos y líneas geométricas másdirección NO-SE, correspondientes a un ruido de baja cortas que corresponden a cadenas de árbolesfrecuencia de 1m de amplitud aproximadamente. cortadoras de viento, típicas en la zona. Los caminosEste ruido lo hemos observado en otras imágenes, son obstáculos reales al escurrimiento superficialincluso conteniendo franjas en el sentido transversal pero no así las regiones arboladas. Figura 2: Porción del DEM donde se observa el ruido de baja amplitud. a) Original, b) Filtrado mediante FourierEn síntesis, para cumplir el objetivo de usar este DEM lagunas.como información de base para modelos hidrológicos En la figura 3 se muestra un diagrama con los pasosde llanura, necesitamos mínimamente: filtrar el ruido seguidos para obtener el DEM adecuado a nuestrade baja amplitud y frecuencia media; filtrar las aplicación.arboledas; y corregir la topografía aledaña a las Figura 3: Diagrama de flujo de los procesos para la obtención del DEM adecuado 39.
VOL. 28 Nº 1 REVISTA38 JUNIO 2009 ADECUACIÓN DEL DEM DEL SRTM DE 90 METROS ... SELPER2-3. Adecuación básica del DEM método del paralelepípedo de zonas con firmasLos píxeles con datos faltantes se completaron con espectrales similares a bosque, en las imágeneslos valores de los vecinos más próximos. Landsat 225-85 y 225-86. Luego, esta clasificaciónPara el filtrado del ruido de baja amplitud y frecuencia fue submuestreada a 90m de píxel.media, primeramente se midió sobre la imagen las Para enfatizar los rasgos salientes de la imagen sefrecuencias en vertical y en horizontal. Luego, sobre trató a la como resultante de la suma de una imagenla imagen transformada de Fourier se buscó las perturbada y otra sin perturbar. La imagen sinregiones correspondientes a la frecuencia medida y perturbar se obtuvo subdividiendo la imagen ensus armónicas y se las eliminó. Antitransformando se cuadrados de 10 píxeles de lado, calculando el valorobtiene la imagen filtrada que se muestra en la figura medio sobre cada uno y luego interpolando2b. cuadráticamente en ambas direcciones para cada punto interior del mismo. La imagen perturbada es2-4. Filtrado de la arboleda. simplemente la diferencia entre la imagen originalPara el filtrado de la arboleda, primeramente se menos la imagen sin perturbación y puedenrealizó una clasificación supervisada mediante el compararse en la figura 4. Figura 4: Porción del DEM para observar el resaltado por medio del método de perturbación. a) DEM original, b) DEM resaltadoPara clasificar los píxeles candidatos a corresponder del DEM sin perturbación y éste valor es inaccesiblea zonas arboladas se determinaron en el DEM lo que se hizo fue repetir 4 veces este proceso desderesaltado aquellos que correspondieran a árboles en la obtención de los DEM perturbados y sin perturbar.la clasificación del Landsat y con alturas mayores a En la figura 5 de muestra en secuencia cómo fue4,5m. Luego en forma iterativa se hizo crecer la actuando este proceso iterativo en la clasificación declasificación a los píxeles contiguos con alturas las arboledas para una porción de la escena.mayores a 2m. Finalmente, se eliminaron píxeles El DEM final filtrando las arboledas resulta de la restadispersos. de los valores de esta perturbación al DEM original.Como en este proceso es importante la estimaciónFigura 5: Secuencia de la clasificación de alturas de árboles. El gris oscuro corresponde a alturas sobre la media entre 2 y 4,5m; el gris claro entre 4,5 y 8m; y el blanco a mayores a 8m 40.
VOL. 28 Nº 1 REVISTA39 JUNIO 2009 ADECUACIÓN DEL DEM DEL SRTM DE 90 METROS ... SELPER2-5. Corrección para zonas bajas. los casos, en trayectos cortos de muy pocos píxeles.Como se dijo en la introducción, la cuantificación de la Por lo tanto, para el agua fluya, la discretización encapacidad de almacenamiento de las zonas bajas es altura debe ser del orden de algunos cm. Para lograruna variable crucial para el entendimiento del balance ésto, se le aplicó al DEM un filtro cuadrático de nuevehidrológico del sistema bajo estudio. Claramente, a puntos en las direcciones vertical, horizontal y lasmenos que en forma fortuita la misión SRTM hubiera diagonales principales. El efecto es que suaviza lascoincidido con una sequía histórica esta información alturas entre los píxeles vecinos pero respetando lasno está en el DEM resultante pues cada bajo alturas mínimas y máximas locales, que en un filtro decorresponde a un cuerpo de agua. media común serían severamente afectados.Como una primera aproximación para corregir estedefecto, se excavó directamente usando la 3. RESULTADOSclasificación de frecuencia de anegamiento sobre elDEM a partir de un valor de frecuencia de Como es claro que las correcciones han sidoanegamiento realizada por Pablo Vázquez (Vázquez realizadas a nivel local y como la visión total del DEMy otros, 2007) en base a una larga serie histórica de no muestra las mejoras realizadas, mostramos en laimágenes Landsat de la región (200 imágenes). El figura 6 una porción del DEM obtenido final,criterio elegido fue restar 10cm al DEM por cada 10 % comparado con los datos originales.de esta frecuencia, a partir del valor de 30 %. La verdadera bondad del trabajo realizado puede verificarse observando los resultados obtenidos por2-6 Filtrado de alta frecuencia. el modelo hidrodinámico que ha usado este DEMLuego de las correcciones anteriores persiste en el (Vázquez y otros, 2006). Como puede observarse enDEM ruido de alta frecuencia que complica su la figuras 7, el agua en superficie obtenida por elutilización para modelos hidrológicos. Téngase mismo modelo puede cambiar sustancialmentepresente que su resolución en altura (no su error) es dependiendo del DEM. Para estas regiones planas sede un metro y que para la región extremadamente observa que las “lomas” transversales del DEM crudoplana que estamos tratando de simular y teniendo en y la altura de la vegetación pueden constituirse encuenta el espaciamiento entre puntos de grilla de 90m, barreras muy importantes para el escurrimiento delel terreno se presenta muy fragmentado de modo que agua.prácticamente el agua sólo podría fluir, en el mejor de Figura 6: Porción del DEM antes (a) y después de los procesamientos (b) Figura 7: Agua en superficie obtenida por el modelo de escurrimiento luego de 10 días. Usando el DEM de la figuras 6a y 6b, respectivamente 41.
VOL. 28 Nº 1 REVISTA40 JUNIO 2009 ADECUACIÓN DEL DEM DEL SRTM DE 90 METROS ... SELPER4. CONCLUSIONES Consecuentemente, pensamos que más convenienteLos problemas señalados el DEM del SRTM sería aplicar esta corrección en forma selectiva,muestran que el mismo no puede usarse aplicándola sobre una máscara de cuerpos de aguadirectamente como base topográfica para un modelo obtenida a partir de una imagen Landsat cercana a lahidrodinámico, en especial para la región de estudio fecha de la misión SRTM (febrero de 2000).extremadamente plana, manteniendo la resolución Por otra parte, algunos de los procesamientosespacial real de 90m de píxel. Para lograr este realizados pueden ser extendidos a otrasobjetivo hay que realizar por lo menos las mejoras y aplicaciones. Por ejemplo, el método usado para elcorrecciones detalladas en este trabajo, e incluso se filtrado de la foresta puede servir para clasificar zonasdeberían hacer con mayor fineza si se espera que el arboladas, dando además una idea de su altura.modelo arroje resultados de mayor precisión.El excavado del DEM para corregir cuerpos de agua Agradecimientos.es una corrección muy importante que debe hacersecon cuidado, pues entender la dinámica de estas Los autores agradecen a las instituciones a las cualesmicrocuencas y sus interrelaciones es importante pertenecen y que posibilitaron la realización de estepara comprender el funcionamiento del sistema trabajo, INTA y FRCU-UTN. Así mismo agradecen alcompleto. La corrección que se hizo ha utilizado la Comodoro Villanueva y al capitán Boolsen del SMNfrecuencia de anegamiento aplicada en forma por la información meteorológica proporcionada, a lageneral y uniforme sobre todo el dominio, lo cual CONAE y a las profesoras Hormaiztegui y Kay, de lapuede estar realizando sobrecorrecciones en varios FRCU-UTN, por la traducción al inglés de loslugares y por lo tanto distorsionando el DEM. resúmenes.BIBLIOGRAFÍA /srtm/SRTM_D31639.pdf)Lee,I.S., H.C. Chang, and L. Ge, 2005: GPS Vazquez, P.M., S. Masuelli, G.R. Platzeck y D. Barrera,Campaigns for Validation of InSAR Derived DEMs. 2006: Modelo hidrológico con uso intensivo deJournal of Global Positioning Systems.Vol. 4, No. 1-2: información satelital para la subcuenca B4 del Río82-87. Salado Bonaerense. XII Simposio Internacional SELPER Colombia.Lehner, B., K. Verdin and A. Jarvis, 2006:HydroSHEDS Technical Documentation. World P. Vázquez, S. Masuelli, G. Platzeck y O. Boolsen,Wildlife Fund US, Washington, DC. Available at 2007: Mapa de riesgo hídrico. XII Congreso dehttp://hydrosheds.cr.usgs.gov. Teledetección de la Asociación Española de Teledetección, Mar del Plata, Argentina.Rodríguez, E., C.S. Morris, J.E. Belz, E.C. Chapin, http://hydrosheds.cr.usgs.gov/datasets.phpJ.M. Martin,W. Daffer, and S. Hensley, 2005: AnAssessment of the SRTM Topographic Products. JPL Sitio web del GLCF de la Universidad de Maryland:D-31639, Jet Propulsion Laboratory, Pasadena, http://glcfapp.umiacs.umd.edu/California, 143 pp (http://www2.jpl.nasa.gov 42.
... y Envisat sigue observando el lugar del próximo simposio (arriba izquierda) ESA - ESRIN Via Galileo Galilei, Casella Postale 64. 00044 Frascati (Roma) ITALIA Tel: (39.06) 94180.409 Fax: (39.06) 94180.280 Att: Sr. Francesco Sarti (América Latina) ESA - Asuntos Internacionales 8-10 rue Mario Nikis, 75738 PARIS Cedex 15, FRANCIA Tel: (33.1) 5369.7731, Fax: (33.1) 5369.7627 Att: Sra. Marie Elisabeth De Vel (América Latina) REVISTA - JOURNAL SELPER VOL. 28, Nº 1, 2009ECUADOR CHILE BRASIL ARGENTINA PERU COLOMBIA MEXICO VENEZUELA ARGENTINA BOLIVIA CHILE COLOMBIA1980 - 1983 1983 - 1986 1986 - 1989 1989 - 1991 1991 - 1993 1993 - 1995 1995 - 1997 1997 - 1999 2000 - 2002 2002 - 2004 2004 - 2006 2006 - 2008 CUBA2008 - 2010 Recommended
Sonny Swokowsky Braun
Impactos de cambios de uso de suelo en el area metropolitana de san salvador
Soza diego.2012.Metodologia de modelación de cuenca hidrológica
Julio Cortes Leyama
Victor Hidalgo Corpancho
PDM Colpa Bélgica
Erosion hidrica cisneros 978-987-688-024-4
06 eae 2011
anaerbol
Foro arauco
Informe final salamanca
HONDURA COURSE - Programa de incentivos economicos ambientales-ACP / Arturo Cer……
Curso gestion de cuencas peru ii
Resumen de consulta publica
Expediente tec. shulcas 31052010 corregido final separ
Erf umafor0901
Guia metodologica docentes
Experiencias de Adaptación al CC en Alta montaña Klaus Schutze Páez
Cdc 17 1 1 versión final
Paper template dinamica colombia 5th_iwfc
milanyela2
Banderillero satelital-completo
Catalogodeobjetos
Ceni space [octubre]
Gv sig 1-1-man-v3-es
Encuesta instructores ap_tabulada.3
270401040 1
270401041 1
270401039 1
270401038 1
270401037 1
270401036 1
170401011 1
170401010 1
Calibracion de combinadas
Gps tecnicana luis joel martinez

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución