Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr^?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.     

Regional patterns of postwildfire streamflow response in the Western United States: The importance of scale‐specific connectivity

Wildfires can impact streamflow by modifying net precipitation, infiltration, evapotranspiration, snowmelt, and hillslope run‐off pathways. Regional differences in fire trends and postwildfire streamflow responses across the conterminous United States have spurred concerns about the impact on streamflow in forests that serve as water resource areas. This is notably the case for the Western United States, where fire activity and burn severity have increased in conjunction with climate change and increased forest density due to human fire suppression. In this review, we discuss the effects of wildfire on hydrological processes with a special focus on regional differences in postwildfire streamflow responses in forests. Postwildfire peak flows and annual water yields are generally higher in regions with a Mediterranean or semi‐arid climate (Southern California and the Southwest) compared to the highlands (Rocky Mountains and the Pacific Northwest), where fire‐induced changes in hydraulic connectivity along the hillslope results in the delivery of more water, more rapidly to streams. No clear streamflow response patterns have been identified in the humid subtropical Southeastern United States, where most fires are prescribed fires with a low burn severity, and more research is needed in that region. Improved assessment of postwildfire streamflow relies on quantitative spatial knowledge of landscape variables such as prestorm soil moisture, burn severity and correlations with soil surface sealing, water repellency, and ash deposition. The latest studies furthermore emphasize that understanding the effects of hydrological processes on postwildfire dynamic hydraulic connectivity, notably at the hillslope and watershed scales, and the relationship between overlapping disturbances including those other than wildfire is necessary for the development of risk assessment tools.     

Variable power‐law scaling of hillslope Hortonian rainfall–runoff processes

Hydrological studies focused on Hortonian rainfall–run‐off scaling have found that the run‐off depth generally declines with the plot length in power‐law scaling. Both the power‐law proportional coefficient and the scaling exponent show great variability for specific conditions, but why and how they vary remain unclear. In the present study, the scaling of hillslope Hortonian rainfall–run‐off processes is investigated for different rainfall, soil infiltration, and hillslope surface characteristics using the physically based cell‐based rainfall‐infiltration‐run‐off model. The results show that both temporally intermittent and steady rainfalls can result in prominent power‐law scaling at the initial stage of run‐off generation. Then, the magnitude of the power‐law scaling decreases gradually due to the decreasing run‐on effect. The power‐law scaling is most sensitive to the rainfall and soil infiltration parameters. When the ratio of rainfall to infiltration exceeds a critical value, the magnitude of the power‐law scaling tends to decrease notably. For different intermittent rainfall patterns, the power‐law exponent varies in the range of ?1.0 to ?0.113, which shows an approximately logarithmic increasing trend for the proportional coefficient as a function of the run‐off coefficient. The scaling is also sensitive to the surface roughness, soil sealing, slope angle, and hillslope geometry because these factors control the run‐off routing and run‐on infiltration processes. These results provide insights into the variable scaling of the Hortonian rainfall–run‐off process, which are expected to benefit modelling of large‐scale hydrological and ecological processes.     

The importance of surface controls on overland flow connectivity in semi‐arid environments: results from a numerical experimental approach

In semi‐arid environments, the characteristics of the land surface determine how rainfall is transformed into surface runoff and influences how this runoff moves from the hillslopes into river channels. Whether or not water reaches the river channel is determined by the hydrological connectivity. This paper uses a numerical experiment‐based approach to systematically assess the effects of slope length, gradient, flow path convergence, infiltration rates and vegetation patterns on the generation and connectivity of runoff. The experiments were performed with the Connectivity of Runoff Model, 2D version distributed, physically based, hydrological model. The experiments presented are set within a semi‐arid environment, characteristic of south‐eastern Spain, which is subject to low frequency high rainfall intensity storm events. As a result, the dominant hydrological processes are infiltration excess runoff generation and surface flow dynamics. The results from the modelling experiments demonstrate that three surface factors are important in determining the form of the discharge hydrograph: the slope length, the slope gradient and the infiltration characteristics at the hillslope‐channel connection. These factors are all related to the time required for generated runoff to reach an efficient flow channel, because once in this channel, the transmission losses significantly decrease. Because these factors are distributed across the landscape, they have a fundamental role in controlling the landscape hydrological response to storm events. Copyright © 2013 John Wiley & Sons, Ltd.     

