Climate models predict amplified warming at high elevations in low latitudes, making tropical glacierized regions some of the most vulnerable hydrological systems in the world. Observations reveal decreasing streamflow due to retreating glaciers in the Andes, which hold 99% of all tropical glaciers. However, the timescales over which meltwater contributes to streamflow and the pathways it takes-surface and subsurface-remain uncertain, hindering our ability to predict how shrinking glaciers will impact water resources. Two major contributors to this uncertainty are the sparsity of hydrologic measurements in tropical glacierized watersheds and the complication of hydrograph separation where there is year-round glacier melt. We address these challenges using a multi-method approach that employs repeat hydrochemical mixing model analysis, hydroclimatic time series analysis, and integrated watershed modeling. Each of these approaches interrogates distinct timescale relationships among meltwater, groundwater, and stream discharge. Our results challenge the commonly held conceptual model that glaciers buffer discharge variability. Instead, in a subhumid watershed on Volcán Chimborazo, Ecuador, glacier melt drives nearly all the variability in discharge (Pearson correlation coefficient of 0.89 in simulations), with glaciers contributing a broad range of 20%-60% or wider of discharge, mostly (86%) through surface runoff on hourly timescales, but also through infiltration that increases annual groundwater contributions by nearly 20%. We further found that rainfall may enhance glacier melt contributions to discharge at timescales that complement glacier melt production, possibly explaining why minimum discharge occurred at the study site during warm but dry El Niño conditions, which typically heighten melt in the Andes. Our findings caution against extrapolations from isolated measurements: stream discharge and glacier melt contributions in tropical glacierized systems can change substantially at hourly to interannual timescales, due to climatic variability and surface to subsurface flow processes.
Bibliographical noteFunding Information:
Acknowledgements. Funding from NSF (EAR-1759071) helped support this work. Rachel T. McLaughlin received travel support from the University of Minnesota’s Walter H. Judd Fellowship. G.-H. Crystal Ng and Andrew D. Wickert received financial support from the University of Minnesota. Jeff La Frenierre received financial support from the National Science Foundation’s Doctoral Dissertation Research Improvement Grant (1103235), the Fulbright Commission of Ecuador, the Geological Society of America, and Gustavus Adolphus College. The authors would like to acknowledge field assistance from Chad Sandell, Casey Decker, Helen Thompson, and Abigail Michels; lab support from Jeff Jeremiason, Chris Harmes, and Scott Alexander; and land-cover mapping assistance from Josh Zoellmer.