GM3.4 | Advances in modelling of erosion processes, sediment dynamics, and landscape evolution
EDI PICO
Advances in modelling of erosion processes, sediment dynamics, and landscape evolution
Convener: Jingtao Lai | Co-conveners: Kimberly Huppert, Boris Gailleton, Fiona Clubb, Jörg Robl
PICO
| Tue, 25 Apr, 08:30–10:15 (CEST)
 
PICO spot 3a
Tue, 08:30
A key goal within geomorphic research is understanding the processes linking topographic form, erosion rates, sediment production, transport and deposition, and external forcings such as tectonics, biotic or climatic. Numerical modelling, by allowing the creation of controlled analogues of natural systems, provides exciting opportunities to explore landscape evolution and generate testable predictions.

In this session, we invite contributions that use numerical modelling to investigate landscape evolution in a broad sense, and over a range of spatial and temporal scales. We welcome studies using models to constrain one or more of: erosion rates and processes, sediment production, transport and deposition, and biotic, climatic or tectonic forcings. We also particularly wish to highlight studies that combine numerical modelling with direct Earth surface process monitoring techniques, such as topographic, field, stratigraphic, geophysical or geochronological data. Contributions using numerical models to unravel the interaction between deep processes, such as mantle dynamics, or biotic processes with topographic patterns are further encouraged. There is no geographical restriction: studies may be focused on mountain environments or sedimentary basins, or they may establish links between the two. Studies beyond planet Earth are welcome too.

PICO: Tue, 25 Apr | PICO spot 3a

Chairpersons: Jingtao Lai, Kimberly Huppert
08:30–08:35
08:35–08:37
|
PICO3a.1
|
EGU23-2864
|
On-site presentation
Michael Li, Yongsheng Wu, Yongxing Ma, and Yolanda Wang

Ocean surface waves and currents can interact to produce strong seabed shear stress and mobilization of sediments that can significantly impact the seabed stability and benthic habitats on continental shelves. Modelled waves, near-bottom tidal current and circulation current data for a 3-year period were used in a widely applied sediment transport module to simulate the seabed shear stresses and the mobilization of observed sediment grain size on the Scotian Shelf of eastern Canada.

The Scotian Shelf is affected by strong waves and tidal currents. These waves, currents and/or their interaction cause maximum mean bed shear velocities of 5 – 10 cm s−1. Observed sediments on the Scotian Shelf can be mobilized by tidal currents at least once during the modelled 3 year period over 28% of the shelf area while waves can mobilize sediments over 60% of the shelf area suggesting much stronger sediment mobilization by waves. Interaction between waves and currents can produce enhanced combined wave-current shear velocity that is capable to mobilize sediments over 74% of the shelf area. The spatial variation of the relative importance of sediment mobilization frequency by component processes was used to classify the Scotian Shelf into six disturbance types. In comparison with previous studies using depth-averaged tidal currents, the present study based on near-bottom tidal currents has resulted in reduced sediment mobilization frequency by tidal currents, smaller extent of high mobility areas and significant changes of the spatial pattern of disturbance type distribution on the Scotian Shelf. The universal Seabed Disturbance Index and Sediment Mobility Index have also been applied to quantify the seabed exposure to physical processes and sediment mobilization on the Scotian Shelf by accounting for both the magnitude and frequency of these processes. The results of this modelling study are important for environmental assessments and for the spatial planning and management of the Scotian Shelf bioregion. 

How to cite: Li, M., Wu, Y., Ma, Y., and Wang, Y.: Modelling seabed shear stresses and sediment mobilization on the Scotian Shelf, eastern Canada, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2864, https://doi.org/10.5194/egusphere-egu23-2864, 2023.

08:37–08:39
|
PICO3a.2
|
EGU23-10870
|
On-site presentation
Christina Tsai and Meng-Jie Wu

