Representación de la hidrología glaciar de la Cuenca del río Maipo a través de un modelo hidrológico simplificado

Abstract

Los Glaciares cumplen un importante rol en la regulación del régimen hidrológico de cuencas montañosas en climas mediterráneos, puesto que son capaces de almacenar un gran volumen de agua dulce y entregarlo en los meses estivales. Gran parte de los glaciares, sin embargo, han ido reduciendo su superficie considerablemente en los últimos años (efecto potenciado por los aumentos de temperatura producto del cambio climático), fenómeno apreciable a lo largo de todo Chile, por lo que es de gran importancia conocer el estado actual y futuro de estas masas de hielo. Se implementó un modelo hidrológico semi-distribuido, desarrollado en el software de modelamiento "WEAP" cuyo objetivo es simular el comportamiento glaciar de las subcuencas de los ríos Olivares, Volcán y Maipo en las Hualtatas, en la cuenca alta del río Maipo. El software incluye un módulo especial para glaciares que permite aplicar ecuaciones de acumulación y derretimiento de hielo, tanto para glaciares cubiertos como descubiertos, para cada una de las URH (Unidad de respuesta Hidrológica) de estudio, que en este caso corresponden a bandas de elevación separadas cada 500 metros de altura, para cada subcuenca. En su formulación original, el modelo no entrega resultados representativos para todas las subcuencas de estudio, entregando una tasa de derretimiento menor a la observada en terreno para estas subcuencas. Una de las hipótesis que explican estos resultados es que las temperaturas ingresadas al modelo no son representativas de la realidad. En un comienzo se obtuvieron las temperaturas promedio de cada una de las URH mediante el método del gradiente térmico, aproximación generalmente aceptada para este tipo de estudios. Sin embargo se omiten muchas variables al aplicar este método, por lo que en el presente estudio se utilizan datos de percepción remota para obtener una representación más confiable de la realidad. Además se revaluaron las ecuaciones que rigen el modelo original, dando énfasis a aquellas que correlacionan la superficie glaciar con su volumen, de manera de corroborar los supuestos que fueron establecidos en un comienzo. Una vez aplicadas las correcciones mencionadas anteriormente se procedió a recalibrar el modelo,incorporando nuevos datos fluviométricos y meteorológicos obtenidos en los últimos años. Como resultado se observó una pérdida aproximada de 20[km2] de superficies glaciares descubiertas y de 6[km2] de superficies glaciares cubiertas durante el período de modelación (2001-2014), equivalentes a una pérdida de 1.53[km2/año] y 0.46[km2/año] respectivamente. Se obtuvo además que para años considerados secos, el aporte proveniente del derretimiento de glaciares descubiertos a la escorrentía total puede llegar a ser del orden de un 70 [ %]. Finalmente se logró obtener un diagnóstico fiable del estado actual de los glaciares de las cuencas de estudio, así como la contribución de estos a los caudales de ríos de montaña.