Effects of forest roads on the hydrological response of a small‐scale mountain watershed in Greece

The effects of land use changes on the ecology and hydrology of natural watersheds have long been debated. However, less attention has been given to the hydrological effects of forest roads. Although less studied, several researchers have claimed that streamflow changes related to forest roads can cause a persistent and pervasive effect on hillslope hydrology and the functioning of the channel system. The main potential direct effects of forest roads on natural watersheds hydrologic response are runoff production on roads surfaces due to reduced infiltration rates, interruption of subsurface flow by road cutslopes and rapid transfer of the produced runoff to the stream network through roadside ditches. The aforementioned effects may significantly modify the total volume and timing of the hillslope flow to the stream network. This study uses detailed field data, spatial data, hydro‐meteorological records, as well as numerical simulation to investigate the effects of forest roads on the hydrological response of a small‐scale mountain experimental watershed, which is situated in the east side of Penteli Mountain, Attica, Greece. The results of this study highlight the possible effects of forest roads on the watersheds hydrological response that may significantly influence direct runoff depths and peak flow rates. It is demonstrated that these effects can be very important in permeable watersheds and that more emphasis should be given on the impact of roads on the watersheds hydrological response. Copyright © 2014 John Wiley & Sons, Ltd.     

Testing of the SIBERIA landscape evolution model using the Tin Camp Creek,Northern Territory,Australia, field catchment

The SIBERIA landscape evolution model was used to simulate the geomorphic development of the Tin Camp Creek natural catchment over geological time. Measured hydrology, erosion and geomorphic data were used to calibrate the SIBERIA model, which was then used to make independent predictions of the landform geomorphology of the study site. The catchment, located in the Northern Territory, Australia is relatively untouched by Europeans so the hydrological and erosion processes that shaped the area can be assumed to be the same today as they have been in the past, subject to the caveats regarding long‐term climate fluctuation. A qualitative, or visual comparison between the natural and simulated catchments indicates that SIBERIA can match hillslope length and hillslope profile of the natural catchments. A comparison of geomorphic and hydrological statistics such as the hypsometric curve, width function, cumulative area distribution and area–slope relationship indicates that SIBERIA can model the geomorphology of the selected Tin Camp Creek catchments. Copyright 2002 © Environmental Research Institute of the Supervising Scientist, Commonwealth of Australia.     

Modelling of sediment dynamics in a laboratory‐scale experimental catchment

The variability of hillslope form and function is examined experimentally using a simple model catchment in which most landscape development parameters are either known or controlled. It is demonstrated that there is considerable variability in sediment output from similar catchments, subjected to the same hydrological processes, and for which the initial hillslope profiles are the same. The results demonstrate that, in the case of catchments with a linear initial hillslope profile, the sediment output is initially high but reduces through time, whereas for a concave initial profile the sediment output was smaller and relatively constant. Concave hillslope profiles also displayed reduced sediment output when compared with linear slopes with the same overall slope. Using this experimental model catchment data, the SIBERIA landscape evolution model was tested for its ability to predict temporal sediment transport. When calibrated for the rainfall and erodible material, SIBERIA is able to simulate mean temporal sediment output for the experimental catchment over a range of hillslope profiles and rainfall intensities. SIBERIA is also able to match the hillslope profile of the experimental catchments. The results of the study provide confidence in the ability of SIBERIA to predict temporal sediment output. The experimental and modelling data also demonstrate that, even with all geomorphic and hydrological variables being known and/or controlled, there is still a need for long‐term stream gauging to obtain reliable assessments of field catchment hydrology and sediment transport. Copyright © 2005 John Wiley & Sons, Ltd.     

The design of post‐mining landscapes using geomorphic principles

Nature can provide analogues for post‐mining landscapes in terms of landscape stability and also in terms of the rehabilitated structure ‘blending in’ with the surrounding undisturbed landscape. In soil‐mantled landscapes, hillslopes typically have a characteristic pro?le that has a convex upper hillslope pro?le with a concave pro?le lower down the slope. In this paper hillslope characteristic form is derived using the area–slope relationship from pre‐mining topography at two sites in Western Australia. Using this relationship, concave hillslope pro?les are constructed and compared to linear hillslopes in terms of sediment loss using the SIBERIA erosion model. It is found that concave hillslopes can reduce sediment loss by up to ?ve times that of linear slopes. Concave slopes can therefore provide an alternative method for the construction of post‐mining landscapes. An understanding of landscape geomorphological properties and the use of erosion models can greatly assist in the design of post‐mining landscapes. Copyright © 2003 John Wiley & Sons, Ltd.     