Studies have categorized turbulent bursting events into outward interactions, ejections, inward interactions, and sweeps. Among these events, ejections and sweeps contribute notably to the time consumption, momentum flux, and sediment flux. Studies have reported that the distribution of turbulent coherent structures is uniform. However, research has revealed that the distribution of bursting events is nonuniform at different bed elevations. Although the nonuniform distribution of turbulent bursting events has been investigated, their influence on sediment transport has yet to be examined. This study established an improved stochastic diffusion particle tracking model (SD-PTM) using the stochastic Lagrangian method to describe sediment particle movement. This model integrates turbulent characteristics determined using a direct numerical simulation data set for comprehensively analyzing the sediment particle motion during turbulent flow. We developed a modified SD-PTM that considers the nonuniform spatial distribution of ejection and sweep events and the particle movement direction during these events. Particle trajectories were obtained using this model, and the anomalous diffusion during sediment transport was analyzed by calculating the variance in the particle trajectories. The performance of the proposed model was evaluated by comparing the flow velocities and sediment concentration profiles obtained using it with those measured in previous studies. Therefore, the sediment particle motion during turbulent flow was comprehensively investigated under extreme flow conditions.

How to cite: Tsai, C. and Wu, M.-J.: Stochastic Modeling for Anomalous Suspended Sediment Transport, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10870, https://doi.org/10.5194/egusphere-egu23-10870, 2023.

08:39–08:41
|
PICO3a.3
|
EGU23-16582
|
On-site presentation
|
Yongpeng Lin, Chenge An, Matthew Czapiga, Marwan Hassan, Gary Parker, and Xudong Fu

Since the 2008 Wenchuan (Ms. 8.0) earthquake, the Shi‑ting River in Sichuan Province, China has suffered massive bed incision, with the largest incision depth being more than 20 m in 7 years. Potential reasons include: breaks in sediment connectivity due to widespread in-channel weirs; the supply of fine sediment after the earthquake; intensive sand mining, etc. In this study, we simulate the combined role of in-channel weirs and fine sediment supply in determining the massive bed incision in the Shi-ting River. A one-dimensional river morphodynamic model is implemented. The simulated results show that the in-channel weirs can lead to bed incision and bed coarsening in the downstream channel. For a weir with a height of 5 m, the maximum incision depth is about 5 m, and the extent of downstream incision is no more than 20 km within 20 years. The supply of fine sediment can enhance the downstream channel incision, as the weir preferentially traps coarse sediment but passes the fine sediment downstream. However, a combination of in-channel weirs and the fine sediment supply cannot explain the dramatic incision (20 m in 7 years) as observed in the Shi-ting River. This suggests that the mining of gravel and sand had a significant role in driving channel degradation.

How to cite: Lin, Y., An, C., Czapiga, M., Hassan, M., Parker, G., and Fu, X.: Morphodynamics of the Shi-ting River after Wenchuan Earthquake: Effects of in-channel weirs and fine sediment supply on river incision, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-16582, https://doi.org/10.5194/egusphere-egu23-16582, 2023.

08:41–08:43
|
PICO3a.4
|
EGU23-7781
|
ECS
|
On-site presentation
Bastian Ringleb and Markus Fuchs

Landscape Evolution Models (LEMs) are powerful tools for simulating erosion and accumulation patterns over a landscape. To accurately capture the complexity of the system, several parameters are used to describe the flow of water on a surface and thereby calculate a change in landscape. However, due to the need for reasonable computing power, simplification is often necessary. The CAESAR-Lisflood LEM simplifies the transition from precipitation to runoff with a cell storage-based system, where water is retained in each cell and released at later iterations, simulating infiltration and forming the hydrograph. This reduced complexity works well for surface processes, but can lead to issues in catchments, where ground water flow and infiltration are unknown parameters, and an unknown proportion of sediment transport is connected to sub-surface processes.

