Search

Since the 1990s, analysis of cosmogenic nuclides, primarily 10Be, in quartz-bearing river sand, has allowed for quantitative determination of erosion rates at a basin scale. Paired measurements of in situ cosmogenic 26Al and 10Be in sediment are less common but offers insight into the history of riverine sediment moving down slopes and through drainage basins. Prolonged sediment burial (>105 years), a violation of assumptions underlying erosion rate calculations, is indicated by higher 26Al-based than 10Be-based erosion rates due to preferential loss of shorter-lived 26Al by decay when quartz is shielded from cosmic rays. Here, we use a global compilation of 26Al and 10Be data generated from quartz-bearing fluvial sediment samples (n = 624, including 121 new measurements) and calculate the discordance between erosion rates derived from each nuclide. We test for correlations between such discordance and topographic metrics for drainage basins, allowing us to infer the likelihood of sediment burial during transport in different geomorphic settings. We find that nearly half of samples (n = 276) exhibit discordance (> 1σ uncertainty) between erosion rates derived from 10Be and 26Al, indicating sediment histories that must include extended burial during residence on hillslopes and/or in the fluvial system after or during initial near-surface exposure. Physical basin parameters such as basin area, slope, and tectonic activity exhibit significant correlation with erosion rate discordance whereas climatic parameters have little correlation. Our analysis suggests that 26Al/10Be erosion rate discordance occurs more regularly in basins larger than 1,000 km2, particularly when such basins have low average slopes and are in tectonically quiescent terrains. Sediment sourced from smaller, steeper basins in tectonically active regions is more likely to have similar 10Be and 26Al erosion rates indicative of limited storage and limited burial during residence in the hillslope and fluvial sediment system. The data and analysis we present demonstrate that paired 26Al and 10Be analyses in detrital fluvial samples can provide a window into watershed processes, elucidating landscape behavior at different spatial scales and allowing a deeper understanding of both sediment routing systems and whether erosion rate assumptions are violated. Large lowland basins are more likely to transport detrital sediment that has experienced prolonged sediment storage and burial either on hillslopes and/or in fluvial networks; thus, erosion rates from such basins are lower limits due to nuclide decay during storage. Conversely, samples from smaller upland basins are more likely to provide reliable erosion rates. [ABSTRACT FROM AUTHOR]

Landslides influence fluvial suspended sediment transport by changing sediment supply and grain size, which alter suspended sediment concentrations and fluxes for a period of time after landsliding. To investigate the duration and scale of altered suspended sediment transport due to landsliding, we analyzed suspended sediment concentration and water discharge measurements at 87 gauging stations across Taiwan over an 11-year period after Typhoon Morakot, which generated nearly 20 000 landslides in 2009. At each gauging station, we computed annual rating curves to quantify changes over time in the sensitivity of suspended sediment concentrations to water discharge. Among the 40 stations in basins that were impacted by landsliding, the discharge-normalized rating curve coefficient ã was higher than that before Morakot by a factor of 5.1±1.1 (mean ± standard error) in 2010, the first year after Morakot. The rating curve exponent b did not decrease at most stations until a year later (2011), when the average b value was lower than that before Morakot by 0.25±0.05. Across the compilation of gauging stations, post-Morakot changes in discharge-normalized sediment concentration (ã) were positively correlated with landslide intensity for 7 years after Morakot, while post-Morakot changes in the exponent of the discharge–concentration relationship (b) were negatively correlated with landslide intensity from 2011 to 2014. This reflects a tendency for larger changes in ã and b to occur in basins with more intense landsliding. At 26 of these 40 stations, elevated values of ã declined after the initial post-Morakot peak, consistent with a gradual return to pre-Morakot suspended sediment transport conditions. Exponential regressions to these ã values reveal a median characteristic decay time of 8.8 years (interquartile range: 5.7–14.8 years). Values of ã increased more and declined faster in basins with more intense landsliding, with a mean characteristic decay time of 6 years in the basins hit hardest by landsliding. Furthermore, changes in ã and b tended to be larger in basins with more intense landsliding. At stations that were not impacted or only minimally impacted by landsliding, neither ã nor b exhibited systematic responses to Morakot. To quantify the effect of landsliding on sediment discharge, we compared the measured sediment discharges after Morakot to the hypothetical sediment discharges that would have occurred if Morakot had induced no landslides, calculated by applying each station's pre-Morakot rating curve to its post-Morakot water discharge history. This analysis suggests that Morakot-induced landsliding increased sediment discharge by as much as > 10-fold in some basins in the 1–2 years after Morakot. Together, these results indicate that the influence of Morakot-induced landsliding on rating curves was large shortly after Morakot but diminished in less than a decade in most of the study rivers and will be imperceptible in another few decades in all of the study rivers. To the extent that these results are applicable to other landscapes, this suggests that periods of elevated sediment transport efficiency after landsliding should persist for years to decades, even if the landslide deposits persist for centuries to millennia. [ABSTRACT FROM AUTHOR]

The Chilean Coastal Range, located in the Mediterranean segment of Chile, is a soil-mantled landscape with the potential to store valuable freshwater supplies and support a biodiverse native forest. Nevertheless, human intervention has been increasing soil erosion for ∼ 200 years, culminating in the intense management of exotic tree plantations throughout the last ∼ 45 years. At the same time, this landscape has been severely affected by a prolonged megadrought. As a result, this combination of stressors complicates disentangling the effects of anthropogenic disturbances and hydroclimatic trends on sediment fluxes at the catchment scale. In this study, we calculate decennial catchment erosion rates from suspended-sediment loads and compare them with a millennial catchment denudation rate estimated from detrital 10Be. We then contrast both of these rates with the effects of discrete anthropogenic-disturbance events and hydroclimatic trends. Erosion and denudation rates are similar in magnitude on decennial and millennial timescales, i.e., 0.018 ± 0.005 and 0.024 ± 0.004 mmyr-1 , respectively. Recent human-made disturbances include logging operations throughout all seasons and a dense network of forestry roads, thereby increasing structural sediment connectivity. Further disturbances include two widespread wildfires (2015 and 2017) and an earthquake with an Mw value of 8.8 in 2010. We observe decreased suspended-sediment loads during the wet seasons for the period 1986–2018, coinciding with declining streamflow, baseflow, and rainfall. The low millennial denudation rate aligns with a landscape dominated by slow diffusive soil creep. However, the low decennial erosion rate and the decrease in suspended sediment disagree with the expected effect of intense anthropogenic disturbances and increased structural (sediment) connectivity. Such a paradox suggests that suspended-sediment loads, and thus respective catchment erosion, are underestimated and that decennial sediment detachment and transport have been masked by decreasing rainfall and streamflow (i.e., weakened hydroclimatic drivers). Our findings indicate that human-made disturbances and hydrologic trends may result in opposite, partially offsetting effects on recent erosion, yet both contribute to landscape degradation. [ABSTRACT FROM AUTHOR]