Coupled numerical–analytical approach to landscape evolution modeling

The Earth's topography is shaped by surface processes that operate on various scales. In particular, river processes control landscape dynamics over large length scales, whereas hillslope processes control the dynamics over smaller length scales. This scale separation challenges numerical treatments of landscape evolution that use space discretization. Large grid spacing cannot account for the dynamics of water divides that control drainage area competition, and erosion rate and slope distribution. Small grid spacing that properly accounts for divide dynamics is computationally inefficient when studying large domains. Here we propose a new approach for landscape evolution modeling that couples irregular grid‐based numerical solutions for the large‐scale fluvial dynamics and continuum‐based analytical solutions for the small‐scale fluvial and hillslope dynamics. The new approach is implemented in the landscape evolution model DAC (divide and capture). The geometrical and topological characteristics of DAC's landscapes show compatibility with those of natural landscapes. A comparative study shows that, even with large grid spacing, DAC predictions fit well an analytical solution for divide migration in the presence of horizontal advection of topography. In addition, DAC is used to study some outstanding problems in landscape evolution. (i) The time to steady‐state is investigated and simulations show that steady‐state requires much more time to achieve than predicted by fixed area calculations, due to divides migration and persistent reorganization of low‐order streams. (ii) Large‐scale stream captures in a strike‐slip environment are studied and show a distinct pattern of erosion rates that can be used to identify recent capture events. (iii) Three tectono‐climatic mechanisms that can lead to asymmetric mountains are studied. Each of the mechanisms produces a distinct morphology and erosion rate distribution. Application to the Southern Alps of New Zealand suggests that tectonic advection, precipitation gradients and non‐uniform tectonic uplift act together to shape the first‐order topography of this mountain range. Copyright © 2014 John Wiley & Sons, Ltd.     

An analysis of the impact of spatial variability in rainfall on runoff and sediment predictions from a distributed model

This paper investigates the effect of introducing spatially varying rainfall fields to a hydrological model simulating runoff and erosion. Pairs of model simulations were run using either spatially uniform (i.e. spatially averaged) or spatially varying rainfall fields on a 500‐m grid. The hydrological model used was a simplified version of Thales which enabled runoff generation processes to be isolated from hillslope averaging processes. Both saturation excess and infiltration excess generation mechanisms were considered, as simplifications of actual hillslope processes. A 5‐year average recurrence interval synthetic rainfall event typical of temperate climates (Melbourne, Australia) was used. The erosion model was based on the WEPP interrill equation, modified to allow nonlinear terms relating the erosion rate to rainfall or runoff‐squared. The model results were extracted at different scales to investigate whether the effects of spatially varying rainfall were scale dependent. A series of statistical metrics were developed to assess the variability due to introducing the spatially varying rainfall field. At the catchment (approximately 150 km²) scale, it was found that particularly for saturation excess runoff, model predictions of runoff were insensitive to the spatial resolution of the rainfall data. Generally, erosion processes at smaller sub‐catchment scales, particularly when the sediment generation equation had non linearity, were more sensitive to spatial rainfall variability. Introducing runon infiltration reduced the total runoff and sediment yield at all scales, and this process was also most sensitive to the rainfall resolution. Copyright © 2011 John Wiley & Sons, Ltd.     

A model for distributed GIUH‐based flow routing on natural and anthropogenic hillslopes

