A spatially distributed model for the dynamic prediction of sediment erosion and transport in mountainous forested watersheds

Erosion and sediment transport in a temperate forested watershed are predicted with a new sediment model that represents the main sources of sediment generation in forested environments (mass wasting, hillslope erosion, and road surface erosion) within the distributed hydrology-soil-vegetation model (DHSVM) environment. The model produces slope failures on the basis of a factor-of-safety analysis with the infinite slope model through use of stochastically generated soil and vegetation parameters. Failed material is routed downslope with a rule-based scheme that determines sediment delivery to streams. Sediment from hillslopes and road surfaces is also transported to the channel network. A simple channel routing scheme is implemented to predict basin sediment yield. We demonstrate through an initial application of this model to the Rainy Creek catchment, a tributary of the Wenatchee River, which drains the east slopes of the Cascade Mountains, that the model produces plausible sediment yield and ratios of landsliding and surface erosion when compared to published rates for similar catchments in the Pacific Northwest. A road removal scenario and a basin-wide fire scenario are both evaluated with the model

The debris flow occurred at ru secco creek, venetian dolomites, on 4 august 2015: Analysis of the phenomenon, its characteristics and reproduction by models

On 4 August 2015, a very high intensity storm, 31.5 mm in 20 min (94.5 mm/h), hit the massif of Mount Antelao on the Venetian Dolomites triggering three stony debris \ufb02ows characterized by high magnitude. Two of them occurred in the historical sites of Rovina di Cancia and Rudan Creek and were stopped by the retaining works upstream the inhabited areas, while the third routed along the Ru Secco Creek and progressively reached the resort area and the village of San Vito di Cadore, causing fatalities and damages. The main triggering factor of the Ru Secco debris \ufb02ow was a large rock collapse on the northern cliffs of Mount Antelao occurred the previous autumn. The fallen debris material deposited on the Vallon d\u2019Antrimoia inclined plateau at the base of the collapsed cliffs and, below it, on the Ru Salvela Creek, covering it from the head to the con\ufb02uence with the Ru Secco Creek. The abundant runoff, caused by the high intensity rainfall on 4 August 2015, entrained about 52,500 m3 of the debris material laying on the Vallon d\u2019Antrimoia forming a debris \ufb02ow surge that hit and eroded the debris deposit covering the downstream Ru Salvela Creek, increasing its volume, about 110,000 m3 of mobilized sediments. This debris \ufb02ow routed downstream the con\ufb02uence, \ufb02ooding the parking of a resort area where three people died, and reached the village downstream damaging some buildings. A geomorphological analysis was initially carried out after surveying the whole basin. All liquid and solid-liquid contributions to the phenomenon were recognized together with the areas subjected to erosion and deposition. The elaboration of pre and post-event topographical surveys provided the map of deposition-erosion depths. Using the rainfall estimated by weather radar and corrected by the nearest rain gauge, about 0.8 km far, we estimated runoff by using a rainfall-runoff model designed for the headwater rocky basins of Dolomites. A triggering model provided the debris \ufb02ow hydrographs in the initiation areas, after using the simulated runoff. The initial solid-liquid surge hydrographs were, then, routed downstream by means of a cell model. The comparison between the simulated and estimated deposition-erosion pattern resulted satisfactory. The results of the simulation captured, in fact, the main features of the occurred phenomenon

A new data assimilation procedure to develop a debris flow run-out model

Abstract Parameter calibration is one of the most problematic phases of numerical modeling since the choice of parameters affects the model\u2019s reliability as far as the physical problems being studied are concerned. In some cases, laboratory tests or physical models evaluating model parameters cannot be completed and other strategies must be adopted; numerical models reproducing debris flow propagation are one of these. Since scale problems affect the reproduction of real debris flows in the laboratory or specific tests used to determine rheological parameters, calibration is usually carried out by comparing in a subjective way only a few parameters, such as the heights of soil deposits calculated for some sections of the debris flows or the distance traveled by the debris flows using the values detected in situ after an event has occurred. Since no automatic or objective procedure has as yet been produced, this paper presents a numerical procedure based on the application of a statistical algorithm, which makes it possible to define, without ambiguities, the best parameter set. The procedure has been applied to a study case for which digital elevation models of both before and after an important event exist, implicating that a good database for applying the method was available. Its application has uncovered insights to better understand debris flows and related phenomena

Simulating the influences of groundwater on regional geomorphology using a distributed, dynamic, landscape evolution modelling platform

