skip to main content


Title: Multi-scale temporal variability in meltwater contributions in a tropical glacierized watershed

Abstract. Climate models predict amplified warming at high elevations in low latitudes,making tropical glacierized regions some of the most vulnerable hydrologicalsystems in the world. Observations reveal decreasing streamflow due toretreating glaciers in the Andes, which hold 99% of all tropicalglaciers. However, the timescales over which meltwater contributes tostreamflow and the pathways it takes – surface and subsurface – remainuncertain, hindering our ability to predict how shrinking glaciers willimpact water resources. Two major contributors to this uncertainty are thesparsity of hydrologic measurements in tropical glacierized watersheds andthe complication of hydrograph separation where there is year-round glaciermelt. We address these challenges using a multi-method approach that employsrepeat hydrochemical mixing model analysis, hydroclimatic time seriesanalysis, and integrated watershed modeling. Each of these approachesinterrogates distinct timescale relationships among meltwater, groundwater,and stream discharge. Our results challenge the commonly held conceptualmodel that glaciers buffer discharge variability. Instead, in a subhumidwatershed on Volcán Chimborazo, Ecuador, glacier melt drives nearly allthe variability in discharge (Pearson correlation coefficient of 0.89 insimulations), with glaciers contributing a broad range of 20%–60%or wider of discharge, mostly (86%) through surface runoff on hourlytimescales, but also through infiltration that increases annual groundwatercontributions by nearly 20%. We further found that rainfall may enhanceglacier melt contributions to discharge at timescales that complement glaciermelt production, possibly explaining why minimum discharge occurred at thestudy site during warm but dry El Niño conditions, which typicallyheighten melt in the Andes. Our findings caution against extrapolations fromisolated measurements: stream discharge and glacier melt contributions intropical glacierized systems can change substantially at hourly tointerannual timescales, due to climatic variability and surface to subsurfaceflow processes.

 
more » « less
Award ID(s):
1759071 1758854 1758795
NSF-PAR ID:
10085310
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Hydrology and Earth System Sciences
Volume:
23
Issue:
1
ISSN:
1607-7938
Page Range / eLocation ID:
405 to 425
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    There is a critical knowledge gap about how glacier retreat in remote and rapidly warming tropical montane watersheds will impact solute export, which has implications for downstream geochemical cycling and ecological function. Because tropical glacierized watersheds are often uniquely characterized by year‐round ablation, upslope vegetation migration, and significant groundwater flow, baseline understanding is needed of how spatiotemporal variables within these watersheds control outlet hydrochemistry. We implemented a recently developed reactive transport watershed model, BioRT‐Flux‐PIHM, for a sub‐humid glacierized watershed in the Ecuadorian Andes with young volcanic soils and fractured bedrock. We found a unique simulated concentration and discharge (C‐Q) pattern that was mostly chemostatic but superimposed by dilution episodes. The chemostatic background was attributed to large simulated contributions of groundwater (subsurface lateral flow) to streamflow, of which a notable fraction (37%) comprised infiltrated ice‐melt. Relatively constant concentrations were further maintained in the model because times and locations of lower mineral surface wetting and dissolution were offset by concentrating effects of greater evapotranspiration. Ice‐melt did not all infiltrate in simulations, especially during large precipitation events, when high surface runoff contributions to discharge triggered dilution episodes. In a model scenario without ice‐melt, major ion concentrations, including Na+, Ca2+, and Mg2+, became more strongly chemostatic and higher, but weathering rates decreased, attenuating export by 23%. We expect this reduction to be exacerbated by higher evapotranspiration and drier conditions with expanded vegetation. This work brings to light the importance of subsurface meltwater flow, ecohydrological variability, and interactions between melt and precipitation for controlling hydrochemical processes in tropical watersheds with rapidly retreating glaciers.

     
    more » « less
  2. Abstract

    Accelerating mountain glacier recession in a warming climate threatens the sustainability of mountain water resources. The extent to which groundwater will provide resilience to these water resources is unknown, in part due to a lack of data and poorly understood interactions between groundwater and surface water. Here we address this knowledge gap by linking climate, glaciers, surface water, and groundwater into an integrated model of the Shullcas Watershed, Peru, in the tropical Andes, the region experiencing the most rapid mountain‐glacier retreat on Earth. For a range of climate scenarios, our model projects that glaciers will disappear by 2100. The loss of glacial meltwater will be buffered by relatively consistent groundwater discharge, which only receives minor recharge (~2%) from glacier melt. However, increasing temperature and associated evapotranspiration, alongside potential decreases in precipitation, will decrease groundwater recharge and streamflow, particularly for the RCP 8.5 emission scenario.

     
    more » « less
  3. Rock glaciers are common geomorphic features in alpine landscapes and comprise a potentially significant but poorly quantified water resource. This project focused on three complementary questions germane to rock glacier hydrology: 1) Does the composition of rock glacier meltwater vary from year to year? 2) How dependent is the composition of rock glacier meltwater on lithology? And 3) How does the presence of rock glaciers in a catchment change stream water chemistry? To address these questions, we deployed automated samplers to collect water from late June through mid-October 2022 in two rock-glacierized mountain ranges in Utah, United States characterized by different lithologies. In the Uinta Mountains of northern Utah, where bedrock is predominantly quartzite, water was collected at springs discharging from two rock glaciers previously shown to release water in late summer sourced from internal ice. In the La Sal Mountains of southeastern Utah, where trachyte bedrock is widespread, water was collected at a rock glacier spring, along the main stream in a watershed containing multiple rock glaciers, and from a stream in a watershed where rock glaciers are absent. Precipitation was also collected, and data loggers for water temperature and electric conductivity were deployed. Water samples were analyzed for stable isotopes with cavity ring-down spectroscopy and hydrochemistry with ICP-MS. Our data show that water discharging from rock glaciers in the Uinta Mountains exhibits a shift from a snowmelt source to an internal ice source over the course of the melt season that is consistent from year to year. We also found that the chemistry of rock glacier water in the two study areas is notably different in ways that can be linked back to their contrasting bedrock types. Finally, in the La Sal Mountains, the properties of water along the main stream in a rock-glacierized basin resemble the properties of water discharging from rock glaciers, and strongly contrast with the water in a catchment lacking rock glaciers. Collectively these results underscore the role of rock glaciers as an agent influencing the hydrochemistry of water in high-elevation stream systems. 
    more » « less
  4. Abstract

    A myriad of downstream communities and industries rely on streams fed by both groundwater discharge and glacier meltwater draining the Cordillera Blanca, Northern Peruvian Andes, which contains the highest density of glaciers in the tropics. During the dry season, approximately half the discharge in the region's proglacial streams comes from groundwater. However, because of the remote and difficult access to the region, there are few field methods that are effective at the reach scale to identify the spatial distribution of groundwater discharge. An energy balance model, Rhodamine WT dye tracing, and high‐definition kite‐borne imagery were used to determine gross and net groundwater inputs to a 4‐km reach of the Quilcay River in Huascaran National Park, Peru. The HFLUX computer programme (http://hydrology.syr.edu/hflux.html) was used to simulate the Quilcay River's energy balance using stream temperature observations, meteorological measurements, and kite‐borne areal photography. Inference from the model indicates 29% of stream discharge at the reach outlet was contributed by groundwater discharge over the study section. Rhodamine WT dye tracing results, coupled with the energy balance, show that approximately 49% of stream water is exchanged (no net gain) with the subsurface as gross gains and losses. The results suggest that gross gains from groundwater are largest in a moraine subreach but because of large gross losses, net gains are larger in the meadow subreaches. These insights into pathways of groundwater–surface water interaction can be applied to improve hydrological modelling in proglacial catchments throughout South America. Copyright © 2016 John Wiley & Sons, Ltd.

     
    more » « less
  5. Abstract. Over the last century, northwestern Canada experienced some of the highest rates of tropospheric warming globally, which caused glaciers in the region to rapidly retreat. Our study seeks to extend the record of glacier fluctuations and assess climate drivers prior to the instrumental record in the Mackenzie and Selwyn mountains of northwestern Canada. We collected 27 10Be surface exposure ages across nine cirque and valley glacier moraines to constrain the timing of their emplacement. Cirque and valley glaciers in this region reached their greatest Holocene extents in the latter half of the Little Ice Age (1600–1850 CE). Four erratic boulders, 10–250 m distal from late Holocene moraines, yielded 10Be exposure ages of 10.9–11.6 ka, demonstrating that by ca. 11 ka, alpine glaciers were no more extensive than during the last several hundred years. Estimated temperature change obtained through reconstruction of equilibrium line altitudes shows that since ca. 1850 CE, mean annual temperatures have risen 0.2–2.3 ∘C. We use our glacier chronology and the Open Global Glacier Model (OGGM) to estimate that from 1000 CE, glaciers in this region reached a maximum total volume of 34–38 km3 between 1765 and 1855 CE and had lost nearly half their ice volume by 2019 CE. OGGM was unable to produce modeled glacier lengths that match the timing or magnitude of the maximum glacier extent indicated by the 10Be chronology. However, when applied to the entire Mackenzie and Selwyn mountain region, past millennium OGGM simulations using the Max Planck Institute Earth System Model (MPI-ESM) and the Community Climate System Model 4 (CCSM4) yield late Holocene glacier volume change temporally consistent with our moraine and remote sensing record, while the Meteorological Research Institute Earth System Model 2 (MRI-ESM2) and the Model for Interdisciplinary Research on Climate (MIROC) fail to produce modeled glacier change consistent with our glacier chronology. Finally, OGGM forced by future climate projections under varying greenhouse gas emission scenarios predicts 85 % to over 97 % glacier volume loss by the end of the 21st century. The loss of glaciers from this region will have profound impacts on local ecosystems and communities that rely on meltwater from glacierized catchments.

     
    more » « less