The recession of Alpine glaciers causes an increase in the extent of proglacial areas and leads to changes in the water discharge and sediment balance (morphodynamics and sediment transport). Although the processes occurring in proglacial areas are relevant not only from a scientific point of view but also for the purpose of climate change adaptation, there is a lack of work on the continuous monitoring and multitemporal characterization of these areas. This study offers a multidisciplinary approach that merges the contributions of different scientific disciplines, such as hydrology, geophysics, geomatics, and water engineering, to characterize the Rutor Glacier and its proglacial area. Since 2020, we have surveyed the glacier and its proglacial area using both uncrewed and crewed aerial surveys (https://doi.org/10.5281/zenodo.8089499 , ; https://doi.org/10.5281/zenodo.10100968 , ; https://doi.org/10.5281/zenodo.10074530 , ; https://doi.org/10.5281/zenodo.10101236 , ; https://doi.org/10.5281/zenodo.7713146 ,). We have determined the bathymetry of the most downstream proglacial lake and the thickness of the sediments deposited on its bottom (https://doi.org/10.5281/zenodo.7682072 ,). The water depth at four different locations within the hydrographic network of the proglacial area (https://doi.org/10.5281/zenodo.7697100 ,) and the bedload at the glacier snout (https://doi.org/10.5281/zenodo.7708800 ,) have also been continuously monitored. The synergy of our approach enables the characterization, monitoring, and understanding of a set of complex and interconnected processes occurring in a proglacial area. [ABSTRACT FROM AUTHOR]

In this paper, an Eulerian two-phase-flow model sedFoam is extended to include an air phase together with the water and sediment phases. The numerical model called sedInterFoam is implemented using the open source library OpenFOAM. SedInterFoam includes the previous features of sedFoam for sediment transport modeling and also solves the air/water interface using the volume of fluid method coupled with the waves2Foam toolbox for free surface wave generation and absorption. Using sedInterFoam, four test cases are successfully reproduced to validate the free-surface evolution algorithm implementation, mass conservation of sediment and fluid phases, predictive capabilities and demonstrate its potential in modelling a broader range of coastal applications with sediment transport dominated by surface waves. [ABSTRACT FROM AUTHOR]

Silt-rich meltwater plume deposits (MPDs) analyzed from marine sediment cores have elucidated relationships that are clearly connected, yet difficult to constrain, between subglacial hydrology, ice-marginal landforms, and grounding-zone retreat patterns for several glacial catchments. Few attempts have been made to infer details of subglacial hydrology, such as flow regime, geometry of drainage pathways, and mode(s) of sediment transport through time, from grain-scale characteristics of MPDs. Using sediment samples from MPD, till, and grounding-zone proximal diamicton collected offshore of six modern and relict glacial catchments in both hemispheres, we examine grain shape distributions and microtextures (collectively, grain micromorphology) of the silt fraction to explore whether grains are measurably altered from their subglacial sources via meltwater action. We find that 75 % of all imaged grains (n = 9400) can be described by 25 % of the full range of measured shape morphometrics, indicating grain shape homogenization through widespread and efficient abrasive processes in subglacial environments. Although silt grains from MPDs exhibit edge rounding more often than silt grains from tills, grain surface textures indicative of fluvial transport (e.g., v-shaped percussions) occur in only a modest number of grains. Furthermore, MPD grain surfaces retain several textures consistent with transport beneath glacial ice (e.g., straight or arcuate steps, (sub)linear fractures) in comparable abundances to till grains. Significant grain shape alteration in MPDs compared to their till sources is observed in sediments from glacial regions where (1) high-magnitude, potentially catastrophic meltwater drainage events are inferred from marine sediment records and (2) submarine landforms suggest supraglacial melt contributed to the subglacial hydrological budget. This implies that quantifiable grain shape alteration in MPDs could reflect a combination of high-energy flow of subglacial meltwater, persistent sediment entrainment, and/or long sediment transport distances through subglacial drainage pathways. Integrating grain micromorphology into analysis of MPDs in site-specific studies could therefore aid in distinguishing periods of persistent, well-connected subglacial discharge from periods of sluggish or disorganized drainage. In the wider context of deglacial marine sedimentary and bathymetric records, a grain micromorphological approach may bolster our ability to characterize ice response to subglacial meltwater transmission through time. This work additionally demonstrates that glacial and fluvial surface textures are retained on silt-sized quartz grains in adequate amounts for microtexture analysis, which has heretofore been conducted exclusively on the sand fraction. Therefore, grain microtextures can be examined on silt-rich glaciogenic deposits that contain little to no sand as a means to evaluate sediment transport processes. [ABSTRACT FROM AUTHOR]

Confluences are dynamic morphological nodes that are found in all river networks. In mountain regions, they are influenced by hydraulic and sedimentary processes that occur in steep channels during extreme events in small watersheds. Sediment transport in the tributary channel and aggradation in the confluence can be massive, potentially causing overbank flooding and sedimentation into adjacent settlement areas. Previous works dealing with confluences have mainly focused on lowland regions, and those that have focused on mountain areas have used sediment concentrations and channel gradients that are largely under-representative of mountain river conditions. The presented work contributes to filling this research gap with 45 experiments that use a large-scale physical model. Geometric model parameters, the applied grain size distribution, and the considered discharges represent the conditions at 135 confluences in South Tyrol (Italy) and Tyrol (Austria). The experimental program allowed for a comprehensive analysis of the effects of (i) the confluence angle, (ii) the tributary gradient, (iii) the channel discharges, and (iv) the tributary sediment concentration. In contrast to most research dealing with confluences, results indicate that, in the presence of an intense tributary sediment supply and a small tributary-to-main-channel discharge ratio (0.1), the confluence angle does not have a decisive effect on confluence morphology. Adjustments to the tributary channel gradient yielded the same results. A reoccurring range of depositional geomorphic units was observed in which a deposition cone transitioned to a bank-attached bar. The confluence morphology and tributary channel gradient rapidly adjusted, tending towards an equilibrium state to accommodate both water discharges and the sediment load from the tributary. Statistical analyses demonstrated that the confluence morphology was controlled by the combined channel discharge and the depositional or erosional extent was controlled by the sediment concentration. Applying conclusions drawn from lowland confluence dynamics could misrepresent depositional and erosional patterns and the related flood hazard at mountain river confluences. [ABSTRACT FROM AUTHOR]

Understanding event-driven sediment transport in coral reef environments is essential to assessing impacts on reef species, habitats, restoration, and mitigation, yet a global knowledge gap remains due to limited quantitative studies. Hurricane Irma made landfall in the Lower Florida Keys with sustained 209 km h -1 winds and waves greater than 8 m on 10 September 2017, directly impacting the Florida Reef Tract (FRT) and providing an opportunity to perform a unique comprehensive, quantitative assessment of its impact on coral reef structure and sediment redistribution. We used lidar and multibeam derived digital elevation models (DEMs) collected before and after the passing of Hurricane Irma over a 15.98 km 2 area along the lower FRT including Looe Key Reef to quantify changes in seafloor elevation, volume, and structure due to storm impacts. Elevation change was calculated at over 4 million point locations across 10 habitat types within this study area for two time periods using data collected (1) approximately 1 year before the passing of Irma and 3 to 6 months following the storm's impact as well as (2) 3 to 6 months after and up to 16.5 months after the storm. Elevation change data were then used to generate triangulated irregular network (TIN) models in ArcMap to calculate changes in seafloor volume during each time period. Our results indicate that Hurricane Irma was primarily a depositional event that increased mean seafloor elevation and volume at this study site by 0.34 m and up to 5.4 Mm 3 , respectively. Sediment was transported primarily west-southwest (WSW) and downslope, modifying geomorphic seafloor features including the migration of sand waves and rubble fields, formation of scour marks in shallow seagrass habitats, and burial of seagrass and coral-dominated habitats. Approximately 16.5 months after Hurricane Irma (during a 13-month period between 2017 and 2019), net erosion was observed across all habitats with mean elevation change of -0.15 m and net volume change up to -2.46 Mm 3. Rates of elevation change during this post-storm period were 1 to 2 orders of magnitude greater than decadal and multi-decadal rates of change in the same location, and changes showed erosion of approximately 50 % of sediment deposited during the storm event as seafloor sediment distribution began to re-equilibrate to non-storm sea-state conditions. Our results suggest that higher-resolution elevation change data collected over seasonal and annual time periods could enhance characterization and understanding of short-term and long-term rates and processes of seafloor change. [ABSTRACT FROM AUTHOR]

In the process of reworking sediments and thus shaping biogeochemical processes, marine bottom-dwelling animals are thought to play a pivotal role in many benthic environments. Bioturbation (particle reworking) includes the downward transport of particles into the sediment as a major process and is sometimes detected as subsurface maxima (peaks) of specific particulate substances (tracers). Here, we document the fact that subsurface peaks, such as those typically attributed to biological particle transport in sediments, may equally be generated by otter boards in bottom-trawling fishery. Boards can generate tracer peaks whereby they scoop sediment from the surface, flip it over, and deposit it onto the adjacent seafloor. These peaks are indistinguishable from those generated by benthic fauna burying surface material at sediment depth. We demonstrate this for the particle tracer chlorophyll a in silty sand from the western Baltic Sea with fauna that generally do not burrow deep in a global comparison. Our inability to distinguish the driving processes generating the peaks indicates limits to our understanding of the magnitude and spatial extent of bioturbation traces in this environment. It also poses a problem for the assessment of fishery resource use and benthic processes. However, based on natural fauna abundance, behavioral information, and fishery intensity data, we identify macrofauna and not otter boards as the dominant cause of peaks at the sites investigated here. [ABSTRACT FROM AUTHOR]

Satellite gravimetry is used to study the global hydrological cycle. It is a key component in the investigation of groundwater depletion on the Indian subcontinent. Terrestrial mass loss caused by river sediment transport is assumed to be below the detection limit in current gravimetric satellites of the Gravity Recovery and Climate Experiment Follow-On mission. Thus, it is not considered in the calculation of terrestrial water storage (TWS) from such satellite data. However, the Ganges and Brahmaputra rivers, which drain the Indian subcontinent, constitute one of the world's most sediment-rich river systems. In this study, we estimate the impact of sediment mass loss within their catchments on local trends in gravity and consequential estimates of TWS trends. We find that for the Ganges–Brahmaputra–Meghna catchment sediment transport accounts for (4 ± 2) % of the gravity decrease currently attributed to groundwater depletion. The sediment is mainly eroded from the Himalayas, where correction for sediment mass loss reduces the decrease in TWS by 0.22 cm of equivalent water height per year (14 %). However, sediment mass loss in the Brahmaputra catchment is more than twice that in the Ganges catchment, and sediment is mainly eroded from mountain regions. Thus, the impact on gravimetric TWS trends within the Indo–Gangetic Plain – the main region identified for groundwater depletion – is found to be comparatively small (< 2 %). [ABSTRACT FROM AUTHOR]

Sediment grains are progressively rounded during their transport down a river. For more than a century, Earth scientists have used the roundness of pebbles within modern sediment, and of clasts within conglomerates, as a key metric to constrain the sediment's transport history and source area(s). However, the current practices of assessment of pebble roundness are mainly qualitative and based on time-consuming manual measurement methods. This qualitative judgement provides the transport history only in a broad sense, such as classifying distance as "near" or "far". In this study, we propose a new model that quantifies the relationship between roundness and the transport distance. We demonstrate that this model can be applied to the clasts of multiple lithologies including modern sediment, as well as conglomerates, deposited by ancient river systems. We present field data from two Himalayan catchments in Nepal. We use the normalized isoperimetric ratio (IRn), which relates a pebble's area (A) to its perimeter (P), to quantify roundness. The maximum analytical value for IRn is 1, and IRn is expected to increase with transport distance. We propose a non-linear roundness model based on our field data, whereby the difference between a grain's IRn and the maximum value of 1 decays exponentially with transport distance, mirroring Sternberg's model of mass loss or size reduction by abrasion. This roundness model predicts an asymptotic behaviour for IRn , and the distance over which IRn approaches the asymptote is controlled by a rounding coefficient. Our field data suggest that the roundness coefficient for granite pebbles is 9 times that of quartzite pebbles. Using this model, we reconstruct the transport history of a Pliocene paleo-river deposit preserved at the base of the Kathmandu intermontane basin. These results, along with other sedimentary evidence, imply that the paleo-river was much longer than the length of the Kathmandu Basin and that it must have lost its headwaters through drainage capture. We further explore the extreme rounding of clasts from Miocene conglomerate of the Siwalik zone and find evidence of sediment recycling. [ABSTRACT FROM AUTHOR]

The construction of check dams is a common practice around the world where the aim is to reduce the damage by flooding events through mountain streams. However, quantifying the effectiveness of such engineering structures has remained very challenging and requires well-selected case studies, since the outcome of such an evaluation depends on site-specific geometric, geologic and climatic conditions. Conventionally, the check dams' effectiveness has been estimated using information about how the bedload sediment flux in the stream changes after the check dams are constructed. A permanent lowering of the bedload flux not only points to a success in reducing the probability of sediment transport occurrence but also implies that the sediment input through the system is likely to decrease. Here, we applied a method for data acquisition and two different equations (Meyer-Peter–Müller and Recking approach) to estimate and compare the sediment transport in a mountain stream in Switzerland under engineered and non-engineered conditions. Whereas the first equation is derived from a classical approach that is based on flume experiment data with a slope of less than 0.02 m m -1 , the second equation (Recking) has been derived based on a bedload field dataset comprising active mountain streams under steeper conditions. We selected the Guerbe (Gürbe) River situated in the Swiss Alps as a case study, which has been engineered since the end of the 19th century. This has resulted in more than 110 check dams along a ca. 5 km reach where sediment has continuously been supplied from adjacent hillslopes, primarily by landsliding. We measured the riverbed grain size, topographic gradients and river widths within selected segments along this reach. Additionally, a gauging station downstream of the reach engineered with check dams yielded information to calibrate the hydroclimatic situation for the study reach, thus offering ideal conditions for our catchment-scale experiment. Using the acquired data and the dataset about historical runoff covering the time interval between 2009 and 2021 and considering the current engineered conditions, we estimated a mean annual volume of transported bedload which ranges from 900 to 6000 m 3 yr -1. We then envisaged possible channel geometries before the check dams were constructed. We inferred (1) higher energy gradients which we averaged over the length of several check dams and which we considered a proxy for the steeper river slope under natural conditions; (2) channel widths that are smaller than those measured today, thereby anticipating that the channel was more confined in the past; and (3) larger grain size percentiles, which we consider to be similar to the values measured from preserved landslides in the region. Using such potential non-engineered scenarios as constraints, the two equations both point towards a larger sediment flux compared to the engineered state, although the results of these equations differed significantly in magnitude. Whereas the Recking approach returned estimates where the bedload sediment flux is ca. 10 times larger in comparison with the current situation, the use of the Meyer-Peter–Müller equation predicts an increase of ca. 100 times in bedload fluxes for a state without check dams. These results suggest that the check dams in the Guerbe River are highly efficient not only in regulating sediment transport by decreasing the probability of high sediment flux occurrence during torrential conditions but also in stabilizing the channel bed by avoiding incision. The most likely consequence is a stabilization of the terrain around such structures by reducing the activation of landslides. [ABSTRACT FROM AUTHOR]

Floods are among the most common natural disasters responsible for severe damages and human losses. Numerically produced data, managed by user-friendly tools for geographically referenced data, have been adopted to increase preparedness and reduce vulnerabilities. This paper describes the locally sensed and numerically produced data that characterize a flood event that occurred in February 2016 in the Portuguese Águeda River, referred to as Agueda.2016Flood for short. The data were managed through the RiverCure portal, a collaborative web platform connected to a validated shallow-water model featuring modelled dynamic bed geometries and sediment transport. The dataset provides a synthesis of topo-bathymetric, hydrometric and numerically produced data from a calibrated hydrodynamic model. Due to the lack of measured hydrometric data near the city, the numerically produced data are crucial for a complete description of the flood event. The Agueda.2016Flood dataset constitutes a complete validation test for flood forecasting models and a tool to better mitigate floods in this river and in similar rivers. Thus, Agueda.2016Flood is a relevant dataset for Águeda River stakeholders as well as for the community of flood modellers, as it provides a well-documented validation event for forecasting tools. The dataset can be accessed at 10.4211/hs.937927473a3a4e66a07a2e2fdd9d581e (Ricardo et al., 2024). [ABSTRACT FROM AUTHOR]

The grain size 2 mm is the conventional border between sand and gravel. This size is used extensively, and generally without much physical justification, to discriminate between such features as sedimentary deposit type (clast-supported versus matrix-supported), river type (gravel bed versus sand bed), and sediment transport relation (gravel versus sand). Here we inquire as to whether this 2 mm boundary is simply a social construct upon which the research community has decided to agree or whether there is some underlying physics. We use dimensionless arguments to show the following for typical conditions on Earth, i.e., natural clasts (e.g., granitic or limestone) in 20 ∘ C water. As grain size ranges from 1 to 5 mm (a narrow band including 2 mm), sediment suspension becomes vanishingly small at normal flood conditions in alluvial rivers. We refer to this range as pea gravel. We further show that bedload movement of a clast in the pea gravel range with, for example, a size of 4 mm moving over a bed of 0.4 mm particles has an enhanced relative mobility compared to a clast with a size of 40 mm moving over a bed of the same 4 mm particles. With this in mind, we use 2 mm here as shorthand for the narrow pea gravel range of 1–5 mm over which transport behavior is distinct from both coarser and finer material. The use of viscosity allows the delineation of a generalized dimensionless bed grain size discriminator between "sand-like" and "gravel-like" rivers. The discriminator is applicable to sediment transport on Titan (ice clasts in flowing methane/ethane liquid at reduced gravity) and Mars (mafic clasts in flowing water at reduced gravity), as well as Earth. [ABSTRACT FROM AUTHOR]

River bifurcations are constituent components of multi-thread fluvial systems, playing a crucial role in their morphodynamic evolution and the partitioning of water and sediment. Although many studies have been directed at exploring bifurcation dynamics, the conditions under which avulsions occur, resulting in the complete abandonment of one branch, are still not well understood. To address this knowledge gap, we develop a novel 1D numerical model based on existing nodal point relations for sediment partitioning, which allows for the simulation of the morphodynamic evolution of a free bifurcation. Model results show that when the discharge asymmetry is so high that the shoaling branch does not transport sediments (partial avulsion conditions) the dominant branch undergoes significant degradation, leading to a higher inlet step between the bifurcates and further amplifying the discharge asymmetry. The degree of asymmetry is found to increase with the length of the downstream channels to the point that when they are sufficiently long, the shoaling branch is completely abandoned (full avulsion conditions). To complement our numerical findings, we also formulate a new analytical model that is able to reproduce the essential characteristics of the partial avulsion equilibrium, which enables us to identify the key parameters that control the transition between different configurations. In summary, this research sheds light on the fundamental processes that drive avulsion through the abandonment of river bifurcations. The insights gained from this study provide a foundation for further investigations and may offer valuable information for the design of sustainable river restoration projects. [ABSTRACT FROM AUTHOR]

The difficulties of measuring bedload transport in gravel-bed rivers have given rise to the morphological method wherein sediment transport can be inferred from changes in riverbed elevation and estimates of the distance traveled by sediment: its path length. Because current methods for estimating path length are time- and labor-intensive, we present a method to estimate a characteristic path length from repeat digital elevation models (DEMs of difference, i.e., DoDs). We propose an automated method to extract the spacing between erosional and depositional sites on the DoD by the application of variational mode decomposition (VMD), a signal processing method, to quantify the spacing as a proxy for path length. We developed this method using flume experiments where bed topography and sediment flux were measured and then applied it to published field data with physical path length measured from tracer measurements. Our sediment transport estimates were not significantly different than the measured sediment flux at lower discharges in the lab. However, we observed an underestimation of sediment flux at the higher discharges in the flume study. We interpret this as a limit of the method in confined settings, where sediment transport becomes decoupled from morphological changes. We also explore how the time between survey acquisitions, the morphological active width relative to the channel width, and DoD thresholding techniques affect the proposed method and the potential issues they pose for the morphological method in general. [ABSTRACT FROM AUTHOR]

A proper understanding of sediment resuspension and transport processes is key to the morphodynamics of shallow tidal environments. However, a complete spatial and temporal coverage of suspended sediment concentration (SSC) to describe these processes is hardly available, preventing the effective representation of depositional dynamics in long-term modelling approaches. Through aiming to couple erosion and deposition dynamics in a unique synthetic theoretical framework, we here investigate SSC dynamics, following a similar approach to that adopted for erosion. The analysis with the peak-over-threshold theory of SSC time series computed using a fully coupled, bi-dimensional model allows us to identify interarrival times, intensities, and durations of over-threshold events and test the hypothesis of modelling SSC dynamics as a Poisson process. The effects of morphological modifications on spatial and temporal SSC patterns are investigated in the Venice Lagoon, for which several historical configurations in the last 4 centuries are available. Our results show that, similar to erosion events, SSC can be modelled as a marked Poisson process in the intertidal flats for all the analysed morphological lagoon configurations because exponentially distributed random variables describe over-threshold events well. Although erosion and resuspension are intimately intertwined, erosion alone does not suffice to describe also SSC because of the non-local dynamics due to advection and dispersion processes. The statistical characterization of SSC events completes the framework introduced for erosion mechanics, and together, they represent a promising tool to generate synthetic, yet realistic, time series of shear stress and SSC for the long-term modelling of tidal environments. [ABSTRACT FROM AUTHOR]

Reliable descriptions of erosion events are foundational to effective frameworks relevant to the fate of tidal landscape evolution. Besides the rhythmic, predictable action of tidal currents, erosion in shallow tidal environments is strongly influenced by the stochastic wave-induced bottom shear stress (BSS), mainly responsible for sediment resuspension on tidal flats. However, the absence of sufficiently long, measured time series of BSS prevents a direct analysis of the combined tide- and wave-driven erosion dynamics and its proper representation in long-term morphodynamic models. Here we test the hypothesis of describing erosion dynamics in shallow tidal environments as a Poisson process by analysing, with the peak-over-threshold theory, the BSS time series computed using a fully coupled, bi-dimensional numerical model. We perform this analysis on the Venice Lagoon, Italy, taking advantage of several historical surveys done in the last 4 centuries, which allow us to investigate the effects of morphological modifications on spatial and temporal erosion patterns. Our analysis suggests that erosion events on intertidal flats can effectively be modelled as a marked Poisson process in different morphological configurations because the interarrival times, durations, and intensities of the over-threshold exceedances are always well described by exponentially distributed random variables. The resulting statistical characterization allows a straightforward computation of morphological indicators, such as the erosion work, and paves the way for a novel synthetic, yet reliable, approach for the long-term morphodynamic modelling of tidal environments. [ABSTRACT FROM AUTHOR]

Based on measurements with the Swiss plate geophone system with a 1 min temporal resolution, bedload transport fluctuations were analysed as a function of the flow and transport conditions in the Swiss Erlenbach stream. The study confirms a finding from an earlier event-based analysis of the same bedload transport data, which showed that the disequilibrium ratio of measured to calculated transport rate (disequilibrium condition) influences the sediment transport behaviour. To analyse the transport conditions, the following elements were examined to characterise bedload transport fluctuations: (i) the autocorrelation coefficient of bedload transport rates as a function of lag time (memory effect), (ii) the critical discharge at the start and end of a transport event, (iii) the variability in the bedload transport rates, and (iv) a hysteresis index as a measure of the strength of bedload transport during the rising and falling limb of the hydrograph. This study underlines that above-average disequilibrium conditions, which are associated with a larger sediment availability on the streambed, generally have a stronger effect on subsequent transport conditions than below-average disequilibrium conditions, which are associated with comparatively less sediment availability on the streambed. The findings highlight the important roles of the sediment availability on the streambed, the disequilibrium ratio, and the hydraulic forcing in view of a better understanding of the bedload transport fluctuations in a steep mountain stream. [ABSTRACT FROM AUTHOR]

Confluences are dynamic morphological nodes in all river networks. In mountain regions, they are influenced by hydraulic and sedimentary processes occurring in steep channels during extreme events in small watersheds. Sediment transport in the tributary channel and aggradation in the confluence can be massive, potentially causing overbank flooding and sedimentation into adjacent settlement areas. Previous works dealing with confluences have been mainly focused on lowland regions, or if focused on mountain areas, the sediment concentrations and channel gradients were largely under-representative of mountain river conditions. The presented work contributes to filling this research gap with 45 experiments using a large-scale physical model. Geometric model parameters, applied grain size distribution, and the considered discharges represent the conditions at 135 confluences in South Tyrol (Italy) and Tyrol (Austria). The experimental program allowed for a comprehensive analysis of the effects of (i) the confluence angle, (ii) the tributary gradient, (iii) the channel discharges, and (iv) the tributary sediment concentration. Results indicate, in contrast to most research dealing with confluences, that in the presence of intense tributary sediment supply and a small tributary to main channel discharge ratio (0.1), the confluence angle does not have a decisive effect on confluence morphology. Adjustments to the tributary channel gradient yielded the same results. A reoccurring range of depositional geomorphic units was observed where a deposition cone transitioned to a bank-attached bar. The confluence morphology and tributary channel gradient rapidly adjusted, tending towards an equilibrium state to accommodate both water discharges and the sediment load from the tributary. Statistical analyses demonstrated that confluence morphology was controlled by the combined channel discharge and the depositional or erosional extents by the sediment concentration. Applying the conclusions drawn from lowland confluence dynamics could misrepresent depositional and erosional patterns and the related flood hazard at mountain river confluences. [ABSTRACT FROM AUTHOR]

Glacial retreat uncovers large bodies of unconsolidated sediment that are prone to erosion. However, our knowledge of overland flow (OF) generation and sediment transport on moraines that have recently become ice-free is still limited. To investigate how the surface characteristics of young moraines affect OF and sediment transport, we installed five bounded runoff plots on two moraines of different ages in a proglacial area of the Swiss Alps. On each plot we conducted three sprinkling experiments to determine OF characteristics (i.e., total OF and peak OF flow rate) and measured sediment transport (turbidity, sediment concentrations, and total sediment yield). To determine and visualize where sediment transport takes place, we used a fluorescent sand tracer with an afterglow as well as ultraviolet (UV) and light-emitting diode (LED) lamps and a high-resolution camera. The results highlight the ability of this field setup to detect sand movement, even for individual fluorescent sand particles (300–500 µm grain size), and to distinguish between the two main mechanisms of sediment transport: OF-driven erosion and splash erosion. The higher rock cover on the younger moraine resulted in longer sediment transport distances and a higher sediment yield. In contrast, the higher vegetation cover on the older moraine promoted infiltration and reduced the length of the sediment transport pathways. Thus, this study demonstrates the potential of the use of fluorescent sand with an afterglow to determine sediment transport pathways as well as the fact that these observations can help to improve our understanding of OF and sediment transport processes on complex natural hillslopes. [ABSTRACT FROM AUTHOR]

The effects of desert dust upon climate and ecosystems depend strongly on its particle size and size-resolved mineralogical composition. However, there is very limited quantitative knowledge on the particle size and composition of the parent sediments along with their variability within dust-source regions, particularly in dust emission hotspots. The lower Drâa Valley, an inland drainage basin and dust hotspot region located in the Moroccan Sahara, was chosen for a comprehensive analysis of sediment particle size and mineralogy. Different sediment type samples (n= 42) were collected, including paleo-sediments, paved surfaces, crusts, and dunes, and analysed for particle-size distribution (minimally and fully dispersed samples) and mineralogy. Furthermore, Fe sequential wet extraction was carried out to characterise the modes of occurrence of Fe, including Fe in Fe (oxyhydr)oxides, mainly from goethite and hematite, which are key to dust radiative effects; the poorly crystalline pool of Fe (readily exchangeable ionic Fe and Fe in nano-Fe oxides), relevant to dust impacts upon ocean biogeochemistry; and structural Fe. Results yield a conceptual model where both particle size and mineralogy are segregated by transport and deposition of sediments during runoff of water across the basin and by the precipitation of salts, which causes a sedimentary fractionation. The proportion of coarser particles enriched in quartz is higher in the highlands, while that of finer particles rich in clay, carbonates, and Fe oxides is higher in the lowland dust emission hotspots. There, when water ponds and evaporates, secondary carbonates and salts precipitate, and the clays are enriched in readily exchangeable ionic Fe, due to sorption of dissolved Fe by illite. The results differ from currently available mineralogical atlases and highlight the need for observationally constrained global high-resolution mineralogical data for mineral-speciated dust modelling. The dataset obtained represents an important resource for future evaluation of surface mineralogy retrievals from spaceborne spectroscopy. [ABSTRACT FROM AUTHOR]

Extreme floods occurred from 4 to 6 August 2023 in Slovenia causing three casualties and causing total direct and indirect damage, including post-disaster needs according to the Post-Disaster Needs Assessment (PDNA), close to EUR 10 billion. The atypical summer weather conditions combined with the high air and sea temperatures in the Mediterranean and the high soil moisture led to the most extreme flood event in Slovenia in recent decades. The return periods of both daily and sub-daily precipitation extremes and peak discharges reached 250–500 years, and the runoff coefficient of a typical torrential and mostly forested mesoscale catchment was around 0.5. In addition, flooding, soil erosion, mass movements and river sediment transport processes caused major damage to buildings (more than 12 000 houses) and diverse infrastructure. [ABSTRACT FROM AUTHOR]

Post-wildfire floods are receiving greater attention as wildland-urban interfaces become more common and as catastrophic wildfires have increased in frequency. Sediment sourcing, transport, and deposition in the post-fire environment receive attention due to the severity of risk caused by debris flows and concentrated sediment flood flows. This study provides a series of sediment model predictions based on MUSLE and the WARSSS suite of models that included: ERMIT, BANCS, and FLOWSED/POWERSED for the 2019 Museum Fire (809 Ha of steep slope Pinus ponderosa forest in the Spruce Wash watershed). A comparison is provided for the internet-based WEPPcloud post-fire sediment model. Empirical evidence from four floods in 2021 indicated 9,900 Mg of sediment yield to city of Flagstaff neighborhoods, the WEPP model estimated 3870 Mg/year, MUSLE predicted 4860 Mg/year (based on the four events), and the WARSSS suite of models predicted 4630 Mg/year. Both WEPP and WARSSS estimated more sediment yield from channels than hillslope (51%/49% and 60%/40% respectively) though the spatial patterns differ between the models. Sediment mitigation structures, or "work areas", are discussed as real-world applications of sediment forecasting for reducing downstream impacts. Continued revisions of sediment forecasts, based on case studies such as this one, can provide managers and policy makers with tools for risk mitigation and emergency management. [ABSTRACT FROM AUTHOR]

The measurement of cosmogenic nuclide (CN) concentrations in riverine sediment has provided breakthroughs in our understanding of landscape evolution. Yet, linking this detrital CN signal and relief evolution is based on hypotheses that are not easy to verify in the field. Models can be used to explore the statistics of CN concentrations in sediment grains. In this work, we present a coupling between the landscape evolution model Cidre and a model of the CN concentration in distinct grains. These grains are exhumed and detached from the bedrock and then transported in the sediment to the catchment outlet with temporary burials and travel according to the erosion–deposition rates calculated spatially in Cidre. The concentrations of various CNs can be tracked in these grains. Because the CN concentrations are calculated in a limited number of grains, they provide an approximation of the whole CN flux. Therefore, this approach is limited by the number of grains that can be handled in a reasonable computing time. Conversely, it becomes possible to record part of the variability in the erosion–deposition processes by tracking the CN concentrations in distinct grains using a Lagrangian approach. We illustrate the robustness and limitations of this approach by deriving the catchment-average erosion rates from the mean 10Be concentration of grains leaving a synthetic catchment and comparing them with the erosion rates calculated from sediment flux, for different uplift scenarios. We show that the catchment-average erosion rates are approximated to within 5 % uncertainty in most of the cases with a limited number of grains. [ABSTRACT FROM AUTHOR]

The prediction of river morphology evolution is very complicated, especially in the case of mountain and Piedmont rivers with complex morphologies, steep slopes, and heterogeneous grain sizes. The Lac des Gaves (LDG) reach, located within the Gave de Pau River in the Hautes-Pyrénées department, France, has precisely the complex morphological characteristics mentioned above. This reach has gone through severe sediment extractions for over 50 years, leading to the construction of two weirs for riverbed stabilisation. Two large floods resulted in changes in the LDG's hydromorphological characteristics as it went from a single channel river section to a braided river reach. In this study, a 2D hydromorphological model is developed with the TELEMAC-MASCARET system to reproduce the evolution of the channel following a flood that occurred in 2018. The model's validity is assessed by comparing the simulated topographic evolution to the observed one. The results reveal the challenge to choose well-fitted sediment transport equations and friction laws that would make it possible to reproduce such complex morphology. Although the exact localisation of the multiple channels forming the braided nature of the LDG was challenging to reproduce, our model was able to provide reliable volumetric predictions as it reproduces the filling of the LDG correctly. The influence of the two weirs on the river's current and future morphology is also studied. The aim is to provide decision-makers with more reliable predictions to design suitable restoration measures for the LDG reach. [ABSTRACT FROM AUTHOR]

Coastal wetland vegetation modulates water flow by exerting drag, which has important implications for sediment transport and geomorphic dynamics. This vegetation effect on flow is commonly represented in hydrodynamic models by approximating the vegetation as an array of vertical cylinders or increased bed roughness. However, this simple approximation may not be valid in the case of Rhizophora mangroves that have complicated three-dimensional root structures. Here, we present a new model to represent the impact of Rhizophora mangroves on flow in hydrodynamic models. The model explicitly accounts for the effects of the three-dimensional root structures on mean flow and turbulence as well as for the effects of two different length scales of vegetation-generated turbulence characterized by stem diameter and root diameter. The model employs an empirical model for the Rhizophora root structures that can be applied using basic vegetation parameters (mean stem diameter and tree density) without rigorous measurements of the root structures. We tested the model against the flows measured by previous studies in a model mangrove forest in the laboratory and an actual mangrove forest in the field, respectively. We show that, compared with the conventional approximation using an array of cylinders or increased bed roughness, the new model significantly improves the predictability of velocity, turbulent kinetic energy, and bed shear stress in Rhizophora mangrove forests. Overall, the presented new model offers a more realistic but feasible framework for simulating flows in Rhizophora mangrove forests with complex root structures using hydrodynamic models. [ABSTRACT FROM AUTHOR]

A field campaign was carried out at a sheltered sandy beach with the aim of gaining new insights into the driving processes behind sheltered beach morphodynamics. Detailed measurements of the local hydrodynamics, bed level changes and sediment composition were collected at the man-made beach on the leeside of the barrier island Texel, bordering the Marsdiep basin that is part of the Dutch Wadden Sea. The data set consists of (1) current, wave and turbidity measurements from a dense cross-shore array and a 3 km alongshore array, (2) sediment composition data from beach surface samples, (3) high-temporal-resolution RTK-GNSS beach profile measurements, (4) a pre-campaign spatially covering bathy-topo map and (5) meteorological data. This manuscript outlines how these measurements were set up and how the data have been processed, stored and can be accessed. The novelty of this data set lies in the detailed approach to resolve forcing conditions on a sheltered beach, where sediment transport is governed by a strong interplay between tidal currents, waves and bed composition, primarily due to the low-energy (near-threshold) forcing. [ABSTRACT FROM AUTHOR]

Reservoirs represent a globally significant source of the greenhouse gas methane (CH4), which is emitted via different emission pathways. In some reservoirs, reservoir flushing is employed as a sediment management strategy to counteract growing sediment deposits that threaten reservoir capacity. Reservoir flushing utilizes the eroding force of water currents during water level drawdown to mobilize and transport sediment deposits through the dam outlet into the downstream river. During this process, CH4 that is stored in the sediment can be released into the water and degas to the atmosphere, resulting in CH4 emissions. Here, we assess the significance of this CH4 emission pathway and compare it to other CH4 emission pathways from reservoirs. We measured seasonal and spatial CH4 concentrations in the sediment of Schwarzenbach Reservoir, providing one of the largest datasets on CH4 pore water concentrations in freshwater systems. Based on this dataset we determined CH4 fluxes from the sediment and estimated potential CH4 emissions due to reservoir flushing. CH4 emissions due to one flushing operation can constitute 7 %–14 % of the typical annual CH4 emissions from Schwarzenbach Reservoir, whereby the amount of released CH4 depends on the seasonal timing of the flushing operation and can differ by a factor of 2. Larger flushing events that mobilize deeper sediment layers lead to non-linear increases in CH4 mobilization. This suggests that regular flushing of smaller sediment layers releases less CH4 than removal of the same sediment volume in fewer flushing events of thicker sediment layers. However, additional indirect CH4 emissions pathways contributing to the total CH4 emissions may vary with the flushing operation. In other reservoirs with higher sediment loadings than Schwarzenbach Reservoir, reservoir flushing could cause substantial CH4 emissions, especially when flushing operations are conducted frequently. Our study recognizes CH4 emissions due to reservoir flushing as an important pathway, identifies potential management strategies to mitigate these CH4 emissions and emphasizes the need for further research. [ABSTRACT FROM AUTHOR]

Reliable forecasts of building damage due to debris flows may provide situational awareness and guide land and emergency management decisions. Application of debris-flow runout models to generate such forecasts requires combining hazard intensity predictions with fragility functions that link hazard intensity with building damage. In this study, we evaluated the performance of building damage forecasts for the 9 January 2018 Montecito postfire debris-flow runout event, in which over 500 buildings were damaged. We constructed forecasts using either peak debris-flow depth or volume flux as the hazard intensity measure and applied each approach using three debris-flow runout models (RAMMS, FLO-2D, and D-Claw). Generated forecasts were based on combining multiple simulations that sampled a range of debris-flow volume and mobility, reflecting typical sources and magnitude of pre-event uncertainty. We found that only forecasts made with volume flux and the D-Claw model could correctly forecast the observed number of damaged buildings and the spatial patterns of building damage. However, the best forecast only predicted 50% of the observed damaged buildings correctly and had coherent spatial patterns of incorrectly forecast building damage (i.e., false positives and false negatives). These results indicate that forecasts made at the building level reliably reflect the spatial pattern of damage, but do not support interpretation at the individual building level. We found the event size strongly influences the number of damaged buildings and the spatial pattern of debrisflow depth and velocity. Consequently, future research on the link between precipitation and the volume of sediment mobilized may have the greatest effect on reducing uncertainty in building damage forecasts. Finally, because we found that both depth and velocity are needed to forecast building damage, comparing debris flow models against spatially distributed observations of building damage is a more stringent test for model fidelity than comparison against the extent of debris-flow runout. [ABSTRACT FROM AUTHOR]

In the process of reworking sediments and thus shaping biogeochemical processes, marine bottom dwelling animals are thought to play a pivotal role in many benthic environments. This aspect of bioturbation (particle reworking) is often partitioned into short distance local and non-local transport acting over relatively longer distances. Here we document that subsurface peaks, such as those typically attributed to biological particle transport in sediments, may equally be generated by otter boards in bottom trawling fishery. Boards can generate tracer peaks whereby they scoop sediment from the surface, flip it over and deposit it onto the adjacent sea floor. These peaks are indistinguishable from those (presumably) generated by benthic fauna in a process whereby. We demonstrate this for the particle tracer chlorophyll a in silty sand from the Western Baltic Sea with fauna that generally does not burrow deep in a global comparison. Our inability to distinguish the driving processes generating the peaks indicates limits to our understanding of magnitude and spatial extend of bioturbation traces in this environment. It also poses a problem to the assessment of fishery resource use and benthic ecosystem services. However, we can clearly identify macrofauna and not otter boards as the cause for peaks at the sites investigated here. [ABSTRACT FROM AUTHOR]

Landscape evolution models (LEMs) are useful for understanding how large-scale processes and perturbations influence the development of the surface of the Earth and other planets. With their increasing sophistication and improvements in computational power, they are finding greater uptake in analyses at finer spatial and temporal scales. For many LEMs, the land surface is represented by a grid of regularly spaced and sized grid cells, or pixels, referred to as a digital elevation model (DEM), yet despite the importance of the DEM to LEM studies, there has been little work to understand the influence of grid cell size (i.e. resolution) on model behaviour. This is despite the choice of grid cell size being arbitrary for many studies, with users needing to balance detail with computational efficiency. Using the Morris method (MM) for global sensitivity analysis, the sensitivity of the CAESAR-Lisflood LEM to the grid cell size is evaluated relative to a set of influential user-defined parameters, showing that it had a similar level of influence as a key hydrological parameter and the choice of sediment transport law. Outputs relating to discharge and sediment yields remained stable across different grid cell sizes until the cells became so large that the representation of the hydrological network degraded. Although total sediment yields remained steady when changing the grid cell sizes, closer analysis revealed that using a coarser grid resulted in it being built up from fewer yet more geomorphically active events, risking outputs that are "the right answer but for the wrong reasons". These results are important considerations for modellers using LEMs and the methodologies detailed provide solutions to understanding the impacts of modelling choices on outputs. [ABSTRACT FROM AUTHOR]

Prediction of bed load sediment transport rates in rivers is a notoriously difficult problem due to inherent variability in river hydraulics and channel morphology. Machine learning (ML) offers a compelling approach to leverage the growing wealth of bed load transport observations towards the development of a data-driven predictive model. We present an artificial neural network (ANN) model for predicting bed load transport rates informed by 8117 measurements from 134 rivers. Inputs to the model were river discharge, flow width, bed slope, and four bed surface sediment sizes. A sensitivity analysis showed that all inputs to the ANN model contributed to a reasonable estimate of bed load flux. At individual sites, the ANN model was able to reproduce observed sediment rating curves with a variety of shapes without site-specific calibration. This ANN model has the potential to be broadly applied to predict bed load fluxes based on discharge and reach properties alone. [ABSTRACT FROM AUTHOR]

Glaciers expel sediment as they melt, in addition to ice and water. As a result, changing glacier dynamics and melt produce changes to glacier erosion and sediment discharge, which can impact the landscape surrounding retreating glaciers, as well as communities and ecosystems downstream. Currently, numerical models that transport subglacial sediment on sub-hourly to decadal scales are one-dimensional, usually along a glacier's flow line. Such models have proven useful in describing the formation of glacial landforms, the impact of sediment transport on glacier dynamics, and the interactions among climate, glacier dynamics, and erosion. However, these models omit the two-dimensional spatial distribution of sediment and its impact on sediment connectivity – the movement of sediment between its detachment in source areas and its deposition in sinks. Here, we present a numerical model that fulfills a need for predictive frameworks that describe subglacial sediment discharge in two spatial dimensions (x and y) over time. SUGSET_2D evolves a two-dimensional subglacial till layer in response to bedrock erosion and changing sediment transport conditions. Numerical experiments performed using an idealized alpine glacier illustrate the heterogeneity in sediment transport and bedrock erosion below the glacier. An increase in sediment discharge follows increased glacier melt, as has been documented in field observations and other numerical experiments. We also apply the model to a real alpine glacier, Griesgletscher in the Swiss Alps, where we compare outputs with annual measurements of sediment discharge. SUGSET_2D accurately reproduces the general quantities of sediment discharge and the year-to-year sediment discharge pattern measured at the glacier terminus. The model's ability to match the measured data depends on the tunable sediment grain size parameter, which controls subglacial sediment transport capacity. Smaller grain sizes allow sediment transport to occur in regions of the bed with reduced water flow and channel size, effectively increasing sediment connectivity into the main channels. The model provides the essential components of modeling subglacial sediment discharge on seasonal to decadal timescales and reveals the importance of including spatial heterogeneities in water discharge and sediment transport in both the x and y dimensions in evaluating sediment discharge. [ABSTRACT FROM AUTHOR]

Submarine channels conveying sediment gravity flows are often topographically confined, but the effect of confinement width on channel morphodynamics is incompletely understood. We use physical experiments and a reduced-complexity model to investigate the effects of confinement width (B) and the inflow-to-sediment discharge ratio (Qin/Qs) on the evolution of submarine braided channels. The results show that a larger confinement width results in increased active braiding intensity (BIA) and that BIA takes longer to stabilize (i.e., a longer critical time; tc). At a fixed confinement width, a higher Qin/Qs slightly decreases the BIA. Digital elevation models of difference (DoD) of the experiments allow measurement of the morphological active width (Wa) of the submarine channels and the bulk morphological change (Vbulk) within an experiment, defined as the sum of total erosion and deposition. We find that Wa and Vbulk are proportional to B. We further confirm that BIA is proportional to both dimensionless sediment–stream power (ω**) and dimensionless stream power (ω*). These trends are consistent for submarine braided channels both with and without confinement width effects. Furthermore, we built a reduced-complexity model (RCM) that can simulate flow bifurcation and confluence of submarine braided channels. The simulated flow distribution provides reliable predictions of flow depth and sediment transport rate in the experiments. Using kernel density estimation (KDE) to forecast the probability and trends of cross-sectional flow distribution and corresponding BIA under extreme events, we find that skewness of the flow distribution decreases as discharge increases. The development of braided submarine channels, shown here to extend to conditions of topographic confinement, suggests that factors not modeled here (e.g., fine sediment) may be necessary to explain the abundance of single-thread submarine channels in nature. [ABSTRACT FROM AUTHOR]

Streams convey coarse-clastic sediments towards coasts, where interactions with deltaic and coastal processes determine their resultant sedimentology and geomorphology. Extracting hydroclimatic signals from such environments is a desired goal, and therefore studies commonly rely on interpreting available paleoclimatic proxy data, but the direct linking of depositional and geomorphic processes with the hydroclimate remains obscure. This is a consequence of the challenge of linking processes that are often studied separately and span across large spatial and temporal scales, including synoptic-scale hydroclimatic forcing, streamflows, water body hydrodynamics, fluvial and coastal sediment transport, and sedimentation. Here, we explore this chain of connected processes in the unique setting of the Dead Sea basin, where present-day hydroclimatology is closely tied with geomorphic evolution and sediment transport of streams and coasts that rapidly respond to lake-level fall. We use a 5-year (2018–2022) rich dataset of (i) high-resolution synoptic-scale circulation patterns, (ii) continuous wind-wave and rain–flood records, and (iii) storm-scale fluvial and coastal sediment transport of "smart" and marked boulders. We show the significance of Mediterranean cyclones in the concurrent activation of fluvial (floods) and coastal (wind waves) sediment conveyors. These synoptic-scale patterns drive the westerlies necessary for (i) delivering the moisture across the Judean desert, which is transformed into floods, and at the same time, (ii) the coeval, topographically funneled winds that turn into surface southerlies (>10 m s -1) along the Dead Sea rift valley. During winter, these mesoscale southerlies generate 10–30 high-amplitude, northward-propagating storm waves, with <4 m wave heights. Such waves transport cobbles for hundreds of meters alongshore, northward and away from the supplying channel mouths. Four to nine times per winter the rainfall generated by these atmospheric patterns is capable of generating floods that reach the stream mouths, delivering poorly sorted, coarse gravel. This usually occurs during the decay of the associated storm waves. This gravel is dispersed alongshore by waves during subsequent storms. As storm waves dominate and are >5 times more frequent than flash floods, coarse-clastic beach berms and fan deltas are deposited preferentially north of the delivering channel mouths. This asymmetric depositional architecture, controlled by the regional hydroclimatology, is identified for both the modern and late Pleistocene coast and delta environments, implying that the dominance of present-day Mediterranean cyclones also persisted in the region during the late Pleistocene when Lake Lisan occupied the basin. [ABSTRACT FROM AUTHOR]

Sediment burial dating using optically stimulated luminescence (OSL) is a well-established tool in geochronology. An important but often inapplicable requirement for its successful use is that the OSL signal is sufficiently reset prior to deposition. However, subaqueous bleaching conditions during fluvial transport are vastly understudied; for example the effect of turbidity and sediment mixing on luminescence bleaching rates is only poorly established. The possibility that slow bleaching rates may dominate under certain transport conditions led to the concept that OSL could be used to derive sediment transport histories. The feasibility of this concept is still to be demonstrated, and experimental set-ups are still to be tested. Our contribution to this scientific challenge involves subaquatic bleaching experiments, in which we suspend saturated coastal sand of Miocene age in a circular flume and illuminate it for discrete time intervals with natural light. We record the in situ energy flux density received by the suspended grains in the UV-NIR frequency range by using a broadband spectrometer with a submersible probe. Our analysis includes pre-profiling of each sample following a polymineral multiple signal (PMS) protocol. Using the PMS, the quartz-dominated, blue-stimulated luminescence signal at 125 ∘ C (BSL-125) decays slower than the K-feldspar-dominated, infrared-stimulated luminescence signal at 25 ∘ C (IR-25) even under subaerial conditions. The BSL-125 from purified quartz shows the opposite behaviour, which renders the PMS unreliable in our case. We find a negative correlation between suspended-sediment concentration and bleaching rate for all the measured signals. For outdoor bleaching experiments we propose to relate the measured luminescence dose to the cumulative received irradiance rather than to the bleaching time. Increases in the sediment concentration lead to a stronger attenuation of the UV–blue compared to the red–NIR wavelength. This attenuation thereby follows an exponential decay that is controlled by the sediment concentration and a wavelength-dependent decay constant, λ. As such λ could potentially be used in numerical models of luminescence signal resetting in turbid suspensions. [ABSTRACT FROM AUTHOR]

Mass movements play an important role in landscape evolution of high mountain areas such as the Himalayas. Yet, establishing numerical age control and reconstructing transport dynamics of past events is challenging. To fill this research gap, we bring Optically Stimulated Luminescence (OSL) dating to the test in an extremely challenging environment: the Pokhara Valley in Nepal. This is challenging for two reasons: i) the OSL sensitivity of quartz, typically the mineral of choice for dating sediments younger than 100 ka, is poor, and ii) highly rapid and turbid conditions during mass movement transport hamper sufficient OSL signal resetting prior to deposition, which eventually results in age overestimation. Here, we first assess the applicability of single-grain feldspar dating of medieval mass movement deposits catastrophically emplaced in the Pokhara Valley. Second, we exploit the poor bleaching mechanisms to get insight into the sediment dynamics of this paleo-mass movement through bleaching proxies. The Pokhara valley is a unique setting for our case-study, considering the availability of an extensive independent radiocarbon dataset as a geochronological benchmark. Single-grain infrared stimulated luminescence signals were measured at 50 °C (IRSL-50) and post-infrared infrared stimulated luminescence signals at 150 °C (pIRIR-150). Our results show that the IRSL-50 signal is better bleached than the pIRIR-150 signal. A bootstrapped Minimum Age Model (bMAM) is applied to retrieve the youngest subpopulation to estimate the paleodose. However, burial ages calculated with this paleodose overestimate the radiocarbon ages by an average factor of ~23 (IRSL-50) and ~72 (pIRIR-150), showing that dating of the Pokhara Formation with a single-grain approach was not successful for most samples. Some samples, however, only slightly overestimate the true emplacement age and thus could be used for a rough age estimation. Large inheritances in combination with the scatter in the single-grain dose distributions show that the sediments have been transported under extremely limited light exposure prior to deposition, which is consistent with the highly turbid nature of the sediment laden flood and debris flows depositing the Pokhara gravels. To investigate the sediment transport dynamics in more detail, we studied three bleaching proxies: the percentage of grains in saturation 2D0 criteria, the percentage of best-bleached grains (2σ range of bMAM-De) and the overdispersion (OD). None of the three bleaching proxies indicate a spatial relationship with run-out distance of the mass movement deposits. We interpret this as evidence for the lack of bleaching during transport, which reflects the catastrophic nature of the event. While not providing reliable burial ages of the Pokhara mass movement deposits, single-grain feldspar dating can potentially be used as an age range finder method. Our approach shows the potential of luminescence techniques to provide insights in sediment transport dynamics of extreme and rare mass movement events in mountainous region. [ABSTRACT FROM AUTHOR]