A dynamic landscape evolution modelling platform (CLiDE) is presented that allows a variety of Earth system interactions to be explored under differing environmental forcing factors. Representation of distributed surface and subsurface hydrology within CLiDE is suited to simulation at sub-annual to centennial time-scales. In this study the hydrological components of CLiDE are evaluated against analytical solutions and recorded datasets. The impact of differing groundwater regimes on sediment discharge is examined for a simple, idealised catchment, Sediment discharge is found to be a function of the evolving catchment morphology. Application of CLiDE to the upper Eden Valley catchment, UK, suggests the addition of baseflow-return from groundwater into the fluvial system modifies the total catchment sediment discharge and the spatio-temporal distribution of sediment fluxes during storm events. The occurrence of a storm following a period of appreciable antecedent rainfall is found to increase simulated sediment fluxes

Analysis of Flood Patterns in Adams County, Pennsylvania Utilizing Drone Technology and Computer Simulations

Drone imagery and photogrammetry models of the Gettysburg College campus and the terrain at Boyer Nurseries and Orchards were utilized to study flood patterns in Adams County, Pennsylvania. Gettysburg College has lower-sloped land and moderately built infrastructure while Boyer Orchards has drastically sloped land with many patches of abundant vegetation. The two locations were selected due to the fact that they have starkly different surface features, while the bedrock geology of the areas are very similar. The terrain of the models was isolated before a 3D carver and 3D printer were used to construct physical models to further analyze potential water flow and speed through virtual, modeled flood simulations. The models were used to compare real world rainfall data and flood events in the investigated areas from the months of June to August in 2018. I hypothesized that the Gettysburg College campus would experience more severe flooding that would take longer to subside in comparison to Boyer Orchards due to the steeper slope of the orchards’ terrain. The research revealed that Boyer Orchards experienced more extreme flooding and rainfall than Gettysburg College but was able to neutralize the effects due to plentiful vegetation and physio-graphic differences. Modeled flood simulations demonstrated less rainfall in comparison to actual rainfall values: there were differences of 0.78 cm and 1.32 cm between the actual and simulated rainfall amounts for Gettysburg and the Boyer Orchards area, respectively

Flood magnitude-frequency and lithologic control on bedrock river incision in post-orogenic terrain

Mixed bedrock-alluvial rivers - bedrock channels lined with a discontinuous alluvial cover - are key agents in the shaping of mountain belt topography by bedrock fluvial incision. Whereas much research focuses upon the erosional dynamics of such rivers in the context of rapidly uplifting orogenic landscapes, the present study investigates river incision processes in a post-orogenic (cratonic) landscape undergoing extremely low rates of incision (> 5 m/Ma). River incision processes are examined as a function of substrate lithology and the magnitude and frequency of formative flows along Sandy Creek gorge, a mixed bedrock-alluvial stream in arid SE-central Australia. Incision is focused along a bedrock channel with a partial alluvial cover arranged into riffle-pool macrobedforms that reflect interactions between rock structure and large-flood hydraulics. Variations in channel width and gradient determine longitudinal trends in mean shear stress (τb) and therefore also patterns of sediment transport and deposition. A steep and narrow, non-propagating knickzone (with 5% alluvial cover) coincides with a resistant quartzite unit that subdivides the gorge into three reaches according to different rock erodibility and channel morphology. The three reaches also separate distinct erosional styles: bedrock plucking (i.e. detachment-limited erosion) prevails along the knickzone, whereas along the upper and lower gorge rock incision is dependent upon large formative floods exceeding critical erosion thresholds (τc) for coarse boulder deposits that line 70% of the channel thalweg (i.e. transport-limited erosion). The mobility of coarse bed materials (up to 2 m diameter) during late Holocene palaeofloods of known magnitude and age is evaluated using step-backwater flow modelling in conjunction with two selective entrainment equations. A new approach for quantifying the formative flood magnitude in mixed bedrock-alluvial rivers is described here based on the mobility of a key coarse fraction of the bed materials; in this case the d84 size fraction. A 350 m3/s formative flood fully mobilises the coarse alluvial cover with τb200-300 N/m2 across the upper and lower gorge riffles, peaking over 500 N/m2 in the knickzone. Such floods have an annual exceedance probability much less than 10- 2 and possibly as low as 10- 3. The role of coarse alluvial cover in the gorge is discussed at two scales: (1) modulation of bedrock exposure at the reach-scale, coupled with adjustment to channel width and gradient, accommodates uniform incision across rocks of different erodibility in steady-state fashion; and (2) at the sub-reach scale where coarse boulder deposits (corresponding to <i>τ</i><sub>b</sub> minima) cap topographic convexities in the rock floor, thereby restricting bedrock incision to rare large floods. While recent studies postulate that decreasing uplift rates during post-orogenic topographic decay might drive a shift to transport-limited conditions in river networks, observations here and elsewhere in post-orogenic settings suggest, to the contrary, that extremely low erosion rates are maintained with substantial bedrock channel exposure. Although bed material mobility is known to be rate-limiting for bedrock river incision under low sediment flux conditions, exactly how a partial alluvial cover might be spatially distributed to either optimise or impede the rate of bedrock incision is open to speculation. Observations here suggest that the small volume of very stable bed materials lining Sandy Creek gorge is distributed so as to minimise the rate of bedrock fluvial incision over time

