Analysis of glacier hydrological model components with applications to Chile's Maipo Basin
Authors: Ahmad Mumtaz , Rahman Mushfiqur 
Published: June 2023

During the dry season, glaciers contribute significantly to the water supply downstream, especially in areas with seasonal rainfall. It is important to have physical-based glacier models, which are more advanced in quantifying the future dynamics of glaciers and melt distribution. Toward this purpose, an investigation of different numerical modeling tools and their approach is conducted and planned further to incorporate a glacier module into a spatially distributed model called mHM (Mesoscale Hydrological Model, mhm-ufz.org). An open-source tool called MATILDA (Modeling Water Resources in Glacierized Catchments) is used beforehand to implement the glacier entities to the application catchment Maipo Basin in Chile. Due to limited time availability and source of data, the glacier melt of the MATILDA model is integrated into the mHM model to see what the implications would be if the glacier module was implemented from scratch in mHM. The analysis of the glacier hydrology resulting from the Maipo-Glacier model was carried out by comparing the modeled glacier flows with the total observed flows. To assess the credibility of the models, the results are compared based on three goodness-fit measures, specifically r^2Coefficient of Determination), NSE (Nash-Sutcliffe Efficiency), and KGE (Kling-Gupta Efficiency). Application of the mHM model in the Maipo basin in Central Chile is performed by using prepared input datasets, showing that the approach is viable as evidenced by a Kling-Gupta Efficiency of 0.80 and Nash-Sutcliffe Efficiency of 0.70 from the starting periods of 1950 to 2020. Overestimation by the mHM model occurs during the late summer season, particularly when temperature decreases. The MATILDA model is used on glacier melt distribution, forcing data is prepared from regional datasets, and a glacier profile is constructed from ice thickness datasets. According to the model, the average KGE result for the years 2012&ndash;14 and 2015&ndash;18 was respectively 0.44 and 0.41. Glacier melt produced by Matilda is in the range of 58-68 % of the total runoff and the glacier area decreased by 22-23% during the simulation period from 2012-14 and 2016-2018. A pilot sub-basin Olivares is used in the coupling mHM model with initializing MATILDA glacier melt. The new model is underestimated during the seasonal transition in the initial days of the simulation. A considerable portion of the ablation period might be missed as a result of the initial storage quality, true glacier extent, and height distribution being understated. As a conclusion, the mHM Model was capable of predicting runoff well in a mountainous region and could also serve as a monitoring tool for watershed management. There are differences in results between the two models because of the different spatial resolutions and methodologies. MATILDA have a glacier retreat routine that gives better result for the glacier characteristics such as. A new glacier dynamics model needs to be implemented in the future to develop intermediate complexity, bridge catchment, and glacier scales after using the mHM, MATILDA, and mHM with MATIDLA model for analysis of glacier and highlighting the model's deficiencies in the initial period of simulative assessments.