The recession of Alpine glaciers causes an increase in the extent of proglacial areas that leads to changes in the water and sediment balance morphodynamics and sediment transport. Although the processes occurring in proglacial areas are relevant not only from a scientific point of view but also for the purpose of climate change adaptation, there is a lack of studies on the continuous monitoring and multitemporal characterization of these areas. This work offers a multidisciplinary approach that merges the contributions of different scientific disciplines such as hydrology, geophysics, geomatics and water engineering to characterise the Rutor glacier and its proglacial area. We surveyed the glacier and its proglacial area since 2020 with both uncrewed (drone) and crewed aerial photogrammetric flights; we determined the bathymetry of the most downstream proglacial lake and the thickness of the sediments deposited on its bottom. Water depth at four different locations within the hydrographic network of the proglacial area and the bedload at the glacier snout were continuously monitored. The synergy of our approach enables the characterisation, monitoring and understanding of a set of complex and interconnected processes occurring in a proglacial area. [ABSTRACT FROM AUTHOR]

Alpine rivers have experienced considerable changes in channel morphology over the last century. Natural factors and human disturbance are the main drivers of changes in channel morphology that modify natural sediment and flow regimes at local, catchment, and regional scales. In glaciated catchments, river sediment loads are likely to increase due to increasing snowmelt and glacier melt runoff, facilitated by climate change. Additionally, channel erosion and depositional dynamics and patterns are influenced by sediment delivery from hillslopes and sediment in the forefields of retreating glaciers. In order to reliably assess the magnitudes of the channel-changing processes and their frequencies due to recent climate change, the investigation period needs to be extended to the last century, ideally back to the end of the Little Ice Age. Moreover, a high temporal resolution is required to account for the history of changes in channel morphology and for better detection and interpretation of related processes. The increasing availability of digitised historical aerial images and advancements in digital photogrammetry provide the basis for reconstructing and assessing the long-term evolution of the surface, in terms of both planimetric mapping and the generation of historical digital elevation models (DEMs). The main issue of current studies is the lack of information over a longer period. Therefore, this study contributes to research on fluvial sediment changes by estimating the sediment balance of a main Alpine river (Fagge) in a glaciated catchment (Kaunertal, Austria) over 19 survey periods from 1953 to 2019. Exploiting the potential of historical multi-temporal DEMs combined with recent topographic data, we quantify 66 years of geomorphic change within the active floodplain, including erosion, deposition, and the amounts of mobilised sediment. Our study focuses on a proglacial river that is undergoing a transition phase, resulting from an extensive glacier retreat of approximately 1.8 km. This has led to the formation of new channel networks and an overall negative cumulative sediment balance for the entire study area. We found that high-magnitude meteorological and hydrological events associated with local glacier retreats have a significant impact on the sediment balance. The gauge record indicates an increase in such events, as well as in runoff and probably in sediment transport capacity. Despite this, the sediment supply has declined in the last decade, which can be attributed to a lower contribution of the lateral moraines coupled to the channel network and less sediment sourced from the melting Gepatsch Glacier as evidenced by roches moutonnées exposed in the current/most recent forefield. Nonetheless, we observed significant erosion in the tributary, leading to the transport of sediment downstream. Overall, this study enhances our understanding of the complexity of sediment dynamics in proglacial rivers across various spatial and temporal scales and their relationship to climate change factors. [ABSTRACT FROM AUTHOR]