Attempts to reduce the number of parameters in distributed rainfall–runoff models have not yet resulted in a model that is accurate for both natural and anthropogenic hillslopes. We take on the challenge by proposing a distributed model for overland flow and channel flow based on a combination of a linear response time distribution and the hillslope geomorphologic instantaneous unit hydrograph (GIUH), which can be calculated with only a digital elevation model and a map with field boundaries and channel network as input. The spatial domain is subdivided into representative elementary hillslopes (REHs) for each of which we define geometric and flow velocity parameters and compute the GIUH. The catchment GIUH is given by the sum of all REH responses. While most distributed models only perform well on natural hillslopes, the advantage of our approach is that it can also be applied to modified hillslopes with for example a rectangular drainage network and terrace cultivation. Tests show that the REH‐GIUH approach performs better than classical routing functions (exponential and gamma). Simulations of four virtual hillslopes suggest that peak flow at the catchment outlet is directly related to drainage density. By combining the distributed flow routing model with a lumped‐parameter infiltration model, we were also able to demonstrate that terrace cultivation delays the response time and reduces peak flow in comparison to the same hillslope, but with a natural stream network. The REH‐GIUH approach is a first step in the process of coupling distributed hydrological models to erosion and water quality models at the REH (associated with agricultural management) and at the catchment scale (associated with the evaluation of the environmental impact of human activities). It furthermore provides a basis for the development of models for large catchments and urban or peri‐urban catchments. Copyright © 2013 John Wiley & Sons, Ltd.     

Precipitation‐runoff processes in Shimen hillslope micro‐catchment of Taihang Mountain,north China

Precipitation runoff is a critical hillslope hydrological process for downslope streamflow and piedmont/floodplain recharge. Shimen hillslope micro‐catchment is strategically located in the central foothill region of Taihang Mountains, where runoff is crucial for water availability in the piedmont corridors and floodplains of north China. This study analyzes precipitation‐runoff processes in the Shimen hillslope micro‐catchment for 2006–2008 using locally designed runoff collection systems. The study shows that slope length is a critical factor, next only to precipitation, in terms of runoff yield. Regression analysis also shows that runoff is related positively to precipitation, and negatively to slope length. Soil mantle in the study area is generally thin and is therefore not as critical a runoff factor as slope length. The study shows a significant difference between overland and subsurface runoff. However, that between the 0–10 and 10–20 cm subsurfaces is insignificant. Runoff hardly occurs under light rains (<10 mm), but is clearly noticeable under moderate‐to‐rainstorm events. In the hillslope catchment, vertical infiltration (accounting for 42–84% of the precipitation) dominates runoff processes in subsurface soils and weathered granite gneiss bedrock. A weak lateral flow (at even the soil/bedrock interface) and the generally small runoff suggest strong infiltration loss via deep percolation. This is critical for groundwater recharge in the downslope piedmont corridors and floodplains. This may enhance water availability, ease water shortage, avert further environmental degradation, and reduce the risk of drought/flood in the event of extreme weather conditions in the catchment and the wider north China Plain. Copyright © 2011 John Wiley & Sons, Ltd.     

Capabilities and limitations of detailed hillslope hydrological modelling

Hillslope hydrological modelling is considered to be of great importance for the understanding and quantification of hydrological processes in hilly or mountainous landscapes. In recent years a few comprehensive hydrological models have been developed at the hillslope scale which have resulted in an advanced representation of hillslope hydrological processes (including their interactions), and in some operational applications, such as in runoff and erosion studies at the field scale or lateral flow simulation in environmental and geotechnical engineering. An overview of the objectives of hillslope hydrological modelling is given, followed by a brief introduction of an exemplary comprehensive hillslope model, which stimulates a series of hydrological processes such as interception, evapotranspiration, infiltration into the soil matrix and into macropores, lateral and vertical subsurface soil water flow both in the matrix and preferential flow paths, surface runoff and channel discharge. Several examples of this model are presented and discussed in order to determine the model's capabilities and limitations. Finally, conclusions about the limitations of detailed hillslope modelling are drawn and an outlook on the future prospects of hydrological models on the hillslope scale is given.The model presented performed reasonable calculations of Hortonian surface runoff and subsequent erosion processes, given detailed information of initial soil water content and soil hydraulic conditions. The vertical and lateral soil moisture dynamics were also represented quite well. However, the given examples of model applications show that quite detailed climatic and soil data are required to obtain satisfactory results. The limitations of detailed hillslope hydrological modelling arise from different points: difficulties in the representations of certain processes (e.g. surface crusting, unsaturated–saturated soil moisture flow, macropore flow), problems of small‐scale variability, a general scarcity of detailed soil data, incomplete process parametrization and problems with the interdependent linkage of several hillslopes and channel–hillslope interactions. Copyright © 1999 John Wiley & Sons, Ltd.     

The effect of urbanization on stream hydrology in hillslope watersheds in central Texas

