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.
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- NSF-PAR ID:
- 10458300
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 11
- Issue:
- 12
- ISSN:
- 1942-2466
- Page Range / eLocation ID:
- p. 4044-4065
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract The seasonal evolution of the ground thermal regime in cold regions influences hydrologic flow paths, soil biogeochemistry, and hillslope geomorphology. In mountain environments, steep topography produces strong gradients in solar insolation, vegetation, and snowpack dynamics that lead to large differences in soil temperature over short distances, suggesting a need for high‐resolution, process‐based models that quantify the influence of topography. We present soil temperature and snow depth results from a coupled thermo‐hydrologic model compared to field observations from Gordon Gulch, a seasonally snow‐covered montane catchment in the Colorado Front Range in the Boulder Creek Critical Zone Observatory. The field site features two instrumented hillslopes with opposing aspects: Despite the persistent snowpack on the north‐facing slope, seasonally frozen ground is more prevalent there than the south‐facing slope, which experiences significantly higher incoming radiation that prevents the persistence of frozen ground. A novel modeling framework is developed by coupling a surface energy balance model incorporating solar radiation and snowpack processes to an existing subsurface model (PFLOTRAN‐ICE). The coupled model is used to reproduce strong aspect‐controlled differences in soil temperature and snow depth evident from observations during water years 2013–2016, including a higher incidence of frozen ground under the north‐facing slope. Representation of the snowpack and its insulating effects significantly improves soil temperature estimates on the north‐facing slope, particularly the duration of soil freezing in the spring, which is underestimated by 1–2 months without including the snowpack.
-
more » « less
Abstract Soil moisture is a key control on runoff generation and biogeochemical processes on hillslopes. Precipitation events can evoke different soil moisture responses with depth through the soil profile, and responses can differ among landscape positions along a hillslope. We sought to elucidate the nature of these responses by estimating changes in water content, response time between peak precipitation and peak soil moisture, and wetting front velocities for 43 storms at 45 locations on three adjacent hillslopes within a headwater catchment of the southern Appalachian Mountains (NC, USA). We used a multivariate modeling approach to quantify the relative influences and the predictability of soil moisture responses by a combination of landscape and storm characteristics. We quantified the lag correlations between hillslope mean soil moisture and catchment runoff to demonstrate how storm properties and hillslope‐scale characteristics may influence runoff at the catchment outlet. Soil moisture responses varied widely, and no consistent patterns were observed among response metrics laterally or vertically along hillslopes. In contrast to other studies, we found that the relative influence of hillslope properties and storm characteristics varied with soil moisture responses and during storms. Antecedent conditions and storm depths influenced the strength of lag correlations between soil moisture and runoff, whereas storm mean intensity was correlated with the lag times. These results highlight the utility of intensive observations for characterizing heterogeneity in soil moisture responses, suggesting, among other things, a need for better representation of the subsurface processes in rainfall‐runoff models. Identifying the relative importance of drivers can be beneficial in building parsimonious hydrological models.
-
more » « less
Abstract Precipitation is the primary driver of hydrological models, and its spatial and temporal variability have a great impact on water partitioning. However, in data‐sparse regions, uncertainty in precipitation estimates is high and the sensitivity of water partitioning to this uncertainty is unknown. This is a particular challenge in drylands (semi‐arid and arid regions) where the water balance is highly sensitive to rainfall, yet there is commonly a lack of in situ rain gauge data. To understand the impact of precipitation uncertainty on the water balance in drylands, here we have performed simulations with a process‐based hydrological model developed to characterize the water balance in arid and semi‐arid regions (DRYP: DRYland water Partitioning model). We performed a series of numerical analyses in the Upper Ewaso Ng'iro basin, Kenya driven by three gridded precipitation datasets with different spatio‐temporal resolutions (IMERG, MSWEP, and ERA5), evaluating simulations against streamflow observations and remotely sensed data products of soil moisture, actual evapotranspiration, and total water storage. We found that despite the great differences in the spatial distribution of rainfall across a climatic gradient within the basin, DRYP shows good performance for representing streamflow (KGE >0.6), soil moisture, actual evapotranspiration, and total water storage (
r >0.5). However, the choice of precipitation datasets greatly influences surface (infiltration, runoff, and transmission losses) and subsurface fluxes (groundwater recharge and discharge) across different climatic zones of the Ewaso Ng'iro basin. Within humid areas, evapotranspiration does not show sensitivity to the choice of precipitation dataset, however, in dry lowland areas it becomes more sensitive to precipitation rates as water‐limited conditions develop. The analysis shows that the highest rates of precipitation produce high rates of diffuse recharge in Ewaso uplands and also propagate into runoff, transmission losses and, ultimately focused recharge, with the latter acting as the main mechanism of groundwater recharge in low dry areas. The results from this modelling exercise suggest that care must be taken in selecting forcing precipitation data to drive hydrological modelling efforts, especially in basins that span a climatic gradient. These results also suggest that more effort is required to reduce uncertainty between different precipitation datasets, which will in turn result in more consistent quantification of the water balance. -
more » « less
Abstract The size distributions of sediment delivered from hillslopes to rivers profoundly influence river morphodynamics, including river incision into bedrock and the quality of aquatic habitat. Yet little is known about the factors that influence size distributions of sediment produced by weathering on hillslopes. We present results of a field study of hillslope sediment size distributions at Inyo Creek, a steep catchment in granitic bedrock of the Sierra Nevada, USA. Particles sampled near the base of hillslopes, adjacent to the trunk stream, show a pronounced decrease in sediment size with decreasing sample elevation across all but the coarsest size classes. Measured size distributions become increasingly bimodal with decreasing elevation, exhibiting a coarse, bouldery mode that does not change with elevation and a more abundant finer mode that shifts from cobbles at the highest elevations to gravel at mid elevations and finally to sand at low elevations. We interpret these altitudinal variations in hillslope sediment size to reflect changes in physical, chemical, and biological weathering that can be explained by the catchment's strong altitudinal gradients in topography, climate, and vegetation cover. Because elevation and travel distance to the outlet are closely coupled, the altitudinal trends in sediment size produce a systematic decrease in sediment size along hillslopes parallel to the trunk stream. We refer to this phenomenon as ‘downvalley fining.’ Forward modeling shows that downvalley fining of hillslope sediment is necessary for downstream fining of the long‐term average flux of coarse sediment in mountain landscapes where hillslopes and channels are coupled and long‐term net sediment deposition is negligible. The model also shows that abrasion plays a secondary role in downstream fining of coarse sediment flux but plays a dominant role in partitioning between the bedload and suspended load. Patterns observed at Inyo Creek may be widespread in mountain ranges around the world. © 2020 John Wiley & Sons Ltd