As humans continue to inhabit and modify river deltas, the natural processes governing material transport through these landscapes are altered. Two common engineering projects undertaken on deltas are the dredging of channels to enable shipping and the construction of embankments to reduce flooding. While the impact of these topographic modifications has been studied at a local level for specific sites, there is a gap in our generalized understanding of how these landscape modifications impact material transport. To narrow this gap, we conduct exploratory numerical modeling to develop deltaic landscapes with different input sediment compositions, modify their topography to mimic dredging and embankments, simulate different flow conditions, and then model the transport of passive particles. We find that human modification of topography lowers hydrological connectivity by reducing the area visited by fluvial inputs. The amount of time particles spend within the delta is reduced by the construction of polders and is lengthened by dredging. Material buoyancy has a greater impact on nourishment areas and exposure times than flow regime or topographic modification, with positively buoyant particles spending longer and visiting a greater area of the delta than neutral and negatively buoyant material. The results of this study can help guide the design of future engineering projects by providing estimates of their likely impact on transport processes. [ABSTRACT FROM AUTHOR]

We develop a phenomenological model of suspended sediment transport on the basis of data acquired in the Capesterre river, which drains a small tropical catchment in Guadeloupe. The model correctly represents the concentration of suspended sediment during floods, provided that the relation between concentration and water level forms a counterclockwise loop. In the model, the properties of the sediment and of the river are all lumped into four parameters: a settling velocity related to the size of the suspended sediment, a threshold water level which acts as a proxy for the threshold shear stress, a characteristic entrainment rate, and a dimensionless exponent. The value of the parameters changes from one flood to the next, probably reflecting changes in the characteristics of the river and the fine sediment. Finally, a test of the model against data acquired in a small catchment in the French Alps suggests that the model is versatile enough to be used in diverse hydrological settings. [ABSTRACT FROM AUTHOR]

