Giordano, Leandro
Uso del índice Normalizado de Agua MODIS para monitoreo, modelación y pronóstico hidrológico en un gran sistema de llanura con morfología fluvial típica [TESIS]. -- Luján : Universidad Nacional de Luján (UNLu), 2018

Trabajo presentado para la obtención del título de Doctor. Carrera de Doctorado en Ciencias Aplicadas, Universidad Nacional de Luján.
Enlace para descargar en texto completo desde el Repositorio Institucional de la Universdidad Nacionalo de Luján URL: http://ri.unlu.edu.ar/xmlui/handle/rediunlu/558

La estimación de variables del ciclo hidrológico mediante el uso de datos de sensores remotos ha experimentado notorios avances en los últimos años. Hoy en día es posible acceder libremente a estimaciones de lluvia, evapotranspiración, humedad del suelo, tipo de cobertura superficial y altimetría, entre otras variables de interés hidrológico. Esto ha impactado en el ámbito de la modelación hidrológica, ya sea con fines de simulación o pronóstico. Así, por ejemplo, las tasas o montos acumulados estimados de lluvia o evapotranspiración se utilizan como forzantes en distintos esquemas de modelación. Además, tanto la información sobre cobertura del suelo como altimétrica suele usarse para estimar distintos parámetros asociados a la respuesta hidrológica (e.g. número de curva) o morfológicos (e.g. relaciones hortonianas), involucrados en la formulación de muchos modelos. El presente trabajo de investigación indaga y explota la asociación entre la dinámica hidrológica de la superficie anegada o saturada y el caudal, en grandes sistemas hídricos de llanura (áreas de aporte >= 2000 km²), de clima predominantemente húmedo (pudiendo presentar marcada variabilidad interanual) y con morfología fluvial típica (red integrada y jerarquizada), postulando una alternativa metodológica para la estimación numérica de la descarga líquida haciendo uso de información obtenida mediante Teledetección, aplicando algoritmos y principios relativamente simples. En primer lugar, para este tipo de sistemas hídricos están caracterizados por suaves pendientes topográficas y procesos de generación de escorrentía vinculados a la saturación del perfil de suelo (excedente hídrico por saturación del perfil de suelo). En consecuencia, se supuso a priori un vínculo significativo entre el valor del área anegada o saturada del sistema de valles y el caudal líquido en la sección de salida, constituyendo una de las hipótesis centrales del trabajo. La hipótesis subyacente es la existencia de un vínculo directo entre el volumen de agua almacenado en el sistema fluvial y dicho valor de área (curva almacenamiento – área anegada). Luego, si existe relación entre este volumen almacenado y el caudal (curva de almacenamiento – caudal) se deduce que esta descarga debe estar asociada al valor del área anegada (curva área anegada – caudal). Asimismo, mediante Teledetección pueden obtenerse estimaciones del valor de área anegada a nivel de sistema hídrico (cuenca vertiente). Ciertamente, debido a la gran variabilidad espacial presente en la dinámica de expansión y contracción de la superficie anegada, propia de los grandes sistemas fluviales de llanura húmeda, tanto como a las dimensiones de estos, es posible utilizar datos de mediana resolución espacial, con resolución temporal adecuada (i.e. diaria o semanal) para generar series temporales con densidad suficiente, a fin de obtener una señal representativa. Así pues, la metodología desarrollada consistió en valerse de datos de reflectancia de superficie, brindados por el sensor MODIS, tomando por caso de estudio a la cuenca vertiente del río Gualeguay en Rosario del Tala (Alto y Medio Gualeguay). Este es un gran sistema de llanura húmeda, con morfología fluvial típica y una dinámica de área anegada acorde a la resolución del sensor MODIS. Para el mismo se dispone de una extensa serie de nivel hidrométrico y se cuenta con datos de aforos líquidos, los cuales son necesarios para establecer una curva de descarga en la sección de control seleccionada (Rosario del Tala). Luego, mediante la aplicación de una técnica de clasificación basada en un valor umbral inferior del 'Índice Normalizado de Agua' (NDWI), elaborado sobre datos MODIS MOD09A1, se obtuvo una serie temporal del valor estimado de área anegada para la cuenca del Alto y Medio Gualeguay. Previamente, se evaluó la eficacia del NDWI como variable asociada a la presencia de agua en superficie, se seleccionó una combinación de bandas apropiada y se estableció un valor umbral inferior. Para esto, se realizó un experimento numérico estudiando la estructura de asociación entre la fracción cubierta por agua en un píxel y el valor de NDWI adquirido, permutando aleatoriamente las combinaciones de coberturas de suelo consideradas (agua, vegetación, suelos y construcciones humanas). Se utilizaron datos de la Biblioteca Espectral ASTER y se asumió un modelo aditivo para la generación de valores sintéticos de la reflectancia del píxel en cada banda MODIS MOD09A1. Particularmente, pudo mostrarse que adoptando la combinación de bandas 3 y 7 del producto MOD09A1, un umbral NDWI inferior fijado en 0 permite discriminar aquellos píxeles completamente anegados o próximos a esta situación. En otras palabras, la discriminación de la superficie de agua mediante esta técnica constituye un método conservador, que bien puede disminuir el error de comisión a expensas de subestimación. Aún así, se postuló que el valor agregado estimado a nivel cuenca vertiente puede estar fuertemente asociado al 'valor real' constituyendo un dato proxy y, de ahí, es posible su utilización para el contraste de la hipótesis de asociación con el caudal. Nuevamente, a fin de sostener esta hipótesis de dato proxy, se evaluó la representatividad de la organización espaciotemporal del índice, en relación al comportamiento esperable de la superficie anegada para un sistema fluvial típico con las características del Alto y Medio Gualeguay. En el mismo sentido, se exploró la relación entre la organización espacial del MODIS NDWI y la dinámica real de la superficie anegada mediante la formulación de un conjunto de consecuencias observacionales propias de la asunción de asociación entre la dinámica de la organización espaciotemporal del campo MODIS NDWI37 y dinámica real (esperable) de la superficie anegada. Dado que la evidencia empírica señaló que esta aproximación era sostenible, se procedió a utilizar la serie temporal obtenida por aplicación de umbral NDWI al conjunto de datos MOD09A1. Así, fue posible analizar la estructura de asociación entre el valor de área anegada y el caudal del Alto y Medio Gualeguay, constatándose un vínculo lineal. Por último, esto fue incorporado en la formulación de un modelo matemático de transformación de lluvia en escorrentía (paso de cálculo diario, agregado espacialmente). La estructura conceptual/procedimental adoptada en el modelo matemático fue diseñada de forma tal que pudiese simular en forma explícita la dinámica del área anegada (reproducir una serie temporal del valor agregado a nivel cuenca vertiente, a paso diario), respetando la forma del vínculo observado con el caudal. A saber, el esquema considera como forzante a la precipitación y los procesos simulados generan información numérica tanto de caudal como de área anegada. Los parámetros del modelo matemático fueron calibrados sobre la base de la información extraída mediante Teledetección, esto es: minimizando el error entre la señal del área anegada simulada al respecto de una muestra de los valores obtenidos mediante clasificación por técnica de umbral MODIS NDWI (considerando a estos valores como observados). Por otro lado, la evaluación (verificación) se realizó contrastando el caudal simulado contra los valores observados disponibles, apreciándose índices de eficiencia satisfactorios (Rnashsutcliffe próximo a 0.7) y con conducta estable (sin tendencia temporal alguna en estos indicadores) para todo el período considerado. También, se evaluó la estabilidad temporal de los valores de los parámetros del ajuste entre las series temporales de área anegada simulada y observada sin observar variaciones significativas. Por tanto, la evidencia empírica muestra que la estructura de la asociación ha permanecido estacionaria y, de ahí, que el vínculo área anegada – caudal puede explotarse. Así, se muestra que en un contexto de escasez de información de campo,es posible obtener caudales confiables mediante la calibración de un modelo de transformación de lluvia en escorrentía sobre la base de datos provistos por Teledetección, de forma simple y rápida, bajo la asunción de un conjunto pequeño de hipótesis. Pudo desarrollarse, evaluarse e implementarse una metodología de modelación hidrológica, apoyada sobre la base de Teledetección, con gran valor para el caso de sistemas hídricos pobremente medidos, tanto en contexto de evaluación del recurso hídrico como en el la gestión de emergencias, particularmente para las tareas de monitoreo y pronóstico hidrológico.

1. AGUA; 2. HIDROLOGIA; 3. CIENCIAS APLICADAS; 4. SUELO; 5. SUELOS-ANALISIS I. Goniadzki, Ora

(1) Inv.: 13.11817 S.T.: 08:556/Cj.Doc.Ap.30