Climate change has the potential to impact headwater streams in the Arctic by thawing permafrost and subsequently altering hydrologic regimes and vegetation distribution, physiognomy and productivity. Permafrost thaw and increased subsurface flow have been inferred from the chemistry of large rivers, but there is limited empirical evidence of the impacts to headwater streams. Here we demonstrate how changing vegetation cover and soil thaw may alter headwater catchment hydrology using water budgets, stream discharge trends, and chemistry across a gradient of ground temperature in northwestern Alaska. Colder, tundra-dominated catchments shed precipitation through stream discharge, whereas in warmer catchments with greater forest extent, evapotranspiration (ET) and infiltration are substantial fluxes. Forest soils thaw earlier, remain thawed longer, and display seasonal water content declines, consistent with greater ET and infiltration. Streambed infiltration and water chemistry indicate that even minor warming can lead to increased infiltration and subsurface flow. Additional warming, permafrost loss, and vegetation shifts in the Arctic will deliver water back to the atmosphere and to subsurface aquifers in many regions, with the potential to substantially reduce discharge in headwater streams, if not compensated by increasing precipitation. Decreasing discharge in headwater streams will have important implications for aquatic and riparian ecosystems.
Upland permafrost regions occupy approximately one third of the Arctic landscape. In upland regions, hydrologic fluxes are influenced by water tracks, curvilinear features on hillslopes that preferentially fill with and route water in response to snowmelt and rainfall when the soil above continuous permafrost thaws in the summer. As continued warming of the Arctic may alter hydrologic cycling leading to increased frequency of extreme hydrologic events like drought and flooding as well as modification of biogeochemical cycling, it is imperative to untangle the interplay between precipitation, runoff, and subsurface flow as water is routed from upland Arctic regions to the Arctic Ocean. This study quantifies how ground surface temperatures affect groundwater discharge from hillslopes with water tracks in the upland Arctic by employing a three‐dimensional, physically based subsurface flow model with variable saturation and freeze and thaw capabilities that is calibrated to field measurements from the Upper Kuparuk River watershed on the North Slope of Alaska, USA. Model analysis indicates that higher ground surface temperatures along water track hillslopes promote increases in groundwater discharge where water tracks act as conduits for large‐recharge events and continue to discharge groundwater into the autumn after the adjacent hillslope has frozen. Simulating the conditions that distinguish water tracks from their hillslope watersheds changes subsurface water storage and ground thermal responses but does not alter the total magnitude of groundwater discharge outside of parameter uncertainty. These findings suggest that water tracks play a complex and critical role in hydrologic cycles of the upland Arctic.
more » « less- NSF-PAR ID:
- 10360094
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Earth Surface
- Volume:
- 125
- Issue:
- 2
- ISSN:
- 2169-9003
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Climate change threatens water resources in snowmelt‐dependent regions by altering the fraction of snow and rain and spurring an earlier snowmelt season. The bulk of hydrological research has focused on forecasting response in streamflow volumes and timing to a shrinking snowpack; however, the degree to which subsurface storage offsets the loss of snow storage in various alpine geologic settings, i.e. the hydrogeologic buffering capacity, is still largely unknown. We address this research need by assessing the affects of climate change on storage and runoff generation for two distinct hydrogeologic settings present in alpine systems: a low storage granitic and a greater storage volcanic hillslope. We use a physically based integrated hydrologic model fully coupled to a land surface model to run a base scenario and then three progressive warming scenarios, and account for the shifts in each component of the water budget. For hillslopes with greater water retention, the larger storage volcanic hillslope buffered streamflow volumes and timing, but at the cost of greater reductions in groundwater storage relative to the low storage granite hillslope. We found that the results were highly sensitive to the unsaturated zone retention parameters, which in the case of alpine systems can be a mix of matrix or fracture flow. The presence of fractures and thus less retention in the unsaturated zone significantly decreased the reduction in recharge and runoff for the volcanic hillslope in climate warming scenarios. This approach highlights the importance of incorporating physically based subsurface flow in to alpine hydrology models, and our findings provide ways forward to arrive at a conceptual model that is both consistent with geology and hydrologic principles. Copyright © 2016 John Wiley & Sons, Ltd.
-
more » « less
Abstract The concept of using representative hillslopes to simulate hydrologically similar areas of a catchment has been incorporated in many hydrologic models but few Earth system models. Here we describe a configuration of the Community Land Model version 5 in which each grid cell is decomposed into one or more multicolumn hillslopes. Within each hillslope, the intercolumn connectivity is specified, and the lateral saturated subsurface flow from each column is passed to its downslope neighbor. We first apply the model to simulate a headwater catchment and assess the results against runoff and evapotranspiration flux measurements. By redistributing soil water within the catchment, the model is able to reproduce the observed difference between evapotranspiration in the upland and lowland portions of the catchment. Next, global simulations based on hypothetical hillslope geomorphic parameters are used to show the model's sensitivity to differences in hillslope shape and discretization. Differences in evapotranspiration between upland and lowland hillslope columns are found to be largest in arid and semiarid regions, while humid tropical and high‐latitude regions show limited evapotranspiration increases in lowlands relative to uplands.
-
more » « less
Abstract Groundwater discharge is an important mechanism through which fresh water and associated solutes are delivered to the ocean. Permafrost environments have traditionally been considered hydrogeologically inactive, yet with accelerated climate change and permafrost thaw, groundwater flow paths are activating and opening subsurface connections to the coastal zone. While warming has the potential to increase land-sea connectivity, sea-level change has the potential to alter land-sea hydraulic gradients and enhance coastal permafrost thaw, resulting in a complex interplay that will govern future groundwater discharge dynamics along Arctic coastlines. Here, we use a recently developed permafrost hydrological model that simulates variable-density groundwater flow and salinity-dependent freeze-thaw to investigate the impacts of sea-level change and land and ocean warming on the magnitude, spatial distribution, and salinity of coastal groundwater discharge. Results project both an increase and decrease in discharge with climate change depending on the rate of warming and sea-level change. Under high warming and low sea-level rise scenarios, results show up to a 58% increase in coastal groundwater discharge by 2100 due to the formation of a supra-permafrost aquifer that enhances freshwater delivery to the coastal zone. With higher rates of sea-level rise, the increase in discharge due to warming is reduced to 21% as sea-level rise decreased land-sea hydraulic gradients. Under lower warming scenarios for which supra-permafrost groundwater flow was not established, discharge decreased by up to 26% between 1980 and 2100 for high sea-level rise scenarios and increased only 8% under low sea-level rise scenarios. Thus, regions with higher warming rates and lower rates of sea-level change (e.g. northern Nunavut, Canada) will experience a greater increase in discharge than regions with lower warming rates and higher rates of sea-level change. The magnitude, location and salinity of discharge have important implications for ecosystem function, water quality, and carbon dynamics in coastal zones.
-
more » « less
Abstract Both hydrological and geophysical data can be used to calibrate hillslope hydrologic models. However, these data often reflect hydrological dynamics occurring at disparate spatial scales. Their use as sole objectives in model calibrations may thus result in different optimum hydraulic parameters and hydrologic model behavior. This is especially true for mountain hillslopes where the subsurface is often heterogeneous and the representative elementary volume can be on the scale of several m3. This study explores differences in hydraulic parameters and hillslope‐scale storage and flux dynamics of models calibrated with different hydrological and geophysical data. Soil water content, groundwater level, and two time‐lapse electrical resistivity tomography (ERT) data sets (transfer resistance and inverted resistivity) from two mountain hillslopes in Wyoming, USA, are used to calibrate physics‐based surface–subsurface hydrologic models of the hillslopes. Calibrations are performed using each data set independently and all data together resulting in five calibrated parameter sets at each site. Model predicted hillslope runoff and internal hydrological dynamics vary significantly depending on the calibration data set. Results indicate that water content calibration data yield models that overestimate near‐surface water storage in mountain hillslopes. Groundwater level calibration data yield models that more reasonably represent hillslope‐scale storage and flux dynamics. Additionally, ERT calibration data yield models with reasonable hillslope runoff predictions but relatively poor predictions of internal hillslope dynamics. These observations highlight the importance of carefully selecting data for hydrologic model calibration in mountain environments. Poor selection of calibration data may yield models with limited predictive capability depending on modeling goals and model complexity.