Floods are among the most common natural disasters responsible for severe damages and human losses. Combining numerical modelling with user-friendly tools for geographically referenced data has been adopted to increase preparedness and reduce vulnerabilities. This paper describes the numerical dataset of hydrometric variables that characterize a flood event occurred in February 2016 in the Portuguese Águeda river, shortly defined as Agueda.2016Flood. The dataset was numerically produced and managed through the RiverCure Portal, a collaborative web platform connected to a validated shallow-water model featuring modelled dynamic bed geometries and sediment transport. The dataset Agueda.2016Flood can be used as a starting point to design other experiments and tools, and to learn and apply the proposed approach by directly using the RiverCure Portal. This dataset includes modelled hydrodynamic data (output data) and the topographic, geometrical, land-use and hydrologic data (input data) necessary to carry out the numerical simulation of the flood event. [ABSTRACT FROM AUTHOR]

The 2008 Wenchuan earthquake triggered rapid local geomorphic changes, shifting abundant material through exogenic processes and generating vast amounts of loose material. The substantial material movement increased the geohazard (flash floods, landslides and debris flows) risks induced by extreme precipitation in the area. Intervention measures such as check dams, levees and vegetated slopes have been constructed in specific locations to reduce sediment transport and thereby mitigate the impact of ensuing geohazards. This study assessed the short–medium-term effects of interventions, including multiple control measures, in a post-earthquake mountainous region. Taking the Xingping valley as an example, we used CAESAR-Lisflood, a two-dimensional landscape evolution model, to simulate three scenarios, unprotected landscape, present protected landscape and enhanced protected landscape, between 2011 and 2013. We defined two indices to assess the intervention effects of the three scenarios by comparing the geomorphic changes and sediment yields. The results show that the mitigation measures are effective, especially the geotechnical engineering efforts in combination with ecological engineering in the upstream area. The spatial patterns of erosion and deposition change considerably due to the intervention measures. Additionally, the effectiveness of each intervention scenario shows a gradual decline over time, mainly due to the reduction in the reservoir storage capacity. The enhanced scenario performs better than the present one, with a more gradual downward trend of effectiveness. The simulation results evaluated the ability and effectiveness of comprehensive control measures and will support optimal mitigation strategies. [ABSTRACT FROM AUTHOR]