Located in the Northern Franconian Jura, Germany, the Weismain River catchment (~125 km2) is a particularly challenging area to model the sediment dynamics, due to its karstic geology. The Weismain river and its tributaries are deeply incised into a limestone plateau forming small, well-defined valleys opening to wider floodplains in the lower parts of the catchment, where sandstone is dominant. To gain a better understanding of the sediment dynamics and the evolution of this catchment, the impact of the karstic environment needs to be evaluated. With CAESAR-Lisflood we are looking at the spatial distribution of alluvial sediment from different model inputs and how differences in discharge calculations affect model results. With this model setup we can gain insight into the sediment dynamics of the catchment and increase confidence in future results.

How to cite: Ringleb, B. and Fuchs, M.: The impact of geology on sediment dynamics. A modelling approach for the Northern Franconian Jura, Germany., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7781, https://doi.org/10.5194/egusphere-egu23-7781, 2023.

08:43–08:45
|
PICO3a.5
|
EGU23-2327
|
On-site presentation
Jürgen Mey, Wolfgang Schwanghart, Angela Landgraf, and Philippe Davy

Repeated alpine glaciations during the Quaternary partly reached far into the foreland and caused profound landscape changes beyond glacial margins in northern Switzerland. Climate-driven glacier growth and decay and commensurate variations in water and sediment delivery caused river systems to aggrade and to incise, leading to the widespread occurrence of glacio-fluvial deposits (Deckenschotter) and associated terraces. Mapping and numerical dating of these depositional complexes increasingly offer insights into the spatial patterns and timing of Quaternary glaciations, and associated changes in (glacio-)fluvial dynamics. An improved understanding of how fluvial systems respond to glacial-interglacial cycles will help to assess the erosion potential around repository sites of nuclear waste over the next one million years. In this study, we contribute to close this research gap using numerical landscape evolution modelling (LEM).

We use EROS, a numerical landscape evolution model, which implements a particle-based approach to simulate water and sediment fluxes that interact with topography through erosive and depositional actions. Unlike LEMs based on the stream-power incision law, the method solves the 2D shallow water equations with both basal and lateral erosion and deposition, which allows for variations in width and lateral mobility of rivers; these variations induce changes in the transport capacity of the sediments that cause specific patterns of deposition and erosion to emerge. We adjusted the model so that it is capable to run over 1 Myrs, and imposed boundary conditions that - informed by estimates on longterm erosion rates – reflect rock uplift and plausible variations in water and sediment fluxes following a 100-kyrs glacial cycle. Our model relies on digital elevation models and sediment thickness data with 60 m spatial resolution and is applied to the Hochrhein river between Stein am Rhein and Basel and the Aare river downstream of the area, where Limmat and Reuss enter.

Our simulations show that the model reproduces widespread aggradation, reworking of the sediments by highly laterally mobile, braided river systems and incision during periods of increased runoff and low sediment availability.

Our model setup and parametrization features several uncertainties. For example, the capacity of rivers to laterally erode strongly determines the thickness and extent of depositional complexes lining the Hochrhein and Aare system. Also, our model is sensitive to temporally varying boundary conditions of water and sediment input about which precise estimates are lacking. Regardless, the more detailed and realistic representation of hydraulic and sediment transport processes by EROS compared to conventionally used landscape evolution models at this spatial and temporal scale provide the opportunity to test different hypotheses using numerical experiments and link the results to field evidence. Further sensitivity analyses and uncertainty quantification will enable us to use our model as simulation tool to hindcast and investigate the behaviour of fluvial system in response to different tectonic and climatic scenarios, thus helping to better understand potential spatial patterns of and sediment assemblages within widespread glacio-fluvial deposits.

How to cite: Mey, J., Schwanghart, W., Landgraf, A., and Davy, P.: Fluvial response to glacial-interglacial cycles - modelling the evolution of the Hochrhein using EROS, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2327, https://doi.org/10.5194/egusphere-egu23-2327, 2023.