Modeling river delta formation

A new model to simulate the time evolution of river delta formation process is presented. It is based on the continuity equation for water and sediment flow and a phenomenological sedimentation/ erosion law. Different delta types are reproduced using different parameters and erosion rules. The structures of the calculated patterns are analyzed in space and time and compared with real data patterns. Furthermore our model is capable to simulate the rich dynamics related to the switching of the mouth of the river delta. The simulation results are then compared with geological records for the Mississippi river

Large Scale Terrain Generation from Tectonic Uplift and Fluvial Erosion

International audienceAt large scale, landscapes result from the combination of two major processes: tectonics which generate the main relief through crust uplift, and weather which accounts for erosion. This paper presents the first method in computer graphics that combines uplift and hydraulic erosion to generate visually plausible terrains. Given a user-painted uplift map, we generate a stream graph over the entire domain embedding elevation information and stream flow. Our approach relies on the stream power equation introduced in geology for hydraulic erosion. By combining crust uplift and stream power erosion we generate large realistic terrains at a low computational cost. Finally, we convert this graph into a digital elevation model by blending landform feature kernels whose parameters are derived from the information in the graph. Our method gives high-level control over the large scale dendritic structures of the resulting river networks, watersheds, and mountains ridges

Evaluating the Differences of Gridding Techniques for Digital Elevation Models Generation and Their Influence on the Modeling of Stony Debris Flows Routing: A Case Study From Rovina di Cancia Basin (North-Eastern Italian Alps)

Debris \ufb02ows are among the most hazardous phenomena in mountain areas. To cope with debris \ufb02ow hazard, it is common to delineate the risk-prone areas through routing models. The most important input to debris \ufb02ow routing models are the topographic data, usually in the form of Digital Elevation Models (DEMs). The quality of DEMs depends on the accuracy, density, and spatial distribution of the sampled points; on the characteristics of the surface; and on the applied gridding methodology. Therefore, the choice of the interpolation method affects the realistic representation of the channel and fan morphology, and thus potentially the debris \ufb02ow routing modeling outcomes. In this paper, we initially investigate the performance of common interpolation methods (i.e., linear triangulation, natural neighbor, nearest neighbor, Inverse Distance to a Power, ANUDEM, Radial Basis Functions, and ordinary kriging) in building DEMs with the complex topography of a debris \ufb02ow channel located in the Venetian Dolomites (North-eastern Italian Alps), by using small footprint full- waveform Light Detection And Ranging (LiDAR) data. The investigation is carried out through a combination of statistical analysis of vertical accuracy, algorithm robustness, and spatial clustering of vertical errors, and multi-criteria shape reliability assessment. After that, we examine the in\ufb02uence of the tested interpolation algorithms on the performance of a Geographic Information System (GIS)-based cell model for simulating stony debris \ufb02ows routing. In detail, we investigate both the correlation between the DEMs heights uncertainty resulting from the gridding procedure and that on the corresponding simulated erosion/deposition depths, both the effect of interpolation algorithms on simulated areas, erosion and deposition volumes, solid-liquid discharges, and channel morphology after the event. The comparison among the tested interpolation methods highlights that the ANUDEM and ordinary kriging algorithms are not suitable for building DEMs with complex topography. Conversely, the linear triangulation, the natural neighbor algorithm, and the thin-plate spline plus tension and completely regularized spline functions ensure the best trade-off among accuracy and shape reliability. Anyway, the evaluation of the effects of gridding techniques on debris \ufb02ow routing modeling reveals that the choice of the interpolation algorithm does not signi\ufb01cantly affect the model outcomes