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.
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
-
more » « less
-
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. There is unambiguous evidence that glaciers have retreated from their 19th century positions, but it is less clear how far glaciers have retreated relative to their long-term Holocene fluctuations. Glaciers in western North America are thought to have advanced from minimum positions in the Early Holocene to maximum positions in the Late Holocene. We assess when four North American glaciers, located between 38–60∘ N, were larger or smaller than their modern (2018–2020 CE) positions during the Holocene. We measured 26 paired cosmogenic in situ 14C and 10Be concentrations in recently exposed proglacial bedrock and applied a Monte Carlo forward model to reconstruct plausible bedrock exposure–burial histories. We find that these glaciers advanced past their modern positions thousands of years apart in the Holocene: a glacier in the Juneau Icefield (BC, Canada) at ∼2 ka, Kokanee Glacier (BC, Canada) at ∼6 ka, and Mammoth Glacier (WY, USA) at ∼1 ka; the fourth glacier, Conness Glacier (CA, USA), was likely larger than its modern position for the duration of the Holocene until present. The disparate Holocene exposure–burial histories are at odds with expectations of similar glacier histories given the presumed shared climate forcings of decreasing Northern Hemisphere summer insolation through the Holocene followed by global greenhouse gas forcing in the industrial era. We hypothesize that the range in histories is the result of unequal amounts of modern retreat relative to each glacier's Holocene maximum position, rather than asynchronous Holocene advance histories. We explore the influence of glacier hypsometry and response time on glacier retreat in the industrial era as a potential cause of the non-uniform burial durations. We also report mean abrasion rates at three of the four glaciers: Juneau Icefield Glacier (0.3±0.3 mm yr−1), Kokanee Glacier (0.04±0.03 mm yr−1), and Mammoth Glacier (0.2±0.2 mm yr−1).
-
Spatiotemporal Drivers of Hydrochemical Variability in a Tropical Glacierized Watershed in the Andesmore » « less
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.
-
Abstract Lemon Creek Glacier, a temperate valley glacier in the Juneau Icefield of Southeast Alaska, is the site of long running (>60 years) glaciological studies. However, the most recent published estimates of its thickness and subglacial topography come from two ~50 years old sources that are not in agreement and do not account for the effects of years of negative mass balance. We collected a 1-km long active-source seismic line on the upper section of the glacier parallel and near to the centerline of the glacier, roughly straddling the equilibrium-line altitude. We used these data to perform joint reflection-refraction velocity modeling and reflection imaging of the glacier bed. We find that this upper section of Lemon Creek Glacier is as much as 150 m (~65%) thicker than previously suggested with a large overdeepening in an area previously believed to have a uniform thickness. Our results lead us to reinterpret the impact of basal motion on ice flow and have a significant impact on expectations of subglacial hydrology. We suggest that further efforts to develop a whole-glacier model of subglacial topography are necessary to support studies that require accurate models of ice thickness and subglacial topography.more » « less