Suspended sediment is an integral part of riverine transport and functioning that has been strongly altered during the Anthropocene due to the overwhelming human pressure on soils, sediments, and the water cycle. Understanding the controls of changing suspended sediment in rivers is therefore vital for effective management strategies. Here we present results from a trend analysis of suspended sediments covering 62 monitoring stations along the German waterways (catchment sizes range between 2000 and 160 000 km 2) with more than 440 000 water samples taken between 1990 and 2010. Based on daily monitoring of suspended sediment concentration (SSC), we found significant declines in mean annual SSC and annual suspended sediment loads (SSLs) at 49 of 62 monitoring stations totaling - 0.92 mg L -1 yr -1. At some stations decreases during the 20 years represent up to 50 % of the long-term average SSC. Significant decreases in SSC are associated with declining SSL despite an increase in sheet and rill erosion by almost 150 % derived from modeling due to changes in land use and management as well as rainfall erosivity. The contemporary suspended sediment loads of the Rhine at the German–Dutch border are approaching the natural base level of ∼1 Mt yr -1 , which was reached by the Rhine during the mid-Holocene when the suspended sediment load was adjusted to the Holocene climatic conditions and before the onset of increased loads due to human-induced land use changes in the Rhine catchment. At this point we can only speculate regarding potential reasons for a decline in sediment transport in larger rivers despite erosion increase. We argue that increased sediment retention in upstream headwaters is presumably the major reason for declining SSC in the large river channels studied. [ABSTRACT FROM AUTHOR]

