skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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. 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
  3. Abstract. The McMurdo Dry Valleys (MDV) are home to a unique microbial ecosystem that is dependent on the availability of freshwater. This is a polar desert and freshwater originates almost entirely from surface and near-surface melt of the cold-based glaciers. Understanding the future evolution of these environments requires the simulation of the full chain of physical processes from net radiative forcing, surface energy balance, melt, runoff and transport of meltwater in stream channels from the glaciers to the terminal lakes where the microbial community resides. To establish a new framework to do this, we present the first application of WRF-Hydro/Glacier in the MDV, which as a fully distributed hydrological model has the capability to resolve the streams from the glaciers to the bare land that surround them. Given that meltwater generation in the MDV is almost entirely dependent on small changes in the energy balance of the glaciers, the aim of this study is to optimize the multi-layer snowpack scheme that is embedded in WRF-Hydro/Glacier to ensure that the feedbacks between albedo, snowfall and melt are fully resolved. To achieve this, WRF-Hydro/Glacier is implemented at a point scale using automatic weather station data on Commonwealth Glacier to physically model the onset, duration and end of melt over a 7-month period (1 August 2021 to 28 February 2022). To resolve the limited energetics controlling melt, it was necessary to (1) limit the percolation of meltwater through the ice layers in the multi-layer snowpack scheme and (2) optimize the parameters controlling the albedo of both snow and ice over the melt season based on observed spectral signatures of albedo. These modifications enabled the variability of broadband albedo over the melt season to be accurately simulated and ensured that modelled surface and near-surface temperatures, surface height change and runoff were fully resolved. By establishing a new framework that couples a detailed snowpack model to a fully distributed hydrological model, this work provides a stepping stone to model the spatial and temporal variability of melt and streamflow in the future, which will enable some of the unknown questions about the hydrological connectivity of the MDV to be answered. 
    more » « less
  4. Abstract. Heterotrophic prokaryotic production (BP) was studied in the western tropical South Pacific (WTSP) using the leucine technique, revealing spatial and temporal variability within the region. Integrated over the euphotic zone, BP ranged from 58 to 120mg Cm−2d−1 within the Melanesian Archipelago, and from 31 to 50mg Cm−2d−1 within the western subtropical gyre. The collapse of a bloom was followed during 6 days in the south of Vanuatu using a Lagrangian sampling strategy. During this period, rapid evolution was observed in the three main parameters influencing the metabolic state: BP, primary production (PP) and bacterial growth efficiency. With N2 fixation being one of the most important fluxes fueling new production, we explored relationships between BP, PP and N2 fixation rates over the WTSP. The contribution of N2 fixation rates to bacterial nitrogen demand ranged from 3 to 81%. BP variability was better explained by the variability of N2 fixation rates than by that of PP in surface waters of the Melanesian Archipelago, which were characterized by N-depleted layers and low DIP turnover times (TDIP<100h). This is consistent with the fact that nitrogen was often one of the main factors controlling BP on short timescales, as shown using enrichment experiments, followed by dissolved inorganic phosphate (DIP) near the surface and labile organic carbon deeper in the euphotic zone. However, BP was more significantly correlated with PP, but not with N2 fixation rates where DIP was more available (TDIP>100h), deeper in the Melanesian Archipelago, or within the entire euphotic zone in the subtropical gyre. The bacterial carbon demand to gross primary production ratio ranged from 0.75 to 3.1. These values are discussed in the framework of various assumptions and conversion factors used to estimate this ratio, including the methodological errors, the daily variability of BP, the bacterial growth efficiency and one bias so far not considered: the ability for Prochlorococcus to assimilate leucine in the dark. 
    more » « less
  5. Abstract. We assessed the influence of the marine diazotrophic cyanobacterium Trichodesmium on the bio-optical properties of western tropical South Pacific (WTSP) waters (18–22°S, 160°E–160°W) during the February–March 2015 OUTPACE cruise. We performed measurements of backscattering and absorption coefficients, irradiance, and radiance in the euphotic zone with a Satlantic MicroPro free-fall profiler and took Underwater Vision Profiler 5 (UPV5) pictures for counting the largest Trichodesmium spp. colonies. Pigment concentrations were determined by fluorimetry and high-performance liquid chromatography and picoplankton abundance by flow cytometry. Trichome concentration was estimated from pigment algorithms and validated by surface visual counts. The abundance of large colonies counted by the UVP5 (maximum 7093coloniesm−3) was well correlated to the trichome concentrations (maximum 2093trichomesL−1) with an aggregation factor of 600. In the Melanesian archipelago, a maximum of 4715trichomesL−1 was enumerated in pump samples (3.2m) at 20°S,16730°E. High Trichodesmium abundance was always associated with absorption peaks of mycosporine-like amino acids (330, 360nm) and high particulate backscattering, but not with high Chl a fluorescence or blue particulate absorption (440nm). Along the west-to-east transect, Trichodesmium together with Prochlorococcus represented the major part of total chlorophyll concentration; the contribution of other groups were relatively small or negligible. The Trichodesmium contribution to total chlorophyll concentration was the highest in the Melanesian archipelago around New Caledonia and Vanuatu (60%), progressively decreased to the vicinity of the islands of Fiji (30%), and reached a minimum in the South Pacific Gyre where Prochlorococcus dominated chlorophyll concentration. The contribution of Trichodesmium to zeaxanthin was respectively 50, 40 and 20% for these regions. During the OUTPACE cruise, the relationship between normalized water-leaving radiance (nLw) in the ultraviolet and visible and chlorophyll concentration was similar to that found during the BIOSOPE cruise in the eastern tropical Pacific. Principal component analysis (PCA) of OUTPACE data showed that nLw at 305, 325, 340, 380, 412 and 440nm was strongly correlated to chlorophyll and zeaxanthin, while nLw at 490 and 565nm exhibited lower correlations. These results, as well as differences in the PCA of BIOSOPE data, indicated that nLw variability in the greenish blue and yellowish green during OUTPACE was influenced by other variables associated with Trichodesmium presence, such as backscattering coefficient, phycoerythrin fluorescence and/or zeaxanthin absorption, suggesting that Trichodesmium detection should involve examination of nLw in this spectral domain. 
    more » « less