Tracking water pathways in steep hillslopes by d18O depth profiles of soil water

Assessing temporal variations in soil water flow is important, especially at the hillslope scale, to identify mechanisms of runoff and flood generation and pathways for nutrients and pollutants in soils. While surface processes are well considered and parameterized, the assessment of subsurface processes at the hillslope scale is still challenging since measurement of hydrological pathways is connected to high efforts in time, money and personnel work. The latter might not even be possible in alpine environments with harsh winter processes. Soil water stable isotope profiles may offer a time-integrating fingerprint of subsurface water pathways. In this study, we investigated the suitability of soil water stable isotope (d18O) depth profiles to identify water flow paths along two transects of steep subalpine hillslopes in the Swiss Alps. We applied a one-dimensional advection–dispersion model using d18O values of precipitation (ranging from _24.7 to _2.9‰) as input data to simulate the d18O profiles of soil water. The variability of d18O values with depth within each soil profile and a comparison of the simulated and measured d18O profiles were used to infer information about subsurface hydrological pathways. The temporal pattern of d18O in precipitation was found in several profiles, ranging from _14.5 to _4.0‰. This suggests that vertical percolation plays an important role even at slope angles of up to 46_. Lateral subsurface flow and/or mixing of soil water at lower slope angles might occur in deeper soil layers and at sites near a small stream. The difference between several observed and simulated d18O profiles revealed spatially highly variable infiltration patterns during the snowmelt periods: The d18O value of snow (_17.7 ± 1.9‰) was absent in several measured d18O profiles but present in the respective simulated d18O profiles. This indicated overland flow and/or preferential flow through the soil profile during the melt period. The applied methods proved to be a fast and promising tool to obtain time-integrated information on soil water flow paths at the hillslope scale in steep subalpine slopes.

Debris flow magnitudes: a global catalog

The episodic occurrence of debris flow events in response to stochastic precipitation and wildfire events makes hazard prediction challenging. Previous work has shown that frequency-magnitude distributions of non-fire-related debris flows follow a power law, but less is known about the distribution of post-fire debris flows. As a first step in parameterizing hazard models, we use frequency-magnitude distributions and cumulative distribution functions to compare volumes of post-fire debris flows to non-fire-related debris flows. Due to the large number of events required to parameterize frequency-magnitude distributions, and the relatively small number of post-fire event magnitudes recorded in the literature, we collected data on 73 recent post-fire events in the field. The resulting catalog of 988 debris flow events is presented as an appendix to this article. We found that the empirical cumulative distribution function of post-fire debris flow volumes is composed of smaller events than that of non-fire-related debris flows. In addition, the slope of the frequency-magnitude distribution of post-fire debris flows is steeper than that of non-fire-related debris flows, evidence that differences in the post-fire environment tend to produce a higher proportion of small events. We propose two possible explanations: 1) post-fire events occur on shorter return intervals than debris flows in similar basins that do not experience fire, causing their distribution to shift toward smaller events due to limitations in sediment supply, or 2) fire causes changes in resisting and driving forces on a package of sediment, such that a smaller perturbation of the system is required in order for a debris flow to occur, resulting in smaller event volumes.

In situ investigation of rapid subsurface flow : Temporal dynamics and catchment-scale implication

Non peer reviewed

Locating the Seattle Fault in Bellevue, WA: Combining geomorphic indicators with borehole data

The Seattle Fault is an active east-west trending reverse fault zone that intersects both Seattle and Bellevue, two highly populated cities in Washington. Rupture along strands of the fault poses a serious threat to infrastructure and thousands of people in the region. Precise locations of fault strands are still poorly constrained in Bellevue due to blind thrusting, urban development, and/or erosion. Seismic reflection and aeromagnetic surveys have shed light on structural geometries of the fault zone in bedrock. However, the fault displaces both bedrock and unconsolidated Quaternary deposits, and seismic data are poor indicators of the locations of fault strands within the unconsolidated strata. Fortunately, evidence of past fault strand ruptures may also be recorded indirectly by fluvial processes and should also be observable in the subsurface. I analyzed hillslope and river geomorphology using LiDAR data and ArcGIS to locate surface fault traces and then compare/correlate these findings to subsurface offsets identified using borehole data. Geotechnical borings were used to locate one fault offset and provide input to a cross section of the fault constructed using Rockworks software. Knickpoints, which may correlate to fault rupture, were found upstream of this newly identified fault offset as well as upstream of a previously known fault segment.