This study examined the effect of urbanization on stream hydrology in hillslope watersheds. Ten streams (seven in hillslope and three in gentle slope watersheds) around Austin, Texas were selected for analysis. For each stream, we compared parameters of transfer function (TF) models estimated from daily rainfall and streamflow data collected in two study periods (October 1988–September 1992 and October 2004–September 2008) representing different degrees of watershed urbanization. As expected, the streams became more intermittent as the watersheds were more urbanized in all the study streams. However, the effect of urbanization on peakflow differs between hillslope and gentle slope watersheds. After watershed urbanization, peakflow increased in gentle slope watersheds, but decreased in hillslope watersheds. Based on the results of the TF models, we found that urbanization made stream not flashier but drier in hillslope watersheds. Overpumpage of aquifer has been recognized as a problem that leads to the stream dryness in the study area. However, the overpumpage alone cannot explain the differences in hydrological changes between the two types of watersheds. We attributed the reduced peakflow and stream dryness in the hillslope watersheds to land grading for construction forming stair‐stepped or terraced landscape. Compared with natural hillslope, a stair‐stepped landscape could infiltrate more stormwater by slowing down surface runoff on tread portions of the stair. Our findings suggest that a watershed management scheme should take into account local hydrogeologic conditions to mitigate the stream dryness resulting from urbanization in hillslope watersheds. Copyright © 2010 John Wiley & Sons, Ltd.     

Hydrology of hillslope‐wetland streams,Polar Bear Pass,Nunavut, Canada

Polar Bear Pass is a large High Arctic low‐gradient wetland (100 km²) bordered by low‐lying hills which are notched by a series of v‐shaped valleys. The spring and summer hydrology of two High Arctic hillslope‐wetland catchments, a first‐order stream, 0·2 km² Landing Strip Creek (LSC) and a larger second‐order basin, 4·2 km² Windy Creek (WC), is described here. A water balance framework was employed in 2008 to examine the movement of water from upland reaches into the low‐lying wetland. Snowcover was low in both basins (<50 mm in water equivalent units), but they both exhibited nival‐type regimes. After the main snowmelt season ended, runoff ceased in the smaller catchment (LSC), but not at the larger basin (WC) which continued to flow throughout the summer. Both basins responded to summer rains in different ways. At LSC, late‐summer continuous streamflow occurred only when rainfall satisfied the large soil moisture deficit in the upper bowl‐shaped zone of the basin. At WC, the presence of thinly thawed, ice‐rich polygonal terrain within the stream channel and in the upper reaches of the catchment likely limited infiltration in these near‐stream zones and enhanced runoff in response to both moderate and high rainfall. Subsequently, seasonal runoff ratios differed between the two sites (0·19 vs 0·68) as did the seasonal storage + residual (+16 vs ?50 mm). This suggests that the post‐snowmelt season runoff response to summer precipitation is very much modified by the unique basin characteristics (soil‐type, vegetation, ground ice) and their location within each stream order type. Copyright © 2010 John Wiley & Sons, Ltd.     

Gully,channel and hillslope erosion – an assessment for a traditionally managed catchment

Understanding landscape features such as gullying and soil erosion is an important issue in the long‐term dynamics and evolution of both natural, agricultural and rehabilitated (i.e. post‐mining) landscapes. Considerable research has been undertaken examining the initiation, movement and overall dynamics of such features. This study reports on a series of 34 gully heads and other erosion features, such as scour holes (five in total), located in channels in a catchment largely undisturbed by European activity in the Northern Territory, Australia over a 5 year period (2002–2007). During this period the erosion features were monitored for their headward advance/retreat, enlargement or in‐filling. The erosion features ranged in depth from 0.2 m to 1.5 m and widths of 0.3 m to 8 m. Hillslope erosion was also monitored using erosion pins. The catchment was subject to a range of rainfall regimes including extreme rainfall and a Category 5 cyclone and also was burnt every second year so that all grass cover was removed according to traditional management practice. The results of this monitoring show that the erosion features have changed little during this 5 year period. A remote sensing assessment found no relationship between erosion feature morphology and hillslope erosion. The monitored gullies heads and scour holes appear to be resilient landscape features, yet have a morphology that suggests they are ready for rapid headward movement and expansion, leading to a destabilisation of the catchment. Hillslope erosion was found to be related to wetness indices derived from a digital elevation model. Significant linkages were found between hillslope erosion and change in erosion feature depth, indicative of a strong hillslope–channel coupling. Copyright © 2010 Commonwealth of Australia     

