An official website of the United States government
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
Spatiotemporal Drivers of Hydrochemical Variability in a Tropical Glacierized Watershed in the Andes
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
- PAR ID:
- 10442534
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 57
- Issue:
- 5
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Terrestrial production and export of dissolved organic and inorganic carbon (DOC and DIC) to streams depends on water flow and biogeochemical processes in and beneath soils. Yet, understanding of these processes in a rapidly changing climate is limited. Using the watershed‐scale reactive‐transport model BioRT‐HBV and stream data from a snow‐dominated catchment in the Rockies, we show deeper groundwater flow averaged about 20% of annual discharge, rising to ∼35% in drier years. DOC and DIC production and export peaked during snowmelt and wet years, driven more by hydrology than temperature. DOC was primarily produced in shallow soils (1.94 ± 1.45 gC/m2/year), stored via sorption, and flushed out during snowmelt. Some DOC was recharged to and further consumed in the deeper subsurface via respiration (−0.27 ± 0.02 gC/m2/year), therefore reducing concentrations in deeper groundwater and stream DOC concentrations at low discharge. Consequently, DOC was primarily exported from the shallow zone (1.62 ± 0.96 gC/m2/year, compared to 0.12 ± 0.02 gC/m2/year from the deeper zone). DIC was produced in both zones but at higher rates in shallow soils (1.34 ± 1.00 gC/m2/year) than in the deep subsurface (0.36 ± 0.02 gC/m2/year). Deep respiration elevated DIC concentrations in the deep zone and stream DIC concentrations at low discharge. In other words, deep respiration is responsible for the commonly‐observed increasing DOC concentrations (flushing) and decreasing DIC concentrations (dilution) with increasing discharge. DIC export from the shallow zone was ~66% of annual export but can drop to ∼53% in drier years. Numerical experiments suggest lower carbon production and export in a warmer, drier future, and a higher proportion from deeper flow and respiration processes. These results underscore the often‐overlooked but growing importance of deeper processes in a warming climate.more » « less
-
Abstract As glaciers around the world rapidly lose mass, the tight coupling between glaciers and downstream ecosystems is resulting in widespread impacts on global hydrologic and biogeochemical cycling. However, a range of challenges make it difficult to conduct research in glacierized systems, and our knowledge of seasonally changing hydrologic processes and solute sources and signatures is limited. This in turn hampers our ability to make predictions on solute composition and flux. We conducted a broad water sampling campaign in order to understand the present‐day partitioning of water sources and associated solutes in Alaska's Wolverine Glacier watershed. We established a relationship between electrical conductivity and streamflow at the watershed outlet to divide the melt season into four hydroclimatic periods. Across hydroclimatic periods, we observed a shift in nonglacial source waters from snowmelt‐dominated overland and shallow subsurface flow paths to deeper groundwater flow paths. We also observed the shift from a low‐ to high‐efficiency subglacial drainage network and the associated flushing of water stored subglacially with higher solute loads. We used calcium, the dominant dissolved ion, from watershed outlet samples to estimate solute fluxes for each hydroclimatic period across two melt seasons. We found between 40% and 55% of Ca2+export occurred during the late season rainy period. This partitioning of the melt season coupled with a characterization of the chemical makeup and magnitude of solute export provides new insight into a rapidly changing watershed and creates a framework to quantify and predict changes to solute fluxes across a melt season.more » « less
-
Abstract Streamflow generation in mountain watersheds is strongly influenced by snow accumulation and melt as well as groundwater connectivity. In mountainous regions with limestone and dolomite geology, bedrock formations can host karst aquifers, which play a significant role in snowmelt–discharge dynamics. However, mapping complex karst features and the resulting surface‐groundwater exchanges at large scales remains infeasible. In this study, timeseries analysis of continuous discharge and specific conductance measurements were combined with gridded snowmelt predictions to characterize seasonal streamflow response and evaluate dominant watershed controls across 12 monitoring sites in a karstified 554 km2watershed in northern Utah, USA. Immense surface water hydrologic variability across subcatchments, years and seasons was linked to geologic controls on groundwater dynamics. Unlike many mountain watersheds, the variability between subcatchments could not be well described by typical watershed properties, including elevation or surficial geology. To fill this gap, a conceptual framework was proposed to characterize subsurface controls on snowmelt–discharge dynamics in karst mountain watersheds in terms of conduit flow direction, aquifer storage capacity and connectivity. This framework requires only readily measured surface water and climatic data from nested monitoring sites and was applied to the study watershed to demonstrate its applicability for evaluating dominant controls and climate sensitivity.more » « less
-
Abstract Biospheric particulate organic carbon (POCbio) burial and rock petrogenic particulate organic carbon (POCpetro) oxidation are opposing long‐term controls on the global carbon cycle, sequestering and releasing carbon, respectively. Here, we examine how watershed glacierization impacts the POC source by assessing the concentration and isotopic composition (δ13C and Δ14C) of POC exported from four watersheds with 0%–49% glacier coverage across a melt season in Southeast Alaska. We used two mixing models (age‐weight percent and dual carbon isotope) to calculate concentrations of POCbioand POCpetrowithin the bulk POC pool. The fraction POCpetrocontribution was highest in the heavily glacierized watershed (age‐weight percent: 0.39 ± 0.05; dual isotope: 0.42 (0.37–0.47)), demonstrating a glacial source of POCpetroto fjords. POCpetrowas mobilized via glacier melt and subglacial flow, while POCbiowas largely flushed from the non‐glacierized landscape by rain. Flow normalized POCbioconcentrations exceeded POCpetroconcentrations for all streams, but surprisingly were highest in the heavily glacierized watershed (mean: 0.70 mgL−1; range 0.16–1.41 mgL−1), suggesting that glacier rivers can contribute substantial POCbioto coastal waters. Further, the most heavily glacierized watershed had the highest sediment concentration (207 mgL−1; 7–708 mgL−1), and thus may facilitate long‐term POCbioprotection via sediment burial in glacier‐dominated fjords. Our results suggest that continuing glacial retreat will decrease POC concentrations and increase POCbio:POCpetroexported from currently glacierized watersheds. Glacier retreat may thus decrease carbon storage in marine sediments and provide a positive feedback mechanism to climate change that is sensitive to future changes in POCpetrooxidation.more » « less