Geology and landslide geomorphology of the Burpee Hills, Skagit County, Washington, USA

Landforms within the Skagit Valley record a complex history of land evolution from Late Pleistocene to the present. Late Pleistocene glacial deposits and subsequent incision by the Skagit River formed the Burpee Hills terrace. The Burpee Hills comprises an approximately 205-m-thick sequence of sediments, including glacio-lacustrine silts and clays, overlain by sandy advance outwash and capped by coarse till, creating a sediment-mantled landscape where mass wasting occurs in the form of debris flows and deep-seated landslides (Heller, 1980; Skagit County, 2014). Landslide probability and location are necessary metrics for informing citizens and policy makers of the frequency of natural hazards. Remote geomorphometric analysis of the site area using airborne LiDAR combined with field investigation provide the information to determine relative ages of landslide deposits, to classify geologic units involved, and to interpret the recent hillslope evolution. Thirty-two percent of the 28-km2 Burpee Hills landform has been mapped as landslide deposits. Eighty-five percent of the south-facing slope is mapped as landslide deposits. The mapped landslides occur predominantly within the advance outwash deposits (Qgav), this glacial unit has a slope angle ranging from 27 to 36 degrees. Quantifying surface roughness as a function of standard deviation of slope provides a relative age of landslide deposits, laying the groundwork for frequency analysis of landslides on the slopes of the Burpee Hills. The south-facing slopes are predominately affected by deep-seated landslides as a result of Skagit River erosion patterns within the floodplain. The slopes eroded at the toe by the Skagit River have the highest roughness coefficients, suggesting that areas with more frequent disturbance at the toe are more prone to sliding or remobilization. Future work including radiocarbon dating and hydrologic-cycle investigations will provide a more accurate timeline of the Burpee Hills hillslope evolution, and better information for emergency management and planners in the future.

Drying front in a sloping aquifer: Nonlinear effects

[1] The profiles for the water table height h(x, t) in a shallow sloping aquifer are reexamined with a solution of the nonlinear Boussinesq equation. We demonstrate that the previous anomaly first reported by Brutsaert [1994] that the point at which the water table h first becomes zero at x = L at time t = t(c) remains fixed at this point for all times t > t(c) is actually a result of the linearization of the Boussinesq equation and not, as previously suggested [Brutsaert, 1994; Verhoest and Troch, 2000], a result of the Dupuit assumption. Rather, by examination of the nonlinear Boussinesq equation the drying front, i.e., the point x(f) at which h is zero for times t greater than or equal to t(c), actually recedes downslope as physically expected. This points out that the linear Boussinesq equation should be used carefully when a zero depth is obtained as the concept of an average'' depth loses meaning at that time.

Similarity solution of axisymmetric flow in porous media

Applications of the axisymmetric Boussinesq equation to groundwater hydrology and reservoir engineering have long been recognised. An archetypal example is invasion by drilling fluid into a permeable bed where there is initially no such fluid present, a circumstance of some importance in the oil industry. It is well known that the governing Boussinesq model can be reduced to a nonlinear ordinary differential equation using a similarity variable, a transformation that is valid for a certain time-dependent flux at the origin. Here, a new analytical approximation is obtained for this case. The new solution,, which has a simple form, is demonstrated to be highly accurate. (c) 2005 Elsevier Ltd. All rights reserved.

A theoretical exploration of catchment-scale sediment delivery

This paper proposes a theoretical explanation of the variations of the sediment delivery ratio (SDR) versus catchment area relationships and the complex patterns in the behavior of sediment transfer processes at catchment scale. Taking into account the effects of erosion source types, deposition, and hydrological controls, we propose a simple conceptual model that consists of two linear stores arranged in series: a hillslope store that addresses transport to the nearest streams and a channel store that addresses sediment routing in the channel network. The model identifies four dimensionless scaling factors, which enable us to analyze a variety of effects on SDR estimation, including (1) interacting processes of erosion sources and deposition, (2) different temporal averaging windows, and (3) catchment runoff response. We show that the interactions between storm duration and hillslope/channel travel times are the major controls of peak-value-based sediment delivery and its spatial variations. The interplay between depositional timescales and the travel/residence times determines the spatial variations of total-volume-based SDR. In practical terms this parsimonious, minimal complexity model could provide a sound physical basis for diagnosing catchment to catchment variability of sediment transport if the proposed scaling factors can be quantified using climatic and catchment properties.

