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.
With a unique biogeophysical signature relative to other freshwater sources, meltwater from glaciers plays a crucial role in the hydrological and ecological regime of high latitude coastal areas. Today, as glaciers worldwide exhibit persistent negative mass balance, glacier runoff is changing in both magnitude and timing, with potential downstream impacts on infrastructure, ecosystems, and ecosystem resources. However, runoff trends may be difficult to detect in coastal systems with large precipitation variability. Here, we use the coupled energy balance and water routing model SnowModel‐HydroFlow to examine changes in timing and magnitude of runoff from the western Juneau Icefield in Southeast Alaska between 1980 and 2016. We find that under sustained glacier mass loss (−0.57 ± 0.12 m w. e. a−1), several hydrological variables related to runoff show increasing trends. This includes annual and spring glacier ice melt volumes (+10% and +16% decade−1) which, because of higher proportions of precipitation, translate to smaller increases in glacier runoff (+3% and +7% decade−1) and total watershed runoff (+1.4% and +3% decade−1). These results suggest that the western Juneau Icefield watersheds are still in an increasing glacier runoff period prior to reaching “peak water.” In terms of timing, we find that maximum glacier ice melt is occurring earlier (2.5 days decade−1), indicating a change in the source and quality of freshwater being delivered downstream in the early summer. Our findings highlight that even in maritime climates with large precipitation variability, high latitude coastal watersheds are experiencing hydrological regime change driven by ongoing glacier mass loss.
more » « less- NSF-PAR ID:
- 10374448
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 57
- Issue:
- 7
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Spatiotemporal Drivers of Hydrochemical Variability in a Tropical Glacierized Watershed in the Andesmore » « less
Abstract -
more » « less
Abstract The west Antarctic Peninsula (WAP) is a region of marked climatic variability, exhibiting strong changes in sea ice extent, retreat of most of its glaciers, and shifts in the amount and form of precipitation. These changes can have significant impacts on the oceanic freshwater budget and marine biogeochemical processes; it is thus important to ascertain the relative balance of the drivers and the spatial scales over which they operate. We present a novel 7‐year summer‐season (October to March; 2011 to 2018) series of oxygen isotopes in seawater (δ18O), augmented with some winter sampling, collected adjacent to Anvers Island at the WAP. These data are used to attribute oceanic freshwater changes to sea ice and meteoric sources, and to deduce information on the spatial scales over which the changes are driven. Sea ice melt shows significant seasonality (∼9% range) and marked interannual changes, with pronounced maxima in seasons 2013/14 and 2016/17. Both of these extrema are driven by anomalous winds, but reflect strongly contrasting dynamic and thermodynamic sea ice responses. Meteoric water also shows seasonality (∼7% range) with interannual variability reflecting changes in the input of accumulated precipitation and glacial melt to the ocean. Unlike sea ice melt, meteoric water extremes are especially pronounced in thin (<10 m) surface layers close to the proximate glacier, associated with enhanced ocean stratification. Isotopic tracers help to deconvolve the complex spatio‐temporal scales inherent in the coastal freshwater budget, and hence improve our knowledge of the separate and cumulative physical and ecological impacts.
-
more » « less
Abstract The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed‐watershed, perennially ice‐covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. We apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice‐covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996–2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2–17 m higher, depending on the lake, relative to the 2020 level.
-
more » « less
Abstract At tidewater glacier termini, ocean‐glacier interactions hinge on two sources of freshwater—submarine melt and subglacial discharge—yet these freshwater fluxes are often unconstrained in their magnitude, seasonality, and relationship. With measurements of ocean velocity, temperature and salinity, fjord budgets can be evaluated to partition the freshwater flux into submarine melt and subglacial discharge. We apply these methods to calculate the freshwater fluxes at LeConte Glacier, Alaska, across a wide range of oceanic and atmospheric conditions during six surveys in 2016–2018. We compare these ocean‐derived fluxes with an estimate of subglacial discharge from a surface mass balance model and with estimates of submarine melt from multibeam sonar and autonomous kayaks, finding relatively good agreement between these independent estimates. Across spring, summer, and fall, the relationship between subglacial discharge and submarine melt follows a scaling law predicted by standard theory (melt ∼ discharge1/3), although the total magnitude of melt is an order of magnitude larger than theoretical estimates. Subglacial discharge is the dominant driver of variability in melt, while the dependence of melt on fjord properties is not discernible. A comparison of oceanic budgets with glacier records indicates that submarine melt removes 33%–49% of the ice flux into the terminus across spring, summer, and fall periods. Thus, melt is a significant component of the glacier's mass balance, and we find that melt correlates with seasonal retreat; however, melt does not appear to directly amplify calving.
-
more » « less
Abstract Dissolved organic matter (DOM) in glacier runoff is aliphatic‐rich, yet studies have proposed that DOM originates mainly from allochthonous, aromatic, and often aged material. Allochthonous organic matter (OM) is exposed to ultraviolet radiation both in atmospheric transport and post‐deposition on the glacier surface. Thus, we evaluate photochemistry as a mechanism to account for the compositional disconnect between allochthonous OM sources and glacier runoff DOM composition. Six endmember OM sources (including soils and diesel particulate matter) were leached and photo‐irradiated for 28 days in a solar simulator, until >90% of initial chromophoric DOM was removed. Ultrahigh‐resolution mass spectrometry was used to compare the molecular composition of endmember leachates pre‐ and post‐irradiation to DOM in supraglacial and bulk runoff from the Greenland Ice Sheet and Juneau Icefield (Alaska), respectively. Photo‐irradiation drove molecular level convergence between the initially aromatic‐rich leachates and aromatic‐poor glacial samples, selectively removing aromatic compounds (−80 ± 19% relative abundance) and producing aliphatics (+75 ± 35% relative abundance). Molecular level glacier runoff DOM composition was statistically indistinguishable to post‐irradiation leachates. Bray‐Curtis analysis showed substantial similarity in the molecular formulae present between glacier samples and post‐irradiation leachates. Post‐irradiation leachates contained 84 ± 7.4% of the molecular formulae, including 72 ± 17% of the aliphatic formulae, detected in glacier samples. Our findings suggest that photodegradation, either in transit to or on glacier surfaces, could provide a mechanistic pathway to account for the disconnect between proposed aromatic, aged sources of OM and the aliphatic‐rich fingerprint of glacial DOM.