08:45–08:47
|
PICO3a.6
|
EGU23-14510
|
On-site presentation
Julien Babault, Lewis A. Owen, Pierre Arroucau, María Charco, Ludovic Bodet, Jean Van Den Driessche, and Marc Caffee

Records of incision history such as topographic data and landform dating can be gathered in inversion schemes to reconstruct base-level fall and uplift history. Here, we use a non-linear inversion scheme and the stream power incision model to study the landscape evolution of a mountainous area to quantify whether it responds to an uplift or a capture-induced local base-level fall. Our inversion is constrained by 10Be cosmogenic nuclide data. The probabilistic approach yields an ensemble of solutions made by a combination of model parameters. We apply our model to the Draa Canyon located in the Anti-Atlas of Morocco at the outlet of the Ouarzazate Basin which was internally drained during the Miocene. We show that incision rates in the Ouarzazate Basin and the southern margin of the High Atlas are compatible with a Pleistocene age for the opening of the Basin. The forcing to drainage integration may be due to capture by regressive erosion in the proto Draa river or tilting to the south of the High Atlas, Ouarzazate basin, and Anti-Atlas as a whole in response to mantle-related continental-scale uplift, or a combination of both. The southern border of the High Atlas in this region displays a transient landscape previously interpreted as evidence for recent shortening and rock uplift. Our results suggest that the rejuvenation of the southern Central High Atlas, in the northern margin of the Ouarzazate Basin, mainly occurred in response to the opening of the Ouarzazate Basin, with only several hundreds of meters of rock uplift localized along the South Atlas front during the late Cenozoic.

How to cite: Babault, J., Owen, L. A., Arroucau, P., Charco, M., Bodet, L., Van Den Driessche, J., and Caffee, M.: Pleistocene opening of the Ouarzazate Basin, and incision rate history of the Draa canyon in the Anti-Atlas of Morocco revealed by 10Be cosmogenic nuclide dating and non-linear river profile inversion, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-14510, https://doi.org/10.5194/egusphere-egu23-14510, 2023.

08:47–08:49
|
PICO3a.7
|
EGU23-6003
|
ECS
|
On-site presentation
|
Marina Ruiz Sánchez-Oro, Simon M Mudd, and Boris Gailleton

Fluvial incision patterns help us understand the role of precipitation in river formation and evolution. Incision is generally thought to be a function of channel gradient and either discharge or drainage area. Many studies use drainage area as a proxy for discharge, but incision could potentially be related to area rather than discharge because area determines sediment fluxes. Assumptions about discharge or drainage area control on incision impact the calculations of chi profiles and steepness indices, which are indicators of the erosional and tectonic history of the landscape. Here we use numerical modelling and topographic analysis to determine if spatially varied precipitation leads to distinctive patterns of channel profiles that can be used to differentiate if a discharge or drainage-area driven incision rule is appropriate for a given landscape.  

We use the modelling framework Fastscape to set up two scenarios: one including an incision rule with drainage area and another with discharge (rainfall). We evolve them to steady state and extract from the resultant topography a disorder metric, which measures the relation between the tributaries and the main river in chi space. To simulate a blind test of the incision conditions, we calculate the disorder metric under the two incision case assumptions. We find that scenarios in which incision is discharge-driven cases are statistically distinguishable from those in which incision is driven by drainage area alone. 

We repeat these analyses with variable rainfall gradients, simulating those found in real landscapes. The results are in agreement: in a sandbox model with homogeneous lithology, it is possible to distinguish the incision case based on the disorder metric under all rainfall conditions. Real landscapes are not sandboxes, however, so we run further simulations to quantify the impact of heterogeneous lithology under different rainfall scenarios. This reveals that it is not possible to distinguish the incision case, as the distortions to channel profiles driven by lithology masks evidence of discharge dominance. 

