Comportamiento de un modelo hidrológico distribuido frente a la incertidumbre de la temperatura del aire

Abstract

La baja cantidad de estaciones de monitoreo es la principal dificultad cuando se estudia la dinámica del manto nival en cuencas de la zona central de Chile. Debido a esto, se propone estimar la temperatura del aire (factor influyente en el balance energético y en la distribución del manto nival) a través de imágenes satelitales MODIS e información de estaciones meteorológicas. Para estudiar los efectos de la metodología propuesta para el cálculo de la temperatura del aire, se calibran siete y validan cuatro años (periodo comprendido entre abril de 2001 marzo 2012), utilizando el modelo hidrológico VIC en las cuencas definidas por las estaciones hidrometeorológicas Aconcagua en Chacabuquito, Juncal en Juncal, y Colorado en Colorado. Los resultados encontrados en la calibración se consideran aceptables (Eficiencia de Nash-Sutcliffe > 0,76), lo mismo pasa en la validación, ya que se mantiene la misma tendencia con valores de Nash-Sutcliffe > 0,65. A continuación, dado que la metodología en el cálculo de la temperatura del aire a través de imágenes satelitales tiene un error asociado, se generan 90 series sintéticas de temperatura del aire que tienen los mismos estadísticos que la serie original. Estos resultados se incluyen dentro del modelo hidrológico durante todo el periodo, representando el caudal modelado como una banda de incertidumbre. Con dicha banda, la eficiencia Nash Sutcliffe aumenta notoriamente, con valores sobre 0,78 para todas las cuencas. Una vez que se obtiene la banda de incertidumbre, se efectúa el pronóstico estacional para la cuenca Aconcagua en Chacabuquito en el período de deshielo comprendido entre el 1 de septiembre de 2012 y el 31 de marzo de 2013, con las respectivas actualizaciones hasta el 1 de octubre y 1 noviembre de 2012. Los tres pronósticos se comparan con el tradicional entregado por Dirección General de Aguas (DGA). Los resultados entregados por VIC y por la DGA se comportan de manera similar y consistente con lo observado hasta cada fecha de actualización. En el caso de VIC, el caudal observado siempre es incluido dentro la banda de incertidumbre. Finalmente, se utilizan las imágenes satelitales MODIS de cobertura nival y se comparan con el equivalente en agua entregado por el modelo VIC. Ambos resultados muestran una dinámica del derretimiento muy similar, incluso se logra presenciar nieve que no se derrite de un año a otro, la que podría estar siendo incorporada al sistema en forma de glaciar.