Detrital 10Be from continental river sands or submarine sediments has been extensively used to determine the average long-term denudation rates of terrestrial catchments, based on the assumption that the rate of cosmogenic nuclide production by the interaction of source rocks with cosmic radiation balances out the loss of these nuclides by surface denudation. However, the 10Be signal recorded in sediments may be affected at the source by the response time of mountainous catchments to high-frequency forcings. In addition, transient sediment storage in piedmonts, alluvial plains and lakes or near the coast may also induce a difference between the erosive signal and its record in the sedimentary sink. Consequently, a significant part of the signal recorded in shallow-water sediments can be lost, as deep marine sediments may simultaneously record a signal coming from newly eroded source rocks along with one coming from the destabilization of previously deposited sediments. In this paper, we use the landscape evolution model Badlands to simulate erosion, deposition and detrital 10Be transfer from a source-to-sink sedimentary system (the Var River catchment, southern French Alps) over the last 100 kyr. We first compare model-based denudation rates with the ones that would be extracted from the 10Be record of local continental sediments (equivalent to river sands) and from sediments deposited offshore over time in order to examine if this record provides an accurate estimate of continental denudation rates. Then, we examine which conditions (precipitation rate, flexure, ice cover) satisfy published measured river incision rates and 10Be concentration in submarine sediments. Our results, based on the Var catchment cosmic ray exposure dating and modelling indicate that, while river sands do accurately estimate the average denudation rate of continental catchments, this is much less the case for deep submarine sediments. We find that deep-sea sediments have a different and often much smoother 10Be signature than continental ones and record a significant time lag with respect to imposed precipitation rate changes, representing the geomorphological response of the margin. A model which allows us to fit both measured 10Be concentration in marine sediments and river incision rates on land involves an increase in precipitation rates from 0.3 to 0.7 m yr -1 after 20 ka, suggesting more intense precipitation starting at the end of the Last Glacial Maximum. [ABSTRACT FROM AUTHOR]