We complete the analysis by testing the hypothesis in 10 real landscapes across the globe, with different lithological and climatic regimes. Some study areas, such as the Pyrénées or the Alburz Mounts, show a weak trend towards discharge dominance, although there is no conclusive statistical evidence of the preference of one incision rule, hinting at the dominance of lithology over rainfall as shown in the models.

How to cite: Ruiz Sánchez-Oro, M., Mudd, S. M., and Gailleton, B.: Determining the correct river incision mechanism under different climatic and lithological conditions using disorder metrics, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6003, https://doi.org/10.5194/egusphere-egu23-6003, 2023.

08:49–08:51
|
PICO3a.8
|
EGU23-14815
|
ECS
|
On-site presentation
|
Fergus McNab, Taylor Schildgen, Jens Turowski, and Andrew Wickert
Alluvial rivers transport sediment from erosional source regions to sedimentary basins downstream. Erosion and precipitation rates control their sediment and water supplies, tectonic processes drive uplift or subsidence along their lengths, while sea or lake levels set their base levels. Changes in these external conditions cause them to adjust their slopes and sediment transport rates. Geomorphic and stratigraphic archives that develop alongside and downstream of alluvial rivers may therefore record information about past climatic and/or tectonic change and their influence on landscapes. Much recent work has aimed to understand precisely how alluvial rivers respond to changes in external forcing, using a range of numerical and laboratory modelling approaches. Many of these conceptual studies have used a simplified, one-dimensional spatial domain in which sediment and water discharge are either held constant or increase continuously downstream (e.g., according to Hack's law). Such studies have emphasised, among other important findings, that the timescales over which a system responds to external forcing is strongly influenced by its length. However, in real rivers, water and sediment accumulate at discrete intervals as tributary streams coalesce. This discrepancy complicates the application to real catchments of concepts developed in one-dimensional modelling studies, since, for example, the 'length' of a river network is not clearly defined. Here, we explore how incorporating realistic network geometries influences the behaviour of an alluvial river model. We use a model describing the long-profile evolution of transport-limited gravel-bed rivers which takes a non-linear diffusive form. We construct networks by linking individual segments, so that their sediment and water supplies are set by segments immediately upstream and their base levels are set by segments immediately downstream. We show that significant complexity can arise locally, so that studies aiming to understand specific segments within a catchment should take into account the geometry of that catchment. However, properties that integrate over the entire catchment, such as its total sediment export, are adequately predicted by a simplified one-dimensional model—provided an appropriate lengthscale is chosen. We conclude that, while care is required in some circumstances, one-dimensional models can provide useful insights into the general behaviour of alluvial river networks.

How to cite: McNab, F., Schildgen, T., Turowski, J., and Wickert, A.: Network geometries influence responses of alluvial river systems to external forcing, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-14815, https://doi.org/10.5194/egusphere-egu23-14815, 2023.

08:51–08:53
|
PICO3a.9
|
EGU23-7602
|
ECS
|
On-site presentation
Anthony Hoskins, Mikael Attal, Simon Mudd, and Miguel Castillo

Any tectonic forcing acting on a mountain range consists of vertical uplift motion and horizontal advection motion. In landscapes with a high advective component, such as those affected by low angle normal faults, advection can form a substantial portion of the total tectonic forcing. Questions therefore remain relating to the role of advection in shaping mountain range topography, surface drainage patterns and erosion rate distributions. We ask, how does advection influence mountain range topography and drainage? How do drainage basins and divides respond to advection? And how does advection influence the spatial distribution of erosion rates? Through numerical modelling with the Channel-Hillslope Integrated Landscape Development (CHILD) model and comparison to a natural landscape, the Sierra de la Laguna (Mexico), we test the extent to which advection affects mountain range evolution. Advection is shown to alter surface drainage patterns by promoting catchment elongation and reducing outlet spacing at the mountain front. The mountain range’s Main Drainage Divide (MDD) exhibits faster erosion on the distal flank’s headwaters; relative to the proximal flank’s headwaters, leading to a migration of the MDD towards the fault. Steady-state is achieved when the MDD migrates towards the fault at a rate equal to the rate at which the footwall is advected away from the fault. This pattern of erosion rates is found across all spatially fixed geomorphic features. Cosmogenic radionuclide analysis for catchment-averaged erosion rates in the Sierra de la Laguna demonstrates a difference in erosion rates across the MDD that is also indicative of a migration of the MDD towards the fault. Topographic, drainage and erosion rate observations for the Sierra de la Laguna are consistent with the numerical modelling when advection is included, suggesting that advection induces continual drainage migration whilst maintaining mountain range flank widths. We suggest that highly elongate catchments are indicative of substantial advection, and we also identify a mechanism through which divide migration can occur without altering the size of adjacent drainage basins.   

How to cite: Hoskins, A., Attal, M., Mudd, S., and Castillo, M.: Drainage Divide Migration in Steady-State Landscapes Experiencing Horizontal Tectonic Advection, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7602, https://doi.org/10.5194/egusphere-egu23-7602, 2023.

08:53–08:55
|
PICO3a.10
|
EGU23-3421
|
ECS
|
On-site presentation
Boris Gailleton, Philippe Steer, Philippe Davy, Wolfgang Schwanghart, and Thomas Bernard

Integrating hydrology in Landscapes Evolution Models (LEMs) is challenging. The drainage-area-based solutions, where drainage area weighted by precipitation rates approximate the amount of water flowing through every location, accumulates drainage-area downstream following the topographic gradient and has been empirically linked to observed discharge. While straightforward and computationally efficient, it implicitly includes hydrology without calculating discharge or water height. A more sophisticated solution consist in the direct calculation of the shallow water equation, which explicitly approximates water height and discharge using physics-based equations. While the latter bears more information about the channels and floodplains dynamics, it is inherently limited to short time scales and is computationally more expensive, with numerical time step typically of the order of the second – making its use for long-term LEMs particularly challenging. Here, we present GraphFlood, a fast iterative method computing river depth and water discharge in 2D on a digital elevation model (DEM). This new method leverages the Directed Acyclic Graph (DAG) nature of water flowing on surface topography to iteratively solve for the 2D shallow water equation without the inertia terms. The main idea of the algorithm is to find the correct water surface height by iteratively balancing discharge input and output. At each iteration, we first use fast DAGs-related algorithms to calculate flow accumulation on the hydraulic surface. We use this as an approximation of the discharge input. Then, the discharge output is calculated using the Manning flow resistance equation, in a manner similar to the Floodos model (Davy et al., 2017). The divergence of the discharges increments the water height. The iterative process is repeated until reaching a stationary state, i.e. a static field of water height and discharge representing an equilibrium state. Note that this method can be slightly modified to solve flood wave propagation by approximating the input discharge function of the immediate upstream neighbours. Water depths obtained with the stationary solution were validated against an analytical solution in the case of a rectangular channel and with the Floodos model for natural DEMs. Compared to previous hydrodynamic models, the main benefits of GraphFlood are its simplicity of implementation, which mainly requires a classical flow routing algorithm, and its efficiency. While case-dependent, our tests suggested a ~10 times speed-up compared to Floodos model (Davy et al., 2017) which was already significantly faster than other hydrodynamic models. Moreover, the computational time scales a little more than linearly with the number of cells, which makes GraphFlood a suitable solution even for DEMs larger than 106 – 107 cells. We demonstrate the suitability of the method for integrating realistic hydrology in a wide range of topographic and morphometric analyses (e.g. channel width assessment, floodplain delineation, Flint-Morisawa metrics) and in LEMs – even for longer timescales. 

How to cite: Gailleton, B., Steer, P., Davy, P., Schwanghart, W., and Bernard, T.: GraphFlood: a fast stationary solution for 2D hydrodynamics in landscapes evolution models, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-3421, https://doi.org/10.5194/egusphere-egu23-3421, 2023.

08:55–08:57
|
PICO3a.11
|
EGU23-10027
|
ECS
|
On-site presentation
|
Shashank Kumar Anand, Matteo B. Bertagni, Arvind Singh, and Amilcare Porporato

Steep landscapes maintain predominantly planar hillslopes over a range of spatial scales. These hillslopes are bounded at a typical angle, beyond which shallow landslides or slope failures remove the excess sediment influx. The evolution of such steep topographies with threshold hillslopes is well approximated by the eikonal equation, which is well-known in problems of geometry, optics, and mechanics. According to this approximation, hillslopes meet upstream to construct a network of sharp ridges and join downstream at the boundary to produce a complementary valley network. We find a good agreement between the proposed approximation and the spatial organization of landscapes with a dominant control of landslide erosion and negligible fluvial erosion. We also show that the eikonal approximation can be utilized to reconstruct the landscapes with threshold hillslopes where fluvial erosion sets the downstream free boundary.

How to cite: Anand, S. K., Bertagni, M. B., Singh, A., and Porporato, A.: Eikonal approximation for landscapes dominated by threshold hillslopes, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10027, https://doi.org/10.5194/egusphere-egu23-10027, 2023.

08:57–08:59
|
PICO3a.12
|
EGU23-10803
|
Highlight
|
On-site presentation
|
Greg Hancock and Jose Martín Duque

New computing and modelling technology allows the reconstruction of post-mining landforms using geomorphic design principles. It is important that such designs be evaluated and if needed, re-designed or reshaped so that soil loss is minimized and geomorphologically and therefore ecologically integrate with the surrounding landscape. One landscape assessment method is to use a computer based Landscape Evolution Models (LEMs). LEMs allow different designs to be input and will highlight where erosion will occur and type of erosion (i.e. sheetwash, rilling, gullying) as well as erosion rate. Here, landscape designs for the Santa Engracia mine in East-Central Spain were designed and constructed using geomorphic principles using the GeoFluv method and Natural Regrade software. These design landscapes were evaluated using the SIBERIA LEM. The results demonstrate that using suitable topsoil and if vegetation can be established the landscapes will have minimal sheet and gully erosion.  Erosion forecasts (5.3 - 15.2 t ha-1 yr-1) are significantly lower than the previous landscape (~350 t ha-1 yr-1) using conventional (terraced) reconstruction methods. The design and assessment methods described here provide procedure that can be used at other sites to highlight strengths and weakness of a landscape design process. The combination of geomorphic design and assessment using a landscape evolution at this project (LIFE RIBERMINE) presents a new standard for mine rehabilitation landform design in Europe.

How to cite: Hancock, G. and Martín Duque, J.: Using geomorphology to design, build and assess a post-mining landscape: a case study of the Santa Engracia mine, Spain, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10803, https://doi.org/10.5194/egusphere-egu23-10803, 2023.

08:59–09:01
|
PICO3a.13
|
EGU23-446
|
ECS
|
Virtual presentation
Nejib Hassen Abdullahi and Zulfequar Ahmad

Most of the alluvial rivers in the globe are suffering and necessitate critical management regarding sand and gravel extraction and other related river ecosystem protection. Sediment mining is the activity of extracting useful materials such as sand, gravel, and aggregate from a river bed, banks, and flood plains. Sediment mining activity causes hydro-morphodynamic variations on the river bed, impacting the river plan and hydraulic structures, and also can harm flora and fauna within the river ecosystem. The period and magnitude of the mining rate with the discharge and channel bed material properties directly impact channel bed mobility and bed-level equilibrium conditions. The flow field, bedload transport, and morphological evolution in the pit's vicinity are varied in space-time, multifaceted, and three-dimensional. Most of the preceding researchers on sediment mining characteristics were concerned with the physical aspects, scouring study, and degradation rate. Very few researches are available on the flow regions and morphological bathymetry of a river under sediment mining. On the bases of this area of gaps, the present study aimed to understand the flow field and bathymetry of the channel under sediment mining. The experimental work is conducted at the Hydraulics Laboratory of the Civil Engineering department, IIT Roorkee, India, on a trapezoidal sediment mining pit constructed of uniform cohesionless bed material in an open channel flow. The average streamwise velocity was measured at five various sections near the sediment mining pit using a three-dimensional acoustic doppler velocimeter. The numerical simulation was conducted in a Flow-3D solver using the Reynolds Averaged Navier–Stokes equations (RANS), standard k-ɛ, the volume of fluid, and FAVOR, and the results were compared to the experimental observations. Both approaches depict an increase in the average longitudinal velocity at the upstream and downstream nick points. The observed longitudinal velocity from the experiment at the downstream nick point is higher than the velocity simulated using CFD at the downstream nick point. Both approaches indicated a significant increase in the longitudinal velocity downstream of the pit center, especially in the flume center. CFD simulations depicted the decreasing velocity at the upstream nick point due to more degradation and increment of flow depth at the section. The detachment and degradation of the upstream nick were observed in the initial stage of the experimental work but not at the downstream nick point. Because the pit is used as a bedload trap, the likely sediment-free water compelled out of the pit causes bed erosion and flattening in the downstream section of the sediment mining pit. The degradation at the upstream nick was 4.6% and 9% for the experimental observations and the CFD simulations, while it was 21.3% and 3.1% for the downstream nick, respectively. The findings of this study can help authorities and experts in the effective maintenance and supervision of river ecosystem balance by supplying cost-effective sediment resources.

How to cite: Abdullahi, N. H. and Ahmad, Z.: Experimental and CFD Studies on the Flow Field and Bed-Morphology in the vicinity of a Sediment Mining Pit, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-446, https://doi.org/10.5194/egusphere-egu23-446, 2023.

09:01–09:03
|
PICO3a.14
|
EGU23-4275
|
Virtual presentation
Daisuke Harada and Shinji Egashira

Recently, several extreme flood hazards with active sediment and driftwood transport have occurred in mountainous areas of Japan, as typically observed in the Akatani river flood disaster in 2017. Therefore, in order to mitigate such hazards, it is required to develop methods to evaluate flood, sediment, and driftwood runoff from a basin during heavy rainfall. This research proposes integrated methods to evaluate flood, sediment and driftwood runoff from a basin during heavy rainfall using numerical methods. The proposed methods are applied to the Terauchi Dam basin, where large amounts of sediment and driftwood were discharged during the 2017 Northern Kyushu heavy rainfall event, to discuss the applicability of the methods.

The methods are composed of a distributed rainfall-runoff model, slope stability analysis, sediment and driftwood transport in the slope based on mass system equation, and sediment and driftwood transport in the river channel by unit channel model. These models are integrated as Rainfall-Sediment Runoff (RSR) model to evaluate sediment and driftwood runoff from a basin. As a result of its application to the Terauchi Dam basin, we found that whether or not the debris flow enters the river channel, i.e., the definition of the upstream end of the river channel, has a significant impact on the results. Therefore, we investigated whether the sediment supplied by debris flow enter the river channel during the actual event. The results show that the upstream end of the river channel in the computation should be defined as approximately 4% to 10% slope, and in case the mesh size is sufficiently fine, the debris flow inflow into the channel is sufficiently evaluated. In conclusion, this research proposes an integrate methods to evaluate flood, sediment and driftwood runoff from a basin, and discusses its applicability to the disasters such as the 2017 Northern Kyushu heavy rainfall event.

How to cite: Harada, D. and Egashira, S.: Method to evaluate sediment-driftwood transport processes with flood runoff in a basin during heavy rainfalls, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-4275, https://doi.org/10.5194/egusphere-egu23-4275, 2023.

09:03–10:15