Modelling sediment delivery ratio over the Murray Darling Basin

This paper presents a scientific and technical description of the modelling framework and the main results of modelling the long-term average sediment delivery at hillslope to medium-scale catchments over the entire Murray Darling Basin (MDB). A theoretical development that relates long-term averaged sediment delivery to the statistics of rainfall and catchment parameters is presented. The derived flood frequency approach was adapted to investigate the problem of regionalization of the sediment delivery ratio (SDR) across the Basin. SDR, a measure of catchment response to the upland erosion rate, was modeled by two lumped linear stores arranged in series: hillslope transport to the nearest streams and flow routing in the channel network. The theory shows that the ratio of catchment sediment residence time (SRT) to average effective rainfall duration is the most important control in the sediment delivery processes. In this study, catchment SRTs were estimated using travel time for overland flow multiplied by an enlargement factor which is a function of particle size. Rainfall intensity and effective duration statistics were regionalized by using long-term measurements from 195 pluviograph sites within and around the Basin. Finally, the model was implemented across the MDB by using spatially distributed soil, vegetation, topographical and land use properties under Geographic Information System (GIs) environment. The results predict strong variations in SDR from close to 0 in floodplains to 70% in the eastern uplands of the Basin. (c) 2005 Elsevier Ltd. All rights reserved.

Spatial organisation of groundwater dynamics and streamflow response from different hydropedological units in a montane catchment

Funding was provided by the Leibniz Association (SAW-2012-IGB 4167) within the International Leibniz Graduate School: Aquatic boundaries and linkages- Aqualink. We would like to thank the NRI staff for their help during field work.

The concept of transport capacity and its use in different geomorphic process domains

The notion of sediment-transport capacity has been engrained in geomorphological and related literature for over 50 years, although its earliest roots date back explicitly to Gilbert in fluvial geomorphology in the 1870s and implicitly to eighteenth to nineteenth century developments in engineering. Despite cross fertilization between different process domains, there seem to have been independent inventions of the idea in aeolian geomorphology by Bagnold in the 1930s and in hillslope studies by Ellison in the 1940s. Here we review the invention and development of the idea of transport capacity in the fluvial, aeolian, coastal, hillslope, débris flow, and glacial process domains. As these various developments have occurred, different definitions have been used, which makes it both a difficult concept to test, and one that may lead to poor communications between those working in different domains of geomorphology. We argue that the original relation between the power of a flow and its ability to transport sediment can be challenged for three reasons. First, as sediment becomes entrained in a flow, the nature of the flow changes and so it is unreasonable to link the capacity of the water or wind only to the ability of the fluid to move sediment. Secondly, environmental sediment transport is complicated, and the range of processes involved in most movements means that simple relationships are unlikely to hold, not least because the movement of sediment often changes the substrate, which in turn affects the flow conditions. Thirdly, the inherently stochastic nature of sediment transport means that any capacity relationships do not scale either in time or in space. Consequently, new theories of sediment transport are needed to improve understanding and prediction and to guide measurement and management of all geomorphic systems.

Modelling the effectiveness of grass buffer strips in managing muddy floods under a changing climate

Muddy floods occur when rainfall generates runoff on agricultural land, detaching and transporting sediment into the surrounding natural and built environment. In the Belgian Loess Belt, muddy floods occur regularly and lead to considerable economic costs associated with damage to property and infrastructure. Mitigation measures designed to manage the problem have been tested in a pilot area within Flanders and were found to be cost-effective within three years. This study assesses whether these mitigation measures will remain effective under a changing climate. To test this, the Water Erosion Prediction Project (WEPP) model was used to examine muddy flooding diagnostics (precipitation, runoff, soil loss and sediment yield) for a case study hillslope in Flanders where grass buffer strips are currently used as a mitigation measure. The model was run for present day conditions and then under 33 future site-specific climate scenarios. These future scenarios were generated from three earth system models driven by four representative concentration pathways and downscaled using quantile mapping and the weather generator CLIGEN. Results reveal that under the majority of future scenarios, muddy flooding diagnostics are projected to increase, mostly as a consequence of large scale precipitation events rather than mean changes. The magnitude of muddy flood events for a given return period is also generally projected to increase. These findings indicate that present day mitigation measures may have a reduced capacity to manage muddy flooding given the changes imposed by a warming climate with an enhanced hydrological cycle. Revisions to the design of existing mitigation measures within existing policy frameworks are considered the most effective way to account for the impacts of climate change in future mitigation planning.