Extracting accurate bathymetric information from clear, shallow waters in complex riverine environments can be challenging, but it is crucial for many applications, such as hydraulic modeling, ecological monitoring, and sediment transport analysis. Multi-view stereo photogrammetry (MVS-SfM) has emerged as a promising technique for acquiring high-resolution bathymetric data from aerial imagery. However, the accuracy of MVS-SfM can be affected by various factors, including water refraction, which can distort the depth measurements. In this study, iterative Dietrich's refraction-correction method is tested for extracting bathymetry from MVS-SfM in a complex riverine environment. Moreover, we proposed a workflow for applying the method using raster data files, which can be more readily available than point clouds. It also compared the obtained results with previous applications and evaluated them using statistical indices and ratios. For this case study, the multicamera refraction method produces bathymetric datasets with accuracies of ~0.019% of the flying height and precisions of ~0.07% of the flying height. This study contributes to increased confidence in exploiting aerial imagery for bathymetric mapping in photogrammetric procedures under field conditions and can facilitate the management and conservation of clear water systems. [ABSTRACT FROM AUTHOR]

The early Cenozoic marine sedimentary record is punctuated by several brief episodes (<200 kyr) of abrupt global warming, called hyperthermals, that have disturbed ocean life and water physicochemistry. Moreover, recent studies of fluvial–deltaic systems, for instance at the Palaeocene–Eocene Thermal Maximum, revealed that these hyperthermals also impacted the hydrologic cycle, triggering an increase in erosion and sediment transport at the Earth's surface. Contrary to the early Cenozoic hyperthermals, the Middle Eocene Climatic Optimum (MECO), lasting from 40.5 to 40.0 Ma, constitutes an event of gradual warming that left a highly variable carbon isotope signature and for which little data exist about its impact on Earth surface systems. In the South Pyrenean foreland basin (SPFB), an episode of prominent deltaic progradation (Belsué–Atarés and Escanilla formations) in the middle Bartonian has been usually associated with increased Pyrenean tectonic activity, but recent magnetostratigraphic data suggest a possible coincidence between the progradation and the MECO warming period. To test this hypothesis, we measured the stable-isotope composition of carbonates (δ13C carb and δ18O carb) and organic matter (δ13C org) of 257 samples in two sections of SPFB fluvial–deltaic successions covering the different phases of the MECO and already dated with magnetostratigraphy. We find a negative shift in δ18O carb and an unclear signal in δ13C carb around the transition from magnetic chron C18r to chron C17r (middle Bartonian). These results allow, by correlation with reference sections in the Atlantic and Tethys, the MECO to be identified and its coincident relationship with the Belsué–Atarès fluvial–deltaic progradation to be documented. Despite its long duration and a more gradual temperature rise, the MECO in the South Pyrenean foreland basin may have led, like lower Cenozoic hyperthermals, to an increase in erosion and sediment transport that is manifested in the sedimentary record. The new data support the hypothesis of a more important hydrological response to the MECO than previously thought in mid-latitude environments, including those around the Tethys. [ABSTRACT FROM AUTHOR]

Following the tradition of modeling fluvial landscape evolution, a novel approach describing glacial erosion based on an empirical stream power law was proposed. This approach differs substantially from well-established process-based models applied to describe glacial erosion in mountain landscapes. Outstanding computational performance but a number of potential limitations compared to process-based models requires extensive testing to evaluate the applicability of this novel approach. In this study, we test the validity of the glacial stream power law and its implementation into a 2-D landform evolution model (OpenLEM) by benchmarking it against a state of the art surface process model based on the integrated second-order shallow-ice approximation (iSOSIA). Despite completely different approaches, OpenLEM and iSOSIA predict similar ice flow patterns and erosion rates for a wide range of climatic conditions without re-adjusting a set of calibrated parameters. This parameter set is valid for full glacial conditions where the entire precipitation is converted to ice but also for an altitude-dependent glacier mass balance as characteristic for most glaciated mountain ranges on Earth. In both models characteristic glacial features, such as overdeepenings, hanging valleys and steps at confluences emerge roughly at the same locations, resulting in a consistent altitude-dependent adjustment of channel slope and relief. Compared to iSOSIA, however, distinctly higher erosion rates occur in OpenLEM at valley flanks during the initial phase of the fluvial to glacial transition. This is mainly due to the simplified description of glacier width and ice surface in OpenLEM. In this respect, we found that the glacial stream power approach cannot replace process-based models such as iSOSIA but is complementary to them by addressing research questions that could not previously be answered due to a lack of computational efficiency. The implementation of the glacial stream power law is primarily suitable for large-scale simulations investigating the evolution of mountain topography in the interplay of tectonics and climate. As coupling glacial and fluvial erosion with sediment transport shows nearly the same computational efficiency as its purely fluvial counterpart, mountain-range-scale simulations at high spatial resolution are not exclusively restricted to the fluvial domain anymore, and a series of exciting research questions can be addressed by this novel approach. [ABSTRACT FROM AUTHOR]

The proportional cover of rubble on reefs is predicted to increase as disturbances increase in intensity and frequency. Unstable rubble can kill coral recruits and impair binding processes that consolidate rubble into a stable substrate for coral recruitment. A clearer understanding of the mechanisms of inhibited coral recovery on rubble requires characterisation of the hydrodynamic conditions that trigger rubble mobilisation. Here, we investigated rubble mobilisation under regular wave conditions in a wave flume and irregular wave conditions in-situ on a coral reef in the Maldives. We examined how changes in near-bed wave orbital velocity influenced the likelihood of rubble motion (e.g., rocking) and transport (by walking, sliding or flipping). Rubble mobilisation was considered as a function of rubble length, branchiness (branched vs. unbranched), and underlying substrate (rubble vs. sand). Rubble was more likely to be transported if ieces were small (4-8 cm) and had no branches, and rubble travelled slightly greater distances (~2 cm) per day on substrates composed of sand than rubble. The effect of near-bed wave orbital velocity on rubble mobilisation was comparable between flume and reef observations. Rubble had a 50% and 90% chance of transport when near-bed wave orbital velocities reached 0.30 m/s and 0.43 m/s, respectively, in the wave flume, and 0.34 m/s and 0.55 m/s, respectively, on the reef. Importantly, the probability of rubble transport per day declined over 3-day deployments in the field, suggesting rubble had settled into more hydrodynamically-stable positions or snagged on the first day of deployment. We expect that settled or snagged rubble may have been mobilised more commonly in locations with higher energy and more variable wave environments. Our results show that rubble beds comprised of small rubble pieces and/or pieces with fewer branches are likely to be more unstable. Such rubble beds are likely to have shorter windows of recovery (stability) between mobilisation events, and thus be good candidates for rubble stabilisation interventions to enhance coral recruitment and binding. [ABSTRACT FROM AUTHOR]