Improvement of surface run‐off in the hydrological model ParFlow by a scale‐consistent river parameterization

We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels.     

Tracing variability of run‐off generation in mountainous permafrost of semi‐arid north‐eastern Mongolia

The headwaters of mountainous, discontinuous permafrost regions in north‐eastern Mongolia are important water resources for the semi‐arid country, but little is known about hydrological processes there. Run‐off generation on south‐facing slopes, which are devoid of permafrost, has so far been neglected and is totally unknown for areas that have been affected by recent forest fires. To fill this knowledge gap, the present study applied artificial tracers on a steppe‐vegetated south‐facing and on two north‐facing slopes, burned and unburned. Combined sprinkling and dye tracer experiments were used to visualize processes of infiltration and water fluxes in the unsaturated zone. On the unburned north‐facing slope, rapid and widespread infiltration through a wet organic layer was observed down to the permafrost. On the burned profile, rapid infiltration occurred through a combusted organic and underlying mineral layer. Stained water seeped out at the bottom of both profiles suggesting a general tendency to subsurface stormflow (SSF). Ongoing SSF could directly be studied 24 h after a high‐intensity rainfall event on a 55‐m hillslope section in the burned forest. Measurements of water temperature proved the role of the permafrost layer as a base horizon for SSF. Repeated tracer injections allowed direct insights into SSF dynamics: A first injection suggested rather slow dispersive subsurface flow paths; whereas 18 h later, a second injection traced a more preferential flow system with 20 times quicker flow velocities. We speculate that these pronounced SSF dynamics are limited to burned slopes where a thermally insulating organic layer is absent. On three south‐facing soil profiles, the applied tracer remained in the uppermost 5 cm of a silt‐rich mineral soil horizon. No signs of preferential infiltration could be found, which suggested reduced biological activity under a harsh, dry and cold climate. Instead, direct observations, distributed tracers and charcoal samples provided evidence for the occurrence of overland flow. Copyright © 2014 John Wiley & Sons, Ltd.     

Relations between rainfall–runoff‐induced erosion and aeolian deposition at archaeological sites in a semi‐arid dam‐controlled river corridor

Process dynamics in fluvial‐based dryland environments are highly complex with fluvial, aeolian, and alluvial processes all contributing to landscape change. When anthropogenic activities such as dam‐building affect fluvial processes, the complexity in local response can be further increased by flood‐ and sediment‐limiting flows. Understanding these complexities is key to predicting landscape behavior in drylands and has important scientific and management implications, including for studies related to paleoclimatology, landscape ecology evolution, and archaeological site context and preservation. Here we use multi‐temporal LiDAR surveys, local weather data, and geomorphological observations to identify trends in site change throughout the 446‐km‐long semi‐arid Colorado River corridor in Grand Canyon, Arizona, USA, where archaeological site degradation related to the effects of upstream dam operation is a concern. Using several site case studies, we show the range of landscape responses that might be expected from concomitant occurrence of dam‐controlled fluvial sand bar deposition, aeolian sand transport, and rainfall‐induced erosion. Empirical rainfall‐erosion threshold analyses coupled with a numerical rainfall–runoff–soil erosion model indicate that infiltration‐excess overland flow and gullying govern large‐scale (centimeter‐ to decimeter‐scale) landscape changes, but that aeolian deposition can in some cases mitigate gully erosion. Whereas threshold analyses identify the normalized rainfall intensity (defined as the ratio of rainfall intensity to hydraulic conductivity) as the primary factor governing hydrologic‐driven erosion, assessment of false positives and false negatives in the dataset highlight topographic slope as the next most important parameter governing site response. Analysis of 4+ years of high resolution (four‐minute) weather data and 75+ years of low resolution (daily) climate records indicates that dryland erosion is dependent on short‐term, storm‐driven rainfall intensity rather than cumulative rainfall, and that erosion can occur outside of wet seasons and even wet years. These results can apply to other similar semi‐arid landscapes where process complexity may not be fully understood. Published 2015. This article is a U.S. Government work and is in the public domain in the USA